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

Detailed Description

Definition at line 93 of file WEquation.h.

#include <WEquation.h>

Public Member Functions

void eval (double const dt, ConstituentDensityData const &constituent_density_data, PrevState< ConstituentDensityData > const &constituent_density_data_prev, PrevState< PureLiquidDensityData > const &rho_W_LR_prev, PureLiquidDensityData const &rho_W_LR, SaturationData const &S_L_data, FW3aData &fW_3a) const
 
void dEval (double const dt, ConstituentDensityData const &constituent_density_data, PhaseTransitionData const &phase_transition_data, PrevState< ConstituentDensityData > const &constituent_density_data_prev, PrevState< PureLiquidDensityData > const &rho_W_LR_prev, PureLiquidDensityData const &rho_W_LR, SaturationData const &S_L_data, SaturationDataDeriv const &dS_L_dp_cap, FW3aDerivativeData &dfW_3a) const
 

Member Function Documentation

◆ dEval()

void ProcessLib::TH2M::ConstitutiveRelations::FW3aModel::dEval ( double const dt,
ConstituentDensityData const & constituent_density_data,
PhaseTransitionData const & phase_transition_data,
PrevState< ConstituentDensityData > const & constituent_density_data_prev,
PrevState< PureLiquidDensityData > const & rho_W_LR_prev,
PureLiquidDensityData const & rho_W_LR,
SaturationData const & S_L_data,
SaturationDataDeriv const & dS_L_dp_cap,
FW3aDerivativeData & dfW_3a ) const

Definition at line 132 of file WEquation.cpp.

142{
143 if (dt == 0.)
144 {
145 dfW_3a.dp_GR = 0.;
146 dfW_3a.dp_cap = 0.;
147 dfW_3a.dT = 0.;
148 return;
149 }
150
151 auto const S_L = S_L_data.S_L;
152 auto const S_G = 1. - S_L;
153
154 double const rho_W_GR_dot = (constituent_density_data.rho_W_GR -
155 constituent_density_data_prev->rho_W_GR) /
156 dt;
157 double const rho_W_LR_dot = (rho_W_LR() - **rho_W_LR_prev) / dt;
158
159 dfW_3a.dp_GR = /*(ds_G_dp_GR = 0) * rho_W_GR_dot +*/
160 S_G * phase_transition_data.drho_W_GR_dp_GR / dt +
161 /*(ds_L_dp_GR = 0) * rho_W_LR_dot +*/
162 S_L * phase_transition_data.drho_W_LR_dp_GR / dt;
163
164 dfW_3a.dp_cap = -dS_L_dp_cap() * rho_W_GR_dot +
165 S_G * phase_transition_data.drho_W_GR_dp_cap / dt +
166 dS_L_dp_cap() * rho_W_LR_dot -
167 S_L * phase_transition_data.drho_W_LR_dp_LR / dt;
168 dfW_3a.dT = S_G * phase_transition_data.drho_W_GR_dT / dt +
169 S_L * phase_transition_data.drho_W_LR_dT / dt;
170}

References ProcessLib::TH2M::ConstitutiveRelations::FW3aDerivativeData::dp_cap, ProcessLib::TH2M::ConstitutiveRelations::FW3aDerivativeData::dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_GR_dp_cap, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_GR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_GR_dT, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_LR_dp_GR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_LR_dp_LR, ProcessLib::TH2M::ConstitutiveRelations::PhaseTransitionData::drho_W_LR_dT, ProcessLib::TH2M::ConstitutiveRelations::FW3aDerivativeData::dT, 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::FW3aModel::eval ( double const dt,
ConstituentDensityData const & constituent_density_data,
PrevState< ConstituentDensityData > const & constituent_density_data_prev,
PrevState< PureLiquidDensityData > const & rho_W_LR_prev,
PureLiquidDensityData const & rho_W_LR,
SaturationData const & S_L_data,
FW3aData & fW_3a ) const

Definition at line 108 of file WEquation.cpp.

116{
117 if (dt == 0.)
118 {
119 fW_3a.a = 0;
120 return;
121 }
122
123 double const rho_W_GR_dot = (constituent_density_data.rho_W_GR -
124 constituent_density_data_prev->rho_W_GR) /
125 dt;
126 double const rho_W_LR_dot = (rho_W_LR() - **rho_W_LR_prev) / dt;
127 auto const S_L = S_L_data.S_L;
128 auto const S_G = 1. - S_L;
129 fW_3a.a = S_G * rho_W_GR_dot + S_L * rho_W_LR_dot;
130}

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