OGS
ProcessLib::TH2M::ConstitutiveRelations::FW2Model Struct Reference

Detailed Description

Definition at line 51 of file WEquation.h.

#include <WEquation.h>

Public Member Functions

void eval (BiotData const biot_data, CapillaryPressureData const pCap, ConstituentDensityData const &constituent_density_data, PorosityData const &porosity_data, PureLiquidDensityData const &rho_W_LR, SaturationData const &S_L_data, SolidCompressibilityData const beta_p_SR, FW2Data &fW_2) const
void dEval (BiotData const &biot_data, CapillaryPressureData const pCap, ConstituentDensityData const &constituent_density_data, PhaseTransitionData const &phase_transition_data, PorosityData const &porosity_data, PorosityDerivativeData const &porosity_d_data, PureLiquidDensityData const &rho_W_LR, SaturationData const &S_L_data, SaturationDataDeriv const &dS_L_dp_cap, SolidCompressibilityData const &beta_p_SR, FW2DerivativeData &dfW_2) const

Member Function Documentation

◆ dEval()

void ProcessLib::TH2M::ConstitutiveRelations::FW2Model::dEval ( BiotData const & biot_data,
CapillaryPressureData const pCap,
ConstituentDensityData const & constituent_density_data,
PhaseTransitionData const & phase_transition_data,
PorosityData const & porosity_data,
PorosityDerivativeData const & porosity_d_data,
PureLiquidDensityData const & rho_W_LR,
SaturationData const & S_L_data,
SaturationDataDeriv const & dS_L_dp_cap,
SolidCompressibilityData const & beta_p_SR,
FW2DerivativeData & dfW_2 ) const

Definition at line 42 of file WEquation.cpp.

53{
54 double const S_L = S_L_data.S_L;
55 double const S_G = 1. - S_L;
56
57 double const drho_C_FR_dp_GR =
58 /*(dS_G_dp_GR = 0) * constituent_density_data.rho_C_GR +*/
59 S_G * phase_transition_data.drho_C_GR_dp_GR +
60 /*(dS_L_dp_GR = 0) * constituent_density_data.rho_C_LR +*/
61 S_L * phase_transition_data.drho_C_LR_dp_GR;
62
63 dfW_2.dp_GR = -porosity_data.phi * phase_transition_data.drho_C_GR_dp_GR -
64 drho_C_FR_dp_GR * pCap.pCap *
65 (biot_data() - porosity_data.phi) * beta_p_SR();
66
67 double const dfW_2a_dp_cap =
68 porosity_data.phi * (-phase_transition_data.drho_W_LR_dp_LR -
69 phase_transition_data.drho_W_GR_dp_cap);
70 double const rho_W_FR =
71 S_G * constituent_density_data.rho_W_GR + S_L * rho_W_LR();
72
73 double const drho_W_FR_dp_cap =
74 -dS_L_dp_cap() * constituent_density_data.rho_W_GR +
75 S_G * phase_transition_data.drho_W_GR_dp_cap +
76 dS_L_dp_cap() * rho_W_LR() -
77 S_L * phase_transition_data.drho_W_LR_dp_LR;
78
79 double const dfW_2b_dp_cap =
80 drho_W_FR_dp_cap * pCap.pCap * (biot_data() - porosity_data.phi) *
81 beta_p_SR() +
82 rho_W_FR * (biot_data() - porosity_data.phi) * beta_p_SR();
83
84 dfW_2.dp_cap = dfW_2a_dp_cap - dfW_2b_dp_cap;
85
86 double const drho_W_FR_dT = S_G * phase_transition_data.drho_W_GR_dT +
87 S_L * phase_transition_data.drho_W_LR_dT;
88
89 double const dfW_2a_dT =
90 porosity_d_data.dphi_dT *
91 (rho_W_LR() - constituent_density_data.rho_W_GR) +
92 porosity_data.phi * (phase_transition_data.drho_W_LR_dT -
93 phase_transition_data.drho_W_GR_dT);
94 double const dfW_2b_dT =
95 drho_W_FR_dT * pCap.pCap * (biot_data() - porosity_data.phi) *
96 beta_p_SR() -
97 rho_W_FR * pCap.pCap * porosity_d_data.dphi_dT * beta_p_SR();
98
99 dfW_2.dT = dfW_2a_dT - dfW_2b_dT;
100}

References ProcessLib::TH2M::ConstitutiveRelations::FW2DerivativeData::dp_cap, ProcessLib::TH2M::ConstitutiveRelations::FW2DerivativeData::dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PorosityDerivativeData::dphi_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_GR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_C_LR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_GR_dp_cap, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_GR_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_LR_dp_LR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_LR_dT, ProcessLib::TH2M::ConstitutiveRelations::FW2DerivativeData::dT, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, ProcessLib::TH2M::ConstitutiveRelations::PorosityData::phi, ProcessLib::TH2M::ConstitutiveRelations::ConstituentDensityData::rho_W_GR, and ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L.

Referenced by ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::updateConstitutiveVariablesDerivatives().

◆ eval()

void ProcessLib::TH2M::ConstitutiveRelations::FW2Model::eval ( BiotData const biot_data,
CapillaryPressureData const pCap,
ConstituentDensityData const & constituent_density_data,
PorosityData const & porosity_data,
PureLiquidDensityData const & rho_W_LR,
SaturationData const & S_L_data,
SolidCompressibilityData const beta_p_SR,
FW2Data & fW_2 ) const

Definition at line 23 of file WEquation.cpp.

31{
32 auto const S_L = S_L_data.S_L;
33 auto const S_G = 1. - S_L;
34 double const rho_W_FR =
35 S_G * constituent_density_data.rho_W_GR + S_L * rho_W_LR();
36
37 fW_2.a =
38 porosity_data.phi * (rho_W_LR() - constituent_density_data.rho_W_GR) -
39 rho_W_FR * pCap.pCap * (biot_data() - porosity_data.phi) * beta_p_SR();
40}

References ProcessLib::TH2M::ConstitutiveRelations::FW2Data::a, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, ProcessLib::TH2M::ConstitutiveRelations::PorosityData::phi, ProcessLib::TH2M::ConstitutiveRelations::ConstituentDensityData::rho_W_GR, and ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L.

Referenced by ProcessLib::TH2M::TH2MLocalAssembler< ShapeFunctionDisplacement, ShapeFunctionPressure, DisplacementDim >::updateConstitutiveVariables().


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