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 Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	Water97FluidProperties (const InputParameters &parameters)

virtual	~Water97FluidProperties ()

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

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

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

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

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

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

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

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

virtual ADReal	rho_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

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

virtual ADReal	v_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	v_from_p_T_template (const T &pressure, const T &temperature) const

virtual void	v_from_p_T (Real pressure, Real temperature, Real &v, Real &dv_dp, Real &dv_dT) const override

virtual void	v_from_p_T (const ADReal &pressure, const ADReal &temperature, ADReal &v, ADReal &dv_dp, ADReal &dv_dT) const override

template<typename T >
void	v_from_p_T_template (const T &pressure, const T &temperature, T &v, T &dv_dp, T &dv_dt) const

template<typename T >
T	rho_from_p_T_template (const T &pressure, const T &temperature) const

template<typename T >
void	rho_from_p_T_template (const T &pressure, const T &temperature, T &rho, T &drho_dp, T &drho_dt) const

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

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

Real	p_from_v_e (Real v, Real e) const override

ADReal	p_from_v_e (const ADReal &v, const ADReal &e) const override

template<typename T >
T	p_from_v_e_template (const T &v, const T &e) const

virtual Real	e_from_p_rho (Real p, Real rho) const override

virtual ADReal	e_from_p_rho (const ADReal &p, const ADReal &rho) const override

template<typename T >
T	e_from_p_rho_template (const T &p, const T &rho) const

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

void	e_from_p_rho (const ADReal &p, const ADReal &rho, ADReal &e, ADReal &de_dp, ADReal &de_drho) const override

template<typename T >
void	e_from_p_rho_template (const T &p, const T &rho, T &e, T &de_dp, T &de_drho) const

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

virtual ADReal	e_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	e_from_p_T_template (const T &pressure, const T &temperature) const

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

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

template<typename T >
void	e_from_p_T_template (const T &pressure, const T &temperature, T &e, T &de_dp, T &de_dT) const

ADReal	e_from_v_h (const ADReal &v, const ADReal &h) const override

Real	T_from_v_e (Real v, Real e) const override

ADReal	T_from_v_e (const ADReal &v, const ADReal &e) const override

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

virtual ADReal	c_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	c_from_p_T_template (const T &pressure, const T &temperature) const

virtual ADReal	c_from_v_e (const ADReal &v, const ADReal &e) const override

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

virtual ADReal	cp_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	cp_from_p_T_template (const T &pressure, const T &temperature) const

virtual ADReal	cp_from_v_e (const ADReal &v, const ADReal &e) const override

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

virtual ADReal	cv_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	cv_from_p_T_template (const T &pressure, const T &temperature) const

virtual ADReal	cv_from_v_e (const ADReal &v, const ADReal &e) const override

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

virtual ADReal	mu_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	mu_from_p_T_template (const T &pressure, const T &temperature) const

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

template<typename T >
T	mu_from_rho_T_template (const T &density, const T &temperature) const

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

ADReal	mu_from_v_e (const ADReal &v, const ADReal &e) const override

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

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

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

virtual ADReal	k_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	k_from_p_T_template (const T &pressure, const T &temperature) const

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

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

template<typename T >
T	k_from_rho_T_template (const T &density, const T &temperature) const

Real	k_from_v_e (Real v, Real e) const override

ADReal	k_from_v_e (const ADReal &v, const ADReal &e) const override

template<typename T >
T	k_from_v_e_template (const T &v, const T &e) const

	propfuncWithDefinitionOverride (s, p, T)

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

virtual ADReal	h_from_p_T (const ADReal &pressure, const ADReal &temperature) const override

template<typename T >
T	h_from_p_T_template (const T &pressure, const T &temperature) const

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

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

template<typename T >
void	h_from_p_T_template (const T &pressure, const T &temperature, T &h, T &dh_dp, T &dh_dT) const

virtual Real	s_from_h_p (Real enthalpy, Real pressure) const override

virtual ADReal	s_from_h_p (const ADReal &enthalpy, const ADReal &pressure) const override

virtual void	s_from_h_p (Real enthalpy, Real pressure, Real &s, Real &ds_dh, Real &ds_dp) const override

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

virtual void	vaporPressure (Real temperature, Real &psat, Real &dpsat_dT) const override

template<typename T >
void	vaporPressureTemplate (const T &temperature, T &psat, T &dpsat_dT) const

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

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

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

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

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

unsigned int	inRegionPH (Real pressure, Real enthalpy) const
	Determines the phase region that the given pressure and enthaply 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...

