Structure to handle the modelling of miscible or immiscible two-phase flows in a porous media. More...
#include <cs_gwf_hydraulic_model.h>
Structure to handle the modelling of miscible or immiscible two-phase flows in a porous media.
The model follows what is depicted in "Finite volume approximation of a diffusion-dissolution model and application to nuclear waste storage" O. Angelini, C. Chavant, E. Chénier, R. Eymard and S. Granet in Mathematics and Computers in Simulation (2011), 81 (10), pp. 2001–2017
Main assumptions are:
Notations are the following :
anderson_param |
Set of parameters to drive the Anderson acceleration (useful if the type of non-linear algorithm is set to the Anderson acceleration).
approx_type |
type of approximation used for the computation of diffusion, unsteady coefficients
b01_w_eqp |
Parameters associated to the (0,1) block in the system of equations. Water conservation w.r.t. the capillarity pressure.
b10_h_eqp |
Parameters associated to the (1,0) block in the system of equations. Conservation of the hydrogen w.r.t. the capillarity pressure.
c_pressure |
Pointer to a cs_field_t structure named "capillarity_pressure". Capillarity pressure
cell_weight |
weight associated to the cell value w.r.t. to the values at vertices when a CS_GWF_TPF_APPROX_PC_CELL_VERTEX_AVERAGE choice is set for 'approx_type'. If the value is < 0 or > 1, then the default settings is kept.
cs_property_t* diff_g_pty |
Property related to the relative permeability in the liquid phase
cs_property_t* diff_hc_pty |
Property related to the relative permeability in the liquid phase
cs_property_t* diff_hl_pty |
Property related to the relative permeability in the liquid phase
cs_property_t* diff_wl_pty |
Property related to the relative permeability in the liquid phase
enforce_pg_positivity |
Apply a rescaling on the gas pressure to avoid a negative value.
g_darcy |
Pointer to a cs_gwf_darcy_flux_t structure. Darcy advective flux in the gas phase
g_pressure |
Pointer to a cs_field_t structure named "gas_pressure". Pressure in the gas phase is denoted by .
g_rho_h |
Pointer to a cs_field_t structure. Property related to the mass density in the gas phase of the component mainly present in the gas phase (e.g. H2). The component present in the gas phase is ruled by the perfect gas law.
g_viscosity |
Viscosity in the gas phase (assumed to be constant) in Pa.s
h_eq |
Equation of mass conservation for (di)hydrogen. Hydrogen can be present in the liquid or in the gaseous phase. This corresponds to the M_10 and M_11 blocks in the system of equations along with the b_h right-hand side
h_molar_mass |
Molar mass of the main component in the gas phase (e.g. hydrogen) in kg.mol^-1
henry_constant |
Value of the Henry constant used in the Henry's law. Setting a very low value for this constant enables the model to degenerate into an immiscible model.
is_miscible |
true or false. Some simplifications can be done if immiscible.
krg_pty |
Property related to the relative permeability in the gas phase
krl_pty |
Property related to the relative permeability in the liquid phase
l_darcy |
Pointer to a cs_gwf_darcy_flux_t structure. Darcy advective flux in the liquid phase
l_diffusivity_h |
Molecular diffusivity of the hydrogen in the liquid phase in m^2.s^-1
l_mass_density |
Mass density in the liquid phase. With the model assumptions, this corresponds to the mass density of the main component in the liquid phase (e.g. water) in kg.m^-3
l_pressure |
Pointer to a cs_field_t structure named "liquid_pressure". Pressure in the liquid phase is denoted by .
l_rho_h |
Pointer to a cs_field_t structure. Property related to the mass density in the liquid phase of the component mainly present in the gas phase (e.g. H2). The component present in the liquid comes from the dissolution process (ruled by the Henry law)
l_saturation |
Pointer to a cs_field_t structure. Liquid saturation at cells. This quantity is denoted by and is defined by the soil model
l_viscosity |
Viscosity in the liquid phase (assumed to be constant) in Pa.s
lcap_pty |
Property related to the liquid capacity (derivative of the liquid saturation w.r.t. the capillarity pressure)
lsat_pre_array |
Values of the liquid saturation at the previous time step when an approximation CS_GWF_TPF_APPROX_VERTEX_SUBCELL has been chosen Size = c2v index size
lsat_pty |
Property related to the liquid saturation
nl_algo |
Structure used to manage the non-linearities
nl_algo_type |
Type of algorithm to solve the non-linearities
nl_algo_verbosity |
Level of verbosity for the algorithm solving the non-linearities
nl_cvg_param |
Set of parameters to drive the convergence of the non-linear solver
ref_temperature |
Reference temperature used in the "perfect gas" law (this is used when no thermal equation is solved). One expects a temperature in Kelvin.
solver_type |
Type of solver considered to solve the system of equations (choice of main unknowns and strategy of resolution (coupled/segregated))
srct_h_array |
Values of the source terms for the hydrogen conservation equation. Size = n_vertrices
srct_w_array |
Values of the source terms for the water conservation equation. Size = n_vertices
system |
System of equations (w_eq, h_eq and the cross-terms defined in the related cs_equation_param_t structures) used for the coupled approach
cs_property_t* time_hc_pty |
Property related to the relative permeability in the liquid phase
cs_property_t* time_hl_pty |
Property related to the relative permeability in the liquid phase
cs_property_t* time_wc_pty |
Property related to the relative permeability in the liquid phase
upwind_weight |
weight associated to the upwind value w.r.t. to the centered value There is no effect when CS_GWF_TPF_APPROX_VERTEX_SUBCELL is chosen to approximate coefficients. If the value is < 0 or > 1, then the default settings is kept.
w_eq |
Equation of mass conservation for water. Only the liquid phase is considered. One assumes no water vapour in the gaseous phase. This corresponds to the M_00 and M_01 blocks in the system of equations and to the b_w right-hand side.