Fluid Properties Module

The Fluid Properties module provides a consistent interface to fluid properties such as density, viscosity, enthalpy and many others, as well as derivatives with respect to the primary variables. The consistent interface allows different fluids to be used in an input file by simply swapping the name of the Fluid Properties UserObject in a plug-and-play manner.

Fluids available

This module provides fluid properties for many different gases, liquids and mixtures.

Available Objects and Subsystems

Fluid Properties App
BrineFluidPropertiesFluid properties for brine
CO2FluidPropertiesFluid properties for carbon dioxide (CO2) using the Span & Wagner EOS
CaloricallyImperfectGasFluid properties for an ideal gas with imperfect caloric behavior.
FlibeFluidPropertiesFluid properties for flibe
FlinakFluidPropertiesFluid properties for flinak
HeliumFluidPropertiesFluid properties for helium
HydrogenFluidPropertiesFluid properties for Hydrogen (H2)
IdealGasFluidPropertiesFluid properties for an ideal gas
IdealRealGasMixtureFluidPropertiesClass for fluid properties of an arbitrary vapor mixture
MethaneFluidPropertiesFluid properties for methane (CH4)
NaClFluidPropertiesFluid properties for NaCl
NitrogenFluidPropertiesFluid properties for Nitrogen (N2)
SimpleFluidPropertiesFluid properties for a simple fluid with a constant bulk density
SodiumPropertiesFluid properties for sodium
SodiumSaturationFluidPropertiesFluid properties for liquid sodium at saturation conditions
StiffenedGasFluidPropertiesFluid properties for a stiffened gas
StiffenedGasTwoPhaseFluidPropertiesTwo-phase stiffened gas fluid properties
TabulatedFluidPropertiesFluid properties using bicubic interpolation on tabulated values provided
TwoPhaseFluidPropertiesIndependent2-phase fluid properties for 2 independent single-phase fluid properties
Water97FluidPropertiesFluid properties for water and steam (H2O) using IAPWS-IF97

Properties available

Each FluidProperties UserObject above provides access to several fluid properties. Two main formulations are available using

(v, e) - specific volume and specific internal energy
(p, T) - pressure and temperature

although some fluid properties may be available in terms of other primary variables. Depending on the actual fluid implementation, any number of the following fluid properties may be available.

Density: rho_from_p_s(pressure, entropy)
Density: rho_from_p_T(pressure, temperature)
Dynamic viscosity: mu_from_v_e(volume, energy)
Dynamic viscosity: mu_from_p_T(pressure, temperature)
Enthalpy: h_from_p_T(pressure, temperature)
Enthalpy: h_from_T_v(temperature, volume)
Fluid name: fluidName()
Gibbs free energy: g_from_v_e(volume, energy)
Internal energy: e_from_v_h(volume, enthalpy)
Internal energy: e_from_p_rho(pressure, density)
Internal energy: e_from_T_v(temperature, volume)
Internal energy: e_from_p_T(pressure, temperature)
Isobaric specific heat: cp_from_v_e(volume, energy)
Isobaric specific heat: cp_from_p_T(pressure, temperature)
Isochoric specific heat: cv_from_v_e(volume, energy)
Isochoric specific heat: cv_from_T_v(temperature, volume)
Isochoric specific heat: cv_from_p_T(pressure, temperature)
Molar mass (kg/mol): molarMass()
Pressure: p_from_v_e(volume, energy)
Pressure: p_from_T_v(temperature, volume)
Pressure: p_from_h_s(pressure, entropy)
Ratio of specific heats: gamma(volume, energy)
Ratio of specific heats: gamma_from_p_T(pressure, temperature)
Specific entropy: s_from_v_e(volume, energy)
Specific entropy: s_from_p_T(pressure, temperature)
Specific entropy: s_from_h_p(enthalpy, pressure)
Specific entropy: s_from_T_v(temperature, volume)
Specific volume: v_from_p_T(pressure, temperature)
Speed of sound: c_from_v_e(volume, energy)
Speed of sound: c_from_p_T(pressure, temperature)
Temperature: T_from_v_e(volume, energy)
Temperature: T_from_p_h(pressure, enthalpy)
Thermal conductivity: k_from_v_e(volume, energy)
Thermal conductivity: k_from_p_T(pressure, temperature)
Thermal expansion coefficient: beta_from_p_T(pressure, temperature)

