Water (H2O) fluid properties as a function of pressure (Pa) and temperature (K) from IAPWS-IF97: Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. More...

#include <Water97FluidProperties.h>

Inheritance diagram for Water97FluidProperties:

[legend]

Public Member Functions
	Water97FluidProperties (const InputParameters &parameters)

virtual	~Water97FluidProperties ()

virtual std::string	fluidName () const override
	Fluid name. More...

virtual Real	molarMass () const override
	Molar mass [kg/mol]. More...

virtual Real	criticalPressure () const override
	Critical pressure. More...

virtual Real	criticalTemperature () const override
	Critical temperature. More...

virtual Real	criticalDensity () const override
	Critical density. More...

virtual Real	triplePointPressure () const override
	Triple point pressure. More...

virtual Real	triplePointTemperature () const override
	Triple point temperature. More...

virtual Real	rho_from_p_T (Real pressure, Real temperature) const override

virtual void	rho_from_p_T (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT) const override

virtual Real	e_from_p_T (Real pressure, Real temperature) const override
	Internal energy from pressure and temperature. More...

virtual void	e_from_p_T (Real pressure, Real temperature, Real &e, Real &de_dp, Real &de_dT) const override
	Internal energy and its derivatives from pressure and temperature. More...

virtual Real	c_from_p_T (Real pressure, Real temperature) const override

virtual Real	cp_from_p_T (Real pressure, Real temperature) const override

virtual Real	cv_from_p_T (Real pressure, Real temperature) const override

virtual Real	mu_from_p_T (Real pressure, Real temperature) const override

virtual void	mu_from_p_T (Real pressure, Real temperature, Real &mu, Real &dmu_dp, Real &dmu_dT) const override

virtual Real	mu_from_rho_T (Real density, Real temperature) const override

virtual void	mu_from_rho_T (Real density, Real temperature, Real ddensity_dT, Real &mu, Real &dmu_drho, Real &dmu_dT) const override
	Dynamic viscosity and its derivatives wrt density and temperature TODO: this shouldn't need 3 input args - AD will assume/call the 2-input version. More...

virtual void	rho_mu_from_p_T (Real pressure, Real temperature, Real &rho, Real &mu) const override

virtual void	rho_mu_from_p_T (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT, Real &mu, Real &dmu_dp, Real &dmu_dT) const override
	Density and viscosity and their derivatives wrt pressure and temperature. More...

virtual Real	k_from_p_T (Real pressure, Real temperature) const override

virtual void	k_from_p_T (Real pressure, Real temperature, Real &k, Real &dk_dp, Real &dk_dT) const override

virtual Real	k_from_rho_T (Real density, Real temperature) const override

virtual Real	s_from_p_T (Real pressure, Real temperature) const override

virtual void	s_from_p_T (Real p, Real T, Real &s, Real &ds_dp, Real &ds_dT) const override

virtual Real	h_from_p_T (Real pressure, Real temperature) const override

virtual void	h_from_p_T (Real pressure, Real temperature, Real &h, Real &dh_dp, Real &dh_dT) const override

virtual Real	vaporPressure (Real temperature) const override
	Vapor pressure. More...

virtual void	vaporPressure (Real temperature, Real &psat, Real &dpsat_dT) const override
	Vapor pressure. More...

Real	vaporTemperature (Real pressure) const
	Saturation temperature as a function of pressure. More...

Real	b23p (Real temperature) const
	Auxillary equation for the boundary between regions 2 and 3. More...

Real	b23T (Real pressure) const
	Auxillary equation for the boundary between regions 2 and 3. More...

unsigned int	inRegion (Real pressure, Real temperature) const
	Determines the phase region that the given pressure and temperature values lie in. More...

unsigned int	subregion3 (Real pressure, Real temperature) const
	Provides the correct subregion index for a (P,T) point in region 3. More...

Real	subregionVolume (Real pi, Real theta, Real a, Real b, Real c, Real d, Real e, unsigned int sid) const
	Specific volume in all subregions of region 3 EXCEPT subregion n (13). More...

Real	densityRegion3 (Real pressure, Real temperature) const
	Density function for Region 3 - supercritical water and steam. More...

