[PL] Re-enable p and T dependent density in TRM

ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.cpp

@@ -18,6 +18,8 @@ namespace ProcessLib::ThermoRichardsMechanics
 template <int DisplacementDim>
 void FluidThermalExpansionModel<DisplacementDim>::eval(
     SpaceTimeData const& x_t, MediaData const& media_data,
+    CapillaryPressureData<DisplacementDim> const& p_cap_data,
+    TemperatureData<DisplacementDim> const& T_data,
     SolidThermalExpansionData<DisplacementDim> const& s_therm_exp_data,
     PorosityData const& poro_data, LiquidDensityData const& rho_L_data,
     BiotData const& biot_data, FluidThermalExpansionData& out) const
@@ -28,6 +30,9 @@ void FluidThermalExpansionModel<DisplacementDim>::eval(
 
     namespace MPL = MaterialPropertyLib;
     MPL::VariableArray variables;
+    variables[static_cast<int>(MPL::Variable::phase_pressure)] =
+        -p_cap_data.p_cap;
+    variables[static_cast<int>(MPL::Variable::temperature)] = -T_data.T;
 
     double const phi = poro_data.phi;
     double const alpha = biot_data.alpha;

ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.h

@@ -27,6 +27,8 @@ struct FluidThermalExpansionModel
 {
     void eval(
         SpaceTimeData const& x_t, MediaData const& media_data,
+        CapillaryPressureData<DisplacementDim> const& p_cap_data,
+        TemperatureData<DisplacementDim> const& T_data,
         SolidThermalExpansionData<DisplacementDim> const& s_therm_exp_data,
         PorosityData const& poro_data, LiquidDensityData const& rho_L_data,
         BiotData const& biot_data, FluidThermalExpansionData& out) const;

ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.cpp

@@ -144,8 +144,9 @@ void ConstitutiveSetting<DisplacementDim>::eval(
                                       dS_L_data, poro_data, p_cap_data, T_data,
                                       vap_data);
 
-    models.f_therm_exp_model.eval(x_t, media_data, s_therm_exp_data, poro_data,
-                                  rho_L_data, biot_data, f_therm_exp_data);
+    models.f_therm_exp_model.eval(x_t, media_data, p_cap_data, T_data,
+                                  s_therm_exp_data, poro_data, rho_L_data,
+                                  biot_data, f_therm_exp_data);
 
     models.storage_model.eval(x_t, biot_data, poro_data, rho_L_data, S_L_data,
                               dS_L_data, prev_state.S_L_data, p_cap_data,