OGS
|
A class to simulate the single phase flow process in porous media described by the governing equation:
\[ \left(\frac{\partial}{\partial p}(\phi \rho) + \rho \beta_s\right) \frac{\partial p}{\partial t} -\nabla \left(\rho\frac{\mathbf K}{\mu}(\nabla p + \rho g \nabla z)\right ) = \rho Q, \]
where
\begin{eqnarray*} &p:& \mbox{pore pressure,}\\ &\phi:& \mbox{porosity,}\\ &\rho:& \mbox{liquid or gas density,}\\ &\beta_s:& \mbox{specific storage,}\\ &{\mathbf K}:& \mbox{permeability,}\\ &\mu:& \mbox{viscosity,}\\ &g:& \mbox{gravitational constant,}\\ &Q:& \mbox{Source/sink term in m}^3/{\text s}.\\ \end{eqnarray*}
This governing equation represents the mass balance.
If the density is assumed constant, for example for a groundwater modelling, the governing equation is scaled with the density, and it becomes volume balanced as:
\[ \left(\frac{1}{\rho}\frac{\partial}{\partial p}(\phi \rho) + \beta_s\right) \frac{\partial p}{\partial t} -\nabla \left(\frac{\mathbf K}{\mu}(\nabla p + \rho g \nabla z)\right ) = Q, \]
An optional input tag equation_balance_type
of this process can be used to select whether to use the volume balanced equation or the mass balanced equation. By default, we assume that volume balanced equation is used.
Be aware that the Neumann condition is
\begin{eqnarray*} & -\frac{\mathbf K}{\mu}(\nabla p + \rho g \nabla z) \cdot \mathbf n = q_v [\text{m/s}]: & \mbox{ for the volume balance equation,}\\ & -\rho\frac{\mathbf K}{\mu}(\nabla p + \rho g \nabla z) \cdot \mathbf n = q_f [\text{kg/m²/s}]: & \mbox{for the mass balance equation,} \end{eqnarray*}
with \( \mathbf n \) the outer normal of the boundary.
Note: This list has been automatically extracted from OGS's benchmark tests (ctests). Therefore it might not be exhaustive, but it should give users a good overview about which properties they can/have to use with this process. Probably most of the properties occurring in this list are mandatory.
The list might contain different property <type>
s for some property <name>
to illustrate different possibilities the users have.
<type>
[case] AqueousLiquid<type>
[case] Gas<type>
[case] Solid<name>
density<type>
[case] Constant<name>
density<type>
[case] IdealGasLaw<name>
density<type>
[case] Linear<name>
molar_mass<type>
[case] Constant<name>
viscosity<type>
[case] Constant<name>
permeability<type>
[case] Constant<name>
permeability<type>
[case] Dupuit<name>
permeability<type>
[case] Parameter<name>
porosity<type>
[case] Constant<name>
reference_temperature<type>
[case] Constant<name>
storage<type>
[case] ConstantNo additional info.