diff --git a/dumux/implicit/1p2c/1p2cfluxvariables.hh b/dumux/implicit/1p2c/1p2cfluxvariables.hh
index 9dd73c15635172005f11af43b8b22105e5cabd7d..eca03684d0d2abb5d61c84759bb105e23b8ad148 100644
--- a/dumux/implicit/1p2c/1p2cfluxvariables.hh
+++ b/dumux/implicit/1p2c/1p2cfluxvariables.hh
@@ -363,9 +363,17 @@ protected:
         }
         else
         {
+            const Element& elementI = *fvGeometry_.neighbors[face().i];
+            FVElementGeometry fvGeometryI;
+            fvGeometryI.subContVol[0].global = elementI.geometry().center();
+            
+            const Element& elementJ = *fvGeometry_.neighbors[face().j];
+            FVElementGeometry fvGeometryJ;
+            fvGeometryJ.subContVol[0].global = elementJ.geometry().center();
+            
             sp.meanK(K_,
-                     sp.elementIntrinsicPermeability(*fvGeometry_.neighbors[face().i]),
-                     sp.elementIntrinsicPermeability(*fvGeometry_.neighbors[face().j]));
+                     sp.intrinsicPermeability(elementI, fvGeometryI, 0),
+                     sp.intrinsicPermeability(elementJ, fvGeometryJ, 0));
         }
     }
 
diff --git a/dumux/implicit/box/boxforchheimerfluxvariables.hh b/dumux/implicit/box/boxforchheimerfluxvariables.hh
index 6fc90ca614fc698fd527181ab6d23290fb7e06ff..184289b2402182ee3041dbf3dd32af58b0d449f2 100644
--- a/dumux/implicit/box/boxforchheimerfluxvariables.hh
+++ b/dumux/implicit/box/boxforchheimerfluxvariables.hh
@@ -145,9 +145,17 @@ protected:
         }
         else
         {
+            const Element& elementI = *this->fvGeometry_.neighbors[this->face().i];
+            FVElementGeometry fvGeometryI;
+            fvGeometryI.subContVol[0].global = elementI.geometry().center();
+            
+            const Element& elementJ = *this->fvGeometry_.neighbors[this->face().j];
+            FVElementGeometry fvGeometryJ;
+            fvGeometryJ.subContVol[0].global = elementJ.geometry().center();
+            
             spatialParams.meanK(K,
-                                spatialParams.elementIntrinsicPermeability(*this->fvGeometry_.neighbors[this->face().i]),
-                                spatialParams.elementIntrinsicPermeability(*this->fvGeometry_.neighbors[this->face().j]));
+                                spatialParams.intrinsicPermeability(elementI, fvGeometryI, 0),
+                                spatialParams.intrinsicPermeability(elementJ, fvGeometryJ, 0));
         }
         
         // obtain the Forchheimer coefficient from the spatial parameters
diff --git a/dumux/implicit/common/implicitdarcyfluxvariables.hh b/dumux/implicit/common/implicitdarcyfluxvariables.hh
index d529fce5e4bc94b2b4115e52efc3f38666543bfe..9507d03a895b2bdea7db71dbabcd84d9b7189d71 100644
--- a/dumux/implicit/common/implicitdarcyfluxvariables.hh
+++ b/dumux/implicit/common/implicitdarcyfluxvariables.hh
@@ -257,9 +257,17 @@ protected:
         }
         else
         {
+            const Element& elementI = *fvGeometry_.neighbors[face().i];
+            FVElementGeometry fvGeometryI;
+            fvGeometryI.subContVol[0].global = elementI.geometry().center();
+            
+            const Element& elementJ = *fvGeometry_.neighbors[face().j];
+            FVElementGeometry fvGeometryJ;
+            fvGeometryJ.subContVol[0].global = elementJ.geometry().center();
+            
             spatialParams.meanK(K,
-                                spatialParams.elementIntrinsicPermeability(*fvGeometry_.neighbors[face().i]),
-                                spatialParams.elementIntrinsicPermeability(*fvGeometry_.neighbors[face().j]));
+                                spatialParams.intrinsicPermeability(elementI, fvGeometryI, 0),
+                                spatialParams.intrinsicPermeability(elementJ, fvGeometryJ, 0));
         }
         
         // loop over all phases
diff --git a/dumux/material/spatialparams/boxspatialparams1p.hh b/dumux/material/spatialparams/boxspatialparams1p.hh
index b1c6e0e7060ac44a42eb260f29110408694dd288..bb73bcfd7b32ac89a5606949faa3e097a8a058b3 100644
--- a/dumux/material/spatialparams/boxspatialparams1p.hh
+++ b/dumux/material/spatialparams/boxspatialparams1p.hh
@@ -107,13 +107,6 @@ public:
                 result[i][j] = harmonicMean(K1[i][j], K2[i][j]);
     }
 
-    Scalar elementIntrinsicPermeability (const Element &element) const
-    {
-        FVElementGeometry fvGeometry;
-        fvGeometry.subContVol[0].global = element.geometry().center();
-        return asImp_().intrinsicPermeability(element, fvGeometry, /*scvIdx=*/0);
-    }
-    
     /*!
      * \brief Function for defining the intrinsic (absolute) permeability.
      *