diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h index bf6372035e13d6e153d59a76095f3371412afa09..dda388621f7e1e2bfb8166c2e736f024e4d96678 100644 --- a/ProcessLib/HT/HTFEM.h +++ b/ProcessLib/HT/HTFEM.h @@ -14,6 +14,8 @@ #include "HTMaterialProperties.h" +#include "MaterialLib/MPL/Medium.h" +#include "MaterialLib/MPL/PropertyType.h" #include "NumLib/DOF/DOFTableUtil.h" #include "NumLib/Extrapolation/ExtrapolatableElement.h" #include "NumLib/Fem/FiniteElement/TemplateIsoparametric.h" @@ -187,8 +189,16 @@ protected: { auto const thermal_conductivity_solid = _material_properties.thermal_conductivity_solid(t, pos)[0]; + + auto const& medium = + *_material_properties.media_map->getMedium(_element.getID()); + auto const& liquid_phase = medium.phase("AqueousLiquid"); auto const thermal_conductivity_fluid = - _material_properties.thermal_conductivity_fluid(t, pos)[0]; + liquid_phase + .property( + MaterialPropertyLib::PropertyType::thermal_conductivity) + .template value<double>(); + double const thermal_conductivity = thermal_conductivity_solid * (1 - porosity) + thermal_conductivity_fluid * porosity;