Derivatives of fluid properties with respect to the primary variables are also available for several of the fluid properties listed above. These can be evaluated using the following notation: rho_from_p_T(p, T, rho, drho_dp, drho_dT) etc.

note

Fluid properties are now available using an interface suitable for use with MOOSE's Automatic Differentiation capability. See example in the next section.

The full list of available methods can be found in either the source code or the Modules Doxygen page for each FluidProperties class.

note

Fluid properties are implemented in GeneralUserObjects that have empty initialize(), execute() and finalize() methods, so do nothing during a simulation. Their purpose is to provide convenient access to fluid properties through the existing UserObject interface.

Usage

All Fluid Properties UserObjects can be accessed in MOOSE objects through the usual UserObject interface. The following example provides a detailed explanation of the steps involved to use the Fluid Properties UserObjects in other MOOSE objects, and the syntax required in the input file.

This example is for a problem that has energy-volume as the primary variables. A material is provided to calculate fluid properties at the quadrature points.

Source

To access the fluid properties defined in the Fluid Properties module in a MOOSE object, the source code of the object must include the following lines of code.

In the header file of the material, a const reference to the base SinglePhaseFluidProperties object is required:

  const SinglePhaseFluidProperties & _fp;

note

A forward declaration to the SinglePhaseFluidProperties class is required at the beginning of the header file.

class SinglePhaseFluidProperties;

In the source file, the SinglePhaseFluidProperties class must be included

#include "SinglePhaseFluidProperties.h"

The Fluid Properties UserObject is passed to this material in the input file by adding a UserObject name parameters in the input parameters:

  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");

The reference to the UserObject is then initialized in the constructor using

    _fp(getUserObject<SinglePhaseFluidProperties>("fp"))

The properties defined in the Fluid Properties UserObject can now be accessed through the reference. In this material, the computeQpProperties method calculates a number of properties at the quadrature points using the values of _v[_qp] and _e[_qp].

FluidPropertiesMaterial::computeQpProperties()
{
  _p[_qp] = _fp.p_from_v_e(_v[_qp], _e[_qp]);
  _T[_qp] = _fp.T_from_v_e(_v[_qp], _e[_qp]);
  _c[_qp] = _fp.c_from_v_e(_v[_qp], _e[_qp]);
  _cp[_qp] = _fp.cp_from_v_e(_v[_qp], _e[_qp]);
  _cv[_qp] = _fp.cv_from_v_e(_v[_qp], _e[_qp]);
  _mu[_qp] = _fp.mu_from_v_e(_v[_qp], _e[_qp]);
  _k[_qp] = _fp.k_from_v_e(_v[_qp], _e[_qp]);
  _g[_qp] = _fp.g_from_v_e(_v[_qp], _e[_qp]);
}

In a similar fashion, fluid properties can be accessed using the Automatic Differentiation interface using the DualReal version which provides both the value and derivatives


DualReal rho = _fp.p_from_T_v(T, v);

