[PL] Pass volumetric strain to permeability model.

ufz · Sep 15, 2023 · 530123c · 530123c
1 parent 26c6f2a
commit 530123c
Show file tree

Hide file tree

Showing 4 changed files with 12 additions and 4 deletions.
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/PermeabilityModel.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/PermeabilityModel.cpp
@@ -23,11 +23,16 @@ void PermeabilityModel<DisplacementDim>::eval(
     LiquidViscosityData const& mu_L_data,
     TransportPorosityData const& transport_poro_data,
     TotalStressData<DisplacementDim> const& total_stress_data,
+    StrainData<DisplacementDim> const& eps_data,
     EquivalentPlasticStrainData const& equiv_plast_strain_data,
     PermeabilityData<DisplacementDim>& out) const
 {
     namespace MPL = MaterialPropertyLib;
 
+    static constexpr int kelvin_vector_size =
+        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
+    using Invariants = MathLib::KelvinVector::Invariants<kelvin_vector_size>;
+
     auto const& medium = media_data.medium;
 
     MPL::VariableArray variables;
@@ -57,6 +62,8 @@ void PermeabilityModel<DisplacementDim>::eval(
     variables.equivalent_plastic_strain =
         equiv_plast_strain_data.equivalent_plastic_strain;
 
+    variables.volumetric_strain = Invariants::trace(eps_data.eps);
+
     auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>(
         medium.property(MPL::PropertyType::permeability)
             .value(variables, x_t.x, x_t.t, x_t.dt));

diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/PermeabilityModel.h b/ProcessLib/ThermoRichardsMechanics/ConstitutiveCommon/PermeabilityModel.h
@@ -28,6 +28,7 @@ struct PermeabilityModel
               LiquidViscosityData const& mu_L_data,
               TransportPorosityData const& transport_poro_data,
               TotalStressData<DisplacementDim> const& total_stress_data,
+              StrainData<DisplacementDim> const& eps_data,
               EquivalentPlasticStrainData const& equiv_plast_strain_data,
               PermeabilityData<DisplacementDim>& out) const;
 };

diff --git a/...sMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ConstitutiveSetting.cpp b/...sMechanics/ConstitutiveStressSaturation_StrainPressureTemperature/ConstitutiveSetting.cpp
@@ -127,8 +127,8 @@ void ConstitutiveSetting<DisplacementDim>::eval(
     assertEvalArgsUnique(models.perm_model);
     models.perm_model.eval(x_t, media_data, S_L_data, p_cap_data, T_data,
                            mu_L_data, state.transport_poro_data,
-                           state.total_stress_data, tmp.equiv_plast_strain_data,
-                           perm_data);
+                           state.total_stress_data, state.eps_data,
+                           tmp.equiv_plast_strain_data, perm_data);
 
     assertEvalArgsUnique(models.th_osmosis_model);
     models.th_osmosis_model.eval(x_t, media_data, T_data, rho_L_data,

diff --git a/...sLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ConstitutiveSetting.cpp b/...sLib/ThermoRichardsMechanics/ConstitutiveStress_StrainTemperature/ConstitutiveSetting.cpp
@@ -149,8 +149,8 @@ void ConstitutiveSetting<DisplacementDim>::eval(
     assertEvalArgsUnique(models.perm_model);
     models.perm_model.eval(x_t, media_data, S_L_data, p_cap_data, T_data,
                            mu_L_data, state.transport_poro_data,
-                           cd.total_stress_data, tmp.equiv_plast_strain_data,
-                           perm_data);
+                           cd.total_stress_data, state.eps_data,
+                           tmp.equiv_plast_strain_data, perm_data);
 
     assertEvalArgsUnique(models.th_osmosis_model);
     models.th_osmosis_model.eval(x_t, media_data, T_data, rho_L_data,