16#ifndef DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_EXPONENTIAL_CUBIC 
   17#define DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_EXPONENTIAL_CUBIC 
   20#include <dune/common/exceptions.hh> 
   21#include <dune/common/float_cmp.hh> 
   35    template<
class Scalar>
 
   39        : a1_(
a1), a2_(
a2), a3_(
a3) {}
 
 
   41        Scalar 
a1()
 const { 
return a1_; }
 
   44        Scalar 
a2()
 const { 
return a2_; }
 
   47        Scalar 
a3()
 const { 
return a3_; }
 
   52            return Dune::FloatCmp::eq(
a1(), p.
a1(), 1e-6)
 
   53                   && Dune::FloatCmp::eq(
a2(), p.
a2(), 1e-6)
 
   54                   && Dune::FloatCmp::eq(
a3(), p.
a3(), 1e-6);
 
 
 
   65    template<
class Scalar = 
double>
 
   85    template<
class Scalar>
 
   89        const Scalar a1 = params.
a1();
 
   90        const Scalar a2 = params.
a2();
 
   91        const Scalar a3 = params.
a3();
 
   94        return a1 * exp( a2 * swe) + a3 * pc * pc * pc ;
 
 
 
Implementation of a exponential function relating specific interfacial area to wetting phase saturati...
Definition exponentialcubic.hh:32
static Params< Scalar > makeParams(const std::string ¶mGroup)
Construct from a subgroup from the global parameter tree.
Definition exponentialcubic.hh:66
static Scalar area(const Scalar swe, const Scalar pc, const Params< Scalar > ¶ms)
The interfacial area.
Definition exponentialcubic.hh:86
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 exponentialcubic.hh:37
void setA2(Scalar a2)
Definition exponentialcubic.hh:45
Params(Scalar a1=0, Scalar a2=0, Scalar a3=0)
Definition exponentialcubic.hh:38
Scalar a1() const
Definition exponentialcubic.hh:41
void setA1(Scalar a1)
Definition exponentialcubic.hh:42
bool operator==(const Params &p) const
Definition exponentialcubic.hh:50
Scalar a2() const
Definition exponentialcubic.hh:44
Scalar a3() const
Definition exponentialcubic.hh:47
void setA3(Scalar a3)
Definition exponentialcubic.hh:48