Root record of JSON input for Flow123d.
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.
If true, the program will wait for key press before it terminates.
The root record of description of particular the problem to solve.
Record with data for a general sequential coupling.\n
Short description of the solved problem.\nIs displayed in the main log, and possibly in other text output files.
Transport of soluted substances, depends on the velocity field from a Flow equation.
Heat transfer, depends on the velocity field from a Flow equation.
Record with mesh related data.
Input file with mesh description.
Parameters of mesh partitioning algorithms.\n
If true, print table of all used regions.
Search algorithm for element intersections.
Maximal distance of observe point from Mesh relative to its size (bounding box). Value is global and it can be rewrite at arbitrary ObservePoint by setting the key search_radius.
Abstract record for Region.
Elementary region declared by ID.Allows to create new region with given id and labelor specify existing region by id which will be renamed.
Label (name) of the region. Has to be unique in one mesh.\n
The ID of the region to which you assign label.
The dim of the region to which you assign label. Value is taken into account only if new region is created.
Gives a new name to an elementary regionwith original name (in the mesh file) given by mesh_label.
New label (name) of the region. Has to be unique in one mesh.
The mesh_label is e.g. physical volume name in GMSH format.
Elementary region declared by a list of elements. The new region is assigned to the list of elements spefied by the keyelement_list.
Label (name) of the region. Has to be unique in one mesh.\n
The ID of the region. If unset a unique ID will be generated automatically.
List of IDs of elements.
Defines region (set) as a union of given two or more regions. Regions can be given by names or IDs or both ways together.
Label (name) of the region. Has to be unique in one mesh.\n
List of region ID numbers that has to be added to the region set.
Defines region as a union of given pair of regions.
Defines region (set) as a difference of given pair of regions.
Label (name) of the region. Has to be unique in one mesh.
List of exactly two regions given by their names.
Defines region (set) as an intersection of given two or more regions.
Label (name) of the region. Has to be unique in one mesh.\n
List of two or more regions given by their names.
Setting for various types of mesh partitioning.
Algorithm for generating graph and its weights from a multidimensional mesh.
Select the partitioning tool to use.
Use PETSc interface to various partitioning tools.
Use direct interface to Metis.
Different algorithms to make the sparse graph with weighted edges\nfrom the multidimensional mesh. Main difference is dealing with \nneighborings of elements of different dimension.
Add edge for any pair of neighboring elements.
Same as before and assign higher weight to cuts of lower dimension in order to make them stick to one face.
Add edge for any pair of neighboring elements of same dimension (bad for matrix multiply).
Use BIH for finding initial candidates, then continue by prolongation.
Use BIH for finding all candidates.
Use bounding boxes for finding initial candidates, then continue by prolongation.
Setting of the simulation time. (can be specific to one equation)
Start time of the simulation.
End time of the simulation. Default value is more then age of universe in seconds.
Initial guess for the time step.\nOnly 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.
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.
Hard upper limit for the time step. Actual length of the time step is also limitedby input and output times.
Darcy flow model. Abstraction of various porous media flow models.
Mixed-Hybrid solver for STEADY saturated Darcy flow.
Vector of the gravity force. Dimensionless.
Parameters of output from MH module.
Time governor setting for the unsteady Darcy flow model.
Number of Schur complements to perform when solving MH system.
Method for coupling Darcy flow between dimensions.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any Flow_Darcy_MH_Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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' except for 'none' and 'seepage'. 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 : '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 for 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.
Abstract for all time-space functions.
R3_to_R[3,3] Field given by a Python script.
Python script given as in place string
Python script given as external file
Function in the given script that returns tuple containing components of the return type.\nFor NxM tensor values: tensor(row,col) = tuple( M*row + col ).
Specify unit of an input value. Evaluation of the unit formula results into a coeficient and a unit in terms of powers of base SI units. The unit must match expected SI unit of the value, while the value provided on the input is multiplied by the coefficient before further processing. The unit formula have a form:\n\n<UnitExpr>;<Variable>=<Number>*<UnitExpr>;...,\n\nwhere <Variable> is a variable name and <UnitExpr> is a units expression which consists of products and divisions of terms.\n\nA term has a form: <Base>^<N>, where <N> is an integer exponent and <Base> is either a base SI unit, a derived unit, or a variable defined in the same unit formula. Example, unit for the pressure head:\n\nMPa/rho/g_; rho = 990*kg*m^-3; g_ = 9.8*m*s^-2
Definition of unit.
R3_to_R[3,3] Field constant in space.
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\n\n - vector of size {$\\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)\n - scalar to enter multiple of the unit matrix.
R3_to_R[3,3] Field given by runtime interpreted formula.
String, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.\nFor vector values, you can use just one string to enter homogeneous vector.\nFor square {$N\\times N$}-matrix values, you can use:\n\n - array of strings of size {$N$} to enter diagonal matrix\n - array of strings of size {$\\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)\n - just one string to enter (spatially variable) multiple of the unit matrix.\nFormula can contain variables x,y,z,t and usual operators and functions.
R3_to_R[3,3] Field piecewise constant on mesh elements.
Input file with ASCII GMSH file format.
The values of the Field are read from the $ElementData section with field name given by this key.
R3_to_R[3,3] Field interpolated from external mesh data and piecewise constant on mesh elements.
Input file with ASCII GMSH file format.
The values of the Field are read from the $ElementData section with field name given by this key.
R3_to_R[3,3] Field time-dependent function in space.
Allow set variable series initialization of Fields.
Value of Field for independent variable.
Independent variable of stamp.
Value of the field in given stamp.
R3_to_R[3,3] Field given by finite element approximation.
GMSH mesh with data. Can be different from actual computational mesh.
Section where to find the field, some sections are specific to file format \npoint_data/node_data, cell_data/element_data, -/element_node_data, native/-.\nIf not given by user we try to find the field in all sections, but report error \nif it is found in more the one section.
The values of the Field are read from the $ElementData section with field name given by this key.
Specify the section in mesh input file where field data is listed.\nSome sections are specific to file format.
point_data (VTK) / node_data (GMSH)
cell_data (VTK) / element_data (GMSH)
element_node_data (only for GMSH)
native_data (only for VTK)
Abstract for all time-space functions.
R3_to_R Field given by a Python script.
Python script given as in place string
Python script given as external file
Function in the given script that returns tuple containing components of the return type.\nFor NxM tensor values: tensor(row,col) = tuple( M*row + col ).
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\n\n - vector of size {$\\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)\n - scalar to enter multiple of the unit matrix.
R3_to_R Field given by runtime interpreted formula.
String, array of strings, or matrix of strings with formulas for individual entries of scalar, vector, or tensor value respectively.\nFor vector values, you can use just one string to enter homogeneous vector.\nFor square {$N\\times N$}-matrix values, you can use:\n\n - array of strings of size {$N$} to enter diagonal matrix\n - array of strings of size {$\\frac12N(N+1)$} to enter symmetric matrix (upper triangle, row by row)\n - just one string to enter (spatially variable) multiple of the unit matrix.\nFormula can contain variables x,y,z,t and usual operators and functions.
R3_to_R Field piecewise constant on mesh elements.
Input file with ASCII GMSH file format.
The values of the Field are read from the $ElementData section with field name given by this key.
R3_to_R Field interpolated from external mesh data and piecewise constant on mesh elements.
Input file with ASCII GMSH file format.
The values of the Field are read from the $ElementData section with field name given by this key.
R3_to_R Field time-dependent function in space.
Allow set variable series initialization of Fields.
Value of Field for independent variable.
Independent variable of stamp.
Value of the field in given stamp.
R3_to_R Field given by finite element approximation.
GMSH mesh with data. Can be different from actual computational mesh.
Section where to find the field, some sections are specific to file format \npoint_data/node_data, cell_data/element_data, -/element_node_data, native/-.\nIf not given by user we try to find the field in all sections, but report error \nif it is found in more the one section.
The values of the Field are read from the $ElementData section with field name given by this key.
Homogeneous Neumann boundary condition. Zero flux
Dirichlet boundary condition. Specify the pressure head through the ''bc_pressure'' field or the piezometric head through the ''bc_piezo_head'' field.
Flux boundary condition (combines Neumann and Robin type). Water inflow equal to {$ \\delta_d(q_d^N + \\sigma_d (h_d^R - h_d) )$}. 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 face boundary condition. Pressure and inflow bounded from above. Boundary with potential seepage flow is described by the pair of inequalities: {$h_d \\le h_d^D$} and {$ -\\boldsymbol q_d\\cdot\\boldsymbol n \\le \\delta q_d^N$}, where the equality holds in at least one of them. Caution. Setting {$q_d^N$} strictly negative may 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_switch_pressure (or bc_switch_piezo_head) and bc_flux respectively.
River boundary condition. For the water level above the bedrock, {$H_d > H_d^S$}, the Robin boundary condition is used with the inflow given by: { $ \\delta_d(q_d^N + \\sigma_d(H_d^D - H_d) )$}. For the water level under the bedrock, constant infiltration is used: { $ \\delta_d(q_d^N + \\sigma_d(H_d^D - H_d^S) )$}. Parameters: bc_pressure, bc_switch_pressure, bc_sigma,bc_flux``.
Parameters to a non-linear solver.
Residual tolerance.
Minimum number of iterations (linear solves) to use. This is usefull if the convergence criteria does not characterize your goal well enough so it converges prematurely possibly without the single linear solve.If greater then 'max_it' the value is set to 'max_it'.
Maximum number of iterations (linear solves) of the non-linear solver.
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.
Linear solver setting.
Interface to PETSc solvers. Convergence criteria is:\n\nnorm( res_n ) < max( norm( res_0 ) * r_tol, a_tol )\n\nwhere 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.
Relative residual tolerance, (to initial error).
Absolute residual tolerance.
Maximum number of outer iterations of the linear solver.
Options passed to PETSC before creating KSP instead of default setting.
Solver setting.
Relative residual tolerance, (to initial error).
Maximum number of outer iterations of the linear solver.
Maximum number of iterations of the linear solver with non-decreasing residual.
Number of levels in the multilevel method (=2 for the standard BDDC).
Use adaptive selection of constraints in BDDCML.
Level of verbosity of the BDDCML library:\n\n - 0 - no output\n - 1 - mild output\n - 2 - detailed output.
Configuration of the spatial output of a single balance equation.
File path to the connected output file.
Output mesh record enables output on a refined mesh [EXPERIMENTAL, VTK only].Sofar refinement is performed only in discontinous sense.Therefore only corner and element data can be written on refined output mesh.Node data are to be transformed to corner data, native data cannot be written.Do not include any node or native data in output fields.
The number of decimal digits used in output of floating point values.\ Default is about 17 decimal digits which is enough to keep double values exect after write-read cycle.
Format of output stream and possible parameters.
Parameters of vtk output format.
Variant of output stream file format.
Parallel or serial version of file format.
ASCII variant of VTK file format
Uncompressed appended binary XML VTK format without usage of base64 encoding of appended data.
Appended binary XML VTK format without usage of base64 encoding of appended data. Compressed with ZLib.
Parameters of gmsh output format.
Equally spaced grid of time points.
The start time of the grid.
The step of the grid. If not specified, the grid consists only of the start time.
The time greater or equal to the last time in the grid.
Parameters of the refined output mesh.
Maximal level of refinement of the output mesh.
Set true for using error_control_field. Set false for global uniform refinement to max_level.
Name of an output field, according to which the output mesh will be refined. The field must be a SCALAR one.
Tolerance for element refinement by error. If tolerance is reached, refinement is stopped.Relative difference between error control field and its linear approximation on element is computedand compared with tolerance.
Specification of the observation point. The actual observe element and the observe point on it is determined as follows:\n\n1. Find an initial element containing the initial point. If no such element exists we report the error.\n2. Use BFS starting from the inital element to find the 'observe element'. The observe element is the closest element 3. Find the closest projection of the inital point on the observe element and snap this projection according to the 'snap_dim'.\n
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.
Initial point for the observe point search.
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.
The region of the initial element for snapping. Without snapping we make a projection to the initial element.
Global value is define in Mesh by the key global_observe_search_radius.
Output 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.
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.
Array of the fields evaluated in the observe points of the associated output stream.
Setting of the field output. The field name, output times, output interpolation (future).
The field name (from selection).
Optional value. Implicit value is given by field and can be changed.
Discrete type of output. Determines type of output data (element, node, native etc).
Node data / point data.
Corner data.
Element data / cell data.
Native data (Flow123D data).
Selection of output fields for the Flow_Darcy_MH model.\n
{$[m]$} Input field:
{$[m]$} Input field:
{$[m]$} Input field:
{$[ms^{-1}]$} Input field:
{$[-]$} Input field:
{$[-]$} Input field:
{$[-]$} Input field: Anisotropy of the conductivity tensor.
{$[m^{3-d}]$} Input field: Complement dimension parameter (cross section for 1D, thickness for 2D).
{$[ms^{-1}]$} Input field: Isotropic conductivity scalar.
{$[-]$} Input field: Transition coefficient between dimensions.
{$[s^{-1}]$} Input field: Water source density.
{$[m]$} Input field: Initial condition for pressure
{$[m^{-1}]$} Input field: Storativity.
{$[m]$} Input field:
{$[ms^{-1}]$} Input field:
{$[s^{-1}]$} Input field:
Specific Darcy flow MH output.
SPECIAL PURPOSE. Computing errors pro non-compatible coupling.
Output file with raw data form MH module.
Balance of a conservative quantity, boundary fluxes and sources.
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.
Compute cumulative balance over time. If true, then balance is calculated at each computational time step, which can slow down the program.
File name for output of balance.
Format of output file for balance.
Legacy format used by previous program versions.
Excel format with tab delimiter.
Format compatible with GnuPlot datafile with fixed column width.
Mortar space: P0 on elements of lower dimension.
Mortar space: P0 on elements of lower dimension.
Mortar space: P1 on intersections, using non-conforming pressures.
Lumped Mixed-Hybrid solver for unsteady saturated Darcy flow.
Vector of the gravity force. Dimensionless.
Parameters of output from MH module.
Time governor setting for the unsteady Darcy flow model.
Number of Schur complements to perform when solving MH system.
Method for coupling Darcy flow between dimensions.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any RichardsLMH_Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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' except for 'none' and 'seepage'. 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 : '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 for pressure {$[m]$}
Storativity. {$[m^{-1}]$}
\nSaturated water content {$ \\theta_s $}.\nrelative volume of the water in a reference volume of a saturated porous media.\n {$[-]$}
\nResidual water content {$ \\theta_r $}.\nRelative volume of the water in a reference volume of an ideally dry porous media.\n {$[-]$}
\nThe van Genuchten pressure head scaling parameter {$ \\alpha $}.\nThe parameter of the van Genuchten's model to scale the pressure head.\nRelated to the inverse of the air entry pressure, i.e. the pressure where the relative water content starts to decrease below 1.\n {$[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.
Setting for the 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.
Fraction of the water content where we cut and rescale the curve.
Van Genuchten soil model with cutting near zero.
Irmay model for conductivity, Van Genuchten model for the water content. Suitable for bentonite.
Abstract advection process. In particular: transport of substances or heat transfer.
Transport by convection and/or diffusion\ncoupled with reaction and adsorption model (ODE per element)\n via operator splitting.
Chemical substance.
Name of the substance.
Molar mass of the substance [kg/mol].
Transport of soluted substances.
Explicit in time finite volume method for advection only solute transport.
Setting of the fields output.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any Solute_Advection_FV:Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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]$}
Selection of output fields for the Solute_Advection_FV model.\n
{$[-]$} Input field: Mobile porosity
{$[-]$} Input field: INTERNAL - water content passed from unsaturated Darcy
{$[m^{-3}kgs^{-1}]$} Input field: Density of concentration sources.
{$[s^{-1}]$} Input field: Concentration flux.
{$[m^{-3}kg]$} Input field: Concentration sources threshold.
{$[m^{-3}kg]$} Input field: Initial concentrations.
{$[m^{-3}kg]$}
{$[-]$} Input field:
{$[-]$} Input field:
DG solver for solute transport.
Density of the solvent [ {$kg.m^(-3)$} ].
Input fields of the equation.
Variant of interior penalty discontinuous Galerkin method.
Polynomial order for finite element in DG method (order 0 is suitable if there is no diffusion/dispersion).
Setting of the field output.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any Solute_AdvectionDiffusion_DG:Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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 in the liquid (for each substance). {$[m]$}
Transversal dispersivity in the liquid (for each substance). {$[m]$}
Molecular diffusivity in the liquid (for each substance). {$[m^{2}s^{-1}]$}
Rock matrix density. {$[m^{-3}kg]$}
Coefficient of linear sorption. {$[m^{3}kg^{-1}]$}
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. {$[-]$}
Types of boundary conditions for advection-diffusion solute transport model.
Default transport boundary condition.\nOn 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 boundary condition {$ c = c_D $}.\nThe prescribed concentration {$c_D$} is specified by the field 'bc_conc'.
Total mass flux boundary condition.\nThe 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 boundary condition.\nThe 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'.
Type of penalty term.
non-symmetric weighted interior penalty DG method
incomplete weighted interior penalty DG method
symmetric weighted interior penalty DG method
Selection of output fields for the Solute_AdvectionDiffusion_DG model.\n
{$[-]$} Input field: Mobile porosity
{$[-]$} Input field: INTERNAL - water content passed from unsaturated Darcy
{$[m^{-3}kgs^{-1}]$} Input field: Density of concentration sources.
{$[s^{-1}]$} Input field: Concentration flux.
{$[m^{-3}kg]$} Input field: Concentration sources threshold.
{$[m^{-3}kg]$} Input field: Initial concentrations.
{$[m]$} Input field: Longitudal dispersivity in the liquid (for each substance).
{$[m]$} Input field: Transversal dispersivity in the liquid (for each substance).
{$[m^{2}s^{-1}]$} Input field: Molecular diffusivity in the liquid (for each substance).
{$[m^{-3}kg]$} Input field: Rock matrix density.
{$[m^{3}kg^{-1}]$} Input field: Coefficient of linear sorption.
{$[m^{-3}kg]$}
{$[-]$} Input field: Coefficient of diffusive transfer through fractures (for each substance).
{$[-]$} Input field: 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.
{$[-]$} Input field:
{$[-]$} Input field:
Abstract equation for a reaction term (dual porosity, sorption, reactions). Can be part of coupling with a transport equation via. operator splitting.
A model of first order chemical reactions (decompositions of a reactant into products).
Numerical solver for the system of first order ordinary differential equations coming from the model.
Describes a single first order chemical reaction.
An array of reactants. Do not use array, reactions with only one reactant (decays) are implemented at the moment!
The reaction rate coefficient of the first order reaction.
A record describing a reactant of a reaction.
The name of the reactant.
A record describing a product of a reaction.
The name of the product.
The branching ratio of the product when there are more products.\nThe value must be positive. Further, the branching ratios of all products are normalized in order to sum to one.\nThe default value 1.0, should only be used in the case of single product.
Record 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.
Polynomial degree of the nominator of Pade approximant.
Polynomial degree of the denominator of Pade approximant
A model of a radioactive decay and possibly of a decay chain.
Numerical solver for the system of first order ordinary differential equations coming from the model.
A model of a radioactive decay.
The name of the parent radionuclide.
The half life of the parent radionuclide in seconds.
An array of the decay products (daughters).
A record describing a product of a radioactive decay.
The name of the product.
Not used at the moment! The released energy in MeV from the decay of the radionuclide into the product.
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).
Sorption model in the reaction term of transport.
Names of the substances that take part in the sorption model.
Density of the solvent.
Number of equidistant substeps, molar mass and isotherm intersections
Specifies solubility limits of all the sorbing species.
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any Sorption:Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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. {$[m^{3}kg^{-1}]$}
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]$}
No sorption considered.
Linear isotherm runs the concentration exchange between liquid and solid.
Langmuir isotherm runs the concentration exchange between liquid and solid.
Freundlich isotherm runs the concentration exchange between liquid and solid.
Abstract equation for a reaction of species in single compartment (e.g. mobile solid).It can be part of: direct operator splitting coupling, dual porosity model, any sorption.
Selection of output fields for the Sorption model.\n
{$[m^{-3}kg]$} Input field: Rock matrix density.
{$[-]$} Input field: 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.
{$[m^{3}kg^{-1}]$} Input field: 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.
{$[-]$} Input field: Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).
{$[kg^{-1}mol]$} Input field: Initial solid concentration of substances. Vector, one value for every substance.
{$[-]$}
Dual porosity model in transport problems.\nProvides computing the concentration of substances in mobile and immobile zone.\n
Containes region specific data necessary to construct dual porosity model.
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.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any DualPorosity:Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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]$}
Abstract equation for a reaction term of the MOBILE pores (sorption, reactions). Is part of dual porosity model.
Sorption model in the mobile zone, following the dual porosity model.
Names of the substances that take part in the sorption model.
Density of the solvent.
Number of equidistant substeps, molar mass and isotherm intersections
Specifies solubility limits of all the sorbing species.
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
Selection of output fields for the SorptionMobile model.\n
{$[m^{-3}kg]$} Input field: Rock matrix density.
{$[-]$} Input field: 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.
{$[m^{3}kg^{-1}]$} Input field: 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.
{$[-]$} Input field: Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).
{$[kg^{-1}mol]$} Input field: Initial solid concentration of substances. Vector, one value for every substance.
{$[-]$}
Abstract equation for a reaction term of the IMMOBILE pores (sorption, reactions). Is part of dual porosity model.
Sorption model in the immobile zone, following the dual porosity model.
Names of the substances that take part in the sorption model.
Density of the solvent.
Number of equidistant substeps, molar mass and isotherm intersections
Specifies solubility limits of all the sorbing species.
Specifies highest aqueous concentration in interpolation table.
Containes region specific data necessary to construct isotherms.
Setting of the fields output.
Selection of output fields for the SorptionImmobile model.\n
{$[m^{-3}kg]$} Input field: Rock matrix density.
{$[-]$} Input field: 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.
{$[m^{3}kg^{-1}]$} Input field: 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.
{$[-]$} Input field: Second parameters (alpha, ...) defining isotherm c_s = omega * (alphac_a)/(1- alphac_a).
{$[kg^{-1}mol]$} Input field: Initial solid concentration of substances. Vector, one value for every substance.
{$[-]$}
Selection of output fields for the DualPorosity model.\n
{$[s^{-1}]$} Input field: Diffusion coefficient of non-equilibrium linear exchange between mobile and immobile zone.
{$[-]$} Input field: Porosity of the immobile zone.
{$[m^{-3}kg]$} Input field: Initial concentration of substances in the immobile zone.
{$[m^{-3}kg]$}
DG solver for heat transfer.
Input fields of the equation.
Variant of interior penalty discontinuous Galerkin method.
Polynomial order for finite element in DG method (order 0 is suitable if there is no diffusion/dispersion).
Setting of the field output.
Record to set fields of the equation.\nThe fields are set only on the domain specified by one of the keys: 'region', 'rid'\nand after the time given by the key 'time'. The field setting can be overridden by\n any Heat_AdvectionDiffusion_DG:Data record that comes later in the boundary data array.
Labels of the regions where to set fields.
ID of the region where to set fields.
Apply field setting in this record after this time.\nThese 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]$}
Density of thermal source in fluid. {$[m^{-1}kgs^{-3}]$}
Density of thermal source in solid. {$[m^{-1}kgs^{-3}]$}
Heat exchange rate of source 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. {$[-]$}
Types of boundary conditions for heat transfer model.
Default heat transfer boundary condition.\nOn 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 boundary condition {$T = T_D $}.\nThe prescribed temperature {$T_D$} is specified by the field 'bc_temperature'.
Total energy flux boundary condition.\nThe 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 boundary condition.\nThe 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'.
Selection of output fields for the Heat_AdvectionDiffusion_DG model.\n
{$[K]$} Input field: Initial temperature.
{$[-]$} Input field: Porosity.
{$[-]$} Input field:
{$[m^{-3}kg]$} Input field: Density of fluid.
{$[m^{2}s^{-2}K^{-1}]$} Input field: Heat capacity of fluid.
{$[mkgs^{-3}K^{-1}]$} Input field: Heat conductivity of fluid.
{$[m^{-3}kg]$} Input field: Density of solid (rock).
{$[m^{2}s^{-2}K^{-1}]$} Input field: Heat capacity of solid (rock).
{$[mkgs^{-3}K^{-1}]$} Input field: Heat conductivity of solid (rock).
{$[m]$} Input field: Longitudal heat dispersivity in fluid.
{$[m]$} Input field: Transversal heat dispersivity in fluid.
{$[m^{-1}kgs^{-3}]$} Input field: Density of thermal source in fluid.
{$[m^{-1}kgs^{-3}]$} Input field: Density of thermal source in solid.
{$[s^{-1}]$} Input field: Heat exchange rate of source in fluid.
{$[s^{-1}]$} Input field: Heat exchange rate of source in solid.
{$[K]$} Input field: Reference temperature of source in fluid.
{$[K]$} Input field: Reference temperature in solid.
{$[K]$}
{$[-]$} Input field: Coefficient of diffusive transfer through fractures (for each substance).
{$[-]$} Input field: 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.
{$[-]$} Input field:
{$[-]$} Input field: