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 : public RichardsSeff1VG
{
public:
  static InputParameters validParams();

  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;
};