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 the 2nd derivative of density wrt pressure
      11             : //
      12             : #include "RichardsDensityPrimePrimeAux.h"
      13             : 
      14             : registerMooseObject("RichardsApp", RichardsDensityPrimePrimeAux);
      15             : 
      16             : InputParameters
      17          25 : RichardsDensityPrimePrimeAux::validParams()
      18             : {
      19          25 :   InputParameters params = AuxKernel::validParams();
      20          50 :   params.addRequiredCoupledVar("pressure_var", "The variable that represents the pressure");
      21          50 :   params.addRequiredParam<UserObjectName>("density_UO",
      22             :                                           "Name of user object that defines the density.");
      23          25 :   params.addClassDescription("auxillary variable which is d^2(density)/dp^2");
      24          25 :   return params;
      25           0 : }
      26             : 
      27          10 : RichardsDensityPrimePrimeAux::RichardsDensityPrimePrimeAux(const InputParameters & parameters)
      28             :   : AuxKernel(parameters),
      29          10 :     _pressure_var(coupledValue("pressure_var")),
      30          20 :     _density_UO(getUserObject<RichardsDensity>("density_UO"))
      31             : {
      32          10 : }
      33             : 
      34             : Real
      35        3030 : RichardsDensityPrimePrimeAux::computeValue()
      36             : {
      37        3030 :   return _density_UO.d2density(_pressure_var[_qp]);
      38             : }