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

virtual Real	molarMass () const override
	Fluid name. 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

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

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

void	mu_from_rho_T (Real rho, Real temperature, Real drho_dT, Real &mu, Real &dmu_drho, Real &dmu_dT) const

virtual void	rho_mu_from_p_T (Real pressure, Real temperature, Real &rho, Real &mu) const override
	Combined methods. More...

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

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

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

virtual void	vaporTemperature (Real pressure, Real &Tsat, Real &dTsat_dp) const override

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	T_from_p_h (Real pressure, Real enthalpy) const override
	Backwards equation T(p, h) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. More...

virtual void	T_from_p_h (Real p, Real h, Real &T, Real &dT_dp, Real &dT_dh) const override

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

Real	henryConstant (Real temperature, const std::vector< Real > &coeffs) 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...

void	henryConstant (Real temperature, const std::vector< Real > &coeffs, Real &Kh, Real &dKh_dT) const

DualReal	henryConstant (const DualReal &temperature, const std::vector< Real > &coeffs) const

virtual Real	criticalInternalEnergy () const
	Critical specific internal energy. 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...

DualReal	vaporPressure (const DualReal &T) const

DualReal	vaporTemperature (const DualReal &p) const

virtual std::vector< Real >	henryCoefficients () const
	Henry's law coefficients for dissolution in water. More...

virtual void	rho_mu_from_p_T (const DualReal &p, const DualReal &T, DualReal &rho, DualReal &mu) const

virtual void	rho_e_from_p_T (Real p, Real T, Real &rho, Real &drho_dp, Real &drho_dT, Real &e, Real &de_dp, Real &de_dT) const

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

Static Public Attributes
static const Real	_R = 8.3144598
	Universal gas constant (J/mol/K) More...

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

FPDualReal	T_from_p_h_ad (const FPDualReal &pressure, const FPDualReal &enthalpy) const
	AD version of backwards equation T(p, h) (used internally) From Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. More...

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

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

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

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

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

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

FPDualReal	vaporTemperature_ad (const FPDualReal &pressure) const
	AD version of saturation temperature as a function of pressure (used internally) More...

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	_T_c2k
	Conversion of temperature from Celsius to Kelvin. More...

const bool	_allow_imperfect_jacobians
	Flag to set unimplemented Jacobian entries to zero. More...

Private Member Functions
template<typename... Args>
void	fluidPropError (Args... args) const

