OGS
ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGModel< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGModel< DisplacementDim >

Definition at line 125 of file CEquation.h.

#include <CEquation.h>

Public Member Functions

void eval (AdvectionData< DisplacementDim > const &advection_data, FluidDensityData const &fluid_density_data, PhaseTransitionData const &phase_transition_data, PorosityData const &porosity_data, SaturationData const &S_L_data, FC4LCpGData< DisplacementDim > &fC_4_LCpG) const
 
void dEval (PermeabilityData< DisplacementDim > const &permeability_data, ViscosityData const &viscosity_data, PhaseTransitionData const &phase_transition_data, AdvectionDerivativeData< DisplacementDim > const &advection_d_data, FC4LCpGDerivativeData< DisplacementDim > &dfC_4_LCpG) const
 

Member Function Documentation

◆ dEval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGModel< DisplacementDim >::dEval ( PermeabilityData< DisplacementDim > const & permeability_data,
ViscosityData const & viscosity_data,
PhaseTransitionData const & phase_transition_data,
AdvectionDerivativeData< DisplacementDim > const & advection_d_data,
FC4LCpGDerivativeData< DisplacementDim > & dfC_4_LCpG ) const

Definition at line 200 of file CEquation.cpp.

206{
207 dfC_4_LCpG.dp_GR = advection_d_data.dadvection_C_dp_GR
208 // + ddiffusion_C_p_dp_GR TODO (naumov)
209 ;
210
211 dfC_4_LCpG.dp_cap = advection_d_data.dadvection_C_dp_cap
212 // + ddiffusion_C_p_dp_cap TODO (naumov)
213 ;
214
215 GlobalDimMatrix<DisplacementDim> const k_over_mu_G =
216 permeability_data.Ki * permeability_data.k_rel_G / viscosity_data.mu_GR;
217 GlobalDimMatrix<DisplacementDim> const k_over_mu_L =
218 permeability_data.Ki * permeability_data.k_rel_L / viscosity_data.mu_LR;
219
220 dfC_4_LCpG.dT = phase_transition_data.drho_C_GR_dT * k_over_mu_G +
221 phase_transition_data.drho_C_LR_dT * k_over_mu_L
222 // + ddiffusion_C_p_dT TODO (naumov)
223 ;
224}
Eigen::Matrix< double, DisplacementDim, DisplacementDim, Eigen::RowMajor > GlobalDimMatrix
Definition Base.h:32

References ProcessLib::TH2M::ConstitutiveRelations::AdvectionDerivativeData< DisplacementDim >::dadvection_C_dp_cap, ProcessLib::TH2M::ConstitutiveRelations::AdvectionDerivativeData< DisplacementDim >::dadvection_C_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGDerivativeData< DisplacementDim >::dp_cap, ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGDerivativeData< DisplacementDim >::dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_GR_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_LR_dT, ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGDerivativeData< DisplacementDim >::dT, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_G, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_L, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::Ki, ProcessLib::TH2M::ConstitutiveRelations::ViscosityData::mu_GR, and ProcessLib::TH2M::ConstitutiveRelations::ViscosityData::mu_LR.

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGModel< DisplacementDim >::eval ( AdvectionData< DisplacementDim > const & advection_data,
FluidDensityData const & fluid_density_data,
PhaseTransitionData const & phase_transition_data,
PorosityData const & porosity_data,
SaturationData const & S_L_data,
FC4LCpGData< DisplacementDim > & fC_4_LCpG ) const

Definition at line 170 of file CEquation.cpp.

177{
178 GlobalDimMatrix<DisplacementDim> const advection_C =
179 advection_data.advection_C_G + advection_data.advection_C_L;
180
181 double const sD_G = phase_transition_data.diffusion_coefficient_vapour;
182 double const sD_L = phase_transition_data.diffusion_coefficient_solute;
183
184 double const phi_G = (1 - S_L_data.S_L) * porosity_data.phi;
185 double const phi_L = S_L_data.S_L * porosity_data.phi;
186
187 double const diffusion_CGpGR = -phi_G * fluid_density_data.rho_GR * sD_G *
188 phase_transition_data.dxmWG_dpGR;
189 double const diffusion_CLpGR = -phi_L * fluid_density_data.rho_LR * sD_L *
190 phase_transition_data.dxmWL_dpGR;
191
192 double const diffusion_C_pGR = diffusion_CGpGR + diffusion_CLpGR;
193
194 auto const I =
195 Eigen::Matrix<double, DisplacementDim, DisplacementDim>::Identity();
196 fC_4_LCpG.L.noalias() = diffusion_C_pGR * I + advection_C;
197}

References ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_C_G, ProcessLib::TH2M::ConstitutiveRelations::AdvectionData< DisplacementDim >::advection_C_L, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::diffusion_coefficient_solute, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::diffusion_coefficient_vapour, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::dxmWG_dpGR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::dxmWL_dpGR, ProcessLib::TH2M::ConstitutiveRelations::FC4LCpGData< DisplacementDim >::L, ProcessLib::TH2M::ConstitutiveRelations::PorosityData::phi, ProcessLib::TH2M::ConstitutiveRelations::FluidDensityData::rho_GR, ProcessLib::TH2M::ConstitutiveRelations::FluidDensityData::rho_LR, and ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L.


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