Low level class for the real description of water. More...

#include <WaterPropsIAPWSphi.h>

Detailed Description

Low level class for the real description of water.

This is a helper class for WaterSSTP and PDSS_Water and does not constitute a complete implementation of a thermo phase by itself (see Thermodynamic Properties and classes WaterSSTP and PDSS_Water).

The reference is Wagner and Pruss [47].

Units Note: This class works with reduced units exclusively.

Definition at line 32 of file WaterPropsIAPWSphi.h.

Public Member Functions
	WaterPropsIAPWSphi ()
	Base constructor. More...

double	phi (double tau, double delta)
	Calculate the Phi function, which is the base function. More...

double	phi_d (double tau, double delta)
	Calculate derivative of phi wrt delta. More...

double	phi_dd (double tau, double delta)
	2nd derivative of phi wrt delta More...

double	phi_t (double tau, double delta)
	First derivative of phi wrt tau. More...

double	phi_tt (double tau, double delta)
	Second derivative of phi wrt tau. More...

double	pressureM_rhoRT (double tau, double delta)
	Calculate the dimensionless pressure at tau and delta;. More...

double	dimdpdrho (double tau, double delta)
	Dimensionless derivative of p wrt rho at constant T. More...

double	dimdpdT (double tau, double delta)
	Dimensionless derivative of p wrt T at constant rho. More...

double	dfind (double p_red, double tau, double deltaGuess)
	This function computes the reduced density, given the reduced pressure and the reduced temperature, tau. More...

double	gibbs_RT () const
	Calculate the dimensionless Gibbs free energy. More...

double	enthalpy_RT () const
	Calculate the dimensionless enthalpy, h/RT. More...

double	entropy_R () const
	Calculate the dimensionless entropy, s/R. More...

double	intEnergy_RT () const
	Calculate the dimensionless internal energy, u/RT. More...

double	cv_R () const
	Calculate the dimensionless constant volume heat capacity, Cv/R. More...

double	cp_R () const
	Calculate the dimensionless constant pressure heat capacity, Cv/R. More...

void	tdpolycalc (double tau, double delta)
	Calculates internal polynomials in tau and delta. More...

double	phiR () const
	Calculate Equation 6.6 for phiR, the residual part of the dimensionless Helmholtz free energy. More...

Protected Member Functions
double	phi0 () const
	Calculate Equation 6.5 for phi0, the ideal gas part of the dimensionless Helmholtz free energy. More...

double	phiR_d () const
	Calculate d_phiR_d(delta), the first derivative of phiR wrt delta. More...

double	phi0_d () const
	Calculate d_phi0_d(delta), the first derivative of phi0 wrt delta. More...

double	phiR_dd () const
	Calculate d2_phiR_dd(delta), the second derivative of phiR wrt delta. More...

double	phi0_dd () const
	Calculate d2_phi0_dd(delta), the second derivative of phi0 wrt delta. More...

double	phi0_t () const
	Calculate d_phi0/d(tau) More...

double	phiR_t () const
	Calculate Equation 6.6 for dphiRdtau, the derivative residual part of the dimensionless Helmholtz free energy wrt temperature. More...

double	phiR_tt () const
	Calculate Equation 6.6 for dphiRdtau, the second derivative residual part of the dimensionless Helmholtz free energy wrt temperature. More...

double	phi0_tt () const
	Calculate d2_phi0/dtau2. More...

double	phiR_dt () const
	Calculate the mixed derivative d2_phiR/(dtau ddelta) More...

double	phi0_dt () const
	Calculate the mixed derivative d2_phi0/(dtau ddelta) More...

Protected Attributes
double	TAUp [52]
	Value of internally calculated polynomials of powers of TAU. More...

double	DELTAp [16]
	Value of internally calculated polynomials of powers of delta. More...

double	TAUsave = -1.0
	Last tau that was used to calculate polynomials. More...

double	TAUsqrt = -1.0
	sqrt of TAU More...

double	DELTAsave = -1.0
	Last delta that was used to calculate polynomials. More...

Constructor & Destructor Documentation

◆ WaterPropsIAPWSphi()

WaterPropsIAPWSphi ( )

Base constructor.

The added constants were calculated so that u = s = 0 for liquid at the triple point.

These where determined by the program testPress. I'm not quite satisfied with the result, but will let it stand for the moment. H didn't turn out to be .611872 J/kg, but .611782 J/kg. There may be a slight error here somehow.

