PorousFlowBroadbridgeWhite.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "MooseTypes.h"
#include "MooseError.h"
#include "libmesh/utility.h"
 
namespace PorousFlowBroadbridgeWhite
{
Real LambertW(Real z);
 
Real effectiveSaturation(Real pc, Real c, Real sn, Real ss, Real las);
 
Real dEffectiveSaturation(Real pc, Real c, Real sn, Real ss, Real las);
 
Real d2EffectiveSaturation(Real pc, Real c, Real sn, Real ss, Real las);
 
template <typename T>
T
relativePermeability(const T & s, Real c, Real sn, Real ss, Real kn, Real ks)
{
  if (MetaPhysicL::raw_value(s) <= sn)
    return kn;
 
  if (MetaPhysicL::raw_value(s) >= ss)
    return ks;
 
  const T coef = (ks - kn) * (c - 1.0);
  const T th = (s - sn) / (ss - sn);
  const T krel = kn + coef * Utility::pow<2>(th) / (c - th);
  return krel;
}
 
Real dRelativePermeability(Real s, Real c, Real sn, Real ss, Real kn, Real ks);
 
Real d2RelativePermeability(Real s, Real c, Real sn, Real ss, Real kn, Real ks);
}