RichardsSeff1VGcut.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
//  "cut" van-Genuchten effective saturation as a function of single pressure, and its derivs wrt to
//  that pressure
//
 
#include "RichardsSeff1VG.h"
 
class RichardsSeff1VGcut;
 
template <>
InputParameters validParams<RichardsSeff1VGcut>();
 
class RichardsSeff1VGcut : public RichardsSeff1VG
{
public:
  RichardsSeff1VGcut(const InputParameters & parameters);
 
  void initialSetup();
 
  Real seff(std::vector<const VariableValue *> p, unsigned int qp) const;
 
  void
  dseff(std::vector<const VariableValue *> p, unsigned int qp, std::vector<Real> & result) const;
 
  void d2seff(std::vector<const VariableValue *> p,
              unsigned int qp,
              std::vector<std::vector<Real>> & result) const;
 
protected:
  Real _al;
 
  Real _m;
 
  Real _p_cut;
 
  Real _s_cut;
 
  Real _ds_cut;
};