where $var element = document.getElementById("moose-equation-3331f02a-28aa-431f-867c-10d3534c62f2");katex.render("T", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ and $var element = document.getElementById("moose-equation-54f0c974-0500-4101-86c4-e4ac17fd5dc3");katex.render("v", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ are DualReal's. The result (density rho in this example) then contains both the value of density and its derivatives with respect to the primary variables T and v.

Input file syntax

The Fluid Properties UserObjects are implemented in an input file in the Modules block. For example, to use the ideal gas formulation for specific volume and energy, the input file syntax would be:

[Modules]
  [./FluidProperties]
    [./ideal_gas]
      type = IdealGasFluidProperties
      gamma = 1.4
      molar_mass = 1.000536678700361
    [../]
  []
[]

In this example, the user has specified a value for gamma (the ratio of isobaric to isochoric specific heat capacites), and R, the universal gas constant.

The fluid properties can then be accessed by other MOOSE objects through the name given in the input file.

[Materials]
  [./fp_mat]
    type = FluidPropertiesMaterial
    e = e
    v = v
    fp = ideal_gas
  [../]
[]

Due to the consistent interface for fluid properties, a different fluid can be substituted in the input file be changing the type of the UserObject. For example, to use a stiffened gas instead of an ideal gas, the only modification required in the input file is

[Modules]
  [./FluidProperties]
    [./sg]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0
      p_inf = 1.e9
      cv = 1816

      mu = 0.9
      k = 0.6
    [../]
  []
[]

Creating additional fluids

New fluids can be added to the Fluid Properties module by inheriting from the base class and overriding the methods that describe the fluid properties. These can then be used in an identical manner as all other Fluid Properties UserObjects.

Utilities

Fluid Properties Interrogator

The FluidPropertiesInterrogator is a user object which can be used to query eligible fluid properties objects.

Additional objects

Several additional objects are provided:

AuxKernels

Available Objects and Subsystems

Fluid Properties App
FluidDensityAuxComputes density from pressure and temperature
PressureAuxComputes pressure given specific volume and specific internal energy
SaturationTemperatureAuxComputes saturation temperature from pressure and 2-phase fluid properties object
SpecificEnthalpyAuxComputes specific enthalpy from pressure and temperature
StagnationPressureAuxComputes stagnation pressure from specific volume, specific internal energy, and velocity
StagnationTemperatureAuxComputes stagnation temperature from specific volume, specific internal energy, and velocity
TemperatureAuxComputes temperature given specific volume and specific internal energy

Functions

Available Objects and Subsystems

Fluid Properties App
SaturationPressureFunctionComputes saturation pressure from temperature function and 2-phase fluid properties object
SaturationTemperatureFunctionComputes saturation temperature from pressure function and 2-phase fluid properties object

Materials

Available Objects and Subsystems

Fluid Properties App
ADSaturationPressureMaterialComputes saturation pressure at some temperature.
ADSaturationTemperatureMaterialComputes saturation temperature at some pressure
ADSurfaceTensionMaterialComputes surface tension at some temperature
FluidPropertiesMaterialComputes fluid properties using (u, v) formulation
FluidPropertiesMaterialPTFluid properties using the (pressure, temperature) formulation
SaturationPressureMaterialComputes saturation pressure at some temperature.
SodiumPropertiesMaterialMaterial properties for liquid sodium sampled from SodiumProperties.

(modules/fluid_properties/include/materials/FluidPropertiesMaterial.h)

// 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

#pragma once

#include "Material.h"

class SinglePhaseFluidProperties;

/**
 * Computes fluid properties using (u, v) formulation
 */
class FluidPropertiesMaterial : public Material
{
public:
  static InputParameters validParams();

  FluidPropertiesMaterial(const InputParameters & parameters);
  virtual ~FluidPropertiesMaterial();

protected:
  virtual void computeQpProperties() override;

  /// Specific internal energy
  const VariableValue & _e;
  /// Specific volume
  const VariableValue & _v;
  /// Pressure
  MaterialProperty<Real> & _p;
  /// Temperature
  MaterialProperty<Real> & _T;
  /// Speed of sound
  MaterialProperty<Real> & _c;
  /// Isobaric specific heat capacity
  MaterialProperty<Real> & _cp;
  /// Isochoric specific heat capacity
  MaterialProperty<Real> & _cv;
  /// Dynamic viscosity
  MaterialProperty<Real> & _mu;
  /// Thermal conductivity
  MaterialProperty<Real> & _k;
  /// Gibbs free energy
  MaterialProperty<Real> & _g;

  /// Fluid properties
  const SinglePhaseFluidProperties & _fp;
};

(modules/fluid_properties/include/materials/FluidPropertiesMaterial.h)

// 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

#pragma once

#include "Material.h"

class SinglePhaseFluidProperties;

/**
 * Computes fluid properties using (u, v) formulation
 */
class FluidPropertiesMaterial : public Material
{
public:
  static InputParameters validParams();

  FluidPropertiesMaterial(const InputParameters & parameters);
  virtual ~FluidPropertiesMaterial();

protected:
  virtual void computeQpProperties() override;

  /// Specific internal energy
  const VariableValue & _e;
  /// Specific volume
  const VariableValue & _v;
  /// Pressure
  MaterialProperty<Real> & _p;
  /// Temperature
  MaterialProperty<Real> & _T;
  /// Speed of sound
  MaterialProperty<Real> & _c;
  /// Isobaric specific heat capacity
  MaterialProperty<Real> & _cp;
  /// Isochoric specific heat capacity
  MaterialProperty<Real> & _cv;
  /// Dynamic viscosity
  MaterialProperty<Real> & _mu;
  /// Thermal conductivity
  MaterialProperty<Real> & _k;
  /// Gibbs free energy
  MaterialProperty<Real> & _g;

  /// Fluid properties
  const SinglePhaseFluidProperties & _fp;
};

(modules/fluid_properties/src/materials/FluidPropertiesMaterial.C)

// 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 "FluidPropertiesMaterial.h"
#include "SinglePhaseFluidProperties.h"

registerMooseObject("FluidPropertiesApp", FluidPropertiesMaterial);

InputParameters
FluidPropertiesMaterial::validParams()
{
  InputParameters params = Material::validParams();
  params.addRequiredCoupledVar("e", "Specific internal energy");
  params.addRequiredCoupledVar("v", "Specific volume");
  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");
  params.addClassDescription("Computes fluid properties using (u, v) formulation");
  return params;
}

FluidPropertiesMaterial::FluidPropertiesMaterial(const InputParameters & parameters)
  : Material(parameters),
    _e(coupledValue("e")),
    _v(coupledValue("v")),

    _p(declareProperty<Real>("pressure")),
    _T(declareProperty<Real>("temperature")),
    _c(declareProperty<Real>("c")),
    _cp(declareProperty<Real>("cp")),
    _cv(declareProperty<Real>("cv")),
    _mu(declareProperty<Real>("mu")),
    _k(declareProperty<Real>("k")),
    _g(declareProperty<Real>("g")),

    _fp(getUserObject<SinglePhaseFluidProperties>("fp"))
{
}

FluidPropertiesMaterial::~FluidPropertiesMaterial() {}

void
FluidPropertiesMaterial::computeQpProperties()
{
  _p[_qp] = _fp.p_from_v_e(_v[_qp], _e[_qp]);
  _T[_qp] = _fp.T_from_v_e(_v[_qp], _e[_qp]);
  _c[_qp] = _fp.c_from_v_e(_v[_qp], _e[_qp]);
  _cp[_qp] = _fp.cp_from_v_e(_v[_qp], _e[_qp]);
  _cv[_qp] = _fp.cv_from_v_e(_v[_qp], _e[_qp]);
  _mu[_qp] = _fp.mu_from_v_e(_v[_qp], _e[_qp]);
  _k[_qp] = _fp.k_from_v_e(_v[_qp], _e[_qp]);
  _g[_qp] = _fp.g_from_v_e(_v[_qp], _e[_qp]);
}

(modules/fluid_properties/src/materials/FluidPropertiesMaterial.C)

// 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 "FluidPropertiesMaterial.h"
#include "SinglePhaseFluidProperties.h"

registerMooseObject("FluidPropertiesApp", FluidPropertiesMaterial);

InputParameters
FluidPropertiesMaterial::validParams()
{
  InputParameters params = Material::validParams();
  params.addRequiredCoupledVar("e", "Specific internal energy");
  params.addRequiredCoupledVar("v", "Specific volume");
  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");
  params.addClassDescription("Computes fluid properties using (u, v) formulation");
  return params;
}

FluidPropertiesMaterial::FluidPropertiesMaterial(const InputParameters & parameters)
  : Material(parameters),
    _e(coupledValue("e")),
    _v(coupledValue("v")),

    _p(declareProperty<Real>("pressure")),
    _T(declareProperty<Real>("temperature")),
    _c(declareProperty<Real>("c")),
    _cp(declareProperty<Real>("cp")),
    _cv(declareProperty<Real>("cv")),
    _mu(declareProperty<Real>("mu")),
    _k(declareProperty<Real>("k")),
    _g(declareProperty<Real>("g")),

    _fp(getUserObject<SinglePhaseFluidProperties>("fp"))
{
}

FluidPropertiesMaterial::~FluidPropertiesMaterial() {}

void
FluidPropertiesMaterial::computeQpProperties()
{
  _p[_qp] = _fp.p_from_v_e(_v[_qp], _e[_qp]);
  _T[_qp] = _fp.T_from_v_e(_v[_qp], _e[_qp]);
  _c[_qp] = _fp.c_from_v_e(_v[_qp], _e[_qp]);
  _cp[_qp] = _fp.cp_from_v_e(_v[_qp], _e[_qp]);
  _cv[_qp] = _fp.cv_from_v_e(_v[_qp], _e[_qp]);
  _mu[_qp] = _fp.mu_from_v_e(_v[_qp], _e[_qp]);
  _k[_qp] = _fp.k_from_v_e(_v[_qp], _e[_qp]);
  _g[_qp] = _fp.g_from_v_e(_v[_qp], _e[_qp]);
}

(modules/fluid_properties/src/materials/FluidPropertiesMaterial.C)

// 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 "FluidPropertiesMaterial.h"
#include "SinglePhaseFluidProperties.h"

registerMooseObject("FluidPropertiesApp", FluidPropertiesMaterial);

InputParameters
FluidPropertiesMaterial::validParams()
{
  InputParameters params = Material::validParams();
  params.addRequiredCoupledVar("e", "Specific internal energy");
  params.addRequiredCoupledVar("v", "Specific volume");
  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");
  params.addClassDescription("Computes fluid properties using (u, v) formulation");
  return params;
}

FluidPropertiesMaterial::FluidPropertiesMaterial(const InputParameters & parameters)
  : Material(parameters),
    _e(coupledValue("e")),
    _v(coupledValue("v")),

    _p(declareProperty<Real>("pressure")),
    _T(declareProperty<Real>("temperature")),
    _c(declareProperty<Real>("c")),
    _cp(declareProperty<Real>("cp")),
    _cv(declareProperty<Real>("cv")),
    _mu(declareProperty<Real>("mu")),
    _k(declareProperty<Real>("k")),
    _g(declareProperty<Real>("g")),

    _fp(getUserObject<SinglePhaseFluidProperties>("fp"))
{
}

FluidPropertiesMaterial::~FluidPropertiesMaterial() {}

void
FluidPropertiesMaterial::computeQpProperties()
{
  _p[_qp] = _fp.p_from_v_e(_v[_qp], _e[_qp]);
  _T[_qp] = _fp.T_from_v_e(_v[_qp], _e[_qp]);
  _c[_qp] = _fp.c_from_v_e(_v[_qp], _e[_qp]);
  _cp[_qp] = _fp.cp_from_v_e(_v[_qp], _e[_qp]);
  _cv[_qp] = _fp.cv_from_v_e(_v[_qp], _e[_qp]);
  _mu[_qp] = _fp.mu_from_v_e(_v[_qp], _e[_qp]);
  _k[_qp] = _fp.k_from_v_e(_v[_qp], _e[_qp]);
  _g[_qp] = _fp.g_from_v_e(_v[_qp], _e[_qp]);
}

(modules/fluid_properties/src/materials/FluidPropertiesMaterial.C)

// 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 "FluidPropertiesMaterial.h"
#include "SinglePhaseFluidProperties.h"

registerMooseObject("FluidPropertiesApp", FluidPropertiesMaterial);

InputParameters
FluidPropertiesMaterial::validParams()
{
  InputParameters params = Material::validParams();
  params.addRequiredCoupledVar("e", "Specific internal energy");
  params.addRequiredCoupledVar("v", "Specific volume");
  params.addRequiredParam<UserObjectName>("fp", "The name of the user object for fluid properties");
  params.addClassDescription("Computes fluid properties using (u, v) formulation");
  return params;
}

FluidPropertiesMaterial::FluidPropertiesMaterial(const InputParameters & parameters)
  : Material(parameters),
    _e(coupledValue("e")),
    _v(coupledValue("v")),

    _p(declareProperty<Real>("pressure")),
    _T(declareProperty<Real>("temperature")),
    _c(declareProperty<Real>("c")),
    _cp(declareProperty<Real>("cp")),
    _cv(declareProperty<Real>("cv")),
    _mu(declareProperty<Real>("mu")),
    _k(declareProperty<Real>("k")),
    _g(declareProperty<Real>("g")),

    _fp(getUserObject<SinglePhaseFluidProperties>("fp"))
{
}

FluidPropertiesMaterial::~FluidPropertiesMaterial() {}

void
FluidPropertiesMaterial::computeQpProperties()
{
  _p[_qp] = _fp.p_from_v_e(_v[_qp], _e[_qp]);
  _T[_qp] = _fp.T_from_v_e(_v[_qp], _e[_qp]);
  _c[_qp] = _fp.c_from_v_e(_v[_qp], _e[_qp]);
  _cp[_qp] = _fp.cp_from_v_e(_v[_qp], _e[_qp]);
  _cv[_qp] = _fp.cv_from_v_e(_v[_qp], _e[_qp]);
  _mu[_qp] = _fp.mu_from_v_e(_v[_qp], _e[_qp]);
  _k[_qp] = _fp.k_from_v_e(_v[_qp], _e[_qp]);
  _g[_qp] = _fp.g_from_v_e(_v[_qp], _e[_qp]);
}

(modules/fluid_properties/test/tests/ideal_gas/test.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 2
  ny = 2
  elem_type = QUAD4
[]

[Functions]
  [./f_fn]
    type = ParsedFunction
    value = -4
  [../]
  [./bc_fn]
    type = ParsedFunction
    value = 'x*x+y*y'
  [../]
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./e]
    initial_condition = 6232.5
  [../]
  [./v]
    initial_condition = 0.02493
  [../]

  [./p]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./T]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cp]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cv]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./c]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./mu]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./k]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./g]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./p]
    type = MaterialRealAux
     variable = p
     property = pressure
  [../]
  [./T]
    type = MaterialRealAux
     variable = T
     property = temperature
  [../]
  [./cp]
    type = MaterialRealAux
     variable = cp
     property = cp
  [../]
  [./cv]
    type = MaterialRealAux
     variable = cv
     property = cv
  [../]
  [./c]
    type = MaterialRealAux
     variable = c
     property = c
  [../]
  [./mu]
    type = MaterialRealAux
     variable = mu
     property = mu
  [../]
  [./k]
    type = MaterialRealAux
     variable = k
     property = k
  [../]
  [./g]
    type = MaterialRealAux
     variable = g
     property = g
  [../]
