OGS
List of all members | Public Member Functions
ProcessLib::TH2M::ConstitutiveRelations::PorosityModelNonConstantSolidPhaseVolumeFraction< DisplacementDim > Struct Template Reference

Detailed Description

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

Definition at line 49 of file Porosity.h.

#include <Porosity.h>

Public Member Functions

void eval (SpaceTimeData const &x_t, MediaData const &media_data, BiotData const &biot, StrainData< DisplacementDim > const &strain_data, SolidThermalExpansionData< DisplacementDim > const &s_therm_exp_data, PorosityData &porosity_data, PorosityDerivativeData &porosity_d_data) const
 

Member Function Documentation

◆ eval()

template<int DisplacementDim>
void ProcessLib::TH2M::ConstitutiveRelations::PorosityModelNonConstantSolidPhaseVolumeFraction< DisplacementDim >::eval ( SpaceTimeData const & x_t,
MediaData const & media_data,
BiotData const & biot,
StrainData< DisplacementDim > const & strain_data,
SolidThermalExpansionData< DisplacementDim > const & s_therm_exp_data,
PorosityData & porosity_data,
PorosityDerivativeData & porosity_d_data ) const

Definition at line 36 of file Porosity.cpp.

44{
45 MaterialPropertyLib::VariableArray variables;
46
47 auto const& mpl_porosity =
48 media_data.medium[MaterialPropertyLib::PropertyType::porosity];
49
50 double const phi_0 =
51 mpl_porosity.template value<double>(variables, x_t.x, x_t.t, x_t.dt);
52
53 static int const KelvinVectorSize =
54 MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
55 using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>;
56 double const div_u = Invariants::trace(strain_data.eps);
57
58 double const phi_S =
59 (1. - phi_0) *
60 (1. + s_therm_exp_data.thermal_volume_strain - biot() * div_u);
61
62 porosity_data.phi = 1. - phi_S;
63
64 auto const dphi_0_dT = mpl_porosity.template dValue<double>(
65 variables, MaterialPropertyLib::Variable::temperature, x_t.x, x_t.t,
66 x_t.dt);
67
68 porosity_d_data.dphi_dT =
69 dphi_0_dT *
70 (1. + s_therm_exp_data.thermal_volume_strain - biot() * div_u) -
71 (1. - phi_0) * s_therm_exp_data.beta_T_SR;
72}
MathLib::KelvinVector::kelvin_vector_dimensions
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
Definition KelvinVector.h:24

References ProcessLib::TH2M::ConstitutiveRelations::SolidThermalExpansionData< DisplacementDim >::beta_T_SR, ProcessLib::TH2M::ConstitutiveRelations::PorosityDerivativeData::dphi_dT, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, MathLib::KelvinVector::kelvin_vector_dimensions(), ProcessLib::TH2M::ConstitutiveRelations::MediaData::medium, ProcessLib::TH2M::ConstitutiveRelations::PorosityData::phi, MaterialPropertyLib::porosity, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MaterialPropertyLib::temperature, ProcessLib::TH2M::ConstitutiveRelations::SolidThermalExpansionData< DisplacementDim >::thermal_volume_strain, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

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