From 6286b8f6a7a7cacfef87895af6338e946846cfd0 Mon Sep 17 00:00:00 2001
From: Philipp Nuske <philipp.nuske@mailbox.org>
Date: Wed, 9 Apr 2014 16:21:24 +0000
Subject: [PATCH] adding the heatpipelaw from devel to stable. This relation
gives capillary pressure as a function of interfacial tension, permeability
and porosity. It was prominently used in Udell (1985), but is also known as
Leverett Function.
reviewed (like the last commit) by Christoph
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12742 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
.../fluidmatrixinteractions/2p/heatpipelaw.hh | 162 ++++++++++++++++++
.../2p/heatpipelawparams.hh | 79 +++++++++
2 files changed, 241 insertions(+)
create mode 100644 dumux/material/fluidmatrixinteractions/2p/heatpipelaw.hh
create mode 100644 dumux/material/fluidmatrixinteractions/2p/heatpipelawparams.hh
diff --git a/dumux/material/fluidmatrixinteractions/2p/heatpipelaw.hh b/dumux/material/fluidmatrixinteractions/2p/heatpipelaw.hh
new file mode 100644
index 0000000000..5197280472
--- /dev/null
+++ b/dumux/material/fluidmatrixinteractions/2p/heatpipelaw.hh
@@ -0,0 +1,162 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 2 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file heatpipelaw.hh
+ *
+ * Implementation of the capillary pressure <-> saturation relation
+ * for the heatpipe problem.
+ */
+#ifndef HEATPIPELAW_HH
+#define HEATPIPELAW_HH
+
+#include "heatpipelawparams.hh"
+
+#include <dumux/common/spline.hh>
+
+#include <algorithm>
+
+#include <math.h>
+#include <assert.h>
+
+namespace Dumux
+{
+/*!
+ * \ingroup material
+ *
+ * \brief Implementation of the capillary pressure <-> saturation
+ * relation for the heatpipe problem. This class bundles the
+ * "raw" curves as static members and doesn't concern itself
+ * converting absolute to effective saturations and vince
+ * versa.
+ */
+template <class ScalarT, class ParamsT = HeatPipeLawParams<ScalarT> >
+class HeatPipeLaw
+{
+public:
+ typedef ParamsT Params;
+ typedef typename Params::Scalar Scalar;
+
+ /*!
+ * \brief The capillary pressure-saturation curve.
+ *
+ * \param Sw Effective saturation of of the wetting phase \f$\overline{S}_w\f$
+ */
+ static Scalar pc(const Params ¶ms, Scalar Sw)
+ {
+ Scalar Sn = 1 - Sw;
+ Scalar p0Gamma = params.p0()*params.gamma();
+
+ // regularization
+ if (Sn >= 1.0) {
+ Scalar y = p0Gamma*( (1.263*1.0 - 2.120)*1.0 + 1.417)*1.0;
+ Scalar m = p0Gamma*((3*1.263*1.0 - 2*2.120)*1.0 + 1.417);
+ return (Sn - 1)*m + y;
+ }
+ else if (Sn <= 0.0) {
+ Scalar y = 0.0;
+ Scalar m = p0Gamma*1.417;
+ return Sn*m + y;
+ }
+
+ return p0Gamma*((1.263*Sn - 2.120)*Sn + 1.417) * Sn;
+ }
+
+ /*!
+ * \brief The saturation-capillary pressure curve.
+ *
+ * \return The effective saturaion of the wetting phase \f$\overline{S}_w\f$
+ */
+ static Scalar Sw(const Params ¶ms, Scalar pC)
+ {
+ DUNE_THROW(Dune::NotImplemented, "HeatPipeLaw::Sw");
+ }
+
+ /*!
+ * \brief Returns the partial derivative of the capillary
+ * pressure to the effective saturation.
+ */
+ static Scalar dpC_dSw(const Params ¶ms, Scalar Sw)
+ {
+ Scalar Sn = 1 - Sw;
+ Scalar p0Gamma = params.p0()*params.gamma();
+ if (Sn > 1.0)
+ Sn = 1.0;
+ else if (Sn <= 0.0) {
+ Scalar m = -p0Gamma*1.417;
+ return m;
+ }
+
+ Scalar m = - p0Gamma*((3*1.263*Sn - 2*2.120)*Sn + 1.417);
+ return m;
+ }
+
+ /*!
+ * \brief Returns the partial derivative of the effective
+ * saturation to the capillary pressure.
+ */
+ static Scalar dSw_dpC(const Params ¶ms, Scalar pC)
+ {
+ DUNE_THROW(Dune::NotImplemented, "HeatPipeLaw::dSw_dpC");
+ }
+
+ /*!
+ * \brief The relative permeability for the wetting phase.
+ *
+ * \param Sw The mobile saturation of the wetting phase.
+ */
+ static Scalar krw(const Params ¶ms, Scalar Sw)
+ {
+ return kr_(Sw);
+ };
+
+ /*!
+ * \brief The relative permeability for the non-wetting phase.
+ *
+ * \param Sw The mobile saturation of the wetting phase.
+ */
+ static Scalar krn(const Params ¶ms, Scalar Sw)
+ {
+ Scalar Sn = 1 - Sw;
+ return kr_(Sn);
+ }
+
+private:
+ static Scalar kr_(Scalar S)
+ {
+ const Scalar eps = 0.95;
+ if (S >= 1)
+ return 1;
+ else if (S <= 0)
+ return 0;
+ else if (S > eps) {
+ // regularize
+ typedef Dumux::Spline<Scalar> Spline;
+ Spline sp(eps, 1.0, // x1, x2
+ eps*eps*eps, 1, // y1, y2
+ 3*eps*eps, 0); // m1, m2
+ return sp.eval(S);
+ }
+
+ return S*S*S;
+ }
+};
+
+}
+
+#endif
diff --git a/dumux/material/fluidmatrixinteractions/2p/heatpipelawparams.hh b/dumux/material/fluidmatrixinteractions/2p/heatpipelawparams.hh
new file mode 100644
index 0000000000..0a2eef559f
--- /dev/null
+++ b/dumux/material/fluidmatrixinteractions/2p/heatpipelawparams.hh
@@ -0,0 +1,79 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 2 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file heatpipelawparams.hh
+ *
+ * Specification of the material params for the heat pipe's capillary
+ * pressure model.
+ */
+#ifndef DUMUX_HEATPIPELAWPARAMS_HH
+#define DUMUX_HEATPIPELAWPARAMS_HH
+
+namespace Dumux
+{
+/*!
+ * \brief Reference implementation of a params for the heat pipe's
+ * material law
+ */
+template<class ScalarT>
+class HeatPipeLawParams
+{
+public:
+ typedef ScalarT Scalar;
+
+ HeatPipeLawParams()
+ {}
+
+ HeatPipeLawParams(Scalar p0, Scalar gamma)
+ {
+ setP0(p0);
+ setGamma(gamma);
+ };
+
+ /*!
+ * \brief Return the \f$\gamma\f$ shape parameter.
+ */
+ Scalar gamma() const
+ { return gamma_; }
+
+ /*!
+ * \brief Set the \f$\gamma\f$ shape parameter.
+ */
+ void setGamma(Scalar v)
+ { gamma_ = v; }
+
+ /*!
+ * \brief Return the entry pressure \f$p_0\f$.
+ */
+ Scalar p0() const
+ { return p0_; }
+
+ /*!
+ * \brief Return the entry pressure \f$p_0\f$.
+ */
+ void setP0(Scalar v)
+ { p0_ = v; }
+
+private:
+ Scalar gamma_;
+ Scalar p0_;
+};
+} // namespace Dumux
+
+#endif
--
GitLab