• String (generic)

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.

• ABSTRACTCoupling_Base

Simulation problem to be solved.

• Bool (generic)

If true, the program will wait for key press before it terminates.

• String (generic)

Short description of the solved problem.
Is displayed in the main log, and possibly in other text output files.

• RECORDMesh

Computational mesh common to all equations.

• RECORDTimeGovernor

Simulation time frame and time step.

• ABSTRACTDarcyFlow

Flow equation, provides the velocity field as a result.

• ABSTRACTAdvectionProcess

Transport of soluted substances, depends on the velocity field from a Flow equation.

• ABSTRACTAdvectionProcess

Heat transfer, depends on the velocity field from a Flow equation.

• file nameinput

Input file with mesh description.

• Array [0, UINT]
• of Region

• RECORDPartition

Parameters of mesh partitioning algorithms.

• Bool (generic)

If true, print table of all used regions.

• SELECTIONTypes of search algorithm for finding intersection candidates.

Search algorithm for element intersections.

• Double [0, +inf)

Maximal snapping distance from Mesh in various search operations. In particular is used in ObservePoint to find closest mesh element and in FieldFormula to find closest surface element in plan view (Z projection).

• file nameoutput

Output file with neighboring data from mesh.

• String (generic)

Name (label) of the region. It has to be unique per single mesh.

• Integer [0, INT]

Id of the region to which you assign the name.

• Integer [0, INT]

Dimension of the region to which you assign the name.
The value is taken into account only if a new region is created.

• String (generic)

New name (label) of the region. It has to be unique per single mesh.

• String (generic)

The original region name in the input file, e.g. a physical volume name in the GMSH format.

• String (generic)

Name (label) of the region. It has to be unique per single mesh.

• Integer [0, INT]

Id of the region. If unset, a unique id will be generated automatically.

• Array [1, UINT]
• of INTEGER

List of ids of elements.

• String (generic)

Name (label) of the new region. It has to be unique per single mesh.

• Array [0, UINT]
• of INTEGER

List of region ids to be added to the new region set.

• Array [0, UINT]
• of STRING

List of region names (labels) to be added to the new region set.

• String (generic)

Name (label) of the new region. It has to be unique per single mesh.

• Array [2, 2]
• of STRING

List of exactly two region (set) names.
Supposing region sets r1, r2, the result includes all regions of r1 that are not in r2.

• String (generic)

Name (label) of the new region. It has to be unique per single mesh.

• Array [2, UINT]
• of STRING

List of two or more region (set) names.

• SELECTIONPartTool

Software package used for partitioning. See corresponding selection.

• SELECTIONGraphType

Algorithm for generating graph and its weights from a multidimensional mesh.

• TUPLETimeValue

Start time of the simulation.

• TUPLETimeValue

End time of the simulation.
The default value is higher than the age of the Universe (given in seconds).

• TUPLETimeValue

Initial guess for the time step.
It applies to equations that use an adaptive time stepping. If set to 0.0, the time step is determined in fully autonomous way, assuming the equation supports it.

• TUPLETimeValue

Soft lower limit for the time step.
Equation using an adaptive time stepping cannot suggest smaller time step. The actual time step can only decrease below the limit in order to match the prescribed input or output times.

• TUPLETimeValue

Hard upper limit for the time step.
The actual time step can only increase above the limit in order to match the prescribed input or output times.

• Array [0, UINT]
• of DtLimits

Allow to set a time dependent changes in min_dt and max_dt limits. This list is processed at individual times overwriting previous values of min_dt/max_dt. Limits equal to 0 are ignored and replaced with min_dt/max_dt values.

• Bool (generic)

Add all times defined in dt_limits table to the list of fixed TimeMarks.

• file nameoutput

Write used time steps to the given file in YAML format corresponding with the format of dt_limits.

• String (generic)

Common time unit of the equation.
This unit will be used for all time inputs and outputs within the equation. Individually, the common time unit can be overwritten for every declared time.
Time units are used in the following cases:
1) Time units of time value keys in: TimeGovernor, FieldDescriptors.
The common time unit can be overwritten for every declared time.
2) Time units in:
a) input fields: FieldFE and FieldTimeFunction
b) time steps definition of OutputTimeSet
Common time unit can be overwritten by one unit value for every whole mesh data file or time function.
3) Time units in output files: observation times, balance times, frame times of VTK and GMSH
Common time unit cannot be overwritten in these cases.

