The physical parameters of the Braginskii model, evaluated at a
reference point. There is a small matlab script
tools/matlab/grillix_parameters.m, which computes these parameters
for a given machine.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
rhos |
NaN |
> 0.0 |
beginner | Drift scale at reference normalised to major radius. (see TODO) |
delta |
NaN |
> 0.0 |
beginner | Major radius normalised to drift scale at reference, used in front of ExB advection terms (see TODO). |
etapar0 |
NaN |
> 0.0 |
beginner | Normalised parallel resistivity at reference (see TODO). |
eta_i0 |
NaN |
> 0.0 |
beginner | Normalised ion resistivity at reference (see TODO). |
tratio |
NaN |
> 0.0 |
beginner | Ion to electron temperature ratio at reference (see TODO). |
beta |
NaN |
> 0.0 |
beginner | Normalized dynamical plasma beta at reference (see TODO). |
mass_ratio_ei |
NaN |
> 0.0 |
beginner | Electron to ion mass ratio (see TODO). |
nu_e0 |
NaN |
> 0.0 |
beginner | Normalised electron collision frequency at reference (see TODO). |
thermal_force_coeff |
0.71 |
> 0.0 |
expert | Thermal force coefficient (see TODO). |
boussinesq_on |
FALSE |
(TRUE, FALSE) |
expert | Switch for Boussinesq approximation. Rather do not swicth on, might be broken. |
Parameters related with parallel heat flux model and parallel viscosity in parallel momentum equation.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
heatflux_model |
'BRAGINSKII_LIM' |
('BRAGINSKII_LIM', 'LANDAU') |
medium | Expressions used for heat flux, either the Braginskii expression limited by free streaming (see TODO) or the Landau model (see TODO). |
chipar0_i |
NaN |
>= 0.0 |
beginner | Normalised ion heat conductivity from Braginskii expression (see TODO). |
chipar0_e |
NaN |
>= 0.0 |
beginner | Normalised ion heat conductivity from Braginskii expression (see TODO). |
heatflux_cutoff_i |
LARGE |
> 0.0 |
medium | Hard cutoff for ion heat flux (see TODO). Only for heatflux_model='BRAGINSKII_LIM'. |
heatflux_cutoff_e |
LARGE |
> 0.0 |
medium | Hard cutoff for electron heat flux (see TODO). Only for heatflux_model='BRAGINSKII_LIM'. |
free_streaming_fraction_e |
LARGE |
> 0.0 |
medium | Free-streaming fraction limiting parallel electron heat conductivity (see TODO). Only for heatflux_model='BRAGINSKII_LIM'. |
free_streaming_fraction_i |
LARGE |
> 0.0 |
medium | Free-streaming fraction limiting parallel ion heat conductivity (see TODO). Only for heatflux_model='BRAGINSKII_LIM'. |
free_streaming_fraction_qfac |
LARGE |
> 0.0 |
medium | Safety factor in the free-streaming heat conductivity limiter expression (see TODO) as well as in neocl. viscosity limiter (see TODO). Take e.g. q95. |
aspect_ratio_inv |
0. |
> 0.0 |
medium | Inverse aspect ratio used in expression for neoclassical viscosity limiter (see TODO). |
landau_numlorentzians_e |
7 |
(3, 7 ,12) |
expert | Number of Lorentzians used for computation of Landau heat flux for electrons (see TODO). Only for heatflux_model='LANDAU'. |
landau_numlorentzians_i |
7 |
(3, 7, 12) |
expert | Number of Lorentzians used for computation of Landau heat flux for ions (see TODO). Only for heatflux_model='LANDAU'. |
landau_flutter_lhs_on |
FALSE |
(TRUE, FALSE) |
expert | Switch if magnetic flutter terms are enabled in left hand side of Landau heat flux computation (see TODO). |
heatflux_timplicit_e |
TRUE |
(TRUE, FALSE) |
medium | Switch if Braginskii electron heat flux term is treated implicit (=TRUE) or explicit (=FALSE) in time. Only for heatflux_model='BRAGINSKII_LIM'. Must be FALSE for heatflux_model='LANDAU'. |
heatflux_timplicit_i |
TRUE |
(TRUE, FALSE) |
medium | Switch if Braginskii ion heat flux term is treated implicit (=TRUE) or explicit (=FALSE) in time. Only for heatflux_model='BRAGINSKII_LIM'. Must be FALSE for heatflux_model='LANDAU'. |
viscosity_timplicit_upar |
TRUE |
(TRUE, FALSE) |
medium | Switch parallel dissipation in parallel momentum equation is treated implicit (=TRUE) or explicit (=FALSE) in time. |
Parameters related with numerical dissipation.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
perpdiss_exclude_buffer |
FALSE |
(TRUE, FALSE) |
beginner | If true, gradually reduces dissipation into buffer region. |
perpdiss_order_ne |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on density (see TODO). |
perpdiss_coeff_ne |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on density (see TODO). |
perpdiss_order_vort |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on vorticity (see TODO). |
perpdiss_coeff_vort |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on vorticity (see TODO). |
perpdiss_order_te |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on electron temperature (see TODO). |
perpdiss_coeff_te |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on electron temperature (see TODO). |
perpdiss_order_pe |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on electron pressure (see TODO). |
perpdiss_coeff_pe |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on electron pressure (see TODO). |
perpdiss_order_ti |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on ion temperature (see TODO). |
perpdiss_coeff_ti |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on ion temperature (see TODO). |
perpdiss_order_pi |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on ion pressure (see TODO). |
perpdiss_coeff_pi |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on ion pressure (see TODO). |
perpdiss_order_upar |
1 |
> 0 |
beginner | Order of (hyper-)viscosity acting on parallel velocity (see TODO). |
perpdiss_coeff_upar |
0.0 |
>= 0.0 |
beginner | Perpendicular (hyper-)viscosity coefficent acting on parallel velocity (see TODO). |
perpdiss_order_ohm |
1 |
> 0 |
medium | Order of (hyper-)viscosity acting in Ohm's law (see TODO). |
perpdiss_coeff_ohm |
0.0 |
>= 0.0 |
medium | Perpendicular (hyper-)viscosity coefficent acting in Ohm's law (see TODO). |
pardiss_coeff_ne |
0.0 |
>= 0 |
beginner | Parallel viscosity coeffcient acting on density (see TODO). |
pardiss_coeff_vort |
0.0 |
>= 0 |
beginner | Parallel viscosity coeffcient acting on vorticity (see TODO). |
pardiss_coeff_ohm |
0.0 |
>= 0 |
medium | Parallel viscosity coeffcient acting in Ohm's law (see TODO). |
Parameters for buffer zone, where additional diffusion smoothes quantities.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
buffer_select |
NONE |
('NONE', 'ZERONULL', 'SINGLENULL', 'SAMPLE_MMS') |
beginner | Type of buffer zone depending on equilibrium. |
buffer_path |
params_braginskii |
String of maxlen=1024 | beginner | Filepath, where to find namelist for corresponding buffer_type parameters, i.e. for bselect='ZERONULL' &buffer_zeronull and for bselect='SINGLENULL' &buffer_singlenull. |
buffer_dbgout |
TRUE |
(TRUE, FALSE) |
beginner | Switch if buffer function written out as netcdf file. |
buffer_coeff_ne |
0.0 |
>= 0 |
beginner | Diffusion coefficient for density inside the buffer zone (see TODO). |
buffer_coeff_te |
0.0 |
>= 0 |
beginner | Diffusion coefficient for electron temperature inside the buffer zone (see TODO). |
buffer_coeff_pe |
0.0 |
>= 0 |
beginner | Diffusion coefficient for electron pressure inside the buffer zone (see TODO). |
buffer_coeff_ti |
0.0 |
>= 0 |
beginner | Diffusion coefficient for ion temperature inside the buffer zone (see TODO). |
buffer_coeff_pi |
0.0 |
>= 0 |
beginner | Diffusion coefficient for ion pressure inside the buffer zone (see TODO). |
buffer_coeff_vort |
0.0 |
>= 0 |
beginner | Diffusion coefficient for vorticity inside the buffer zone (see TODO). |
buffer_coeff_upar |
0.0 |
>= 0 |
beginner | Diffusion coefficient for parallel velocity inside the buffer zone (see TODO). |
buffer_coeff_ohm |
0.0 |
>= 0 |
beginner | Diffusion coefficient in Ohm's law inside the buffer zone (see TODO). |
Floors are present to avoid density and temperatures to become negative or zero.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
floor_ne |
EPS |
real | beginner | Minumum allowed density. |
floor_te |
EPS |
real | beginner | Minumum allowed electron temperature. |
floor_ti |
EPS |
real | beginner | Minumum allowed ion temperature. |
Switches for each individual term of the Braginskii model. This is an exhaustive list and for the individual parameters and default values the reader is referred to the corresponding source code file params_brag_switches_m.f90.
Parameters for selection of external source types.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
nesrc_select |
NONE |
('SOURCE_ZONAL', 'SOURCE_CONSTANTRATE') |
beginner | Type of density source (see TODO). |
nesrc_path |
'' |
String of maxlen=1024 | beginner | Filepath, where to find namelist for corresponding density source parameters, i.e. for nesrc_select='SOURCE_ZONAL' &source_zonal and for nesrc_select='SOURCE_CONSTANTRATE' source_constantrate(see TODO). |
tesrc_select |
NONE |
('SOURCE_ZONAL', 'SOURCE_CONSTANTRATE') |
beginner | Type of electron temperature source respectively electron thermal power source if modify_te_source_to_powersource=TRUE (see TODO). |
tesrc_path |
'' |
String of maxlen=1024 | beginner | Filepath, where to find namelist for corresponding electron thermal source parameters, i.e. for tesrc_select='SOURCE_ZONAL' &source_zonal and for tesrc_select='SOURCE_CONSTANTRATE' source_constantrate(see TODO). |
is_powersource |
FALSE |
(TRUE, FALSE) |
expert | Switch that modifies temperature source to power sources. |
tisrc_select |
NONE |
('SOURCE_ZONAL', 'SOURCE_CONSTANTRATE') |
beginner | Type of ion temperature source (see TODO). |
tisrc_path |
'' |
String of maxlen=1024 | beginner | Filepath, where to find namelist for corresponding ion thermal source parameters, i.e. for tisrc_select='SOURCE_ZONAL' &source_zonal and for tisrc_select='SOURCE_CONSTANTRATE' source_constantrate(see TODO). |
Parameters for perpendicular boundary conditions.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
bndtype_<ne,te,ti,jpar,upar>_<core,wall,dome,out> |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN') |
beginner | Type of boudbary condition at corresponding boundary section for corresponding quantity. |
bndval_<ne,te,ti,jpar,upar>_<core,wall,dome,out> |
'0.0' |
any real | beginner | Value (for BND_TYPE_DIRICHLET_ZERO) or gradient (for BND_TYPE_NEUMANN) at corresponding boundary section for corresponding quantity. |
bndtype_apar_<core,wall,dome,out> |
'BND_TYPE_DIRICHLET_ZERO' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN') |
beginner | Type of boudbary condition at corresponding boundary section for parallel electromagnetic potential. At least one boundary section must have Dirichlet boundary condition. |
bndval_apar_<core,wall,dome,out> |
'0.0' |
any real | beginner | Value (for BND_TYPE_DIRICHLET_ZERO) or gradient (for BND_TYPE_NEUMANN) at corresponding boundary section for parallel electromagnetic potential. |
bndtype_pot_core |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN', 'BND_BRAGTYPE_ZONAL_NEUMANN') |
medium | Type of boudbary condition at core boundary for electrostatic potential. |
bndval_pot_core |
'0.0' |
any real | medium | Value (for BND_TYPE_DIRICHLET_ZERO) or gradient (for BND_TYPE_NEUMANN, BND_BRAGTYPE_ZONAL_NEUMANN) at core boundary for electrostatic potential. |
bndtype_pot_<wall,dome,out> |
'BND_TYPE_DIRICHLET_ZERO' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN', 'BND_BRAGTYPE_FLOATING_POTENTIAL_LOCAL') |
medium | Type of boudbary condition at corresponding boundary section for electrostatic potential. BND_BRAGTYPE_FLOATING_POTENTIAL_LOCAL sets floating potential (pot=lambda*Te) irrespective of val_pot<wall,dome,out> (see TODO). |
bndval_pot_<wall,dome,out> |
'0.0' |
any real | medium | Value (for BND_TYPE_DIRICHLET_ZERO) or gradient (for BND_TYPE_NEUMANN) for electrostatic potential. |
bndtype_<qe,qi>_<core,wall,dome,out> |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN') |
expert | Type of boundary condition at corresponding boundary section for corresponding quantity. Currently only used for the boundary condition of the 3D solver within the Landau-fluid model. |
bndval_<qe,qi>_<core,wall,dome,out> |
'0.0' |
any real | expert | Value (for BND_TYPE_DIRICHLET_ZERO) or gradient (for BND_TYPE_NEUMANN) at corresponding boundary section for corresponding quantity. |
Parameters for parallel boundary conditions treated via penalisation.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
epsinv |
'GP_NAN' |
>= 0 |
beginner | Value for epsilon in penalisation source term (see TODO). |
bndtype_ne_pen |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN') |
beginner | Type of boundary for density in penalisation region (see TODO). |
val_ne_pen |
1.0 |
> 0 |
expert | Value of density in penalisation region, where magnetic field is directed into boundary. Only for bndtype_ne='BND_TYPE_DIRICHLET_ZERO'. |
bndtype_te_pen |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN','BND_BRAGTYPE_PARALLEL_HEAT_TRANSMISSION') |
beginner | Type of boundary for electron temperature in penalisation region (see TODO). |
val_te_pen |
1.0 |
> 0 |
expert | Value of electron temperature in penalisation region. Only for bndtype_te='BND_TYPE_DIRICHLET_ZERO'. |
sheath_heattransfac_e |
2.5 |
>= 0 |
medium | Electron heat sheath transmissiion factor. Only for bndtype_te='BND_BRAGTYPE_PARALLEL_HEAT_TRANSMISSION' (see TODO). |
bndtype_ti_pen |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN','BND_BRAGTYPE_PARALLEL_HEAT_TRANSMISSION') |
beginner | Type of boundary for ion temperature in penalisation region (see TODO). |
val_ti_pen |
1.0 |
> 0 |
expert | Value of ion temperature in penalisation region. Only for bndtype_ti='BND_TYPE_DIRICHLET_ZERO'. |
sheath_heattransfac_i |
0.1 |
>= 0 |
medium | Ion heat sheath transmissiion factor. Only for bndtype_ti='BND_BRAGTYPE_PARALLEL_HEAT_TRANSMISSION' (see TODO). |
bndtype_upar_pen |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN','BND_BRAGTYPE_SONIC_LOCAL', 'BND_BRAGTYPE_BOHM_LOCAL', 'BND_BRAGTYPE_SONIC_DRIFT_LOCAL', 'BND_BRAGTYPE_BOHM_DRIFT_LOCAL') |
beginner | Type of boundary for parallel velocity. Sonic refers to that parallel velocity is set equal to and Bohm to larger or equal than sound speed. Drift applies correction due to ExB and diamagnetic drift (see TODO). |
bohm_frac_ti |
1.0 |
> 0 |
medium | Fraction/switch of the ion temperature included in the calculation of the sound speed at the Bohm boundary conditions. Only for bndtype_upar='..._<SONIC,BOHM>_...' |
bndtype_vort_pen |
'BND_TYPE_NEUMANN' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN') |
beginner | Type of boundary for density in penalisation region (see TODO). |
val_vort_pen |
0.0 |
any real | expert | Value of vorticity in penalisation region. Only for bndtype_vort='BND_TYPE_DIRICHLET_ZERO'. |
bndtype_pot_pen |
'BND_TYPE_DIRICHLET_ZERO' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NONE' 'BND_BRAGTYPE_FLOATING_POTENTIAL_LOCAL') |
beginner | Type of boundary for electrostatic potential in penalisation region (see TODO). |
sheath_potential_lambda_sh |
2.69 |
>= 0 |
beginner | Value for potential drop across the sheath. Only for bndtype_pot='BND_BRAGTYPE_FLOATING_POTENTIAL_LOCAL' (see TODO). |
bndtype_ohm_pen |
'BND_BRAGTYPE_JPAR_SMOOTHPENZERO' |
('BND_TYPE_DIRICHLET_ZERO', 'BND_TYPE_NEUMANN' 'BND_BRAGTYPE_JPAR_SMOOTHPENZERO') |
medium | Type of boundary for modified parallel electromagnetci potential in penalisation region (see TODO). |
Parameters for selection of initial state.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
brag_iselect |
'NONE' |
('CONTINUE', 'RANDOM', 'BLOB', 'MMS', 'PROFILE') |
beginner | Type of initial state. |
brag_init_path |
'params_braginskii.in' |
String of maxlen=1024 | beginner | Filepath, where further parameters for specific initial state are read from. |
For heatflux_model='BRAGINSKII_LIM' and implicit_e=TRUE in namelist
&pardiss_model, this namelist contains parameters related with the
aligned 3D GMRES solver for the implicit treatment of the parallel heat
conduction term in the electron tempereature equation (see TODO). For
heatflux_model='LANDAU' in namelist &pardiss_model, it contains
parameters related with the aligned 3D GMRES solver for the parallel
electron heat conduction (see TODO).
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
solver |
PIM |
('PIM', ('CERFACS')) |
expert | Type of solver. 'CERFACS' only available if compiled with -DENABLE_CERFACS=ON |
krylov_method |
RGMRES |
('CG', 'RGMRES', 'BICGSTAB') |
medium | Krylov method to be used as solver |
resmax |
1.0E-8 |
> 0.0 |
medium | Requested tolerance of residuum of solution. |
maxiter |
10 |
> 0, < 50 |
medium | Maximum number of Krylov outer iterations. |
nrestart |
10 |
> 0, < 50 |
medium | Number of inner iterations before its restart (for 'RGMRES'). |
precond_type |
'JAC' |
('NONE', 'JAC') |
expert | Type of Preconditioner, either 'NONE' or Jacobi ('JAC') are available. |
dbgout |
0 |
(0, 1) |
beginner | Prints out debug information of solver (1) or not (0). |
For heatflux_model='BRAGINSKII_LIM' and implicit_i=TRUE in namelist
&pardiss_model, this namelist contains parameters related with the
aligned 3D GMRES solver for the implicit treatment of the parallel heat
conduction term in the ion tempereature equation (see TODO). For
heatflux_model='LANDAU' in namelist &pardiss_model, it contains
parameters related with the aligned 3D GMRES solver for the parallel
ion heat conduction (see TODO). Same as namelist &parsolver_te.
For implicit_upar=TRUE in namelist &pardiss_model, this namelist
contains parameters related with the aligned 3D GMRES solver for the
implicit treatment of the parallel dissipation term in the parallel
momentum equation (see TODO). Same as namelist &parsolver_te.
Parameters for selecting linout diagnostics data.
| Name | Default value | Valid values | Complexity level | Description |
|---|---|---|---|---|
lineout_select |
'NONE' |
('NONE', 'ON) |
medium | Switch for lineout diagnostcis. |
lineout_path |
'' |
String | medium | Filepath, where to find specialised input data for lineout diagnostics (see TODO). |