13#ifndef DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_EXPONENTIAL 
   14#define DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_EXPONENTIAL 
   17#include <dune/common/exceptions.hh> 
   18#include <dune/common/float_cmp.hh> 
   33    template<
class Scalar>
 
   39        Scalar 
swr()
 const { 
return swr_; }
 
   42        Scalar 
a1()
 const { 
return a1_; }
 
   45        Scalar 
a2()
 const { 
return a2_; }
 
   48        Scalar 
a3()
 const { 
return a3_; }
 
   53            return Dune::FloatCmp::eq(
swr(), p.
swr(), 1e-6)
 
   54                   && Dune::FloatCmp::eq(
a1(), p.
a1(), 1e-6)
 
   55                   && Dune::FloatCmp::eq(
a2(), p.
a2(), 1e-6)
 
   56                   && Dune::FloatCmp::eq(
a3(), p.
a3(), 1e-6);
 
 
   60        Scalar swr_, a1_, a2_, a3_;
 
 
   67    template<
class Scalar = 
double>
 
   74        return {swr, a1, a2, a3};
 
 
   89    template<
class Scalar>
 
   92        const Scalar a1 = params.
a1();
 
   93        const Scalar a2 = params.
a2();
 
   94        const Scalar a3 = params.
a3();
 
   95        const Scalar swr = params.
swr();
 
   97        return a1 * (swr-swe) * (1-swe) + a2 * (swr-swe) * (1-swe) * exp( a3 * pc) ;
 
 
  106    template<
class Scalar>
 
  109        const Scalar a2 = params.
a2();
 
  110        const Scalar a3 = params.
a3();
 
  111        const Scalar swr = params.
swr();
 
  113        return a2 * a3 * (swr-swe) * (1-swe) * exp(a3*pc);
 
 
  122    template<
class Scalar>
 
  125        Scalar a1 = params.
a1();
 
  126        Scalar a2 = params.
a2();
 
  127        Scalar a3 = params.
a3();
 
  128        Scalar swr = params.
swr();
 
  130        return - a1 *(1+swr-2*swe) - a2 * exp(a3*pc) * (1+swr-2*swe);
 
 
 
Implementation of the exponential function relating specific interfacial area to wetting phase satura...
Definition exponential.hh:30
static Scalar darea_dsw(const Scalar swe, const Scalar pc, const Params< Scalar > ¶ms)
the derivative of specific interfacial area function w.r.t. saturation
Definition exponential.hh:123
static Params< Scalar > makeParams(const std::string ¶mGroup)
Construct from a subgroup from the global parameter tree.
Definition exponential.hh:68
static Scalar darea_dpc(const Scalar swe, const Scalar pc, const Params< Scalar > ¶ms)
the derivative of specific interfacial area function w.r.t. capillary pressure
Definition exponential.hh:107
static Scalar area(const Scalar swe, const Scalar pc, const Params< Scalar > ¶ms)
The interfacial area.
Definition exponential.hh:90
T getParamFromGroup(Args &&... args)
A free function to get a parameter from the parameter tree singleton with a model group.
Definition parameters.hh:149
Definition brookscorey.hh:23
The infrastructure to retrieve run-time parameters from Dune::ParameterTrees.
Definition exponential.hh:35
void setSwr(Scalar swr)
Definition exponential.hh:40
Scalar a3() const
Definition exponential.hh:48
Params(Scalar swr=0, Scalar a1=0, Scalar a2=0, Scalar a3=0)
Definition exponential.hh:36
Scalar a2() const
Definition exponential.hh:45
void setA1(Scalar a1)
Definition exponential.hh:43
void setA2(Scalar a2)
Definition exponential.hh:46
Scalar a1() const
Definition exponential.hh:42
Scalar swr() const
Definition exponential.hh:39
bool operator==(const Params &p) const
Definition exponential.hh:51
void setA3(Scalar a3)
Definition exponential.hh:49