You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
<tdclass="markdownTableBodyRight"><code>adap_dt</code></td><tdclass="markdownTableBodyCenter">Logical </td><tdclass="markdownTableBodyLeft">Strang splitting scheme with adaptive time stepping </td></tr>
485
485
<trclass="markdownTableRowOdd">
486
-
<tdclass="markdownTableBodyRight"><code>adap_dt_tol</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Tolerance for adaptive time stepping in Strang splitting scheme</td></tr>
<tdclass="markdownTableBodyRight"><code>adap_dt_max_iters</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">Max iteration for adaptive time stepping in Strang splitting scheme </td></tr>
488
+
<tdclass="markdownTableBodyRight"><code>adap_dt_tol</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">Tolerance for adaptive time stepping in Strang splitting scheme </td></tr>
489
489
<trclass="markdownTableRowOdd">
490
-
<tdclass="markdownTableBodyRight"><code>weno_order</code></td><tdclass="markdownTableBodyCenter">Integer </td><tdclass="markdownTableBodyLeft">WENO order [1,3,5]</td></tr>
490
+
<tdclass="markdownTableBodyRight"><code>adap_dt_max_iters</code></td><tdclass="markdownTableBodyCenter">Integer </td><tdclass="markdownTableBodyLeft">Max iteration for adaptive time stepping in Strang splitting scheme</td></tr>
491
491
<trclass="markdownTableRowEven">
492
-
<tdclass="markdownTableBodyRight"><code>weno_eps</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">WENO perturbation (avoid division by zero)</td></tr>
492
+
<tdclass="markdownTableBodyRight"><code>weno_order</code></td><tdclass="markdownTableBodyCenter">Integer</td><tdclass="markdownTableBodyLeft">WENO order [1,3,5]</td></tr>
493
493
<trclass="markdownTableRowOdd">
494
-
<tdclass="markdownTableBodyRight"><code>mapped_weno</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">WENO-M (WENO with mapping of nonlinear weights) </td></tr>
494
+
<tdclass="markdownTableBodyRight"><code>weno_eps</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">WENO perturbation (avoid division by zero) </td></tr>
<tdclass="markdownTableBodyRight"><code>mapped_weno</code></td><tdclass="markdownTableBodyCenter">Logical </td><tdclass="markdownTableBodyLeft">WENO-M (WENO with mapping of nonlinear weights)</td></tr>
497
497
<trclass="markdownTableRowOdd">
498
-
<tdclass="markdownTableBodyRight"><code>wenoz_q</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">WENO-Z power parameter q (only for WENO7)</td></tr>
<tdclass="markdownTableBodyRight"><code>wenoz_q</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">WENO-Z power parameter q (only for WENO7) </td></tr>
501
501
<trclass="markdownTableRowOdd">
502
-
<tdclass="markdownTableBodyRight"><code>teno_CT</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">TENO threshold for smoothness detection</td></tr>
<tdclass="markdownTableBodyRight"><code>null_weights</code></td><tdclass="markdownTableBodyCenter">Logical</td><tdclass="markdownTableBodyLeft">Null WENO weights at boundaries</td></tr>
504
+
<tdclass="markdownTableBodyRight"><code>teno_CT</code></td><tdclass="markdownTableBodyCenter">Real</td><tdclass="markdownTableBodyLeft">TENO threshold for smoothness detection</td></tr>
505
505
<trclass="markdownTableRowOdd">
506
+
<tdclass="markdownTableBodyRight"><code>null_weights</code></td><tdclass="markdownTableBodyCenter">Logical </td><tdclass="markdownTableBodyLeft">Null WENO weights at boundaries </td></tr>
<tdclass="markdownTableBodyRight"><code>muscl_order</code></td><tdclass="markdownTableBodyCenter">Integer </td><tdclass="markdownTableBodyLeft">MUSCL order [1,2] </td></tr>
511
+
<trclass="markdownTableRowEven">
512
+
<tdclass="markdownTableBodyRight"><code>muscl_lim</code></td><tdclass="markdownTableBodyCenter">Integer </td><tdclass="markdownTableBodyLeft">MUSCL Slope Limiter: [1] minmod; [2] monotonized central; [3] Van Albada; [4] Van Leer; [5] SUPERBEE </td></tr>
<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>
586
598
<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>
587
599
<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 (<aclass="el" href="md_references.html">Balsara and Shu, 2000</a>; <aclass="el" href="md_references.html">Suresh and Huynh, 1997</a>).</li>
588
-
<li><code>riemann_solver</code> specifies the choice of the Riemann solver that is used in simulation by an integer from 1 through 3. <code>riemann_solver = 1</code>, <code>2</code>, and <code>3</code> correspond to HLL, HLLC, and Exact Riemann solver, respectively (<aclass="el" href="md_references.html">Toro, 2013</a>). <code>riemann_solver = 4</code> is only for MHD simulations. It resolves 5 of the full seven-wave structure of the MHD equations (<aclass="el" href="md_references.html">Miyoshi and Kusano, 2005</a>).</li>
600
+
<li><code>muscl_order</code> specifies the order of the MUSCL scheme that is used for spatial reconstruction of variables by an integer of 1, or 2, that corresponds to the 1st, and 2nd order respectively. When using <code>muscl_order = 2</code>, <code>muscl_lim</code> must be defined.</li>
601
+
<li><code>muscl_lim</code> specifies the slope limiter that is used in 2nd order MUSCL Reconstruction by an integer from 1 through 5. <code>muscl_lim = 1</code>, <code>2</code>, <code>3</code>, <code>4</code>, and <code>5</code> correspond to minmod, monotonized central, Van Albada, Van Leer, and SUPERBEE, respectively.</li>
602
+
<li><code>int_comp</code> activates interface compression using THINC used in MUSCL Reconstruction, with control parameters (<code>ic_eps</code>, and <code>ic_beta</code>).</li>
603
+
<li><code>riemann_solver</code> specifies the choice of the Riemann solver that is used in simulation by an integer from 1 through 4. <code>riemann_solver = 1</code>, <code>2</code>, and <code>3</code> correspond to HLL, HLLC, and Exact Riemann solver, respectively (<aclass="el" href="md_references.html">Toro, 2013</a>). <code>riemann_solver = 4</code> is only for MHD simulations. It resolves 5 of the full seven-wave structure of the MHD equations (<aclass="el" href="md_references.html">Miyoshi and Kusano, 2005</a>).</li>
589
604
<li><code>low_Mach</code> specifies the choice of the low Mach number correction scheme for the HLLC Riemann solver. <code>low_Mach = 0</code> is default value and does not apply any correction scheme. <code>low_Mach = 1</code> and <code>2</code> apply the anti-dissipation pressure correction method (<aclass="el" href="md_references.html">Chen et al., 2022</a>) and the improved velocity reconstruction method (<aclass="el" href="md_references.html">Thornber et al., 2008</a>). This feature requires <code>model_eqns = 2</code> or <code>3</code>. <code>low_Mach = 1</code> works for <code>riemann_solver = 1</code> and <code>2</code>, but <code>low_Mach = 2</code> only works for <code>riemann_solver = 2</code>.</li>
590
605
<li><code>avg_state</code> specifies the choice of the method to compute averaged variables at the cell-boundaries from the left and the right states in the Riemann solver by an integer of 1 or 2. <code>avg_state = 1</code> and <code>2</code> correspond to Roe- and arithmetic averages, respectively.</li>
591
606
<li><code>wave_speeds</code> specifies the choice of the method to compute the left, right, and middle wave speeds in the Riemann solver by an integer of 1 and 2. <code>wave_speeds = 1</code> and <code>2</code> correspond to the direct method (<aclass="el" href="md_references.html">Batten et al., 1997</a>), and indirect method that approximates the pressures and velocity (<aclass="el" href="md_references.html">Toro, 2013</a>), respectively.</li>
<trclass="memdesc:ad82707f55f487d6868f76ea0074a8882"><tdclass="mdescLeft"> </td><tdclass="mdescRight">Checks constraints regarding WENO order. Called by s_check_inputs_common for all three stages. <br/></td></tr>
<trclass="memdesc:a799f031bb6612816625d841ba4e19e00"><tdclass="mdescLeft"> </td><tdclass="mdescRight">Check constraints regarding MUSCL order Called by s_check_inputs_common for all three stages. <br/></td></tr>
<trclass="memdesc:ae5e684a5330dc4903b6cb1234e0b4f79"><tdclass="mdescLeft"> </td><tdclass="mdescRight">Checks constraints on the boundary conditions in the x-direction. Called by s_check_inputs_common for all three stages. <br/></td></tr>
0 commit comments