e e e e p h T T T T T v v v s h	propfuncWithDefault (beta, p, T) propfuncWithDefault(v

e e e e p h T T T T T v v v s h T	propfuncWithDefault (e, p, T) propfuncWithDefault(gamma

	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(T

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

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

e e e e p h T	propfunc (mu, rho, T) propfunc(k

e e e e p h T T	propfunc (c, p, T) propfunc(cp

e e e e p h T T T	propfunc (cv, p, T) propfunc(mu

e e e e p h T T T T	propfunc (k, p, T) propfunc(rho

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

e e e e p h T T T T T v	propfunc (p, T, v) propfunc(h

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

e e e e p h T T T T T v v v	propfunc (h, p, T) propfunc(p

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

	v

e	v

e e	v

e e e	v

e e e e p	v

e e e e p h T T T T T v v v s h T	v

e e e e	h

e e e e p h T T T T T v v v	h

e e e e p h	p

e e e e p h T T	p

e e e e p h T T T	p

e e e e p h T T T T	p

e e e e p h T T T T T v v v s	p

e e e e p h T T T T T v v v s h	p

e e e e p h T	rho

e e e e p h T T T T T	T

e e e e p h T T T T T v	T

e e e e p h T T T T T v v	T

e e e e p h T T T T T v v v s h T e	propfuncWithDefault (gamma, p, 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 42 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 494 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 25 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 37 of file Water97FluidProperties.C.

37 {}

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 884 of file Water97FluidProperties.C.

 {
   // 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)
     mooseException(name(),
                    ": b23p(): Temperature ",
                    temperature,
                    " is outside range of 623.15 K <= T <= 863.15 K");
  
   return (_n23[0] + _n23[1] * temperature + _n23[2] * temperature * temperature) * 1.e6;
 }

Referenced by b3ab(), and inRegion().

◆ 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 898 of file Water97FluidProperties.C.

 {
   // 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)
     mooseException(
         name(), ": b23T(): Pressure ", pressure, "is outside range 16.529 MPa <= p <= 100 MPa");
  
   return _n23[3] + std::sqrt((pressure / 1.e6 - _n23[4]) / _n23[2]);
 }

Referenced by inRegionPH().

◆ 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 1663 of file Water97FluidProperties.C.

 {
   // Check whether the input pressure is within the region of validity of this equation.
   if (pressure < 6.5467e6 || pressure > 100.0e6)
     mooseException(
         name(), ": b2bc(): Pressure ", 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;
 }

Referenced by subregion2ph().

◆ 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 1823 of file Water97FluidProperties.C.

 {
   // Check whether the input pressure is within the region of validity of this equation.
   if (pressure < b23p(623.15) || pressure > 100.0e6)
     mooseException(
         name(), ": b3ab(): Pressure ", 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;
 }

Referenced by subregion3ph().

◆ c_from_p_T()

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

overridevirtual

Definition at line 310 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 368 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 64 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 IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 126 of file SinglePhaseFluidProperties.C.

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

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e(), and IdealRealGasMixtureFluidProperties::T_from_p_v().

◆ criticalPressure()

Real Water97FluidProperties::criticalPressure ( ) const

overridevirtual

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 52 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 58 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 417 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 974 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1007 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 996 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1050 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1098 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1081 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1208 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1256 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1239 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1133 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1165 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1154 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1534 of file Water97FluidProperties.C.

 {
   // 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;
 }

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().

◆ 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 963 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 985 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1036 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1065 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1194 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1223 of file Water97FluidProperties.C.

 {
   // 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;
 }

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

◆ 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 1123 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 1144 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ e_from_p_T() [1/2]

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

overridevirtual

Definition at line 187 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by e_from_p_T().

◆ e_from_p_T() [2/2]

void Water97FluidProperties::e_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

overridevirtual

Definition at line 234 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 IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 286 of file SinglePhaseFluidProperties.C.

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

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e(), and IdealRealGasMixtureFluidProperties::T_from_p_v().

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Definition at line 34 of file FluidProperties.h.

34 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 36 of file FluidProperties.h.

36 {}

◆ fluidName()

std::string Water97FluidProperties::fluidName ( ) const

overridevirtual

Definition at line 40 of file Water97FluidProperties.C.

 {
   return "water";
 }

◆ fluidPropError()

template<typename... Args>

void SinglePhaseFluidProperties::fluidPropError ( Args... args ) const

inlineprivateinherited

Definition at line 326 of file SinglePhaseFluidProperties.h.

   {
     if (_allow_imperfect_jacobians)
       mooseDoOnce(mooseWarning(std::forward<Args>(args)...));
     else
       mooseError(std::forward<Args>(args)...);
   }

Referenced by SinglePhaseFluidProperties::vaporPressure(), and SinglePhaseFluidProperties::vaporTemperature().

◆ 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 953 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by s_from_p_T().

◆ 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 1018 of file Water97FluidProperties.C.

 {
   // 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;
 }

Referenced by s_from_p_T().

◆ 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 1176 of file Water97FluidProperties.C.

 {
   // 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;
 }

Referenced by s_from_p_T().

◆ h_from_p_T() [1/2]

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

overridevirtual

Definition at line 682 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by inRegionPH().

◆ 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 726 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;
 }

◆ henryCoefficients()

std::vector< Real > SinglePhaseFluidProperties::henryCoefficients ( ) const

virtualinherited

Henry's law coefficients for dissolution in water.

Returns: Henry's constant coefficients

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

Definition at line 183 of file SinglePhaseFluidProperties.C.

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

Referenced by TabulatedFluidProperties::henryCoefficients().

◆ henryConstant() [1/3]

DualReal Water97FluidProperties::henryConstant	(	const DualReal &	temperature,
		const std::vector< Real > &	coeffs
	)		const

Definition at line 1935 of file Water97FluidProperties.C.

 {
   const Real T = temperature.value();
   Real Kh_real = 0.0;
   Real dKh_dT_real = 0.0;
   henryConstant(T, coeffs, Kh_real, dKh_dT_real);
  
   DualReal Kh = Kh_real;
   Kh.derivatives() = temperature.derivatives() * dKh_dT_real;
  
   return Kh;
 }

◆ henryConstant() [2/3]

Real Water97FluidProperties::henryConstant	(	Real	temperature,
		const std::vector< Real > &	coeffs
	)		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)

Parameters

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

Definition at line 1866 of file Water97FluidProperties.C.

 {
   const Real A = coeffs[0];
   const Real B = coeffs[1];
   const Real C = coeffs[2];
  
   const Real Tr = temperature / 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);
 }

Referenced by PorousFlowWaterNCG::enthalpyOfDissolution(), PorousFlowWaterNCG::equilibriumMassFractions(), BrineFluidProperties::henryConstant(), and henryConstant().

◆ henryConstant() [3/3]

void Water97FluidProperties::henryConstant	(	Real	temperature,
		const std::vector< Real > &	coeffs,
		Real &	Kh,
		Real &	dKh_dT
	)		const

Definition at line 1891 of file Water97FluidProperties.C.

 {
   const Real A = coeffs[0];
   const Real B = coeffs[1];
   const Real C = coeffs[2];
  
   const Real pc = 22.064e6;
   const Real Tc = 647.096;
  
   const Real Tr = temperature / 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);
 }

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Definition at line 35 of file FluidProperties.h.

35 {}

◆ 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 910 of file Water97FluidProperties.C.

 {
   // 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)
       mooseException("Pressure ", pressure, " is out of range in ", name(), ": inRegion()");
   }
   else if (temperature > 1073.15 && temperature <= 2273.15)
   {
     if (pressure < 0.0 || pressure > 50.0e6)
       mooseException("Pressure ", pressure, " is out of range in ", name(), ": inRegion()");
   }
   else
     mooseException("Temperature ", 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;
 }

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().

◆ 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 1574 of file Water97FluidProperties.C.

 {
   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
       mooseException("Enthalpy ", 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
       mooseException("Enthalpy ", 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
       mooseException("Enthalpy ", enthalpy, " is out of range in ", name(), ": inRegionPH()");
   }
   else
     mooseException("Pressure ", pressure, " is out of range in ", name(), ": inRegionPH()");
  
   return region;
 }

Referenced by T_from_p_h_ad().

◆ k_from_p_T() [1/2]

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

overridevirtual

Definition at line 587 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 594 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 600 of file Water97FluidProperties.C.

 {
   // 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;
 }

Referenced by k_from_p_T().

◆ molarMass()

Real Water97FluidProperties::molarMass ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name Molar mass [kg/mol]; molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 46 of file Water97FluidProperties.C.

 {
   return _Mh2o;
 }

◆ mu_from_p_T() [1/2]

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

overridevirtual

Definition at line 469 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 476 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 487 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by mu_from_p_T(), and rho_mu_from_p_T().

◆ mu_from_rho_T() [2/2]

void Water97FluidProperties::mu_from_rho_T	(	Real	rho,
		Real	temperature,
		Real	drho_dT,
		Real &	mu,
		Real &	dmu_drho,
		Real &	dmu_dT
	)		const

Definition at line 513 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 1113 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by s_from_p_T().

◆ propfunc() [1/16]

e e e e p h T T SinglePhaseFluidProperties::propfunc	(	c	,
		p	,
		T
	)

inherited

◆ propfunc() [2/16]

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

inherited

◆ propfunc() [3/16]

e e e e p h T T T SinglePhaseFluidProperties::propfunc	(	cv	,
		p	,
		T
	)

inherited

◆ propfunc() [4/16]

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

inherited

◆ propfunc() [5/16]

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

inherited

◆ propfunc() [6/16]

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

inherited

◆ propfunc() [7/16]

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

inherited

◆ propfunc() [8/16]

e e e e p h T T T T SinglePhaseFluidProperties::propfunc	(	k	,
		p	,
		T
	)

inherited

◆ propfunc() [9/16]

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

inherited

◆ propfunc() [10/16]

e e e e p h T SinglePhaseFluidProperties::propfunc	(	mu	,
		rho	,
		T
	)

inherited

◆ propfunc() [11/16]

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

inherited

◆ propfunc() [12/16]

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:

@begincode 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] gamma Adiabatic ratio (cp/cv) [-]

As an example:

@begincode // 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) and x_from_a_b(DualReal a, DualReal b, DualReal x, DualReal dx_da, DualReal dx_db) versions of the functions where a and b must be ADReal/DualNumber's calculated using all AD-supporting values:

@begincode 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() [13/16]

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

inherited

◆ propfunc() [14/16]

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

inherited

◆ propfunc() [15/16]

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

inherited

◆ propfunc() [16/16]

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

inherited

◆ propfuncWithDefault() [1/2]

e e e e p h T T T T T v v v s h SinglePhaseFluidProperties::propfuncWithDefault	(	beta	,
		p	,
		T
	)

inherited

◆ propfuncWithDefault() [2/2]

e e e e p h T T T T T v v v s h T SinglePhaseFluidProperties::propfuncWithDefault	(	e	,
		p	,
		T
	)

inherited

◆ rho_e_from_p_T()

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

virtualinherited

Definition at line 242 of file SinglePhaseFluidProperties.C.

 {
   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 82 of file Water97FluidProperties.C.

 {
   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;
 }

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

◆ 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 120 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_from_p_T() [1/3]

void SinglePhaseFluidProperties::rho_mu_from_p_T	(	const DualReal &	p,
		const DualReal &	T,
		DualReal &	rho,
		DualReal &	mu
	)		const

virtualinherited

Definition at line 277 of file SinglePhaseFluidProperties.C.

 {
   rho = rho_from_p_T(p, T);
   mu = mu_from_p_T(p, T);
 }

◆ rho_mu_from_p_T() [2/3]

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 571 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;
 }

◆ rho_mu_from_p_T() [3/3]

void Water97FluidProperties::rho_mu_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	mu
	)		const

overridevirtual

Combined methods.

These methods are particularly useful for the PorousFlow module, where density and viscosity are typically both computed everywhere. The combined methods allow the most efficient means of calculating both properties, especially where rho(p, T) and mu(rho, T). In this case, an extra density calculation would be required to calculate mu(p, T). All propery names are described above.

Reimplemented from SinglePhaseFluidProperties.

Definition at line 561 of file Water97FluidProperties.C.

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

◆ s_from_p_T() [1/2]

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

overridevirtual

Definition at line 676 of file Water97FluidProperties.C.

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

◆ s_from_p_T() [2/2]

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

overridevirtual

Definition at line 635 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;
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 39 of file FluidProperties.h.

39 {}

◆ 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 1630 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by T_from_p_h_ad().

◆ 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 1271 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by densityRegion3().

◆ 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 1650 of file Water97FluidProperties.C.

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

Referenced by T_from_p_h_ad().

◆ 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 1521 of file Water97FluidProperties.C.

 {
   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);
 }

Referenced by densityRegion3().

◆ T_from_p_h() [1/2]

void Water97FluidProperties::T_from_p_h	(	Real	p,
		Real	h,
		Real &	T,
		Real &	dT_dp,
		Real &	dT_dh
	)		const

overridevirtual

Definition at line 1685 of file Water97FluidProperties.C.

 {
   FPDualReal p = pressure;
   Moose::derivInsert(p.derivatives(), 0, 1.0);
   FPDualReal h = enthalpy;
   Moose::derivInsert(h.derivatives(), 1, 1.0);
  
   const FPDualReal T = T_from_p_h_ad(p, h);
  
   temperature = T.value();
   dT_dp = T.derivatives()[0];
   dT_dh = T.derivatives()[1];
 }

◆ T_from_p_h() [2/2]

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

overridevirtual

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 1676 of file Water97FluidProperties.C.

 {
   const FPDualReal p = pressure;
   const FPDualReal h = enthalpy;
  
   return T_from_p_h_ad(p, h).value();
 }

◆ T_from_p_h_ad()

FPDualReal Water97FluidProperties::T_from_p_h_ad	(	const FPDualReal &	pressure,
		const FPDualReal &	enthalpy
	)		const

protected

AD version of backwards equation T(p, h) (used internally) 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 1701 of file Water97FluidProperties.C.

 {
   FPDualReal temperature = 0.0;
  
   // Determine which region the point is in
   const unsigned int region = inRegionPH(pressure.value(), enthalpy.value());
  
   switch (region)
   {
     case 1:
       temperature = temperature_from_ph1(pressure, enthalpy);
       break;
  
     case 2:
     {
       // First, determine which subregion the point is in:
       const unsigned int subregion = subregion2ph(pressure.value(), enthalpy.value());
  
       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:
       const unsigned int subregion = subregion3ph(pressure.value(), enthalpy.value());
  
       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;
 }

Referenced by T_from_p_h().

◆ temperature_from_ph1()

FPDualReal Water97FluidProperties::temperature_from_ph1	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1753 of file Water97FluidProperties.C.

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

Referenced by T_from_p_h_ad().

◆ temperature_from_ph2a()

FPDualReal Water97FluidProperties::temperature_from_ph2a	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1767 of file Water97FluidProperties.C.

 {
   const FPDualReal pi = pressure / 1.0e6;
   const FPDualReal eta = enthalpy / 2000.0e3;
   FPDualReal sum = 0.0;
  
   // Factor out the negative in std::pow(eta - 2.1, _Jph2a[i]) to avoid fpe in dbg (see #13163)
   const Real sgn = MathUtils::sign(eta.value() - 2.1);
  
   for (std::size_t i = 0; i < _nph2a.size(); ++i)
     sum += _nph2a[i] * std::pow(pi, _Iph2a[i]) * std::pow(std::abs(eta - 2.1), _Jph2a[i]) *
            std::pow(sgn, _Jph2a[i]);
  
   return sum;
 }

Referenced by T_from_p_h_ad().

◆ temperature_from_ph2b()

FPDualReal Water97FluidProperties::temperature_from_ph2b	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1785 of file Water97FluidProperties.C.

 {
   const FPDualReal pi = pressure / 1.0e6;
   const FPDualReal eta = enthalpy / 2000.0e3;
   FPDualReal sum = 0.0;
  
   // Factor out the negatives in std::pow(pi - 2.0, _Iph2b[i])* std::pow(eta - 2.6, _Jph2b[i])
   // to avoid fpe in dbg (see #13163)
   const Real sgn0 = MathUtils::sign(pi.value() - 2.0);
   const Real sgn1 = MathUtils::sign(eta.value() - 2.6);
  
   for (std::size_t i = 0; i < _nph2b.size(); ++i)
     sum += _nph2b[i] * std::pow(std::abs(pi - 2.0), _Iph2b[i]) * std::pow(sgn0, _Iph2b[i]) *
            std::pow(std::abs(eta - 2.6), _Jph2b[i]) * std::pow(sgn1, _Jph2b[i]);
  
   return sum;
 }

Referenced by T_from_p_h_ad().

◆ temperature_from_ph2c()

FPDualReal Water97FluidProperties::temperature_from_ph2c	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1805 of file Water97FluidProperties.C.

 {
   const FPDualReal pi = pressure / 1.0e6;
   const FPDualReal eta = enthalpy / 2000.0e3;
   FPDualReal sum = 0.0;
  
   // Factor out the negative in std::pow(eta - 1.8, _Jph2c[i]) to avoid fpe in dbg (see #13163)
   const Real sgn = MathUtils::sign(eta.value() - 1.8);
  
   for (std::size_t i = 0; i < _nph2c.size(); ++i)
     sum += _nph2c[i] * std::pow(pi + 25.0, _Iph2c[i]) * std::pow(std::abs(eta - 1.8), _Jph2c[i]) *
            std::pow(sgn, _Jph2c[i]);
  
   return sum;
 }

Referenced by T_from_p_h_ad().

◆ temperature_from_ph3a()

FPDualReal Water97FluidProperties::temperature_from_ph3a	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1838 of file Water97FluidProperties.C.

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

Referenced by T_from_p_h_ad().

◆ temperature_from_ph3b()

FPDualReal Water97FluidProperties::temperature_from_ph3b	(	const FPDualReal &	pressure,
		const FPDualReal &	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 1852 of file Water97FluidProperties.C.

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

Referenced by T_from_p_h_ad().

◆ 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 1460 of file Water97FluidProperties.C.

 {
   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;
 }

Referenced by subregion3().

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Definition at line 38 of file FluidProperties.h.

38 {}

◆ triplePointPressure()

Real Water97FluidProperties::triplePointPressure ( ) const

overridevirtual

Triple point pressure.

Returns: triple point pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 70 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 76 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 IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 292 of file SinglePhaseFluidProperties.C.

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

Referenced by IdealRealGasMixtureFluidProperties::v_from_p_T().

◆ vaporPressure() [1/3]

DualReal SinglePhaseFluidProperties::vaporPressure ( const DualReal & T ) const

inherited

Definition at line 198 of file SinglePhaseFluidProperties.C.

 {
   Real p = 0.0;
   Real temperature = T.value();
   Real dpdT = 0.0;
  
   vaporPressure(temperature, p, dpdT);
  
   DualReal result = p;
   result.derivatives() = T.derivatives() * dpdT;
  
   return result;
 }

◆ vaporPressure() [2/3]

Real Water97FluidProperties::vaporPressure ( Real T ) const

overridevirtual

Vapor pressure.

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

Parameters

	T	fluid temperature (K)
[out]	saturation	pressure (Pa)
[out]	derivative	of saturation pressure wrt temperature (Pa/K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 786 of file Water97FluidProperties.C.

 {
   // 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)
     mooseException(name(),
                    ": vaporPressure(): Temperature ",
                    temperature,
                    " is outside range 273.15 K <= T "
                    "<= 647.096 K");
  
   Real theta, theta2, a, b, c, p;
   theta = temperature + _n4[8] / (temperature - _n4[9]);
   theta2 = theta * theta;
   a = theta2 + _n4[0] * theta + _n4[1];
   b = _n4[2] * theta2 + _n4[3] * theta + _n4[4];
   c = _n4[5] * theta2 + _n4[6] * theta + _n4[7];
   p = Utility::pow<4>(2.0 * c / (-b + std::sqrt(b * b - 4.0 * a * c)));
  
   return p * 1.e6;
 }

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

◆ vaporPressure() [3/3]

void Water97FluidProperties::vaporPressure	(	Real	temperature,
		Real &	psat,
		Real &	dpsat_dT
	)		const

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 809 of file Water97FluidProperties.C.

 {
   // 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)
     mooseException(name(),
                    ": vaporPressure(): Temperature is outside range 273.15 K <= T <= 647.096 K");
  
   Real theta, dtheta_dT, theta2, a, b, c, da_dtheta, db_dtheta, dc_dtheta;
   theta = temperature + _n4[8] / (temperature - _n4[9]);
   dtheta_dT = 1.0 - _n4[8] / (temperature - _n4[9]) / (temperature - _n4[9]);
   theta2 = theta * theta;
  
   a = theta2 + _n4[0] * theta + _n4[1];
   b = _n4[2] * theta2 + _n4[3] * theta + _n4[4];
   c = _n4[5] * theta2 + _n4[6] * theta + _n4[7];
  
   da_dtheta = 2.0 * theta + _n4[0];
   db_dtheta = 2.0 * _n4[2] * theta + _n4[3];
   dc_dtheta = 2.0 * _n4[5] * theta + _n4[6];
  
   Real denominator = -b + std::sqrt(b * b - 4.0 * a * c);
  
   psat = Utility::pow<4>(2.0 * c / denominator) * 1.0e6;
  
   // The derivative wrt temperature is given by the chain rule
   Real dpsat = 4.0 * Utility::pow<3>(2.0 * c / denominator);
   dpsat *= (2.0 * dc_dtheta / denominator -
             2.0 * c / denominator / denominator *
                 (-db_dtheta + std::pow(b * b - 4.0 * a * c, -0.5) *
                                   (b * db_dtheta - 2.0 * da_dtheta * c - 2.0 * a * dc_dtheta)));
   dpsat_dT = dpsat * dtheta_dT * 1.0e6;
 }

◆ vaporTemperature() [1/3]

DualReal SinglePhaseFluidProperties::vaporTemperature ( const DualReal & p ) const

inherited

Definition at line 227 of file SinglePhaseFluidProperties.C.

 {
   Real T = 0.0;
   Real pressure = p.value();
   Real dTdp = 0.0;
  
   vaporTemperature(pressure, T, dTdp);
  
   DualReal result = T;
   result.derivatives() = p.derivatives() * dTdp;
  
   return result;
 }

◆ vaporTemperature() [2/3]

Real Water97FluidProperties::vaporTemperature ( Real pressure ) const

overridevirtual

Saturation temperature as a function of pressure.

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

Valid for 611.213 Pa <= p <= 22.064 MPa

Parameters

pressure water pressure (Pa)

Returns: saturation temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 864 of file Water97FluidProperties.C.

 {
   const FPDualReal p = pressure;
  
   return vaporTemperature_ad(p).value();
 }

Referenced by inRegionPH(), and subregion3().

◆ vaporTemperature() [3/3]

void Water97FluidProperties::vaporTemperature	(	Real	pressure,
		Real &	Tsat,
		Real &	dTsat_dp
	)		const

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 872 of file Water97FluidProperties.C.

 {
   FPDualReal p = pressure;
   Moose::derivInsert(p.derivatives(), 0, 1.0);
  
   const FPDualReal T = vaporTemperature_ad(p);
  
   Tsat = T.value();
   dTsat_dp = T.derivatives()[0];
 }

◆ vaporTemperature_ad()

FPDualReal Water97FluidProperties::vaporTemperature_ad ( const FPDualReal & pressure ) const

protected

AD version of saturation temperature as a function of pressure (used internally)

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

Valid for 611.213 Pa <= p <= 22.064 MPa

Parameters

pressure water pressure (Pa)

Returns: saturation temperature (K)

Definition at line 844 of file Water97FluidProperties.C.

 {
   // Check whether the input pressure is within the region of validity of this equation.
   // Valid for 611.213 Pa <= p <= 22.064 MPa
   if (pressure.value() < 611.23 || pressure.value() > _p_critical)
     mooseException(name() + ": vaporTemperature(): Pressure ",
                    pressure.value(),
                    " is outside range 611.213 Pa <= p <= 22.064 MPa");
  
   const FPDualReal beta = std::pow(pressure / 1.e6, 0.25);
   const FPDualReal beta2 = beta * beta;
   const FPDualReal e = beta2 + _n4[2] * beta + _n4[5];
   const FPDualReal f = _n4[0] * beta2 + _n4[3] * beta + _n4[6];
   const FPDualReal g = _n4[1] * beta2 + _n4[4] * beta + _n4[7];
   const FPDualReal d = 2.0 * g / (-f - std::sqrt(f * f - 4.0 * e * g));
  
   return (_n4[9] + d - std::sqrt((_n4[9] + d) * (_n4[9] + d) - 4.0 * (_n4[8] + _n4[9] * d))) / 2.0;
 }

Referenced by vaporTemperature().

Member Data Documentation

◆ _allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians

protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 48 of file FluidProperties.h.

Referenced by SinglePhaseFluidProperties::fluidPropError().

◆ _I1

const std::array<int, 34> Water97FluidProperties::_I1

protected

Initial value:

{{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2,

2, 2, 3, 3, 3, 4, 4, 4, 5, 8, 8, 21, 23, 29, 30, 31, 32}}

Definition at line 692 of file Water97FluidProperties.h.

Referenced by d2gamma1_dpi2(), d2gamma1_dpitau(), d2gamma1_dtau2(), dgamma1_dpi(), dgamma1_dtau(), and gamma1().

◆ _I2

const std::array<int, 43> Water97FluidProperties::_I2

protected

Initial value:

{{1,  1,  1,  1,  1,  2,  2,  2,  2,  2,  3,  3,  3, 3,  3,
                                 4,  4,  4,  5,  6,  6,  6,  7,  7,  7,  8,  8,  9, 10, 10,
                                 10, 16, 16, 18, 20, 20, 20, 21, 22, 23, 24, 24, 24}}

Definition at line 737 of file Water97FluidProperties.h.

Referenced by d2gamma2_dpi2(), d2gamma2_dpitau(), d2gamma2_dtau2(), dgamma2_dpi(), dgamma2_dtau(), and gamma2().

◆ _I3

const std::array<int, 40> Water97FluidProperties::_I3

protected

Initial value:

{{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3,

3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 8, 9, 9, 10, 10, 11}}

Definition at line 816 of file Water97FluidProperties.h.

Referenced by d2phi3_ddelta2(), d2phi3_ddeltatau(), d2phi3_dtau2(), dphi3_ddelta(), dphi3_dtau(), and phi3().

◆ _I3s

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

protected

Definition at line 1102 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _I5

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

protected

Definition at line 1261 of file Water97FluidProperties.h.

Referenced by d2gamma5_dpi2(), d2gamma5_dpitau(), d2gamma5_dtau2(), dgamma5_dpi(), dgamma5_dtau(), and gamma5().

◆ _Iph1

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

protected

Definition at line 706 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _Iph2a

const std::array<int, 36> Water97FluidProperties::_Iph2a

protected

Initial value:

{{0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2,

2, 2, 2, 2, 2, 2, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 7}}

Definition at line 756 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _Iph2b

const std::array<int, 38> Water97FluidProperties::_Iph2b

protected

Initial value:

{{0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2,

2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 6, 7, 7, 9, 9}}

Definition at line 775 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _Iph2c

const std::array<int, 23> Water97FluidProperties::_Iph2c

protected

Initial value:

{

{-7, -7, -6, -6, -5, -5, -2, -2, -1, -1, 0, 0, 1, 1, 2, 6, 6, 6, 6, 6, 6, 6, 6}}

Definition at line 790 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _Iph3a

const std::array<int, 31> Water97FluidProperties::_Iph3a

protected

Initial value:

{{-12, -12, -12, -12, -12, -12, -12, -12, -10, -10, -10,
                                    -8,  -8,  -8,  -8,  -5,  -3,  -2,  -2,  -2,  -1,  -1,
                                    0,   0,   1,   3,   3,   4,   4,   10,  12}}

Definition at line 1215 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _Iph3b

const std::array<int, 33> Water97FluidProperties::_Iph3b

protected

Initial value:

{{-12, -12, -10, -10, -10, -10, -10, -8, -8, -8, -8,
                                    -8,  -6,  -6,  -6,  -4,  -4,  -3,  -2, -2, -1, -1,
                                    -1,  -1,  -1,  -1,  0,   0,   1,   3,  5,  6,  8}}

Definition at line 1233 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _J02

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

protected

Definition at line 735 of file Water97FluidProperties.h.

Referenced by d2gamma2_dtau2(), dgamma2_dtau(), and gamma2().

◆ _J05

const std::array<int, 6> Water97FluidProperties::_J05 {{0, 1, -3, -2, -1, 2}}

protected

Constants for region 5.

Definition at line 1253 of file Water97FluidProperties.h.

Referenced by d2gamma5_dtau2(), dgamma5_dtau(), and gamma5().

◆ _J1

const std::array<int, 34> Water97FluidProperties::_J1

protected

Initial value:

{{-2, -1, 0,   1,  2,   3,   4,   5,   -9,  -7, -1, 0,
                                 1,  3,  -3,  0,  1,   3,   17,  -4,  0,   6,  -5, -2,
                                 10, -8, -11, -6, -29, -31, -38, -39, -40, -41}}

Definition at line 695 of file Water97FluidProperties.h.

Referenced by d2gamma1_dpi2(), d2gamma1_dpitau(), d2gamma1_dtau2(), dgamma1_dpi(), dgamma1_dtau(), and gamma1().

◆ _J2

const std::array<int, 43> Water97FluidProperties::_J2

protected

Initial value:

{{0,  1,  2,  3,  6,  1,  2,  4,  7,  36, 0,  1,  3,  6, 35,
                                 1,  2,  3,  7,  3,  16, 35, 0,  11, 25, 8,  36, 13, 4, 10,
                                 14, 29, 50, 57, 20, 35, 48, 21, 53, 39, 26, 40, 58}}

Definition at line 741 of file Water97FluidProperties.h.

Referenced by d2gamma2_dpi2(), d2gamma2_dpitau(), d2gamma2_dtau2(), dgamma2_dpi(), dgamma2_dtau(), and gamma2().

◆ _J3

const std::array<int, 40> Water97FluidProperties::_J3

protected

Initial value:

{{0, 0,  1,  2,  7,  10, 12, 23, 2,  6, 15, 17, 0,  2,
                                 6, 7,  22, 26, 0,  2,  4,  16, 26, 0, 2,  4,  26, 1,
                                 3, 26, 0,  2,  26, 2,  26, 2,  26, 0, 1,  26}}

Definition at line 819 of file Water97FluidProperties.h.

Referenced by d2phi3_ddelta2(), d2phi3_ddeltatau(), d2phi3_dtau2(), dphi3_ddelta(), dphi3_dtau(), and phi3().

◆ _J3s

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

protected

Definition at line 1154 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _J5

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

protected

Definition at line 1262 of file Water97FluidProperties.h.

Referenced by d2gamma5_dpi2(), d2gamma5_dpitau(), d2gamma5_dtau2(), dgamma5_dpi(), dgamma5_dtau(), and gamma5().

◆ _Jph1

const std::array<int, 20> Water97FluidProperties::_Jph1

protected

Initial value:

{

{0, 1, 2, 6, 22, 32, 0, 1, 2, 3, 4, 10, 32, 10, 32, 10, 32, 32, 32, 32}}

Definition at line 708 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _Jph2a

const std::array<int, 36> Water97FluidProperties::_Jph2a

protected

Initial value:

{{0,  1,  2,  3,  7,  20, 0,  1,  2,  3,  7,  9,
                                    11, 18, 44, 0,  2,  7,  36, 38, 40, 42, 44, 24,
                                    44, 12, 32, 44, 32, 36, 42, 34, 44, 28}}

Definition at line 759 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _Jph2b

const std::array<int, 38> Water97FluidProperties::_Jph2b

protected

Initial value:

{{0,  1,  2,  12, 18, 24, 28, 40, 0, 2,  6,  12, 18,
                                    24, 28, 40, 2,  8,  18, 40, 1,  2, 12, 24, 2,  12,
                                    18, 24, 28, 40, 18, 24, 40, 28, 2, 28, 1,  40}}

Definition at line 778 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _Jph2c

const std::array<int, 23> Water97FluidProperties::_Jph2c

protected

Initial value:

{

{0, 4, 0, 2, 0, 2, 0, 1, 0, 2, 0, 1, 4, 8, 4, 0, 1, 4, 10, 12, 16, 20, 22}}

Definition at line 793 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _Jph3a

const std::array<int, 31> Water97FluidProperties::_Jph3a

protected

Initial value:

{{0, 1, 2, 6, 14, 16, 20, 22, 1, 5, 12, 0, 2, 4, 10, 2,

0, 1, 3, 4, 0, 2, 0, 1, 1, 0, 1, 0, 3, 4, 5}}

Definition at line 1219 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _Jph3b

const std::array<int, 33> Water97FluidProperties::_Jph3b

protected

Initial value:

{{0, 1, 0, 1, 5, 10, 12, 0, 1, 2, 4, 10, 0, 1, 2, 0, 1,

5, 0, 4, 2, 4, 6, 10, 14, 16, 0, 2, 1, 1, 1, 1, 1}}

Definition at line 1237 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _k_a

std::array<Real, 4> Water97FluidProperties::_k_a {{0.0102811, 0.0299621, 0.0156146, -0.00422464}}

protected

Constants for thermal conductivity.

Definition at line 1331 of file Water97FluidProperties.h.

Referenced by k_from_rho_T().

◆ _Mh2o

const Real Water97FluidProperties::_Mh2o

protected

Water molar mass (kg/mol)

Definition at line 658 of file Water97FluidProperties.h.

Referenced by molarMass().

◆ _mu_H0

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

protected

Definition at line 1324 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T().

◆ _mu_H1

const std::array<Real, 21> Water97FluidProperties::_mu_H1

protected

Initial value:

{
      {5.20094e-1, 8.50895e-2, -1.08374,    -2.89555e-1, 2.22531e-1,  9.99115e-1,  1.88797,
       1.26613,    1.20573e-1, -2.81378e-1, -9.06851e-1, -7.72479e-1, -4.89837e-1, -2.57040e-1,
       1.61913e-1, 2.57399e-1, -3.25372e-2, 6.98452e-2,  8.72102e-3,  -4.35673e-3, -5.93264e-4}}

Definition at line 1325 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T().

◆ _mu_I

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

protected

Constants from Release on the IAPWS Formulation 2008 for the Viscosity of Ordinary Water Substance.

Definition at line 1322 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T().

◆ _mu_J

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

protected

Definition at line 1323 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T().

◆ _n02

const std::array<Real, 9> Water97FluidProperties::_n02

protected

Initial value:

{{-0.96927686500217e1,
                                  0.10086655968018e2,
                                  -0.56087911283020e-2,
                                  0.71452738081455e-1,
                                  -0.40710498223928e0,
                                  0.14240819171444e1,
                                  -0.43839511319450e1,
                                  -0.28408632460772e0,
                                  0.21268463753307e-1}}

Constants for region 2.

Definition at line 712 of file Water97FluidProperties.h.

Referenced by d2gamma2_dtau2(), dgamma2_dtau(), and gamma2().

◆ _n05

const std::array<Real, 6> Water97FluidProperties::_n05

protected

Initial value:

{{-0.13179983674201e2,
                                  0.68540841634434e1,
                                  -0.24805148933466e-1,
                                  0.36901534980333,
                                  -0.31161318213925e1,
                                  -0.32961626538917}}

Definition at line 1254 of file Water97FluidProperties.h.

Referenced by d2gamma5_dtau2(), dgamma5_dtau(), and gamma5().

◆ _n1

const std::array<Real, 34> Water97FluidProperties::_n1

protected

Initial value:

{
      {0.14632971213167e0,    -0.84548187169114e0,  -0.37563603672040e1,   0.33855169168385e1,
       -0.95791963387872e0,   0.15772038513228e0,   -0.16616417199501e-1,  0.81214629983568e-3,
       0.28319080123804e-3,   -0.60706301565874e-3, -0.18990068218419e-1,  -0.32529748770505e-1,
       -0.21841717175414e-1,  -0.52838357969930e-4, -0.47184321073267e-3,  -0.30001780793026e-3,
       0.47661393906987e-4,   -0.44141845330846e-5, -0.72694996297594e-15, -0.31679644845054e-4,
       -0.28270797985312e-5,  -0.85205128120103e-9, -0.22425281908000e-5,  -0.65171222895601e-6,
       -0.14341729937924e-12, -0.40516996860117e-6, -0.12734301741641e-8,  -0.17424871230634e-9,
       -0.68762131295531e-18, 0.14478307828521e-19, 0.26335781662795e-22,  -0.11947622640071e-22,
       0.18228094581404e-23,  -0.93537087292458e-25}}

Reference constants used in to calculate thermophysical properties of water.

Taken from Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam and 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 Constants for region 1

Definition at line 681 of file Water97FluidProperties.h.

Referenced by d2gamma1_dpi2(), d2gamma1_dpitau(), d2gamma1_dtau2(), dgamma1_dpi(), dgamma1_dtau(), and gamma1().

◆ _n2

const std::array<Real, 43> Water97FluidProperties::_n2

protected

Initial value:

{
      {-0.17731742473213e-2, -0.17834862292358e-1,  -0.45996013696365e-1,  -0.57581259083432e-1,
       -0.50325278727930e-1, -0.33032641670203e-4,  -0.18948987516315e-3,  -0.39392777243355e-2,
       -0.43797295650573e-1, -0.26674547914087e-4,  0.20481737692309e-7,   0.43870667284435e-6,
       -0.32277677238570e-4, -0.15033924542148e-2,  -0.40668253562649e-1,  -0.78847309559367e-9,
       0.12790717852285e-7,  0.48225372718507e-6,   0.22922076337661e-5,   -0.16714766451061e-10,
       -0.21171472321355e-2, -0.23895741934104e2,   -0.59059564324270e-17, -0.12621808899101e-5,
       -0.38946842435739e-1, 0.11256211360459e-10,  -0.82311340897998e1,   0.19809712802088e-7,
       0.10406965210174e-18, -0.10234747095929e-12, -0.10018179379511e-8,  -0.80882908646985e-10,
       0.10693031879409e0,   -0.33662250574171e0,   0.89185845355421e-24,  0.30629316876232e-12,
       -0.42002467698208e-5, -0.59056029685639e-25, 0.37826947613457e-5,   -0.12768608934681e-14,
       0.73087610595061e-28, 0.55414715350778e-16,  -0.94369707241210e-6}}

Definition at line 722 of file Water97FluidProperties.h.

Referenced by d2gamma2_dpi2(), d2gamma2_dpitau(), d2gamma2_dtau2(), dgamma2_dpi(), dgamma2_dtau(), and gamma2().

◆ _n23

const std::array<Real, 5> Water97FluidProperties::_n23

protected

Initial value:

{{0.34805185628969e3,
                                  -0.11671859879975e1,
                                  0.10192970039326e-2,
                                  0.57254459862746e3,
                                  0.13918839778870e2}}

Constants for the boundary between regions 2 and 3.

Definition at line 797 of file Water97FluidProperties.h.

Referenced by b23p(), and b23T().

◆ _n3

const std::array<Real, 40> Water97FluidProperties::_n3

protected

Initial value:

{
      {0.10658070028513e1,  -0.15732845290239e2,  0.20944396974307e2,   -0.76867707878716e1,
       0.26185947787954e1,  -0.28080781148620e1,  0.12053369696517e1,   -0.84566812812502e-2,
       -0.12654315477714e1, -0.11524407806681e1,  0.88521043984318e0,   -0.64207765181607e0,
       0.38493460186671e0,  -0.85214708824206e0,  0.48972281541877e1,   -0.30502617256965e1,
       0.39420536879154e-1, 0.12558408424308e0,   -0.27999329698710e0,  0.13899799569460e1,
       -0.20189915023570e1, -0.82147637173963e-2, -0.47596035734923e0,  0.43984074473500e-1,
       -0.44476435428739e0, 0.90572070719733e0,   0.70522450087967e0,   0.10770512626332e0,
       -0.32913623258954e0, -0.50871062041158e0,  -0.22175400873096e-1, 0.94260751665092e-1,
       0.16436278447961e0,  -0.13503372241348e-1, -0.14834345352472e-1, 0.57922953628084e-3,
       0.32308904703711e-2, 0.80964802996215e-4,  -0.16557679795037e-3, -0.44923899061815e-4}}

Constants for region 3.

Definition at line 804 of file Water97FluidProperties.h.

Referenced by d2phi3_ddelta2(), d2phi3_ddeltatau(), d2phi3_dtau2(), dphi3_ddelta(), dphi3_dtau(), and phi3().

◆ _n3s

const std::vector<std::vector<Real> > Water97FluidProperties::_n3s

protected

Constants for all 26 subregions in region 3.

Definition at line 855 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _n4

const std::array<Real, 10> Water97FluidProperties::_n4

protected

Initial value:

{{0.11670521452767e4,
                                  -0.72421316703206e6,
                                  -0.17073846940092e2,
                                  0.12020824702470e5,
                                  -0.32325550322333e7,
                                  0.14915108613530e2,
                                  -0.48232657361591e4,
                                  0.40511340542057e6,
                                  -0.238555575678490,
                                  0.65017534844798e3}}

Constants for region 4 (the saturation curve up to the critical point)

Definition at line 1241 of file Water97FluidProperties.h.

Referenced by vaporPressure(), and vaporTemperature_ad().

◆ _n5

const std::array<Real, 6> Water97FluidProperties::_n5

protected

Initial value:

{{0.15736404855259e-2,
                                 0.90153761673944e-3,
                                 -0.50270077677648e-2,
                                 0.22440037409485e-5,
                                 -0.41163275453471e-5,
                                 0.37919454822955e-7}}

Definition at line 1263 of file Water97FluidProperties.h.

Referenced by d2gamma5_dpi2(), d2gamma5_dpitau(), d2gamma5_dtau2(), dgamma5_dpi(), dgamma5_dtau(), and gamma5().

◆ _nph1

const std::array<Real, 20> Water97FluidProperties::_nph1

protected

Initial value:

{
      {-0.23872489924521e3,  0.40421188637945e3,    0.11349746881718e3,   -0.58457616048039e1,
       -0.15285482413140e-3, -0.10866707695377e-5,  -0.13391744872602e2,  0.43211039183559e2,
       -0.54010067170506e2,  0.30535892203916e2,    -0.65964749423638e1,  0.93965400878363e-2,
       0.11573647505340e-6,  -0.25858641282073e-4,  -0.40644363084799e-8, 0.66456186191635e-7,
       0.80670734103027e-10, -0.93477771213947e-12, 0.58265442020601e-14, -0.15020185953503e-16}}

Definition at line 699 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _nph2a

const std::array<Real, 36> Water97FluidProperties::_nph2a

protected

Initial value:

{
      {0.10898952318288e4,  0.84951654495535e3,   -0.10781748091826e3, 0.33153654801263e2,
       -0.74232016790248e1, 0.11765048724356e2,   0.18445749355790e1,  -0.41792700549624e1,
       0.62478196935812e1,  -0.17344563108114e2,  -0.20058176862096e3, 0.27196065473796e3,
       -0.45511318285818e3, 0.30919688604755e4,   0.25226640357872e6,  -0.61707422868339e-2,
       -0.31078046629583,   0.11670873077107e2,   0.12812798404046e9,  -0.98554909623276e9,
       0.28224546973002e10, -0.35948971410703e10, 0.17227349913197e10, -0.13551334240775e5,
       0.12848734664650e8,  0.13865724283226e1,   0.23598832556514e6,  -0.13105236545054e8,
       0.73999835474766e4,  -0.55196697030060e6,  0.37154085996233e7,  0.19127729239660e5,
       -0.41535164835634e6, -0.62459855192507e2}}

Definition at line 745 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _nph2b

const std::array<Real, 38> Water97FluidProperties::_nph2b

protected

Initial value:

{
      {0.14895041079516e4,    0.74307798314034e3,   -0.97708318797837e2,   0.24742464705674e1,
       -0.63281320016026,     0.11385952129658e1,   -0.47811863648625,     0.85208123431544e-2,
       0.93747147377932,      0.33593118604916e1,   0.33809355601454e1,    0.16844539671904,
       0.73875745236695,      -0.47128737436186,    0.15020273139707,      -0.21764114219750e-2,
       -0.21810755324761e-1,  -0.10829784403677,    -0.46333324635812e-1,  0.71280351959551e-4,
       0.11032831789999e-3,   0.18955248387902e-3,  0.30891541160537e-2,   0.13555504554949e-2,
       0.28640237477456e-6,   -0.10779857357512e-4, -0.76462712454814e-4,  0.14052392818316e-4,
       -0.31083814331434e-4,  -0.10302738212103e-5, 0.28217281635040e-6,   0.12704902271945e-5,
       0.73803353468292e-7,   -0.11030139238909e-7, -0.81456365207833e-13, -0.25180545682962e-10,
       -0.17565233969407e-17, 0.86934156344163e-14}}

Definition at line 763 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _nph2c

const std::array<Real, 23> Water97FluidProperties::_nph2c

protected

Initial value:

{
      {-0.32368398555242e13, 0.73263350902181e13,   0.35825089945447e12,  -0.58340131851590e12,
       -0.10783068217470e11, 0.20825544563171e11,   0.61074783564516e6,   0.85977722535580e6,
       -0.25745723604170e5,  0.31081088422714e5,    0.12082315865936e4,   0.48219755109255e3,
       0.37966001272486e1,   -0.10842984880077e2,   -0.45364172676660e-1, 0.14559115658698e-12,
       0.11261597407230e-11, -0.17804982240686e-10, 0.12324579690832e-6,  -0.11606921130984e-5,
       0.27846367088554e-4,  -0.59270038474176e-3,  0.12918582991878e-2}}

Definition at line 782 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _nph3a

const std::array<Real, 31> Water97FluidProperties::_nph3a

protected

Initial value:

{
      {-0.133645667811215e-6, 0.455912656802978e-5,  -0.146294640700979e-4, 0.639341312970080e-2,
       0.372783927268847e3,   -0.718654377460447e4,  0.573494752103400e6,   -0.267569329111439e7,
       -0.334066283302614e-4, -0.245479214069597e-1, 0.478087847764996e2,   0.764664131818904e-5,
       0.128350627676972e-2,  0.171219081377331e-1,  -0.851007304583213e1,  -0.136513461629781e-1,
       -0.384460997596657e-5, 0.337423807911655e-2,  -0.551624873066791,    0.729202277107470,
       -0.992522757376041e-2, -.119308831407288,     .793929190615421,      .454270731799386,
       .20999859125991,       -0.642109823904738e-2, -0.235155868604540e-1, 0.252233108341612e-2,
       -0.764885133368119e-2, 0.136176427574291e-1,  -0.133027883575669e-1}}

Definition at line 1205 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _nph3b

const std::array<Real, 33> Water97FluidProperties::_nph3b

protected

Initial value:

{
      {0.323254573644920e-4,  -0.127575556587181e-3, -0.475851877356068e-3, 0.156183014181602e-2,
       0.105724860113781,     -0.858514221132534e2,  0.724140095480911e3,   0.296475810273257e-2,
       -0.592721983365988e-2, -0.126305422818666e-1, -0.115716196364853,    0.849000969739595e2,
       -0.108602260086615e-1, 0.154304475328851e-1,  0.750455441524466e-1,  0.252520973612982e-1,
       -0.602507901232996e-1, -0.307622221350501e1,  -0.574011959864879e-1, 0.503471360939849e1,
       -0.925081888584834,    0.391733882917546e1,   -0.773146007130190e2,  0.949308762098587e4,
       -0.141043719679409e7,  0.849166230819026e7,   0.861095729446704,     0.323346442811720,
       0.873281936020439,     -0.436653048526683,    0.286596714529479,     -0.131778331276228,
       0.676682064330275e-2}}

Definition at line 1222 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _p_critical

const Real Water97FluidProperties::_p_critical

protected

Critical pressure (Pa)

Definition at line 662 of file Water97FluidProperties.h.

Referenced by criticalPressure(), tempXY(), and vaporTemperature_ad().

◆ _p_star

const std::array<Real, 5> Water97FluidProperties::_p_star {{16.53e6, 1.0e6, 1.0e6, 1.0e6, 1.0e6}}

protected

Pressure scale for each region.

Definition at line 1336 of file Water97FluidProperties.h.

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().

◆ _p_triple

const Real Water97FluidProperties::_p_triple

protected

Triple point pressure (Pa)

Definition at line 668 of file Water97FluidProperties.h.

Referenced by triplePointPressure().

◆ _par3

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

protected

Initial value:

{
      {{{0.0024, 100.0, 760.0, 0.085, 0.817, 1.0, 1.0, 1.0}},
       {{0.0041, 100.0, 860.0, 0.280, 0.779, 1.0, 1.0, 1.0}},
       {{0.0022, 40.0, 690.0, 0.259, 0.903, 1.0, 1.0, 1.0}},
       {{0.0029, 40.0, 690.0, 0.559, 0.939, 1.0, 1.0, 4.0}},
       {{0.0032, 40.0, 710.0, 0.587, 0.918, 1.0, 1.0, 1.0}},
       {{0.0064, 40.0, 730.0, 0.587, 0.891, 0.5, 1.0, 4.0}},
       {{0.0027, 25.0, 660.0, 0.872, 0.971, 1.0, 1.0, 4.0}},
       {{0.0032, 25.0, 660.0, 0.898, 0.983, 1.0, 1.0, 4.0}},
       {{0.0041, 25.0, 660.0, 0.910, 0.984, 0.5, 1.0, 4.0}},
       {{0.0054, 25.0, 670.0, 0.875, 0.964, 0.5, 1.0, 4.0}},
       {{0.0077, 25.0, 680.0, 0.802, 0.935, 1.0, 1.0, 1.0}},
       {{0.0026, 24.0, 650.0, 0.908, 0.989, 1.0, 1.0, 4.0}},
       {{0.0028, 23.0, 650.0, 1.000, 0.997, 1.0, 0.25, 1.0}},
       {{0.0031, 23.0, 650.0, 0.976, 0.997, 0.0, 0.0, 0.0}},
       {{0.0034, 23.0, 650.0, 0.974, 0.996, 0.5, 1.0, 1.0}},
       {{0.0041, 23.0, 650.0, 0.972, 0.997, 0.5, 1.0, 1.0}},
       {{0.0022, 23.0, 650.0, 0.848, 0.983, 1.0, 1.0, 4.0}},
       {{0.0054, 23.0, 650.0, 0.874, 0.982, 1.0, 1.0, 1.0}},
       {{0.0022, 21.0, 640.0, 0.886, 0.990, 1.0, 1.0, 4.0}},
       {{0.0088, 20.0, 650.0, 0.803, 1.020, 1.0, 1.0, 1.0}},
       {{0.0026, 23.0, 650.0, 0.902, 0.988, 1.0, 1.0, 1.0}},
       {{0.0031, 23.0, 650.0, 0.960, 0.995, 1.0, 1.0, 1.0}},
       {{0.0039, 23.0, 650.0, 0.959, 0.995, 1.0, 1.0, 4.0}},
       {{0.0049, 23.0, 650.0, 0.910, 0.988, 1.0, 1.0, 1.0}},
       {{0.0031, 22.0, 650.0, 0.996, 0.994, 1.0, 1.0, 4.0}},
       {{0.0038, 22.0, 650.0, 0.993, 0.994, 1.0, 1.0, 4.0}}}}

Definition at line 823 of file Water97FluidProperties.h.

Referenced by densityRegion3().

◆ _par3N

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

protected

Initial value:

{{30, 32, 35, 38, 29, 42, 38, 29, 42, 29, 34, 43, 40,

39, 24, 27, 24, 27, 29, 33, 38, 39, 35, 36, 20, 23}}

Definition at line 851 of file Water97FluidProperties.h.

Referenced by densityRegion3().

◆ _R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 42 of file FluidProperties.h.

Referenced by HelmholtzFluidProperties::c_from_p_T(), IdealGasFluidProperties::c_from_p_T(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), HelmholtzFluidProperties::p_from_rho_T(), HelmholtzFluidProperties::rho_from_p_T(), IdealGasFluidProperties::rho_from_p_T(), and HelmholtzFluidProperties::s_from_p_T().

◆ _rho_critical

const Real Water97FluidProperties::_rho_critical

protected

Critical density (kg/m^3)

Definition at line 666 of file Water97FluidProperties.h.

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

◆ _Rw

const Real Water97FluidProperties::_Rw

protected

Specific gas constant for H2O (universal gas constant / molar mass of water - J/kg/K)

Definition at line 660 of file Water97FluidProperties.h.

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().

◆ _T_c2k

const Real FluidProperties::_T_c2k

protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 46 of file FluidProperties.h.

Referenced by NaClFluidProperties::cp_from_p_T(), BrineFluidProperties::cp_from_p_T_X(), NaClFluidProperties::h_from_p_T(), BrineFluidProperties::h_from_p_T_X(), BrineFluidProperties::haliteSolubility(), NaClFluidProperties::k_from_p_T(), BrineFluidProperties::k_from_p_T_X(), BrineFluidProperties::mu_from_p_T_X(), CO2FluidProperties::partialDensity(), NaClFluidProperties::rho_from_p_T(), and BrineFluidProperties::rho_from_p_T_X().

◆ _T_critical

const Real Water97FluidProperties::_T_critical

protected

Critical temperature (K)

Definition at line 664 of file Water97FluidProperties.h.

Referenced by criticalTemperature(), mu_from_rho_T(), tempXY(), and vaporPressure().

◆ _T_star

const std::array<Real, 5> Water97FluidProperties::_T_star {{1386.0, 540.0, _T_critical, 1.0, 1000.0}}

protected

Temperature scale for each region.

Definition at line 1334 of file Water97FluidProperties.h.

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().

◆ _T_triple

const Real Water97FluidProperties::_T_triple

protected

Triple point temperature (K)

Definition at line 670 of file Water97FluidProperties.h.

Referenced by triplePointTemperature().

◆ _tempXY_I

const std::vector<std::vector<int> > Water97FluidProperties::_tempXY_I

protected

Initial value:

{{0, 1, 2, -1, -2},
                                                {0, 1, 2, 3},
                                                {0, 1, 2, 3, 4},
                                                {0, 1, 2, 3, 4},
                                                {0, 1, 2, 3, 4},
                                                {0, 1, 2, 3},
                                                {0, 1, 2, -1, -2},
                                                {0, 1, 2, 3},
                                                {0, 1, 2, 3},
                                                {0, 1, 2, 3, 4},
                                                {0, 1, 2, -1, -2}}

Constnats for the tempXY() method.

Definition at line 1271 of file Water97FluidProperties.h.

Referenced by tempXY().

◆ _tempXY_n

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

protected

Definition at line 1283 of file Water97FluidProperties.h.

Referenced by tempXY().

◆ h [1/2]

e e e e SinglePhaseFluidProperties::h

inherited

Definition at line 163 of file SinglePhaseFluidProperties.h.

◆ h [2/2]

e e e e p h T T T T T v v v SinglePhaseFluidProperties::h

inherited

Definition at line 183 of file SinglePhaseFluidProperties.h.

◆ p [1/6]

e e e e p h SinglePhaseFluidProperties::p

inherited

Definition at line 167 of file SinglePhaseFluidProperties.h.

◆ p [2/6]

e e e e p h T T SinglePhaseFluidProperties::p

inherited

Definition at line 171 of file SinglePhaseFluidProperties.h.

◆ p [3/6]

e e e e p h T T T SinglePhaseFluidProperties::p

inherited

Definition at line 173 of file SinglePhaseFluidProperties.h.

◆ p [4/6]

e e e e p h T T T T SinglePhaseFluidProperties::p

inherited

Definition at line 175 of file SinglePhaseFluidProperties.h.

◆ p [5/6]

e e e e p h T T T T T v v v s SinglePhaseFluidProperties::p

inherited

Definition at line 185 of file SinglePhaseFluidProperties.h.

◆ p [6/6]

e e e e p h T T T T T v v v s h SinglePhaseFluidProperties::p

inherited

Definition at line 187 of file SinglePhaseFluidProperties.h.

◆ propfuncWithDefault

e e e e p h T T T T T v v v s h T e SinglePhaseFluidProperties::propfuncWithDefault(gamma, p, T)

inherited

Definition at line 190 of file SinglePhaseFluidProperties.h.

◆ rho

e e e e p h T SinglePhaseFluidProperties::rho

inherited

Definition at line 169 of file SinglePhaseFluidProperties.h.

◆ T [1/3]

e e e e p h T T T T T SinglePhaseFluidProperties::T

inherited

Definition at line 177 of file SinglePhaseFluidProperties.h.

◆ T [2/3]

e e e e p h T T T T T v SinglePhaseFluidProperties::T

inherited

Definition at line 179 of file SinglePhaseFluidProperties.h.

◆ T [3/3]

e e e e p h T T T T T v v SinglePhaseFluidProperties::T

inherited

Definition at line 181 of file SinglePhaseFluidProperties.h.

◆ v [1/6]

SinglePhaseFluidProperties::v

inherited

Definition at line 155 of file SinglePhaseFluidProperties.h.

◆ v [2/6]

e SinglePhaseFluidProperties::v

inherited

Definition at line 157 of file SinglePhaseFluidProperties.h.

◆ v [3/6]

e e SinglePhaseFluidProperties::v

inherited

Definition at line 159 of file SinglePhaseFluidProperties.h.

◆ v [4/6]

e e e SinglePhaseFluidProperties::v

inherited

Definition at line 161 of file SinglePhaseFluidProperties.h.

◆ v [5/6]

e e e e p SinglePhaseFluidProperties::v

inherited

Definition at line 165 of file SinglePhaseFluidProperties.h.

◆ v [6/6]

e e e e p h T T T T T v v v s h T SinglePhaseFluidProperties::v

inherited

Definition at line 189 of file SinglePhaseFluidProperties.h.

The documentation for this class was generated from the following files:

fluid_properties/include/userobjects/Water97FluidProperties.h
fluid_properties/src/userobjects/Water97FluidProperties.C

Public Member Functions

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Enumeration Documentation

◆ subregionEnum

Constructor & Destructor Documentation

◆ Water97FluidProperties()

◆ ~Water97FluidProperties()

Member Function Documentation

◆ b23p()

◆ b23T()

◆ b2bc()

◆ b3ab()

◆ 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_from_p_T() [1/2]

◆ e_from_p_T() [2/2]

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ fluidPropError()

◆ gamma1()

◆ gamma2()

◆ gamma5()

◆ h_from_p_T() [1/2]

◆ h_from_p_T() [2/2]

◆ henryCoefficients()

◆ henryConstant() [1/3]

◆ henryConstant() [2/3]

◆ henryConstant() [3/3]

◆ initialize()

◆ inRegion()

◆ inRegionPH()

◆ k_from_p_T() [1/2]

◆ k_from_p_T() [2/2]

◆ k_from_rho_T()

◆ molarMass()

◆ mu_from_p_T() [1/2]

◆ mu_from_p_T() [2/2]

◆ mu_from_rho_T() [1/2]

◆ mu_from_rho_T() [2/2]

◆ phi3()

◆ propfunc() [1/16]

◆ propfunc() [2/16]

◆ propfunc() [3/16]

◆ propfunc() [4/16]

◆ propfunc() [5/16]

◆ propfunc() [6/16]

◆ propfunc() [7/16]