diff --git a/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh index 7e046bef6a6516baab074a17c565bb5a2a7ed871..59e616cde51470c69ae55cdc2e6518255977b69b 100644 --- a/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh +++ b/dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh @@ -178,12 +178,12 @@ public: { if (Sw <= 0) return 0; - else if (Sw > 0) + else if (Sw >= 1) return 1.0; - else if (Sw >= 1 - 0.01) { - Scalar m1 = BrooksCorey::dkrw_dSw(params, 1.0 - 0.01); - Scalar y1 = BrooksCorey::krw(params, 1.0 - 0.01); - Dumux::Spline<Scalar> sp(1 - 0.01, 1.0, + else if (Sw >= 1 - 0.05) { + Scalar m1 = BrooksCorey::dkrw_dSw(params, 1.0 - 0.05); + Scalar y1 = BrooksCorey::krw(params, 1.0 - 0.05); + Dumux::Spline<Scalar> sp(1 - 0.05, 1.0, y1, 1.0, m1, 0); return sp.eval(Sw); @@ -206,10 +206,10 @@ public: // check if we need to regularize the relative permeability else if (Sw <= 0) return 1.0; - else if (Sw < 0.01) { - Scalar m1 = BrooksCorey::dkrn_dSw(params, 0.01); - Scalar y1 = BrooksCorey::krn(params, 0.01); - Dumux::Spline<Scalar> sp(0.0, 0.01, + else if (Sw < 0.05) { + Scalar m1 = BrooksCorey::dkrn_dSw(params, 0.05); + Scalar y1 = BrooksCorey::krn(params, 0.05); + Dumux::Spline<Scalar> sp(0.0, 0.05, 1.0, y1, 0, m1); return sp.eval(Sw);