template<typename T >
T	subregionVolume (const T &pi, const T &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...

template<typename T >
T	densityRegion3 (const T &pressure, const T &temperature) const
	Density function for Region 3 - supercritical water and steam. More...

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

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

virtual ADReal	T_from_p_h (const ADReal &pressure, const ADReal &enthalpy) const override

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

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

Real	henryConstant (Real temperature, const std::vector< Real > &coeffs) const
	IAPWS formulation of Henry's law constant for dissolution in water From Guidelines on the Henry's constant and vapour liquid distribution constant for gases in H20 and D20 at high temperatures, IAPWS (2004) More...

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

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

template<typename T >
std::pair< T, T >	p_T_from_v_e (const T &v, const T &e) const
	Computes the pressure (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific internal energy. More...

template<typename T >
std::pair< T, T >	rho_T_from_v_e (const T &v, const T &e) const
	Computes the density (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific internal energy. More...

template<typename T >
std::pair< T, T >	p_T_from_v_h (const T &v, const T &h) const
	Computes the pressure (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific enthalpy. More...

template<typename T >
void	p_T_from_v_h (const T &v, const T &h, Real p0, Real T0, T &pressure, T &temperature, bool &conversion_succeeded) const
	Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

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

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

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

virtual ADReal	vaporPressure (const ADReal &T) const

virtual ADReal	vaporTemperature (const ADReal &p) const

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

template<typename CppType >
void	v_e_from_p_T (const CppType &p, const CppType &T, CppType &v, CppType &e) const

template<typename CppType >
void	v_e_from_p_T (const CppType &p, const CppType &T, CppType &v, CppType &dv_dp, CppType &dv_dT, CppType &e, CppType &de_dp, CppType &de_dT) const

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

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

template<typename CppType >
void	p_T_from_v_e (const CppType &v, const CppType &e, Real p0, Real T0, CppType &p, CppType &T, bool &conversion_succeeded) const
	Determines (p,T) from (v,e) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<typename T >
void	p_T_from_v_h (const T &v, const T &h, Real p0, Real T0, T &pressure, T &temperature, bool &conversion_succeeded) const
	Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<typename T >
void	p_T_from_h_s (const T &h, const T &s, Real p0, Real T0, T &pressure, T &temperature, bool &conversion_succeeded) const
	Determines (p,T) from (h,s) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<>
std::pair< Real, Real >	makeZeroAndOne (const Real &)

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

bool	needThreadedCopy () const override final

virtual Real	spatialValue (const Point &) const

virtual const std::vector< Point >	spatialPoints () const

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

std::set< UserObjectName >	getDependObjects () const

const std::set< std::string > &	getRequestedItems () override

const std::set< std::string > &	getSuppliedItems () override

unsigned int	systemNumber () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	haveParameter (const std::string &name) const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &name) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

std::string	messagePrefix (const bool hit_prefix=true) const

std::string	errorPrefix (const std::string &) const

void	mooseError (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecatedNoTrace (Args &&... args) const

void	mooseInfo (Args &&... args) const

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

UserObjectName	getUserObjectName (const std::string &param_name) const

const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const

const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const

const UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObjectBase &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const

const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const

const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name_in)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > >	getKokkosBlockMaterialProperty (const MaterialPropertyName &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

virtual bool	getMaterialPropertyCalled () const

virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const

virtual void	resolveOptionalProperties ()

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

bool	isImplicit ()

Moose::StateArg	determineState () const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (const FunctionName &name) const

bool	hasFunction (const std::string &param_name) const

bool	hasFunctionByName (const FunctionName &name) const

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

const T &	getKokkosFunction (const std::string &name) const

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

bool	hasKokkosFunction (const std::string &param_name) const

bool	hasKokkosFunctionByName (const FunctionName &name) const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const

T &	getSampler (const std::string &name)

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshChanged ()

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header, NameFunc &&name_func)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

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

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

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
template<typename T >
T	gamma1 (const T &pi, const T &tau) const
	Gibbs free energy in Region 1 - single phase liquid region. More...

template<typename T >
T	dgamma1_dpi (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 1 wrt pi. More...

template<typename T >
T	d2gamma1_dpi2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 1 wrt pi. More...

template<typename T >
T	dgamma1_dtau (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 1 wrt tau. More...

template<typename T >
T	d2gamma1_dtau2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 1 wrt tau. More...

template<typename T >
T	d2gamma1_dpitau (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 1 wrt pi and tau. More...

template<typename T >
T	gamma2 (const T &pi, const T &tau) const
	Gibbs free energy in Region 2 - superheated steam. More...

template<typename T >
T	dgamma2_dpi (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 2 wrt pi. More...

template<typename T >
T	d2gamma2_dpi2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 2 wrt pi. More...

template<typename T >
T	dgamma2_dtau (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 2 wrt tau. More...

template<typename T >
T	d2gamma2_dtau2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 2 wrt tau. More...

template<typename T >
T	d2gamma2_dpitau (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 2 wrt pi and tau. More...

template<typename T >
T	phi3 (const T &delta, const T &tau) const
	Helmholtz free energy in Region 3. More...

template<typename T >
T	dphi3_ddelta (const T &delta, const T &tau) const
	Derivative of Helmholtz free energy in Region 3 wrt delta. More...

template<typename T >
T	d2phi3_ddelta2 (const T &delta, const T &tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt delta. More...

template<typename T >
T	dphi3_dtau (const T &delta, const T &tau) const
	Derivative of Helmholtz free energy in Region 3 wrt tau. More...

template<typename T >
T	d2phi3_dtau2 (const T &delta, const T &tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt tau. More...

template<typename T >
T	d2phi3_ddeltatau (const T &delta, const T &tau) const
	Second derivative of Helmholtz free energy in Region 3 wrt delta and tau. More...

template<typename T >
T	gamma5 (const T &pi, const T &tau) const
	Gibbs free energy in Region 5. More...

template<typename T >
T	dgamma5_dpi (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 5 wrt pi. More...

template<typename T >
T	d2gamma5_dpi2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 5 wrt pi. More...

template<typename T >
T	dgamma5_dtau (const T &pi, const T &tau) const
	Derivative of Gibbs free energy in Region 5 wrt tau. More...

template<typename T >
T	d2gamma5_dtau2 (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 5 wrt tau. More...

template<typename T >
T	d2gamma5_dpitau (const T &pi, const T &tau) const
	Second derivative of Gibbs free energy in Region 5 wrt pi and tau. More...

template<typename T >
T	tempXY (const T &pressure, subregionEnum xy) const
	Boundaries between subregions in region 3. 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...

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

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

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

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

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

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

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

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

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObjectBase &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

const ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const

const ReporterName &	getReporterName (const std::string &param_name) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const

unsigned int	coupledScalarComponents (const std::string &var_name) const

unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const

libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const

const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const

virtual void	checkMaterialProperty (const std::string &name, const unsigned int state)

virtual void	getKokkosMaterialPropertyHook (const std::string &, const unsigned int)

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) const

void	checkExecutionStage ()

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

Static Protected Member Functions
template<typename T , typename Functor >
static void	xyDerivatives (const T x, const T &y, T &z, T &dz_dx, T &dz_dy, const Functor &z_from_x_y)
	Computes the dependent variable z and its derivatives with respect to the independent variables x and y using the simple two parameter `z_from_x_y` functor. More...

template<typename T >
static std::pair< T, T >	makeZeroAndOne (const T &)
	Given a type example, this method returns zero and unity representations of that type (first and second members of returned pair respectively) More...

static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

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< Real, 4 >	_mu_H0 {{1.67752, 2.20462, 0.6366564, -0.241605}}
	Constants from Release on the IAPWS Formulation 2008 for the Viscosity of Ordinary Water Substance. More...

const std::array< std::array< Real, 7 >, 6 >	_mu_Hij

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	_tolerance
	Newton's method may be used to convert between variable sets. More...

const Real	_T_initial_guess
	Initial guess for temperature (or temperature used to compute the initial guess) More...

const Real	_p_initial_guess
	Initial guess for pressure (or pressure used to compute the initial guess) More...

const unsigned int	_max_newton_its
	Maximum number of iterations for the variable conversion newton solves. More...

const bool	_verbose_newton
	Whether to output information about newton solves to console. More...

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

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

const Moose::CoordinateSystemType &	_coord_sys

const THREAD_ID	_tid

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

Assembly &	_assembly

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_mi_params

const std::string	_mi_name

const MooseObjectName	_mi_moose_object_name

FEProblemBase &	_mi_feproblem

SubProblem &	_mi_subproblem

const THREAD_ID	_mi_tid

const bool	_is_kokkos_object

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

const Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Member Functions
template<typename T >
std::pair< T, T >	T_drhodT_from_p_rho (const T &p, const T &rho) const
	Computes the temperature (first member of the pair) and the derivative of density (second member of the pair) with respect to temperature as a function of pressure and density. More...

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

	v

e	v

e e	v

e e e	v

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

e e e e	p

e e e e s	p

e e e e s T T	p

e e e e s T T T	p

e e e e s T T T T	p

e e e e s T T T T T	p

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

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

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

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

e e e e s T	rho

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

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

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

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 38 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 667 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 24 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 36 of file Water97FluidProperties.C.

36 {}

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

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

◆ 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 705 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)
     mooseException(
         name(), ": b2bc(): Pressure ", pressure, " is outside range of 6.5467 MPa <= p <= 100 MPa");
 
   Real pi = pressure / 1.0e6;
 
   return (0.26526571908428e4 + std::sqrt((pi - 0.45257578905948e1) / 0.12809002730136e-3)) * 1.0e3;
 }

◆ 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 903 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)
     mooseException(
         name(), ": b3ab(): Pressure ", pressure, "is outside range of 16.529 MPa <= p <= 100 MPa");
 
   Real pi = pressure / 1.0e6;
   Real eta = 0.201464004206875e4 + 0.374696550136983e1 * pi - 0.219921901054187e-1 * pi * pi +
              0.87513168600995e-4 * pi * pi * pi;
 
   return eta * 1.0e3;
 }

◆ c_from_p_T() [1/2]

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

overridevirtual

Definition at line 232 of file Water97FluidProperties.C.

Referenced by c_from_v_e().

 {
   return c_from_p_T_template(p, T);
 }

◆ c_from_p_T() [2/2]

ADReal Water97FluidProperties::c_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 238 of file Water97FluidProperties.C.

 {
   return c_from_p_T_template(p, T);
 }

◆ c_from_p_T_template()

template<typename T >

T Water97FluidProperties::c_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2222 of file Water97FluidProperties.h.

Referenced by c_from_p_T().

 {
   T speed2, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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
       T 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("inRegion() has given an incorrect region");
   }
 
   using std::sqrt;
   return sqrt(speed2);
 }

◆ c_from_v_e()

ADReal Water97FluidProperties::c_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 225 of file Water97FluidProperties.C.

 {
   const auto [p, T] = p_T_from_v_e(v, e);
   return c_from_p_T(p, T);
 }

◆ cp_from_p_T() [1/2]

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

overridevirtual

Definition at line 244 of file Water97FluidProperties.C.

Referenced by cp_from_v_e().

 {
   return cp_from_p_T_template(pressure, temperature);
 }

◆ cp_from_p_T() [2/2]

ADReal Water97FluidProperties::cp_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 250 of file Water97FluidProperties.C.

 {
   return cp_from_p_T_template(pressure, temperature);
 }

◆ cp_from_p_T_template()

template<typename T >

T Water97FluidProperties::cp_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2283 of file Water97FluidProperties.h.

Referenced by cp_from_p_T().

 {
   T specific_heat, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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
       T 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("inRegion() has given an incorrect region");
   }
   return specific_heat;
 }

◆ cp_from_v_e()

ADReal Water97FluidProperties::cp_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 256 of file Water97FluidProperties.C.

 {
   const auto [p, T] = p_T_from_v_e(v, e);
   return cp_from_p_T(p, T);
 }

◆ criticalDensity()

Real Water97FluidProperties::criticalDensity ( ) const

overridevirtual

Critical density.

Returns: critical density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 63 of file Water97FluidProperties.C.

 {
   return _rho_critical;
 }

◆ criticalInternalEnergy()

Real SinglePhaseFluidProperties::criticalInternalEnergy ( ) const

virtualinherited

Critical specific internal energy.

Returns: specific internal energy (J/kg)

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 320 of file SinglePhaseFluidProperties.C.

 {
   return e_from_p_rho(criticalPressure(), criticalDensity());
 }

◆ criticalPressure()

Real Water97FluidProperties::criticalPressure ( ) const

overridevirtual

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

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

 {
   return _T_critical;
 }

◆ cv_from_p_T() [1/2]

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

overridevirtual

Definition at line 263 of file Water97FluidProperties.C.

Referenced by cv_from_v_e().

 {
   return cv_from_p_T_template(pressure, temperature);
 }

◆ cv_from_p_T() [2/2]

ADReal Water97FluidProperties::cv_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 269 of file Water97FluidProperties.C.

 {
   return cv_from_p_T_template(pressure, temperature);
 }

◆ cv_from_p_T_template()

template<typename T >

T Water97FluidProperties::cv_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2334 of file Water97FluidProperties.h.

Referenced by cv_from_p_T().

 {
   T specific_heat, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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
       T 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("inRegion() has given an incorrect region");
   }
   return specific_heat;
 }

◆ cv_from_v_e()

ADReal Water97FluidProperties::cv_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 275 of file Water97FluidProperties.C.

 {
   const auto [p, T] = p_T_from_v_e(v, e);
   return cv_from_p_T(p, T);
 }

◆ d2gamma1_dpi2()

template<typename T >

T Water97FluidProperties::d2gamma1_dpi2	(	const T &	pi,
		const T &	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 1599 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::d2gamma1_dpitau	(	const T &	pi,
		const T &	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 1635 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::d2gamma1_dtau2	(	const T &	pi,
		const T &	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 1623 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::d2gamma2_dpi2	(	const T &	pi,
		const T &	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 1668 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = -1.0 / pi / pi;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::d2gamma2_dpitau	(	const T &	pi,
		const T &	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 1719 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = 0.0;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::d2gamma2_dtau2	(	const T &	pi,
		const T &	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 1701 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T 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
   T 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()

template<typename T >

T Water97FluidProperties::d2gamma5_dpi2	(	const T &	pi,
		const T &	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 1842 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = -1.0 / pi / pi;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::d2gamma5_dpitau	(	const T &	pi,
		const T &	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 1893 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = 0.0;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::d2gamma5_dtau2	(	const T &	pi,
		const T &	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 1875 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T 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
   T 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()

template<typename T >

T Water97FluidProperties::d2phi3_ddelta2	(	const T &	delta,
		const T &	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 1759 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::d2phi3_ddeltatau	(	const T &	delta,
		const T &	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 1794 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::d2phi3_dtau2	(	const T &	delta,
		const T &	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 1782 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), and e_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::densityRegion3	(	const T &	pressure,
		const T &	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 1543 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), and rho_from_p_T_template().

 {
   using std::exp;
 
   // 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(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(temperature));
 
   T vstar, pi, theta;
   Real 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];
 
   T sum = 0.0;
   T 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 * 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()

template<typename T >

T Water97FluidProperties::dgamma1_dpi	(	const T &	pi,
		const T &	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 1516 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), and rho_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::dgamma1_dtau	(	const T &	pi,
		const T &	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 1611 of file Water97FluidProperties.h.

Referenced by e_from_p_T_template(), h_from_p_T_template(), and inRegionPH().

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

template<typename T >

T Water97FluidProperties::dgamma2_dpi	(	const T &	pi,
		const T &	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 1528 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), and rho_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = 1.0 / pi;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::dgamma2_dtau	(	const T &	pi,
		const T &	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 1684 of file Water97FluidProperties.h.

Referenced by e_from_p_T_template(), h_from_p_T_template(), and inRegionPH().

 {
   // Ideal gas part of the Gibbs free energy
   T 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
   T 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()

template<typename T >

T Water97FluidProperties::dgamma5_dpi	(	const T &	pi,
		const T &	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 1827 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), and rho_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T dg0 = 1.0 / pi;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::dgamma5_dtau	(	const T &	pi,
		const T &	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 1858 of file Water97FluidProperties.h.

Referenced by e_from_p_T_template(), and h_from_p_T_template().

 {
   // Ideal gas part of the Gibbs free energy
   T 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
   T 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()

template<typename T >

T Water97FluidProperties::dphi3_ddelta	(	const T &	delta,
		const T &	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 1748 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), e_from_p_T_template(), and h_from_p_T_template().

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

template<typename T >

T Water97FluidProperties::dphi3_dtau	(	const T &	delta,
		const T &	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 1771 of file Water97FluidProperties.h.

Referenced by e_from_p_T_template(), and h_from_p_T_template().

 {
   T 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_from_p_rho() [1/4]

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

overridevirtual

Definition at line 151 of file Water97FluidProperties.C.

Referenced by p_from_v_e_template().

 {
   return e_from_p_rho_template(p, rho);
 }

◆ e_from_p_rho() [2/4]

ADReal Water97FluidProperties::e_from_p_rho	(	const ADReal &	p,
		const ADReal &	rho
	)		const

overridevirtual

Definition at line 157 of file Water97FluidProperties.C.

 {
   return e_from_p_rho_template(p, rho);
 }

◆ e_from_p_rho() [3/4]

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

override

Definition at line 163 of file Water97FluidProperties.C.

 {
   e_from_p_rho_template(p, rho, e, de_dp, de_drho);
 }

◆ e_from_p_rho() [4/4]

void Water97FluidProperties::e_from_p_rho	(	const ADReal &	p,
		const ADReal &	rho,
		ADReal &	e,
		ADReal &	de_dp,
		ADReal &	de_drho
	)		const

override

Definition at line 170 of file Water97FluidProperties.C.

 {
   e_from_p_rho_template(p, rho, e, de_dp, de_drho);
 }

◆ e_from_p_rho_template() [1/2]

template<typename T >

T Water97FluidProperties::e_from_p_rho_template	(	const T &	p,
		const T &	rho
	)		const

Definition at line 2456 of file Water97FluidProperties.h.

Referenced by e_from_p_rho(), and e_from_p_rho_template().

 {
   const auto temperature = T_drhodT_from_p_rho(p, rho).first;
   return e_from_p_T_template(p, temperature);
 }

◆ e_from_p_rho_template() [2/2]

template<typename T >

void Water97FluidProperties::e_from_p_rho_template	(	const T &	p,
		const T &	rho,
		T &	e,
		T &	de_dp,
		T &	de_drho
	)		const

Definition at line 2464 of file Water97FluidProperties.h.

 {
   auto functor = [this](const auto & p, const auto & rho)
   { return this->e_from_p_rho_template(p, rho); };
 
   xyDerivatives(p, rho, e, de_dp, de_drho, functor);
 }

◆ e_from_p_T() [1/4]

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

overridevirtual

Definition at line 177 of file Water97FluidProperties.C.

Referenced by e_from_p_T_template(), and e_from_v_h().

 {
   return e_from_p_T_template(pressure, temperature);
 }

◆ e_from_p_T() [2/4]

ADReal Water97FluidProperties::e_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 183 of file Water97FluidProperties.C.

 {
   return e_from_p_T_template(pressure, temperature);
 }

◆ e_from_p_T() [3/4]

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

overridevirtual

Definition at line 189 of file Water97FluidProperties.C.

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

◆ e_from_p_T() [4/4]

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

overridevirtual

Definition at line 196 of file Water97FluidProperties.C.

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

◆ e_from_p_T_template() [1/2]

template<typename T >

T Water97FluidProperties::e_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2046 of file Water97FluidProperties.h.

Referenced by e_from_p_rho_template(), and e_from_p_T().

 {
   T internal_energy, pi, tau;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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
       T density3 = densityRegion3(pressure, temperature);
       T 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("inRegion() has given an incorrect region");
   }
   // Output in J/kg
   return internal_energy;
 }

◆ e_from_p_T_template() [2/2]

template<typename T >

void Water97FluidProperties::e_from_p_T_template	(	const T &	pressure,
		const T &	temperature,
		T &	e,
		T &	de_dp,
		T &	de_dT
	)		const

Definition at line 2095 of file Water97FluidProperties.h.

 {
   T pi, tau, dinternal_energy_dp, dinternal_energy_dT;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(temperature));
   switch (region)
   {
     case 1:
     {
       pi = pressure / _p_star[0];
       tau = _T_star[0] / temperature;
       T dgdp = dgamma1_dpi(pi, tau);
       T 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;
       T dgdp = dgamma2_dpi(pi, tau);
       T 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
       T density3 = densityRegion3(pressure, temperature);
       T delta = density3 / _rho_critical;
       tau = _T_star[2] / temperature;
       T dpdd = dphi3_ddelta(delta, tau);
       T d2pddt = d2phi3_ddeltatau(delta, tau);
       T 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;
       T dgdp = dgamma5_dpi(pi, tau);
       T 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("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_v_h()

ADReal Water97FluidProperties::e_from_v_h	(	const ADReal &	v,
		const ADReal &	h
	)		const

override

Definition at line 206 of file Water97FluidProperties.C.

 {
   const auto [p, T] = p_T_from_v_h(v, h);
   return e_from_p_T(p, T);
 }

◆ e_spndl_from_v()

Real SinglePhaseFluidProperties::e_spndl_from_v ( Real v ) const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 483 of file SinglePhaseFluidProperties.C.

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

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 33 of file FluidProperties.h.

33 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 35 of file FluidProperties.h.

35 {}

◆ fluidName()

std::string Water97FluidProperties::fluidName ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name

Reimplemented from SinglePhaseFluidProperties.

Definition at line 39 of file Water97FluidProperties.C.

 {
   return "water";
 }

◆ gamma1()

template<typename T >

T Water97FluidProperties::gamma1	(	const T &	pi,
		const T &	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 1588 of file Water97FluidProperties.h.

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

template<typename T >

T Water97FluidProperties::gamma2	(	const T &	pi,
		const T &	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 1647 of file Water97FluidProperties.h.

 {
   using std::log;
 
   // Ideal gas part of the Gibbs free energy
   T sum0 = 0.0;
   for (std::size_t i = 0; i < _n02.size(); ++i)
     sum0 += _n02[i] * MathUtils::pow(tau, _J02[i]);
 
   T g0 = log(pi) + sum0;
 
   // Residual part of the Gibbs free energy
   T 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()

template<typename T >

T Water97FluidProperties::gamma5	(	const T &	pi,
		const T &	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 1806 of file Water97FluidProperties.h.

 {
   using std::log;
 
   // Ideal gas part of the Gibbs free energy
   T sum0 = 0.0;
   for (std::size_t i = 0; i < _n05.size(); ++i)
     sum0 += _n05[i] * MathUtils::pow(tau, _J05[i]);
 
   T g0 = log(pi) + sum0;
 
   // Residual part of the Gibbs free energy
   T 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;
 }

◆ h_from_p_T() [1/4]

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

overridevirtual

Definition at line 492 of file Water97FluidProperties.C.

Referenced by inRegionPH(), and T_from_p_h_ad().

 {
   return h_from_p_T_template(pressure, temperature);
 }

◆ h_from_p_T() [2/4]

ADReal Water97FluidProperties::h_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 498 of file Water97FluidProperties.C.

 {
   return h_from_p_T_template(pressure, temperature);
 }

◆ h_from_p_T() [3/4]

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

overridevirtual

Definition at line 504 of file Water97FluidProperties.C.

 {
   h_from_p_T_template(pressure, temperature, h, dh_dp, dh_dT);
 }

◆ h_from_p_T() [4/4]

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

overridevirtual

Definition at line 511 of file Water97FluidProperties.C.

 {
   h_from_p_T_template(pressure, temperature, h, dh_dp, dh_dT);
 }

◆ h_from_p_T_template() [1/2]

template<typename T >

T Water97FluidProperties::h_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2518 of file Water97FluidProperties.h.

Referenced by h_from_p_T(), and h_from_p_T_template().

 {
   T enthalpy, pi, tau, delta;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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
       T 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_template() [2/2]

template<typename T >

void Water97FluidProperties::h_from_p_T_template	(	const T &	pressure,
		const T &	temperature,
		T &	h,
		T &	dh_dp,
		T &	dh_dT
	)		const

Definition at line 2564 of file Water97FluidProperties.h.

 {
   auto functor = [this](const auto & pressure, const auto & temperature)
   { return this->h_from_p_T_template(pressure, temperature); };
 
   xyDerivatives(pressure, temperature, h, dh_dp, dh_dT, functor);
 }

◆ henryCoefficients()

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

virtualinherited

Henry's law coefficients for dissolution in water.

Returns: Henry's constant coefficients

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

Definition at line 378 of file SinglePhaseFluidProperties.C.

Referenced by TabulatedFluidProperties::henryCoefficients().

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

◆ henryConstant() [1/3]

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

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

Parameters

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

Definition at line 950 of file Water97FluidProperties.C.

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

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

◆ henryConstant() [2/3]

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

Definition at line 975 of file Water97FluidProperties.C.

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

◆ henryConstant() [3/3]

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

Definition at line 1019 of file Water97FluidProperties.C.

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

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 34 of file FluidProperties.h.

34 {}

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

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), and rho_from_p_T_template().

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

◆ inRegionPH()

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

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

Referenced by T_from_p_h_ad().

 {
   unsigned int region;
 
   // Need to calculate enthalpies at the boundaries to delineate regions
   Real p273 = vaporPressure(273.15);
   Real p623 = vaporPressure(623.15);
 
   if (enthalpy < h_from_p_T(pressure, 273.15))
     mooseException("Enthalpy ", enthalpy, " is out of range in ", name(), ": inRegionPH()");
 
   if (pressure >= p273 && pressure <= p623)
   {
     // Use region 1 definition of h_from_p_T to get the lower bound
     if (enthalpy <=
         _Rw * _T_star[0] *
             dgamma1_dtau(pressure / _p_star[0], _T_star[0] / vaporTemperature(pressure)))
       region = 1;
     // Use region 2 definition of h_from_p_T to get the upper bound
     else if (enthalpy <=
              _Rw * _T_star[1] *
                  dgamma2_dtau(pressure / _p_star[1], _T_star[1] / vaporTemperature(pressure)))
       region = 4;
     else if (enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else if (enthalpy <= h_from_p_T(pressure, 2273.15))
       region = 5;
     else
       mooseException("Enthalpy ", enthalpy, " is out of range in ", name(), ": inRegionPH()");
   }
   else if (pressure > p623 && pressure <= 50.0e6)
   {
     if (enthalpy <= h_from_p_T(pressure, 623.15))
       region = 1;
     else if (enthalpy <= h_from_p_T(pressure, b23T(pressure)))
       region = 3;
     else if (enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else if (enthalpy <= h_from_p_T(pressure, 2273.15))
       region = 5;
     else
       mooseException("Enthalpy ", enthalpy, " is out of range in ", name(), ": inRegionPH()");
   }
   else if (pressure > 50.0e6 && pressure <= 100.0e6)
   {
     if (enthalpy <= h_from_p_T(pressure, 623.15))
       region = 1;
     else if (enthalpy <= h_from_p_T(pressure, b23T(pressure)))
       region = 3;
     else if (enthalpy <= h_from_p_T(pressure, 1073.15))
       region = 2;
     else
       mooseException("Enthalpy ", enthalpy, " is out of range in ", name(), ": inRegionPH()");
   }
   else
     mooseException("Pressure ", pressure, " is out of range in ", name(), ": inRegionPH()");
 
   return region;
 }

◆ k_from_p_T() [1/3]

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

overridevirtual

Definition at line 395 of file Water97FluidProperties.C.

 {
   return k_from_p_T_template(pressure, temperature);
 }

◆ k_from_p_T() [2/3]

ADReal Water97FluidProperties::k_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 401 of file Water97FluidProperties.C.

 {
   return k_from_p_T_template(pressure, temperature);
 }

◆ k_from_p_T() [3/3]

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

overridevirtual

Definition at line 407 of file Water97FluidProperties.C.

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

◆ k_from_p_T_template()

template<typename T >

T Water97FluidProperties::k_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2424 of file Water97FluidProperties.h.

Referenced by k_from_p_T(), and k_from_v_e_template().

 {
   T rho = this->rho_from_p_T(pressure, temperature);
   return this->k_from_rho_T_template(rho, temperature);
 }

◆ k_from_rho_T()

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

overridevirtual

Definition at line 413 of file Water97FluidProperties.C.

 {
   return k_from_rho_T_template(density, temperature);
 }

◆ k_from_rho_T_template()

template<typename T >

T Water97FluidProperties::k_from_rho_T_template	(	const T &	density,
		const T &	temperature
	)		const

Definition at line 2388 of file Water97FluidProperties.h.

Referenced by k_from_p_T_template(), and k_from_rho_T().

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

◆ k_from_v_e() [1/2]

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

override

Definition at line 419 of file Water97FluidProperties.C.

 {
   return k_from_v_e_template(v, e);
 }

◆ k_from_v_e() [2/2]

ADReal Water97FluidProperties::k_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

override

Definition at line 425 of file Water97FluidProperties.C.

 {
   return k_from_v_e_template(v, e);
 }

◆ k_from_v_e_template()

template<typename T >

T Water97FluidProperties::k_from_v_e_template	(	const T &	v,
		const T &	e
	)		const

Definition at line 2432 of file Water97FluidProperties.h.

Referenced by k_from_v_e().

 {
   const auto [p, temperature] = p_T_from_v_e(v, e);
   return k_from_p_T_template(p, temperature);
 }

◆ makeZeroAndOne() [1/2]

template<typename T >

std::pair< T, T > SinglePhaseFluidProperties::makeZeroAndOne ( const T & )

staticprotectedinherited

Given a type example, this method returns zero and unity representations of that type (first and second members of returned pair respectively)

Definition at line 478 of file SinglePhaseFluidProperties.h.

Referenced by SinglePhaseFluidProperties::xyDerivatives().

 {
   return {T{0, 0}, T{1, 0}};
 }

◆ makeZeroAndOne() [2/2]

template<>

std::pair<Real, Real> SinglePhaseFluidProperties::makeZeroAndOne ( const Real & )

inlineinherited

Definition at line 485 of file SinglePhaseFluidProperties.h.

 {
   return {Real{0}, Real{1}};
 }

◆ molarMass()

Real Water97FluidProperties::molarMass ( ) const

overridevirtual

Molar mass [kg/mol].

Returns: molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 45 of file Water97FluidProperties.C.

 {
   return _Mh2o;
 }

◆ mu_from_p_T() [1/3]

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

overridevirtual

Definition at line 282 of file Water97FluidProperties.C.

 {
   return mu_from_p_T_template(pressure, temperature);
 }

◆ mu_from_p_T() [2/3]

ADReal Water97FluidProperties::mu_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 288 of file Water97FluidProperties.C.

 {
   return mu_from_p_T_template(pressure, temperature);
 }

◆ mu_from_p_T() [3/3]

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

overridevirtual

Definition at line 294 of file Water97FluidProperties.C.

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

◆ mu_from_p_T_template()

template<typename T >

T Water97FluidProperties::mu_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2615 of file Water97FluidProperties.h.

Referenced by mu_from_p_T().

 {
   T rho = this->rho_from_p_T_template(pressure, temperature);
   return this->mu_from_rho_T_template(rho, temperature);
 }

◆ mu_from_rho_T() [1/2]

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

overridevirtual

Definition at line 305 of file Water97FluidProperties.C.

Referenced by mu_from_p_T(), and rho_mu_from_p_T().

 {
   return mu_from_rho_T_template(density, temperature);
 }

◆ mu_from_rho_T() [2/2]

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

Definition at line 311 of file Water97FluidProperties.C.

 {
   const Real mu_star = 1.0e-6;
   const Real rhobar = density / _rho_critical;
   const Real Tbar = temperature / _T_critical;
   const Real drhobar_drho = 1.0 / _rho_critical;
   const 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);
   }
 
   const Real mu0 = 100.0 * std::sqrt(Tbar) / sum0;
   const 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_drhobar = 0.0, dsum1_dTbar = 0.0;
   for (unsigned int i = 0; i < 6; ++i)
   {
     const Real fact = MathUtils::pow(1.0 / Tbar - 1.0, i);
     const Real dfact_dTbar = i * MathUtils::pow(1.0 / Tbar - 1.0, i - 1) / Tbar / Tbar;
 
     for (unsigned int j = 0; j < 7; ++j)
     {
       sum1 += fact * _mu_Hij[i][j] * MathUtils::pow(rhobar - 1.0, j);
       dsum1_dTbar -= dfact_dTbar * _mu_Hij[i][j] * MathUtils::pow(rhobar - 1.0, j);
       dsum1_drhobar += j * fact * _mu_Hij[i][j] * MathUtils::pow(rhobar - 1.0, j - 1);
     }
   }
 
   const Real mu1 = std::exp(rhobar * sum1);
   const Real dmu1_drhobar = (sum1 + rhobar * dsum1_drhobar) * mu1;
   const Real dmu1_dTbar = (rhobar * dsum1_dTbar) * mu1;
 
   // Viscosity and its derivatives are then
   mu = mu_star * mu0 * mu1;
   dmu_drho = mu_star * mu0 * dmu1_drhobar * drhobar_drho;
   dmu_dT = mu_star * (dmu0_dTbar * mu1 + mu0 * dmu1_dTbar) * dTbar_dT + dmu_drho * ddensity_dT;
 }

◆ mu_from_rho_T_template()

template<typename T >

T Water97FluidProperties::mu_from_rho_T_template	(	const T &	density,
		const T &	temperature
	)		const

Definition at line 2475 of file Water97FluidProperties.h.

Referenced by mu_from_p_T_template(), mu_from_rho_T(), and mu_from_v_e().

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

◆ mu_from_v_e()

ADReal Water97FluidProperties::mu_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

override

Definition at line 362 of file Water97FluidProperties.C.

 {
   const auto [rho, T] = rho_T_from_v_e(v, e);
   return mu_from_rho_T_template(rho, T);
 }

◆ p_from_v_e() [1/2]

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

override

Definition at line 139 of file Water97FluidProperties.C.

Referenced by p_T_from_v_e(), and rho_T_from_v_e().

 {
   return p_from_v_e_template(v, e);
 }

◆ p_from_v_e() [2/2]

ADReal Water97FluidProperties::p_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

override

Definition at line 145 of file Water97FluidProperties.C.

 {
   return p_from_v_e_template(v, e);
 }

◆ p_from_v_e_template()

template<typename T >

T Water97FluidProperties::p_from_v_e_template	(	const T &	v,
		const T &	e
	)		const

Definition at line 2440 of file Water97FluidProperties.h.

Referenced by p_from_v_e().

 {
   const auto rho = 1 / v;
   // For NewtonSolve
   // y = e
   // x = rho
   // z = p
   auto lambda = [this](const T & rho, const T & p, T & e, T & de_drho, T & de_dp)
   { e_from_p_rho(p, rho, e, de_dp, de_drho); };
   return FluidPropertiesUtils::NewtonSolve(
              rho, e, _p_initial_guess, _tolerance, lambda, name() + "::p_from_v_e")
       .first;
 }

◆ p_T_from_h_s()

template<typename T >

void SinglePhaseFluidProperties::p_T_from_h_s	(	const T &	h,
		const T &	s,
		Real	p0,
		Real	T0,
		T &	pressure,
		T &	temperature,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (h,s) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	h	specific enthalpy (J / kg)
[in]	s	specific entropy (J/K*kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 599 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::T_from_h_s().

 {
   auto h_lambda = [&](const T & pressure, const T & temperature, T & new_h, T & dh_dp, T & dh_dT)
   { h_from_p_T(pressure, temperature, new_h, dh_dp, dh_dT); };
   auto s_lambda = [&](const T & pressure, const T & temperature, T & new_s, T & ds_dp, T & ds_dT)
   { s_from_p_T(pressure, temperature, new_s, ds_dp, ds_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(h,
                                         s,
                                         p0,
                                         T0,
                                         pressure,
                                         temperature,
                                         _tolerance,
                                         _tolerance,
                                         h_lambda,
                                         s_lambda,
                                         "p_T_from_h_s",
                                         _max_newton_its,
                                         _verbose_newton);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (h, s)=(", h, ", ", s, ") to (p, T) failed"));
 }

◆ p_T_from_v_e() [1/2]

template<typename CppType >

void SinglePhaseFluidProperties::p_T_from_v_e	(	const CppType &	v,
		const CppType &	e,
		Real	p0,
		Real	T0,
		CppType &	p,
		CppType &	T,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (v,e) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	v	specific volume (m^3 / kg)
[in]	e	specific internal energy (J / kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 513 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::g_from_v_e(), TemperaturePressureFunctionFluidProperties::k_from_v_e(), TemperaturePressureFunctionFluidProperties::mu_from_v_e(), and TemperaturePressureFunctionFluidProperties::T_from_v_e().

 {
   auto v_lambda = [&](const CppType & pressure,
                       const CppType & temperature,
                       CppType & new_v,
                       CppType & dv_dp,
                       CppType & dv_dT) { v_from_p_T(pressure, temperature, new_v, dv_dp, dv_dT); };
   auto e_lambda = [&](const CppType & pressure,
                       const CppType & temperature,
                       CppType & new_e,
                       CppType & de_dp,
                       CppType & de_dT) { e_from_p_T(pressure, temperature, new_e, de_dp, de_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(v,
                                         e,
                                         p0,
                                         T0,
                                         p,
                                         T,
                                         _tolerance,
                                         _tolerance,
                                         v_lambda,
                                         e_lambda,
                                         "p_T_from_v_e",
                                         _max_newton_its,
                                         _verbose_newton);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (v, e)=(", v, ", ", e, ") to (p, T) failed"));
 }

◆ p_T_from_v_e() [2/2]

template<typename T >

std::pair< T, T > Water97FluidProperties::p_T_from_v_e	(	const T &	v,
		const T &	e
	)		const

Computes the pressure (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific internal energy.

Definition at line 2204 of file Water97FluidProperties.h.

Referenced by c_from_v_e(), cp_from_v_e(), cv_from_v_e(), k_from_v_e_template(), and T_from_v_e().

 {
   const auto rho = 1 / v;
   const auto p = p_from_v_e(v, e);
   return {p, T_drhodT_from_p_rho(p, rho).first};
 }

◆ p_T_from_v_h() [1/3]

template<typename T >

void SinglePhaseFluidProperties::p_T_from_v_h	(	const T &	v,
		const T &	h,
		Real	p0,
		Real	T0,
		T &	pressure,
		T &	temperature,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	v	specific volume (m^3 / kg)
[in]	h	specific enthalpy (J / kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 559 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedBicubicFluidProperties::constructInterpolation().

 {
   auto v_lambda = [&](const T & pressure, const T & temperature, T & new_v, T & dv_dp, T & dv_dT)
   { v_from_p_T(pressure, temperature, new_v, dv_dp, dv_dT); };
   auto h_lambda = [&](const T & pressure, const T & temperature, T & new_h, T & dh_dp, T & dh_dT)
   { h_from_p_T(pressure, temperature, new_h, dh_dp, dh_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(v,
                                         h,
                                         p0,
                                         T0,
                                         pressure,
                                         temperature,
                                         _tolerance,
                                         _tolerance,
                                         v_lambda,
                                         h_lambda,
                                         "p_T_from_v_h",
                                         _max_newton_its,
                                         _verbose_newton);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (v, h)=(", v, ", ", h, ") to (p, T) failed"));
 }

◆ p_T_from_v_h() [2/3]

template<typename T >

std::pair< T, T > Water97FluidProperties::p_T_from_v_h	(	const T &	v,
		const T &	h
	)		const

Computes the pressure (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific enthalpy.

Definition at line 2507 of file Water97FluidProperties.h.

Referenced by e_from_v_h().

 {
   T pressure, temperature;
   bool conversion_succeeded;
   p_T_from_v_h(
       v, h, _p_initial_guess, _T_initial_guess, pressure, temperature, conversion_succeeded);
   return {std::move(pressure), std::move(temperature)};
 }

◆ p_T_from_v_h() [3/3]

template<typename T >

void SinglePhaseFluidProperties::p_T_from_v_h ( typename T )

Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	v	specific volume (m^3 / kg)
[in]	h	specific enthalpy (J / kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 559 of file SinglePhaseFluidProperties.h.

 {
   auto v_lambda = [&](const T & pressure, const T & temperature, T & new_v, T & dv_dp, T & dv_dT)
   { v_from_p_T(pressure, temperature, new_v, dv_dp, dv_dT); };
   auto h_lambda = [&](const T & pressure, const T & temperature, T & new_h, T & dh_dp, T & dh_dT)
   { h_from_p_T(pressure, temperature, new_h, dh_dp, dh_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(v,
                                         h,
                                         p0,
                                         T0,
                                         pressure,
                                         temperature,
                                         _tolerance,
                                         _tolerance,
                                         v_lambda,
                                         h_lambda,
                                         "p_T_from_v_h",
                                         _max_newton_its,
                                         _verbose_newton);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (v, h)=(", v, ", ", h, ") to (p, T) failed"));
 }

◆ phi3()

template<typename T >

T Water97FluidProperties::phi3	(	const T &	delta,
		const T &	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 1735 of file Water97FluidProperties.h.

 {
   using std::log;
 
   T 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] * log(delta) + sum;
 }

◆ propfunc() [1/12]

SinglePhaseFluidProperties::propfunc	(	p	,
		v	,
		e
	)

inherited

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

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

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

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

As an example:

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

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

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

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

◆ propfunc() [2/12]

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

inherited

◆ propfunc() [3/12]

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

inherited

◆ propfunc() [4/12]

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

inherited

◆ propfunc() [5/12]

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

inherited

◆ propfunc() [6/12]

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

inherited

◆ propfunc() [7/12]

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

inherited

◆ propfunc() [8/12]

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

inherited

◆ propfunc() [9/12]

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

inherited

◆ propfunc() [10/12]

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

inherited

◆ propfunc() [11/12]

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

inherited

◆ propfunc() [12/12]

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

inherited

◆ propfuncWithDefault() [1/7]

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

inherited

◆ propfuncWithDefault() [2/7]

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

inherited

◆ propfuncWithDefault() [3/7]

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

inherited

◆ propfuncWithDefault() [4/7]

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

inherited

◆ propfuncWithDefault() [5/7]

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

inherited

◆ propfuncWithDefault() [6/7]

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

inherited

◆ propfuncWithDefault() [7/7]

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

inherited

◆ propfuncWithDefinitionOverride()

Water97FluidProperties::propfuncWithDefinitionOverride	(	s	,
		p	,
		T
	)

◆ rho_e_from_p_T()

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

virtualinherited

Definition at line 438 of file SinglePhaseFluidProperties.C.

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

◆ rho_from_p_T() [1/4]

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

overridevirtual

Definition at line 81 of file Water97FluidProperties.C.

Referenced by k_from_p_T_template(), mu_from_p_T(), and rho_mu_from_p_T().

 {
   return rho_from_p_T_template(pressure, temperature);
 }

◆ rho_from_p_T() [2/4]

ADReal Water97FluidProperties::rho_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 87 of file Water97FluidProperties.C.

 {
   return rho_from_p_T_template(pressure, temperature);
 }

◆ rho_from_p_T() [3/4]

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

overridevirtual

Definition at line 93 of file Water97FluidProperties.C.

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

◆ rho_from_p_T() [4/4]

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

overridevirtual

Definition at line 100 of file Water97FluidProperties.C.

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

◆ rho_from_p_T_template() [1/2]

template<typename T >

T Water97FluidProperties::rho_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2575 of file Water97FluidProperties.h.

Referenced by mu_from_p_T_template(), rho_from_p_T(), rho_from_p_T_template(), T_drhodT_from_p_rho(), and v_from_p_T_template().

 {
   T density, pi, tau;
 
   // Determine which region the point is in
   unsigned int region =
       inRegion(MetaPhysicL::raw_value(pressure), MetaPhysicL::raw_value(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("inRegion() has given an incorrect region");
   }
   return density;
 }

◆ rho_from_p_T_template() [2/2]

template<typename T >

void Water97FluidProperties::rho_from_p_T_template	(	const T &	pressure,
		const T &	temperature,
		T &	rho,
		T &	drho_dp,
		T &	drho_dt
	)		const

Definition at line 2017 of file Water97FluidProperties.h.

 {
   auto functor = [this](const auto & pressure, const auto & temperature)
   { return this->rho_from_p_T_template(pressure, temperature); };
 
   xyDerivatives(pressure, temperature, rho, drho_dp, drho_dT, functor);
 }

◆ rho_mu_from_p_T() [1/3]

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

overridevirtual

Combined methods.

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 369 of file Water97FluidProperties.C.

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

◆ rho_mu_from_p_T() [2/3]

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

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 379 of file Water97FluidProperties.C.

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

◆ rho_mu_from_p_T() [3/3]

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

virtualinherited

Definition at line 473 of file SinglePhaseFluidProperties.C.

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

◆ rho_T_from_v_e()

template<typename T >

std::pair< T, T > Water97FluidProperties::rho_T_from_v_e	(	const T &	v,
		const T &	e
	)		const

Computes the density (first member of the pair) and temperature (second member of the pair) as functions of specific volume and specific internal energy.

Definition at line 2213 of file Water97FluidProperties.h.

Referenced by mu_from_v_e().

 {
   const auto rho = 1 / v;
   const auto p = p_from_v_e(v, e);
   return {rho, T_drhodT_from_p_rho(p, rho).first};
 }

◆ s_from_h_p() [1/3]

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

overridevirtual

Definition at line 450 of file Water97FluidProperties.C.

 {
   Real T = T_from_p_h(pressure, enthalpy);
   return s_from_p_T(pressure, T);
 }

◆ s_from_h_p() [2/3]

ADReal Water97FluidProperties::s_from_h_p	(	const ADReal &	enthalpy,
		const ADReal &	pressure
	)		const

overridevirtual

Definition at line 457 of file Water97FluidProperties.C.

 {
   ADReal temperature = T_from_p_h_ad(pressure, enthalpy);
   return s_from_p_T_template(pressure, temperature);
 }

◆ s_from_h_p() [3/3]

void Water97FluidProperties::s_from_h_p	(	Real	enthalpy,
		Real	pressure,
		Real &	s,
		Real &	ds_dh,
		Real &	ds_dp
	)		const

overridevirtual

Definition at line 464 of file Water97FluidProperties.C.

 {
   ADReal p = pressure;
   Moose::derivInsert(p.derivatives(), 0, 1.0);
 
   ADReal h = enthalpy;
   Moose::derivInsert(h.derivatives(), 1, 1.0);
 
   ADReal T = T_from_p_h_ad(p, h);
   ADReal entropy = s_from_p_T_template(p, T);
 
   ds_dh = entropy.derivatives()[1];
   ds_dp = entropy.derivatives()[0];
   s = entropy.value();
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 38 of file FluidProperties.h.

38 {}

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

Referenced by T_from_p_h_ad().

 {
   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 1034 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("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 692 of file Water97FluidProperties.C.

Referenced by T_from_p_h_ad().

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

◆ subregionVolume()

template<typename T >

T Water97FluidProperties::subregionVolume	(	const T &	pi,
		const T &	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 2173 of file Water97FluidProperties.h.

Referenced by densityRegion3().

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

◆ T_drhodT_from_p_rho()

template<typename T >

std::pair< T, T > Water97FluidProperties::T_drhodT_from_p_rho	(	const T &	p,
		const T &	rho
	)		const

private

Computes the temperature (first member of the pair) and the derivative of density (second member of the pair) with respect to temperature as a function of pressure and density.

Definition at line 2189 of file Water97FluidProperties.h.

Referenced by e_from_p_rho_template(), p_T_from_v_e(), and rho_T_from_v_e().

 {
   // For NewtonSolve
   // y = rho
   // x = p
   // z = T
   auto lambda =
       [this](const T & p, const T & temperature, T & rho, T & drho_dp, T & drho_dtemperature)
   { rho_from_p_T_template(p, temperature, rho, drho_dp, drho_dtemperature); };
   return FluidPropertiesUtils::NewtonSolve(
       p, rho, _T_initial_guess, _tolerance, lambda, name() + "::T_drhodT_from_p_rho");
 }

◆ T_from_p_h() [1/3]

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

overridevirtual

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 718 of file Water97FluidProperties.C.

Referenced by s_from_h_p().

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

◆ T_from_p_h() [2/3]

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

overridevirtual

Definition at line 727 of file Water97FluidProperties.C.

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

◆ T_from_p_h() [3/3]

ADReal Water97FluidProperties::T_from_p_h	(	const ADReal &	pressure,
		const ADReal &	enthalpy
	)		const

overridevirtual

Definition at line 743 of file Water97FluidProperties.C.

 {
   return T_from_p_h_ad(pressure, enthalpy);
 }

◆ T_from_p_h_ad()

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

protected

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

Parameters

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

Returns: temperature water temperature (K)

Definition at line 749 of file Water97FluidProperties.C.

Referenced by s_from_h_p(), and T_from_p_h().

 {
   ADReal temperature = 0.0;
 
   // Determine which region the point is in
   const unsigned int region = inRegionPH(pressure.value(), enthalpy.value());
 
   switch (region)
   {
     case 1:
       temperature = temperature_from_ph1(pressure, enthalpy);
       break;
 
     case 2:
     {
       // First, determine which subregion the point is in:
       const unsigned int subregion = subregion2ph(pressure.value(), enthalpy.value());
 
       if (subregion == 1)
         temperature = temperature_from_ph2a(pressure, enthalpy);
       else if (subregion == 2)
         temperature = temperature_from_ph2b(pressure, enthalpy);
       else
         temperature = temperature_from_ph2c(pressure, enthalpy);
       break;
     }
 
     case 3:
     {
       // First, determine which subregion the point is in:
       const unsigned int subregion = subregion3ph(pressure.value(), enthalpy.value());
 
       if (subregion == 1)
         temperature = temperature_from_ph3a(pressure, enthalpy);
       else
         temperature = temperature_from_ph3b(pressure, enthalpy);
       break;
     }
 
     case 4:
       return vaporTemperature_ad(pressure);
       break;
 
     case 5:
     {
       Real T_min = 1073.15;
       Real T_max = 2273.15;
       Real T_find = 0.;
       Real T_error = 1.0;
 
       // Bisection solve
       while (T_error > libMesh::TOLERANCE)
       {
         T_find = 0.5 * (T_min + T_max);
         Real h_find = h_from_p_T(pressure.value(), T_find);
 
         if (h_find > enthalpy.value())
           T_max = T_find;
         else
           T_min = T_find;
 
         T_error = std::abs((T_max - T_min) / T_find);
       }
       if (!_allow_imperfect_jacobians)
         mooseDoOnce(
             mooseWarning("The derivatives for T_from_p_h_ad are currently neglected, set to 0."));
       temperature = T_find;
 
       break;
     }
     default:
       mooseError("inRegionPH() has given an incorrect region");
   }
 
   return temperature;
 }

◆ T_from_v_e() [1/2]

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

override

Definition at line 213 of file Water97FluidProperties.C.

 {
   return p_T_from_v_e(v, e).second;
 }

◆ T_from_v_e() [2/2]

ADReal Water97FluidProperties::T_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

override

Definition at line 219 of file Water97FluidProperties.C.

 {
   return p_T_from_v_e(v, e).second;
 }

◆ temperature_from_ph1()

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

Referenced by T_from_p_h_ad().

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

◆ temperature_from_ph2a()

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

Referenced by T_from_p_h_ad().

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

◆ temperature_from_ph2b()

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

Referenced by T_from_p_h_ad().

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

◆ temperature_from_ph2c()

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

Referenced by T_from_p_h_ad().

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

◆ temperature_from_ph3a()

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

Referenced by T_from_p_h_ad().

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

◆ temperature_from_ph3b()

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

Referenced by T_from_p_h_ad().

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

◆ tempXY()

template<typename T >

T Water97FluidProperties::tempXY	(	const T &	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 1909 of file Water97FluidProperties.h.

Referenced by subregion3().

 {
   T 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;
   }
 
   T sum = 0.0;
 
   using std::log;
 
   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(log(pi), _tempXY_I[row][i]);
   else if (xy == EF)
     sum += 3.727888004 * (pi - _p_critical / 1.0e6) + _T_critical;
   else
     for (std::size_t i = 0; i < _tempXY_n[row].size(); ++i)
       sum += _tempXY_n[row][i] * MathUtils::pow(pi, _tempXY_I[row][i]);
 
   return sum;
 }

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 37 of file FluidProperties.h.

37 {}

◆ triplePointPressure()

Real Water97FluidProperties::triplePointPressure ( ) const

overridevirtual

Triple point pressure.

Returns: triple point pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 69 of file Water97FluidProperties.C.

 {
   return _p_triple;
 }

◆ triplePointTemperature()

Real Water97FluidProperties::triplePointTemperature ( ) const

overridevirtual

Triple point temperature.

Returns: triple point temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 75 of file Water97FluidProperties.C.

 {
   return _T_triple;
 }

◆ v_e_from_p_T() [1/2]

template<typename CppType >

void SinglePhaseFluidProperties::v_e_from_p_T	(	const CppType &	p,
		const CppType &	T,
		CppType &	v,
		CppType &	e
	)		const

inherited

Definition at line 639 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::h_from_p_T(), and TabulatedFluidProperties::s_from_p_T().

 {
   const CppType rho = rho_from_p_T(p, T);
   v = 1.0 / rho;
   try
   {
     // more likely to not involve a Newton search
     e = e_from_p_T(p, T);
   }
   catch (...)
   {
     e = e_from_p_rho(p, rho);
   }
 }

◆ v_e_from_p_T() [2/2]

template<typename CppType >

void SinglePhaseFluidProperties::v_e_from_p_T	(	const CppType &	p,
		const CppType &	T,
		CppType &	v,
		CppType &	dv_dp,
		CppType &	dv_dT,
		CppType &	e,
		CppType &	de_dp,
		CppType &	de_dT
	)		const

inherited

Definition at line 659 of file SinglePhaseFluidProperties.h.

 {
   CppType rho, drho_dp, drho_dT;
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
 
   v = 1.0 / rho;
   const CppType dv_drho = -1.0 / (rho * rho);
   dv_dp = dv_drho * drho_dp;
   dv_dT = dv_drho * drho_dT;
 
   CppType de_dp_partial, de_drho;
   e_from_p_rho(p, rho, e, de_dp_partial, de_drho);
   de_dp = de_dp_partial + de_drho * drho_dp;
   de_dT = de_drho * drho_dT;
 }

◆ v_e_spndl_from_T()

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

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 489 of file SinglePhaseFluidProperties.C.

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

◆ v_from_p_T() [1/4]

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

overridevirtual

Definition at line 110 of file Water97FluidProperties.C.

 {
   return v_from_p_T_template(pressure, temperature);
 }

◆ v_from_p_T() [2/4]

ADReal Water97FluidProperties::v_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature
	)		const

overridevirtual

Definition at line 116 of file Water97FluidProperties.C.

 {
   return v_from_p_T_template(pressure, temperature);
 }

◆ v_from_p_T() [3/4]

void Water97FluidProperties::v_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	v,
		Real &	dv_dp,
		Real &	dv_dT
	)		const

overridevirtual

Definition at line 122 of file Water97FluidProperties.C.

 {
   v_from_p_T_template(pressure, temperature, v, dv_dp, dv_dT);
 }

◆ v_from_p_T() [4/4]

void Water97FluidProperties::v_from_p_T	(	const ADReal &	pressure,
		const ADReal &	temperature,
		ADReal &	v,
		ADReal &	dv_dp,
		ADReal &	dv_dT
	)		const

overridevirtual

Definition at line 129 of file Water97FluidProperties.C.

 {
   v_from_p_T_template(pressure, temperature, v, dv_dp, dv_dT);
 }

◆ v_from_p_T_template() [1/2]

template<typename T >

T Water97FluidProperties::v_from_p_T_template	(	const T &	pressure,
		const T &	temperature
	)		const

Definition at line 2028 of file Water97FluidProperties.h.

Referenced by v_from_p_T(), and v_from_p_T_template().

 {
   return 1 / rho_from_p_T_template(pressure, temperature);
 }

◆ v_from_p_T_template() [2/2]

template<typename T >

void Water97FluidProperties::v_from_p_T_template	(	const T &	pressure,
		const T &	temperature,
		T &	v,
		T &	dv_dp,
		T &	dv_dt
	)		const

Definition at line 2035 of file Water97FluidProperties.h.

 {
   auto functor = [this](const auto & pressure, const auto & temperature)
   { return this->v_from_p_T_template(pressure, temperature); };
 
   xyDerivatives(pressure, temperature, v, dv_dp, dv_dT, functor);
 }

◆ validParams()

InputParameters Water97FluidProperties::validParams ( )

static

Definition at line 17 of file Water97FluidProperties.C.

 {
   InputParameters params = SinglePhaseFluidProperties::validParams();
   params.addClassDescription("Fluid properties for water and steam (H2O) using IAPWS-IF97");
   return params;
 }

◆ vaporPressure() [1/3]

Real Water97FluidProperties::vaporPressure ( Real T ) const

overridevirtual

Vapor pressure.

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

Parameters

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

Reimplemented from SinglePhaseFluidProperties.

Definition at line 521 of file Water97FluidProperties.C.

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

 {
   // Check whether the input temperature is within the region of validity of this equation.
   // Valid for 273.15 K <= t <= 647.096 K
   if (temperature < 273.15 || temperature > _T_critical)
     mooseException(name(),
                    ": vaporPressure(): Temperature ",
                    temperature,
                    " is outside range 273.15 K <= T "
                    "<= 647.096 K");
 
   Real theta, theta2, a, b, c, p;
   theta = temperature + _n4[8] / (temperature - _n4[9]);
   theta2 = theta * theta;
   a = theta2 + _n4[0] * theta + _n4[1];
   b = _n4[2] * theta2 + _n4[3] * theta + _n4[4];
   c = _n4[5] * theta2 + _n4[6] * theta + _n4[7];
   p = Utility::pow<4>(2.0 * c / (-b + std::sqrt(b * b - 4.0 * a * c)));
 
   return p * 1.e6;
 }

◆ vaporPressure() [2/3]

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

inlineoverridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 2010 of file Water97FluidProperties.h.

 {
   vaporPressureTemplate(temperature, psat, dpsat_dT);
 }

◆ vaporPressure() [3/3]

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

virtualinherited

Definition at line 393 of file SinglePhaseFluidProperties.C.

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

◆ vaporPressureTemplate()

template<typename T >

void Water97FluidProperties::vaporPressureTemplate	(	const T &	temperature,
		T &	psat,
		T &	dpsat_dT
	)		const

Definition at line 1973 of file Water97FluidProperties.h.

Referenced by vaporPressure().

 {
   using std::sqrt, std::pow;
 
   // Check whether the input temperature is within the region of validity of this equation.
   // Valid for 273.15 K <= t <= 647.096 K
   if (temperature < 273.15 || temperature > _T_critical)
     mooseException(name(),
                    ": vaporPressure(): Temperature is outside range 273.15 K <= T <= 647.096 K");
 
   T theta, dtheta_dT, theta2, a, b, c, da_dtheta, db_dtheta, dc_dtheta;
   theta = temperature + _n4[8] / (temperature - _n4[9]);
   dtheta_dT = 1.0 - _n4[8] / (temperature - _n4[9]) / (temperature - _n4[9]);
   theta2 = theta * theta;
 
   a = theta2 + _n4[0] * theta + _n4[1];
   b = _n4[2] * theta2 + _n4[3] * theta + _n4[4];
   c = _n4[5] * theta2 + _n4[6] * theta + _n4[7];
 
   da_dtheta = 2.0 * theta + _n4[0];
   db_dtheta = 2.0 * _n4[2] * theta + _n4[3];
   dc_dtheta = 2.0 * _n4[5] * theta + _n4[6];
 
   T denominator = -b + sqrt(b * b - 4.0 * a * c);
 
   psat = Utility::pow<4>(2.0 * c / denominator) * 1.0e6;
 
   // The derivative wrt temperature is given by the chain rule
   T dpsat = 4.0 * Utility::pow<3>(2.0 * c / denominator);
   dpsat *= (2.0 * dc_dtheta / denominator -
             2.0 * c / denominator / denominator *
                 (-db_dtheta + pow(b * b - 4.0 * a * c, -0.5) *
                                   (b * db_dtheta - 2.0 * da_dtheta * c - 2.0 * a * dc_dtheta)));
   dpsat_dT = dpsat * dtheta_dT * 1.0e6;
 }

◆ vaporTemperature() [1/3]

Real Water97FluidProperties::vaporTemperature ( Real pressure ) const

overridevirtual

Saturation temperature as a function of pressure.

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

Valid for 611.213 Pa <= p <= 22.064 MPa

Parameters

pressure water pressure (Pa)

Returns: saturation temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 565 of file Water97FluidProperties.C.

Referenced by inRegionPH(), and subregion3().

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

◆ vaporTemperature() [2/3]

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

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 573 of file Water97FluidProperties.C.

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

◆ vaporTemperature() [3/3]

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

virtualinherited

Definition at line 423 of file SinglePhaseFluidProperties.C.

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

◆ vaporTemperature_ad()

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

protected

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

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

Valid for 611.213 Pa <= p <= 22.064 MPa

Parameters

pressure water pressure (Pa)

Returns: saturation temperature (K)

Definition at line 544 of file Water97FluidProperties.C.

Referenced by T_from_p_h_ad(), and vaporTemperature().

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

◆ xyDerivatives()

template<typename T , typename Functor >

void SinglePhaseFluidProperties::xyDerivatives	(	const T	x,
		const T &	y,
		T &	z,
		T &	dz_dx,
		T &	dz_dy,
		const Functor &	z_from_x_y
	)

staticprotectedinherited

Computes the dependent variable z and its derivatives with respect to the independent variables x and y using the simple two parameter z_from_x_y functor.

The derivatives are computed using a compound automatic differentiation type

Definition at line 492 of file SinglePhaseFluidProperties.h.

Referenced by e_from_p_rho_template(), h_from_p_T_template(), rho_from_p_T_template(), and v_from_p_T_template().

 {
   typedef MetaPhysicL::DualNumber<T, MetaPhysicL::NumberArray<2, T>> CompoundType;
   const auto [zero, one] = makeZeroAndOne(x);
 
   CompoundType x_c(x, zero);
   auto & x_cd = x_c.derivatives();
   x_cd[0] = one;
   CompoundType y_c(y, zero);
   auto & y_cd = y_c.derivatives();
   y_cd[1] = one;
 
   const auto z_c = z_from_x_y(x_c, y_c);
   z = z_c.value();
   dz_dx = z_c.derivatives()[0];
   dz_dy = z_c.derivatives()[1];
 }

Member Data Documentation

◆ _allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians

protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 47 of file FluidProperties.h.

Referenced by VaporMixtureFluidProperties::imperfectJacobianMessage(), T_from_p_h_ad(), and SinglePhaseFluidProperties::unimplementedDerivativeMethod().

◆ _I1

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

protected

Initial value:

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

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

Definition at line 856 of file Water97FluidProperties.h.

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

◆ _I2

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

protected

Initial value:

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

Definition at line 901 of file Water97FluidProperties.h.

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

◆ _I3

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

protected

Initial value:

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

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

Definition at line 980 of file Water97FluidProperties.h.

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

◆ _I3s

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

protected

Definition at line 1266 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _I5

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

protected

Definition at line 1425 of file Water97FluidProperties.h.

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

◆ _Iph1

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

protected

Definition at line 870 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _Iph2a

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

protected

Initial value:

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

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

Definition at line 920 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _Iph2b

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

protected

Initial value:

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

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

Definition at line 939 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _Iph2c

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

protected

Initial value:

{

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

Definition at line 954 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _Iph3a

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

protected

Initial value:

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

Definition at line 1379 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _Iph3b

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

protected

Initial value:

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

Definition at line 1397 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _J02

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

protected

Definition at line 899 of file Water97FluidProperties.h.

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

◆ _J05

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

protected

Constants for region 5.

Definition at line 1417 of file Water97FluidProperties.h.

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

◆ _J1

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

protected

Initial value:

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

Definition at line 859 of file Water97FluidProperties.h.

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

◆ _J2

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

protected

Initial value:

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

Definition at line 905 of file Water97FluidProperties.h.

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

◆ _J3

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

protected

Initial value:

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

Definition at line 983 of file Water97FluidProperties.h.

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

◆ _J3s

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

protected

Definition at line 1318 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _J5

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

protected

Definition at line 1426 of file Water97FluidProperties.h.

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

◆ _Jph1

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

protected

Initial value:

{

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

Definition at line 872 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _Jph2a

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

protected

Initial value:

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

Definition at line 923 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _Jph2b

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

protected

Initial value:

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

Definition at line 942 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _Jph2c

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

protected

Initial value:

{

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

Definition at line 957 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _Jph3a

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

protected

Initial value:

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

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

Definition at line 1383 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _Jph3b

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

protected

Initial value:

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

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

Definition at line 1401 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _k_a

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

protected

Constants for thermal conductivity.

Definition at line 1496 of file Water97FluidProperties.h.

Referenced by k_from_rho_T_template().

◆ _max_newton_its

const unsigned int SinglePhaseFluidProperties::_max_newton_its

protectedinherited

Maximum number of iterations for the variable conversion newton solves.

Definition at line 449 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::e_from_p_rho(), TabulatedFluidProperties::e_from_p_T(), TabulatedFluidProperties::e_from_v_h(), SinglePhaseFluidProperties::p_T_from_h_s(), SinglePhaseFluidProperties::p_T_from_v_e(), SinglePhaseFluidProperties::p_T_from_v_h(), TemperaturePressureFunctionFluidProperties::T_from_p_h(), SimpleFluidProperties::T_from_p_h(), TabulatedFluidProperties::T_from_p_h(), NaKFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_s(), and TabulatedFluidProperties::v_from_p_T().

◆ _Mh2o

const Real Water97FluidProperties::_Mh2o

protected

Water molar mass (kg/mol)

Definition at line 822 of file Water97FluidProperties.h.

Referenced by molarMass().

◆ _mu_H0

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

protected

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

Definition at line 1486 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T(), and mu_from_rho_T_template().

◆ _mu_Hij

const std::array<std::array<Real, 7>, 6> Water97FluidProperties::_mu_Hij

protected

Initial value:

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

Definition at line 1487 of file Water97FluidProperties.h.

Referenced by mu_from_rho_T(), and mu_from_rho_T_template().

◆ _n02

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

protected

Initial value:

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

Constants for region 2.

Definition at line 876 of file Water97FluidProperties.h.

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

◆ _n05

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

protected

Initial value:

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

Definition at line 1418 of file Water97FluidProperties.h.

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

◆ _n1

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

protected

Initial value:

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

Reference constants used in to calculate thermophysical properties of water.

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

Definition at line 845 of file Water97FluidProperties.h.

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

◆ _n2

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

protected

Initial value:

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

Definition at line 886 of file Water97FluidProperties.h.

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

◆ _n23

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

protected

Initial value:

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

Constants for the boundary between regions 2 and 3.

Definition at line 961 of file Water97FluidProperties.h.

Referenced by b23p(), and b23T().

◆ _n3

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

protected

Initial value:

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

Constants for region 3.

Definition at line 968 of file Water97FluidProperties.h.

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

◆ _n3s

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

protected

Constants for all 26 subregions in region 3.

Definition at line 1019 of file Water97FluidProperties.h.

Referenced by densityRegion3(), and subregionVolume().

◆ _n4

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

protected

Initial value:

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

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

Definition at line 1405 of file Water97FluidProperties.h.

Referenced by vaporPressure(), vaporPressureTemplate(), and vaporTemperature_ad().

◆ _n5

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

protected

Initial value:

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

Definition at line 1427 of file Water97FluidProperties.h.

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

◆ _nph1

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

protected

Initial value:

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

Definition at line 863 of file Water97FluidProperties.h.

Referenced by temperature_from_ph1().

◆ _nph2a

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

protected

Initial value:

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

Definition at line 909 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2a().

◆ _nph2b

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

protected

Initial value:

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

Definition at line 927 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2b().

◆ _nph2c

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

protected

Initial value:

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

Definition at line 946 of file Water97FluidProperties.h.

Referenced by temperature_from_ph2c().

◆ _nph3a

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

protected

Initial value:

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

Definition at line 1369 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3a().

◆ _nph3b

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

protected

Initial value:

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

Definition at line 1386 of file Water97FluidProperties.h.

Referenced by temperature_from_ph3b().

◆ _p_critical

const Real Water97FluidProperties::_p_critical

protected

Critical pressure (Pa)

Definition at line 826 of file Water97FluidProperties.h.

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

◆ _p_initial_guess

const Real SinglePhaseFluidProperties::_p_initial_guess

protectedinherited

Initial guess for pressure (or pressure used to compute the initial guess)

Definition at line 447 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::checkInitialGuess(), TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::e_from_v_h(), TabulatedFluidProperties::g_from_v_e(), TemperaturePressureFunctionFluidProperties::k_from_v_e(), TemperaturePressureFunctionFluidProperties::mu_from_v_e(), TemperaturePressureFunctionFluidProperties::p_from_v_e(), p_from_v_e_template(), p_T_from_v_h(), TabulatedFluidProperties::T_from_h_s(), and TemperaturePressureFunctionFluidProperties::T_from_v_e().

◆ _p_star

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

protected

Pressure scale for each region.

Definition at line 1501 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), inRegionPH(), and rho_from_p_T_template().

◆ _p_triple

const Real Water97FluidProperties::_p_triple

protected

Triple point pressure (Pa)

Definition at line 832 of file Water97FluidProperties.h.

Referenced by triplePointPressure().

◆ _par3

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

protected

Initial value:

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

Definition at line 987 of file Water97FluidProperties.h.

Referenced by densityRegion3().

◆ _par3N

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

protected

Initial value:

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

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

Definition at line 1015 of file Water97FluidProperties.h.

Referenced by densityRegion3().

◆ _R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 41 of file FluidProperties.h.

Referenced by HelmholtzFluidProperties::c_from_p_T(), IdealGasFluidProperties::c_from_p_T(), MassDiffusionEnergyGasMixDGKernel::computeComponentTotalEnthalpy(), BinaryDiffusionCoefMaterial::computeQpProperties(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), LeadLithiumFluidProperties::mu_from_p_T(), LeadLithiumFluidProperties::mu_from_v_e(), HelmholtzFluidProperties::p_from_rho_T(), HelmholtzFluidProperties::rho_from_p_T(), CaloricallyImperfectGas::rho_from_p_T(), IdealGasFluidProperties::rho_from_p_T(), HelmholtzFluidProperties::s_from_p_T(), and TEST_F().

◆ _rho_critical

const Real Water97FluidProperties::_rho_critical

protected

Critical density (kg/m^3)

Definition at line 830 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), criticalDensity(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), mu_from_rho_T(), and mu_from_rho_T_template().

◆ _Rw

const Real Water97FluidProperties::_Rw

protected

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

Definition at line 824 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), inRegionPH(), and rho_from_p_T_template().

◆ _T_c2k

const Real FluidProperties::_T_c2k

protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 45 of file FluidProperties.h.

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

◆ _T_critical

const Real Water97FluidProperties::_T_critical

protected

Critical temperature (K)

Definition at line 828 of file Water97FluidProperties.h.

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

◆ _T_initial_guess

const Real SinglePhaseFluidProperties::_T_initial_guess

protectedinherited

Initial guess for temperature (or temperature used to compute the initial guess)

Definition at line 445 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::checkInitialGuess(), TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::e_from_p_rho(), TabulatedFluidProperties::g_from_v_e(), TemperaturePressureFunctionFluidProperties::k_from_v_e(), TemperaturePressureFunctionFluidProperties::mu_from_v_e(), p_T_from_v_h(), T_drhodT_from_p_rho(), TabulatedFluidProperties::T_from_h_s(), HelmholtzFluidProperties::T_from_p_h(), TemperaturePressureFunctionFluidProperties::T_from_p_h(), TabulatedFluidProperties::T_from_p_h(), NaKFluidProperties::T_from_p_rho(), TemperaturePressureFunctionFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_s(), and TemperaturePressureFunctionFluidProperties::T_from_v_e().

◆ _T_star

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

protected

Temperature scale for each region.

Definition at line 1499 of file Water97FluidProperties.h.

Referenced by c_from_p_T_template(), cp_from_p_T_template(), cv_from_p_T_template(), e_from_p_T_template(), h_from_p_T_template(), inRegionPH(), and rho_from_p_T_template().

◆ _T_triple

const Real Water97FluidProperties::_T_triple

protected

Triple point temperature (K)

Definition at line 834 of file Water97FluidProperties.h.

Referenced by triplePointTemperature().

◆ _tempXY_I

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

protected

Initial value:

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

Constnats for the tempXY() method.

Definition at line 1435 of file Water97FluidProperties.h.

Referenced by tempXY().

◆ _tempXY_n

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

protected

Definition at line 1447 of file Water97FluidProperties.h.

Referenced by tempXY().

◆ _tolerance

const Real SinglePhaseFluidProperties::_tolerance

protectedinherited

Newton's method may be used to convert between variable sets.

Relative tolerance of the solves

Definition at line 443 of file SinglePhaseFluidProperties.h.

◆ _verbose_newton

const bool SinglePhaseFluidProperties::_verbose_newton

protectedinherited

Whether to output information about newton solves to console.

Definition at line 451 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::e_from_p_rho(), TabulatedFluidProperties::e_from_p_T(), TabulatedFluidProperties::e_from_v_h(), SinglePhaseFluidProperties::p_T_from_h_s(), SinglePhaseFluidProperties::p_T_from_v_e(), SinglePhaseFluidProperties::p_T_from_v_h(), TabulatedFluidProperties::T_from_p_h(), TabulatedFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_s(), and TabulatedFluidProperties::v_from_p_T().

◆ h

e e e e s T T T T T rho v v T e SinglePhaseFluidProperties::h

inherited

Definition at line 212 of file SinglePhaseFluidProperties.h.

◆ p [1/10]

e e e e SinglePhaseFluidProperties::p

inherited

Definition at line 190 of file SinglePhaseFluidProperties.h.

◆ p [2/10]

e e e e s SinglePhaseFluidProperties::p

inherited

Definition at line 192 of file SinglePhaseFluidProperties.h.

◆ p [3/10]

e e e e s T T SinglePhaseFluidProperties::p

inherited

Definition at line 196 of file SinglePhaseFluidProperties.h.

◆ p [4/10]

e e e e s T T T SinglePhaseFluidProperties::p

inherited

Definition at line 198 of file SinglePhaseFluidProperties.h.

◆ p [5/10]

e e e e s T T T T SinglePhaseFluidProperties::p

inherited

Definition at line 200 of file SinglePhaseFluidProperties.h.

◆ p [6/10]

e e e e s T T T T T SinglePhaseFluidProperties::p

inherited

Definition at line 202 of file SinglePhaseFluidProperties.h.

◆ p [7/10]

e e e e s T T T T T rho v v SinglePhaseFluidProperties::p

inherited

Definition at line 208 of file SinglePhaseFluidProperties.h.

◆ p [8/10]

e e e e s T T T T T rho v v T e p SinglePhaseFluidProperties::p

inherited

Definition at line 214 of file SinglePhaseFluidProperties.h.

◆ p [9/10]

e e e e s T T T T T rho v v T e p T SinglePhaseFluidProperties::p

inherited

Definition at line 216 of file SinglePhaseFluidProperties.h.

◆ p [10/10]

e e e e s T T T T T rho v v T e p T T SinglePhaseFluidProperties::p

inherited

Definition at line 218 of file SinglePhaseFluidProperties.h.

◆ rho

e e e e s T SinglePhaseFluidProperties::rho

inherited

Definition at line 194 of file SinglePhaseFluidProperties.h.

◆ T [1/2]

e e e e s T T T T T rho SinglePhaseFluidProperties::T

inherited

Definition at line 204 of file SinglePhaseFluidProperties.h.

◆ T [2/2]

e e e e s T T T T T rho v SinglePhaseFluidProperties::T

inherited

Definition at line 206 of file SinglePhaseFluidProperties.h.

◆ v [1/5]

SinglePhaseFluidProperties::v

inherited

Definition at line 182 of file SinglePhaseFluidProperties.h.

◆ v [2/5]

e SinglePhaseFluidProperties::v

inherited

Definition at line 184 of file SinglePhaseFluidProperties.h.

◆ v [3/5]

e e SinglePhaseFluidProperties::v

inherited

Definition at line 186 of file SinglePhaseFluidProperties.h.

◆ v [4/5]

e e e SinglePhaseFluidProperties::v

inherited

Definition at line 188 of file SinglePhaseFluidProperties.h.

◆ v [5/5]

e e e e s T T T T T rho v v T SinglePhaseFluidProperties::v

inherited

Definition at line 210 of file SinglePhaseFluidProperties.h.

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

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

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Detailed Description

Member Enumeration Documentation

◆ subregionEnum

Constructor & Destructor Documentation

◆ Water97FluidProperties()

◆ ~Water97FluidProperties()

Member Function Documentation

◆ b23p()

◆ b23T()

◆ b2bc()

◆ b3ab()

◆ c_from_p_T() [1/2]

◆ c_from_p_T() [2/2]

◆ c_from_p_T_template()

◆ c_from_v_e()

◆ cp_from_p_T() [1/2]

◆ cp_from_p_T() [2/2]

◆ cp_from_p_T_template()

◆ cp_from_v_e()

◆ criticalDensity()

◆ criticalInternalEnergy()

◆ criticalPressure()

◆ criticalTemperature()

◆ cv_from_p_T() [1/2]

◆ cv_from_p_T() [2/2]

◆ cv_from_p_T_template()

◆ 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_from_p_rho() [1/4]

◆ e_from_p_rho() [2/4]

◆ e_from_p_rho() [3/4]

◆ e_from_p_rho() [4/4]

◆ e_from_p_rho_template() [1/2]

◆ e_from_p_rho_template() [2/2]

◆ e_from_p_T() [1/4]

◆ e_from_p_T() [2/4]

◆ e_from_p_T() [3/4]

◆ e_from_p_T() [4/4]

◆ e_from_p_T_template() [1/2]

◆ e_from_p_T_template() [2/2]

◆ e_from_v_h()

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ gamma1()

◆ gamma2()

◆ gamma5()

◆ h_from_p_T() [1/4]