OGS
List of all members | Public Member Functions
ProcessLib::TH2M::ConstitutiveRelations::FW2Model Struct Reference

Detailed Description

Definition at line 57 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 48 of file WEquation.cpp.

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

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::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 29 of file WEquation.cpp.

37{
38 auto const S_L = S_L_data.S_L;
39 auto const S_G = 1. - S_L;
40 double const rho_W_FR =
41 S_G * constituent_density_data.rho_W_GR + S_L * rho_W_LR();
42
43 fW_2.a =
44 porosity_data.phi * (rho_W_LR() - constituent_density_data.rho_W_GR) -
45 rho_W_FR * pCap() * (biot_data() - porosity_data.phi) * beta_p_SR();
46}

References ProcessLib::TH2M::ConstitutiveRelations::FW2Data::a, 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: