Docs @ e362dc9

Author: MFC Action

documentation/md_case.html

@@ -238,7 +238,7 @@ <h2><a class="anchor" id="autotoc_md6"></a>
 <p>$$ x \in \left[ x \_ domain \% beg, x \_ domain \% end \right], y \in \left[ y \_ domain \% beg, y \_ domain \% end \right], z \in \left[ z \_ domain \% beg, z \_ domain \% end \right] $$</p>
 <ul>
 <li>$m$, $n$, and $p$ define the number of finite volume cells that uniformly discretize the domain along the $x$, $y$, and $z$ axes, respectively. Note that the actual number of cells in each coordinate axis is given as $[m,n,p]+1$. For example, $m=n=p=499$ discretizes the domain into $500^3$ cells. When the simulation is 2D/axi-symmetric or 1D, it requires that $p=0$ or $p=n=0$, respectively.</li>
-<li><code>stretch_[x,y,z]</code> activates grid stretching in the $[x,y,z]$ directions. The grid is gradually stretched such that the domain boundaries are pushed away from the origin along a specified axis. WENO7 does not support grid stretching.</li>
+<li><code>stretch_[x,y,z]</code> activates grid stretching in the $[x,y,z]$ directions. The grid is gradually stretched such that the domain boundaries are pushed away from the origin along a specified axis.</li>
 <li><code>a_[x,y,z]</code>, <code>[x,y,z]_a</code>, and <code>[x,y,z]_b</code> are parameters that define the grid stretching function. When grid stretching along the $x$ axis is considered, the stretching function is given as:</li>
 </ul>
 <p>$$ x_{cb,stretch} = x_{cb} + \frac{x_{cb}}{a_x} \Bigg[ \mathrm{log}\left[\mathrm{cosh} \left( \frac{a_x(x_{cb}-x_a)}{L} \right) \right] + \mathrm{log}\left[\mathrm{cosh} \left( \frac{a_x(x_{cb}-x_b)}{L} \right) \right] -2 \mathrm{log}\left[\mathrm{cosh} \left( \frac{a_x(x_b-x_a)}{2L} \right) \right] \Bigg] $$</p>
@@ -523,12 +523,12 @@ <h2><a class="anchor" id="autotoc_md14"></a>
 <li><code>mixture_err</code> activates correction of solutions to avoid imaginary speed of sound at each grid cell.</li>
 <li><code>time_stepper</code> specifies the order of the Runge-Kutta (RK) time integration scheme that is used for temporal integration in simulation, from the 1st to 5th order by corresponding integer. Note that <code>time_stepper = 3</code> specifies the total variation diminishing (TVD), third order RK scheme (<a href="references.md#Gottlieb98">Gottlieb and Shu, 1998</a>).</li>
 <li><code>adap_dt</code> activates the Strang operator splitting scheme which splits flux and source terms in time marching, and an adaptive time stepping strategy is implemented for the source term. It requires <code>bubbles = 'T'</code>, <code>polytropic = 'T'</code>, <code>adv_n = 'T'</code> and <code>time_stepper = 3</code>.</li>
-<li><code>weno_order</code> specifies the order of WENO scheme that is used for spatial reconstruction of variables by an integer of 1, 3, 5, and 7, that correspond to the 1st, 3rd, 5th, and 7th order, respectively. WENO7 does not support grid stretching.</li>
+<li><code>weno_order</code> specifies the order of WENO scheme that is used for spatial reconstruction of variables by an integer of 1, 3, 5, and 7, that correspond to the 1st, 3rd, 5th, and 7th order, respectively.</li>
 <li><code>weno_eps</code> specifies the lower bound of the WENO nonlinear weights. It is recommended to set <code>weno_eps</code> to $10^{-6}$ for WENO-JS, and to $10^{-40}$ for other WENO variants.</li>
 <li><code>mapped_weno</code> activates the WENO-M scheme in place of the default WENO-JS scheme (<a href="references.md#Henrick05">Henrick et al., 2005</a>). WENO-M a variant of the WENO scheme that remaps the nonlinear WENO-JS weights by assigning larger weights to non-smooth stencils, reducing dissipation compared to the default WENO-JS scheme, at the expense of higher computational cost. Only one of <code>mapped_weno</code>, <code>wenoz</code>, and <code>teno</code> can be activated.</li>
 <li><code>wenoz</code> activates the WENO-Z scheme in place of the default WENO-JS scheme (<a href="references.md#Borges08">Borges et al., 2008</a>). WENO-Z is a variant of the WENO scheme that further reduces the dissipation compared to the WENO-M scheme. It has similar computational cost to the WENO-JS scheme.</li>
 <li><code>wenoz_q</code> specifies the power parameter <code>q</code> used in the WENO-Z scheme. It controls how aggressively the smoothness coefficients scale the weights. A higher value of <code>wenoz_q</code> increases the sensitivity to smoothness, improving stability but worsening numerical dissipation. For WENO3 and WENO5, <code>q=1</code> is fixed, so <code>wenoz_q</code> must not be set. For WENO7, <code>wenoz_q</code> can be set to 2, 3, or 4.</li>
-<li><code>teno</code> activates the TENO scheme in place of the default WENO-JS scheme (<a href="references.md#Fu16">Fu et al., 2016</a>). TENO is a variant of the ENO scheme that is the least dissipative, but could be less robust for extreme cases. It uses a threshold to identify smooth and non-smooth stencils, and applies optimal weights to the smooth stencils. Only available for <code>weno_order = 5</code> and <code>7</code>. Requires <code>teno_CT</code> to be set.</li>
+<li><code>teno</code> activates the TENO scheme in place of the default WENO-JS scheme (<a href="references.md#Fu16">Fu et al., 2016</a>). TENO is a variant of the ENO scheme that is the least dissipative, but could be less robust for extreme cases. It uses a threshold to identify smooth and non-smooth stencils, and applies optimal weights to the smooth stencils. Only available for <code>weno_order = 5</code> and <code>7</code>. Requires <code>teno_CT</code> to be set. Does not support grid stretching.</li>
 <li><code>teno_CT</code> specifies the threshold for the TENO scheme. This dimensionless constant, also known as $C_T$, sets a threshold to identify smooth and non-smooth stencils. Larger values make the scheme more robust but also more dissipative. A recommended value for teno_CT is <code>1e-6</code>. When adjusting this parameter, it is recommended to try values like <code>1e-5</code> or <code>1e-7</code> for TENO5. A smaller value can be used for TENO7.</li>
 <li><code>null_weights</code> activates nullification of the nonlinear WENO weights at the buffer regions outside the domain boundaries when the Riemann extrapolation boundary condition is specified (<code>bc_[x,y,z]%beg[end]} = -4</code>).</li>
 <li><code>mp_weno</code> activates monotonicity preservation in the WENO reconstruction (MPWENO) such that the values of reconstructed variables do not reside outside the range spanned by WENO stencil (<a href="references.md#Balsara00">Balsara and Shu, 2000</a>; <a href="references.md#Suresh97">Suresh and Huynh, 1997</a>).</li>