Update 3.1 Code basics authored by Udo Ziegler's avatar Udo Ziegler
...@@ -428,15 +428,15 @@ criteria which are ...@@ -428,15 +428,15 @@ criteria which are
![amr_deriv](uploads/25c608e1a335aca4242b34afb9e4ffeb/amr_deriv.png) ![amr_deriv](uploads/25c608e1a335aca4242b34afb9e4ffeb/amr_deriv.png)
The criterion is applied to primary (or related) physical variables The criterion is applied to primary (or related) physical variables
*U* ( = {𝜚, **m**, *e*, **B**, *C*<sub>*c*</sub>}). Undivided first *U* ( = {𝜚, **m**, *e*, **B**, *C*<sub>*c*</sub>}). Undivided first
(* *U*) and second (*δ* <sup>2</sup>*U*) differences are (*δU*) and second (*δ*<sup>2</sup>*U*) differences are
computed along various spatial directions. The switch *α*\[0, 1\] computed along various spatial directions. The switch *α*\[0, 1\]
(code parameter: `_C.amr_d1`) selects between a purely (code parameter: `_C.amr_d1`) selects between a purely
gradient-based criterion (*α* = 1) and a purely gradient-based criterion (*α*=1) and a purely
second-derivatives-based criterion (*α* = 0). second-derivatives-based criterion (*α* = 0).
*U*<sub>*r**e**f*</sub> are reference values to avoid division by *U*<sub>*r**e**f*</sub> are reference values to avoid division by
zero for indefinite variables. The exponent *ξ*\[0, 2\] (code zero for indefinite variables. The exponent *ξ*\[0, 2\] (code
parameter: `_C.amr_exp`) introduces a power law functional on the parameter: `_C.amr_exp`) introduces a power law functional on the
grid spacing ratio * *x*<sup>(*l*)</sup>/* *x*<sup>(0)</sup> grid spacing ratio *δx*<sup>(*l*)</sup>/*δx*<sup>(0)</sup>
which allows to control the behavior with progressive mesh which allows to control the behavior with progressive mesh
refinement. The macro `FILTER` (preset to 10<sup> − 2</sup>) is a refinement. The macro `FILTER` (preset to 10<sup> − 2</sup>) is a
filter to suppress refinement at small-scale solution wiggles. The filter to suppress refinement at small-scale solution wiggles. The
...@@ -453,13 +453,11 @@ The criterion is applied to primary (or related) physical variables ...@@ -453,13 +453,11 @@ The criterion is applied to primary (or related) physical variables
flow features right in time.* flow features right in time.*
- Jeans-length-based (important in simulations involving selfgravity) - Jeans-length-based (important in simulations involving selfgravity)
![amr_Jeans](uploads/f4f4b19589de242e49db03a088acdf8f/amr_Jeans.png)
![amr_Jeans](uploads/88131f0e5135bb8e7b795ae51651b0bd/amr_Jeans.png)
![amr_Jeans](uploads/7cc51d2e1126f09b32c9187ea6f87490/amr_Jeans.png) ![amr_Jeans](uploads/7cc51d2e1126f09b32c9187ea6f87490/amr_Jeans.png)
where the first factor is the local Jeans length with where the first factor is the local Jeans length with
*c*<sub>*s*</sub> the sound speed and *G* the gravitational *c*<sub>*s*</sub> the sound speed and *G* the gravitational
constant, and where constant, and where
* *s* = min {* *x*, *h*<sub>*y*</sub>* *y*, *h*<sub>*z*</sub>* *z*}. *δs* = min {*δx*, *h*<sub>*y*</sub>*δy*, *h*<sub>*z*</sub>*δz*}.
<sub>*J**e**a**n**s*</sub> is a user-specific threshold with the <sub>*J**e**a**n**s*</sub> is a user-specific threshold with the
reciprocal value giving the minimum number of grid cells the local reciprocal value giving the minimum number of grid cells the local
Jeans length should be resolved. Jeans length should be resolved.
... ...
......