HydrostaticPressureFunction.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
 
#include "HydrostaticPressureFunction.h"
#include "SinglePhaseFluidProperties.h"
 
registerMooseObject("FluidPropertiesApp", HydrostaticPressureFunction);
 
InputParameters
HydrostaticPressureFunction::validParams()
{
  InputParameters params = Function::validParams();
  params += GravityVectorInterface::validParams();
 
  params.addClassDescription("Computes hydrostatic pressure from a reference point.");
 
  params.addRequiredParam<Real>("reference_pressure", "Pressure at the reference point [Pa]");
  params.addRequiredParam<Real>("reference_temperature",
                                "Reference temperature for density evaluation [K]");
  params.addRequiredParam<Point>("reference_point",
                                 "Reference point where pressure is specified [m]");
 
  params.addRequiredParam<UserObjectName>("fluid_properties",
                                          "Name of the SinglePhaseFluidProperties user object");
 
  return params;
}
 
HydrostaticPressureFunction::HydrostaticPressureFunction(const InputParameters & parameters)
  : Function(parameters),
    GravityVectorInterface(this),
    _p_ref(getParam<Real>("reference_pressure")),
    _T_ref(getParam<Real>("reference_temperature")),
    _r_ref(getParam<Point>("reference_point"))
{
}
 
void
HydrostaticPressureFunction::initialSetup()
{
  Function::initialSetup();
 
  const auto & fp = getUserObject<SinglePhaseFluidProperties>("fluid_properties");
  _rho = fp.rho_from_p_T(_p_ref, _T_ref);
}
 
Real
HydrostaticPressureFunction::value(Real /*t*/, const Point & r) const
{
  return _p_ref - _rho * gravityVector() * (_r_ref - r);
}