Structure to handle the modelling of miscible two-phase flows in a porous media. More...
#include <cs_gwf_priv.h>
Data Fields | |
Equations | |
cs_equation_t * | w_eq |
cs_equation_t * | h_eq |
Darcy (advection) fields | |
cs_gwf_darcy_flux_t * | l_darcy |
cs_gwf_darcy_flux_t * | g_darcy |
Properties related to the model | |
cs_property_t * | time_w_eq_pty |
cs_property_t * | diff_w_eq_pty |
cs_property_t * | time_h_eq_pty |
cs_property_t * | diff_hl_eq_pty |
cs_property_t * | diff_hg_eq_pty |
Additional fields | |
cs_field_t * | l_saturation |
cs_field_t * | c_pressure |
cs_field_t * | l_pressure |
cs_field_t * | g_pressure |
Additional arrays | |
cs_real_t * | time_w_eq_array |
cs_real_t * | diff_w_eq_array |
cs_real_t * | time_h_eq_array |
cs_real_t * | diff_hl_eq_array |
cs_real_t * | diff_hg_eq_array |
cs_real_t * | l_rel_permeability |
Model parameters | |
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 | w_molar_mass |
cs_real_t | h_molar_mass |
cs_real_t | ref_temperature |
cs_real_t | henry_constant |
Structure to handle the modelling of miscible 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:
The two primitive variables are the liquid and gas pressures with a specific treatment in the saturated case to handle the gaz pressure (cf. the cited article or Angelini's PhD thesis
Notations are the following :
c_pressure |
Pointer to a cs_field_t structure named "capillarity_pressure". Capillarity pressure P_c = P_g - P_l
diff_hg_eq_array |
Values in each cell of the coefficient appearing in the diffusion term for the gas phase in the hydrogen conservation equation. This array is linked to the diff_hl_eq_pty (size = n_cells)
diff_hg_eq_pty |
Property related to the diffusion term of the hydrogen conservation equation (part related to the gas phase)
diff_hl_eq_array |
Values in each cell of the coefficient appearing in the diffusion term for the liquid phase in the hydrogen conservation equation. This array is linked to the diff_hl_eq_pty (size = n_cells)
diff_hl_eq_pty |
Property related to the diffusion term of the hydrogen conservation equation (part related to the liquid phase)
diff_w_eq_array |
Values in each cell of the coefficient appearing in the diffusion term in the water conservation equation. This array is linked to the diff_w_eq_pty (size = n_cells)
diff_w_eq_pty |
Property related to the diffusion term of the water conservation equation
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 P_g.
g_viscosity |
Viscosity in the gas phase (assumed to be constant) in Pa.s
h_eq |
Equation of conservation for the (di)hydrogen. Hydrogen can be present in the liquid or in the gas phase.
h_molar_mass |
Molar mass of the main component in the gas phase (e.g. hydrogen) in kg.mol^-1
cs_real_t henry_constant |
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 P_l.
l_rel_permeability |
Values in each cell of the relative permeability in the liquid phase. This quantity is used either in the water conservation or in the hydrogen conservation. This enables also to recover the (full) permeability since permeability = abs_permeability * rel_l_permeability
l_saturation |
Pointer to a cs_field_t structure. Liquid saturation at cells.
l_viscosity |
Viscosity in the liquid phase (assumed to be constant) in Pa.s
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.
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.
time_h_eq_array |
Values in each cell of the coefficient appearing in front of the unsteady term in the hydrogen conservation equation. This array is linked to the time_h_eq_pty (size = n_cells)
time_h_eq_pty |
Property related to the unsteady term of the hydrogen conservation equation
time_w_eq_array |
Values in each cell of the coefficient appearing in front of the unsteady term in the water conservation equation. This array is linked to the time_w_eq_pty (size = n_cells)
time_w_eq_pty |
Property related to the unsteady term of the water conservation equation
w_eq |
Equation of conservation for the water component. Only the liquid phase is considered. One assumes to wapter vapour in the gas phase.
w_molar_mass |
Molar mass of the main component in the liquid phase (e.g. water) in kg.mol^-1