Commit 85dbcbc7 authored by Timo Koch's avatar Timo Koch
Browse files

Merge branch 'feature/brooks-corey-regularized-kr-derivatives' into 'master'

[brookscorey] add regularization for rel. permeability derivatives

See merge request !1927
parents 4a2d3b7f a83f78dd
......@@ -263,6 +263,25 @@ public:
return BrooksCorey::krw(params, swe);
}
/*!
* \brief A regularized version of the derivative of the relative
* permeability for the wetting phase in regard to the wetting
* saturation of the medium implied by the Brooks-Corey parameterization.
*
* \copydetails BrooksCorey::dkrw_dswe()
*/
static Scalar dkrw_dswe(const Params &params, Scalar swe)
{
// derivative of the regularization
// the slope is zero below sw=0.0 and above sw=1.0
if (swe <= 0.0)
return 0.0;
else if (swe >= 1.0)
return 0.0;
return BrooksCorey::dkrw_dswe(params, swe);
}
/*!
* \brief Regularized version of the relative permeability
* for the non-wetting phase of
......@@ -286,6 +305,25 @@ public:
return BrooksCorey::krn(params, swe);
}
/*!
* \brief A regularized version of the derivative of the relative permeability
* for the non-wetting phase in regard to the wetting saturation of
* the medium as implied by the Brooks-Corey parameterization.
*
* \copydetails BrooksCorey::dkrn_dswe()
*/
static Scalar dkrn_dswe(const Params &params, Scalar swe)
{
// derivative of the regularization
// the slope is zero below sw=0.0 and above sw=1.0
if (swe <= 0)
return 0.0;
else if (swe >= 1)
return 0.0;
return BrooksCorey::dkrn_dswe(params, swe);
}
};
} // end namespace Dumux
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment