SaturationTemperatureFunction.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 "SaturationTemperatureFunction.h"
#include "TwoPhaseFluidProperties.h"

registerMooseObject("FluidPropertiesApp", SaturationTemperatureFunction);

InputParameters
SaturationTemperatureFunction::validParams()
{
  InputParameters params = Function::validParams();

  params.addRequiredParam<FunctionName>("p", "Pressure function");
  params.addRequiredParam<UserObjectName>("fp_2phase", "2-phase fluid properties");

  params.addClassDescription(
      "Computes saturation temperature from pressure function and 2-phase fluid properties object");

  return params;
}

SaturationTemperatureFunction::SaturationTemperatureFunction(const InputParameters & parameters)
  : Function(parameters),
    FunctionInterface(this),

    _p_fn(getFunction("p"))
{
}

void
SaturationTemperatureFunction::initialSetup()
{
  _fp_2phase = &getUserObject<TwoPhaseFluidProperties>("fp_2phase");
}

Real
SaturationTemperatureFunction::value(Real t, const Point & point) const
{
  return _fp_2phase->T_sat(_p_fn.value(t, point));
}

RealVectorValue
SaturationTemperatureFunction::gradient(Real t, const Point & point) const
{
  return _fp_2phase->dT_sat_dp(_p_fn.value(t, point)) * _p_fn.gradient(t, point);
}