Definition at line 351 of file WaterPropsIAPWSphi.cpp.

Member Function Documentation

◆ phi()

double phi	(	double	tau,
		double	delta
	)

Calculate the Phi function, which is the base function.

The phi function is basically the Helmholtz free energy Eqn. (6.4) All internal polynomials are recalculated.

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 440 of file WaterPropsIAPWSphi.cpp.

◆ phi_d()

double phi_d	(	double	tau,
		double	delta
	)

Calculate derivative of phi wrt delta.

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 512 of file WaterPropsIAPWSphi.cpp.

◆ phi_dd()

double phi_dd	(	double	tau,
		double	delta
	)

2nd derivative of phi wrt delta

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 616 of file WaterPropsIAPWSphi.cpp.

◆ phi_t()

double phi_t	(	double	tau,
		double	delta
	)

First derivative of phi wrt tau.

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 702 of file WaterPropsIAPWSphi.cpp.

◆ phi_tt()

double phi_tt	(	double	tau,
		double	delta
	)

Second derivative of phi wrt tau.

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 783 of file WaterPropsIAPWSphi.cpp.

◆ pressureM_rhoRT()

double pressureM_rhoRT	(	double	tau,
		double	delta
	)

Calculate the dimensionless pressure at tau and delta;.

  pM/(rhoRT) = delta * phi_d() = 1.0 + delta phiR_d()

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

note: this is done so much, we have a separate routine.

Definition at line 520 of file WaterPropsIAPWSphi.cpp.

◆ dimdpdrho()

double dimdpdrho	(	double	tau,
		double	delta
	)

Dimensionless derivative of p wrt rho at constant T.

dp/drho * 1/RT = (2. delta phi_d() + delta**2 phi_dd()) (1.0 + 2. delta phiR_d() + delta**2 phiR_dd())

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 624 of file WaterPropsIAPWSphi.cpp.

◆ dimdpdT()

double dimdpdT	(	double	tau,
		double	delta
	)

Dimensionless derivative of p wrt T at constant rho.

