PorousFlowBroadbridgeWhite.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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);
}