virtual Real	temperature_from_ph (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. More...

Real	b2bc (Real pressure) const
	Boundary between subregions b and c in region 2. More...

Real	b3ab (Real pressure) const
	Boundary between subregions a and b in region 3. More...

virtual Real	rho (Real p, Real T) const
	Density from pressure and temperature. More...

e e e e s T T T T T rho v v T virtual e Real	s (Real pressure, Real temperature) const

virtual Real	beta_from_p_T (Real, Real) const

virtual void	beta_from_p_T (Real, Real, Real &, Real &, Real &) const

virtual void	rho_dpT (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT) const
	Density and its derivatives from pressure and temperature. More...

virtual Real	v_from_p_T (Real p, Real T) const
	Specific volume from pressure and temperature. More...

virtual void	v_from_p_T (Real p, Real T, Real &v, Real &dv_dp, Real &dv_dT) const
	Specific volume and its derivatives from pressure and temperature. More...

virtual Real	e_spndl_from_v (Real v) const
	Specific internal energy from temperature and specific volume. More...

virtual void	v_e_spndl_from_T (Real T, Real &v, Real &e) const
	Specific internal energy from temperature and specific volume. More...

virtual Real	h (Real p, Real T) const
	Specific enthalpy from pressure and temperature. More...

virtual void	h_dpT (Real pressure, Real temperature, Real &h, Real &dh_dp, Real &dh_dT) const
	Specific enthalpy and its derivatives from pressure and temperature. More...

virtual Real	e (Real pressure, Real temperature) const

virtual void	e_dpT (Real pressure, Real temperature, Real &e, Real &de_dp, Real &de_dT) const

virtual Real	beta (Real pressure, Real temperature) const

virtual Real	T_from_p_h (Real pressure, Real enthalpy) const
	Temperature from pressure and specific enthalpy. More...

virtual Real	criticalInternalEnergy () const
	Critical specific internal energy. More...

virtual void	rho_e_from_p_T (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT, Real &e, Real &de_dp, Real &de_dT) const
	Density and internal energy and their derivatives wrt pressure and temperature. More...

virtual void	rho_e_dpT (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT, Real &e, Real &de_dp, Real &de_dT) const

virtual Real	c (Real pressure, Real temperature) const

virtual Real	gamma_from_p_T (Real pressure, Real temperature) const
	Adiabatic index - ratio of specific heats. More...

virtual Real	mu (Real pressure, Real temperature) const
	Dynamic viscosity. More...

virtual void	mu_dpT (Real pressure, Real temperature, Real &mu, Real &dmu_dp, Real &dmu_dT) const
	Dynamic viscosity and its derivatives wrt pressure and temperature. More...

virtual void	mu_drhoT_from_rho_T (Real density, Real temperature, Real ddensity_dT, Real &mu, Real &dmu_drho, Real &dmu_dT) const

virtual void	rho_mu (Real pressure, Real temperature, Real &rho, Real &mu) const
	Density and viscosity. More...

virtual void	rho_mu_dpT (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT, Real &mu, Real &dmu_dp, Real &dmu_dT) const

virtual Real	k (Real pressure, Real temperature) const
	Thermal conductivity. More...

virtual void	k_dpT (Real pressure, Real temperature, Real &k, Real &dk_dp, Real &dk_dT) const
	Thermal conductivity and its derivatives wrt pressure and temperature. More...

virtual Real	henryConstant (Real temperature) const
	Henry's law constant for dissolution in water. More...

virtual void	henryConstant (Real temperature, Real &Kh, Real &dKh_dT) const
	Henry's law constant for dissolution in water and derivative wrt temperature. More...

virtual void	henryConstant_dT (Real temperature, Real &Kh, Real &dKh_dT) const

virtual void	vaporPressure_dT (Real temperature, Real &psat, Real &dpsat_dT) const

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

Protected Types
enum	subregionEnum { AB, CD, GH, IJ, JK, MN, OP, QU, RX, UV, WX, EF }
	Enum of subregion ids for region 3. More...

Protected Member Functions
Real	gamma1 (Real pi, Real tau) const
	Gibbs free energy in Region 1 - single phase liquid region. More...

Real	dgamma1_dpi (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 1 wrt pi. More...

Real	d2gamma1_dpi2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 1 wrt pi. More...

Real	dgamma1_dtau (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 1 wrt tau. More...

Real	d2gamma1_dtau2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 1 wrt tau. More...

Real	d2gamma1_dpitau (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 1 wrt pi and tau. More...

Real	gamma2 (Real pi, Real tau) const
	Gibbs free energy in Region 2 - superheated steam. More...

Real	dgamma2_dpi (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 2 wrt pi. More...

Real	d2gamma2_dpi2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 2 wrt pi. More...

Real	dgamma2_dtau (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 2 wrt tau. More...

Real	d2gamma2_dtau2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 2 wrt tau. More...

Real	d2gamma2_dpitau (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 2 wrt pi and tau. More...

Real	phi3 (Real delta, Real tau) const
	Helmholtz free energy in Region 3. More...

Real	dphi3_ddelta (Real delta, Real tau) const
	Derivative of Helmholtz free energy in Region 3 wrt delta. More...

Real	d2phi3_ddelta2 (Real delta, Real tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt delta. More...

Real	dphi3_dtau (Real delta, Real tau) const
	Derivative of Helmholtz free energy in Region 3 wrt tau. More...

Real	d2phi3_dtau2 (Real delta, Real tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt tau. More...

Real	d2phi3_ddeltatau (Real delta, Real tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt delta and tau. More...

Real	gamma5 (Real pi, Real tau) const
	Gibbs free energy in Region 5. More...

Real	dgamma5_dpi (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 5 wrt pi. More...

Real	d2gamma5_dpi2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 5 wrt pi. More...

Real	dgamma5_dtau (Real pi, Real tau) const
	Derivative of Gibbs free energy in Region 5 wrt tau. More...

Real	d2gamma5_dtau2 (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 5 wrt tau. More...

Real	d2gamma5_dpitau (Real pi, Real tau) const
	Second derivative of Gibbs free energy in Region 5 wrt pi and tau. More...

Real	tempXY (Real pressure, subregionEnum xy) const
	Boundaries between subregions in region 3. More...

unsigned int	inRegionPH (Real pressure, Real enthalpy) const
	Determines the phase region that the given pressure and enthaply values lie in. More...

unsigned int	subregion2ph (Real pressure, Real enthalpy) const
	Provides the correct subregion index for a (P,h) point in region 2. More...

unsigned int	subregion3ph (Real pressure, Real enthalpy) const
	Provides the correct subregion index for a (P,h) point in region 3. More...

Real	temperature_from_ph1 (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 1 Eq. More...

Real	temperature_from_ph2a (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 2a Eq. More...

Real	temperature_from_ph2b (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 2b Eq. More...

Real	temperature_from_ph2c (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 2c Eq. More...

Real	temperature_from_ph3a (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 3a Eq. More...

Real	temperature_from_ph3b (Real pressure, Real enthalpy) const
	Backwards equation T(p, h) in Region 3b Eq. More...

virtual Real	henryConstantIAPWS (Real temperature, Real A, Real B, Real C) const
	IAPWS formulation of Henry's law constant for dissolution in water From Guidelines on the Henry's constant and vapour liquid distribution constant for gases in H20 and D20 at high temperatures, IAPWS (2004) More...

virtual void	henryConstantIAPWS (Real temperature, Real &Kh, Real &dKh_dT, Real A, Real B, Real C) const

virtual void	henryConstantIAPWS_dT (Real temperature, Real &Kh, Real &dKh_dT, Real A, Real B, Real C) const

Protected Attributes
const Real	_Mh2o
	Water molar mass (kg/mol) More...

const Real	_Rw
	Specific gas constant for H2O (universal gas constant / molar mass of water - J/kg/K) More...

const Real	_p_critical
	Critical pressure (Pa) More...

const Real	_T_critical
	Critical temperature (K) More...

const Real	_rho_critical
	Critical density (kg/m^3) More...

const Real	_p_triple
	Triple point pressure (Pa) More...

const Real	_T_triple
	Triple point temperature (K) More...

const std::array< Real, 34 >	_n1
	Reference constants used in to calculate thermophysical properties of water. More...

const std::array< int, 34 >	_I1

const std::array< int, 34 >	_J1

const std::array< Real, 20 >	_nph1

const std::array< int, 20 >	_Iph1 {{0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6}}

const std::array< int, 20 >	_Jph1

const std::array< Real, 9 >	_n02
	Constants for region 2. More...

const std::array< Real, 43 >	_n2

const std::array< int, 9 >	_J02 {{0, 1, -5, -4, -3, -2, -1, 2, 3}}

const std::array< int, 43 >	_I2

const std::array< int, 43 >	_J2

const std::array< Real, 36 >	_nph2a

const std::array< int, 36 >	_Iph2a

const std::array< int, 36 >	_Jph2a

const std::array< Real, 38 >	_nph2b

const std::array< int, 38 >	_Iph2b

const std::array< int, 38 >	_Jph2b

const std::array< Real, 23 >	_nph2c

const std::array< int, 23 >	_Iph2c

const std::array< int, 23 >	_Jph2c

const std::array< Real, 5 >	_n23
	Constants for the boundary between regions 2 and 3. More...

const std::array< Real, 40 >	_n3
	Constants for region 3. More...

const std::array< int, 40 >	_I3

const std::array< int, 40 >	_J3

const std::array< std::array< Real, 8 >, 26 >	_par3

const std::array< unsigned int, 26 >	_par3N

const std::vector< std::vector< Real > >	_n3s
	Constants for all 26 subregions in region 3. More...

const std::vector< std::vector< int > >	_I3s

const std::vector< std::vector< int > >	_J3s

const std::array< Real, 31 >	_nph3a

const std::array< int, 31 >	_Iph3a

const std::array< int, 31 >	_Jph3a

const std::array< Real, 33 >	_nph3b

const std::array< int, 33 >	_Iph3b

const std::array< int, 33 >	_Jph3b

const std::array< Real, 10 >	_n4
	Constants for region 4 (the saturation curve up to the critical point) More...

const std::array< int, 6 >	_J05 {{0, 1, -3, -2, -1, 2}}
	Constants for region 5. More...

const std::array< Real, 6 >	_n05

const std::array< int, 6 >	_I5 {{1, 1, 1, 2, 2, 3}}

const std::array< int, 6 >	_J5 {{1, 2, 3, 3, 9, 7}}

const std::array< Real, 6 >	_n5

const std::vector< std::vector< int > >	_tempXY_I
	Constnats for the tempXY() method. More...

const std::vector< std::vector< Real > >	_tempXY_n

const std::array< int, 21 >	_mu_I {{0, 1, 2, 3, 0, 1, 2, 3, 5, 0, 1, 2, 3, 4, 0, 1, 0, 3, 4, 3, 5}}
	Constants from Release on the IAPWS Formulation 2008 for the Viscosity of Ordinary Water Substance. More...

const std::array< int, 21 >	_mu_J {{0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 4, 4, 5, 6, 6}}

const std::array< Real, 4 >	_mu_H0 {{1.67752, 2.20462, 0.6366564, -0.241605}}

const std::array< Real, 21 >	_mu_H1

std::array< Real, 4 >	_k_a {{0.0102811, 0.0299621, 0.0156146, -0.00422464}}
	Constants for thermal conductivity. More...

const std::array< Real, 5 >	_T_star {{1386.0, 540.0, _T_critical, 1.0, 1000.0}}
	Temperature scale for each region. More...

const std::array< Real, 5 >	_p_star {{16.53e6, 1.0e6, 1.0e6, 1.0e6, 1.0e6}}
	Pressure scale for each region. More...

const Real	_R
	Universal gas constant (J/mol/K) More...

const Real	_T_c2k
	Conversion of temperature from Celsius to Kelvin. More...

	propfunc (p, v, e) propfunc(T
	Compute a fluid property given for the state defined by two given properties. More...

e	propfunc (c, v, e) propfunc(cp

e e	propfunc (cv, v, e) propfunc(mu

e e e	propfunc (k, v, e) propfunc(s

e e e e	propfunc (s, h, p) propfunc(rho

e e e e s	propfunc (e, v, h) propfunc(s

e e e e s T	propfunc (pp_sat, p, T) propfunc(mu

e e e e s T T	propfunc (k, rho, T) propfunc(c

e e e e s T T T	propfunc (cp, p, T) propfunc(cv

e e e e s T T T T	propfunc (mu, p, T) propfunc(k

e e e e s T T T T T	propfunc (rho, p, T) propfunc(e

e e e e s T T T T T rho	propfunc (e, T, v) propfunc(p

e e e e s T T T T T rho v	propfunc (h, T, v) propfunc(s

e e e e s T T T T T rho v v	propfunc (cv, T, v) propfunc(h

e e e e s T T T T T rho v v T	propfunc (p, h, s) propfunc(g

	v

e	v

e e	v

e e e	v

e e e e s T T T T T rho v v T	v

e e e e	p

e e e e s	p

e e e e s T T	p

e e e e s T T T	p

e e e e s T T T T	p

e e e e s T T T T T	p

e e e e s T T T T T rho v v	p

e e e e s T	rho

e e e e s T T T T T rho	T

e e e e s T T T T T rho v	T

Detailed Description

Water (H2O) fluid properties as a function of pressure (Pa) and temperature (K) from IAPWS-IF97: Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam.

To avoid iteration in Region 3, the backwards equations from: Revised Supplementary Release on Backward Equations for Specific Volume as a Function of Pressure and Temperature v(p,T) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam are implemented.

Water viscosity as a function of density (kg/m^3) and temperature (K) from: Release on the IAPWS Formulation 2008 for the Viscosity of Ordinary Water Substance.

Thermal conductivity as a function of density (kg/m^3) and temperature (K) from: Revised Release on the IAPS Formulation 1985 for the Thermal Conductivity of Ordinary Water Substance

Definition at line 43 of file Water97FluidProperties.h.

Member Enumeration Documentation

◆ subregionEnum

enum Water97FluidProperties::subregionEnum

protected

Enum of subregion ids for region 3.

Enumerator
AB
CD
GH
IJ
JK
MN
OP
QU
RX
UV
WX
EF

Definition at line 482 of file Water97FluidProperties.h.

   {
     AB,
     CD,
     GH,
     IJ,
     JK,
     MN,
     OP,
     QU,
     RX,
     UV,
     WX,
     EF
   };

Constructor & Destructor Documentation

◆ Water97FluidProperties()

Water97FluidProperties::Water97FluidProperties ( const InputParameters & parameters )

Definition at line 26 of file Water97FluidProperties.C.

   : SinglePhaseFluidProperties(parameters),
     _Mh2o(18.015e-3),
     _Rw(461.526),
     _p_critical(22.064e6),
     _T_critical(647.096),
     _rho_critical(322.0),
     _p_triple(611.657),
     _T_triple(273.16)
 {
 }

◆ ~Water97FluidProperties()

Water97FluidProperties::~Water97FluidProperties ( )

virtual

Definition at line 38 of file Water97FluidProperties.C.

38 {}

Member Function Documentation

◆ b23p()

Real Water97FluidProperties::b23p ( Real temperature ) const

Auxillary equation for the boundary between regions 2 and 3.

Eq. (5) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

temperature water temperature (K)

Returns: pressure (Pa) on the boundary between region 2 and 3

Definition at line 862 of file Water97FluidProperties.C.

Referenced by b3ab(), and inRegion().

 {
   // Check whether the input temperature is within the region of validity of this equation.
   // Valid for 623.15 K <= t <= 863.15 K
   if (temperature < 623.15 || temperature > 863.15)
     throw MooseException(name() +
                          ": b23p(): Temperature is outside range of 623.15 K <= T <= 863.15 K");
 
   return (_n23[0] + _n23[1] * temperature + _n23[2] * temperature * temperature) * 1.e6;
 }

◆ b23T()

Real Water97FluidProperties::b23T ( Real pressure ) const

Auxillary equation for the boundary between regions 2 and 3.

Eq. (6) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure water pressure (Pa)

Returns: temperature (K) on the boundary between regions 2 and 3

Definition at line 874 of file Water97FluidProperties.C.

Referenced by inRegionPH().

 {
   // Check whether the input pressure is within the region of validity of this equation.
   // Valid for 16.529 MPa <= p <= 100 MPa
   if (pressure < 16.529e6 || pressure > 100.0e6)
     throw MooseException(name() + ": b23T(): Pressure is outside range 16.529 MPa <= p <= 100 MPa");
 
   return _n23[3] + std::sqrt((pressure / 1.e6 - _n23[4]) / _n23[2]);
 }

◆ b2bc()

Real Water97FluidProperties::b2bc ( Real pressure ) const

Boundary between subregions b and c in region 2.

Equation 21 from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam.

Parameters

pressure water pressure (Pa)

Returns: enthalpy at boundary (J/kg)

Definition at line 1645 of file Water97FluidProperties.C.

Referenced by subregion2ph().

 {
   // Check whether the input pressure is within the region of validity of this equation.
   if (pressure < 6.5467e6 || pressure > 100.0e6)
     throw MooseException(name() +
                          ": b2bc(): Pressure is outside range of 6.5467 MPa <= p <= 100 MPa");
 
   Real pi = pressure / 1.0e6;
 
   return (0.26526571908428e4 + std::sqrt((pi - 0.45257578905948e1) / 0.12809002730136e-3)) * 1.0e3;
 }

◆ b3ab()

Real Water97FluidProperties::b3ab ( Real pressure ) const

Boundary between subregions a and b in region 3.

Equation 1 from Revised Supplementary Release on Backward Equations for the Functions T(p,h), v(p,h) and T(p,s), v(p,s) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure water pressure (Pa)

Returns: enthalpy at boundary (J/kg)

Definition at line 1761 of file Water97FluidProperties.C.

Referenced by subregion3ph().

 {
   // Check whether the input pressure is within the region of validity of this equation.
   if (pressure < b23p(623.15) || pressure > 100.0e6)
     throw MooseException(name() +
                          ": b3ab(): Pressure is outside range of 16.529 MPa <= p <= 100 MPa");
 
   Real pi = pressure / 1.0e6;
   Real eta = 0.201464004206875e4 + 0.374696550136983e1 * pi - 0.219921901054187e-1 * pi * pi +
              0.87513168600995e-4 * pi * pi * pi;
 
   return eta * 1.0e3;
 }

◆ beta()

Real SinglePhaseFluidProperties::beta	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Definition at line 155 of file SinglePhaseFluidProperties.C.

Referenced by SimpleFluidProperties::beta_from_p_T(), and vaporTemperature().

 {
   return beta_from_p_T(p, T);
 }

◆ beta_from_p_T() [1/2]

Real SinglePhaseFluidProperties::beta_from_p_T	(	Real	p,
		Real	T
	)		const

virtualinherited

Reimplemented in SimpleFluidProperties.

Definition at line 86 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::beta().

 {
   // The volumetric thermal expansion coefficient is defined as
   //   1/v dv/dT)_p
   // It is the fractional change rate of volume with respect to temperature change
   // at constant pressure. Here it is coded as
   //   - 1/rho drho/dT)_p
   // using chain rule with v = v(rho)
 
   Real rho, drho_dp, drho_dT;
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
   return -drho_dT / rho;
 }

◆ beta_from_p_T() [2/2]

void SinglePhaseFluidProperties::beta_from_p_T	(	Real	,
		Real	,
		Real &	,
		Real &	,
		Real &
	)		const

virtualinherited

Reimplemented in SimpleFluidProperties.

Definition at line 80 of file SinglePhaseFluidProperties.C.

 {
   mooseError(name(), ": ", __PRETTY_FUNCTION__, " is not implemented.");
 }

◆ c()

Real SinglePhaseFluidProperties::c	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Definition at line 337 of file SinglePhaseFluidProperties.C.

Referenced by IdealGasFluidProperties::c_from_v_e(), StiffenedGasFluidProperties::c_from_v_e(), densityRegion3(), subregionVolume(), and vaporPressure().

 {
   mooseDeprecated(name(), ": c() is deprecated. Use c_from_p_T() instead");
 
   return c_from_p_T(p, T);
 }

◆ c_from_p_T()

Real Water97FluidProperties::c_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 311 of file Water97FluidProperties.C.

 {
   Real speed2, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       speed2 = _Rw * temperature * Utility::pow<2>(dgamma1_dpi(pi, tau)) /
                (Utility::pow<2>(dgamma1_dpi(pi, tau) - tau * d2gamma1_dpitau(pi, tau)) /
                     (tau * tau * d2gamma1_dtau2(pi, tau)) -
                 d2gamma1_dpi2(pi, tau));
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       speed2 = _Rw * temperature * Utility::pow<2>(pi * dgamma2_dpi(pi, tau)) /
                ((-pi * pi * d2gamma2_dpi2(pi, tau)) +
                 Utility::pow<2>(pi * dgamma2_dpi(pi, tau) - tau * pi * d2gamma2_dpitau(pi, tau)) /
                     (tau * tau * d2gamma2_dtau2(pi, tau)));
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       speed2 =
           _Rw * temperature *
           (2.0 * delta * dphi3_ddelta(delta, tau) + delta * delta * d2phi3_ddelta2(delta, tau) -
            Utility::pow<2>(delta * dphi3_ddelta(delta, tau) -
                            delta * tau * d2phi3_ddeltatau(delta, tau)) /
                (tau * tau * d2phi3_dtau2(delta, tau)));
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       speed2 = _Rw * temperature * Utility::pow<2>(pi * dgamma5_dpi(pi, tau)) /
                ((-pi * pi * d2gamma5_dpi2(pi, tau)) +
                 Utility::pow<2>(pi * dgamma5_dpi(pi, tau) - tau * pi * d2gamma5_dpitau(pi, tau)) /
                     (tau * tau * d2gamma5_dtau2(pi, tau)));
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
 
   return std::sqrt(speed2);
 }

◆ cp_from_p_T()

Real Water97FluidProperties::cp_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 369 of file Water97FluidProperties.C.

 {
   Real specific_heat, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       specific_heat = -_Rw * tau * tau * d2gamma1_dtau2(pi, tau);
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       specific_heat = -_Rw * tau * tau * d2gamma2_dtau2(pi, tau);
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       specific_heat =
           _Rw *
           (-tau * tau * d2phi3_dtau2(delta, tau) +
            (delta * dphi3_ddelta(delta, tau) - delta * tau * d2phi3_ddeltatau(delta, tau)) *
                (delta * dphi3_ddelta(delta, tau) - delta * tau * d2phi3_ddeltatau(delta, tau)) /
                (2.0 * delta * dphi3_ddelta(delta, tau) +
                 delta * delta * d2phi3_ddelta2(delta, tau)));
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       specific_heat = -_Rw * tau * tau * d2gamma5_dtau2(pi, tau);
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
   return specific_heat;
 }

◆ criticalDensity()

Real Water97FluidProperties::criticalDensity ( ) const

overridevirtual

Critical density.

Returns: critical density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 65 of file Water97FluidProperties.C.

 {
   return _rho_critical;
 }

◆ criticalInternalEnergy()

Real SinglePhaseFluidProperties::criticalInternalEnergy ( ) const

virtualinherited

Critical specific internal energy.

Returns: specific internal energy (J/kg)

Reimplemented in StiffenedGasFluidProperties.

Definition at line 131 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e().

 {
   mooseError(name(), ": criticalInternalEnergy() is not implemented");
 }

◆ criticalPressure()

Real Water97FluidProperties::criticalPressure ( ) const

overridevirtual

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 53 of file Water97FluidProperties.C.

 {
   return _p_critical;
 }

◆ criticalTemperature()

Real Water97FluidProperties::criticalTemperature ( ) const

overridevirtual

Critical temperature.

Returns: critical temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 59 of file Water97FluidProperties.C.

 {
   return _T_critical;
 }

◆ cv_from_p_T()

Real Water97FluidProperties::cv_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 418 of file Water97FluidProperties.C.

 {
   Real specific_heat, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       specific_heat =
           _Rw * (-tau * tau * d2gamma1_dtau2(pi, tau) +
                  Utility::pow<2>(dgamma1_dpi(pi, tau) - tau * d2gamma1_dpitau(pi, tau)) /
                      d2gamma1_dpi2(pi, tau));
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       specific_heat =
           _Rw * (-tau * tau * d2gamma2_dtau2(pi, tau) +
                  Utility::pow<2>(dgamma2_dpi(pi, tau) - tau * d2gamma2_dpitau(pi, tau)) /
                      d2gamma2_dpi2(pi, tau));
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       specific_heat = _Rw * (-tau * tau * d2phi3_dtau2(delta, tau));
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       specific_heat =
           _Rw * (-tau * tau * d2gamma5_dtau2(pi, tau) +
                  Utility::pow<2>(dgamma5_dpi(pi, tau) - tau * d2gamma5_dpitau(pi, tau)) /
                      d2gamma5_dpi2(pi, tau));
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
   return specific_heat;
 }

◆ d2gamma1_dpi2()

Real Water97FluidProperties::d2gamma1_dpi2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 1 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi (-)

Definition at line 952 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     sum += _n1[i] * _I1[i] * (_I1[i] - 1) * MathUtils::pow(7.1 - pi, _I1[i] - 2) *
            MathUtils::pow(tau - 1.222, _J1[i]);
 
   return sum;
 }

◆ d2gamma1_dpitau()

Real Water97FluidProperties::d2gamma1_dpitau	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 1 wrt pi and tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi and tau (-)

Definition at line 985 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   Real dg = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     dg += -_n1[i] * _I1[i] * _J1[i] * MathUtils::pow(7.1 - pi, _I1[i] - 1) *
           MathUtils::pow(tau - 1.222, _J1[i] - 1);
 
   return dg;
 }

◆ d2gamma1_dtau2()

Real Water97FluidProperties::d2gamma1_dtau2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 1 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt tau (-)

Definition at line 974 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), and h_from_p_T().

 {
   Real dg = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     dg += _n1[i] * _J1[i] * (_J1[i] - 1) * MathUtils::pow(7.1 - pi, _I1[i]) *
           MathUtils::pow(tau - 1.222, _J1[i] - 2);
 
   return dg;
 }

◆ d2gamma2_dpi2()

Real Water97FluidProperties::d2gamma2_dpi2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 2 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi (-)

Definition at line 1028 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = -1.0 / pi / pi;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     dgr += _n2[i] * _I2[i] * (_I2[i] - 1) * MathUtils::pow(pi, _I2[i] - 2) *
            MathUtils::pow(tau - 0.5, _J2[i]);
 
   return dg0 + dgr;
 }

◆ d2gamma2_dpitau()

Real Water97FluidProperties::d2gamma2_dpitau	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 2 wrt pi and tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi and tau (-)

Definition at line 1076 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     dgr += _n2[i] * _I2[i] * _J2[i] * MathUtils::pow(pi, _I2[i] - 1) *
            MathUtils::pow(tau - 0.5, _J2[i] - 1);
 
   return dg0 + dgr;
 }

◆ d2gamma2_dtau2()

Real Water97FluidProperties::d2gamma2_dtau2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 2 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt tau (-)

Definition at line 1059 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), and h_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
   for (std::size_t i = 0; i < _n02.size(); ++i)
     dg0 += _n02[i] * _J02[i] * (_J02[i] - 1) * MathUtils::pow(tau, _J02[i] - 2);
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     dgr += _n2[i] * _J2[i] * (_J2[i] - 1) * MathUtils::pow(pi, _I2[i]) *
            MathUtils::pow(tau - 0.5, _J2[i] - 2);
 
   return dg0 + dgr;
 }

◆ d2gamma5_dpi2()

Real Water97FluidProperties::d2gamma5_dpi2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 5 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi (-)

Definition at line 1186 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = -1.0 / pi / pi;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     dgr += _n5[i] * _I5[i] * (_I5[i] - 1) * MathUtils::pow(pi, _I5[i] - 2) *
            MathUtils::pow(tau, _J5[i]);
 
   return dg0 + dgr;
 }

◆ d2gamma5_dpitau()

Real Water97FluidProperties::d2gamma5_dpitau	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 5 wrt pi and tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt pi and tau (-)

Definition at line 1234 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     dgr +=
         _n5[i] * _I5[i] * _J5[i] * MathUtils::pow(pi, _I5[i] - 1) * MathUtils::pow(tau, _J5[i] - 1);
 
   return dg0 + dgr;
 }

◆ d2gamma5_dtau2()

Real Water97FluidProperties::d2gamma5_dtau2	(	Real	pi,
		Real	tau
	)		const

protected

Second derivative of Gibbs free energy in Region 5 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: second derivative of Gibbs free energy wrt tau (-)

Definition at line 1217 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), and h_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
   for (std::size_t i = 0; i < _n05.size(); ++i)
     dg0 += _n05[i] * _J05[i] * (_J05[i] - 1) * MathUtils::pow(tau, _J05[i] - 2);
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     dgr += _n5[i] * _J5[i] * (_J5[i] - 1) * MathUtils::pow(pi, _I5[i]) *
            MathUtils::pow(tau, _J5[i] - 2);
 
   return dg0 + dgr;
 }

◆ d2phi3_ddelta2()

Real Water97FluidProperties::d2phi3_ddelta2	(	Real	delta,
		Real	tau
	)		const

protected

Second derivative of Helmholtz free energy in Region 3 wrt delta.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: second derivative of Helmholtz free energy wrt delta (-)

Definition at line 1111 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * _I3[i] * (_I3[i] - 1) * MathUtils::pow(delta, _I3[i] - 2) *
            MathUtils::pow(tau, _J3[i]);
 
   return -_n3[0] / delta / delta + sum;
 }

◆ d2phi3_ddeltatau()

Real Water97FluidProperties::d2phi3_ddeltatau	(	Real	delta,
		Real	tau
	)		const

protected

Second derivative of Helmholtz free energy in Region 3 wrt delta and tau.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: second derivative of Helmholtz free energy wrt delta and tau (-)

Definition at line 1143 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * _I3[i] * _J3[i] * MathUtils::pow(delta, _I3[i] - 1) *
            MathUtils::pow(tau, _J3[i] - 1);
 
   return sum;
 }

◆ d2phi3_dtau2()

Real Water97FluidProperties::d2phi3_dtau2	(	Real	delta,
		Real	tau
	)		const

protected

Second derivative of Helmholtz free energy in Region 3 wrt tau.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: second derivative of Helmholtz free energy wrt tau (-)

Definition at line 1132 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), and h_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * _J3[i] * (_J3[i] - 1) * MathUtils::pow(delta, _I3[i]) *
            MathUtils::pow(tau, _J3[i] - 2);
 
   return sum;
 }

◆ densityRegion3()

Real Water97FluidProperties::densityRegion3	(	Real	pressure,
		Real	temperature
	)		const

Density function for Region 3 - supercritical water and steam.

To avoid iteration, use the backwards equations for region 3 from Revised Supplementary Release on Backward Equations for Specific Volume as a Function of Pressure and Temperature v(p,T) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam.

Parameters

pressure	water pressure (Pa)
temperature	water temperature (K)

Returns: density (kg/m^3) in region 3

Definition at line 1512 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), h_from_p_T(), rho_from_p_T(), and s_from_p_T().

 {
   // Region 3 is subdivided into 26 subregions, each with a given backwards equation
   // to directly calculate density from pressure and temperature without the need for
   // expensive iterations. Find the subregion id that the point is in:
   unsigned int sid = subregion3(pressure, temperature);
 
   Real vstar, pi, theta, a, b, c, d, e;
   unsigned int N;
 
   vstar = _par3[sid][0];
   pi = pressure / _par3[sid][1] / 1.0e6;
   theta = temperature / _par3[sid][2];
   a = _par3[sid][3];
   b = _par3[sid][4];
   c = _par3[sid][5];
   d = _par3[sid][6];
   e = _par3[sid][7];
   N = _par3N[sid];
 
   Real sum = 0.0;
   Real volume = 0.0;
 
   // Note that subregion 13 is the only different formulation
   if (sid == 13)
   {
     for (std::size_t i = 0; i < N; ++i)
       sum += _n3s[sid][i] * MathUtils::pow(pi - a, _I3s[sid][i]) *
              MathUtils::pow(theta - b, _J3s[sid][i]);
 
     volume = vstar * std::exp(sum);
   }
   else
     volume = vstar * subregionVolume(pi, theta, a, b, c, d, e, sid);
 
   // Density is the inverse of volume
   return 1.0 / volume;
 }

◆ dgamma1_dpi()

Real Water97FluidProperties::dgamma1_dpi	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 1 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt pi (-)

Definition at line 941 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     sum += -_n1[i] * _I1[i] * MathUtils::pow(7.1 - pi, _I1[i] - 1) *
            MathUtils::pow(tau - 1.222, _J1[i]);
 
   return sum;
 }

◆ dgamma1_dtau()

Real Water97FluidProperties::dgamma1_dtau	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 1 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt tau (-)

Definition at line 963 of file Water97FluidProperties.C.

Referenced by e_from_p_T(), h_from_p_T(), and s_from_p_T().

 {
   Real g = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     g += _n1[i] * _J1[i] * MathUtils::pow(7.1 - pi, _I1[i]) *
          MathUtils::pow(tau - 1.222, _J1[i] - 1);
 
   return g;
 }

◆ dgamma2_dpi()

Real Water97FluidProperties::dgamma2_dpi	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 2 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt pi (-)

Definition at line 1014 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 1.0 / pi;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     dgr += _n2[i] * _I2[i] * MathUtils::pow(pi, _I2[i] - 1) * MathUtils::pow(tau - 0.5, _J2[i]);
 
   return dg0 + dgr;
 }

◆ dgamma2_dtau()

Real Water97FluidProperties::dgamma2_dtau	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 2 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt tau (-)

Definition at line 1043 of file Water97FluidProperties.C.

Referenced by e_from_p_T(), h_from_p_T(), and s_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
   for (std::size_t i = 0; i < _n02.size(); ++i)
     dg0 += _n02[i] * _J02[i] * MathUtils::pow(tau, _J02[i] - 1);
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     dgr += _n2[i] * _J2[i] * MathUtils::pow(pi, _I2[i]) * MathUtils::pow(tau - 0.5, _J2[i] - 1);
 
   return dg0 + dgr;
 }

◆ dgamma5_dpi()

Real Water97FluidProperties::dgamma5_dpi	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 5 wrt pi.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt pi (-)

Definition at line 1172 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cv_from_p_T(), e_from_p_T(), and rho_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 1.0 / pi;
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     dgr += _n5[i] * _I5[i] * MathUtils::pow(pi, _I5[i] - 1) * MathUtils::pow(tau, _J5[i]);
 
   return dg0 + dgr;
 }

◆ dgamma5_dtau()

Real Water97FluidProperties::dgamma5_dtau	(	Real	pi,
		Real	tau
	)		const

protected

Derivative of Gibbs free energy in Region 5 wrt tau.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: derivative of Gibbs free energy wrt tau (-)

Definition at line 1201 of file Water97FluidProperties.C.

Referenced by e_from_p_T(), h_from_p_T(), and s_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real dg0 = 0.0;
   for (std::size_t i = 0; i < _n05.size(); ++i)
     dg0 += _n05[i] * _J05[i] * MathUtils::pow(tau, _J05[i] - 1);
 
   // Residual part of the Gibbs free energy
   Real dgr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     dgr += _n5[i] * _J5[i] * MathUtils::pow(pi, _I5[i]) * MathUtils::pow(tau, _J5[i] - 1);
 
   return dg0 + dgr;
 }

◆ dphi3_ddelta()

Real Water97FluidProperties::dphi3_ddelta	(	Real	delta,
		Real	tau
	)		const

protected

Derivative of Helmholtz free energy in Region 3 wrt delta.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: derivative of Helmholtz free energy wrt delta (-)

Definition at line 1101 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), e_from_p_T(), h_from_p_T(), and rho_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * _I3[i] * MathUtils::pow(delta, _I3[i] - 1) * MathUtils::pow(tau, _J3[i]);
 
   return _n3[0] / delta + sum;
 }

◆ dphi3_dtau()

Real Water97FluidProperties::dphi3_dtau	(	Real	delta,
		Real	tau
	)		const

protected

Derivative of Helmholtz free energy in Region 3 wrt tau.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: derivative of Helmholtz free energy wrt tau (-)

Definition at line 1122 of file Water97FluidProperties.C.

Referenced by e_from_p_T(), h_from_p_T(), and s_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * _J3[i] * MathUtils::pow(delta, _I3[i]) * MathUtils::pow(tau, _J3[i] - 1);
 
   return sum;
 }

◆ e()

Real SinglePhaseFluidProperties::e	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Definition at line 292 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": e() is deprecated. Use e_from_p_T() instead");
 
   return e_from_p_T(p, T);
 }

◆ e_dpT()

void SinglePhaseFluidProperties::e_dpT	(	Real	pressure,
		Real	temperature,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

virtualinherited

Definition at line 300 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": e_dpT() is deprecated. Use e_from_p_T() instead");
 
   e_from_p_T(p, T, e, de_dp, de_dT);
 }

◆ e_from_p_T() [1/2]

Real Water97FluidProperties::e_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Internal energy from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)

Returns: internal energy (J/kg)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 188 of file Water97FluidProperties.C.

Referenced by e_from_p_T().

 {
   Real internal_energy, pi, tau;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       internal_energy =
           _Rw * temperature * (tau * dgamma1_dtau(pi, tau) - pi * dgamma1_dpi(pi, tau));
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       internal_energy =
           _Rw * temperature * (tau * dgamma2_dtau(pi, tau) - pi * dgamma2_dpi(pi, tau));
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       Real delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       internal_energy = _Rw * temperature * tau * dphi3_dtau(delta, tau);
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       internal_energy =
           _Rw * temperature * (tau * dgamma5_dtau(pi, tau) - pi * dgamma5_dpi(pi, tau));
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
   // Output in J/kg
   return internal_energy;
 }

◆ e_from_p_T() [2/2]

void Water97FluidProperties::e_from_p_T	(	Real	p,
		Real	T,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

overridevirtual

Internal energy and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)
[out]	e	internal energy (J/kg)
[out]	de_dp	derivative of internal energy w.r.t. pressure
[out]	de_dT	derivative of internal energy w.r.t. temperature

Reimplemented from SinglePhaseFluidProperties.

Definition at line 235 of file Water97FluidProperties.C.

 {
   Real pi, tau, dinternal_energy_dp, dinternal_energy_dT;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
     {
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       Real dgdp = dgamma1_dpi(pi, tau);
       Real d2gdpt = d2gamma1_dpitau(pi, tau);
       dinternal_energy_dp =
           _Rw * temperature * (tau * d2gdpt - dgdp - pi * d2gamma1_dpi2(pi, tau)) / _p_star[0];
       dinternal_energy_dT =
           _Rw * (pi * tau * d2gdpt - tau * tau * d2gamma1_dtau2(pi, tau) - pi * dgdp);
       break;
     }
 
     case 2:
     {
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       Real dgdp = dgamma2_dpi(pi, tau);
       Real d2gdpt = d2gamma2_dpitau(pi, tau);
       dinternal_energy_dp =
           _Rw * temperature * (tau * d2gdpt - dgdp - pi * d2gamma2_dpi2(pi, tau)) / _p_star[1];
       dinternal_energy_dT =
           _Rw * (pi * tau * d2gdpt - tau * tau * d2gamma2_dtau2(pi, tau) - pi * dgdp);
       break;
     }
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       Real delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       Real dpdd = dphi3_ddelta(delta, tau);
       Real d2pddt = d2phi3_ddeltatau(delta, tau);
       Real d2pdd2 = d2phi3_ddelta2(delta, tau);
       dinternal_energy_dp =
           _T_star[2] * d2pddt / _rho_critical /
           (2.0 * temperature * delta * dpdd + temperature * delta * delta * d2pdd2);
       dinternal_energy_dT =
           -_Rw * (delta * tau * d2pddt * (dpdd - tau * d2pddt) / (2.0 * dpdd + delta * d2pdd2) +
                   tau * tau * d2phi3_dtau2(delta, tau));
       break;
     }
 
     case 5:
     {
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       Real dgdp = dgamma5_dpi(pi, tau);
       Real d2gdpt = d2gamma5_dpitau(pi, tau);
       dinternal_energy_dp =
           _Rw * temperature * (tau * d2gdpt - dgdp - pi * d2gamma5_dpi2(pi, tau)) / _p_star[4];
       dinternal_energy_dT =
           _Rw * (pi * tau * d2gdpt - tau * tau * d2gamma5_dtau2(pi, tau) - pi * dgdp);
       break;
     }
 
     default:
       mooseError(name(), ": inRegion has given an incorrect region");
   }
 
   e = this->e_from_p_T(pressure, temperature);
   de_dp = dinternal_energy_dp;
   de_dT = dinternal_energy_dT;
 }

◆ e_spndl_from_v()

Real SinglePhaseFluidProperties::e_spndl_from_v ( Real v ) const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in StiffenedGasFluidProperties.

Definition at line 452 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e().

 {
   mooseError(name(), ": ", __PRETTY_FUNCTION__, " not implemented.");
 }

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Definition at line 27 of file FluidProperties.h.

27 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 29 of file FluidProperties.h.

29 {}

◆ fluidName()

std::string Water97FluidProperties::fluidName ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name

Reimplemented from SinglePhaseFluidProperties.

Definition at line 41 of file Water97FluidProperties.C.

 {
   return "water";
 }

◆ gamma1()

Real Water97FluidProperties::gamma1	(	Real	pi,
		Real	tau
	)		const

protected

Gibbs free energy in Region 1 - single phase liquid region.

From Eq. (7) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, IAPWS 2007.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: Gibbs free energy (-)

Definition at line 931 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 0; i < _n1.size(); ++i)
     sum += _n1[i] * MathUtils::pow(7.1 - pi, _I1[i]) * MathUtils::pow(tau - 1.222, _J1[i]);
 
   return sum;
 }

