RichardsRelPermPrimePrimeAux.C

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
 
//  This post processor returns d^2(relperm)/d(Seff)^2
//
#include "RichardsRelPermPrimePrimeAux.h"
 
registerMooseObject("RichardsApp", RichardsRelPermPrimePrimeAux);
 
template <>
InputParameters
validParams<RichardsRelPermPrimePrimeAux>()
{
  InputParameters params = validParams<AuxKernel>();
  params.addRequiredCoupledVar("seff_var", "The variable that represents the effective saturation");
  params.addRequiredParam<UserObjectName>(
      "relperm_UO", "Name of user object that defines the relative permeability.");
  params.addClassDescription("auxillary variable which is d^2(relative permeability)/dSeff^2");
  return params;
}
 
RichardsRelPermPrimePrimeAux::RichardsRelPermPrimePrimeAux(const InputParameters & parameters)
  : AuxKernel(parameters),
    _seff_var(coupledValue("seff_var")),
    _relperm_UO(getUserObject<RichardsRelPerm>("relperm_UO"))
{
}
 
Real
RichardsRelPermPrimePrimeAux::computeValue()
{
  return _relperm_UO.d2relperm(_seff_var[_qp]);
}