Loading [MathJax]/jax/output/HTML-CSS/config.js
OGS
ProcessLib::TH2M::ConstitutiveRelations::PorosityModel< DisplacementDim > Struct Template Reference

Detailed Description

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

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

74{
75 static int const KelvinVectorSize =
78
81
82 variables.grain_compressibility = solid_compressibility();
83
84 variables.liquid_saturation = S_L_data.S_L;
85 variables_prev.liquid_saturation = S_L_prev_data->S_L;
86
87 variables.effective_pore_pressure =
88 (1 - chi_S_L.chi_S_L) * p_GR.pG +
89 chi_S_L.chi_S_L * (p_GR.pG - p_cap.pCap);
90
91 // Used in MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
92 // and MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
93 variables_prev.effective_pore_pressure =
94 (1 - chi_S_L_prev->chi_S_L) * p_GR.pG_prev +
95 chi_S_L_prev->chi_S_L * (p_GR.pG_prev - p_cap.pCap_prev);
96
97 variables.volumetric_strain = Invariants::trace(eps_data.eps);
98 variables_prev.volumetric_strain = Invariants::trace(eps_prev);
99
100 variables_prev.porosity = porosity_prev_data->phi;
101
102 auto const& mpl_porosity =
104
105 porosity_d_data.dphi_dT = mpl_porosity.template dValue<double>(
106 variables, variables_prev, MaterialPropertyLib::Variable::temperature,
107 x_t.x, x_t.t, x_t.dt);
108
109 porosity_d_data.dphi_L_dp_cap = mpl_porosity.template dValue<double>(
110 variables, variables_prev,
112 x_t.dt);
113}
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.

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::MediaData::medium, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap_prev, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG_prev, MaterialPropertyLib::porosity, MaterialPropertyLib::VariableArray::porosity, 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 18 of file Porosity.cpp.

29{
30 static int const KelvinVectorSize =
33
36
37 variables.grain_compressibility = solid_compressibility();
38
39 variables.liquid_saturation = S_L_data.S_L;
40 variables_prev.liquid_saturation = S_L_prev_data->S_L;
41
42 variables.effective_pore_pressure =
43 (1 - chi_S_L.chi_S_L) * p_GR.pG +
44 chi_S_L.chi_S_L * (p_GR.pG - p_cap.pCap);
45
46 // Used in MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
47 // and MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
48 variables_prev.effective_pore_pressure =
49 (1 - chi_S_L_prev->chi_S_L) * p_GR.pG_prev +
50 chi_S_L_prev->chi_S_L * (p_GR.pG_prev - p_cap.pCap_prev);
51
52 variables.volumetric_strain = Invariants::trace(eps_data.eps);
53 variables_prev.volumetric_strain = Invariants::trace(eps_prev);
54
55 variables_prev.porosity = porosity_prev_data->phi;
56 porosity_data.phi =
57 media_data.medium.property(MaterialPropertyLib::PropertyType::porosity)
58 .template value<double>(variables, variables_prev, x_t.x, x_t.t,
59 x_t.dt);
60}

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::MediaData::medium, 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::porosity, MaterialPropertyLib::VariableArray::porosity, MaterialPropertyLib::Medium::property(), 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: