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

#include <Water97FluidProperties.h>

Inheritance diagram for Water97FluidProperties:

[legend]

Public Member Functions
	Water97FluidProperties (const InputParameters &parameters)

virtual	~Water97FluidProperties ()

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

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

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

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

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

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

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

virtual Real	rho_from_p_T (Real pressure, Real temperature) const override
	Density from pressure and temperature. More...

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

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

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

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

virtual Real	c_from_p_T (Real pressure, Real temperature) const override
	Speed of sound. More...

virtual Real	cp_from_p_T (Real pressure, Real temperature) const override
	Isobaric specific heat capacity. More...

virtual Real	cv_from_p_T (Real pressure, Real temperature) const override
	Isochoric specific heat. More...

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
	Dynamic viscosity as a function of density and temperature. More...

virtual void	mu_drhoT_from_rho_T (Real density, Real temperature, Real ddensity_dT, Real &mu, Real &dmu_drho, Real &dmu_dT) const override
	Dynamic viscosity and its derivatives wrt density and temperature. More...

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

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

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

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

virtual Real	k_from_rho_T (Real density, Real temperature) const override
	Thermal conductivity as a function of density and temperature. More...

virtual Real	s_from_p_T (Real pressure, Real temperature) const override
	Specific entropy from pressure and temperature. More...

virtual void	s_from_p_T (Real p, Real T, Real &s, Real &ds_dp, Real &ds_dT) const override
	Specific entropy and its derivatives from pressure and temperature. More...

virtual Real	h_from_p_T (Real pressure, Real temperature) const override
	Specific enthalpy from pressure and temperature. More...

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual Real	p_from_v_e (Real v, Real e) const
	Pressure from specific volume and specific internal energy. More...

virtual void	p_from_v_e (Real v, Real e, Real &p, Real &dp_dv, Real &dp_de) const
	Pressure and its derivatives from specific volume and specific internal energy. More...

virtual Real	T_from_v_e (Real v, Real e) const
	Temperature from specific volume and specific internal energy. More...

virtual void	T_from_v_e (Real v, Real e, Real &T, Real &dT_dv, Real &dT_de) const
	Temperature and its derivatives from specific volume and specific internal energy. More...

virtual Real	c_from_v_e (Real v, Real e) const
	Sound speed from specific volume and specific internal energy. More...

virtual void	c_from_v_e (Real v, Real e, Real &c, Real &dc_dv, Real &dc_de) const
	Sound speed and its derivatives from specific volume and specific internal energy. More...

virtual Real	cp_from_v_e (Real v, Real e) const
	Isobaric (constant-pressure) specific heat from specific volume and specific internal energy. More...

virtual Real	cv_from_v_e (Real v, Real e) const
	Isochoric (constant-volume) specific heat from specific volume and specific internal energy. More...

virtual Real	mu_from_v_e (Real v, Real e) const
	Dynamic viscosity from specific volume and specific internal energy. More...

virtual Real	k_from_v_e (Real v, Real e) const
	Thermal conductivity from specific volume and specific internal energy. More...

virtual Real	s_from_v_e (Real v, Real e) const
	Specific entropy from specific volume and specific internal energy. More...

virtual void	s_from_v_e (Real v, Real e, Real &s, Real &ds_dv, Real &ds_de) const
	Specific entropy and its derivatives from specific volume and specific internal energy. More...

virtual Real	s (Real pressure, Real temperature) const

virtual Real	s_from_h_p (Real h, Real p) const
	Specific entropy from specific enthalpy and pressure. More...

virtual void	s_from_h_p (Real h, Real p, Real &s, Real &ds_dh, Real &ds_dp) const
	Specific entropy and its derivatives from specific enthalpy and pressure. More...

virtual Real	rho_from_p_s (Real p, Real s) const
	Density from pressure and specific entropy. More...

virtual void	rho_from_p_s (Real p, Real s, Real &rho, Real &drho_dp, Real &drho_ds) const
	Density and its derivatives from pressure and specific entropy. More...

virtual Real	e_from_v_h (Real v, Real h) const
	Specific internal energy as a function of specific volume and specific enthalpy. More...

virtual void	e_from_v_h (Real v, Real h, Real &e, Real &de_dv, Real &de_dh) const
	Specific internal energy and derivatives as a function of specific volume and specific enthalpy. More...

virtual Real	rho (Real p, Real T) const

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

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

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

virtual Real	e_from_p_rho (Real p, Real rho) const
	Specific internal energy from pressure and density. More...

virtual void	e_from_p_rho (Real p, Real rho, Real &e, Real &de_dp, Real &de_drho) const
	Specific internal energy and its derivatives from pressure and density. More...

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

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

virtual Real	e_from_T_v (Real T, Real v) const
	Specific volume and specific internal energy from temerature at the vapor spinodal. More...

virtual void	e_from_T_v (Real T, Real v, Real &e, Real &de_dT, Real &de_dv) const
	Specific internal energy and its derivatives from temperature and specific volume. More...

virtual Real	p_from_T_v (Real T, Real v) const
	Pressure from temperature and specific volume. More...

virtual void	p_from_T_v (Real T, Real v, Real &p, Real &dp_dT, Real &dp_dv) const
	Pressure and its derivatives from temperature and specific volume. More...

virtual Real	h_from_T_v (Real T, Real v) const
	Specific enthalpy from temperature and specific volume. More...

virtual void	h_from_T_v (Real T, Real v, Real &h, Real &dh_dT, Real &dh_dv) const
	Specific enthalpy and its derivatives from temperature and specific volume. More...