◆ gamma2()

Real Water97FluidProperties::gamma2	(	Real	pi,
		Real	tau
	)		const

protected

Gibbs free energy in Region 2 - superheated steam.

From Eq. (15) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, IAPWS 2007.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: Gibbs free energy (-)

Definition at line 996 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real sum0 = 0.0;
   for (std::size_t i = 0; i < _n02.size(); ++i)
     sum0 += _n02[i] * MathUtils::pow(tau, _J02[i]);
 
   Real g0 = std::log(pi) + sum0;
 
   // Residual part of the Gibbs free energy
   Real gr = 0.0;
   for (std::size_t i = 0; i < _n2.size(); ++i)
     gr += _n2[i] * MathUtils::pow(pi, _I2[i]) * MathUtils::pow(tau - 0.5, _J2[i]);
 
   return g0 + gr;
 }

◆ gamma5()

Real Water97FluidProperties::gamma5	(	Real	pi,
		Real	tau
	)		const

protected

Gibbs free energy in Region 5.

From Eq. (32) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, IAPWS 2007.

Parameters

pi	reduced pressure (-)
tau	reduced temperature (-)

Returns: Gibbs free energy (-)

Definition at line 1154 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

 {
   // Ideal gas part of the Gibbs free energy
   Real sum0 = 0.0;
   for (std::size_t i = 0; i < _n05.size(); ++i)
     sum0 += _n05[i] * MathUtils::pow(tau, _J05[i]);
 
   Real g0 = std::log(pi) + sum0;
 
   // Residual part of the Gibbs free energy
   Real gr = 0.0;
   for (std::size_t i = 0; i < _n5.size(); ++i)
     gr += _n5[i] * MathUtils::pow(pi, _I5[i]) * MathUtils::pow(tau, _J5[i]);
 
   return g0 + gr;
 }

