RichardsDensityIdeal.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

//  Fluid density ideal gas
//
#include "RichardsDensityIdeal.h"

registerMooseObject("RichardsApp", RichardsDensityIdeal);

InputParameters
RichardsDensityIdeal::validParams()
{
  InputParameters params = RichardsDensity::validParams();
  params.addRequiredParam<Real>("slope", "Density = slope*(p - p0)");
  params.addRequiredParam<Real>("p0", "Density = slope*(p - p0)");
  params.addClassDescription("Fluid density of ideal gas.  Density = slope*(p - p0)");
  return params;
}

RichardsDensityIdeal::RichardsDensityIdeal(const InputParameters & parameters)
  : RichardsDensity(parameters), _slope(getParam<Real>("slope")), _p0(getParam<Real>("p0"))
{
}

Real
RichardsDensityIdeal::density(Real p) const
{
  return _slope * (p - _p0);
}

Real RichardsDensityIdeal::ddensity(Real /*p*/) const { return _slope; }

Real RichardsDensityIdeal::d2density(Real /*p*/) const { return 0.0; }