version 4.2 authored by Udo Ziegler's avatar Udo Ziegler
...@@ -156,8 +156,7 @@ possible values or a numeric range. ...@@ -156,8 +156,7 @@ possible values or a numeric range.
physical evolution time `_C.time_max` is reached the simulation physical evolution time `_C.time_max` is reached the simulation
stops ahead of schedule. stops ahead of schedule.
- `03` (`_C.freq_log`, `_C.freq_nir`, `_C.freq_ana`, - `03` (`_C.freq_log`, `_C.freq_nir`, `_C.freq_ana`, `_C.freq_walltime`)
`_C.freq_walltime`)
- `_C.freq_log`: interval in units of timesteps at which the - `_C.freq_log`: interval in units of timesteps at which the
NIRVANA log file `nirvana.log` and monitoring file `nirvana.mon` NIRVANA log file `nirvana.log` and monitoring file `nirvana.mon`
...@@ -181,8 +180,10 @@ possible values or a numeric range. ...@@ -181,8 +180,10 @@ possible values or a numeric range.
- `01` (`_C.geometry`, `_C.omega[0-2]`) - `01` (`_C.geometry`, `_C.omega[0-2]`)
- `_C.geometry` ({CART,CYL,SPH}): choice of coordinate system - `_C.geometry` ({CART,CYL,SPH}): choice of coordinate system where
where CART=Cartesian, CYL=cylindrical and SPH=spherical. - CART: Cartesian
- CYL: cylindrical
- SPH: spherical
- `_C.omega[0-2]`: components of the angular velocity vector of a - `_C.omega[0-2]`: components of the angular velocity vector of a
rotating frame of reference with respect to the inertial frame rotating frame of reference with respect to the inertial frame
...@@ -202,32 +203,32 @@ possible values or a numeric range. ...@@ -202,32 +203,32 @@ possible values or a numeric range.
- `01` (`_C.lo[0]`, `_C.up[0]`, `_C.dim[0]`) - `01` (`_C.lo[0]`, `_C.up[0]`, `_C.dim[0]`)
- `_C.lo[0]`,\_C.up\[0\]: lower,upper $x$-coordinate of the - `_C.lo[0]`,\_C.up\[0\]: lower,upper **x**-coordinate of the
computational domain. computational domain.
- `_C.dim[0]`: number of *base-level* grid points in - `_C.dim[0]`: number of *base-level* grid points in
$x$-direction. `_C.dim[0]` must be an integral factor of 4, and **x**-direction. `_C.dim[0]` must be an integral factor of 4, and
excludes ghost cells which are automatically added by the code. excludes ghost cells which are automatically added by the code.
- `02` (`_C.lo[1]`, `_C.up[1]`, `_C.dim[1]`) - `02` (`_C.lo[1]`, `_C.up[1]`, `_C.dim[1]`)
- `_C.lo[1]`,`_C.up[1]`: lower,upper $y$-coordinate of the - `_C.lo[1]`,`_C.up[1]`: lower,upper **y**-coordinate of the
computational domain. In case of spherical geometry computational domain. In case of spherical geometry
($y\equiv \theta$) `_C.lo[1]`,`_C.up[1]` have to be specified in ($y\equiv \theta$) `_C.lo[1]`,`_C.up[1]` have to be specified in
units of $\pi$. units of $\pi$.
- `_C.dim[1]`: number of *base-level* grid points in - `_C.dim[1]`: number of *base-level* grid points in
$y$-direction. `_C.dim[1]` must be a multiple factor of 4. **y**-direction. `_C.dim[1]` must be a multiple factor of 4.
- `03` (`_C.lo[2]`, `_C.up[2]`, `_C.dim[2]`) - `03` (`_C.lo[2]`, `_C.up[2]`, `_C.dim[2]`)
- `_C.lo[2]`,`_C.up[2]`: lower,upper $z$-coordinate of the - `_C.lo[2]`,`_C.up[2]`: lower,upper **z**-coordinate of the
computational domain. In case of cylindrical- or spherical computational domain. In case of cylindrical- or spherical
geometry ($z\equiv \phi$) `_C.lo[2]`,`_C.up[2]` have to be geometry ($z\equiv \phi$) `_C.lo[2]`,`_C.up[2]` have to be
specified in units of $\pi$. specified in units of $\pi$.
- `_C.dim[2]`: number of *base-level* grid points in - `_C.dim[2]`: number of *base-level* grid points in
$z$-direction. `_C.dim[2]` must be a multiple factor of 4. If **z**-direction. `_C.dim[2]` must be a multiple factor of 4. If
`_C.dim[2]`=0 the simulation is assumed 2D, i.e., axisymmetric `_C.dim[2]`=0 the simulation is assumed 2D, i.e., axisymmetric
in case of cylindrical- or spherical coordinates. in case of cylindrical- or spherical coordinates.
...@@ -241,7 +242,7 @@ possible values or a numeric range. ...@@ -241,7 +242,7 @@ possible values or a numeric range.
SFC-decomposition is automatically used instead. SFC-decomposition is automatically used instead.
- `_C.bnx`,`_C.bny`,`_C.bnz`: number of domain subdivisions in - `_C.bnx`,`_C.bny`,`_C.bnz`: number of domain subdivisions in
$x,y,z$-direction in case \_C.partitioning_type=BLOCK. Numbers **x,y,z**-direction in case \_C.partitioning_type=BLOCK. Numbers
must be chosen such that the grid dimension of subdomains is a must be chosen such that the grid dimension of subdomains is a
multiple factor of 4 in each coordinate direction. Moreover, the multiple factor of 4 in each coordinate direction. Moreover, the
total number of subdomains must equal the number of MPI threads, total number of subdomains must equal the number of MPI threads,
...@@ -341,14 +342,14 @@ possible values or a numeric range. ...@@ -341,14 +342,14 @@ possible values or a numeric range.
- `01` (`_C.imr`, `_C.amr`) - `01` (`_C.imr`, `_C.amr`)
- `_C.imr` ($\le$`MAXLEVEL`): maximum refinement level for an - `_C.imr` (<=`MAXLEVEL`): maximum refinement level for an
initially refined mesh by the user (cf. [User-defined initial initially refined mesh by the user (cf. [User-defined initial
mesh refinement and refinement mesh refinement and refinement
control](#user-defined-initial-mesh-refinement-and-refinement-control)). control](#user-defined-initial-mesh-refinement-and-refinement-control)).
`_C.imr` cannot be larger than the macro `MAXLEVEL` defined in `_C.imr` cannot be larger than the macro `MAXLEVEL` defined in
the header file `nirvanaUser.h`. the header file `nirvanaUser.h`.
- `_C.amr` ($\le$`MAXLEVEL`): allowed maximum mesh refinement - `_C.amr` (<=`MAXLEVEL`): allowed maximum mesh refinement
level in a AMR simulations. `_C.amr` cannot be larger than the level in a AMR simulations. `_C.amr` cannot be larger than the
macro `MAXLEVEL` defined in the header file `nirvanaUser.h`. macro `MAXLEVEL` defined in the header file `nirvanaUser.h`.
...@@ -411,14 +412,13 @@ possible values or a numeric range. ...@@ -411,14 +412,13 @@ possible values or a numeric range.
grid cell. A zero or negative value means that the grid cell. A zero or negative value means that the
Jeans-length-based criterion is disabled. Jeans-length-based criterion is disabled.
- `_C.amr_dJeans` ($\ge 0$): tuning parameter for the - `_C.amr_dJeans` (>=0): tuning parameter for the
Jeans-length-based mesh refinement criterion allowing a Jeans-length-based mesh refinement criterion allowing a
systematic reduction of the Jeans threshold with increasing systematic reduction of the Jeans threshold with increasing
refinement level $l$ according to the expression refinement level l according to the expression
`_C.amr_Jeans` $-l*$ `_C.amr_dJeans` `_C.amr_Jeans`-l*`_C.amr_dJeans` i.e., the local
i.e., the local
Jeans length becomes gradually higher resolved with increasing Jeans length becomes gradually higher resolved with increasing
$l$. `_C.amr_dJeans` must be positiv. l. `_C.amr_dJeans` must be positiv.
**Important**: `_C.amr_dJeans` must be chosen with care such **Important**: `_C.amr_dJeans` must be chosen with care such
that the actual Jeans threshold never becomes too small or that the actual Jeans threshold never becomes too small or
...@@ -497,7 +497,7 @@ possible values or a numeric range. ...@@ -497,7 +497,7 @@ possible values or a numeric range.
numerical technique (equiv. to genuinely 2D-HLL in the numerical technique (equiv. to genuinely 2D-HLL in the
current implementation). current implementation).
- `_C.mhd_courant` (typical value: $<0.5$): CFL number in the - `_C.mhd_courant` (typical value: <0.5): CFL number in the
MHD timestep. MHD timestep.
- `02` (`_C.viscosity_solver`, `_C.viscosity_courant`) - `02` (`_C.viscosity_solver`, `_C.viscosity_courant`)
...@@ -512,7 +512,7 @@ possible values or a numeric range. ...@@ -512,7 +512,7 @@ possible values or a numeric range.
Coupling to MHD integrator is via Strang-type splitting. Coupling to MHD integrator is via Strang-type splitting.
- `_C.viscosity_courant`: CFL-like number in the fluid viscosity - `_C.viscosity_courant`: CFL-like number in the fluid viscosity
timestep. A typical value in case of STD is $<0.4$. Values much timestep. A typical value in case of STD is <0.4. Values much
larger than 1 are possible in the case of RKL. larger than 1 are possible in the case of RKL.
- `03` (`_C.diffusion_solver`, `_C.diffusion_courant`) - `03` (`_C.diffusion_solver`, `_C.diffusion_courant`)
...@@ -529,13 +529,13 @@ possible values or a numeric range. ...@@ -529,13 +529,13 @@ possible values or a numeric range.
- `06` (`_C.heatloss_max_change`) - `06` (`_C.heatloss_max_change`)
- `_C.heatloss_max_change` (typical value: $<0.1$): allowed - `_C.heatloss_max_change` (typical value: <0.1): allowed
maximal relative change in the temperature due to the heatloss maximal relative change in the temperature due to the heatloss
source term. source term.
- `07` (`_C.reactions_max_changeX`, `_C.reactions_max_changeT`) - `07` (`_C.reactions_max_changeX`, `_C.reactions_max_changeT`)
- `_C.reactions_max_changeX` (typical value: $<0.1$): allowed - `_C.reactions_max_changeX` (typical value: <0.1): allowed
maximal relative change of species number densities (or total maximal relative change of species number densities (or total
number density) in the time integration of the chemo-thermal number density) in the time integration of the chemo-thermal
rate equations. rate equations.
...@@ -546,7 +546,7 @@ possible values or a numeric range. ...@@ -546,7 +546,7 @@ possible values or a numeric range.
number density (`SXN`=1 uses individual number densities; number density (`SXN`=1 uses individual number densities;
`SXN`=0 uses the total number density). `SXN`=0 uses the total number density).
- `_C.reactions_max_changeT` (typical value: $<0.1$): allowed - `_C.reactions_max_changeT` (typical value: <0.1): allowed
maximal relative change in the temperature in the time maximal relative change in the temperature in the time
integration of the chemo-thermal rate equations. integration of the chemo-thermal rate equations.
...@@ -635,7 +635,7 @@ possible values or a numeric range. ...@@ -635,7 +635,7 @@ possible values or a numeric range.
perpendicular to the magnetic field (meaningless in isotropic perpendicular to the magnetic field (meaningless in isotropic
conduction). conduction).
- `_C.conduction_coeff_sat` (typical value: $0.3$): $\Psi$ is the - `_C.conduction_coeff_sat` (typical value: 0.3): $\Psi$ is the
parameter in the saturation heat flux model of \[[CM77](#references)\] (cf. parameter in the saturation heat flux model of \[[CM77](#references)\] (cf.
[physics [physics
guide](https://gitlab.aip.de/ziegler/NIRVANA/-/tree/master/doc/pdf/PhysicsGuide.pdf)). guide](https://gitlab.aip.de/ziegler/NIRVANA/-/tree/master/doc/pdf/PhysicsGuide.pdf)).
...@@ -675,7 +675,7 @@ possible values or a numeric range. ...@@ -675,7 +675,7 @@ possible values or a numeric range.
**Important:** Disabling the energy equation is not compatible **Important:** Disabling the energy equation is not compatible
with the use of an adiabatic EOS, for instance. with the use of an adiabatic EOS, for instance.
- `_C.energy_dual_sw` ($[0,1]$; typical value: $0.01$): threshold - `_C.energy_dual_sw` (\[0,1\]; typical value: <=0.01): threshold
value for the thermal-to-total energy density ratio in the dual value for the thermal-to-total energy density ratio in the dual
energy formalism (when `_C.energy`=DUAL). energy formalism (when `_C.energy`=DUAL).
... ...
......