From cccc36c5abcc714cb234dad81f8a288c67709537 Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Tue, 31 Mar 2020 10:52:59 +0200
Subject: [PATCH] [efftoabs] Add sneToSn for symmetry reasons

---
 .../fluidmatrixinteractions/2p/efftoabslaw.hh      | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh b/dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh
index 2afe1fb7ab..ba4cad3a1d 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 &params, 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 &params, Scalar sne)
+    {
+        return sne*(1. - params.swr() - params.snr()) + params.snr();
+    }
+
     /*!
      * \brief Derivative of the effective saturation w.r.t. the absolute saturation.
      *
-- 
GitLab