32 auto const& mpl_thermal_conductivity =
33 media_data.
medium[MPL::PropertyType::thermal_conductivity];
35 thermal_conductivity_data.
lambda = MPL::formEigenTensor<DisplacementDim>(
36 mpl_thermal_conductivity.value(variables, x_t.
x, x_t.
t, x_t.
dt));
42 media_data.
gas.
hasProperty(MPL::PropertyType::thermal_conductivity)
43 ? MPL::formEigenTensor<DisplacementDim>(
44 media_data.
gas[MPL::PropertyType::thermal_conductivity].value(
45 variables, x_t.
x, x_t.
t, x_t.
dt))
46 : MPL::formEigenTensor<DisplacementDim>(0.);
48 auto const dlambda_GR_dT =
49 media_data.
gas.
hasProperty(MPL::PropertyType::thermal_conductivity)
50 ? MPL::formEigenTensor<DisplacementDim>(
51 media_data.
gas[MPL::PropertyType::thermal_conductivity]
52 .dValue(variables, MPL::Variable::temperature, x_t.
x,
54 : MPL::formEigenTensor<DisplacementDim>(0.);
58 ? MPL::formEigenTensor<DisplacementDim>(
59 media_data.
liquid[MPL::PropertyType::thermal_conductivity]
60 .value(variables, x_t.
x, x_t.
t, x_t.
dt))
61 : MPL::formEigenTensor<DisplacementDim>(0.);
63 auto const dlambda_LR_dT =
65 ? MPL::formEigenTensor<DisplacementDim>(
66 media_data.
liquid[MPL::PropertyType::thermal_conductivity]
67 .dValue(variables, MPL::Variable::temperature, x_t.
x,
69 : MPL::formEigenTensor<DisplacementDim>(0.);
73 ? MPL::formEigenTensor<DisplacementDim>(
74 media_data.
solid[MPL::PropertyType::thermal_conductivity]
75 .value(variables, x_t.
x, x_t.
t, x_t.
dt))
76 : MPL::formEigenTensor<DisplacementDim>(0.);
78 auto const dlambda_SR_dT =
80 ? MPL::formEigenTensor<DisplacementDim>(
81 media_data.
solid[MPL::PropertyType::thermal_conductivity]
82 .dValue(variables, MPL::Variable::temperature, x_t.
x,
84 : MPL::formEigenTensor<DisplacementDim>(0.);
87 double const dphi_G_dp_cap = -dS_L_dp_cap() * porosity_data.
phi;
89 double const dphi_L_dp_cap = dS_L_dp_cap() * porosity_data.
phi;
92 dphi_G_dp_cap * lambdaGR + dphi_L_dp_cap * lambdaLR;
94 double const phi_L = S_L_data.
S_L * porosity_data.
phi;
95 double const phi_G = (1. - S_L_data.
S_L) * porosity_data.
phi;
96 double const phi_S = 1. - porosity_data.
phi;
103 (1 - S_L_data.
S_L) * porosity_d_data.
dphi_dT * lambdaGR +
104 phi_G * dlambda_GR_dT +
105 S_L_data.
S_L * porosity_d_data.
dphi_dT * lambdaLR +
106 +phi_L * dlambda_LR_dT - porosity_d_data.
dphi_dT * lambdaSR +
107 phi_S * dlambda_SR_dT;
void eval(SpaceTimeData const &x_t, MediaData const &media_data, TemperatureData const &T_data, PorosityData const &porosity_data, PorosityDerivativeData const &porosity_d_data, SaturationData const &S_L_data, SaturationDataDeriv const &dS_L_dp_cap, ThermalConductivityData< DisplacementDim > &thermal_conductivity_data) const