[]

[Modules]
  [./FluidProperties]
    [./ideal_gas]
      type = IdealGasFluidProperties
      gamma = 1.4
      molar_mass = 1.000536678700361
    [../]
  []
[]

[Materials]
  [./fp_mat]
    type = FluidPropertiesMaterial
    e = e
    v = v
    fp = ideal_gas
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
  [./ffn]
    type = BodyForce
    variable = u
    function = f_fn
  [../]
[]

[BCs]
  [./all]
    type = FunctionDirichletBC
    variable = u
    boundary = 'left right top bottom'
    function = bc_fn
  [../]
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
[]

[Outputs]
  exodus = true
[]

(modules/fluid_properties/test/tests/ideal_gas/test.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 2
  ny = 2
  elem_type = QUAD4
[]

[Functions]
  [./f_fn]
    type = ParsedFunction
    value = -4
  [../]
  [./bc_fn]
    type = ParsedFunction
    value = 'x*x+y*y'
  [../]
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./e]
    initial_condition = 6232.5
  [../]
  [./v]
    initial_condition = 0.02493
  [../]

  [./p]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./T]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cp]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cv]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./c]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./mu]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./k]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./g]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./p]
    type = MaterialRealAux
     variable = p
     property = pressure
  [../]
  [./T]
    type = MaterialRealAux
     variable = T
     property = temperature
  [../]
  [./cp]
    type = MaterialRealAux
     variable = cp
     property = cp
  [../]
  [./cv]
    type = MaterialRealAux
     variable = cv
     property = cv
  [../]
  [./c]
    type = MaterialRealAux
     variable = c
     property = c
  [../]
  [./mu]
    type = MaterialRealAux
     variable = mu
     property = mu
  [../]
  [./k]
    type = MaterialRealAux
     variable = k
     property = k
  [../]
  [./g]
    type = MaterialRealAux
     variable = g
     property = g
  [../]
