RichardsDensityPrimePrimeAux.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 the 2nd derivative of density wrt pressure
//
#include "RichardsDensityPrimePrimeAux.h"

registerMooseObject("RichardsApp", RichardsDensityPrimePrimeAux);

InputParameters
RichardsDensityPrimePrimeAux::validParams()
{
  InputParameters params = AuxKernel::validParams();
  params.addRequiredCoupledVar("pressure_var", "The variable that represents the pressure");
  params.addRequiredParam<UserObjectName>("density_UO",
                                          "Name of user object that defines the density.");
  params.addClassDescription("auxillary variable which is d^2(density)/dp^2");
  return params;
}

RichardsDensityPrimePrimeAux::RichardsDensityPrimePrimeAux(const InputParameters & parameters)
  : AuxKernel(parameters),
    _pressure_var(coupledValue("pressure_var")),
    _density_UO(getUserObject<RichardsDensity>("density_UO"))
{
}

Real
RichardsDensityPrimePrimeAux::computeValue()
{
  return _density_UO.d2density(_pressure_var[_qp]);
}