◆ gamma_from_p_T()

Real SinglePhaseFluidProperties::gamma_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Adiabatic index - ratio of specific heats.

Parameters

pressure	fluid pressure (Pa)
temperature	fluid temperature (K)

Returns: gamma (-)

Definition at line 149 of file SinglePhaseFluidProperties.C.

 {
   return cp_from_p_T(p, T) / cv_from_p_T(p, T);
 }

◆ h()

Real SinglePhaseFluidProperties::h	(	Real	p,
		Real	T
	)		const

virtualinherited

Specific enthalpy from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)

Returns: h (J/kg)

Definition at line 429 of file SinglePhaseFluidProperties.C.

Referenced by NaClFluidProperties::e_from_p_T(), IdealGasFluidProperties::e_from_v_h(), StiffenedGasFluidProperties::e_from_v_h(), SinglePhaseFluidProperties::h_dpT(), HelmholtzFluidProperties::h_from_p_T(), IdealGasFluidProperties::h_from_p_T(), StiffenedGasFluidProperties::h_from_p_T(), NaClFluidProperties::h_from_p_T(), IdealGasFluidPropertiesPT::h_from_p_T(), SimpleFluidProperties::h_from_p_T(), TabulatedFluidProperties::h_from_p_T(), h_from_p_T(), StiffenedGasFluidProperties::h_from_T_v(), IdealGasFluidProperties::h_from_T_v(), IdealGasFluidProperties::p_from_h_s(), StiffenedGasFluidProperties::p_from_h_s(), StiffenedGasFluidProperties::s_from_h_p(), IdealGasFluidProperties::s_from_h_p(), IdealGasFluidProperties::T_from_p_h(), and SinglePhaseFluidProperties::T_from_p_h().

 {
   mooseDeprecated(name(), ": h() is deprecated. Use h_from_p_T() instead");
 
   return h_from_p_T(p, T);
 }

◆ h_dpT()

void SinglePhaseFluidProperties::h_dpT	(	Real	pressure,
		Real	temperature,
		Real &	h,
		Real &	dh_dp,
		Real &	dh_dT
	)		const

virtualinherited

Specific enthalpy and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)
[out]	h	specific enthalpy (J/kg)
[out]	dh_dp	derivative of specific enthalpy w.r.t. pressure
[out]	dh_dT	derivative of specific enthalpy w.r.t. temperature

Definition at line 437 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": h_dpT() is deprecated. Use h_from_p_T() instead");
 
   h_from_p_T(p, T, h, dh_dp, dh_dT);
 }

◆ h_from_p_T() [1/2]

Real Water97FluidProperties::h_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 683 of file Water97FluidProperties.C.

Referenced by inRegionPH().

 {
   Real enthalpy, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       enthalpy = _Rw * _T_star[0] * dgamma1_dtau(pi, tau);
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       enthalpy = _Rw * _T_star[1] * dgamma2_dtau(pi, tau);
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       enthalpy =
           _Rw * temperature * (tau * dphi3_dtau(delta, tau) + delta * dphi3_ddelta(delta, tau));
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       enthalpy = _Rw * _T_star[4] * dgamma5_dtau(pi, tau);
       break;
 
     default:
       mooseError("Water97FluidProperties::inRegion has given an incorrect region");
   }
   return enthalpy;
 }

◆ h_from_p_T() [2/2]

void Water97FluidProperties::h_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	h,
		Real &	dh_dp,
		Real &	dh_dT
	)		const

overridevirtual

Definition at line 727 of file Water97FluidProperties.C.

 {
   Real enthalpy, pi, tau, delta, denthalpy_dp, denthalpy_dT;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       enthalpy = _Rw * _T_star[0] * dgamma1_dtau(pi, tau);
       denthalpy_dp = _Rw * _T_star[0] * d2gamma1_dpitau(pi, tau) / _p_star[0];
       denthalpy_dT = -_Rw * tau * tau * d2gamma1_dtau2(pi, tau);
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       enthalpy = _Rw * _T_star[1] * dgamma2_dtau(pi, tau);
       denthalpy_dp = _Rw * _T_star[1] * d2gamma2_dpitau(pi, tau) / _p_star[1];
       denthalpy_dT = -_Rw * tau * tau * d2gamma2_dtau2(pi, tau);
       break;
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       Real dpdd = dphi3_ddelta(delta, tau);
       Real d2pddt = d2phi3_ddeltatau(delta, tau);
       Real d2pdd2 = d2phi3_ddelta2(delta, tau);
       enthalpy = _Rw * temperature * (tau * dphi3_dtau(delta, tau) + delta * dpdd);
       denthalpy_dp = (d2pddt + dpdd + delta * d2pdd2) / _rho_critical /
                      (2.0 * delta * dpdd + delta * delta * d2pdd2);
       denthalpy_dT = _Rw * delta * dpdd * (1.0 - tau * d2pddt / dpdd) *
                          (1.0 - tau * d2pddt / dpdd) / (2.0 + delta * d2pdd2 / dpdd) -
                      _Rw * tau * tau * d2phi3_dtau2(delta, tau);
       break;
     }
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       enthalpy = _Rw * _T_star[4] * dgamma5_dtau(pi, tau);
       denthalpy_dp = _Rw * _T_star[4] * d2gamma5_dpitau(pi, tau) / _p_star[4];
       denthalpy_dT = -_Rw * tau * tau * d2gamma5_dtau2(pi, tau);
       break;
 
     default:
       mooseError("Water97FluidProperties::inRegion has given an incorrect region");
   }
   h = enthalpy;
   dh_dp = denthalpy_dp;
   dh_dT = denthalpy_dT;
 }

◆ henryConstant() [1/2]

Real SinglePhaseFluidProperties::henryConstant ( Real temperature ) const

virtualinherited

Henry's law constant for dissolution in water.

Parameters

temperature fluid temperature (K)

Returns: Henry's constant

Reimplemented in CO2FluidProperties, TabulatedFluidProperties, SimpleFluidProperties, IdealGasFluidPropertiesPT, HydrogenFluidProperties, NitrogenFluidProperties, and MethaneFluidProperties.

Definition at line 244 of file SinglePhaseFluidProperties.C.

Referenced by PorousFlowWaterNCG::enthalpyOfDissolution(), PorousFlowWaterNCG::equilibriumMassFractions(), TabulatedFluidProperties::henryConstant(), PorousFlowBrineCO2::henryConstant(), and SinglePhaseFluidProperties::henryConstant_dT().

 {
   mooseError(name(), ": henryConstant() is not implemented");
 }

◆ henryConstant() [2/2]

void SinglePhaseFluidProperties::henryConstant	(	Real	temperature,
		Real &	Kh,
		Real &	dKh_dT
	)		const

virtualinherited

Henry's law constant for dissolution in water and derivative wrt temperature.

Parameters

	temperature	fluid temperature (K)
[out]	Kh	Henry's constant
[out]	dKh_dT	derivative of Kh wrt temperature

Reimplemented in CO2FluidProperties, TabulatedFluidProperties, SimpleFluidProperties, IdealGasFluidPropertiesPT, HydrogenFluidProperties, NitrogenFluidProperties, and MethaneFluidProperties.

Definition at line 250 of file SinglePhaseFluidProperties.C.

 {
   mooseError(name(), ": henryConstant() is not implemented");
 }

◆ henryConstant_dT()

void SinglePhaseFluidProperties::henryConstant_dT	(	Real	temperature,
		Real &	Kh,
		Real &	dKh_dT
	)		const

virtualinherited