• Array [3, 3]
• of DOUBLE

Vector of the gravity force. Dimensionless.

• Array [0, UINT]
• of Flow_Darcy_MH_Data

Input data for Darcy flow model.

• RECORDNonlinearSolver

Non-linear solver for MH problem.

• RECORDOutputStream

Output stream settings.
Specify file format, precision etc.

• instance of RECORD
• generic typeEquationOutput
• parameter output_field_selection = Flow_Darcy_MH:OutputFields

Specification of output fields and output times.

• instance of RECORD
• generic typeOutput_DarcyMHSpecific
• parameter output_field_selection = Flow_Darcy_MH_specific:OutputFields

Output settings specific to Darcy flow model.
Includes raw output and some experimental functionality.

• RECORDBalance

Settings for computing mass balance.

• RECORDTimeGovernor

Time governor settings for the unsteady Darcy flow model.

• Integer [0, 2]

Number of Schur complements to perform when solving MH system.

• SELECTIONMH_MortarMethod

Method for coupling Darcy flow between dimensions on incompatible meshes. [Experimental]

• Array [1, UINT]
• of STRING

Labels of the regions where to set fields.

• Integer [0, INT]

ID of the region where to set fields.

• TUPLETimeValue

Apply field setting in this record after this time.
These times have to form an increasing sequence.

• instance of ABSTRACT
• generic typeField_R3_to_R[3,3]
• parameter element_input_type = DOUBLE

