31 static constexpr int kelvin_vector_size =
35 auto const& medium = media_data.
medium;
43 out.
k_rel = medium.property(MPL::PropertyType::relative_permeability)
44 .template value<double>(variables, x_t.
x, x_t.
t, x_t.
dt);
46 medium.property(MPL::PropertyType::relative_permeability)
47 .template dValue<double>(variables,
48 MPL::Variable::liquid_saturation,
54 using SymmetricTensor =
66 auto const K_intrinsic = MPL::formEigenTensor<DisplacementDim>(
67 medium.property(MPL::PropertyType::permeability)
68 .value(variables, x_t.
x, x_t.
t, x_t.
dt));
void eval(SpaceTimeData const &x_t, MediaData const &media_data, SaturationData const &S_L_data, CapillaryPressureData< DisplacementDim > const &p_cap_data, TemperatureData< DisplacementDim > const &T_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