 |
OGS
|
|
OGS
|
Global assembler for the monolithic scheme of the non-isothermal Richards flow.
Governing equations without vapor diffusion
The energy balance equation is given by
\[ (\rho c_p)^{eff}\dot T - \nabla (\mathbf{k}_T^{eff} \nabla T)+\rho^l c_p^l \nabla T \cdot \mathbf{v}^l = Q_T \]
with \(T\) the temperature, \((\rho c_p)^{eff}\) the effective volumetric heat capacity, \(\mathbf{k}_T^{eff} \) the effective thermal conductivity, \(\rho^l\) the density of liquid, \(c_p^l\) the specific heat capacity of liquid, \(\mathbf{v}^l\) the liquid velocity, and \(Q_T\) the point heat source. The effective volumetric heat can be considered as a composite of the contributions of solid phase and the liquid phase as
\[ (\rho c_p)^{eff} = (1-\phi) \rho^s c_p^s + S^l \phi \rho^l c_p^l \]
with \(\phi\) the porosity, \(S^l\) the liquid saturation, \(\rho^s \) the solid density, and \(c_p^s\) the specific heat capacity of solid. Similarly, the effective thermal conductivity is given by
\[ \mathbf{k}_T^{eff} = (1-\phi) \mathbf{k}_T^s + S^l \phi k_T^l \mathbf I \]
where \(\mathbf{k}_T^s\) is the thermal conductivity tensor of solid, \( k_T^l\) is the thermal conductivity of liquid, and \(\mathbf I\) is the identity tensor.
The mass balance equation is given by
\begin{eqnarray*} \left(S^l\beta - \phi\frac{\partial S}{\partial p_c}\right) \rho^l\dot p - S \left( \frac{\partial \rho^l}{\partial T} +\rho^l(\alpha_B -S) \alpha_T^s \right)\dot T\\ +\nabla (\rho^l \mathbf{v}^l) + S \alpha_B \rho^l \nabla \cdot \dot {\mathbf u}= Q_H \end{eqnarray*}
where \(p\) is the pore pressure, \(p_c\) is the capillary pressure, which is \(-p\) under the single phase assumption, \(\beta\) is a composite coefficient by the liquid compressibility and solid compressibility, \(\alpha_B\) is the Biot's constant, \(\alpha_T^s\) is the linear thermal expansivity of solid, \(Q_H\) is the point source or sink term, \(H(S-1)\) is the Heaviside function, and \( \mathbf u\) is the displacement. While this process does not contain a fully mechanical coupling, simplfied expressions can be given to approximate the latter term under certain stress conditions. The liquid velocity \(\mathbf{v}^l\) is described by the Darcy's law as
\[ \mathbf{v}^l=-\frac{{\mathbf k} k_{ref}}{\mu} (\nabla p - \rho^l \mathbf g) \]
with \({\mathbf k}\) the intrinsic permeability, \(k_{ref}\) the relative permeability, \(\mathbf g\) the gravitational force.
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] Function<name> density<type> [case] Linear<name> density<type> [case] WaterVapourDensity<name> diffusion<type> [case] VapourDiffusionFEBEX<name> poissons_ratio<type> [case] Constant<name> poissons_ratio<type> [case] Parameter<name> specific_heat_capacity<type> [case] Constant<name> specific_latent_heat<type> [case] LinearWaterVapourLatentHeat<name> storage<type> [case] Constant<name> storage_contribution<type> [case] Constant<name> thermal_conductivity<type> [case] Constant<name> thermal_diffusion_enhancement_factor <type> [case] Constant<name> thermal_expansivity<type> [case] Constant<name> thermal_expansivity_contribution<type> [case] Constant<name> thermal_osmosis_coefficient<type> [case] Constant<name> viscosity<type> [case] Constant<name> viscosity<type> [case] Curve<name> youngs_modulus<type> [case] Constant<name> youngs_modulus<type> [case] Parameter<name> biot_coefficient<type> [case] Constant<name> bishops_effective_stress<type> [case] BishopsPowerLaw<name> bishops_effective_stress<type> [case] BishopsSaturationCutoff<name> permeability<type> [case] Constant<name> permeability<type> [case] Parameter<name> porosity<type> [case] Constant<name> relative_permeability<type> [case] Constant<name> relative_permeability<type> [case] Curve<name> relative_permeability<type> [case] RelativePermeabilityVanGenuchten<name> saturation<type> [case] Constant<name> saturation<type> [case] Curve<name> saturation<type> [case] SaturationVanGenuchten<name> thermal_conductivity<type> [case] Curve<name> thermal_conductivity<type> [case] EffectiveThermalConductivityPorosityMixing<name> thermal_conductivity<type> [case] Parameter<name> thermal_longitudinal_dispersivity<type> [case] Constant<name> thermal_transversal_dispersivity<type> [case] Constant<name> tortuosity<type> [case] ConstantNo additional info.
1.11.0