dp/dT * M/(Rho R) = (1.0 + delta phiR_d()

tau delta (phiR_dt())

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 632 of file WaterPropsIAPWSphi.cpp.

◆ dfind()

double dfind	(	double	p_red,
		double	tau,
		double	deltaGuess
	)

This function computes the reduced density, given the reduced pressure and the reduced temperature, tau.

It takes an initial guess, deltaGuess. DeltaGuess is important as this is a multivalued function below the critical point.

Parameters

p_red	Value of the dimensionless pressure
tau	Dimensionless temperature = T_c/T
deltaGuess	Initial guess for the dimensionless density

Returns: the dimensionless density.

Definition at line 864 of file WaterPropsIAPWSphi.cpp.

◆ gibbs_RT()

double gibbs_RT ( ) const

Calculate the dimensionless Gibbs free energy.

Definition at line 939 of file WaterPropsIAPWSphi.cpp.

◆ enthalpy_RT()

double enthalpy_RT ( ) const

Calculate the dimensionless enthalpy, h/RT.

Definition at line 946 of file WaterPropsIAPWSphi.cpp.

◆ entropy_R()

double entropy_R ( ) const

Calculate the dimensionless entropy, s/R.

Definition at line 956 of file WaterPropsIAPWSphi.cpp.

◆ intEnergy_RT()

double intEnergy_RT ( ) const

Calculate the dimensionless internal energy, u/RT.

Definition at line 966 of file WaterPropsIAPWSphi.cpp.

◆ cv_R()

double cv_R ( ) const

Calculate the dimensionless constant volume heat capacity, Cv/R.

Definition at line 974 of file WaterPropsIAPWSphi.cpp.

◆ cp_R()

double cp_R ( ) const

Calculate the dimensionless constant pressure heat capacity, Cv/R.

Definition at line 982 of file WaterPropsIAPWSphi.cpp.

◆ tdpolycalc()

void tdpolycalc	(	double	tau,
		double	delta
	)

Calculates internal polynomials in tau and delta.

This routine is used to store the internal state of tau and delta for later use by the other routines in the class.

Parameters

tau	Dimensionless temperature = T_c/T
delta	Dimensionless density = delta = rho / Rho_c

Definition at line 361 of file WaterPropsIAPWSphi.cpp.

◆ phiR()

double phiR ( ) const

Calculate Equation 6.6 for phiR, the residual part of the dimensionless Helmholtz free energy.

Definition at line 394 of file WaterPropsIAPWSphi.cpp.

◆ phi0()

double phi0 ( ) const

protected

Calculate Equation 6.5 for phi0, the ideal gas part of the dimensionless Helmholtz free energy.

Definition at line 380 of file WaterPropsIAPWSphi.cpp.

◆ phiR_d()

double phiR_d ( ) const

protected

Calculate d_phiR_d(delta), the first derivative of phiR wrt delta.

Definition at line 448 of file WaterPropsIAPWSphi.cpp.

◆ phi0_d()

double phi0_d ( ) const

protected

Calculate d_phi0_d(delta), the first derivative of phi0 wrt delta.

Definition at line 506 of file WaterPropsIAPWSphi.cpp.

◆ phiR_dd()

double phiR_dd ( ) const

protected

Calculate d2_phiR_dd(delta), the second derivative of phiR wrt delta.

Definition at line 527 of file WaterPropsIAPWSphi.cpp.

◆ phi0_dd()

double phi0_dd ( ) const

protected

Calculate d2_phi0_dd(delta), the second derivative of phi0 wrt delta.

Definition at line 610 of file WaterPropsIAPWSphi.cpp.

◆ phi0_t()

double phi0_t ( ) const

protected

Calculate d_phi0/d(tau)

Definition at line 640 of file WaterPropsIAPWSphi.cpp.

◆ phiR_t()

double phiR_t ( ) const

protected

Calculate Equation 6.6 for dphiRdtau, the derivative residual part of the dimensionless Helmholtz free energy wrt temperature.

Definition at line 652 of file WaterPropsIAPWSphi.cpp.

◆ phiR_tt()

double phiR_tt ( ) const

protected

Calculate Equation 6.6 for dphiRdtau, the second derivative residual part of the dimensionless Helmholtz free energy wrt temperature.

Definition at line 723 of file WaterPropsIAPWSphi.cpp.

◆ phi0_tt()

double phi0_tt ( ) const

protected

Calculate d2_phi0/dtau2.

Definition at line 710 of file WaterPropsIAPWSphi.cpp.

◆ phiR_dt()

double phiR_dt ( ) const

protected

Calculate the mixed derivative d2_phiR/(dtau ddelta)

Definition at line 796 of file WaterPropsIAPWSphi.cpp.

◆ phi0_dt()

double phi0_dt ( ) const

protected

Calculate the mixed derivative d2_phi0/(dtau ddelta)

Definition at line 791 of file WaterPropsIAPWSphi.cpp.

Member Data Documentation

◆ TAUp

double TAUp[52]

protected

Value of internally calculated polynomials of powers of TAU.

Definition at line 183 of file WaterPropsIAPWSphi.h.

◆ DELTAp

double DELTAp[16]

protected

Value of internally calculated polynomials of powers of delta.

Definition at line 186 of file WaterPropsIAPWSphi.h.

◆ TAUsave

double TAUsave = -1.0

protected

Last tau that was used to calculate polynomials.

Definition at line 189 of file WaterPropsIAPWSphi.h.

◆ TAUsqrt

double TAUsqrt = -1.0

protected

sqrt of TAU

Definition at line 192 of file WaterPropsIAPWSphi.h.

◆ DELTAsave

double DELTAsave = -1.0

protected

Last delta that was used to calculate polynomials.

Definition at line 195 of file WaterPropsIAPWSphi.h.

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

Detailed Description

Public Member Functions

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ WaterPropsIAPWSphi()

Member Function Documentation

◆ phi()

◆ phi_d()

◆ phi_dd()

◆ phi_t()

◆ phi_tt()

◆ pressureM_rhoRT()

◆ dimdpdrho()

◆ dimdpdT()

◆ dfind()

◆ gibbs_RT()

◆ enthalpy_RT()

◆ entropy_R()

◆ intEnergy_RT()

◆ cv_R()

◆ cp_R()

◆ tdpolycalc()

◆ phiR()

◆ phi0()

◆ phiR_d()

◆ phi0_d()

◆ phiR_dd()

◆ phi0_dd()

◆ phi0_t()

◆ phiR_t()

◆ phiR_tt()

◆ phi0_tt()

◆ phiR_dt()

◆ phi0_dt()

Member Data Documentation

◆ TAUp

◆ DELTAp

◆ TAUsave

◆ TAUsqrt

◆ DELTAsave