Definition at line 256 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": henryConstant_dT() is deprecated. Use henryConstant() instead");
 
   henryConstant(T, Kh, dKh_dT);
 }

◆ henryConstantIAPWS() [1/2]

Real SinglePhaseFluidProperties::henryConstantIAPWS	(	Real	temperature,
		Real	A,
		Real	B,
		Real	C
	)		const

protectedvirtualinherited

IAPWS formulation of Henry's law constant for dissolution in water From Guidelines on the Henry's constant and vapour liquid distribution constant for gases in H20 and D20 at high temperatures, IAPWS (2004)

Definition at line 161 of file SinglePhaseFluidProperties.C.

Referenced by MethaneFluidProperties::henryConstant(), NitrogenFluidProperties::henryConstant(), HydrogenFluidProperties::henryConstant(), CO2FluidProperties::henryConstant(), and SinglePhaseFluidProperties::henryConstantIAPWS_dT().

 {
   const Real Tr = T / 647.096;
   const Real tau = 1.0 - Tr;
 
   const Real lnkh =
       A / Tr + B * std::pow(tau, 0.355) / Tr + C * std::pow(Tr, -0.41) * std::exp(tau);
 
   // The vapor pressure used in this formulation
   const std::vector<Real> a{
       -7.85951783, 1.84408259, -11.7866497, 22.6807411, -15.9618719, 1.80122502};
   const std::vector<Real> b{1.0, 1.5, 3.0, 3.5, 4.0, 7.5};
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < a.size(); ++i)
     sum += a[i] * std::pow(tau, b[i]);
 
   return 22.064e6 * std::exp(sum / Tr) * std::exp(lnkh);
 }

◆ henryConstantIAPWS() [2/2]

void SinglePhaseFluidProperties::henryConstantIAPWS	(	Real	temperature,
		Real &	Kh,
		Real &	dKh_dT,
		Real	A,
		Real	B,
		Real	C
	)		const

protectedvirtualinherited

Definition at line 182 of file SinglePhaseFluidProperties.C.

 {
   const Real pc = 22.064e6;
   const Real Tc = 647.096;
 
   const Real Tr = T / Tc;
   const Real tau = 1.0 - Tr;
 
   const Real lnkh =
       A / Tr + B * std::pow(tau, 0.355) / Tr + C * std::pow(Tr, -0.41) * std::exp(tau);
   const Real dlnkh_dT =
       (-A / Tr / Tr - B * std::pow(tau, 0.355) / Tr / Tr - 0.355 * B * std::pow(tau, -0.645) / Tr -
        0.41 * C * std::pow(Tr, -1.41) * std::exp(tau) - C * std::pow(Tr, -0.41) * std::exp(tau)) /
       Tc;
 
   // The vapor pressure used in this formulation
   const std::vector<Real> a{
       -7.85951783, 1.84408259, -11.7866497, 22.6807411, -15.9618719, 1.80122502};
   const std::vector<Real> b{1.0, 1.5, 3.0, 3.5, 4.0, 7.5};
   Real sum = 0.0;
   Real dsum = 0.0;
 
   for (std::size_t i = 0; i < a.size(); ++i)
   {
     sum += a[i] * std::pow(tau, b[i]);
     dsum += a[i] * b[i] * std::pow(tau, b[i] - 1.0);
   }
 
   const Real p = pc * std::exp(sum / Tr);
   const Real dp_dT = -p / Tc / Tr * (sum / Tr + dsum);
 
   // Henry's constant and its derivative wrt temperature
   Kh = p * std::exp(lnkh);
   dKh_dT = (p * dlnkh_dT + dp_dT) * std::exp(lnkh);
 }

◆ henryConstantIAPWS_dT()

void SinglePhaseFluidProperties::henryConstantIAPWS_dT	(	Real	temperature,
		Real &	Kh,
		Real &	dKh_dT,
		Real	A,
		Real	B,
		Real	C
	)		const

protectedvirtualinherited

Definition at line 220 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(),
                   ":henryConstantIAPWS_dT() is deprecated. Use henryConstantIAPWS() instead");
 
   henryConstantIAPWS(T, Kh, dKh_dT, A, B, C);
 }

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Definition at line 28 of file FluidProperties.h.

28 {}

◆ inRegion()

unsigned int Water97FluidProperties::inRegion	(	Real	pressure,
		Real	temperature
	)		const

Determines the phase region that the given pressure and temperature values lie in.

Parameters

pressure	water pressure (Pa)
temperature	water temperature (K)

Returns: region phase region index

Definition at line 885 of file Water97FluidProperties.C.

Referenced by c_from_p_T(), cp_from_p_T(), cv_from_p_T(), e_from_p_T(), h_from_p_T(), rho_from_p_T(), and s_from_p_T().

 {
   // Valid for 273.15 K <= T <= 1073.15 K, p <= 100 MPa
   //          1073.15 K <= T <= 2273.15 K, p <= 50 Mpa
   if (temperature >= 273.15 && temperature <= 1073.15)
   {
     if (pressure < vaporPressure(273.15) || pressure > 100.0e6)
       throw MooseException("Pressure " + Moose::stringify(pressure) + " is out of range in " +
                            name() + ": inRegion()");
   }
   else if (temperature > 1073.15 && temperature <= 2273.15)
   {
     if (pressure < 0.0 || pressure > 50.0e6)
       throw MooseException("Pressure " + Moose::stringify(pressure) + " is out of range in " +
                            name() + ": inRegion()");
   }
   else
     throw MooseException("Temperature " + Moose::stringify(temperature) + " is out of range in " +
                          name() + ": inRegion()");
 
   // Determine the phase region that the (P, T) point lies in
   unsigned int region;
 
   if (temperature >= 273.15 && temperature <= 623.15)
   {
     if (pressure > vaporPressure(temperature) && pressure <= 100.0e6)
       region = 1;
     else
       region = 2;
   }
   else if (temperature > 623.15 && temperature <= 863.15)
   {
     if (pressure <= b23p(temperature))
       region = 2;
     else
       region = 3;
   }
   else if (temperature > 863.15 && temperature <= 1073.15)
     region = 2;
   else
     region = 5;
 
   return region;
 }

◆ inRegionPH()

unsigned int Water97FluidProperties::inRegionPH	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Determines the phase region that the given pressure and enthaply values lie in.

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: region phase region index

Definition at line 1552 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   unsigned int region;
 
   // Need to calculate enthalpies at the boundaries to delineate regions
   Real p273 = vaporPressure(273.15);
   Real p623 = vaporPressure(623.15);
 
   if (pressure >= p273 && pressure <= p623)
   {
     if (enthalpy >= h_from_p_T(pressure, 273.15) &&
         enthalpy <= h_from_p_T(pressure, vaporTemperature(pressure)))
       region = 1;
     else if (enthalpy > h_from_p_T(pressure, vaporTemperature(pressure)) &&
              enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else if (enthalpy > h_from_p_T(pressure, 1073.15) && enthalpy <= h_from_p_T(pressure, 2273.15))
       region = 5;
     else
       throw MooseException("Enthalpy " + Moose::stringify(enthalpy) + " is out of range in " +
                            name() + ": inRegionPH()");
   }
   else if (pressure > p623 && pressure <= 50.0e6)
   {
     if (enthalpy >= h_from_p_T(pressure, 273.15) && enthalpy <= h_from_p_T(pressure, 623.15))
       region = 1;
     else if (enthalpy > h_from_p_T(pressure, 623.15) &&
              enthalpy <= h_from_p_T(pressure, b23T(pressure)))
       region = 3;
     else if (enthalpy > h_from_p_T(pressure, b23T(pressure)) &&
              enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else if (enthalpy > h_from_p_T(pressure, 1073.15) && enthalpy <= h_from_p_T(pressure, 2273.15))
       region = 5;
     else
       throw MooseException("Enthalpy " + Moose::stringify(enthalpy) + " is out of range in " +
                            name() + ": inRegionPH()");
   }
   else if (pressure > 50.0e6 && pressure <= 100.0e6)
   {
     if (enthalpy >= h_from_p_T(pressure, 273.15) && enthalpy <= h_from_p_T(pressure, 623.15))
       region = 1;
     else if (enthalpy > h_from_p_T(pressure, 623.15) &&
              enthalpy <= h_from_p_T(pressure, b23T(pressure)))
       region = 3;
     else if (enthalpy > h_from_p_T(pressure, b23T(pressure)) &&
              enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else
       throw MooseException("Enthalpy " + Moose::stringify(enthalpy) + " is out of range in " +
                            name() + ": inRegionPH()");
   }
   else
     throw MooseException("Pressure " + Moose::stringify(pressure) + " is out of range in " +
                          name() + ": inRegionPH()");
 
   return region;
 }

◆ k()

Real SinglePhaseFluidProperties::k	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Thermal conductivity.

Parameters

pressure	fluid pressure (Pa)
temperature	fluid temperature (K)

Returns: thermal conductivity (W/m/K)

Definition at line 405 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::k_dpT(), IdealGasFluidPropertiesPT::k_from_p_T(), MethaneFluidProperties::k_from_p_T(), SimpleFluidProperties::k_from_p_T(), NitrogenFluidProperties::k_from_p_T(), HydrogenFluidProperties::k_from_p_T(), NaClFluidProperties::k_from_p_T(), IdealGasFluidProperties::k_from_p_T(), StiffenedGasFluidProperties::k_from_p_T(), TabulatedFluidProperties::k_from_p_T(), and CO2FluidProperties::k_from_p_T().

 {
   mooseDeprecated(name(), ": k() is deprecated. Use k_from_p_T() instead");
 
   return k_from_p_T(p, T);
 }

◆ k_dpT()

void SinglePhaseFluidProperties::k_dpT	(	Real	pressure,
		Real	temperature,
		Real &	k,
		Real &	dk_dp,
		Real &	dk_dT
	)		const

virtualinherited

Thermal conductivity and its derivatives wrt pressure and temperature.

Parameters

	pressure	fluid pressure (Pa)
	temperature	fluid temperature (K)
[out]	thermal	conductivity (W/m/K)
[out]	derivative	of thermal conductivity wrt pressure
[out]	derivative	of thermal conductivity wrt temperature

Definition at line 413 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": k_dpT() is deprecated. Use k_from_p_T() instead");
 
   k_from_p_T(p, T, k, dk_dp, dk_dT);
 }

◆ k_from_p_T() [1/2]

Real Water97FluidProperties::k_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 588 of file Water97FluidProperties.C.

 {
   Real rho = this->rho_from_p_T(pressure, temperature);
   return this->k_from_rho_T(rho, temperature);
 }

◆ k_from_p_T() [2/2]

void Water97FluidProperties::k_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	k,
		Real &	dk_dp,
		Real &	dk_dT
	)		const

overridevirtual

Definition at line 595 of file Water97FluidProperties.C.

 {
   mooseError(name(), ": k_from_p_T() is not implemented.");
 }

◆ k_from_rho_T()

Real Water97FluidProperties::k_from_rho_T	(	Real	density,
		Real	temperature
	)		const

overridevirtual

Definition at line 601 of file Water97FluidProperties.C.

Referenced by k_from_p_T().

 {
   // Scale the density and temperature. Note that the scales are slightly
   // different to the critical values used in IAPWS-IF97
   Real Tbar = temperature / 647.26;
   Real rhobar = density / 317.7;
 
   // Ideal gas component
   Real sum0 = 0.0;
 
   for (std::size_t i = 0; i < _k_a.size(); ++i)
     sum0 += _k_a[i] * MathUtils::pow(Tbar, i);
 
   Real lambda0 = std::sqrt(Tbar) * sum0;
 
   // The contribution due to finite density
   Real lambda1 = -0.39707 + 0.400302 * rhobar +
                  1.06 * std::exp(-0.171587 * Utility::pow<2>(rhobar + 2.392190));
 
   // Critical enhancement
   Real DeltaT = std::abs(Tbar - 1.0) + 0.00308976;
   Real Q = 2.0 + 0.0822994 / std::pow(DeltaT, 0.6);
   Real S = (Tbar >= 1.0 ? 1.0 / DeltaT : 10.0932 / std::pow(DeltaT, 0.6));
 
   Real lambda2 =
       (0.0701309 / Utility::pow<10>(Tbar) + 0.011852) * std::pow(rhobar, 1.8) *
           std::exp(0.642857 * (1.0 - std::pow(rhobar, 2.8))) +
       0.00169937 * S * std::pow(rhobar, Q) *
           std::exp((Q / (1.0 + Q)) * (1.0 - std::pow(rhobar, 1.0 + Q))) -
       1.02 * std::exp(-4.11717 * std::pow(Tbar, 1.5) - 6.17937 / Utility::pow<5>(rhobar));
 
   return lambda0 + lambda1 + lambda2;
 }

◆ molarMass()

Real Water97FluidProperties::molarMass ( ) const

overridevirtual

Molar mass [kg/mol].

Returns: molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 47 of file Water97FluidProperties.C.

 {
   return _Mh2o;
 }

◆ mu()

Real SinglePhaseFluidProperties::mu	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Dynamic viscosity.

Parameters

pressure	fluid pressure (Pa)
temperature	fluid temperature (K)

Returns: viscosity (Pa.s)

Definition at line 345 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": mu() is deprecated. Use mu_from_p_T() instead");
 
   return mu_from_p_T(p, T);
 }

◆ mu_dpT()

void SinglePhaseFluidProperties::mu_dpT	(	Real	pressure,
		Real	temperature,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

virtualinherited

Dynamic viscosity and its derivatives wrt pressure and temperature.

Parameters

	pressure	fluid pressure (Pa)
	temperature	fluid temperature (K)
[out]	mu	viscosity (Pa.s)
[out]	dmu_dp	derivative of viscosity wrt pressure
[out]	dmu_dT	derivative of viscosity wrt temperature

Definition at line 353 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": mu_dpT() is deprecated. Use mu_from_p_T() instead");
 
   mu_from_p_T(p, T, mu, dmu_dp, dmu_dT);
 }