virtual Real	s_from_T_v (Real T, Real v) const
	Specific entropy from temperature and specific volume. More...

virtual void	s_from_T_v (Real T, Real v, Real &s, Real &ds_dT, Real &ds_dv) const
	Specific entropy and its derivatives from temperature and specific volume. More...

virtual Real	cv_from_T_v (Real T, Real v) const
	Specific isochoric heat capacity from temperature and specific volume. More...

virtual Real	h (Real p, Real T) const

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

virtual Real	e (Real pressure, Real temperature) const

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

virtual Real	p_from_h_s (Real h, Real s) const
	Pressure from specific enthalpy and specific entropy. More...

virtual void	p_from_h_s (Real h, Real s, Real &p, Real &dp_dh, Real &dp_ds) const
	Pressure and its derivatives from specific enthalpy and specific entropy. More...

virtual Real	g_from_v_e (Real v, Real e) const
	Gibbs free energy from specific volume and specific internal energy. More...

virtual Real	beta_from_p_T (Real p, Real T) const
	Thermal expansion coefficient from pressure and temperature. More...

virtual Real	beta (Real pressure, Real temperature) const

virtual Real	pp_sat_from_p_T (Real p, Real T) const
	Partial pressure at saturation in a gas mixture. More...

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

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

virtual Real	c (Real pressure, Real temperature) const

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

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

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

virtual Real	k (Real pressure, Real temperature) const

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

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

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

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

virtual void	henryConstantIAPWS_dT (Real temperature, Real &Kh, Real &dKh_dT, Real A, Real B, Real C) const
	IAPWS formulation of Henry's law constant for dissolution in water and derivative wrt temperature. 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	_R
	Universal gas constant (J/mol/K) More...

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

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 40 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 485 of file Water97FluidProperties.h.

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

Constructor & Destructor Documentation

◆ Water97FluidProperties()

Water97FluidProperties::Water97FluidProperties ( const InputParameters & parameters )

Definition at line 26 of file Water97FluidProperties.C.

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

◆ ~Water97FluidProperties()

Water97FluidProperties::~Water97FluidProperties ( )

virtual

Definition at line 38 of file Water97FluidProperties.C.

38 {}

Member Function Documentation

◆ b23p()

Real Water97FluidProperties::b23p ( Real temperature ) const

Auxillary equation for the boundary between regions 2 and 3.

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

Parameters

temperature water temperature (K)

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

Definition at line 874 of file Water97FluidProperties.C.

Referenced by b3ab(), and inRegion().

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

◆ b23T()

Real Water97FluidProperties::b23T ( Real pressure ) const

Auxillary equation for the boundary between regions 2 and 3.

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

Parameters

pressure water pressure (Pa)

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

Definition at line 886 of file Water97FluidProperties.C.

Referenced by inRegionPH().

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

◆ b2bc()

Real Water97FluidProperties::b2bc ( Real pressure ) const

Boundary between subregions b and c in region 2.

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

Parameters

pressure water pressure (Pa)

Returns: enthalpy at boundary (J/kg)

Definition at line 1657 of file Water97FluidProperties.C.

Referenced by subregion2ph().

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

◆ b3ab()

Real Water97FluidProperties::b3ab ( Real pressure ) const

Boundary between subregions a and b in region 3.

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

Parameters

pressure water pressure (Pa)

Returns: enthalpy at boundary (J/kg)

Definition at line 1773 of file Water97FluidProperties.C.

Referenced by subregion3ph().

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

◆ beta()

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

virtualinherited

Definition at line 256 of file SinglePhaseFluidProperties.C.

Referenced by vaporTemperature().

 {
   return beta_from_p_T(pressure, temperature);
 }

◆ beta_from_p_T()

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

virtualinherited

Thermal expansion coefficient from pressure and temperature.

Parameters

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

Returns: beta (1/K)

Reimplemented in SimpleFluidProperties.

Definition at line 171 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::beta().

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

◆ c()

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

virtualinherited

Definition at line 391 of file SinglePhaseFluidProperties.C.

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

 {
   return c_from_p_T(pressure, temperature);
 }

◆ c_from_p_T()

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

overridevirtual

Speed of sound.

Parameters

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

Returns: speed of sound (m/s)

Reimplemented from SinglePhaseFluidProperties.

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

◆ c_from_v_e() [1/2]

Real SinglePhaseFluidProperties::c_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Sound speed from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 50 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), and NSMachAux::computeValue().

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

◆ c_from_v_e() [2/2]

void SinglePhaseFluidProperties::c_from_v_e	(	Real	v,
		Real	e,
		Real &	c,
		Real &	dc_dv,
		Real &	dc_de
	)		const

virtualinherited

Sound speed and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[out]	dc_dv	derivative of sound speed w.r.t. specific volume
[out]	dc_de	derivative of sound speed w.r.t. specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 56 of file SinglePhaseFluidProperties.C.

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

◆ cp_from_p_T()

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

overridevirtual

Isobaric specific heat capacity.

Parameters

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

Returns: cp (J/kg/K)

Reimplemented from SinglePhaseFluidProperties.

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

◆ cp_from_v_e()

Real SinglePhaseFluidProperties::cp_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Isobaric (constant-pressure) specific heat from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 61 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), and GeneralVaporMixtureFluidProperties::cp_from_p_T().

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

◆ criticalDensity()

Real Water97FluidProperties::criticalDensity ( ) const

overridevirtual

Critical density.

Returns: critical density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 65 of file Water97FluidProperties.C.

 {
   return _rho_critical;
 }

◆ criticalInternalEnergy()

Real SinglePhaseFluidProperties::criticalInternalEnergy ( ) const

virtualinherited

Critical specific internal energy.

Returns: specific internal energy (J/kg)

Reimplemented in StiffenedGasFluidProperties.

Definition at line 232 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e().

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

◆ criticalPressure()

Real Water97FluidProperties::criticalPressure ( ) const

overridevirtual

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 53 of file Water97FluidProperties.C.

 {
   return _p_critical;
 }

◆ criticalTemperature()

Real Water97FluidProperties::criticalTemperature ( ) const

overridevirtual

Critical temperature.

Returns: critical temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 59 of file Water97FluidProperties.C.

 {
   return _T_critical;
 }

◆ cv_from_p_T()

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

overridevirtual

Isochoric specific heat.

Parameters

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

Returns: cv (J/kg/K)

Reimplemented from SinglePhaseFluidProperties.

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

◆ cv_from_T_v()

Real SinglePhaseFluidProperties::cv_from_T_v	(	Real	T,
		Real	v
	)		const

virtualinherited

Specific isochoric heat capacity from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 568 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::cv_from_p_T(), and IdealRealGasMixtureFluidProperties::cv_from_T_v().

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

◆ cv_from_v_e()

Real SinglePhaseFluidProperties::cv_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Isochoric (constant-volume) specific heat from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 66 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), and GeneralVaporMixtureFluidProperties::cv_from_p_T().

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

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

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

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

◆ d2gamma1_dpitau()

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

protected

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

Parameters

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

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

Definition at line 997 of file Water97FluidProperties.C.

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

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

◆ d2gamma1_dtau2()

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

protected

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

Parameters

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

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

Definition at line 986 of file Water97FluidProperties.C.

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

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

◆ d2gamma2_dpi2()

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

protected

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

Parameters

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

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

Definition at line 1040 of file Water97FluidProperties.C.

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

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

◆ d2gamma2_dpitau()

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

protected

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

Parameters

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

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

Definition at line 1088 of file Water97FluidProperties.C.

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

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

◆ d2gamma2_dtau2()

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

protected

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

Parameters

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

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

Definition at line 1071 of file Water97FluidProperties.C.

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

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

◆ d2gamma5_dpi2()

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

protected

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

Parameters

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

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

Definition at line 1198 of file Water97FluidProperties.C.

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

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

◆ d2gamma5_dpitau()

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

protected

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

Parameters

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

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

Definition at line 1246 of file Water97FluidProperties.C.

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

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

◆ d2gamma5_dtau2()

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

protected

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

Parameters

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

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

Definition at line 1229 of file Water97FluidProperties.C.

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

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

◆ d2phi3_ddelta2()

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

protected

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

Parameters

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

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

Definition at line 1123 of file Water97FluidProperties.C.

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

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

◆ d2phi3_ddeltatau()

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

protected

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

Parameters

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

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

Definition at line 1155 of file Water97FluidProperties.C.

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

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

◆ d2phi3_dtau2()

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

protected

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

Parameters

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

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

Definition at line 1144 of file Water97FluidProperties.C.

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

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

◆ densityRegion3()

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

Density function for Region 3 - supercritical water and steam.

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

Parameters

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

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

Definition at line 1524 of file Water97FluidProperties.C.

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

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

◆ dgamma1_dpi()

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

protected

Derivative of Gibbs free energy in Region 1 wrt pi.

Parameters

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

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

Definition at line 953 of file Water97FluidProperties.C.

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

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

◆ dgamma1_dtau()

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

protected

Derivative of Gibbs free energy in Region 1 wrt tau.

Parameters

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

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

Definition at line 975 of file Water97FluidProperties.C.

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

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

◆ dgamma2_dpi()

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

protected

Derivative of Gibbs free energy in Region 2 wrt pi.

Parameters

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

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

Definition at line 1026 of file Water97FluidProperties.C.

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

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

◆ dgamma2_dtau()

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

protected

Derivative of Gibbs free energy in Region 2 wrt tau.

Parameters

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

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

Definition at line 1055 of file Water97FluidProperties.C.

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

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

◆ dgamma5_dpi()

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

protected

Derivative of Gibbs free energy in Region 5 wrt pi.

Parameters

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

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

Definition at line 1184 of file Water97FluidProperties.C.

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

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

◆ dgamma5_dtau()

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

protected

Derivative of Gibbs free energy in Region 5 wrt tau.

Parameters

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

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

Definition at line 1213 of file Water97FluidProperties.C.

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

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

◆ dphi3_ddelta()

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

protected

Derivative of Helmholtz free energy in Region 3 wrt delta.

Parameters

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

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

Definition at line 1113 of file Water97FluidProperties.C.

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

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

◆ dphi3_dtau()

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

protected

Derivative of Helmholtz free energy in Region 3 wrt tau.

Parameters

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

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

Definition at line 1134 of file Water97FluidProperties.C.

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

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

◆ e()

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

virtualinherited

Definition at line 366 of file SinglePhaseFluidProperties.C.

 {
   return e_from_p_T(p, T);
 }

◆ e_dpT()

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

virtualinherited

Definition at line 372 of file SinglePhaseFluidProperties.C.

 {
   e_from_p_T(pressure, temperature, e, de_dp, de_dT);
 }

◆ e_from_p_rho() [1/2]

Real SinglePhaseFluidProperties::e_from_p_rho	(	Real	p,
		Real	rho
	)		const

virtualinherited

Specific internal energy from pressure and density.

Parameters

