Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : //  This post processor returns d^2(relperm)/d(Seff)^2
      11             : //
      12             : #include "RichardsRelPermPrimePrimeAux.h"
      13             : 
      14             : registerMooseObject("RichardsApp", RichardsRelPermPrimePrimeAux);
      15             : 
      16             : InputParameters
      17          72 : RichardsRelPermPrimePrimeAux::validParams()
      18             : {
      19          72 :   InputParameters params = AuxKernel::validParams();
      20         144 :   params.addRequiredCoupledVar("seff_var", "The variable that represents the effective saturation");
      21         144 :   params.addRequiredParam<UserObjectName>(
      22             :       "relperm_UO", "Name of user object that defines the relative permeability.");
      23          72 :   params.addClassDescription("auxillary variable which is d^2(relative permeability)/dSeff^2");
      24          72 :   return params;
      25           0 : }
      26             : 
      27          18 : RichardsRelPermPrimePrimeAux::RichardsRelPermPrimePrimeAux(const InputParameters & parameters)
      28             :   : AuxKernel(parameters),
      29          18 :     _seff_var(coupledValue("seff_var")),
      30          36 :     _relperm_UO(getUserObject<RichardsRelPerm>("relperm_UO"))
      31             : {
      32          18 : }
      33             : 
      34             : Real
      35        5454 : RichardsRelPermPrimePrimeAux::computeValue()
      36             : {
      37        5454 :   return _relperm_UO.d2relperm(_seff_var[_qp]);
      38             : }