◆ mu_drhoT_from_rho_T()

void SinglePhaseFluidProperties::mu_drhoT_from_rho_T	(	Real	density,
		Real	temperature,
		Real	ddensity_dT,
		Real &	mu,
		Real &	dmu_drho,
		Real &	dmu_dT
	)		const

virtualinherited

Definition at line 236 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ":mu_drhoT_from_rho_T() is deprecated. Use mu_from_rho_T() instead");
 
   mu_from_rho_T(rho, T, drho_dT, mu, dmu_drho, dmu_dT);
 }

◆ mu_from_p_T() [1/2]

Real Water97FluidProperties::mu_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 470 of file Water97FluidProperties.C.

 {
   Real rho = this->rho_from_p_T(pressure, temperature);
   return this->mu_from_rho_T(rho, temperature);
 }

◆ mu_from_p_T() [2/2]

void Water97FluidProperties::mu_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

overridevirtual

Definition at line 477 of file Water97FluidProperties.C.

 {
   Real rho, drho_dp, drho_dT;
   this->rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
   Real dmu_drho;
   this->mu_from_rho_T(rho, temperature, drho_dT, mu, dmu_drho, dmu_dT);
   dmu_dp = dmu_drho * drho_dp;
 }

◆ mu_from_rho_T() [1/2]

Real Water97FluidProperties::mu_from_rho_T	(	Real	density,
		Real	temperature
	)		const

overridevirtual

Definition at line 488 of file Water97FluidProperties.C.

Referenced by mu_from_p_T(), and rho_mu_from_p_T().

 {
   Real mu_star = 1.e-6;
   Real rhobar = density / _rho_critical;
   Real Tbar = temperature / _T_critical;
 
   // Viscosity in limit of zero density
   Real sum0 = 0.0;
   for (std::size_t i = 0; i < _mu_H0.size(); ++i)
     sum0 += _mu_H0[i] / MathUtils::pow(Tbar, i);
 
   Real mu0 = 100.0 * std::sqrt(Tbar) / sum0;
 
   // Residual component due to finite density
   Real sum1 = 0.0;
   for (std::size_t i = 0; i < _mu_H1.size(); ++i)
     sum1 += MathUtils::pow(1.0 / Tbar - 1.0, _mu_I[i]) * _mu_H1[i] *
             MathUtils::pow(rhobar - 1.0, _mu_J[i]);
 
   Real mu1 = std::exp(rhobar * sum1);
 
   // The water viscosity (in Pa.s) is then given by
   return mu_star * mu0 * mu1;
 }

◆ mu_from_rho_T() [2/2]

void Water97FluidProperties::mu_from_rho_T	(	Real	density,
		Real	temperature,
		Real	ddensity_dT,
		Real &	mu,
		Real &	dmu_drho,
		Real &	dmu_dT
	)		const

overridevirtual

Dynamic viscosity and its derivatives wrt density and temperature TODO: this shouldn't need 3 input args - AD will assume/call the 2-input version.

Parameters

	density	fluid density (kg/m^3)
	temperature	fluid temperature (K)
	ddensity_dT	derivative of density wrt temperature
[out]	mu	viscosity (Pa.s)
[out]	dmu_drho	derivative of viscosity wrt density
[out]	dmu_dT	derivative of viscosity wrt temperature

Reimplemented from SinglePhaseFluidProperties.

Definition at line 514 of file Water97FluidProperties.C.

 {
   Real mu_star = 1.0e-6;
   Real rhobar = density / _rho_critical;
   Real Tbar = temperature / _T_critical;
   Real drhobar_drho = 1.0 / _rho_critical;
   Real dTbar_dT = 1.0 / _T_critical;
 
   // Limit of zero density. Derivative wrt rho is 0
   Real sum0 = 0.0, dsum0_dTbar = 0.0;
   for (std::size_t i = 0; i < _mu_H0.size(); ++i)
   {
     sum0 += _mu_H0[i] / MathUtils::pow(Tbar, i);
     dsum0_dTbar -= i * _mu_H0[i] / MathUtils::pow(Tbar, i + 1);
   }
 
   Real mu0 = 100.0 * std::sqrt(Tbar) / sum0;
   Real dmu0_dTbar =
       50.0 / std::sqrt(Tbar) / sum0 - 100.0 * std::sqrt(Tbar) * dsum0_dTbar / sum0 / sum0;
 
   // Residual component due to finite density
   Real sum1 = 0.0, dsum1_drho = 0.0, dsum1_dTbar = 0.0;
   for (std::size_t i = 0; i < _mu_H1.size(); ++i)
   {
     sum1 += MathUtils::pow(1.0 / Tbar - 1.0, _mu_I[i]) * _mu_H1[i] *
             MathUtils::pow(rhobar - 1.0, _mu_J[i]);
     dsum1_drho += MathUtils::pow(1.0 / Tbar - 1.0, _mu_I[i]) * _mu_H1[i] * _mu_J[i] *
                   MathUtils::pow(rhobar - 1.0, _mu_J[i] - 1);
     dsum1_dTbar -= _mu_I[i] * MathUtils::pow(1.0 / Tbar - 1.0, _mu_I[i] - 1) * _mu_H1[i] *
                    MathUtils::pow(rhobar - 1.0, _mu_J[i]) / Tbar / Tbar;
   }
 
   Real mu1 = std::exp(rhobar * sum1);
   Real dmu1_drho = (sum1 + rhobar * dsum1_drho) * mu1;
   Real dmu1_dTbar = (rhobar * dsum1_dTbar) * mu1;
 
   // Viscosity and its derivatives are then
   mu = mu_star * mu0 * mu1;
   dmu_drho = mu_star * mu0 * dmu1_drho * drhobar_drho;
   dmu_dT = mu_star * (dmu0_dTbar * mu1 + mu0 * dmu1_dTbar) * dTbar_dT + dmu_drho * ddensity_dT;
 }

◆ phi3()

Real Water97FluidProperties::phi3	(	Real	delta,
		Real	tau
	)		const

protected

Helmholtz free energy in Region 3.

From Eq. (28) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, IAPWS 2007.

Parameters

delta	reduced density (-)
tau	reduced temperature (-)

Returns: Helmholtz free energy (-)

Definition at line 1091 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

 {
   Real sum = 0.0;
   for (std::size_t i = 1; i < _n3.size(); ++i)
     sum += _n3[i] * MathUtils::pow(delta, _I3[i]) * MathUtils::pow(tau, _J3[i]);
 
   return _n3[0] * std::log(delta) + sum;
 }

◆ propfunc() [1/15]

SinglePhaseFluidProperties::propfunc	(	p	,
		v	,
		e
	)

inherited

Compute a fluid property given for the state defined by two given properties.

For all functions, the first two arguments are the given properties that define the fluid state. For the two-argument variants, the desired property is the return value. The five-argument variants also provide partial derivatives dx/da and dx/db where x is the desired property being computed, a is the first given property, and b is the second given property. The desired property, dx/da, and dx/db are stored into the 3rd, 4th, and 5th arguments respectively.

Properties/parameters used in these function are listed below with their units:

p pressure [Pa] T temperature [K] e specific internal energy [J/kg] v specific volume [m^3/kg] rho density [kg/m^3] h specific enthalpy [J/kg] s specific entropy [J/(kg*K)] mu viscosity [Pa*s] k thermal conductivity [W/(m*K)] c speed of sound [m/s] cp constant-pressure specific heat [J/K] cv constant-volume specific heat [J/K] beta volumetric thermal expansion coefficient [1/K] g Gibbs free energy [J] pp_sat partial pressure at saturation [Pa]

As an example:

// calculate pressure given specific vol and energy: auto pressure = your_fluid_properties_object.p_from_v_e(specific_vol, specific_energy);

// or use the derivative variant: Real dp_dv = 0; // derivative will be stored into here Real dp_de = 0; // derivative will be stored into here your_fluid_properties_object.p_from_v_e(specific_vol, specific_energy, pressure, dp_dv, dp_de);

Automatic differentiation (AD) support is provided through x_from_a_b(DualReal a, DualReal b) versions of the functions where a and b must be ADReal/DualNumber's calculated using all AD-supporting values:

auto v = 1/rho; // rho must be an AD non-linear variable. auto e = rhoE/rho - vel_energy; // rhoE and vel_energy must be AD variables/numbers also. auto pressure = your_fluid_properties_object.p_from_v_e(v, e); // pressure now contains partial derivatives w.r.t. all degrees of freedom

◆ propfunc() [2/15]

e SinglePhaseFluidProperties::propfunc	(	c	,
		v	,
		e
	)

inherited

◆ propfunc() [3/15]

e e SinglePhaseFluidProperties::propfunc	(	cv	,
		v	,
		e
	)

inherited

◆ propfunc() [4/15]

e e e SinglePhaseFluidProperties::propfunc	(	k	,
		v	,
		e
	)

inherited

◆ propfunc() [5/15]

e e e e SinglePhaseFluidProperties::propfunc	(	s	,
		h	,
		p
	)

inherited

◆ propfunc() [6/15]

e e e e s SinglePhaseFluidProperties::propfunc	(	e	,
		v	,
		h
	)

inherited

◆ propfunc() [7/15]

e e e e s T SinglePhaseFluidProperties::propfunc	(	pp_sat	,
		p	,
		T
	)

inherited

◆ propfunc() [8/15]

e e e e s T T SinglePhaseFluidProperties::propfunc	(	k	,
		rho	,
		T
	)

inherited

◆ propfunc() [9/15]

e e e e s T T T SinglePhaseFluidProperties::propfunc	(	cp	,
		p	,
		T
	)

inherited

◆ propfunc() [10/15]

e e e e s T T T T SinglePhaseFluidProperties::propfunc	(	mu	,
		p	,
		T
	)

inherited

◆ propfunc() [11/15]

e e e e s T T T T T SinglePhaseFluidProperties::propfunc	(	rho	,
		p	,
		T
	)

inherited

◆ propfunc() [12/15]

e e e e s T T T T T rho SinglePhaseFluidProperties::propfunc	(	e	,
		T	,
		v
	)

inherited

◆ propfunc() [13/15]

e e e e s T T T T T rho v SinglePhaseFluidProperties::propfunc	(	h	,
		T	,
		v
	)

inherited

◆ propfunc() [14/15]

e e e e s T T T T T rho v v SinglePhaseFluidProperties::propfunc	(	cv	,
		T	,
		v
	)

inherited

◆ propfunc() [15/15]

e e e e s T T T T T rho v v T SinglePhaseFluidProperties::propfunc	(	p	,
		h	,
		s
	)

inherited

◆ rho()

virtual Real SinglePhaseFluidProperties::rho	(	Real	p,
		Real	T
	)		const

virtualinherited

Density from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)

Returns: density (kg/m^3)

◆ rho_dpT()

void SinglePhaseFluidProperties::rho_dpT	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT
	)		const

virtualinherited

Density and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure (Pa)
[in]	T	temperature (K)
[out]	rho	density (kg/m^3)
[out]	drho_dp	derivative of density w.r.t. pressure
[out]	drho_dT	derivative of density w.r.t. temperature

Definition at line 283 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": rho_dpT() is deprecated. Use rho_from_p_T() instead");
 
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
 }

◆ rho_e_dpT()

void SinglePhaseFluidProperties::rho_e_dpT	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

virtualinherited

Reimplemented in IdealGasFluidPropertiesPT.

Definition at line 322 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": rho_e_dpT() is deprecated. Use rho_e_from_p_T() instead");
 
   rho_e_from_p_T(p, T, rho, drho_dp, drho_dT, e, de_dp, de_dT);
 }

◆ rho_e_from_p_T()

void SinglePhaseFluidProperties::rho_e_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

virtualinherited

Density and internal energy and their derivatives wrt pressure and temperature.

Parameters

	pressure	fluid pressure (Pa)
	temperature	fluid temperature (K)
[out]	rho	density (kg/m^3)
[out]	drho_dp	derivative of density wrt pressure
[out]	drho_dT	derivative of density wrt temperature
[out]	e	internal energy (J/kg)
[out]	de_dp	derivative of internal energy wrt pressure
[out]	de_dT	derivative of internal energy wrt temperature

Definition at line 308 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::rho_e_dpT().

 {
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
   e_from_p_T(p, T, e, de_dp, de_dT);
 }

◆ rho_from_p_T() [1/2]

Real Water97FluidProperties::rho_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 83 of file Water97FluidProperties.C.

Referenced by k_from_p_T(), mu_from_p_T(), rho_from_p_T(), and rho_mu_from_p_T().

 {
   Real density, pi, tau;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
 
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       density = pressure / (pi * _Rw * temperature * dgamma1_dpi(pi, tau));
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       density = pressure / (pi * _Rw * temperature * dgamma2_dpi(pi, tau));
       break;
 
     case 3:
       density = densityRegion3(pressure, temperature);
       break;
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       density = pressure / (pi * _Rw * temperature * dgamma5_dpi(pi, tau));
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
   return density;
 }

◆ rho_from_p_T() [2/2]

void Water97FluidProperties::rho_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT
	)		const

overridevirtual

