diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h
index d92822dd710290f2cb2972d3abbc28837d1bc53c..804836e85fa7c41dfa58683d3027ebc4d74fc3b0 100644
--- a/ProcessLib/HT/HTFEM.h
+++ b/ProcessLib/HT/HTFEM.h
@@ -130,11 +130,9 @@ public:
         auto const& medium =
             *_process_data.media_map->getMedium(_element.getID());
         auto const& liquid_phase = medium.phase("AqueousLiquid");
-        auto const& solid_phase = medium.phase("Solid");
 
         auto const K = MaterialPropertyLib::formEigenTensor<GlobalDim>(
-            solid_phase
-                .property(MaterialPropertyLib::PropertyType::permeability)
+            medium.property(MaterialPropertyLib::PropertyType::permeability)
                 .value(vars, pos, t));
 
         auto const mu =
@@ -273,7 +271,6 @@ protected:
         auto const& medium =
             *_process_data.media_map->getMedium(_element.getID());
         auto const& liquid_phase = medium.phase("AqueousLiquid");
-        auto const& solid_phase = medium.phase("Solid");
 
         for (unsigned ip = 0; ip < n_integration_points; ++ip)
         {
@@ -294,8 +291,7 @@ protected:
                 MaterialPropertyLib::Variable::phase_pressure)] = p_int_pt;
 
             auto const K = MaterialPropertyLib::formEigenTensor<GlobalDim>(
-                solid_phase
-                    .property(MaterialPropertyLib::PropertyType::permeability)
+                medium.property(MaterialPropertyLib::PropertyType::permeability)
                     .value(vars, pos, t));
 
             auto const mu =
diff --git a/ProcessLib/HT/MonolithicHTFEM.h b/ProcessLib/HT/MonolithicHTFEM.h
index ce0e7da613d154554137d755719dce42b486afdc..18a9a47a59eb54ac094bfec54ed0a17fd741b1b4 100644
--- a/ProcessLib/HT/MonolithicHTFEM.h
+++ b/ProcessLib/HT/MonolithicHTFEM.h
@@ -143,7 +143,7 @@ public:
 
             auto const intrinsic_permeability =
                 MaterialPropertyLib::formEigenTensor<GlobalDim>(
-                    solid_phase
+                    medium
                         .property(
                             MaterialPropertyLib::PropertyType::permeability)
                         .value(vars, pos, t));
diff --git a/ProcessLib/HT/StaggeredHTFEM-impl.h b/ProcessLib/HT/StaggeredHTFEM-impl.h
index 2d94a53c4bf1590f2667f987e1f811c4713253cb..ddba8cbc16040ca6b04d5bb49e0de7cb1b09de2c 100644
--- a/ProcessLib/HT/StaggeredHTFEM-impl.h
+++ b/ProcessLib/HT/StaggeredHTFEM-impl.h
@@ -129,8 +129,7 @@ void StaggeredHTFEM<ShapeFunction, IntegrationMethod, GlobalDim>::
 
         auto const intrinsic_permeability =
             MaterialPropertyLib::formEigenTensor<GlobalDim>(
-                solid_phase
-                    .property(MaterialPropertyLib::PropertyType::permeability)
+                medium.property(MaterialPropertyLib::PropertyType::permeability)
                     .value(vars, pos, t));
         GlobalDimMatrixType const K_over_mu =
             intrinsic_permeability / viscosity;
@@ -267,7 +266,7 @@ void StaggeredHTFEM<ShapeFunction, IntegrationMethod, GlobalDim>::
 
         auto const intrinsic_permeability =
             MaterialPropertyLib::formEigenTensor<GlobalDim>(
-                solid_phase
+                medium
                     .property(MaterialPropertyLib::PropertyType::permeability)
                     .value(vars, pos, t));