[in]	p	pressure
[in]	rho	density

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 502 of file SinglePhaseFluidProperties.C.

Referenced by StagnationTemperatureAux::computeValue(), InternalEnergyAux::computeValue(), and SinglePhaseFluidProperties::e_from_p_T().

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

◆ e_from_p_rho() [2/2]

void SinglePhaseFluidProperties::e_from_p_rho	(	Real	p,
		Real	rho,
		Real &	e,
		Real &	de_dp,
		Real &	de_drho
	)		const

virtualinherited

Specific internal energy and its derivatives from pressure and density.

Parameters

[in]	p	pressure
[in]	rho	density
[out]	e	specific internal energy
[out]	de_dp	derivative of specific internal energy w.r.t. pressure
[out]	de_drho	derivative of specific internal energy w.r.t. density

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 508 of file SinglePhaseFluidProperties.C.

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

◆ e_from_p_T() [1/2]

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

overridevirtual

Internal energy from pressure and temperature.

Parameters

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

Returns: internal energy (J/kg)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 188 of file Water97FluidProperties.C.

Referenced by e_from_p_T(), and rho_e_dpT().

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

◆ e_from_p_T() [2/2]

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

overridevirtual

Internal energy and its derivatives from pressure and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 235 of file Water97FluidProperties.C.

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

◆ e_from_T_v() [1/2]

Real SinglePhaseFluidProperties::e_from_T_v	(	Real	T,
		Real	v
	)		const

virtualinherited

Specific volume and specific internal energy from temerature at the vapor spinodal.

Parameters

[in] T temerature

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 513 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::e_from_T_v(), IdealRealGasMixtureFluidProperties::k_from_p_T(), IdealRealGasMixtureFluidProperties::k_from_T_v(), IdealRealGasMixtureFluidProperties::mu_from_p_T(), and IdealRealGasMixtureFluidProperties::mu_from_T_v().

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

◆ e_from_T_v() [2/2]

void SinglePhaseFluidProperties::e_from_T_v	(	Real	T,
		Real	v,
		Real &	e,
		Real &	de_dT,
		Real &	de_dv
	)		const

virtualinherited

Specific internal energy and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[out]	e	specific internal energy (J/kg)
[out]	de_dT	derivative of specific internal energy w.r.t. temperature
[out]	de_dv	derivative of specific internal energy w.r.t. specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 530 of file SinglePhaseFluidProperties.C.

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

◆ e_from_v_h() [1/2]

Real SinglePhaseFluidProperties::e_from_v_h	(	Real	v,
		Real	h
	)		const

virtualinherited

Specific internal energy as a function of specific volume and specific enthalpy.

Parameters

[in]	v	specific volume
[in]	h	specific enthalpy

Reimplemented in StiffenedGasFluidProperties, and IdealGasFluidProperties.

Definition at line 114 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::T_from_p_h().

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

◆ e_from_v_h() [2/2]

void SinglePhaseFluidProperties::e_from_v_h	(	Real	v,
		Real	h,
		Real &	e,
		Real &	de_dv,
		Real &	de_dh
	)		const

virtualinherited

Specific internal energy and derivatives as a function of specific volume and specific enthalpy.

Parameters

[in]	v	specific volume
[in]	h	specific enthalpy
[out]	de_dv	derivative of specific internal energy w.r.t. specific volume
[out]	de_dh	derivative of specific internal energy w.r.t. specific enthalpy

Reimplemented in StiffenedGasFluidProperties, and IdealGasFluidProperties.

Definition at line 120 of file SinglePhaseFluidProperties.C.

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

◆ e_spndl_from_v()

Real SinglePhaseFluidProperties::e_spndl_from_v ( Real v ) const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in StiffenedGasFluidProperties.

Definition at line 518 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::p_T_from_v_e().

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

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Definition at line 27 of file FluidProperties.h.

27 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 29 of file FluidProperties.h.

29 {}

◆ fluidName()

std::string Water97FluidProperties::fluidName ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name

Reimplemented from SinglePhaseFluidProperties.

Definition at line 41 of file Water97FluidProperties.C.

 {
   return "water";
 }

◆ g_from_v_e()

Real SinglePhaseFluidProperties::g_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Gibbs free energy from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in StiffenedGasFluidProperties, and IdealGasFluidProperties.

Definition at line 595 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties().

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

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

Referenced by s_from_p_T().

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

◆ gamma2()

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

protected

Gibbs free energy in Region 2 - superheated steam.

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

Parameters

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

Returns: Gibbs free energy (-)

Definition at line 1008 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

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

◆ gamma5()

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

protected

Gibbs free energy in Region 5.

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

Parameters

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

Returns: Gibbs free energy (-)

Definition at line 1166 of file Water97FluidProperties.C.

Referenced by s_from_p_T().

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

◆ gamma_from_p_T()

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

virtualinherited

Adiabatic index - ratio of specific heats.

Parameters

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

Returns: gamma (-)

Definition at line 250 of file SinglePhaseFluidProperties.C.

 {
   return cp_from_p_T(pressure, temperature) / cv_from_p_T(pressure, temperature);
 }

◆ h()

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

virtualinherited

Definition at line 474 of file SinglePhaseFluidProperties.C.

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

 {
   return h_from_p_T(p, T);
 }

◆ h_dpT()

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

virtualinherited

Definition at line 480 of file SinglePhaseFluidProperties.C.

 {
   h_from_p_T(p, T, h, dh_dp, dh_dT);
 }

◆ h_from_p_T() [1/2]

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

overridevirtual

Specific enthalpy from pressure and temperature.

Parameters

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

