diff --git a/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh
index f10982966c054fc2b4f2eb82b2f2bd3d6531b429..4d1516a76a6769c1f6d6b405bb645cba289d1040 100644
--- a/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh
+++ b/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh
@@ -203,6 +203,12 @@ public:
     {
         const Scalar Sthres = params.thresholdSw();
 
+        //instead of return value = inf, return a very large number
+        if (params.pe() == 0.0)
+        {
+            return 1e100;
+        }
+
         // calculate the saturation which corresponds to the
         // saturation in the non-regularized version of the
         // Brooks-Corey law