Anisotropy of the conductivity tensor. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Isotropic conductivity scalar. {$[ms^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Transition coefficient between dimensions. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Water source density. {$[s^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = Flow_Darcy_BC_Type

Boundary condition type. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Prescribed pressure value on the boundary. Used for all values of bc_type except none and seepage. See documentation of bc_type for exact meaning of bc_pressure in individual boundary condition types. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Incoming water boundary flux. Used for bc_types : total_flux, seepage, river. {$[ms^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Conductivity coefficient in the total_flux or the river boundary condition type. {$[s^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Critical switch pressure for seepage and river boundary conditions. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial condition for pressure in time dependent problems. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Storativity (in time dependent problems). {$[m^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Boundary piezometric head for BC types: dirichlet, robin, and river.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Boundary switch piezometric head for BC types: seepage, river.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial condition for the pressure given as the piezometric head.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• String (generic)

Python script given as in place string

• file nameinput

Python script given as external file

• String (generic)

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 ).

• String (generic)

Definition of unit.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• Array [1, UINT]
• of ARRAY

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• Array [1, UINT]
• of ARRAY

• String (generic)

The vector used to project evaluation point onto the surface.

• String (generic)

The name of region set considered as the surface. You have to set surface region if you want to use formula variable d.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• RECORDTableFunction

Values of time series initialization of Field.

• String (generic)

Definition of unit of all times defined in FieldTimeFunction.

• Array [2, UINT]
• of IndependentValue

Initizaliation values of Field.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• file nameinput

GMSH mesh with data. Can be different from actual computational mesh.

• SELECTIONFE_discretization

Section where to find the field.
Some sections are specific to file format: point_data/node_data, cell_data/element_data, -/element_node_data, native/-.
If not given by a user, we try to find the field in all sections, but we report an error if it is found in more than one section.

• String (generic)

The values of the Field are read from the $ElementData section with field name given by this key.

• Double (-inf, +inf)

Default value is set on elements which values have not been listed in the mesh data file.

• String (generic)

Definition of the unit of all times defined in the mesh data file.

• TUPLETimeValue

This key allows reading field data from the mesh data file shifted in time. Considering the time 't', field descriptor with time 'T', time shift 'S', then if 't > T', we read the time frame 't + S'.

• SELECTIONinterpolation

Type of interpolation applied to the input spatial data.
The default value 'equivalent_mesh' assumes the data being constant on elements living on the same mesh as the computational mesh, but possibly with different numbering. In the case of the same numbering, the user can set 'identical_mesh' to omit algorithm for guessing node and element renumbering. Alternatively, in case of different input mesh, several interpolation algorithms are available.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• String (generic)

Python script given as in place string

• file nameinput

Python script given as external file

• String (generic)

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 ).

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• String (generic)

• String (generic)

The vector used to project evaluation point onto the surface.

• String (generic)

The name of region set considered as the surface. You have to set surface region if you want to use formula variable d.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• RECORDTableFunction

Values of time series initialization of Field.

• String (generic)

Definition of unit of all times defined in FieldTimeFunction.

• Array [2, UINT]
• of IndependentValue

Initizaliation values of Field.

• RECORDUnit

Unit of the field values provided in the main input file, in the external file, or by a function (FieldPython).

• file nameinput

GMSH mesh with data. Can be different from actual computational mesh.

• SELECTIONFE_discretization

Section where to find the field.
Some sections are specific to file format: point_data/node_data, cell_data/element_data, -/element_node_data, native/-.
If not given by a user, we try to find the field in all sections, but we report an error if it is found in more than one section.

• String (generic)

The values of the Field are read from the $ElementData section with field name given by this key.

• Double (-inf, +inf)

Default value is set on elements which values have not been listed in the mesh data file.

• String (generic)

Definition of the unit of all times defined in the mesh data file.

• TUPLETimeValue

This key allows reading field data from the mesh data file shifted in time. Considering the time 't', field descriptor with time 'T', time shift 'S', then if 't > T', we read the time frame 't + S'.

• SELECTIONinterpolation

Type of interpolation applied to the input spatial data.
The default value 'equivalent_mesh' assumes the data being constant on elements living on the same mesh as the computational mesh, but possibly with different numbering. In the case of the same numbering, the user can set 'identical_mesh' to omit algorithm for guessing node and element renumbering. Alternatively, in case of different input mesh, several interpolation algorithms are available.

• ABSTRACTLinSys

Linear solver for MH problem.

• Double [0, +inf)

Residual tolerance.

• Integer [0, INT]

Minimum number of iterations (linear solutions) to use.
This is usefull if the convergence criteria does not characterize your goal well enough so it converges prematurely, possibly even without a single linear solution.If greater then 'max_it' the value is set to 'max_it'.

• Integer [0, INT]

Maximum number of iterations (linear solutions) of the non-linear solver.

• Bool (generic)

If a stagnation of the nonlinear solver is detected the solver stops. A divergence is reported by default, forcing the end of the simulation. By setting this flag to 'true', the solver ends with convergence success on stagnation, but it reports warning about it.

• Double [0, 1]

Residual tolerance relative to the initial error.

• Double [0, +inf)

Absolute residual tolerance.

• Integer [0, INT]

Maximum number of outer iterations of the linear solver.

• String (generic)

This options is passed to PETSC to create a particular KSP (Krylov space method).
If the string is left empty (by default), the internal default options is used.

• Double [0, 1]

Residual tolerance relative to the initial error.

• Integer [0, INT]

Maximum number of outer iterations of the linear solver.

• Integer [0, INT]

Maximum number of iterations of the linear solver with non-decreasing residual.

• Integer [0, INT]

Number of levels in the multilevel method (=2 for the standard BDDC).

• Bool (generic)

Use adaptive selection of constraints in BDDCML.

• Integer [0, 2]

• file nameoutput

File path to the connected output file.

• ABSTRACTOutputTime

File format of the output stream and possible parameters.

• Array [0, UINT]
• of TimeGrid

Output times used for fields that do not have their own output times defined.

• RECORDOutputMesh

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.

• Integer [0, INT]

The number of decimal digits used in output of floating point values.
Default is 17 decimal digits which are necessary to reproduce double values exactly after write-read cycle.

• Array [0, UINT]
• of ObservePoint

Array of observe points.

• SELECTIONVTK variant (ascii or binary)

Variant of output stream file format.

• Bool (generic)

Parallel or serial version of file format.

• Double [0, +inf)

The start time of the grid.

• Double [0, +inf)

The step of the grid. If not specified, the grid consists of the single time given by the begin key.

• Double [0, +inf)

The time greater or equal to the last time in the grid.

• String (generic)

Definition of unit of all defined times (begin, step and end time).

• Integer [1, 20]

Maximal level of refinement of the output mesh.

• Bool (generic)

Set true for using error_control_field. Set false for global uniform refinement to max_level.

• String (generic)

Name of an output field, according to which the output mesh will be refined. The field must be a SCALAR one.

• Double [0, +inf)

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.

• String (generic)

Optional point name, which has to be unique.
Any string that is a valid YAML key in record without any quoting can be used, however, using just alpha-numerical characters, and underscore instead of the space, is recommended.

• Array [3, 3]
• of DOUBLE

Initial point for the observe point search.

• Integer [0, 4]

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.

• String (generic)

The region of the initial element for snapping. Without snapping we make a projection to the initial element.

• Double [0, +inf)

Global value is define in Mesh by the key global_snap_radius.

• Array [0, UINT]
• of TimeGrid

Output times used for the output fields without is own time series specification.

• Bool (generic)

Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.

• Array [0, UINT]
• of FieldOutputSetting

Array of output fields and their individual output settings.

• Array [0, UINT]
• of PARAMETER

Array of the fields evaluated in the observe points of the associated output stream.

The field name (from selection).

• Array [0, UINT]
• of TimeGrid

Output times specific to particular field.

• SELECTIONDiscrete_output

Optional value. Implicit value is given by field and can be changed.

• Array [0, UINT]
• of TimeGrid

Output times used for the output fields without is own time series specification.

• Bool (generic)

Add all input time points of the equation, mentioned in the 'input_fields' list, also as the output points.

• Array [0, UINT]
• of FieldOutputSetting

Array of output fields and their individual output settings.

• Array [0, UINT]
• of PARAMETER

Array of the fields evaluated in the observe points of the associated output stream.

• Bool (generic)

SPECIAL PURPOSE. Computes error norms of the solution, particulary suited for non-compatible coupling models.

• file nameoutput

Output file with raw data from MH module.

• Array [0, UINT]
• of TimeGrid

• Bool (generic)

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.

• SELECTIONBalance_output_format

Format of output file.

• Bool (generic)

Compute cumulative balance over time. If true, then balance is calculated at each computational time step, which can slow down the program.

• file nameoutput

File name for output of balance.

• Array [3, 3]
• of DOUBLE

Vector of the gravity force. Dimensionless.

• Array [0, UINT]
• of RichardsLMH_Data

Input data for Darcy flow model.

• RECORDNonlinearSolver

Non-linear solver for MH problem.

• RECORDOutputStream

Output stream settings.
Specify file format, precision etc.

• instance of RECORD
• generic typeEquationOutput
• parameter output_field_selection = Flow_Darcy_MH:OutputFields

Specification of output fields and output times.

• instance of RECORD
• generic typeOutput_DarcyMHSpecific
• parameter output_field_selection = Flow_Darcy_MH_specific:OutputFields

Output settings specific to Darcy flow model.
Includes raw output and some experimental functionality.

• RECORDBalance

Settings for computing mass balance.

• RECORDTimeGovernor

Time governor settings for the unsteady Darcy flow model.

• Integer [0, 2]

Number of Schur complements to perform when solving MH system.

• SELECTIONMH_MortarMethod

Method for coupling Darcy flow between dimensions on incompatible meshes. [Experimental]

• RECORDSoilModel

Soil model settings.

• Array [1, UINT]
• of STRING

Labels of the regions where to set fields.

• Integer [0, INT]

ID of the region where to set fields.

• TUPLETimeValue

Apply field setting in this record after this time.
These times have to form an increasing sequence.

• instance of ABSTRACT
• generic typeField_R3_to_R[3,3]
• parameter element_input_type = DOUBLE

Anisotropy of the conductivity tensor. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Complement dimension parameter (cross section for 1D, thickness for 2D). {$[m^{3-d}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Isotropic conductivity scalar. {$[ms^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Transition coefficient between dimensions. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Water source density. {$[s^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = Flow_Darcy_BC_Type

Boundary condition type. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Prescribed pressure value on the boundary. Used for all values of bc_type except none and seepage. See documentation of bc_type for exact meaning of bc_pressure in individual boundary condition types. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Incoming water boundary flux. Used for bc_types : total_flux, seepage, river. {$[ms^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Conductivity coefficient in the total_flux or the river boundary condition type. {$[s^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Critical switch pressure for seepage and river boundary conditions. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial condition for pressure in time dependent problems. {$[m]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Storativity (in time dependent problems). {$[m^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Saturated water content {$ \theta_s $}.
Relative volume of water in a reference volume of a saturated porous media. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Residual water content {$ \theta_r $}.
Relative volume of water in a reference volume of an ideally dry porous media. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

The van Genuchten pressure head scaling parameter {$ \alpha $}.
It is 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}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

The van Genuchten exponent parameter {$ n $}. {$[-]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Boundary piezometric head for BC types: dirichlet, robin, and river.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Boundary switch piezometric head for BC types: seepage, river.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial condition for the pressure given as the piezometric head.

• SELECTIONSoil_Model_Type

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.

• Double [0, 1]

Fraction of the water content where we cut and rescale the curve.

• RECORDTimeGovernor

Time governor settings for the transport equation.

• RECORDBalance

Settings for computing mass balance.

• RECORDOutputStream

Output stream settings.
Specify file format, precision etc.

• Array [1, UINT]
• of Substance

Specification of transported substances.

• ABSTRACTSolute

Type of the numerical method for the transport equation.

• ABSTRACTReactionTerm

Reaction model involved in transport.

• String (generic)

Name of the substance.

• Double [0, +inf)

Molar mass of the substance [kg/mol].

• Array [0, UINT]
• of Solute_Advection_FV:Data

• instance of RECORD
• generic typeEquationOutput
• parameter output_field_selection = Solute_Advection_FV:OutputFields

Specification of output fields and output times.

• Array [1, UINT]
• of STRING

Labels of the regions where to set fields.

• Integer [0, INT]

ID of the region where to set fields.

• TUPLETimeValue

Apply field setting in this record after this time.
These times have to form an increasing sequence.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Porosity of the mobile phase. {$[-]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Density of concentration sources. {$[m^{-3}kgs^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Concentration flux. {$[s^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Concentration sources threshold. {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Boundary condition for concentration of substances. {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial values for concentration of substances. {$[m^{-3}kg]$}

• Double [0, +inf)

Density of the solvent [ {$kg.m^{-3}$} ].

• RECORDPetsc

Solver for the linear system.

• Array [0, UINT]
• of Solute_AdvectionDiffusion_DG:Data

Input fields of the equation.

• SELECTIONDG_variant

Variant of the interior penalty discontinuous Galerkin method.

• Integer [0, 3]

Polynomial order for the finite element in DG method (order 0 is suitable if there is no diffusion/dispersion).

• instance of RECORD
• generic typeEquationOutput
• parameter output_field_selection = Solute_AdvectionDiffusion_DG:OutputFields

Specification of output fields and output times.

• Array [1, UINT]
• of STRING

Labels of the regions where to set fields.

• Integer [0, INT]

ID of the region where to set fields.

• TUPLETimeValue

Apply field setting in this record after this time.
These times have to form an increasing sequence.

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Porosity of the mobile phase. {$[-]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Density of concentration sources. {$[m^{-3}kgs^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Concentration flux. {$[s^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Concentration sources threshold. {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = Solute_AdvectionDiffusion_BC_Type

Type of boundary condition. {$[-]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Dirichlet boundary condition (for each substance). {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Flux in Neumann boundary condition. {$[m^{1-d}kgs^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Conductivity coefficient in Robin boundary condition. {$[m^{4-d}s^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Initial values for concentration of substances. {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Longitudinal dispersivity in the liquid (for each substance). {$[m]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Transverse dispersivity in the liquid (for each substance). {$[m]$}

• Array [1, UINT]
• generic typeField_R3_to_R[3,3]
• parameter element_input_type = DOUBLE

Molecular diffusivity in the liquid (for each substance). {$[m^{2}s^{-1}]$}

• instance of ABSTRACT
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Rock matrix density. {$[m^{-3}kg]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Coefficient of linear sorption. {$[m^{3}kg^{-1}]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

Coefficient of diffusive transfer through fractures (for each substance). {$[-]$}

• Array [1, UINT]
• generic typeField_R3_to_R
• parameter element_input_type = DOUBLE

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. {$[-]$}