Returns: h (J/kg)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 695 of file Water97FluidProperties.C.

Referenced by inRegionPH().

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

◆ h_from_p_T() [2/2]

void Water97FluidProperties::h_from_p_T	(	Real	p,
		Real	T,
		Real &	h,
		Real &	dh_dp,
		Real &	dh_dT
	)		const

overridevirtual

Specific enthalpy and its derivatives from pressure and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

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

◆ h_from_T_v() [1/2]

Real SinglePhaseFluidProperties::h_from_T_v	(	Real	T,
		Real	v
	)		const

virtualinherited

Specific enthalpy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 546 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::cp_from_p_T(), and IdealRealGasMixtureFluidProperties::cp_from_T_v().

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

◆ h_from_T_v() [2/2]

void SinglePhaseFluidProperties::h_from_T_v	(	Real	T,
		Real	v,
		Real &	h,
		Real &	dh_dT,
		Real &	dh_dv
	)		const

virtualinherited

Specific enthalpy and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[out]	h	specific enthalpy (J/kg)
[out]	dh_dT	derivative of specific enthalpy w.r.t. temperature
[out]	dh_dv	derivative of specific enthalpy w.r.t. specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 552 of file SinglePhaseFluidProperties.C.

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

◆ henryConstant()

Real SinglePhaseFluidProperties::henryConstant ( Real temperature ) const

virtualinherited

Henry's law constant for dissolution in water.

Parameters

temperature fluid temperature (K)

Returns: Henry's constant

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

Definition at line 332 of file SinglePhaseFluidProperties.C.

Referenced by TabulatedFluidProperties::henryConstant().

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

◆ henryConstant_dT()

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

virtualinherited

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

Parameters

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

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

Definition at line 338 of file SinglePhaseFluidProperties.C.

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

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

◆ henryConstantIAPWS()

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

protectedvirtualinherited

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

Definition at line 262 of file SinglePhaseFluidProperties.C.

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

 {
   Real Tr = temperature / 647.096;
   Real tau = 1.0 - Tr;
 
   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
   std::vector<Real> a{-7.85951783, 1.84408259, -11.7866497, 22.6807411, -15.9618719, 1.80122502};
   std::vector<Real> b{1.0, 1.5, 3.0, 3.5, 4.0, 7.5};
   Real sum = 0.0;
 
   for (std::size_t i = 0; i < a.size(); ++i)
     sum += a[i] * std::pow(tau, b[i]);
 
   return 22.064e6 * std::exp(sum / Tr) * std::exp(lnkh);
 }

◆ henryConstantIAPWS_dT()

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

protectedvirtualinherited

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

Definition at line 281 of file SinglePhaseFluidProperties.C.

Referenced by MethaneFluidProperties::henryConstant_dT(), NitrogenFluidProperties::henryConstant_dT(), HydrogenFluidProperties::henryConstant_dT(), and CO2FluidProperties::henryConstant_dT().

 {
   Real pc = 22.064e6;
   Real Tc = 647.096;
 
   Real Tr = temperature / Tc;
   Real tau = 1.0 - Tr;
 
   Real lnkh = A / Tr + B * std::pow(tau, 0.355) / Tr + C * std::pow(Tr, -0.41) * std::exp(tau);
   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
   std::vector<Real> a{-7.85951783, 1.84408259, -11.7866497, 22.6807411, -15.9618719, 1.80122502};
   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);
   }
 
   Real p = pc * std::exp(sum / Tr);
   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 28 of file FluidProperties.h.

28 {}

◆ inRegion()

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

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

Parameters

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

Returns: region phase region index

Definition at line 897 of file Water97FluidProperties.C.

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

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

◆ inRegionPH()

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

protected

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

Parameters

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

Returns: region phase region index

Definition at line 1564 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ k()

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

virtualinherited

Definition at line 449 of file SinglePhaseFluidProperties.C.

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

 {
   return k_from_p_T(p, T);
 }

◆ k_dpT()

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

virtualinherited

Definition at line 461 of file SinglePhaseFluidProperties.C.

 {
   k_from_p_T(pressure, temperature, k, dk_dp, dk_dT);
 }

◆ k_from_p_T() [1/2]

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

overridevirtual

Thermal conductivity.

Parameters

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