Definition at line 121 of file Water97FluidProperties.C.

 {
   Real pi, tau, ddensity_dp, ddensity_dT;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
 
   switch (region)
   {
     case 1:
     {
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       Real dgdp = dgamma1_dpi(pi, tau);
       ddensity_dp = -d2gamma1_dpi2(pi, tau) / (_Rw * temperature * dgdp * dgdp);
       ddensity_dT = -pressure * (dgdp - tau * d2gamma1_dpitau(pi, tau)) /
                     (_Rw * pi * temperature * temperature * dgdp * dgdp);
       break;
     }
 
     case 2:
     {
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       Real dgdp = dgamma2_dpi(pi, tau);
       ddensity_dp = -d2gamma2_dpi2(pi, tau) / (_Rw * temperature * dgdp * dgdp);
       ddensity_dT = -pressure * (dgdp - tau * d2gamma2_dpitau(pi, tau)) /
                     (_Rw * pi * temperature * temperature * dgdp * dgdp);
       break;
     }
 
     case 3:
     {
       // Calculate density first, then use that in Helmholtz free energy
       Real density = densityRegion3(pressure, temperature);
       Real delta = density / _rho_critical;
       tau = _T_star[2] / temperature;
       Real dpdd = dphi3_ddelta(delta, tau);
       Real d2pdd2 = d2phi3_ddelta2(delta, tau);
       ddensity_dp = 1.0 / (_Rw * temperature * delta * (2.0 * dpdd + delta * d2pdd2));
       ddensity_dT = density * (tau * d2phi3_ddeltatau(delta, tau) - dpdd) / temperature /
                     (2.0 * dpdd + delta * d2pdd2);
       break;
     }
 
     case 5:
     {
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       Real dgdp = dgamma5_dpi(pi, tau);
       ddensity_dp = -d2gamma5_dpi2(pi, tau) / (_Rw * temperature * dgdp * dgdp);
       ddensity_dT = -pressure * (dgdp - tau * d2gamma5_dpitau(pi, tau)) /
                     (_Rw * pi * temperature * temperature * dgdp * dgdp);
       break;
     }
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
 
   rho = this->rho_from_p_T(pressure, temperature);
   drho_dp = ddensity_dp;
   drho_dT = ddensity_dT;
 }

◆ rho_mu()

void SinglePhaseFluidProperties::rho_mu	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	mu
	)		const

virtualinherited

Density and viscosity.

Parameters

	pressure	fluid pressure (Pa)
	temperature	fluid temperature (K)
[out]	rho	density (kg/m^3)
[out]	mu	viscosity (Pa.s)

Definition at line 361 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": rho_mu() is deprecated. Use rho_mu_from_p_T() instead");
 
   rho_mu_from_p_T(p, T, rho, mu);
 }

◆ rho_mu_dpT()

void SinglePhaseFluidProperties::rho_mu_dpT	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

virtualinherited

Definition at line 376 of file SinglePhaseFluidProperties.C.

 {
   mooseDeprecated(name(), ": rho_mu_dpT() is deprecated. Use rho_mu_from_p_T() instead");
 
   rho_mu_from_p_T(p, T, rho, drho_dp, drho_dT, mu, dmu_dp, dmu_dT);
 }

◆ rho_mu_from_p_T() [1/2]

void Water97FluidProperties::rho_mu_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	mu
	)		const

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 562 of file Water97FluidProperties.C.

 {
   rho = this->rho_from_p_T(pressure, temperature);
   mu = this->mu_from_rho_T(rho, temperature);
 }

◆ rho_mu_from_p_T() [2/2]

void Water97FluidProperties::rho_mu_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

overridevirtual

Density and viscosity and their derivatives wrt pressure and temperature.

Parameters

	pressure	fluid pressure (Pa)
	temperature	fluid temperature (K)
[out]	rho	density (kg/m^3)
[out]	drho_dp	derivative of density wrt pressure
[out]	drho_dT	derivative of density wrt temperature
[out]	mu	viscosity (Pa.s)
[out]	dmu_dp	derivative of viscosity wrt pressure
[out]	dmu_dT	derivative of viscosity wrt temperature

Reimplemented from SinglePhaseFluidProperties.

Definition at line 572 of file Water97FluidProperties.C.

 {
   this->rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
   Real dmu_drho;
   this->mu_from_rho_T(rho, temperature, drho_dT, mu, dmu_drho, dmu_dT);
   dmu_dp = dmu_drho * drho_dp;
 }

◆ s()

Real SinglePhaseFluidProperties::s	(	Real	pressure,
		Real	temperature
	)		const

virtualinherited

Definition at line 421 of file SinglePhaseFluidProperties.C.

Referenced by IdealGasFluidProperties::p_from_h_s(), StiffenedGasFluidProperties::p_from_h_s(), IdealGasFluidProperties::rho_from_p_s(), StiffenedGasFluidProperties::rho_from_p_s(), StiffenedGasFluidProperties::s_from_h_p(), IdealGasFluidProperties::s_from_h_p(), IdealGasFluidProperties::s_from_p_T(), StiffenedGasFluidProperties::s_from_p_T(), IdealGasFluidPropertiesPT::s_from_p_T(), HelmholtzFluidProperties::s_from_p_T(), SimpleFluidProperties::s_from_p_T(), s_from_p_T(), TabulatedFluidProperties::s_from_p_T(), StiffenedGasFluidProperties::s_from_T_v(), IdealGasFluidProperties::s_from_T_v(), IdealGasFluidProperties::s_from_v_e(), StiffenedGasFluidProperties::s_from_v_e(), and SinglePhaseFluidProperties::T_from_p_h().

 {
   mooseDeprecated(name(), ": s() is deprecated. Use s_from_p_T() instead");
 
   return s_from_p_T(p, T);
 }

◆ s_from_p_T() [1/2]

Real Water97FluidProperties::s_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 636 of file Water97FluidProperties.C.

 {
   Real entropy, pi, tau, density3, delta;
 
   // Determine which region the point is in
   unsigned int region = inRegion(pressure, temperature);
   switch (region)
   {
     case 1:
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       entropy = _Rw * (tau * dgamma1_dtau(pi, tau) - gamma1(pi, tau));
       break;
 
     case 2:
       pi = pressure / _p_star[1];
       tau = _T_star[1] / temperature;
       entropy = _Rw * (tau * dgamma2_dtau(pi, tau) - gamma2(pi, tau));
       break;
 
     case 3:
       // Calculate density first, then use that in Helmholtz free energy
       density3 = densityRegion3(pressure, temperature);
       delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       entropy = _Rw * (tau * dphi3_dtau(delta, tau) - phi3(delta, tau));
       break;
 
     case 5:
       pi = pressure / _p_star[4];
       tau = _T_star[4] / temperature;
       entropy = _Rw * (tau * dgamma5_dtau(pi, tau) - gamma5(pi, tau));
       break;
 
     default:
       mooseError(name(), ": inRegion() has given an incorrect region");
   }
   return entropy;
 }

◆ s_from_p_T() [2/2]

void Water97FluidProperties::s_from_p_T	(	Real	p,
		Real	T,
		Real &	s,
		Real &	ds_dp,
		Real &	ds_dT
	)		const

overridevirtual

Definition at line 677 of file Water97FluidProperties.C.

 {
   SinglePhaseFluidProperties::s_from_p_T(p, T, s, ds_dp, ds_dT);
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 32 of file FluidProperties.h.

32 {}

◆ subregion2ph()

unsigned int Water97FluidProperties::subregion2ph	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Provides the correct subregion index for a (P,h) point in region 2.

From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: subregion index

Definition at line 1612 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   unsigned int subregion;
 
   if (pressure <= 4.0e6)
     subregion = 1;
   else if (pressure > 4.0e6 && pressure < 6.5467e6)
     subregion = 2;
   else
   {
     if (enthalpy >= b2bc(pressure))
       subregion = 2;
     else
       subregion = 3;
   }
 
   return subregion;
 }

◆ subregion3()

unsigned int Water97FluidProperties::subregion3	(	Real	pressure,
		Real	temperature
	)		const

Provides the correct subregion index for a (P,T) point in region 3.

From Revised Supplementary Release on Backward Equations for Specific Volume as a Function of Pressure and Temperature v(p,T) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
temperature	water temperature (K)

Returns: subregion index

Definition at line 1249 of file Water97FluidProperties.C.

Referenced by densityRegion3().

 {
   Real pMPa = pressure / 1.0e6;
   const Real P3cd = 19.00881189173929;
   unsigned int subregion = 0;
 
   if (pMPa > 40.0 && pMPa <= 100.0)
   {
     if (temperature <= tempXY(pressure, AB))
       subregion = 0;
     else // (temperature > tempXY(pressure, AB))
       subregion = 1;
   }
   else if (pMPa > 25.0 && pMPa <= 40.0)
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= tempXY(pressure, AB))
       subregion = 3;
     else if (temperature > tempXY(pressure, AB) && temperature <= tempXY(pressure, EF))
       subregion = 4;
     else // (temperature > tempXY(pressure, EF))
       subregion = 5;
   }
   else if (pMPa > 23.5 && pMPa <= 25.0)
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= tempXY(pressure, GH))
       subregion = 6;
     else if (temperature > tempXY(pressure, GH) && temperature <= tempXY(pressure, EF))
       subregion = 7;
     else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, IJ))
       subregion = 8;
     else if (temperature > tempXY(pressure, IJ) && temperature <= tempXY(pressure, JK))
       subregion = 9;
     else // (temperature > tempXY(pressure, JK))
       subregion = 10;
   }
   else if (pMPa > 23.0 && pMPa <= 23.5)
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= tempXY(pressure, GH))
       subregion = 11;
     else if (temperature > tempXY(pressure, GH) && temperature <= tempXY(pressure, EF))
       subregion = 7;
     else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, IJ))
       subregion = 8;
     else if (temperature > tempXY(pressure, IJ) && temperature <= tempXY(pressure, JK))
       subregion = 9;
     else // (temperature > tempXY(pressure, JK))
       subregion = 10;
   }
   else if (pMPa > 22.5 && pMPa <= 23.0)
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= tempXY(pressure, GH))
       subregion = 11;
     else if (temperature > tempXY(pressure, GH) && temperature <= tempXY(pressure, MN))
       subregion = 12;
     else if (temperature > tempXY(pressure, MN) && temperature <= tempXY(pressure, EF))
       subregion = 13;
     else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, OP))
       subregion = 14;
     else if (temperature > tempXY(pressure, OP) && temperature <= tempXY(pressure, IJ))
       subregion = 15;
     else if (temperature > tempXY(pressure, IJ) && temperature <= tempXY(pressure, JK))
       subregion = 9;
     else // (temperature > tempXY(pressure, JK))
       subregion = 10;
   }
   else if (pMPa > vaporPressure(643.15) * 1.0e-6 &&
            pMPa <= 22.5) // vaporPressure(643.15) = 21.04 MPa
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= tempXY(pressure, QU))
       subregion = 16;
     else if (temperature > tempXY(pressure, QU) && temperature <= tempXY(pressure, RX))
     {
       if (pMPa > 22.11 && pMPa <= 22.5)
       {
         if (temperature <= tempXY(pressure, UV))
           subregion = 20;
         else if (temperature > tempXY(pressure, UV) && temperature <= tempXY(pressure, EF))
           subregion = 21;
         else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, WX))
           subregion = 22;
         else // (temperature > tempXY(pressure, WX) && temperature <= tempXY(pressure, RX))
           subregion = 23;
       }
       else if (pMPa > 22.064 && pMPa <= 22.11)
       {
         if (temperature <= tempXY(pressure, UV))
           subregion = 20;
         else if (temperature > tempXY(pressure, UV) && temperature <= tempXY(pressure, EF))
           subregion = 24;
         else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, WX))
           subregion = 25;
         else // (temperature > tempXY(pressure, WX) && temperature <= tempXY(pressure, RX))
           subregion = 23;
       }
       else if (temperature <= vaporTemperature(pressure))
       {
         if (pMPa > 21.93161551 && pMPa <= 22.064)
           if (temperature > tempXY(pressure, QU) && temperature <= tempXY(pressure, UV))
             subregion = 20;
           else
             subregion = 24;
         else // (pMPa > vaporPressure(643.15) * 1.0e-6 && pMPa <= 21.93161551)
           subregion = 20;
       }
       else if (temperature > vaporTemperature(pressure))
       {
         if (pMPa > 21.90096265 && pMPa <= 22.064)
         {
           if (temperature <= tempXY(pressure, WX))
             subregion = 25;
           else
             subregion = 23;
         }
         else
           subregion = 23;
       }
     }
     else if (temperature > tempXY(pressure, RX) && temperature <= tempXY(pressure, JK))
       subregion = 17;
     else
       subregion = 10;
   }
   else if (pMPa > 20.5 &&
            pMPa <= vaporPressure(643.15) * 1.0e-6) // vaporPressure(643.15) = 21.04 MPa
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= vaporTemperature(pressure))
       subregion = 16;
     else if (temperature > vaporTemperature(pressure) && temperature <= tempXY(pressure, JK))
       subregion = 17;
     else // (temperature > tempXY(pressure, JK))
       subregion = 10;
   }
   else if (pMPa > P3cd && pMPa <= 20.5) // P3cd  = 19.00881189173929
   {
     if (temperature <= tempXY(pressure, CD))
       subregion = 2;
     else if (temperature > tempXY(pressure, CD) && temperature <= vaporTemperature(pressure))
       subregion = 18;
     else
       subregion = 19;
   }
   else if (pMPa > vaporPressure(623.15) * 1.0e-6 && pMPa <= P3cd)
   {
     if (temperature < vaporTemperature(pressure))
       subregion = 2;
     else
       subregion = 19;
   }
   else if (pMPa > 22.11 && pMPa <= 22.5)
   {
     if (temperature > tempXY(pressure, QU) && temperature <= tempXY(pressure, UV))
       subregion = 20;
     else if (temperature > tempXY(pressure, UV) && temperature <= tempXY(pressure, EF))
       subregion = 21;
     else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, WX))
       subregion = 22;
     else // (temperature > tempXY(pressure, WX) && temperature <= tempXY(pressure, RX))
       subregion = 23;
   }
   else if (pMPa > 22.064 && pMPa <= 22.11)
   {
     if (temperature > tempXY(pressure, QU) && temperature <= tempXY(pressure, UV))
       subregion = 20;
     else if (temperature > tempXY(pressure, UV) && temperature <= tempXY(pressure, EF))
       subregion = 24;
     else if (temperature > tempXY(pressure, EF) && temperature <= tempXY(pressure, WX))
       subregion = 25;
     else // (temperature > tempXY(pressure, WX) && temperature <= tempXY(pressure, RX))
       subregion = 23;
   }
   else
     mooseError(name(), ": subregion3(): Shouldn't have got here!");
 
   return subregion;
 }

◆ subregion3ph()

unsigned int Water97FluidProperties::subregion3ph	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Provides the correct subregion index for a (P,h) point in region 3.

From Revised Supplementary Release on Backward Equations for the Functions T(p,h), v(p,h) and T(p,s), v(p,s) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: subregion index

Definition at line 1632 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   unsigned int subregion;
 
   if (enthalpy <= b3ab(pressure))
     subregion = 1;
   else
     subregion = 2;
 
   return subregion;
 }

◆ subregionVolume()

Real Water97FluidProperties::subregionVolume	(	Real	pi,
		Real	theta,
		Real	a,
		Real	b,
		Real	c,
		Real	d,
		Real	e,
		unsigned int	sid
	)		const

Specific volume in all subregions of region 3 EXCEPT subregion n (13).

Parameters

pi	scaled water pressure
theta	scaled water temperature
a	to e constants
sid	subregion ID of the subregion

Returns: volume water specific volume (m^3/kg)

Definition at line 1499 of file Water97FluidProperties.C.

Referenced by densityRegion3().

 {
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _n3s[sid].size(); ++i)
     sum += _n3s[sid][i] * MathUtils::pow(std::pow(pi - a, c), _I3s[sid][i]) *
            MathUtils::pow(std::pow(theta - b, d), _J3s[sid][i]);
 
   return std::pow(sum, e);
 }

◆ T_from_p_h()

Real SinglePhaseFluidProperties::T_from_p_h	(	Real	pressure,
		Real	enthalpy
	)		const

virtualinherited

Temperature from pressure and specific enthalpy.

Parameters

[in]	pressure	pressure (Pa)
[in]	enthalpy	enthalpy (J/kg)

Returns: Temperature (K)

Reimplemented in IdealGasFluidProperties.

Definition at line 464 of file SinglePhaseFluidProperties.C.

 {
   const Real s = s_from_h_p(h, p);
   const Real rho = rho_from_p_s(p, s);
   const Real v = 1. / rho;
   const Real e = e_from_v_h(v, h);
   return T_from_v_e(v, e);
 }

◆ temperature_from_ph()

Real Water97FluidProperties::temperature_from_ph	(	Real	pressure,
		Real	enthalpy
	)		const

virtual

Backwards equation T(p, h) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam.

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1658 of file Water97FluidProperties.C.

 {
   Real temperature = 0.0;
 
   // Determine which region the point is in
   unsigned int region = inRegionPH(pressure, enthalpy);
 
   switch (region)
   {
     case 1:
       temperature = temperature_from_ph1(pressure, enthalpy);
       break;
 
     case 2:
     {
       // First, determine which subregion the point is in:
       unsigned int subregion = subregion2ph(pressure, enthalpy);
 
       if (subregion == 1)
         temperature = temperature_from_ph2a(pressure, enthalpy);
       else if (subregion == 2)
         temperature = temperature_from_ph2b(pressure, enthalpy);
       else
         temperature = temperature_from_ph2c(pressure, enthalpy);
       break;
     }
 
     case 3:
     {
       // First, determine which subregion the point is in:
       unsigned int subregion = subregion3ph(pressure, enthalpy);
 
       if (subregion == 1)
         temperature = temperature_from_ph3a(pressure, enthalpy);
       else
         temperature = temperature_from_ph3b(pressure, enthalpy);
       break;
     }
 
     case 5:
       mooseError(name(), ": temperature_from_ph() not implemented for region 5");
       break;
 
     default:
       mooseError(name(), ": inRegionPH() has given an incorrect region");
   }
 
   return temperature;
 }

◆ temperature_from_ph1()

Real Water97FluidProperties::temperature_from_ph1	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 1 Eq.

(11) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1709 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 1.0e6;
   Real eta = enthalpy / 2500.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph1.size(); ++i)
     sum += _nph1[i] * MathUtils::pow(pi, _Iph1[i]) * MathUtils::pow(eta + 1.0, _Jph1[i]);
 
   return sum;
 }

◆ temperature_from_ph2a()

Real Water97FluidProperties::temperature_from_ph2a	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 2a Eq.

(22) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1722 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 1.0e6;
   Real eta = enthalpy / 2000.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph2a.size(); ++i)
     sum += _nph2a[i] * MathUtils::pow(pi, _Iph2a[i]) * MathUtils::pow(eta - 2.1, _Jph2a[i]);
 
   return sum;
 }

◆ temperature_from_ph2b()

Real Water97FluidProperties::temperature_from_ph2b	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 2b Eq.

(23) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1735 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 1.0e6;
   Real eta = enthalpy / 2000.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph2b.size(); ++i)
     sum += _nph2b[i] * MathUtils::pow(pi - 2.0, _Iph2b[i]) * MathUtils::pow(eta - 2.6, _Jph2b[i]);
 
   return sum;
 }

◆ temperature_from_ph2c()

Real Water97FluidProperties::temperature_from_ph2c	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 2c Eq.

(24) from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1748 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 1.0e6;
   Real eta = enthalpy / 2000.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph2c.size(); ++i)
     sum += _nph2c[i] * MathUtils::pow(pi + 25.0, _Iph2c[i]) * MathUtils::pow(eta - 1.8, _Jph2c[i]);
 
   return sum;
 }

◆ temperature_from_ph3a()

Real Water97FluidProperties::temperature_from_ph3a	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 3a Eq.

(2) from Revised Supplementary Release on Backward Equations for the Functions T(p,h), v(p,h) and T(p,s), v(p,s) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1776 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 100.0e6;
   Real eta = enthalpy / 2300.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph3a.size(); ++i)
     sum +=
         _nph3a[i] * MathUtils::pow(pi + 0.24, _Iph3a[i]) * MathUtils::pow(eta - 0.615, _Jph3a[i]);
 
   return sum * 760.0;
 }

◆ temperature_from_ph3b()

Real Water97FluidProperties::temperature_from_ph3b	(	Real	pressure,
		Real	enthalpy
	)		const

protected

Backwards equation T(p, h) in Region 3b Eq.

(3) from Revised Supplementary Release on Backward Equations for the Functions T(p,h), v(p,h) and T(p,s), v(p,s) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
enthalpy	water enthalpy (J/kg)

Returns: temperature water temperature (K)

Definition at line 1790 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

 {
   Real pi = pressure / 100.0e6;
   Real eta = enthalpy / 2800.0e3;
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < _nph3b.size(); ++i)
     sum +=
         _nph3b[i] * MathUtils::pow(pi + 0.298, _Iph3b[i]) * MathUtils::pow(eta - 0.72, _Jph3b[i]);
 
   return sum * 860.0;
 }

◆ tempXY()

Real Water97FluidProperties::tempXY	(	Real	pressure,
		subregionEnum	xy
	)		const

protected

Boundaries between subregions in region 3.

From Revised Supplementary Release on Backward Equations for Specific Volume as a Function of Pressure and Temperature v(p,T) for Region 3 of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam

Parameters

pressure	water pressure (Pa)
xy	string to select the boundary between two subregions

Returns: temperature (K) along the boundary

Definition at line 1438 of file Water97FluidProperties.C.

Referenced by subregion3().

 {
   Real pi = pressure / 1.0e6;
 
   // Choose the constants based on the string xy
   unsigned int row;
 
   switch (xy)
   {
     case AB:
       row = 0;
       break;
     case CD:
       row = 1;
       break;
     case GH:
       row = 2;
       break;
     case IJ:
       row = 3;
       break;
     case JK:
       row = 4;
       break;
     case MN:
       row = 5;
       break;
     case OP:
       row = 6;
       break;
     case QU:
       row = 7;
       break;
     case RX:
       row = 8;
       break;
     case UV:
       row = 9;
       break;
     case WX:
       row = 10;
       break;
     default:
       row = 0;
   }
 
   Real sum = 0.0;
 
   if (xy == AB || xy == OP || xy == WX)
     for (std::size_t i = 0; i < _tempXY_n[row].size(); ++i)
       sum += _tempXY_n[row][i] * MathUtils::pow(std::log(pi), _tempXY_I[row][i]);
   else if (xy == EF)
     sum += 3.727888004 * (pi - _p_critical / 1.0e6) + _T_critical;
   else
     for (std::size_t i = 0; i < _tempXY_n[row].size(); ++i)
       sum += _tempXY_n[row][i] * MathUtils::pow(pi, _tempXY_I[row][i]);
 
   return sum;
 }

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Definition at line 31 of file FluidProperties.h.

31 {}

◆ triplePointPressure()

Real Water97FluidProperties::triplePointPressure ( ) const

overridevirtual

Triple point pressure.

Returns: triple point pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 71 of file Water97FluidProperties.C.

 {
   return _p_triple;
 }

◆ triplePointTemperature()

Real Water97FluidProperties::triplePointTemperature ( ) const

overridevirtual

Triple point temperature.

Returns: triple point temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 77 of file Water97FluidProperties.C.

 {
   return _T_triple;
 }

◆ v_e_spndl_from_T()

void SinglePhaseFluidProperties::v_e_spndl_from_T	(	Real	T,
		Real &	v,
		Real &	e
	)		const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in StiffenedGasFluidProperties.

Definition at line 458 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::v_from_p_T().

 {
   mooseError(name(), ": ", __PRETTY_FUNCTION__, " not implemented.");
 }

◆ v_from_p_T() [1/2]

Real SinglePhaseFluidProperties::v_from_p_T	(	Real	p,
		Real	T
	)		const

virtualinherited

Specific volume from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature

Definition at line 60 of file SinglePhaseFluidProperties.C.

Referenced by GeneralVaporMixtureFluidProperties::c_from_p_T(), GeneralVaporMixtureFluidProperties::cp_from_p_T(), GeneralVaporMixtureFluidProperties::cv_from_p_T(), GeneralVaporMixtureFluidProperties::k_from_p_T(), GeneralVaporMixtureFluidProperties::mu_from_p_T(), GeneralVaporMixtureFluidProperties::v_from_p_T(), and IdealRealGasMixtureFluidProperties::xs_prim_from_p_T().

 {
   const Real rho = rho_from_p_T(p, T);
   return 1.0 / rho;
 }

◆ v_from_p_T() [2/2]

void SinglePhaseFluidProperties::v_from_p_T	(	Real	p,
		Real	T,
		Real &	v,
		Real &	dv_dp,
		Real &	dv_dT
	)		const

virtualinherited

Specific volume and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[out]	v	specific volume
[out]	dv_dp	derivative of specific volume w.r.t. pressure
[out]	dv_dT	derivative of specific volume w.r.t. temperature

Definition at line 67 of file SinglePhaseFluidProperties.C.

 {
   Real rho, drho_dp, drho_dT;
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
 
   v = 1.0 / rho;
   const Real dv_drho = -1.0 / (rho * rho);
 
   dv_dp = dv_drho * drho_dp;
   dv_dT = dv_drho * drho_dT;
 }

◆ vaporPressure() [1/2]

Real Water97FluidProperties::vaporPressure ( Real temperature ) const

overridevirtual

Vapor pressure.

Used to delineate liquid and gas phases. Valid for temperatures between the triple point temperature and the critical temperature

Parameters

temperature water temperature (K)

Returns: saturation pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 787 of file Water97FluidProperties.C.

Referenced by inRegion(), inRegionPH(), and subregion3().

 {
   // Check whether the input temperature is within the region of validity of this equation.
   // Valid for 273.15 K <= t <= 647.096 K
   if (temperature < 273.15 || temperature > _T_critical)
     throw MooseException(name() + ": vaporPressure(): Temperature is outside range 273.15 K <= T "
                                   "<= 647.096 K");

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ subregionEnum

Constructor & Destructor Documentation

◆ Water97FluidProperties()

◆ ~Water97FluidProperties()

Member Function Documentation

◆ b23p()

◆ b23T()

◆ b2bc()

◆ b3ab()

◆ beta()

◆ beta_from_p_T() [1/2]

◆ beta_from_p_T() [2/2]

◆ c()

◆ c_from_p_T()

◆ cp_from_p_T()

◆ criticalDensity()

◆ criticalInternalEnergy()

◆ criticalPressure()

◆ criticalTemperature()

◆ cv_from_p_T()

◆ d2gamma1_dpi2()

◆ d2gamma1_dpitau()

◆ d2gamma1_dtau2()

◆ d2gamma2_dpi2()

◆ d2gamma2_dpitau()

◆ d2gamma2_dtau2()

◆ d2gamma5_dpi2()

◆ d2gamma5_dpitau()

◆ d2gamma5_dtau2()

◆ d2phi3_ddelta2()

◆ d2phi3_ddeltatau()

◆ d2phi3_dtau2()

◆ densityRegion3()

◆ dgamma1_dpi()

◆ dgamma1_dtau()

◆ dgamma2_dpi()

◆ dgamma2_dtau()

◆ dgamma5_dpi()

◆ dgamma5_dtau()

◆ dphi3_ddelta()

◆ dphi3_dtau()

◆ e()

◆ e_dpT()

◆ e_from_p_T() [1/2]

◆ e_from_p_T() [2/2]

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ gamma1()

◆ gamma2()

◆ gamma5()

◆ gamma_from_p_T()

◆ h()

◆ h_dpT()

◆ h_from_p_T() [1/2]

◆ h_from_p_T() [2/2]

◆ henryConstant() [1/2]

◆ henryConstant() [2/2]

◆ henryConstant_dT()

◆ henryConstantIAPWS() [1/2]

◆ henryConstantIAPWS() [2/2]

◆ henryConstantIAPWS_dT()

◆ initialize()

◆ inRegion()

◆ inRegionPH()

◆ k()

◆ k_dpT()

◆ k_from_p_T() [1/2]

◆ k_from_p_T() [2/2]

◆ k_from_rho_T()

◆ molarMass()

◆ mu()

◆ mu_dpT()