8.3
general documentation
cs_gwf_tpf_t Struct Reference

Structure to handle the modelling of miscible or immiscible two-phase flows in a porous media. More...

#include <cs_gwf_hydraulic_model.h>

+ Collaboration diagram for cs_gwf_tpf_t:

Data Fields

Equations and system of equations
cs_equation_tw_eq
 
cs_equation_th_eq
 
cs_equation_param_tb01_w_eqp
 
cs_equation_param_tb10_h_eqp
 
cs_equation_system_tsystem
 
Darcy (advection) fields
cs_gwf_darcy_flux_tl_darcy
 
cs_gwf_darcy_flux_tg_darcy
 
Properties related to the model
cs_property_tkrl_pty
 
cs_property_tkrg_pty
 
cs_property_tlsat_pty
 
cs_property_tlcap_pty
 
cs_property_ttime_wc_pty
 
cs_property_tdiff_wl_pty
 
cs_property_ttime_hc_pty
 
cs_property_tdiff_hc_pty
 
cs_property_ttime_hl_pty
 
cs_property_tdiff_hl_pty
 
cs_property_tdiff_g_pty
 
Additional fields
cs_field_tc_pressure
 
cs_field_tl_pressure
 
cs_field_tg_pressure
 
cs_field_tl_saturation
 
cs_field_tl_rho_h
 
cs_field_tg_rho_h
 
Additional arrays
cs_real_tsrct_w_array
 
cs_real_tsrct_h_array
 
cs_real_tlsat_pre_array
 
Physical model parameters
bool is_miscible
 
cs_real_t l_mass_density
 
cs_real_t l_viscosity
 
cs_real_t g_viscosity
 
cs_real_t l_diffusivity_h
 
cs_real_t h_molar_mass
 
cs_real_t ref_temperature
 
cs_real_t henry_constant
 
Numerical parameters
cs_gwf_tpf_solver_type_t solver_type
 Type of solver considered to solve the system of equations (choice of main unknowns and strategy of resolution (coupled/segregated)) More...
 
bool enforce_pg_positivity
 Apply a rescaling on the gas pressure to avoid a negative value. More...
 
cs_gwf_tpf_approx_type_t approx_type
 type of approximation used for the computation of diffusion, unsteady coefficients More...
 
double 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. More...
 
double 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. More...
 
cs_param_nl_algo_t nl_algo_type
 
int nl_algo_verbosity
 
cs_param_convergence_t nl_cvg_param
 
cs_iter_algo_param_aac_t anderson_param
 
cs_iter_algo_tnl_algo
 

Detailed Description

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:

  • No water in the gaseous phase
  • Incompressibility of the liquid phase
  • Hydrogen pressure is given by the "perfect gas" law in the gas phase and the Henry's law in the liquid phase (when a miscible model is used) The two primitive variables are the capillarity and liquid pressures with a specific treatment in the saturated case (cf. the cited article or Angelini's PhD thesis)

Notations are the following :

  • Two phases: Liquid phase denoted by "l" and gaseous phase denoted by "g"
  • indice "c" refers to the capillarity pressure
  • Two components: water denoted by "w" and a gaseous component (let's say hydrogen) denoted by "h". The gaseous component is present in the two phases whereas water is only considered in the liquid phase.

Field Documentation

◆ anderson_param

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

approx_type

type of approximation used for the computation of diffusion, unsteady coefficients

◆ b01_w_eqp

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

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

c_pressure

Pointer to a cs_field_t structure named "capillarity_pressure". Capillarity pressure $ P_c = P_g - P_l $

◆ cell_weight

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.

◆ diff_g_pty

cs_property_t* diff_g_pty

Property related to the relative permeability in the liquid phase

◆ diff_hc_pty

cs_property_t* diff_hc_pty

Property related to the relative permeability in the liquid phase

◆ diff_hl_pty

cs_property_t* diff_hl_pty

Property related to the relative permeability in the liquid phase

◆ diff_wl_pty

cs_property_t* diff_wl_pty

Property related to the relative permeability in the liquid phase

◆ enforce_pg_positivity

enforce_pg_positivity

Apply a rescaling on the gas pressure to avoid a negative value.

◆ g_darcy

g_darcy

Pointer to a cs_gwf_darcy_flux_t structure. Darcy advective flux in the gas phase

◆ g_pressure

g_pressure

Pointer to a cs_field_t structure named "gas_pressure". Pressure in the gas phase is denoted by $ P_g $.

◆ g_rho_h

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

g_viscosity

Viscosity in the gas phase (assumed to be constant) in Pa.s

◆ h_eq

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

h_molar_mass

Molar mass of the main component in the gas phase (e.g. hydrogen) in kg.mol^-1

◆ henry_constant

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

is_miscible

true or false. Some simplifications can be done if immiscible.

◆ krg_pty

krg_pty

Property related to the relative permeability in the gas phase

◆ krl_pty

krl_pty

Property related to the relative permeability in the liquid phase

◆ l_darcy

l_darcy

Pointer to a cs_gwf_darcy_flux_t structure. Darcy advective flux in the liquid phase

◆ l_diffusivity_h

l_diffusivity_h

Molecular diffusivity of the hydrogen in the liquid phase in m^2.s^-1

◆ l_mass_density

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

l_pressure

Pointer to a cs_field_t structure named "liquid_pressure". Pressure in the liquid phase is denoted by $ P_l $.

◆ l_rho_h

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

l_saturation

Pointer to a cs_field_t structure. Liquid saturation at cells. This quantity is denoted by $ S_l $ and is defined by the soil model

◆ l_viscosity

l_viscosity

Viscosity in the liquid phase (assumed to be constant) in Pa.s

◆ lcap_pty

lcap_pty

Property related to the liquid capacity (derivative of the liquid saturation w.r.t. the capillarity pressure)

◆ lsat_pre_array

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

lsat_pty

Property related to the liquid saturation

◆ nl_algo

nl_algo

Structure used to manage the non-linearities

◆ nl_algo_type

nl_algo_type

Type of algorithm to solve the non-linearities

◆ nl_algo_verbosity

nl_algo_verbosity

Level of verbosity for the algorithm solving the non-linearities

◆ nl_cvg_param

nl_cvg_param

Set of parameters to drive the convergence of the non-linear solver

◆ ref_temperature

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

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

srct_h_array

Values of the source terms for the hydrogen conservation equation. Size = n_vertrices

◆ srct_w_array

srct_w_array

Values of the source terms for the water conservation equation. Size = n_vertices

◆ system

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

◆ time_hc_pty

cs_property_t* time_hc_pty

Property related to the relative permeability in the liquid phase

◆ time_hl_pty

cs_property_t* time_hl_pty

Property related to the relative permeability in the liquid phase

◆ time_wc_pty

cs_property_t* time_wc_pty

Property related to the relative permeability in the liquid phase

◆ upwind_weight

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

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.


The documentation for this struct was generated from the following file: