16#ifndef DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_PC_MAX 
   17#define DUMUX_MATERIAL_FLUIDMATRIX_TWO_P_INTERFACIAL_AREA_PC_MAX 
   20#include <dune/common/exceptions.hh> 
   21#include <dune/common/float_cmp.hh> 
   35    template<
class Scalar>
 
   41        Scalar 
pcMax()
 const { 
return pcMax_; }
 
   44        Scalar 
a1()
 const { 
return a1_; }
 
   47        Scalar 
a2()
 const { 
return a2_; }
 
   50        Scalar 
a3()
 const { 
return a3_; }
 
   55            return Dune::FloatCmp::eq(
pcMax(), p.
pcMax(), 1e-6)
 
   56                   && Dune::FloatCmp::eq(
a1(), p.
a1(), 1e-6)
 
   57                   && Dune::FloatCmp::eq(
a2(), p.
a2(), 1e-6)
 
   58                   && Dune::FloatCmp::eq(
a3(), p.
a3(), 1e-6);
 
 
   62        Scalar pcMax_, a1_, a2_, a3_;
 
 
   69    template<
class Scalar = 
double>
 
   76        return {pcMax, a1, a2, a3};
 
 
   90    template<
class Scalar>
 
   93      const Scalar a1 = params.
a1();
 
   94      const Scalar a2 = params.
a2();
 
   95      const Scalar a3 = params.
a3();
 
   96      const Scalar pcMax = params.
pcMax();
 
   97      return a1 * (pcMax-pc) * (1.0-swe) +  a2*(pcMax-pc)*(pcMax-pc) * (1.-swe) +  a3 * (pcMax-pc)*(1-swe)*(1-swe);
 
 
 
Implementation of the polynomial of second order relating specific interfacial area to wetting phase ...
Definition pcmax.hh:32
static Scalar area(const Scalar swe, const Scalar pc, const Params< Scalar > ¶ms)
The interfacial area.
Definition pcmax.hh:91
static Params< Scalar > makeParams(const std::string ¶mGroup)
Construct from a subgroup from the global parameter tree.
Definition pcmax.hh:70
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.
Scalar a2() const
Definition pcmax.hh:47
void setA1(Scalar a1)
Definition pcmax.hh:45
Params(Scalar pcMax=0, Scalar a1=0, Scalar a2=0, Scalar a3=0)
Definition pcmax.hh:38
void setA3(Scalar a3)
Definition pcmax.hh:51
Scalar a3() const
Definition pcmax.hh:50
bool operator==(const Params &p) const
Definition pcmax.hh:53
Scalar a1() const
Definition pcmax.hh:44
void setA2(Scalar a2)
Definition pcmax.hh:48
void setPcMax(Scalar pcMax)
Definition pcmax.hh:42
Scalar pcMax() const
Definition pcmax.hh:41