|
45 | 45 | <ul class="current"> |
46 | 46 | <li class="toctree-l1"><a class="reference internal" href="1_requirements.html">Requirements and Quickstart</a></li> |
47 | 47 | <li class="toctree-l1"><a class="reference internal" href="2_theory.html">Theory of FTLE and LCS</a></li> |
48 | | -<li class="toctree-l1 current"><a class="current reference internal" href="#">Numerical Methods</a></li> |
49 | | -<li class="toctree-l1"><a class="reference internal" href="#computational-density-and-comparison">Computational Density and Comparison</a></li> |
| 48 | +<li class="toctree-l1 current"><a class="current reference internal" href="#">Numerical Methods</a><ul> |
| 49 | +<li class="toctree-l2"><a class="reference internal" href="#wall-treatment">Wall Treatment</a></li> |
| 50 | +<li class="toctree-l2"><a class="reference internal" href="#lagrangian-advection">Lagrangian Advection</a></li> |
| 51 | +<li class="toctree-l2"><a class="reference internal" href="#velocity-interpolation">Velocity Interpolation</a></li> |
| 52 | +<li class="toctree-l2"><a class="reference internal" href="#gradient-discretization">Gradient Discretization</a></li> |
| 53 | +<li class="toctree-l2"><a class="reference internal" href="#eigenvalue-solver">Eigenvalue Solver</a></li> |
| 54 | +</ul> |
| 55 | +</li> |
| 56 | +<li class="toctree-l1"><a class="reference internal" href="#computational-density-and-comparison">Computational Density and Comparison</a><ul> |
| 57 | +<li class="toctree-l2"><a class="reference internal" href="#general-tips">General Tips</a></li> |
| 58 | +</ul> |
| 59 | +</li> |
50 | 60 | <li class="toctree-l1"><a class="reference internal" href="4_input.html">Input Data Files</a></li> |
51 | 61 | <li class="toctree-l1"><a class="reference internal" href="5_gui.html">Graphical User Interface (GUI)</a></li> |
52 | 62 | <li class="toctree-l1"><a class="reference internal" href="6_inputdeck.html">Input Deck for Solver Parameters</a></li> |
|
81 | 91 |
|
82 | 92 | <section id="numerical-methods"> |
83 | 93 | <span id="numerical"></span><h1>Numerical Methods<a class="headerlink" href="#numerical-methods" title="Link to this heading"></a></h1> |
| 94 | +<section id="wall-treatment"> |
| 95 | +<span id="wall"></span><h2>Wall Treatment<a class="headerlink" href="#wall-treatment" title="Link to this heading"></a></h2> |
| 96 | +</section> |
| 97 | +<section id="lagrangian-advection"> |
| 98 | +<span id="advc"></span><h2>Lagrangian Advection<a class="headerlink" href="#lagrangian-advection" title="Link to this heading"></a></h2> |
| 99 | +</section> |
| 100 | +<section id="velocity-interpolation"> |
| 101 | +<span id="intp"></span><h2>Velocity Interpolation<a class="headerlink" href="#velocity-interpolation" title="Link to this heading"></a></h2> |
| 102 | +</section> |
| 103 | +<section id="gradient-discretization"> |
| 104 | +<span id="grad"></span><h2>Gradient Discretization<a class="headerlink" href="#gradient-discretization" title="Link to this heading"></a></h2> |
| 105 | +</section> |
| 106 | +<section id="eigenvalue-solver"> |
| 107 | +<span id="eigen"></span><h2>Eigenvalue Solver<a class="headerlink" href="#eigenvalue-solver" title="Link to this heading"></a></h2> |
| 108 | +</section> |
84 | 109 | </section> |
85 | 110 | <section id="computational-density-and-comparison"> |
86 | | -<h1>Computational Density and Comparison<a class="headerlink" href="#computational-density-and-comparison" title="Link to this heading"></a></h1> |
| 111 | +<span id="numcompare"></span><h1>Computational Density and Comparison<a class="headerlink" href="#computational-density-and-comparison" title="Link to this heading"></a></h1> |
| 112 | +<section id="general-tips"> |
| 113 | +<span id="numtips"></span><h2>General Tips<a class="headerlink" href="#general-tips" title="Link to this heading"></a></h2> |
| 114 | +<p>As for your reference, and configured as defaults, the <em>Berkeley LCS Tutorials</em> used <code class="docutils literal notranslate"><span class="pre">RK4</span></code> for advection. |
| 115 | +The velocity fields were interpolated with <code class="docutils literal notranslate"><span class="pre">tricubic-FL</span></code> by them, originating from <span id="id1">[Lekien2005]</span>, which has higher performance by solving a 64*64 linear system using the function values, gradients, and mixed partial derivatives at its eight corners, which is in future development plan for <code class="docutils literal notranslate"><span class="pre">Py3DFTLE</span></code> with high priority. |
| 116 | +Although not detailed, <code class="docutils literal notranslate"><span class="pre">grad_order=2</span></code> was employed by them from the equation, supposing the mesh is sufficiently refined.</p> |
| 117 | +<p>Please always notice that, although providing much better numerical precision and looks cool in papers, high-order methods could be resource-consuming, even several hundred times.</p> |
| 118 | +</section> |
87 | 119 | </section> |
88 | 120 |
|
89 | 121 |
|
|
0 commit comments