[]

[Modules]
  [./FluidProperties]
    [./ideal_gas]
      type = IdealGasFluidProperties
      gamma = 1.4
      molar_mass = 1.000536678700361
    [../]
  []
[]

[Materials]
  [./fp_mat]
    type = FluidPropertiesMaterial
    e = e
    v = v
    fp = ideal_gas
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
  [./ffn]
    type = BodyForce
    variable = u
    function = f_fn
  [../]
[]

[BCs]
  [./all]
    type = FunctionDirichletBC
    variable = u
    boundary = 'left right top bottom'
    function = bc_fn
  [../]
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
[]

[Outputs]
  exodus = true
[]

(modules/fluid_properties/test/tests/stiffened_gas/test.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 2
  ny = 2
  elem_type = QUAD4
[]

[Functions]
  [./f_fn]
    type = ParsedFunction
    value = -4
  [../]
  [./bc_fn]
    type = ParsedFunction
    value = 'x*x+y*y'
  [../]
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./e]
    initial_condition = 113206.45935406466
  [../]
  [./v]
    initial_condition = 0.0007354064593540647
  [../]

  [./p]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./T]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cp]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cv]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./c]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./mu]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./k]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./g]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./p]
    type = MaterialRealAux
     variable = p
     property = pressure
  [../]
  [./T]
    type = MaterialRealAux
     variable = T
     property = temperature
  [../]
  [./cp]
    type = MaterialRealAux
     variable = cp
     property = cp
  [../]
  [./cv]
    type = MaterialRealAux
     variable = cv
     property = cv
  [../]
  [./c]
    type = MaterialRealAux
     variable = c
     property = c
  [../]
  [./mu]
    type = MaterialRealAux
     variable = mu
     property = mu
  [../]
  [./k]
    type = MaterialRealAux
     variable = k
     property = k
  [../]
  [./g]
    type = MaterialRealAux
     variable = g
     property = g
  [../]
[]

[Modules]
  [./FluidProperties]
    [./sg]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0
      p_inf = 1.e9
      cv = 1816

      mu = 0.9
      k = 0.6
    [../]
  []
[]

[Materials]
  [./fp_mat]
    type = FluidPropertiesMaterial
    e = e
    v = v
    fp = sg
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
  [./ffn]
    type = BodyForce
    variable = u
    function = f_fn
  [../]
[]

[BCs]
  [./all]
    type = FunctionDirichletBC
    variable = u
    boundary = 'left right top bottom'
    function = bc_fn
  [../]
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
[]

[Outputs]
  exodus = true
[]

Fluids available
Properties available
Usage
Creating additional fluids
Utilities
Additional objects

Install MOOSE

New Users

Examples and Tutorials

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools

INL Applications and Remote Access

Fluid Properties Module

Fluids available

Available Objects and Subsystems

Properties available

Usage

Source

Input file syntax

Creating additional fluids

Utilities

Fluid Properties Interrogator

Additional objects

AuxKernels

Available Objects and Subsystems

Functions

Available Objects and Subsystems

Materials

Available Objects and Subsystems