Returns: thermal conductivity (W/m/K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 599 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

Thermal conductivity and its derivatives wrt pressure and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 606 of file Water97FluidProperties.C.

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

◆ k_from_rho_T()

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

overridevirtual

Thermal conductivity as a function of density and temperature.

Parameters

density	fluid density (kg/m^3)
temperature	fluid temperature (K)

Returns: thermal conductivity (W/m/K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 613 of file Water97FluidProperties.C.

Referenced by k_from_p_T().

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

◆ k_from_v_e()

Real SinglePhaseFluidProperties::k_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Thermal conductivity from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 76 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), GeneralVaporMixtureFluidProperties::k_from_p_T(), IdealRealGasMixtureFluidProperties::k_from_p_T(), and IdealRealGasMixtureFluidProperties::k_from_T_v().

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

◆ molarMass()

Real Water97FluidProperties::molarMass ( ) const

overridevirtual

Molar mass [kg/mol].

Returns: molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 47 of file Water97FluidProperties.C.

 {
   return _Mh2o;
 }

◆ mu()

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

virtualinherited

Dynamic viscosity.

Parameters

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

Returns: viscosity (Pa.s)

Definition at line 407 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::mu_dpT(), NitrogenFluidProperties::mu_drhoT_from_rho_T(), HydrogenFluidProperties::mu_drhoT_from_rho_T(), CO2FluidProperties::mu_drhoT_from_rho_T(), mu_drhoT_from_rho_T(), MethaneFluidProperties::mu_from_p_T(), CO2FluidProperties::mu_from_p_T(), NitrogenFluidProperties::mu_from_p_T(), HydrogenFluidProperties::mu_from_p_T(), IdealGasFluidPropertiesPT::mu_from_p_T(), mu_from_p_T(), SimpleFluidProperties::mu_from_p_T(), TabulatedFluidProperties::mu_from_p_T(), HydrogenFluidProperties::mu_from_rho_T(), HelmholtzFluidProperties::rho_mu(), IdealGasFluidPropertiesPT::rho_mu(), SimpleFluidProperties::rho_mu(), rho_mu(), TabulatedFluidProperties::rho_mu(), HelmholtzFluidProperties::rho_mu_dpT(), IdealGasFluidPropertiesPT::rho_mu_dpT(), SimpleFluidProperties::rho_mu_dpT(), rho_mu_dpT(), and TabulatedFluidProperties::rho_mu_dpT().

 {
   return mu_from_p_T(pressure, temperature);
 }

◆ mu_dpT()

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

virtualinherited

Dynamic viscosity and its derivatives wrt pressure and temperature.

Parameters

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

Definition at line 418 of file SinglePhaseFluidProperties.C.

 {
   mu_from_p_T(pressure, temperature, mu, dmu_dp, dmu_dT);
 }

◆ mu_drhoT_from_rho_T()

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

overridevirtual

Dynamic viscosity and its derivatives wrt density and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 528 of file Water97FluidProperties.C.

Referenced by mu_from_p_T(), and rho_mu_dpT().

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

◆ mu_from_p_T() [1/2]

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

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 484 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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 491 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_drhoT_from_rho_T(rho, temperature, drho_dT, mu, dmu_drho, dmu_dT);
   dmu_dp = dmu_drho * drho_dp;
 }

◆ mu_from_rho_T()

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

overridevirtual

Dynamic viscosity as a function of density and temperature.

Parameters

density	fluid density (kg/m^3)
temperature	fluid temperature (K)

Returns: viscosity (Pa.s)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 502 of file Water97FluidProperties.C.

Referenced by mu_from_p_T(), and rho_mu().

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

◆ mu_from_v_e()

Real SinglePhaseFluidProperties::mu_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Dynamic viscosity from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 71 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), GeneralVaporMixtureFluidProperties::mu_from_p_T(), IdealRealGasMixtureFluidProperties::mu_from_p_T(), and IdealRealGasMixtureFluidProperties::mu_from_T_v().

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

◆ p_from_h_s() [1/2]

Real SinglePhaseFluidProperties::p_from_h_s	(	Real	h,
		Real	s
	)		const

virtualinherited

Pressure from specific enthalpy and specific entropy.

Parameters

[in]	h	specific enthalpy
[in]	s	specific entropy

Reimplemented in StiffenedGasFluidProperties, and IdealGasFluidProperties.

Definition at line 584 of file SinglePhaseFluidProperties.C.

Referenced by StagnationPressureAux::computeValue(), and StagnationTemperatureAux::computeValue().

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

◆ p_from_h_s() [2/2]

void SinglePhaseFluidProperties::p_from_h_s	(	Real	h,
		Real	s,
		Real &	p,
		Real &	dp_dh,
		Real &	dp_ds
	)		const

virtualinherited

Pressure and its derivatives from specific enthalpy and specific entropy.

Parameters

[in]	h	specific enthalpy
[in]	s	specific entropy

Reimplemented in StiffenedGasFluidProperties, and IdealGasFluidProperties.

Definition at line 590 of file SinglePhaseFluidProperties.C.

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

◆ p_from_T_v() [1/2]

Real SinglePhaseFluidProperties::p_from_T_v	(	Real	T,
		Real	v
	)		const

virtualinherited

Pressure from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 535 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::p_from_T_v().

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

◆ p_from_T_v() [2/2]

void SinglePhaseFluidProperties::p_from_T_v	(	Real	T,
		Real	v,
		Real &	p,
		Real &	dp_dT,
		Real &	dp_dv
	)		const

virtualinherited

Pressure and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[out]	p	pressure (Pa)
[out]	dp_dT	derivative of pressure w.r.t. temperature
[out]	dp_dv	derivative of pressure w.r.t. specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 541 of file SinglePhaseFluidProperties.C.

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

◆ p_from_v_e() [1/2]

Real SinglePhaseFluidProperties::p_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Pressure from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, StiffenedGasFluidProperties, and NaNInterfaceTestFluidProperties.

Definition at line 28 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), PressureAux::computeValue(), StagnationPressureAux::computeValue(), and StagnationTemperatureAux::computeValue().

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

◆ p_from_v_e() [2/2]

void SinglePhaseFluidProperties::p_from_v_e	(	Real	v,
		Real	e,
		Real &	p,
		Real &	dp_dv,
		Real &	dp_de
	)		const

virtualinherited

Pressure and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[out]	p	pressure
[out]	dp_dv	derivative of pressure w.r.t. specific volume
[out]	dp_de	derivative of pressure w.r.t. specific internal energy

Reimplemented in IdealGasFluidProperties, StiffenedGasFluidProperties, and NaNInterfaceTestFluidProperties.

Definition at line 34 of file SinglePhaseFluidProperties.C.

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

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

Referenced by s_from_p_T().

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

◆ pp_sat_from_p_T()

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

virtualinherited

Partial pressure at saturation in a gas mixture.

