diff --git a/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh b/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh index 2afe1fb7ab1ca5f2f709997fd5d4a9d678268aeb..ba4cad3a1dad70e440f6aa8fdd944858f47cbe18 100644 --- a/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh +++ b/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh @@ -250,6 +250,20 @@ public: static Scalar sweToSw_(const Params ¶ms, Scalar swe) { return sweToSw(params, swe); } + /*! + * \brief Convert an effective non-wetting saturation to an absolute one. + * + * \param sne Effective saturation of the non-wetting phase \f$\mathrm{[{S}_n]}\f$. + * \param params A container object that is populated with the appropriate coefficients for the respective law. + * Therefore, in the (problem specific) spatialParameters first, the material law is chosen, + * and then the params container is constructed accordingly. Afterwards the values are set there, too. + * \return Absolute saturation of the non-wetting phase. + */ + static Scalar sneToSn(const Params ¶ms, Scalar sne) + { + return sne*(1. - params.swr() - params.snr()) + params.snr(); + } + /*! * \brief Derivative of the effective saturation w.r.t. the absolute saturation. *