MultiComponentFluidPropertiesMaterialPT.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "MultiComponentFluidPropertiesMaterialPT.h"
 
registerMooseObject("FluidPropertiesTestApp", MultiComponentFluidPropertiesMaterialPT);
 
template <>
InputParameters
validParams<MultiComponentFluidPropertiesMaterialPT>()
{
  InputParameters params = validParams<Material>();
  params.addRequiredCoupledVar("pressure", "pressure (Pa)");
  params.addRequiredCoupledVar("temperature", "temperature (K)");
  params.addRequiredCoupledVar("xmass", "Solute mass fraction (-)");
  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");
  params.addClassDescription(
      "Fluid properties of a multicomponent fluid using the (pressure, temperature) formulation");
  return params;
}
 
MultiComponentFluidPropertiesMaterialPT::MultiComponentFluidPropertiesMaterialPT(
    const InputParameters & parameters)
  : Material(parameters),
    _pressure(coupledValue("pressure")),
    _temperature(coupledValue("temperature")),
    _xmass(coupledValue("xmass")),
 
    _rho(declareProperty<Real>("density")),
    _h(declareProperty<Real>("enthalpy")),
    _cp(declareProperty<Real>("cp")),
    _e(declareProperty<Real>("e")),
 
    _fp(getUserObject<MultiComponentFluidProperties>("fp"))
{
}
 
MultiComponentFluidPropertiesMaterialPT::~MultiComponentFluidPropertiesMaterialPT() {}
 
void
MultiComponentFluidPropertiesMaterialPT::computeQpProperties()
{
  _rho[_qp] = _fp.rho_from_p_T_X(_pressure[_qp], _temperature[_qp], _xmass[_qp]);
  _h[_qp] = _fp.h_from_p_T_X(_pressure[_qp], _temperature[_qp], _xmass[_qp]);
  _cp[_qp] = _fp.cp_from_p_T_X(_pressure[_qp], _temperature[_qp], _xmass[_qp]);
  _e[_qp] = _fp.e_from_p_T_X(_pressure[_qp], _temperature[_qp], _xmass[_qp]);
}