Parameters

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

Returns: pp_sat (Pa)

Reimplemented in StiffenedGasFluidProperties.

Definition at line 196 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::xs_prim_from_p_T().

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

◆ rho()

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

virtualinherited

Definition at line 491 of file SinglePhaseFluidProperties.C.

 {
   return rho_from_p_T(p, T);
 }

◆ rho_dpT()

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

virtualinherited

Definition at line 359 of file SinglePhaseFluidProperties.C.

 {
   rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
 }

◆ rho_e_dpT()

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

overridevirtual

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

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 311 of file Water97FluidProperties.C.

 {
   rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
   e_from_p_T(pressure, temperature, e, de_dp, de_dT);
 }

◆ rho_from_p_s() [1/2]

Real SinglePhaseFluidProperties::rho_from_p_s	(	Real	p,
		Real	s
	)		const

virtualinherited

Density from pressure and specific entropy.

Parameters

[in]	p	pressure
[in]	s	specific entropy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 103 of file SinglePhaseFluidProperties.C.

Referenced by StagnationTemperatureAux::computeValue(), and SinglePhaseFluidProperties::T_from_p_h().

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

◆ rho_from_p_s() [2/2]

void SinglePhaseFluidProperties::rho_from_p_s	(	Real	p,
		Real	s,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_ds
	)		const

virtualinherited

Density and its derivatives from pressure and specific entropy.

Parameters

[in]	p	pressure
[in]	s	specific entropy
[out]	rho	density
[out]	drho_dp	derivative of density w.r.t. pressure
[out]	drho_ds	derivative of density w.r.t. specific entropy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 109 of file SinglePhaseFluidProperties.C.

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

◆ rho_from_p_T() [1/2]

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

overridevirtual

Density from pressure and temperature.

Parameters

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

Returns: density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 83 of file Water97FluidProperties.C.

Referenced by k_from_p_T(), mu_from_p_T(), rho_e_dpT(), rho_from_p_T(), rho_mu(), and rho_mu_dpT().

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

◆ rho_from_p_T() [2/2]

void Water97FluidProperties::rho_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT
	)		const

overridevirtual

Density and its derivatives from pressure and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 121 of file Water97FluidProperties.C.

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

◆ rho_mu()

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

overridevirtual

Density and viscosity.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 576 of file Water97FluidProperties.C.

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

◆ rho_mu_dpT()

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

overridevirtual

Density and viscosity and their derivatives wrt pressure and temperature.

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 583 of file Water97FluidProperties.C.

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

◆ s()

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

virtualinherited

Definition at line 468 of file SinglePhaseFluidProperties.C.

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

 {
   return s_from_p_T(p, T);
 }

◆ s_from_h_p() [1/2]

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

virtualinherited

Specific entropy from specific enthalpy and pressure.

Parameters

[in]	h	specific enthalpy
[in]	p	pressure

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 92 of file SinglePhaseFluidProperties.C.

Referenced by SinglePhaseFluidProperties::T_from_p_h().

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

◆ s_from_h_p() [2/2]

void SinglePhaseFluidProperties::s_from_h_p	(	Real	h,
		Real	p,
		Real &	s,
		Real &	ds_dh,
		Real &	ds_dp
	)		const

virtualinherited

Specific entropy and its derivatives from specific enthalpy and pressure.

Parameters

[in]	h	specific enthalpy
[in]	p	pressure
[out]	s	specific entropy
[out]	ds_dh	derivative of specific entropy w.r.t. specific enthalpy
[out]	ds_dp	derivative of specific entropy w.r.t. pressure

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 98 of file SinglePhaseFluidProperties.C.

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

◆ s_from_p_T() [1/2]

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

overridevirtual

Specific entropy from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature

Reimplemented from SinglePhaseFluidProperties.

Definition at line 648 of file Water97FluidProperties.C.

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

◆ s_from_p_T() [2/2]

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

overridevirtual

Specific entropy and its derivatives from pressure and temperature.

Parameters

[in]	p	pressure
[in]	T	temperature
[out]	s	specific entropy
[out]	ds_dp	derivative of specific entropy w.r.t. pressure
[out]	ds_dT	derivative of specific entropy w.r.t. temperature

Reimplemented from SinglePhaseFluidProperties.

Definition at line 689 of file Water97FluidProperties.C.

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

◆ s_from_T_v() [1/2]

Real SinglePhaseFluidProperties::s_from_T_v	(	Real	T,
		Real	v
	)		const

virtualinherited

Specific entropy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 557 of file SinglePhaseFluidProperties.C.

Referenced by IdealRealGasMixtureFluidProperties::s_from_T_v().

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

◆ s_from_T_v() [2/2]

void SinglePhaseFluidProperties::s_from_T_v	(	Real	T,
		Real	v,
		Real &	s,
		Real &	ds_dT,
		Real &	ds_dv
	)		const

virtualinherited

Specific entropy and its derivatives from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume
[out]	s	specific entropy (J/kg)
[out]	ds_dT	derivative of specific entropy w.r.t. temperature
[out]	ds_dv	derivative of specific entropy w.r.t. specific volume

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 563 of file SinglePhaseFluidProperties.C.

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

◆ s_from_v_e() [1/2]

Real SinglePhaseFluidProperties::s_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Specific entropy from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 81 of file SinglePhaseFluidProperties.C.

Referenced by StagnationPressureAux::computeValue(), and StagnationTemperatureAux::computeValue().

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

◆ s_from_v_e() [2/2]

