Skip to content

Commit 72b1aaa

Browse files
committed
minor
1 parent cf030ec commit 72b1aaa

File tree

2 files changed

+12
-8
lines changed

2 files changed

+12
-8
lines changed

source/3_numerical.rst

Lines changed: 10 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -176,6 +176,10 @@ The continuous velocity field is reconstructed by trilinear interpolation of the
176176

177177

178178

179+
**Tricubic by F. Lekien** ``I.P.``
180+
181+
Marked as ``tricubicFL``, this variation of tricubic interpolator is still under development by the author.
182+
179183
**Hermite**
180184

181185
The ``hermite``
@@ -185,13 +189,10 @@ The ``hermite``
185189
The weighted essentially non-oscillatory ``WENO`` used here is a fifth-order WENO reconstruction (WENO-5). It is suggested to be used in research with intermittent capture need, e.g., high-speed flows and shock capture.
186190
It shows relatively poor performance in general cases, and comsuming more wall time.
187191
The method originates from [Jiang1996]_ and expanded to three-dimensional computation, and [Shu2009]_ gave a review on the WEMO method.
188-
189192
The process is given as follows.
190193

191-
**Fifth-Order WENO Reconstruction (WENO-5)**
192-
193194
The WENO-5 method reconstructs a non-oscillatory, fifth-order-accurate approximation of a function value at an arbitrary location :math:`x = x_{i+1/2} + t\,\Delta x`, where :math:`t \in [0,1)` and :math:`x_{i+1/2} = x_i + \tfrac{1}{2}\,\Delta x` on a uniform grid with :math:`\Delta x = 1`.
194-
A five-point stencil ``{f_{i-2}``, ``f_{i-1}``, ``f_i``, ``f_{i+1}``, ``f_{i+2}}`` is used.
195+
A five-point stencil ``f_{i-2}``, ``f_{i-1}``, ``f_i``, ``f_{i+1}``, ``f_{i+2}`` is used.
195196

196197
.. math::
197198
@@ -266,7 +267,6 @@ Finally, reconstruct at :math:`x = x_{i+1/2} + t\,\Delta x` by combining:
266267
267268
f_{\mathrm{WENO5}}(x_{i+1/2} + t\,\Delta x) = \omega_{0}\,p_{0}(t) + \omega_{1}\,p_{1}(t) + \omega_{2}\,p_{2}(t).
268269
269-
The code computes :math:`t^{2}`, forms coefficients :math:`C_{\ell,k}`, evaluates :math:`p_{\ell}(t)`, computes :math:`\beta_{\ell}`, then :math:`\tilde{\alpha}_{\ell}`, :math:`\omega_{\ell}`, and returns :math:`\omega_{0}\,p_{0} + \omega_{1}\,p_{1} + \omega_{2}\,p_{2}` exactly as above.
270270
271271
272272
@@ -305,6 +305,11 @@ Eigenvalue Solver
305305
~~~~~~~~~~~~~~~~~~~~
306306

307307

308+
.. _numdyn:
309+
310+
Windowing for Dynamic LCS
311+
----------------------------------------------
312+
308313

309314

310315
.. _numcompare:

source/9_references.rst

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -12,10 +12,9 @@ Cited and Suggested by ``Py3DFTLE``
1212
1313
.. [Butcher1964] Butcher J. C. "On Runge-Kutta processes of high order." *Journal of the Australian Mathematical Society* **4(2)** (1964). `link <https://doi.org/10.1017/S1446788700023387>`__
1414
15+
.. [Jiang1996] G.-S. Jiang and C.-W. Shu. "Efficient Implementation of Weighted ENO Schemes." *Journal of Computational Physics* **126(1)** (1996). `link <https://doi.org/10.1006/jcph.1996.0130>`__
1516
16-
.. [Jiang1996] G.-S. Jiang and C.-W. Shu. "Efficient Implementation of Weighted ENO Schemes." *Journal of Computational Physics* **126(1)** (1996). `link <https://doi.org/10.1006/jcph.1996.0045>`__
17-
18-
.. [Shu2009] C.-W. Shu. "High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems." *SIAM Review* **51(1)** (2009). `link <https://doi.org/10.1137/080716034>`__
17+
.. [Shu2009] C.-W. Shu. "High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems." *SIAM Review* **51(1)** (2009). `link <https://doi.org/10.1137/070679065>`__
1918
2019
2120

0 commit comments

Comments
 (0)