Keystimes
no description provided
add_output_times
Add all output times of the balanced equation to the balance output times set. Note that this is not the time set of the output stream.
cumulative
Compute cumulative balance over time. If true, then balance is calculated at each computational time step, which can slow down the program.
file
File name for output of balance.
Implements abstract type: DescriptionSolver setting.
KeysTYPE
Sub-record Selection.
r_tol
Relative residual tolerance, (to initial error).
max_it
Maximum number of outer iterations of the linear solver.
max_nondecr_it
Maximum number of iterations of the linear solver with non-decreasing residual.
number_of_levels
Number of levels in the multilevel method (=2 for the standard BDDC).
bddcml_verbosity_level
Level of verbosity of the BDDCML library:
- 0 - no output
- 1 - mild output
- 2 - detailed output.
KeysTYPE
Sub-record Selection.
time
Time governor setting for the secondary equation.
balance
Settings for computing balance.
substances
Specification of transported substances.
transport
Type of numerical method for solute transport.
KeysTYPE
Sub-record Selection.
description
Short description of the solved problem.
Is displayed in the main log, and possibly in other text output files.
mesh
Computational mesh common to all equations.
time
Simulation time frame and time step.
flow_equation
Flow equation, provides the velocity field as a result.
solute_equation
Transport of soluted substances, depends on the velocity field from a Flow equation.
heat_equation
Heat transfer, depends on the velocity field from a Flow equation.
Keyshalf_life
The half life of the parent radionuclide in seconds.
products
An array of the decay products (daughters).
KeysTYPE
Sub-record Selection.
name
Label (name) of the region. Has to be unique in one mesh.
regions
Defines region as a difference of given pair of regions.
KeysTYPE
Sub-record Selection.
Containes region specific data necessary to construct dual porosity model.
scheme_tolerance
Tolerance according to which the explicit Euler scheme is used or not.Set 0.0 to use analytic formula only (can be slower).
Setting of the fields output.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any DualPorosity_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Diffusion coefficient of non-equilibrium linear exchange between mobile and immobile zone. {$[s^{-1}]$}
Porosity of the immobile zone. {$[-]$}
Initial concentration of substances in the immobile zone. {$[m^{-3}kg]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any DualPorosity_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Descriptionno description provided
DescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
Array of output fields and their individual output settings.
observe_fields
Array of the fields evaluated in the observe points of the associated output stream.
ParametersDescriptionOutput of the equation's fields.The output is done through the output stream of the associated balance law equation.The stream defines output format for the full space information in selected times and observe points for the full time information. The key 'fields' select the fields for the full spatial output.The set of output times may be specified per field otherwise common time set 'times' is used. If even this is not providedthe time set of the output_stream is used. The initial time of the equation is automatically added to the time set of every selected field. The end time of the equation is automatically added to the common output time set.
Keystimes
Output times used for the output fields without is own time series specification.
Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.
fields
Array of output fields and their individual output settings.
observe_fields
- Array [0, UINT]
- of PARAMETER
Array of the fields evaluated in the observe points of the associated output stream.
Generic parametersKeysTYPE
Sub-record Selection.
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
element_input_type
DOUBLE type of Double
KeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
Parameterselement_input_type
DOUBLE type of Double
KeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
ParametersKeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
ParametersKeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
ParametersKeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
ParametersKeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
KeysTYPE
Sub-record Selection.
value
- Array [1, UINT]
- of PARAMETER
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
KeysTYPE
Sub-record Selection.
value
- Array [1, 3]
- of PARAMETER
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
KeysTYPE
Sub-record Selection.
value
- Array [1, 2]
- of PARAMETER
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
KeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
KeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
Generic parametersKeysTYPE
Sub-record Selection.
value
Value of the constant field.
For vector values, you can use scalar value to enter constant vector.
For square {$N\times N$}-matrix values, you can use:
- vector of size {$N$} to enter diagonal matrix
- vector of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- scalar to enter multiple of the unit matrix.
Generic parametersKeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysString, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
KeysString, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
KeysString, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
Keys- Array [1, UINT]
- of STRING
String, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
Keys- Array [1, UINT]
- of STRING
String, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
KeysString, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.
For vector values, you can use just one string to enter homogeneous vector.
For square {$N\times N$}-matrix values, you can use:
- array of strings of size {$N$} to enter diagonal matrix
- array of strings of size {$\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)
- just one string to enter (spatially variable) multiple of the unit matrix.
Formula can contain variables x,y,z,t and usual operators and functions.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
KeysTYPE
Sub-record Selection.
gmsh_file
Input file with ASCII GMSH file format.
field_name
The values of the Field are read from the $ElementData section with field name given by this key.
Keysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
ParametersKeysfield
The field name (from selection).
times
Output times specific to particular field.
Generic parametersKeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
function
Function in the given script that returns tuple containing components of the return type.
For NxM tensor values: tensor(row,col) = tuple( M*row + col ).
KeysTYPE
Sub-record Selection.
reactions
An array of first order chemical reactions.
ode_solver
Numerical solver for the system of first order ordinary differential equations coming from the model.
Keysname
The name of the product.
branching_ratio
The branching ratio of the product when there are more products.
The value must be positive. Further, the branching ratios of all products are normalized in order to sum to one.
The default value 1.0, should only be used in the case of single product.
Keysname
The name of the reactant.
KeysTYPE
Sub-record Selection.
gravity
Vector of the gravitational acceleration (divided by the acceleration). Dimensionless, magnitude one for the Earth conditions.
Input data for Darcy flow model.
Parameters of output from MH module.
balance
Settings for computing mass balance.
time
Time governor setting for the unsteady Darcy flow model.
n_schurs
Number of Schur complements to perform when solving MH system.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Flow_Darcy_MH_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Anisotropy of the conductivity tensor. {$[-]$}
Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}
Isotropic conductivity scalar. {$[ms^{-1}]$}
Transition coefficient between dimensions. {$[-]$}
Water source density. {$[s^{-1}]$}
Boundary condition type, possible values: {$[-]$}
Prescribed pressure value on the boundary. Used for all values of 'bc_type' save the bc_type='none'.See documentation of 'bc_type' for exact meaning of 'bc_pressure' in individual boundary condition types. {$[m]$}
Incoming water boundary flux. Used for bc_types : 'none', 'total_flux', 'seepage', 'river'. {$[m^{4-d}s^{-1}]$}
Conductivity coefficient in the 'total_flux' or the 'river' boundary condition type. {$[m^{3-d}s^{-1}]$}
Critical switch pressure for 'seepage' and 'river' boundary conditions. {$[m]$}
Initial condition as pressure {$[m]$}
Storativity. {$[m^{-1}]$}
Boundary piezometric head for BC types: dirichlet, robin, and river.
Boundary switch piezometric head for BC types: seepage, river.
Initial condition for the pressure given as the piezometric head.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Flow_Darcy_MH_Data_aux_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Anisotropy of the conductivity tensor. {$[-]$}
Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}
Isotropic conductivity scalar. {$[ms^{-1}]$}
Transition coefficient between dimensions. {$[-]$}
Water source density. {$[s^{-1}]$}
Boundary condition type, possible values: {$[-]$}
Prescribed pressure value on the boundary. Used for all values of 'bc_type' save the bc_type='none'.See documentation of 'bc_type' for exact meaning of 'bc_pressure' in individual boundary condition types. {$[m]$}
Incoming water boundary flux. Used for bc_types : 'none', 'total_flux', 'seepage', 'river'. {$[m^{4-d}s^{-1}]$}
Conductivity coefficient in the 'total_flux' or the 'river' boundary condition type. {$[m^{3-d}s^{-1}]$}
Critical switch pressure for 'seepage' and 'river' boundary conditions. {$[m]$}
Initial condition as pressure {$[m]$}
Storativity. {$[m^{-1}]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Flow_Darcy_MH_Data_aux_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
KeysTYPE
Sub-record Selection.
gravity
Vector of the gravitational acceleration (divided by the acceleration). Dimensionless, magnitude one for the Earth conditions.
Input data for Darcy flow model.
Parameters of output from MH module.
balance
Settings for computing mass balance.
time
Time governor setting for the unsteady Darcy flow model.
n_schurs
Number of Schur complements to perform when solving MH system.
soil_model
Selection of the globally applied soil model. In future we replace this key by a field for selection of the model.That will allow usage of different soil model in a single simulation.
KeysTYPE
Sub-record Selection.
name
Label (name) of the region. Has to be unique in one mesh.
id
The ID of the region to which you assign label.
If new region is created and ID is not set, unique ID will be generated automatically.
element_list
- Array [1, UINT]
- of INTEGER
Specification of the region by the list of elements.
KeysTYPE
Sub-record Selection.
name
Label (name) of the region. Has to be unique in one mesh.
id
The ID of the region to which you assign label.
dim
The dim of the region to which you assign label. Value is taken into account only if new region is created.
KeysTYPE
Sub-record Selection.
name
New label (name) of the region. Has to be unique in one mesh.
mesh_label
The mesh_label is e.g. physical volume name in GMSH format.
KeysTYPE
Sub-record Selection.
time
Time governor setting for the secondary equation.
balance
Settings for computing balance.
solver
Linear solver for MH problem.
Input fields of the equation.
dg_variant
Variant of interior penalty discontinuous Galerkin method.
dg_order
Polynomial order for finite element in DG method (order 0 is suitable if there is no diffusion/dispersion).
Setting of the field output.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Heat_AdvectionDiffusion_DG_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Type of boundary condition. {$[-]$}
Boundary value of temperature. {$[K]$}
Flux in Neumann boundary condition. {$[m^{1-d}kgs^{-1}]$}
Conductivity coefficient in Robin boundary condition. {$[m^{4-d}s^{-1}]$}
Initial temperature. {$[K]$}
Porosity. {$[-]$}
Density of fluid. {$[m^{-3}kg]$}
Heat capacity of fluid. {$[m^{2}s^{-2}K^{-1}]$}
Heat conductivity of fluid. {$[mkgs^{-3}K^{-1}]$}
Density of solid (rock). {$[m^{-3}kg]$}
Heat capacity of solid (rock). {$[m^{2}s^{-2}K^{-1}]$}
Heat conductivity of solid (rock). {$[mkgs^{-3}K^{-1}]$}
Longitudal heat dispersivity in fluid. {$[m]$}
Transversal heat dispersivity in fluid. {$[m]$}
Thermal source density in fluid. {$[m^{-1}kgs^{-3}]$}
Thermal source density in solid. {$[m^{-1}kgs^{-3}]$}
Heat exchange rate in fluid. {$[s^{-1}]$}
Heat exchange rate of source in solid. {$[s^{-1}]$}
Reference temperature of source in fluid. {$[K]$}
Reference temperature in solid. {$[K]$}
Coefficient of diffusive transfer through fractures (for each substance). {$[-]$}
Penalty parameter influencing the discontinuity of the solution (for each substance). Its default value 1 is sufficient in most cases. Higher value diminishes the inter-element jumps. {$[-]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Heat_AdvectionDiffusion_DG_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
KeysTYPE
Sub-record Selection.
name
Label (name) of the region. Has to be unique in one mesh.
regions
- Array [2, UINT]
- of STRING
Defines region as an intersection of given pair of regions.
KeysTYPE
Sub-record Selection.
Keysregions
List of additional region and region set definitions not contained in the mesh.
There are three region sets implicitly defined:
- ALL (all regions of the mesh)
- .BOUNDARY (all boundary regions)
- and BULK (all bulk regions)
Keysmax_it
Maximal number of iterations (linear solves) of the non-linear solver.
converge_on_stagnation
If a stagnation of the nonlinear solver is detected the solver stops. A divergence is reported by default forcing the end of the simulation. Setting this flag to 'true', the solverends with convergence success on stagnation, but report warning about it.
Keysname
Optional point name. Has to be unique. Any string that is valid YAML key in record without any quoting can be used howeverusing just alpha-numerical characters and underscore instead of the space is recommended.
point
Initial point for the observe point search.
snap_dim
The dimension of the sub-element to which center we snap. For value 4 no snapping is done. For values 0 up to 3 the element containing the initial point is found and then the observepoint is snapped to the nearest center of the sub-element of the given dimension. E.g. for dimension 2 we snap to the nearest center of the face of the initial element.
snap_region
The region of the initial element for snapping. Without snapping we make a projection to the initial element.
n_search_levels
Maximum number of levels of the breadth first search used to find the observe element from the initial element. Value zero means to search only the initial element itself.
DescriptionParameters of output.
Keysfile
File path to the connected output file.
Format of output stream and possible parameters.
times
Output times used for equations without is own output times key.
output_mesh
Output mesh record enables output on a refined mesh.
precision
The number of decimal digits used in output of floating point values.
DescriptionParameters of output.
Keysmax_level
Maximal level of refinement of the output mesh.
refine_by_error
Set true for using error_control_field. Set false for global uniform refinement to max_level.
error_control_field
Name of an output field, according to which the output mesh will be refined. The field must be a SCALAR one.
Keyscompute_errors
SPECIAL PURPOSE. Computing errors pro non-compatible coupling.
DescriptionRecord with an information about pade approximant parameters.Note that stable method is guaranteed only if d-n=1 or d-n=2, where d=degree of denominator and n=degree of nominator. In those cases the Pade approximant corresponds to an implicit Runge-Kutta method which is both A- and L-stable. The default values n=2, d=3 yield relatively good precision while keeping the order moderately low.
KeysKeysSoftware package used for partitioning. See corresponding selection.
graph_type
Algorithm for generating graph and its weights from a multidimensional mesh.
Implements abstract type: DescriptionInterface to PETSc solvers. Convergence criteria is:
norm( res_n ) < max( norm( res_0 ) * r_tol, a_tol )
where res_i is the residuum vector after i-th iteration of the solver and res_0 is an estimate of the norm of initial residual.
If the initial guess of the solution is provided (usually only for transient equations) the residual of this estimate is used,
otherwise the norm of preconditioned RHS is used.
The default norm is L2 norm of preconditioned residual: {$ P^{-1}(Ax-b)$}, usage of other norm may be prescribed using the 'option' key.
See also PETSc documentation for KSPSetNormType.
KeysTYPE
Sub-record Selection.
r_tol
Relative residual tolerance, (to initial error).
a_tol
Absolute residual tolerance.
max_it
Maximum number of outer iterations of the linear solver.
options
Options passed to PETSC before creating KSP instead of default setting.
KeysTYPE
Sub-record Selection.
decays
An array of radioactive decays.
ode_solver
Numerical solver for the system of first order ordinary differential equations coming from the model.
Keysname
The name of the product.
energy
Not used at the moment! The released energy in MeV from the decay of the radionuclide into the product.
branching_ratio
The branching ratio of the product when there is more than one.Considering only one product, the default ratio 1.0 is used.Its value must be positive. Further, the branching ratios of all products are normalizedby their sum, so the sum then gives 1.0 (this also resolves possible rounding errors).
Keysreactants
An array of reactants. Do not use array, reactions with only one reactant (decays) are implemented at the moment!
reaction_rate
The reaction rate coefficient of the first order reaction.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any RichardsLMH_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Anisotropy of the conductivity tensor. {$[-]$}
Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}
Isotropic conductivity scalar. {$[ms^{-1}]$}
Transition coefficient between dimensions. {$[-]$}
Water source density. {$[s^{-1}]$}
Boundary condition type, possible values: {$[-]$}
Prescribed pressure value on the boundary. Used for all values of 'bc_type' save the bc_type='none'.See documentation of 'bc_type' for exact meaning of 'bc_pressure' in individual boundary condition types. {$[m]$}
Incoming water boundary flux. Used for bc_types : 'none', 'total_flux', 'seepage', 'river'. {$[m^{4-d}s^{-1}]$}
Conductivity coefficient in the 'total_flux' or the 'river' boundary condition type. {$[m^{3-d}s^{-1}]$}
Critical switch pressure for 'seepage' and 'river' boundary conditions. {$[m]$}
Initial condition as pressure {$[m]$}
Storativity. {$[m^{-1}]$}
Saturated water content {$ \theta_s $}.
relative volume of the water in a reference volume of a saturated porous media.
{$[-]$}
Residual water content {$ \theta_r $}.
Relative volume of the water in a reference volume of an ideally dry porous media.
{$[-]$}
The van Genuchten pressure head scaling parameter {$ \alpha $}.
The parameter of the van Genuchten's model to scale the pressure head.
Related to the inverse of the air entry pressure, i.e. the pressure where the relative water content starts to decrease below 1.
{$[m^{-1}]$}
The van Genuchten exponent parameter {$ n $}. {$[-]$}
Boundary piezometric head for BC types: dirichlet, robin, and river.
Boundary switch piezometric head for BC types: seepage, river.
Initial condition for the pressure given as the piezometric head.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any RichardsLMH_Data_aux_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Anisotropy of the conductivity tensor. {$[-]$}
Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}
Isotropic conductivity scalar. {$[ms^{-1}]$}
Transition coefficient between dimensions. {$[-]$}
Water source density. {$[s^{-1}]$}
Boundary condition type, possible values: {$[-]$}
Prescribed pressure value on the boundary. Used for all values of 'bc_type' save the bc_type='none'.See documentation of 'bc_type' for exact meaning of 'bc_pressure' in individual boundary condition types. {$[m]$}
Incoming water boundary flux. Used for bc_types : 'none', 'total_flux', 'seepage', 'river'. {$[m^{4-d}s^{-1}]$}
Conductivity coefficient in the 'total_flux' or the 'river' boundary condition type. {$[m^{3-d}s^{-1}]$}
Critical switch pressure for 'seepage' and 'river' boundary conditions. {$[m]$}
Initial condition as pressure {$[m]$}
Storativity. {$[m^{-1}]$}
Saturated water content {$ \theta_s $}.
relative volume of the water in a reference volume of a saturated porous media.
{$[-]$}
Residual water content {$ \theta_r $}.
Relative volume of the water in a reference volume of an ideally dry porous media.
{$[-]$}
The van Genuchten pressure head scaling parameter {$ \alpha $}.
The parameter of the van Genuchten's model to scale the pressure head.
Related to the inverse of the air entry pressure, i.e. the pressure where the relative water content starts to decrease below 1.
{$[m^{-1}]$}
The van Genuchten exponent parameter {$ n $}. {$[-]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any RichardsLMH_Data_aux_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Keysflow123d_version
Version of Flow123d for which the input file was created.Flow123d only warn about version incompatibility. However, external tools may use this information to provide conversion of the input file to the structure required by another version of Flow123d.
problem
Simulation problem to be solved.
pause_after_run
If true, the program will wait for key press before it terminates.
KeysTYPE
Sub-record Selection.
solver
Linear solver for MH problem.
Input fields of the equation.
dg_variant
Variant of interior penalty discontinuous Galerkin method.
dg_order
Polynomial order for finite element in DG method (order 0 is suitable if there is no diffusion/dispersion).
Setting of the field output.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Solute_AdvectionDiffusion_DG_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Mobile porosity {$[-]$}
Density of concentration sources. {$[m^{-3}kgs^{-1}]$}
Concentration flux. {$[s^{-1}]$}
Concentration sources threshold. {$[m^{-3}kg]$}
Type of boundary condition. {$[-]$}
Dirichlet boundary condition (for each substance). {$[m^{-3}kg]$}
Flux in Neumann boundary condition. {$[m^{1-d}kgs^{-1}]$}
Conductivity coefficient in Robin boundary condition. {$[m^{4-d}s^{-1}]$}
Initial concentrations. {$[m^{-3}kg]$}
Longitudal dispersivity (for each substance). {$[m]$}
Transversal dispersivity (for each substance). {$[m]$}
Molecular diffusivity (for each substance). {$[m^{2}s^{-1}]$}
Rock matrix density. {$[m^{-3}kg]$}
Coefficient of linear sorption. {$[kg^{-1}mol]$}
Coefficient of diffusive transfer through fractures (for each substance). {$[-]$}
Penalty parameter influencing the discontinuity of the solution (for each substance). Its default value 1 is sufficient in most cases. Higher value diminishes the inter-element jumps. {$[-]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Solute_AdvectionDiffusion_DG_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
KeysTYPE
Sub-record Selection.
Setting of the fields output.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Solute_Advection_FV_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Mobile porosity {$[-]$}
Density of concentration sources. {$[m^{-3}kgs^{-1}]$}
Concentration flux. {$[s^{-1}]$}
Concentration sources threshold. {$[m^{-3}kg]$}
Boundary conditions for concentrations. {$[m^{-3}kg]$}
Initial concentrations. {$[m^{-3}kg]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Solute_Advection_FV_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
KeysTYPE
Sub-record Selection.
substances
- Array [1, UINT]
- of STRING
Names of the substances that take part in the sorption model.
substeps
Number of equidistant substeps, molar mass and isotherm intersections
solubility
- Array [0, UINT]
- of DOUBLE
Specifies solubility limits of all the sorbing species.
table_limits
- Array [0, UINT]
- of DOUBLE
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
Keyssubstances
- Array [1, UINT]
- of STRING
Names of the substances that take part in the sorption model.
substeps
Number of equidistant substeps, molar mass and isotherm intersections
solubility
- Array [0, UINT]
- of DOUBLE
Specifies solubility limits of all the sorbing species.
table_limits
- Array [0, UINT]
- of DOUBLE
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
KeysTYPE
Sub-record Selection.
substances
- Array [1, UINT]
- of STRING
Names of the substances that take part in the sorption model.
substeps
Number of equidistant substeps, molar mass and isotherm intersections
solubility
- Array [0, UINT]
- of DOUBLE
Specifies solubility limits of all the sorbing species.
table_limits
- Array [0, UINT]
- of DOUBLE
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
KeysTYPE
Sub-record Selection.
substances
- Array [1, UINT]
- of STRING
Names of the substances that take part in the sorption model.
substeps
Number of equidistant substeps, molar mass and isotherm intersections
solubility
- Array [0, UINT]
- of DOUBLE
Specifies solubility limits of all the sorbing species.
table_limits
- Array [0, UINT]
- of DOUBLE
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Sorption_Data record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Rock matrix density. {$[m^{-3}kg]$}
Considered sorption is described by selected isotherm. If porosity on an element is equal or even higher than 1.0 (meaning no sorbing surface), then type 'none' will be selected automatically. {$[-]$}
Multiplication parameters (k, omega) in either Langmuir c_s = omega * (alphac_a)/(1- alphac_a) or in linear c_s = k * c_a isothermal description. {$[kg^{-1}mol]$}
Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a). {$[-]$}
Initial solid concentration of substances. Vector, one value for every substance. {$[kg^{-1}mol]$}
DescriptionRecord to set fields of the equation.
The fields are set only on the domain specified by one of the keys: 'region', 'rid'
and after the time given by the key 'time'. The field setting can be overridden by
any Sorption_Data_aux record that comes later in the boundary data array.
Keysregion
- Array [1, UINT]
- of STRING
Labels of the regions where to set fields.
rid
ID of the region where to set fields.
time
Apply field setting in this record after this time.
These times have to form an increasing sequence.
Constructible from key: name DescriptionChemical substance.
Keysname
Name of the substance.
Keysinit_dt
Initial guess for the time step.
Only useful for equations that use adaptive time stepping.If set to 0.0, the time step is determined in fully autonomous way if the equation supports it.
min_dt
Soft lower limit for the time step. Equation using adaptive time stepping can notsuggest smaller time step, but actual time step could be smaller in order to match prescribed input or output times.
max_dt
Hard upper limit for the time step. Actual length of the time step is also limitedby input and output times.
Keysbegin
The start time of the grid.
step
The step of the grid. If not specified, the grid consists only of the start time.
end
The time greater or equal to the last time in the grid.
KeysTYPE
Sub-record Selection.
name
Label (name) of the region. Has to be unique in one mesh.
region_ids
- Array [0, UINT]
- of INTEGER
List of region ID numbers that has to be added to the region set.
regions
- Array [0, UINT]
- of STRING
Defines region as a union of given pair of regions.
KeysTYPE
Sub-record Selection.
KeysTYPE
Sub-record Selection.
variant
Variant of output stream file format.
parallel
Parallel or serial version of file format.
compression
Compression used in output stream file format.
Implementations- DG solver for heat transfer.
- Transport by convection and/or diffusion
coupled with reaction and adsorption model (ODE per element)
via operator splitting.
Implementations- Record with data for a general sequential coupling.
Implementations- Mixed-Hybrid solver for STEADY saturated Darcy flow.
- Lumped Mixed-Hybrid solver for unsteady saturated Darcy flow.
Descriptionno description provided
Implementations- R2 -> R Field constant in space.
- R2 -> R Field given by a Python script.
- R2 -> R Field given by runtime interpreted formula.
- R2 -> R Field constant in space.
- R2 -> R Field constant in space.
Implementations- R2 -> R[2,2] Field constant in space.
- R2 -> R[2,2] Field given by a Python script.
- R2 -> R[2,2] Field given by runtime interpreted formula.
- R2 -> R[2,2] Field constant in space.
- R2 -> R[2,2] Field constant in space.
Implementations- R2 -> R[2] Field constant in space.
- R2 -> R[2] Field given by a Python script.
- R2 -> R[2] Field given by runtime interpreted formula.
- R2 -> R[2] Field constant in space.
- R2 -> R[2] Field constant in space.
Implementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
Parameterselement_input_type
DOUBLE type of Double
Implementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
ParametersImplementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
ParametersImplementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
ParametersImplementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
ParametersImplementations- R3 -> R Field constant in space.
- R3 -> R Field given by a Python script.
- R3 -> R Field given by runtime interpreted formula.
- R3 -> R Field constant in space.
- R3 -> R Field constant in space.
Generic parametersImplementations- R3 -> R[3,3] Field constant in space.
- R3 -> R[3,3] Field given by a Python script.
- R3 -> R[3,3] Field given by runtime interpreted formula.
- R3 -> R[3,3] Field constant in space.
- R3 -> R[3,3] Field constant in space.
Parameterselement_input_type
DOUBLE type of Double
Implementations- R3 -> R[3,3] Field constant in space.
- R3 -> R[3,3] Field given by a Python script.
- R3 -> R[3,3] Field given by runtime interpreted formula.
- R3 -> R[3,3] Field constant in space.
- R3 -> R[3,3] Field constant in space.
Generic parametersImplementations- R3 -> R[3] Field constant in space.
- R3 -> R[3] Field given by a Python script.
- R3 -> R[3] Field given by runtime interpreted formula.
- R3 -> R[3] Field constant in space.
- R3 -> R[3] Field constant in space.
Implementations- R3 -> R[n] Field constant in space.
- R3 -> R[n] Field given by a Python script.
DescriptionLinear solver setting.
Implementations- Interface to PETSc solvers. Convergence criteria is:
```norm( res_n ) < max( norm( res_0 ) * r_tol, a_tol )```
where res_i is the residuum vector after i-th iteration of the solver and res_0 is an estimate of the norm of initial residual.
If the initial guess of the solution is provided (usually only for transient equations) the residual of this estimate is used,
otherwise the norm of preconditioned RHS is used.
The default norm is L2 norm of preconditioned residual: (($ P^{-1}(Ax-b)$)), usage of other norm may be prescribed using the 'option' key.
See also PETSc documentation for KSPSetNormType.
Implementations- Evaluate analytic solution of the system of ODEs.
Implementations- Parameters of gmsh output format.
- Parameters of vtk output format.
Implementations- Dual porosity model in transport problems.
Provides computing the concentration of substances in mobile and immobile zone.
- Sorption model in the reaction term of transport.
- Sorption model in the mobile zone, following the dual porosity model.
- Sorption model in the immobile zone, following the dual porosity model.
- A model of a radioactive decay and possibly of a decay chain.
- A model of first order chemical reactions (decompositions of a reactant into products).
Implementations- Region declared by id and name.
Allows to create new region with given id and label
or specify existing region by id which will be renamed.
- Allows to rename existing region specified by mesh_label.
- Region declared by name, ID and enum of elements.
Allows to create new region and assign elements to its.
Elements are specified by ids.
- Defines region as a union of given two or more regions.
Regions can be given by names or IDs or both ways together.
- Defines region as a difference of given pair of regions.
- Defines region as an intersection of given two or more regions.
Implementations- DG solver for solute transport.
- Explicit in time finite volume method for advection only solute transport.
ValuesLegacy format used by previous program versions.
Excel format with tab delimiter.
Format compatible with GnuPlot datafile with fixed column width.
Valuesnon-symmetric
non-symmetric weighted interior penalty DG method
incomplete
incomplete weighted interior penalty DG method
symmetric
symmetric weighted interior penalty DG method
Valuesdiffusion_rate_immobile
Output of the field diffusion_rate_immobile {$[s^{-1}]$} (Diffusion coefficient of non-equilibrium linear exchange between mobile and immobile zone.).
porosity_immobile
Output of the field porosity_immobile {$[-]$} (Porosity of the immobile zone.).
init_conc_immobile
Output of the field init_conc_immobile {$[m^{-3}kg]$} (Initial concentration of substances in the immobile zone.).
conc_immobile
Output of the field conc_immobile {$[m^{-3}kg]$}.
Valuesnone
Homogeneous Neumann boundary condition. Zero flux
dirichlet
Dirichlet boundary condition. Specify the pressure head through the 'bc_pressure' field or the piezometric head through the 'bc_piezo_head' field.
total_flux
Flux boundary condition (combines Neumann and Robin type). Water inflow equal to {$q^N + \sigma (h^R - h)$}. Specify the water inflow by the 'bc_flux' field, the transition coefficient by 'bc_robin_sigma' and the reference pressure head or pieozmetric head through 'bc_pressure' or 'bc_piezo_head' respectively.
seepage
Seepage face boundary condition. Pressure and inflow bounded from above. Boundary with potential seepage flow is described by the pair of inequalities:{$h \le h_d^D$} and {$ q \le q_d^N$}, where the equality holds in at least one of them. Caution! Setting $q_d^N$ strictly negativemay lead to an ill posed problem since a positive outflow is enforced.Parameters {$h_d^D$} and {$q_d^N$} are given by fields bc_pressure (or bc_piezo_head) and bc_flux respectively.
river
River boundary condition. For the water level above the bedrock, {$H > H^S$}, the Robin boundary condition is used with the inflow given by: { $q^N + \sigma(H^D - H)$}. For the water level under the bedrock, constant infiltration is used { $q^N + \sigma(H^D - H^S)$}. Parameters: bc_pressure, bc_switch_pressure, bc_sigma,bc_flux``.
Valuesvan_genuchten
Van Genuchten soil model with cutting near zero.
irmay
Irmay model for conductivity, Van Genuchten model for the water content. Suitable for bentonite.
Valuespressure_p0
Output of the field pressure_p0 {$[m]$}.
pressure_p1
Output of the field pressure_p1 {$[m]$}.
piezo_head_p0
Output of the field piezo_head_p0 {$[m]$}.
velocity_p0
Output of the field velocity_p0 {$[ms^{-1}]$}.
subdomain
Output of the field subdomain {$[-]$}.
region_id
Output of the field region_id {$[-]$}.
anisotropy
Output of the field anisotropy {$[-]$} (Anisotropy of the conductivity tensor.).
cross_section
Output of the field cross_section {$[m^{3-d}]$} (Complement dimension parameter (cross section for 1D, thickness for 2D).).
conductivity
Output of the field conductivity {$[ms^{-1}]$} (Isotropic conductivity scalar.).
sigma
Output of the field sigma {$[-]$} (Transition coefficient between dimensions.).
water_source_density
Output of the field water_source_density {$[s^{-1}]$} (Water source density.).
init_pressure
Output of the field init_pressure {$[m]$} (Initial condition as pressure).
storativity
Output of the field storativity {$[m^{-1}]$} (Storativity.).
pressure_diff
Output of the field pressure_diff {$[m]$}.
velocity_diff
Output of the field velocity_diff {$[ms^{-1}]$}.
div_diff
Output of the field div_diff {$[s^{-1}]$}.
Valuesinit_temperature
Output of the field init_temperature {$[K]$} (Initial temperature.).
porosity
Output of the field porosity {$[-]$} (Porosity.).
water_content
Output of the field water_content {$[-]$}.
fluid_density
Output of the field fluid_density {$[m^{-3}kg]$} (Density of fluid.).
fluid_heat_capacity
Output of the field fluid_heat_capacity {$[m^{2}s^{-2}K^{-1}]$} (Heat capacity of fluid.).
fluid_heat_conductivity
Output of the field fluid_heat_conductivity {$[mkgs^{-3}K^{-1}]$} (Heat conductivity of fluid.).
solid_density
Output of the field solid_density {$[m^{-3}kg]$} (Density of solid (rock).).
solid_heat_capacity
Output of the field solid_heat_capacity {$[m^{2}s^{-2}K^{-1}]$} (Heat capacity of solid (rock).).
solid_heat_conductivity
Output of the field solid_heat_conductivity {$[mkgs^{-3}K^{-1}]$} (Heat conductivity of solid (rock).).
disp_l
Output of the field disp_l {$[m]$} (Longitudal heat dispersivity in fluid.).
disp_t
Output of the field disp_t {$[m]$} (Transversal heat dispersivity in fluid.).
fluid_thermal_source
Output of the field fluid_thermal_source {$[m^{-1}kgs^{-3}]$} (Thermal source density in fluid.).
solid_thermal_source
Output of the field solid_thermal_source {$[m^{-1}kgs^{-3}]$} (Thermal source density in solid.).
fluid_heat_exchange_rate
Output of the field fluid_heat_exchange_rate {$[s^{-1}]$} (Heat exchange rate in fluid.).
solid_heat_exchange_rate
Output of the field solid_heat_exchange_rate {$[s^{-1}]$} (Heat exchange rate of source in solid.).
fluid_ref_temperature
Output of the field fluid_ref_temperature {$[K]$} (Reference temperature of source in fluid.).
solid_ref_temperature
Output of the field solid_ref_temperature {$[K]$} (Reference temperature in solid.).
temperature
Output of the field temperature {$[K]$}.
fracture_sigma
Output of the field fracture_sigma {$[-]$} (Coefficient of diffusive transfer through fractures (for each substance).).
dg_penalty
Output of the field dg_penalty {$[-]$} (Penalty parameter influencing the discontinuity of the solution (for each substance). Its default value 1 is sufficient in most cases. Higher value diminishes the inter-element jumps.).
region_id
Output of the field region_id {$[-]$}.
Valuesinflow
Default heat transfer boundary condition.
On water inflow {$(q_w \le 0)$}, total energy flux is given by the reference temperature 'bc_temperature'. On water outflow we prescribe zero diffusive flux, i.e. the energy flows out only due to advection.
dirichlet
Dirichlet boundary condition {$T = T_D $}.
The prescribed temperature {$T_D$} is specified by the field 'bc_temperature'.
total_flux
Total energy flux boundary condition.
The prescribed incoming total flux can have the general form {$\delta(f_N+\sigma_R(T_R-T) )$}, where the absolute flux {$f_N$} is specified by the field 'bc_flux', the transition parameter {$\sigma_R$} by 'bc_robin_sigma', and the reference temperature {$T_R$} by 'bc_temperature'.
diffusive_flux
Diffusive flux boundary condition.
The prescribed incoming energy flux due to diffusion can have the general form {$\delta(f_N+\sigma_R(T_R-T) )$}, where the absolute flux {$f_N$} is specified by the field 'bc_flux', the transition parameter {$\sigma_R$} by 'bc_robin_sigma', and the reference temperature {$T_R$} by 'bc_temperature'.
ValuesNone
Mortar space: P0 on elements of lower dimension.
P0
Mortar space: P0 on elements of lower dimension.
P1
Mortar space: P1 on intersections, using non-conforming pressures.
Valuesinflow
Default transport boundary condition.
On water inflow {$(q_w \le 0)$}, total flux is given by the reference concentration 'bc_conc'. On water outflow we prescribe zero diffusive flux, i.e. the mass flows out only due to advection.
dirichlet
Dirichlet boundary condition {$ c = c_D $}.
The prescribed concentration {$c_D$} is specified by the field 'bc_conc'.
total_flux
Total mass flux boundary condition.
The prescribed total incoming flux can have the general form {$\delta(f_N+\sigma_R(c_R-c) )$}, where the absolute flux {$f_N$} is specified by the field 'bc_flux', the transition parameter {$\sigma_R$} by 'bc_robin_sigma', and the reference concentration {$c_R$} by 'bc_conc'.
diffusive_flux
Diffusive flux boundary condition.
The prescribed incoming mass flux due to diffusion can have the general form {$\delta(f_N+\sigma_R(c_R-c) )$}, where the absolute flux {$f_N$} is specified by the field 'bc_flux', the transition parameter {$\sigma_R$} by 'bc_robin_sigma', and the reference concentration {$c_R$} by 'bc_conc'.
Valuesporosity
Output of the field porosity {$[-]$} (Mobile porosity).
water_content
Output of the field water_content {$[-]$} (INTERNAL - water content passed from unsaturated Darcy).
sources_density
Output of the field sources_density {$[m^{-3}kgs^{-1}]$} (Density of concentration sources.).
sources_sigma
Output of the field sources_sigma {$[s^{-1}]$} (Concentration flux.).
sources_conc
Output of the field sources_conc {$[m^{-3}kg]$} (Concentration sources threshold.).
init_conc
Output of the field init_conc {$[m^{-3}kg]$} (Initial concentrations.).
disp_l
Output of the field disp_l {$[m]$} (Longitudal dispersivity (for each substance).).
disp_t
Output of the field disp_t {$[m]$} (Transversal dispersivity (for each substance).).
diff_m
Output of the field diff_m {$[m^{2}s^{-1}]$} (Molecular diffusivity (for each substance).).
rock_density
Output of the field rock_density {$[m^{-3}kg]$} (Rock matrix density.).
sorption_mult
Output of the field sorption_mult {$[kg^{-1}mol]$} (Coefficient of linear sorption.).
conc
Output of the field conc {$[m^{-3}kg]$}.
fracture_sigma
Output of the field fracture_sigma {$[-]$} (Coefficient of diffusive transfer through fractures (for each substance).).
dg_penalty
Output of the field dg_penalty {$[-]$} (Penalty parameter influencing the discontinuity of the solution (for each substance). Its default value 1 is sufficient in most cases. Higher value diminishes the inter-element jumps.).
region_id
Output of the field region_id {$[-]$}.
Valuesporosity
Output of the field porosity {$[-]$} (Mobile porosity).
water_content
Output of the field water_content {$[-]$} (INTERNAL - water content passed from unsaturated Darcy).
sources_density
Output of the field sources_density {$[m^{-3}kgs^{-1}]$} (Density of concentration sources.).
sources_sigma
Output of the field sources_sigma {$[s^{-1}]$} (Concentration flux.).
sources_conc
Output of the field sources_conc {$[m^{-3}kg]$} (Concentration sources threshold.).
init_conc
Output of the field init_conc {$[m^{-3}kg]$} (Initial concentrations.).
conc
Output of the field conc {$[m^{-3}kg]$}.
region_id
Output of the field region_id {$[-]$}.
Valuesrock_density
Output of the field rock_density {$[m^{-3}kg]$} (Rock matrix density.).
sorption_type
Output of the field sorption_type {$[-]$} (Considered sorption is described by selected isotherm. If porosity on an element is equal or even higher than 1.0 (meaning no sorbing surface), then type 'none' will be selected automatically.).
isotherm_mult
Output of the field isotherm_mult {$[kg^{-1}mol]$} (Multiplication parameters (k, omega) in either Langmuir c_s = omega * (alphac_a)/(1- alphac_a) or in linear c_s = k * c_a isothermal description.).
isotherm_other
Output of the field isotherm_other {$[-]$} (Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).).
init_conc_solid
Output of the field init_conc_solid {$[kg^{-1}mol]$} (Initial solid concentration of substances. Vector, one value for every substance.).
conc_immobile_solid
Output of the field conc_immobile_solid {$[m^{-3}kg]$}.
Valuesrock_density
Output of the field rock_density {$[m^{-3}kg]$} (Rock matrix density.).
sorption_type
Output of the field sorption_type {$[-]$} (Considered sorption is described by selected isotherm. If porosity on an element is equal or even higher than 1.0 (meaning no sorbing surface), then type 'none' will be selected automatically.).
isotherm_mult
Output of the field isotherm_mult {$[kg^{-1}mol]$} (Multiplication parameters (k, omega) in either Langmuir c_s = omega * (alphac_a)/(1- alphac_a) or in linear c_s = k * c_a isothermal description.).
isotherm_other
Output of the field isotherm_other {$[-]$} (Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).).
init_conc_solid
Output of the field init_conc_solid {$[kg^{-1}mol]$} (Initial solid concentration of substances. Vector, one value for every substance.).
conc_solid
Output of the field conc_solid {$[m^{-3}kg]$}.
Valuesnone
No sorption considered.
linear
Linear isotherm runs the concentration exchange between liquid and solid.
langmuir
Langmuir isotherm runs the concentration exchange between liquid and solid.
freundlich
Freundlich isotherm runs the concentration exchange between liquid and solid.
Valuesrock_density
Output of the field rock_density {$[m^{-3}kg]$} (Rock matrix density.).
sorption_type
Output of the field sorption_type {$[-]$} (Considered sorption is described by selected isotherm. If porosity on an element is equal or even higher than 1.0 (meaning no sorbing surface), then type 'none' will be selected automatically.).
isotherm_mult
Output of the field isotherm_mult {$[kg^{-1}mol]$} (Multiplication parameters (k, omega) in either Langmuir c_s = omega * (alphac_a)/(1- alphac_a) or in linear c_s = k * c_a isothermal description.).
isotherm_other
Output of the field isotherm_other {$[-]$} (Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).).
init_conc_solid
Output of the field init_conc_solid {$[kg^{-1}mol]$} (Initial solid concentration of substances. Vector, one value for every substance.).
conc_solid
Output of the field conc_solid {$[m^{-3}kg]$}.