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

Detailed Description

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

Definition at line 42 of file TH2M/ConstitutiveRelations/Porosity.h.

#include <Porosity.h>

Public Member Functions

void eval (SpaceTimeData const &x_t, MediaData const &media_data, SaturationData const &S_L_data, PrevState< SaturationData > const &S_L_prev_data, CapillaryPressureData const &p_cap, GasPressureData const &p_GR, BishopsData const &chi_S_L, PrevState< BishopsData > const &chi_S_L_prev, SolidCompressibilityData const &solid_compressibility, StrainData< DisplacementDim > const &eps_data, MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const &&eps_prev, PrevState< PorosityData > const &porosity_prev_data, PorosityData &porosity_data) const
void dEval (SpaceTimeData const &x_t, MediaData const &media_data, SaturationData const &S_L_data, PrevState< SaturationData > const &S_L_prev_data, CapillaryPressureData const &p_cap, GasPressureData const &p_GR, BishopsData const &chi_S_L, PrevState< BishopsData > const &chi_S_L_prev, SolidCompressibilityData const &solid_compressibility, StrainData< DisplacementDim > const &eps_data, MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const &&eps_prev, PrevState< PorosityData > const &porosity_prev_data, PorosityDerivativeData &porosity_d_data) const

Member Function Documentation

◆ dEval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::PorosityModel< DisplacementDim >::dEval ( SpaceTimeData const & x_t,
MediaData const & media_data,
SaturationData const & S_L_data,
PrevState< SaturationData > const & S_L_prev_data,
CapillaryPressureData const & p_cap,
GasPressureData const & p_GR,
BishopsData const & chi_S_L,
PrevState< BishopsData > const & chi_S_L_prev,
SolidCompressibilityData const & solid_compressibility,
StrainData< DisplacementDim > const & eps_data,
MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const && eps_prev,
PrevState< PorosityData > const & porosity_prev_data,
PorosityDerivativeData & porosity_d_data ) const

Definition at line 55 of file TH2M/ConstitutiveRelations/Porosity.cpp.

66{
67 static int const KelvinVectorSize =
68 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
69 using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>;
70
71 MaterialPropertyLib::VariableArray variables;
72 MaterialPropertyLib::VariableArray variables_prev;
73
74 variables.grain_compressibility = solid_compressibility();
75
76 variables.liquid_saturation = S_L_data.S_L;
77 variables_prev.liquid_saturation = S_L_prev_data->S_L;
78
79 variables.effective_pore_pressure =
80 (1 - chi_S_L.chi_S_L) * p_GR.pG +
81 chi_S_L.chi_S_L * (p_GR.pG - p_cap.pCap);
82
83 // Used in MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
84 // and MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
85 variables_prev.effective_pore_pressure =
86 (1 - chi_S_L_prev->chi_S_L) * p_GR.pG_prev +
87 chi_S_L_prev->chi_S_L * (p_GR.pG_prev - p_cap.pCap_prev);
88
89 variables.volumetric_strain = Invariants::trace(eps_data.eps);
90 variables_prev.volumetric_strain = Invariants::trace(eps_prev);
91
92 variables_prev.porosity = porosity_prev_data->phi;
93
94 auto const& mpl_porosity = media_data.porosity_prop;
95
96 porosity_d_data.dphi_dT = mpl_porosity.template dValue<double>(
97 variables, variables_prev, MaterialPropertyLib::Variable::temperature,
98 x_t.x, x_t.t, x_t.dt);
99
100 porosity_d_data.dphi_L_dp_cap = mpl_porosity.template dValue<double>(
101 variables, variables_prev,
102 MaterialPropertyLib::Variable::capillary_pressure, x_t.x, x_t.t,
103 x_t.dt);
104}
MathLib::KelvinVector::kelvin_vector_dimensions
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Definition KelvinVector.h:18

References MaterialPropertyLib::capillary_pressure, ProcessLib::TH2M::ConstitutiveRelations::BishopsData::chi_S_L, ProcessLib::TH2M::ConstitutiveRelations::PorosityDerivativeData::dphi_dT, ProcessLib::TH2M::ConstitutiveRelations::PorosityDerivativeData::dphi_L_dp_cap, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, MaterialPropertyLib::VariableArray::effective_pore_pressure, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, MaterialPropertyLib::VariableArray::grain_compressibility, MathLib::KelvinVector::kelvin_vector_dimensions(), MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap_prev, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG_prev, MaterialPropertyLib::VariableArray::porosity, ProcessLib::TH2M::ConstitutiveRelations::MediaData::porosity_prop, ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MaterialPropertyLib::temperature, MaterialPropertyLib::VariableArray::volumetric_strain, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::PorosityModel< DisplacementDim >::eval ( SpaceTimeData const & x_t,
MediaData const & media_data,
SaturationData const & S_L_data,
PrevState< SaturationData > const & S_L_prev_data,
CapillaryPressureData const & p_cap,
GasPressureData const & p_GR,
BishopsData const & chi_S_L,
PrevState< BishopsData > const & chi_S_L_prev,
SolidCompressibilityData const & solid_compressibility,
StrainData< DisplacementDim > const & eps_data,
MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const && eps_prev,
PrevState< PorosityData > const & porosity_prev_data,
PorosityData & porosity_data ) const

Definition at line 12 of file TH2M/ConstitutiveRelations/Porosity.cpp.

23{
24 static int const KelvinVectorSize =
25 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
26 using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>;
27
28 MaterialPropertyLib::VariableArray variables;
29 MaterialPropertyLib::VariableArray variables_prev;
30
31 variables.grain_compressibility = solid_compressibility();
32
33 variables.liquid_saturation = S_L_data.S_L;
34 variables_prev.liquid_saturation = S_L_prev_data->S_L;
35
36 variables.effective_pore_pressure =
37 (1 - chi_S_L.chi_S_L) * p_GR.pG +
38 chi_S_L.chi_S_L * (p_GR.pG - p_cap.pCap);
39
40 // Used in MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
41 // and MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
42 variables_prev.effective_pore_pressure =
43 (1 - chi_S_L_prev->chi_S_L) * p_GR.pG_prev +
44 chi_S_L_prev->chi_S_L * (p_GR.pG_prev - p_cap.pCap_prev);
45
46 variables.volumetric_strain = Invariants::trace(eps_data.eps);
47 variables_prev.volumetric_strain = Invariants::trace(eps_prev);
48
49 variables_prev.porosity = porosity_prev_data->phi;
50 porosity_data.phi = media_data.porosity_prop.template value<double>(
51 variables, variables_prev, x_t.x, x_t.t, x_t.dt);
52}

References ProcessLib::TH2M::ConstitutiveRelations::BishopsData::chi_S_L, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, MaterialPropertyLib::VariableArray::effective_pore_pressure, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, MaterialPropertyLib::VariableArray::grain_compressibility, MathLib::KelvinVector::kelvin_vector_dimensions(), MaterialPropertyLib::VariableArray::liquid_saturation, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap_prev, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG_prev, ProcessLib::TH2M::ConstitutiveRelations::PorosityData::phi, MaterialPropertyLib::VariableArray::porosity, ProcessLib::TH2M::ConstitutiveRelations::MediaData::porosity_prop, ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MaterialPropertyLib::VariableArray::volumetric_strain, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

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