void SinglePhaseFluidProperties::s_from_v_e	(	Real	v,
		Real	e,
		Real &	s,
		Real &	ds_dv,
		Real &	ds_de
	)		const

virtualinherited

Specific entropy and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[out]	s	specific entropy
[out]	ds_dv	derivative of specific entropy w.r.t. specific volume
[out]	ds_de	derivative of specific entropy w.r.t. specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 87 of file SinglePhaseFluidProperties.C.

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

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 32 of file FluidProperties.h.

32 {}

◆ subregion2ph()

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

protected

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

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

Parameters

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

Returns: subregion index

Definition at line 1624 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ subregion3()

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

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

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

Parameters

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

Returns: subregion index

Definition at line 1261 of file Water97FluidProperties.C.

Referenced by densityRegion3().

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

◆ subregion3ph()

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

protected

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

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

Parameters

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

Returns: subregion index

Definition at line 1644 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ subregionVolume()

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

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

Parameters

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

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

Definition at line 1511 of file Water97FluidProperties.C.

Referenced by densityRegion3().

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

◆ T_from_p_h()

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

virtualinherited

Temperature from pressure and specific enthalpy.

Parameters

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

Returns: Temperature (K)

Reimplemented in IdealGasFluidProperties.

Definition at line 601 of file SinglePhaseFluidProperties.C.

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

◆ T_from_v_e() [1/2]

Real SinglePhaseFluidProperties::T_from_v_e	(	Real	v,
		Real	e
	)		const

virtualinherited

Temperature from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy

Returns: sound speed

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 39 of file SinglePhaseFluidProperties.C.

Referenced by FluidPropertiesMaterial::computeQpProperties(), StagnationTemperatureAux::computeValue(), TemperatureAux::computeValue(), IdealRealGasMixtureFluidProperties::p_T_from_v_e(), and SinglePhaseFluidProperties::T_from_p_h().

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

◆ T_from_v_e() [2/2]

void SinglePhaseFluidProperties::T_from_v_e	(	Real	v,
		Real	e,
		Real &	T,
		Real &	dT_dv,
		Real &	dT_de
	)		const

virtualinherited

Temperature and its derivatives from specific volume and specific internal energy.

Parameters

[in]	v	specific volume
[in]	e	specific internal energy
[out]	T	temperature
[out]	dT_dv	derivative of temperature w.r.t. specific volume
[out]	dT_de	derivative of temperature w.r.t. specific internal energy

Reimplemented in IdealGasFluidProperties, and StiffenedGasFluidProperties.

Definition at line 45 of file SinglePhaseFluidProperties.C.

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

◆ temperature_from_ph()

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

virtual

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1670 of file Water97FluidProperties.C.

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

◆ temperature_from_ph1()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1721 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ temperature_from_ph2a()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1734 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ temperature_from_ph2b()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1747 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ temperature_from_ph2c()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1760 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ temperature_from_ph3a()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1788 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ temperature_from_ph3b()

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

protected

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

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 1802 of file Water97FluidProperties.C.

Referenced by temperature_from_ph().

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

◆ tempXY()

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

protected

Boundaries between subregions in region 3.

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

Parameters

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

Returns: temperature (K) along the boundary

Definition at line 1450 of file Water97FluidProperties.C.

Referenced by subregion3().

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

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ subregionEnum

Constructor & Destructor Documentation

◆ Water97FluidProperties()

◆ ~Water97FluidProperties()

Member Function Documentation

◆ b23p()

◆ b23T()

◆ b2bc()

◆ b3ab()

◆ beta()

◆ beta_from_p_T()

◆ c()

◆ c_from_p_T()

◆ c_from_v_e() [1/2]

◆ c_from_v_e() [2/2]

◆ cp_from_p_T()

◆ cp_from_v_e()

◆ criticalDensity()

◆ criticalInternalEnergy()

◆ criticalPressure()

◆ criticalTemperature()

◆ cv_from_p_T()

◆ cv_from_T_v()

◆ cv_from_v_e()

◆ d2gamma1_dpi2()

◆ d2gamma1_dpitau()

◆ d2gamma1_dtau2()

◆ d2gamma2_dpi2()

◆ d2gamma2_dpitau()

◆ d2gamma2_dtau2()

◆ d2gamma5_dpi2()

◆ d2gamma5_dpitau()

◆ d2gamma5_dtau2()

◆ d2phi3_ddelta2()

◆ d2phi3_ddeltatau()

◆ d2phi3_dtau2()

◆ densityRegion3()

◆ dgamma1_dpi()

◆ dgamma1_dtau()

◆ dgamma2_dpi()

◆ dgamma2_dtau()

◆ dgamma5_dpi()

◆ dgamma5_dtau()

◆ dphi3_ddelta()

◆ dphi3_dtau()

◆ e()

◆ e_dpT()

◆ e_from_p_rho() [1/2]

◆ e_from_p_rho() [2/2]

◆ e_from_p_T() [1/2]

◆ e_from_p_T() [2/2]

◆ e_from_T_v() [1/2]

◆ e_from_T_v() [2/2]

◆ e_from_v_h() [1/2]

◆ e_from_v_h() [2/2]

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ g_from_v_e()

◆ gamma1()

◆ gamma2()

◆ gamma5()

◆ gamma_from_p_T()

◆ h()

◆ h_dpT()

◆ h_from_p_T() [1/2]

◆ h_from_p_T() [2/2]

◆ h_from_T_v() [1/2]

◆ h_from_T_v() [2/2]

◆ henryConstant()

◆ henryConstant_dT()

◆ henryConstantIAPWS()

◆ henryConstantIAPWS_dT()