WaterSteamPro functions list

Category: General

Recommended to use functions.

1.      Specific isobaric heat capacity at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspCPEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

2.      Specific isobaric heat capacity at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspCPEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

1.      Specific isobaric heat capacity [J/(kgK)] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspCPHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

2.      Specific isobaric heat capacity [J/(kgK)] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspCPPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspCPxPT or wspCPSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

3.      Specific isobaric heat capacity [J/(kgK)] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspCPPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspCPxPT or wspCPSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

4.      Specific isobaric heat capacity [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCPPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspCPxPT) is used.

5.      Specific isobaric heat capacity [J/(kgK)] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspCPPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspCPxPT or wspCPSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

6.      Specific isochoric heat capacity at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspCVEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

7.      Specific isochoric heat capacity at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspCVEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

8.      Specific isochoric heat capacity [J/(kgK)] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspCVHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

9.      Specific isochoric heat capacity [J/(kgK)] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspCVPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspCVxPT or wspCVSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

10.  Specific isochoric heat capacity [J/(kgK)] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspCVPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspCVxPT or wspCVSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

11.  Specific isochoric heat capacity [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCVPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspCVxPT) is used.

12.  Specific isochoric heat capacity [J/(kgK)] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspCVPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspCVxPT or wspCVSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

13.  Static dielectric constant of ordinary water substance [-] as function of pressure p [Pa], temperature t [K]:

wspDCPT(p, t)

where:

Function is based upon the function wspDCRT(r, t). The latest in it's turn is based upon the "Release on the Static Dielectric Constant of Ordinary Water Substance for Temperatures from 238K to 873K and Pressures up to 1000 MPa", 1997 from IAPWS. The range of validity is from 238 to 273K in the metastable liquid at atmospheric pressure (0.101325 MPa); from 273 to 323 K at pressures up to the lower of the ice IV melting pressure or 1000 MPa; above 323 K at pressures up to 600 MPa. The formulation also extrapolates smoothly up to at least 1200 K and 1200 MPa.

14.  Density [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspDENSPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspDENSxPT) is used.

15.  Density [kg/m3] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspDENSPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspDENSxPT or wspDENSSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

16.  Dynamic viscosity at the end of expansion/compression process [Pasec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspDYNVISEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

17.  Dynamic viscosity at the end of expansion/compression process [Pasec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspDYNVISEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

18.  Dynamic viscosity [Pasec] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspDYNVISHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

19.  Dynamic viscosity [Pasec] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspDYNVISPH(p, h)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspDYNVISPT or wspDYNVISSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

20.  Dynamic viscosity [Pasec] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspDYNVISPS(p, s)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspDYNVISPT or wspDYNVISSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

21.  Dynamic viscosity [Pasec] as function of pressure p [Pa], temperature t [K]:

wspDYNVISPT(p, t)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspDYNVISRT is used for calculation the argument of which is density defined by reverse value of the function wspVPT.

22.  Dynamic viscosity [Pasec] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspDYNVISPTX(p, t, x)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspDYNVISPT or wspDYNVISSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

23.  Specific enthalpy at the end of expansion/compression process [J/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspHEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

24.  Specific enthalpy at the end of expansion/compression process [J/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspHEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

25.  Specific enthalpy [J/kg] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspHPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspHxPT or wspHSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

26.  Specific enthalpy [J/kg] as function of pressure p [Pa], temperature t [K]:

wspHPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspHxPT) is used.

27.  Specific enthalpy [J/kg] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspHPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspHxPT or wspHSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

28.  Joule-Thomson coefficient at the end of expansion/compression process [K/Pa] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspJOULETHOMPSONEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

29.  Joule-Thomson coefficient at the end of expansion/compression process [K/Pa] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspJOULETHOMPSONEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

30.  Joule-Thomson coefficient [K/Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspJOULETHOMPSONHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

31.  Joule-Thomson coefficient [K/Pa] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspJOULETHOMPSONPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspJOULETHOMPSONPT or wspJOULETHOMPSONSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

32.  Joule-Thomson coefficient [K/Pa] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspJOULETHOMPSONPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspJOULETHOMPSONPT or wspJOULETHOMPSONSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

33.  Joule-Thomson coefficient [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSONPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspJOULETHOMPSONxPT) is used.

34.  Joule-Thomson coefficient [K/Pa] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspJOULETHOMPSONPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspJOULETHOMPSONPT or wspJOULETHOMPSONSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

35.  Isoentropic exponent at the end of expansion/compression process [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspKEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

36.  Isoentropic exponent [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspKEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

37.  Isoentropic exponent [-] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspKHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

38.  Kinematic viscosity at the end of expansion/compression process [m2/sec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspKINVISEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

39.  Kinematic viscosity at the end of expansion/compression process [m2/sec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspKINVISEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

40.  Kinematic viscosity [m2/sec] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspKINVISHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

41.  Kinematic viscosity [m2/sec] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspKINVISPH(p, h)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspKINVISPT or wspKINVISSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

42.  Kinematic viscosity [m2/sec] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspKINVISPS(p, s)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspKINVISPT or wspKINVISSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

43.  Kinematic viscosity [m2/sec] as function of pressure p [Pa], temperature t [K]:

wspKINVISPT(p, t)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The equation is: KINVIS = DYNVIS V is used for calculation, where DYNVIS defined by the function wspDYNVISPT and V - by wspVPT.

44.  Kinematic viscosity [m2/sec] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspKINVISPTX(p, t, x)

where:

The range of the validity is within that described in IF-97 and in formulation for dynamic viscosity (see function wspDYNVISRT). The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspKINVISPT or wspKINVISSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

45.  Isoentropic exponent [-] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspKPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspKPT or wspKSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

46.  Isoentropic exponent [-] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspKPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspKPT or wspKSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

47.  Isoentropic exponent [-] as function of pressure p [Pa], temperature t [K]:

wspKPT(p, t)

where:

The range of the validity is within that described in IF-97. Used equation is: K = W W / (P V) is used for calculation, where W (sound velocity) is defined by the function wspWPT, P - pressure and V (specific volume) - by wspVPT.

48.  Isoentropic exponent [-] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspKPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspKPT or wspKSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

49.  Ion product of water substance [mole2/kg2] as function of pressure p [Pa], temperature t [K]:

wspKWPT(p, t)

where:

Function is based upon the function wspKWRT(r, t). The latest in it's turn based upon the "Release on the Ion Product of Water Substance, May 1980" from IAPWS. The range of validity is for pressure from 1 to 10000 bar, temperature is from 0 to 1000C. Also is recommended do not use this formulation for densities less than 0.45 g/cm3. Warning! Function return ion product constant in SI base units (mole/kg)^2 (molality^2) but usually used the dimension (mole/dm^3)^2 (molarity^2). To convert from (mole/kg)^2 to (mole/dm^3)^2 you must to multiply value by squared density (r * r). Density must be in kg/dm^3.

50.  Pressure [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

51.  Prandtl number at the end of expansion/compression process [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspPRANDTLEEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

52.  Prandtl number at the end of expansion/compression process [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspPRANDTLEEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

53.  Prandtl number [-] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPRANDTLEHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

54.  Prandtl number [-] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspPRANDTLEPH(p, h)

where:

The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspPRANDTLEPT or wspPRANDTLESTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

55.  Prandtl number [-] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspPRANDTLEPS(p, s)

where:

The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspPRANDTLEPT or wspPRANDLTESTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

56.  Prandtl number [-] as function of pressure p [Pa], temperature t [K]:

wspPRANDTLEPT(p, t)

where:

The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity. The equation is: Pr = DYNVIS CP / THERMCOND is used for calculation, where DYNVIS is calculated by the function wspDYNVISPT, CP - by wspCPPT and THERMCOND - by wspTHERMCONDPT.

57.  Prandtl number [-] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspPRANDTLEPTX(p, t, x)

where:

The range of the validity is within that described in IF-97 and in formulations for the properties included in the equation to calculate sought quantity (see function wspPRANDTLEPT, wspTHERMCONDRT, wspDYNVISRT). The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspPRANDTLEPPT or wspPRANDTLEPSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

58.  Properties calculation result as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPTHS(h, s, *p, *t)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantities are calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

59.  Specific entropy at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspSEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

60.  Specific entropy at the end of expansion/compression process [J/(kgK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspSEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

61.  Specific entropy [J/(kgK)] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspSPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspSPTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

62.  Specific entropy [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspSPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspSxPT) is used.

63.  Specific entropy [J/(kgK)] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspSPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspSxPT or wspSSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

64.  Surface tension [N/m] as function of temperature t [K]:

wspSURFTENT(t)

where:

The function is based on the IAPWS Release on The Surface Tension of Ordinary Water Substance 1995. The range of validity is from triple point (0.01C) to 647.096 K.

65.  Temperature at the end of expansion/compression process [K] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspTEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

66.  Temperature at the end of expansion/compression process [K] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspTEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

67.  Thermal conductivity coefficient at the end of expansion/compression process [W/(mK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspTHERMCONDEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

68.  Thermal conductivity coefficient at the end of expansion/compression process [W/(mK)] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspTHERMCONDEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

69.  Thermal conductivity coefficient [W/(mK)] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspTHERMCONDHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

70.  Thermal conductivity coefficient [W/(mK)] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspTHERMCONDPH(p, h)

where:

The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspTHERMCONDPT or wspTHERMCONDSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

71.  Thermal conductivity coefficient [W/(mK)] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspTHERMCONDPS(p, s)

where:

The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspTHERMCONDPT or wspTHERMCONDSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

72.  Thermal conductivity coefficient [W/(mK)] as function of pressure p [Pa], temperature t [K]:

wspTHERMCONDPT(p, t)

where:

The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). The function wspTHERMCONDRT is used for calculation the argument of which is density defined by the reverse value of function wspVPT.

73.  Thermal conductivity coefficient [W/(mK)] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspTHERMCONDPTX(p, t, x)

where:

The range of the validity is within that described in IF-97 and in formulation for thermal conductivity (see function wspTHERMCONDRT). The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspTHERMCONDPT or wspTHERMCONDSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

74.  Temperature [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspTHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

75.  Temperature [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspTPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

76.  Temperature [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspTPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

77.  Specific internal energy at the end of expansion/compression process [J/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspUEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

78.  Specific internal energy at the end of expansion/compression process [J/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspUEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), w- specific enthalpy in initial point, h1 - specific enthalpy in final point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in final point with isoentropic expansion (s = Const). For the calculating the process of compression you must to use instead the efficiency eff the inverse number (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

79.  Specific internal energy [J/kg] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspUHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

80.  Specific internal energy [J/kg] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspUPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspUxPT or wspUSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

81.  Specific internal energy [J/kg] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspUPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspUxPT or wspUSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

82.  Specific internal energy [J/kg] as function of pressure p [Pa], temperature t [K]:

wspUPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspUxPT) is used.

83.  Specific internal energy [J/kg] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspUPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspUxPT or wspUSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

84.  Specific volume at the end of expansion/compression process [m3/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspVEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

85.  Specific volume at the end of expansion/compression process [m3/kg] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspVEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

86.  Specific volume [m3/kg] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspVHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

87.  Specific volume [m3/kg] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspVPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspVPTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

88.  Specific volume [m3/kg] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspVPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspVxPT or wspVPTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

89.  Specific volume [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspVPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspVxPT) is used.

90.  Specific volume [m3/kg] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspVPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspVxPT or wspVSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

91.  Properties calculation result as function of pressure p [Pa], temperature t [K]:

wspVUSHCVWDERPTPT(p, t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspVUSHDERPTxPT) is used to return the set of properties and so speed up calculation of several values in one point.

Note: In Mathcad the function have only parameters (p, t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (p, t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

92.  Sound velocity at the end of expansion/compression process [m/sec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspWEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

93.  Sound velocity at the end of expansion/compression process [m/sec] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspWEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

94.  Speed of sound [m/sec] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspWHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

95.  Speed of sound [m/sec] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspWPH(p, h)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPH is used to determine the IF-97 region. Then the necessary function (wspTxPH) is used. Finally, the necessary function (wspWxPT or wspWSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

96.  Speed of sound [m/sec] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspWPS(p, s)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREAPS is used to determine the IF-97 region. Then the necessary function (wspTxPS) is used. Finally, the necessary function (wspWxPT or wspWSTX) is called. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

97.  Speed of sound [m/sec] as function of pressure p [Pa], temperature t [K]:

wspWPT(p, t)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA2 is used to determine the IF-97 region. Then the necessary function (wspWxPT) is used.

98.  Speed of sound [m/sec] as function of pressure p [Pa], temperature t [K], vapor fraction x [-]:

wspWPTX(p, t, x)

where:

The range of the validity is within that described in IF-97. The function wspWATERSTATEAREA is used to determine the IF-97 region. Then the necessary function (wspWxPT or wspWSTX) is used. If the point is out of the double-phase area the value of vapor fraction is ignored.

99.  Vapor fraction at the end of expansion/compression process [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspXEXPANSIONPTPEFF(p0, t0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

100.                     Vapor fraction at the end of expansion/compression process [-] as function of pressure at initial point p0 [Pa], temperature at initial point t0 [K], vapor fraction at initial point x0 [-], pressure at final point p1 [Pa], internal efficiency of process eff [-]:

wspXEXPANSIONPTXPEFF(p0, t0, x0, p1, eff)

where:

It is a function. The range of the validity is within that described in IF-97. Function returns a parameter at the end of the expansion process from initial point with arguments p0, t0, x0 to pressure p1 with relative efficiency eff. The process is calculated by the formula: h1 = h0 - eff * (h0 - h1isoent), where h0 - specific enthalpy in initial point, h1 - specific enthalpy in finite point of the expansion, eff - efficiency, h1isoentr - specific enthalpy in finite point at isoentropic equilibrium expansion (s = Const). For compression process calculation it is necessary to use the value inverse to the efficiency eff (i.e. 1/eff). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

101.                     Vapor fraction [-] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspXHS(h, s)

where:

It is a function. The range of the validity is within that described in IF-97. Function works as follows. For first the IF-97 region is determined by function wspWATERSTATEAREAHS. Then the original variables are defined for basic equation of IF-97 region at this point (by functions wspPTxHS and wspRT3HS). Then the sought quantity is calculated on these data. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

102.                     Vapor fraction [-] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspXPH(p, h)

where:

It is a function. The range of the validity is within that described in IF-97. The IF-97 region is determined for given p and h (by function wspWATERSTATEAREAPH). If the region is double-phase the sought parameter is calculated. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

103.                     Vapor fraction [-] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspXPS(p, s)

where:

It is a function. The range of the validity is within that described in IF-97. The IF-97 region is determined for given p and s (by function wspWATERSTATEAREAPS). If the region is double-phase the sought parameter is calculated. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Category: Saturation Line

Functions for calculation properties at saturation line for water and vapor.

104.                     Specific isobaric heat capacity of steam at saturation line [J/(kgK)] as function of temperature t [K]:

wspCPSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

105.                     Specific isobaric heat capacity of water at saturation line [J/(kgK)] as function of temperature t [K]:

wspCPSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

106.                     Specific isochoric heat capacity of steam at saturation line from the double-phase region [J/(kgK)] as function of temperature t [K]:

wspCVDPSST(t)

where:

Function returns the sum of specific isochoric heat capacity of steam from one-phase region (function wspCVSST) and additional increment calculated by corresponding thermodynamic formulas. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

107.                     Specific isochoric heat capacity of water at saturation line from the double-phase region [J/(kgK)] as function of temperature t [K]:

wspCVDPSWT(t)

where:

Function returns the sum of specific isochoric heat capacity of water from one-phase region (function wspCVSWT) and additional increment calculated by corresponding thermodynamic formulas. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

108.                     Specific isochoric heat capacity of steam at saturation line from the one-phase region [J/(kgK)] as function of temperature t [K]:

wspCVSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

109.                     Specific isochoric heat capacity of water at saturation line from the one-phase region [J/(kgK)] as function of temperature t [K]:

wspCVSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

110.                     Density of steam at saturation line [kg/m3] as function of temperature t [K]:

wspDENSSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

111.                     Density of water at saturation line [kg/m3] as function of temperature t [K]:

wspDENSSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

112.                     Derivative of saturation pressure on saturation temperature [Pa/K] as function of temperature t [K]:

wspDPDTST(t)

where:

Function is based upon "Supplementary Release on Saturation Properties of Ordinary Water Substance" from International Association for the Properties of Water and Steam. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

113.                     Dynamic viscosity of steam at saturation line [Pasec] as function of temperature t [K]:

wspDYNVISSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

114.                     Dynamic viscosity of water at saturation line [Pasec] as function of temperature t [K]:

wspDYNVISSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

115.                     Specific enthalpy of steam at saturation line [J/kg] as function of temperature t [K]:

wspHSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

116.                     Specific enthalpy of water at saturation line [J/kg] as function of temperature t [K]:

wspHSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

117.                     Joule-Thomson coefficient of steam at saturation line [K/Pa] as function of temperature t [K]:

wspJOULETHOMPSONSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

118.                     Joule-Thomson coefficient of water at saturation line [K/Pa] as function of temperature t [K]:

wspJOULETHOMPSONSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

119.                     Kinematic viscosity of steam at saturation line [m2/sec] as function of temperature t [K]:

wspKINVISSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

120.                     Kinematic viscosity of water at saturation line [m2/sec] as function of temperature t [K]:

wspKINVISSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

121.                     Isoentropic exponent of steam at saturation line [-] as function of temperature t [K]:

wspKSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

122.                     Isoentropic exponent of water at saturation line [-] as function of temperature t [K]:

wspKSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

123.                     Prandtl number of steam at saturation line [-] as function of temperature t [K]:

wspPRANDTLESST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

124.                     Prandtl number of water at saturation line [-] as function of temperature t [K]:

wspPRANDTLESWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

125.                     Pressure at saturation line [Pa] as function of temperature t [K]:

wspPST(t)

where:

Function is based upon the Supplementary Release on Saturation Properties of Ordinary Water Substance from IAPWS. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

126.                     Specific enthalpy of steam at saturation line [J/kg] as function of specific entropy s [J/(kgK)]:

wspROUGHHSSS(s)

where:

Function uses the special equation from Supplementory Release for IF-97 Formulation. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

127.                     Specific enthalpy of water at saturation line [J/kg] as function of specific entropy s [J/(kgK)]:

wspROUGHHSWS(s)

where:

Function uses the special equation from Supplementory Release for IF-97 Formulation. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

128.                     Rough value of density of steam at saturation line [kg/m3] as function of temperature t [K]:

wspROUGHRSST(t)

where:

This function allows to calculate quickly rough value of steam density at saturation line. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

129.                     Rough value of density of water at saturation line [kg/m3] as function of temperature t [K]:

wspROUGHRSWT(t)

where:

This function allows to calculate quickly rough value of steam density at saturation line. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

130.                     Specific evaporation heat [J/kg] as function of temperature t [K]:

wspRST(t)

where:

The result of function is calculated by formula r = (hs - hw), where hs - specific enthalpy of steam at saturation line, hw - specific enthalpy of water at saturation line. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

131.                     Specific entropy of steam at saturation line [J/(kgK)] as function of temperature t [K]:

wspSSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

132.                     Specific entropy of water at saturation line [J/(kgK)] as function of temperature t [K]:

wspSSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

133.                     Thermal conductivity coefficient of steam at saturation line [W/(mK)] as function of temperature t [K]:

wspTHERMCONDSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

134.                     Thermal conductivity coefficient of water at saturation line [W/(mK)] as function of temperature t [K]:

wspTHERMCONDSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

135.                     Temperature at saturation line [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspTSHS(h, s)

where:

Function calculates saturation temperature on a base of thermodunamic formulas. Newton method is used to determine the root. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

136.                     Temperature at saturation line [K] as function of pressure p [Pa]:

wspTSP(p)

where:

Function is based upon the Supplementary Release on Saturation Properties of Ordinary Water Substance from IAPWS. The range of validity is from triple point of water (0.01C or 611.657 Pa) to critical point (647.096K or 373.946C or 22.064 MPa).

137.                     Specific internal energy of steam at saturation line [J/kg] as function of temperature t [K]:

wspUSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

138.                     Specific internal energy of water at saturation line [J/kg] as function of temperature t [K]:

wspUSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

139.                     Specific volume of steam at saturation line [m3/kg] as function of temperature t [K]:

wspVSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

140.                     Specific volume of water at saturation line [m3/kg] as function of temperature t [K]:

wspVSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

141.                     Properties calculation result for steam at saturation line as function of temperature t [K]:

wspVUSHCVWDERPTSST(t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

The function is based upon IF-97 formulation. It returns the property set that speed up calculation considerably. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

Note: In Mathcad the function have only parameters (t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

142.                     Properties calculation result for water at saturation line as function of temperature t [K]:

wspVUSHCVWDERPTSWT(t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

The function is based upon IF-97 formulation. It returns the property set that speed up calculation considerably. The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

Note: In Mathcad the function have only parameters (t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

143.                     Sound velocity in steam at saturation line [m/sec] as function of temperature t [K]:

wspWSST(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

144.                     Sound velocity in water at saturation line [m/sec] as function of temperature t [K]:

wspWSWT(t)

where:

The range of validity is from triple point of water (0.01C) to critical point (647.096K or 373.946C).

Category: Double-phase area

Functions for calculation properties in double-phase area: water-vapor.

145.                     Specific isobaric heat capacity in double-phase area [J/(kgK)] as function of temperature t [K], vapor fraction x [-]:

wspCPSTX(t, x)

where:

Function uses functions wspCPSST and wspCPSWT which return specific heat capacities at constant pressure (Cp) of steam and water at saturation line. The function uses following formula: CPx = (1 - X)CPw + XCPs, where CPs = wspCPSST, CPw = wspCPSWT and X - vapor fraction. The range of calculation is double-phase (water-vapor) area.

146.                     Specific isochoric heat capacity in double-phase area [J/(kgK)] as function of temperature t [K], vapor fraction x [-]:

wspCVSTX(t, x)

where:

Function uses functions wspCVDPSST and wspCVDPSWT which return specific heat capacities at constant volume (Cv) of steam and water at saturation line in the double-phase region area. The function uses following formula: CVx = (1 - X)CVDPw + XCVDPs, where CVDPs = wspCVDPSST, CVDPw = wspCVDPSWT and X - vapor fraction. This function used functions wspCVSST and wspCVSWT until version 5.6 which is not conform to thermodynamic, because it neglects the jump of properties at saturation line. The range of validity is double-phase (water-vapor) area.

147.                     Density in double-phase area [kg/m3] as function of temperature t [K], vapor fraction x [-]:

wspDENSSTX(t, x)

where:

Function uses the function named wspVSTX which return specific volume in double-phase area. The function uses following formula: DENSx = 1 / Vx, where Vx = wspVSTX(t, x). The range of validity is double-phase (water-vapor) area.

148.                     Dynamic viscosity in double-phase area [Pasec] as function of temperature t [K], vapor fraction x [-]:

wspDYNVISSTX(t, x)

where:

Function uses functions wspDYNVISSST and wspDYNVISSWT which return dynamic viscosity of steam and water at saturation line. The function uses following formula: DYNVISx = (1 - X)DYNVISw + XDYNVISs, where DYNVISs = wspDYNVISSST, DYNVISw = wspDYNVISSWT and X - vapor fraction. The range of calculation is double-phase (water-vapor) area.

149.                     Specific enthalpy in double-phase area [J/kg] as function of temperature t [K], vapor fraction x [-]:

wspHSTX(t, x)

where:

Function uses functions wspHSST and wspHSWT which return specific enthalpies of steam and water at saturation line. The function uses following formula: Hx = (1 - X)Hw + XHs, where Hs = wspHSST, Hw = wspHSWT and X - vapor fraction. The range of validity is double-phase (water-vapor) area.

150.                     Joule-Thomson coefficient in double-phase area [K/Pa] as function of temperature t [K], vapor fraction x [-]:

wspJOULETHOMPSONSTX(t, x)

where:

Function calculate the Joule-Thomson coefficient in double-phase area by equation at saturation line: JT = (dT/dP)h. Until version 5.6 the function used wspJOULETHOMPSONSST and wspJOULETHOMPSONSWT and returned value from following formula: JTx = (1 - X)JOULETHOMPSONw + XJOULETHOMPSONs, where JOULETHOMPSONs = wspJOULETHOMPSONSST, JOULETHOMPSONw = wspJOULETHOMPSONSWT and X - vapor vapor fraction that did not conform to definition of Joule-Thomson coefficient. Now it calculates correctly. The range of validity is double-phase (water-vapor) area.

151.                     Kinematic viscosity in double-phase area [m2/sec] as function of temperature t [K], vapor fraction x [-]:

wspKINVISSTX(t, x)

where:

Function uses functions wspKINVISSST and wspKINVISSWT which return kinematic viscosities of steam and water at saturation line. The function uses next formula: KINVISx = (1 - X)KINVISw + XKINVISs, where KINVISs = wspKINVISSST, KINVISw = wspKINVISSWT and X - vapor fraction. The range of calculation is double-phase (water-vapor) area.

152.                     Isoentropic exponent in double-phase area [-] as function of temperature t [K], vapor fraction x [-]:

wspKSTX(t, x)

where:

Function calculate isoentropic exponent in double-phase area considering property jump at saturation line. Until version 5.6 this function uses following formula: Kx = (1 - X)Kw + XKs, where Ks = wspKSST, Kw = wspKSWT and X - vapor fraction. But that is not agree with thermodynamic. Now it is calculated properly. The range of validity is double-phase (water-vapor) area.

153.                     Prandtl number in double-phase area [-] as function of temperature t [K], vapor fraction x [-]:

wspPRANDTLESTX(t, x)

where:

Function uses functions wspPRANDTLESST and wspPRANDTLESWT which returns Prandtl numbers of steam and water at saturation line. The function uses following formula: PRANDTLEx = (1 - X)PRANDTLEw + XPRANDTLEs, where PRANDTLEs = wspPRANDTLESST, PRANDTLEw = wspPRANDTLESWT and X - vapor fraction. The range of calculation is double-phase (water-vapor) area.

154.                     Specific entropy in double-phase area [J/(kgK)] as function of temperature t [K], vapor fraction x [-]:

wspSSTX(t, x)

where:

Function uses functions wspSSST and wspSSWT which return specific entropies of steam and water at saturation line. The function uses following formula: Sx = (1 - X)Sw + XSs, where Ss = wspSSST, Sw = wspSSWT and X - vapor fraction. The range of validity is double-phase (water-vapor) area.

155.                     Thermal conductivity coefficient in double-phase area [W/(mK)] as function of temperature t [K], vapor fraction x [-]:

wspTHERMCONDSTX(t, x)

where:

Function uses functions wspTHERMCONDSST and wspTHERMCONDSWT which return thermal conductivity of steam and water at saturation line. The function uses following formula: THERMCONDx = (1 - X)THERMCONDw + XTHERMCONDs, where THERMCONDs = wspTHERMCONDSST, THERMCONDw = wspTHERMCONDSWT and X - vapor fraction. The range of calculation is double-phase (water-vapor) area.

156.                     Properties calculation result in double-phase area as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspTXSHS(h, s, *t, *x)

where:

Function calculates saturation temperature and vapor fraction on the base of thermodunamic formulas. Newton method is used to determine the root. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE). The range of validity is double-phase (water-vapor) area.

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

157.                     Specific internal energy in double-phase area [J/kg] as function of temperature t [K], vapor fraction x [-]:

wspUSTX(t, x)

where:

Function uses functions wspUSST and wspUSWT which return specific internal energies of steam and water at saturation line. The function uses following formula: Ux = (1 - X)Uw + XUs, where Us = wspUSST, Uw = wspUSWT and X - vapor fraction. The range of validity is double-phase (water-vapor) area.

158.                     Specific volume in double-phase area [m3/kg] as function of temperature t [K], vapor fraction x [-]:

wspVSTX(t, x)

where:

Function uses the functions wspVSST and wspVSWT which return specific volumes of steam and water at saturation line. The function uses following formula: Vx = (1 - X)Vw + XVs, where Vs = wspVSST, Vw = wspVSWT and X - vapor fraction. The range of validity is double-phase (water-vapor) area.

159.                     Sound velocity in double-phase area [m/sec] as function of temperature t [K], vapor fraction x [-]:

wspWSTX(t, x)

where:

Function uses function wspKSTX to calculate sound velocity in double-phase area conforming with thermodynamic. The function neglects the flow structure. Until version 5.6 this function used following formula: Wx = (1 - X)Ww + XWs, where Ws = wspWSST, Ww = wspWSWT and X - vapor fraction, which disagrees with thermodynamic because it neglects the property jump at saturation line. The range of validity is double-phase (water-vapor) area.

160.                     Vapor fraction [-] as function of temperature t [K], specific isobaric heat capacity Cp [J/(kgK)]:

wspXSTCP(t, Cp)

where:

This function uses formula: X = (CP - CPw)/(CPs - CPw), where CPw = wspCPSWT, CPs = wspCPSST. The range of calculation is double-phase (water-vapor) area.

161.                     Vapor fraction [-] as function of temperature t [K], specific isochoric heat capacity Cv [J/(kgK)]:

wspXSTCV(t, Cv)

where:

This function uses formula: X = (CV - CVDPw)/(CVDPs - CVDPw), where CVDPw = wspCVDPSWT, CVDPs = wspCVDPSST. The range of validity is double-phase (water-vapor) area.

162.                     Vapor fraction [-] as function of temperature t [K], density r [kg/m3]:

wspXSTDENS(t, r)

where:

This function uses formula: X = rs * (rw - r) / (r * (rw - rs)), where rw = wspDENSSWT, rs = wspDENSSST. The range of validity is double-phase (water-vapor) area.

163.                     Vapor fraction [-] as function of temperature t [K], dynamic viscosity dv [Pasec]:

wspXSTDYNVIS(t, dv)

where:

This function uses formula: X = (DV - DVw)/(DVs - DVw), where DVw = wspDYNVISSWT, DVs = wspDYNVISSST. The range of calculation is double-phase (water-vapor) area.

164.                     Vapor fraction [-] as function of temperature t [K], specific enthalpy h [J/kg]:

wspXSTH(t, h)

where:

This function uses formula: X = (H - Hw)/(Hs - Hw), where Hw = wspHSWT, Hs = wspHSST. The range of validity is double-phase (water-vapor) area.

165.                     Vapor fraction [-] as function of temperature t [K], Joule-Thomson coefficient jt [K/Pa]:

wspXSTJOULETHOMPSON(t, jt)

where:

This function uses formula: X = (JOULETHOMPSON - JOULETHOMPSONw)/(JOULETHOMPSONs - JOULETHOMPSONw), where JOULETHOMPSONw = wspJOULETHOMPSONSWT, JOULETHOMPSONs = wspJOULETHOMPSONSST. The range of validity is double-phase (water-vapor) area.

166.                     Vapor fraction [-] as function of temperature t [K], isoentropic exponent k [-]:

wspXSTK(t, k)

where:

This function uses formula: X = (K - Kw)/(Ks - Kw), where Kw = wspKSWT, Ks = wspKSST. The range of validity is double-phase (water-vapor) area.

167.                     Vapor fraction [-] as function of temperature t [K], kinematic viscosity kv [m2/sec]:

wspXSTKINVIS(t, kv)

where:

This function uses formula: X = (KV - KVw)/(KVs - KVw), where KVw = wspKINVISSWT, KVs = wspKINVISSST. The range of calculation is double-phase (water-vapor) area.

168.                     Vapor fraction [-] as function of temperature t [K], Prandtl number pr [-]:

wspXSTPRANDTLE(t, pr)

where:

This function uses formula: X = (PR - PRw)/(PRs - PRw), where PRw = wspPRANDTLESWT, PRs = wspPRANDTLESST. The range of calculation is double-phase (water-vapor) area.

169.                     Vapor fraction [-] as function of temperature t [K], specific entropy s [J/(kgK)]:

wspXSTS(t, s)

where:

This function uses formula: X = (S - Sw)/(Ss - Sw), where Sw = wspSSWT, Ss = wspSSST. The range of validity is double-phase (water-vapor) area.

170.                     Vapor fraction [-] as function of temperature t [K], thermal conductivity coefficient tc [W/(mK)]:

wspXSTTHERMCOND(t, tc)

where:

This function uses formula: X = (TC - TCw)/(TCs - TCw), where TCw = wspTHERMCONDSWT, TCs = wspTHERMCONDSST. The range of calculation is double-phase (water-vapor) area.

171.                     Vapor fraction [-] as function of temperature t [K], specific internal energy u [J/kg]:

wspXSTU(t, u)

where:

This function uses formula: X = (U - Uw)/(Us - Uw), where Uw = wspUSWT, Us = wspUSST. The range of validity is double-phase (water-vapor) area.

172.                     Vapor fraction [-] as function of temperature t [K], specific volume v [m3/kg]:

wspXSTV(t, v)

where:

This function uses formula: X = (V - Vw)/(Vs - Vw), where Vw = wspVSWT, Vs = wspVSST. The range of validity is double-phase (water-vapor) area.

173.                     Vapor fraction [-] as function of temperature t [K], speed of sound w [m/sec]:

wspXSTW(t, w)

where:

From version 5.6 this function calculates value of sound velocity in double-phase region from isoentropic exponent that is conform to thermodynamic. Until version 5.6 this function used formula: X = (W - Ww)/(Ws - Ww), where Ww = wspWSWT, Ws = wspWSST which was not right. Now it is calculated correctly. The range of validity is double-phase (water-vapor) area.

Category: Sublimation and Melting Lines

Functions for calculation properties at sublimation and melting lines.

174.                     Pressure at melting line of ice I [Pa] as function of temperature t [K]:

wspPMELTIT(t)

where:

Function is based on "Release on the Pressure along the Melting and the Sublimation Curves of Ordinary Water Substance" IAPWS (Milan, Italy, September 1993).

175.                     Pressure at sublimation line [Pa] as function of temperature t [K]:

wspPSUBT(t)

where:

Function is based on "Release on the Pressure along the Melting and the Sublimation Curves of Ordinary Water Substance" IAPWS (Milan, Italy, September 1993).

Category: MetaStable

Functions for calculation properties of meta-stable supercooled steam

176.                     Specific isobaric heat capacity of meta-stable supercooled steam [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCPMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

177.                     Specific isochoric heat capacity of meta-stable supercooled steam [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCVMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

178.                     Dynamic viscosity of meta-stable supercooled steam [Pasec] as function of pressure p [Pa], temperature t [K]:

wspDYNVISMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

179.                     Specific enthalpy of meta-stable supercooled steam [J/kg] as function of pressure p [Pa], temperature t [K]:

wspHMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

180.                     Joule-Thomson coefficient of meta-stable supercooled steam [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSONMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

181.                     Kinematic viscosity of meta-stable supercooled steam [m2/sec] as function of pressure p [Pa], temperature t [K]:

wspKINVISMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

182.                     Isoentropic exponent of meta-stable supercooled steam [-] as function of pressure p [Pa], temperature t [K]:

wspKMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

183.                     Prandtl number of meta-stable supercooled steam [-] as function of pressure p [Pa], temperature t [K]:

wspPRANDTLEMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

184.                     Specific entropy of meta-stable supercooled steam [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspSMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

185.                     Thermal conductivity coefficient of meta-stable supercooled steam [W/(mK)] as function of pressure p [Pa], temperature t [K]:

wspTHERMCONDMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

186.                     Specific internal energy of meta-stable supercooled steam [J/kg] as function of pressure p [Pa], temperature t [K]:

wspUMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

187.                     Specific volume of meta-stable supercooled steam [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspVMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

188.                     Properties calculation result of super-cooled steam (meta-stable region) as function of pressure p [Pa], temperature t [K]:

wspVUSHCVWDERPTMSPT(p, t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

Note: In Mathcad the function have only parameters (p, t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (p, t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

189.                     Sound velocity of meta-stable supercooled steam [m/sec] as function of pressure p [Pa], temperature t [K]:

wspWMSPT(p, t)

where:

This function calculates the properties of super-cooled steam (in meta-stable region). The parameters range: pressure up to 10 MPa, vapor fraction above 95%. Special equation for meta-stable vapor region is used for calculation which is based upon IAPWS IF-97.

Category: Source

Functions from IAPWS IF-97 and other formulations.

190.                     Specific isobaric heat capacity in IF-97 region 1 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCP1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

191.                     Specific isobaric heat capacity in IF-97 region 2 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCP2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

192.                     Specific isobaric heat capacity in IF-97 region 3 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCP3PT(p, t)

where:

Function is based upon functions wspR3PT and wspU3RT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

193.                     Specific isobaric heat capacity in IF-97 region 3 [J/(kgK)] as function of density r [kg/m3], temperature t [K]:

wspCP3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

194.                     Specific isobaric heat capacity in IF-97 region 5 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCP5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

195.                     Specific isochoric heat capacity in IF-97 region 1 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCV1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

196.                     Specific isochoric heat capacity in IF-97 region 2 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCV2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

197.                     Specific isochoric heat capacity in IF-97 region 3 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCV3PT(p, t)

where:

Function is based upon functions wspR3PT and wspU3RT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

198.                     Specific isochoric heat capacity in IF-97 region 3 [J/(kgK)] as function of density r [kg/m3], temperature t [K]:

wspCV3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

199.                     Specific isochoric heat capacity in IF-97 region 5 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspCV5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

200.                     Static dielectric constant of ordinary water substance [-] as function of density r [kg/m3], temperature t [K]:

wspDCRT(r, t)

where:

Function is based upon the "Release on the Static Dielectric Constant of Ordinary Water Substance for Temperatures from 238K to 873K and Pressures up to 1000 MPa", 1997 from IAPWS. The range of validity is from 238 to 273K in the metastable liquid at atmospheric pressure (0.101325 MPa); from 273 to 323 K at pressures up to the lower of the ice IV melting pressure or 1000 MPa; above 323 K at pressures up to 600 MPa. The formulation also extrapolates smoothly up to at least 1200 K and 1200 MPa.

201.                     Density in IF-97 region 1 [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspDENS1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

202.                     Density in IF-97 region 2 [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspDENS2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

203.                     Density in IF-97 region 3 [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspDENS3PT(p, t)

where:

Function calculates the density in region 3 by calling the function wspR3PTR0 with the corresponding initial estimate for density. Used for unification of the functions in all IF-97 regions. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

204.                     Density in IF-97 region 5 [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspDENS5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

205.                     Dynamic viscosity [Pasec] as function of density r [kg/m3], temperature t [K]:

wspDYNVISRT(r, t)

where:

Function is based upon the Revised Release on the IAPS Formulation 1985 for the Viscosity of Ordinary Water Substance with correction for ITS-90 (International Temperature Scale). The range of validity is for pressure up to 500 MPa temperature is from 0 to 150C, for pressure up to 350 MPa temperature is from 150 to 600C, for pressure up to 300 MPa temperature is from 600 to 900C.

206.                     Specific enthalpy in IF-97 region 1 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspH1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

207.                     Specific enthalpy at line between areas 2b and 2c [J/kg] as function of pressure p [Pa]:

wspH2B2CP(p)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

208.                     Specific enthalpy in IF-97 region 2 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspH2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

209.                     Specific enthalpy in IF-97 region 3 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspH3PT(p, t)

where:

Function is based upon functions wspR3PT and wspU3RT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

210.                     Specific enthalpy in IF-97 region 3 [J/kg] as function of density r [kg/m3], temperature t [K]:

wspH3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

211.                     Specific enthalpy in IF-97 region 5 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspH5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

212.                     Specific enthalpy at boundary line between areas 1 and 3 [J/kg] as function of specific entropy s [J/(kgK)]:

wspHB13S(s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004).

213.                     Joule-Thomson coefficient in IF-97 region 1 [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSON1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The equation: JT = (T*dV/dT - V)/Cp is used, where T - temperature, dV/dT - differential quotient (dV/dT)p, V - volume, Cp - isobaric heat capacity (Cp).

214.                     Joule-Thomson coefficient in IF-97 region 2 [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSON2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The equation: JT = (T*dV/dT - V)/Cp is used, where T - temperature, dV/dT - differential quotient (dV/dT)p, V - volume, Cp - isobaric heat capacity (Cp).

215.                     Joule-Thomson coefficient in IF-97 region 3 [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSON3PT(p, t)

where:

Function calculates Joule-Thomson coefficient for region 3 of IF-97 Formulation. This function uses function wspR3PT(p, t) to calculate density and then returns the value from the function wspJOULETHOMPSON3RT(r, t). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

216.                     Joule-Thomson coefficient in IF-97 region 3 [K/Pa] as function of density r [kg/m3], temperature t [K]:

wspJOULETHOMPSON3RT(r, t)

where:

Function calculates the Joule-Thomson coefficient for region 3 of IF-97 Formulation.

217.                     Joule-Thomson coefficient in IF-97 region 5 [K/Pa] as function of pressure p [Pa], temperature t [K]:

wspJOULETHOMPSON5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The equation: JT = (T*dV/dT - V)/Cp is used, where T - temperature, dV/dT - differential quotient (dV/dT)p, V - volume, Cp - isobaric heat capacity (Cp).

218.                     Ion product of water substance [mole2/kg2] as function of density r [kg/m3], temperature t [K]:

wspKWRT(r, t)

where:

Function is based upon the "Release on the Ion Product of Water Substance, May 1980" from IAPWS. The range of validity is for pressure from 1 to 10000 bar, temperature is from 0 to 1000C. Also is recommended do not use this formulation for densities less than 0.45 g/cm3. Warning! Function return ion product constant in SI base units (mole/kg)^2 (molality^2) but usually used the dimension (mole/dm^3)^2 (molarity^2). To convert from (mole/kg)^2 to (mole/dm^3)^2 you must to multiply value by squared density (r * r). Density must be in kg/dm^3.

219.                     Pressure in IF-97 region 1 [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspP1HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for Pressure as a Function of Enthalpy and Entropy p(h,s) to the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2001). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

220.                     Pressure at line between areas 2 and 3 [Pa] as function of temperature t [K]:

wspP23T(t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. Used in function wspWATERSTATEAREA when region of IF-97 is determined.

221.                     Pressure at line between areas 2b and 2c [Pa] as function of specific enthalpy h [J/kg]:

wspP2B2CH(h)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

222.                     Pressure in IF-97 region 2 [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspP2HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for Pressure as a Function of Enthalpy and Entropy p(h,s) to the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2001) and used the functions: wspP2AHS(h, s), wspP2BHS(h, s) and wspP2CHS(h, s). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

223.                     Pressure in IF-97 region 3 [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspP3HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

224.                     Pressure in IF-97 region 3 [Pa] as function of density r [kg/m3], temperature t [K]:

wspP3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

225.                     Pressure in IF-97 region 5 [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspP5HS(h, s)

where:

Function is used to find the pressure in region 5 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

226.                     Pressure at boundary line between areas 2 and 3 [Pa] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPB23HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004).

227.                     Area of phase state as function of pressure p [Pa], temperature t [K]:

wspPHASESTATEPT(p, t)

where:

Function returns code of phase state in point with given p, t. If pressure p or temperature t are above critical (Prk or Tkr), the result is "3" (phase state - supercritical). If the point is situated above the saturation line (in P-T diagramm) the result is "1" (phase state - liquid). If the point is situated below the saturation line (in P-T diagram) the result is "2" (phase state - steam). If point is situated below triple point, the return is error value ("-1").

228.                     Properties calculation result in IF-97 region 1 as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPT1HS(h, s, *p, *t)

where:

Function is used to find parameters in region 1 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

229.                     Properties calculation result in IF-97 region 1 as function of density r [kg/m3], specific enthalpy h [J/kg]:

wspPT1RH(r, h, *p, *t)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (r, h) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (r, h) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

230.                     Properties calculation result in IF-97 region 2 as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPT2HS(h, s, *p, *t)

where:

Function is used to find parameters in region 2 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

231.                     Properties calculation result in IF-97 region 2 as function of density r [kg/m3], specific enthalpy h [J/kg]:

wspPT2RH(r, h, *p, *t)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (r, h) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (r, h) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

232.                     Properties calculation result in IF-97 region 3 as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPT3HS(h, s, *p, *t)

where:

Function is based on function wspRT3HS. After calling wspRT3HS the standard equations are used to determine the pressure. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

233.                     Properties calculation result in IF-97 region 5 as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspPT5HS(h, s, *p, *t)

where:

Function is used to find parameters in region 5 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

234.                     Properties calculation result in IF-97 region 5 as function of density r [kg/m3], specific enthalpy h [J/kg]:

wspPT5RH(r, h, *p, *t)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (r, h) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (r, h) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

235.                     Density in IF-97 region 3 [kg/m3] as function of pressure p [Pa], temperature t [K]:

wspR3PT(p, t)

where:

Function calculates the density in region 3 by calling the function wspR3PTR0 with the corresponding initial estimate for density. Used for unification of the functions in all IF-97 regions. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

236.                     Density in IF-97 region 3 [kg/m3] as function of pressure p [Pa], temperature t [K], initial density r0 [kg/m3]:

wspR3PTR0(p, t, r0)

where:

Function used the Newton method with initial estimate for density to calculate the density. Used for unification of the functions in all IF-97 regions. Since at some stages the function use iterations, to speed up calculations you may vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

237.                     Properties calculation result in IF-97 region 3 as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspRT3HS(h, s, *r, *t)

where:

Function is used to find parameters in region 3 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

Note: In Mathcad the function have only parameters (h, s) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (h, s) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

238.                     Specific entropy in IF-97 region 1 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspS1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

239.                     Specific entropy in IF-97 region 2 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspS2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

240.                     Specific entropy in IF-97 region 3 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspS3PT(p, t)

where:

Function is based upon functions wspR3PT and wspU3RT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

241.                     Specific entropy in IF-97 region 3 [J/(kgK)] as function of density r [kg/m3], temperature t [K]:

wspS3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

242.                     Specific entropy in IF-97 region 5 [J/(kgK)] as function of pressure p [Pa], temperature t [K]:

wspS5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

243.                     Temperature in IF-97 region 1 [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspT1HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for Pressure as a Function of Enthalpy and Entropy p(h,s) to the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2001). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

244.                     Temperature in IF-97 region 1 [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT1PH(p, h)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

245.                     Temperature in IF-97 region 1 [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT1PS(p, s)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

246.                     Temperature at line between areas 2 and 3 [K] as function of pressure p [Pa]:

wspT23P(p)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. Used in function wspWATERSTATEAREA when region of IF-97 is determined.

247.                     Temperature in IF-97 region 2a [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT2APH(p, h)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

248.                     Temperature in IF-97 region 2a [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT2APS(p, s)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

249.                     Temperature in IF-97 region 2b [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT2BPH(p, h)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

250.                     Temperature in IF-97 region 2b [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT2BPS(p, s)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

251.                     Temperature in IF-97 region 2c [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT2CPH(p, h)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

252.                     Temperature in IF-97 region 2c [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT2CPS(p, s)

where:

Function is based upon additional equations of the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

253.                     Temperature in IF-97 region 2 [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspT2HS(h, s)

where:

Function is based upon the "Supplementary Release on Backward Equations for Pressure as a Function of Enthalpy and Entropy p(h,s) to the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2001). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

254.                     Temperature in IF-97 region 2 [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT2PH(p, h)

where:

Function is based upon functions wspT2APH, wspT2BPH and wspT2CPH. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

255.                     Temperature in IF-97 region 2 [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT2PS(p, s)

where:

Function is based upon functions wspT2APS, wspT2BPS and wspT2CPS. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

256.                     Temperature in IF-97 region 3 [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspT3HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

257.                     Temperature in IF-97 region 3 [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT3PH(p, h)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

258.                     Temperature in IF-97 region 3 [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT3PS(p, s)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

259.                     Temperature in IF-97 region 3 [K] as function of density r [kg/m3], specific enthalpy h [J/kg]:

wspT3RH(r, h)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

260.                     Temperature in IF-97 region 5 [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspT5HS(h, s)

where:

Function is used to find the temperature in region 5 by Newton method. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

261.                     Temperature in IF-97 region 5 [K] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspT5PH(p, h)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

262.                     Temperature in IF-97 region 5 [K] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspT5PS(p, s)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. It uses the Newton method to determine the roots of the function with two arguments. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

263.                     Temperature at boundary line between areas 2 and 3 [K] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspTB23HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004).

264.                     Thermal conductivity coefficient [W/(mK)] as function of density r [kg/m3], temperature t [K]:

wspTHERMCONDRT(r, t)

where:

Function is based upon the Revised Release on the IAPS Formulation 1985 for the Thermal Conductivity of Ordinary Water Substance with correction for ITS-90 (International Temperature Scale). The range of validity is: temperature from 0.01 to 800C for pressure up to 40 MPa, from 0.01 to 650C for pressure from 40 to 70 MPa, from 0.01 to 500 C for pressure from 70 to 100 MPa.

265.                     Specific internal energy in IF-97 region 1 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspU1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

266.                     Specific internal energy in IF-97 region 2 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspU2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

267.                     Specific internal energy in IF-97 region 3 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspU3PT(p, t)

where:

Function is based upon functions wspR3PT and wspU3RT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

268.                     Specific internal energy in IF-97 region 3 [J/kg] as function of density r [kg/m3], temperature t [K]:

wspU3RT(r, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

269.                     Specific internal energy in IF-97 region 5 [J/kg] as function of pressure p [Pa], temperature t [K]:

wspU5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

270.                     Specific volume in IF-97 region 1 [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspV1PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

271.                     Specific volume in IF-97 region 2 [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspV2PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

272.                     Specific volume in IF-97 region 3 [m3/kg] as function of specific enthalpy h [J/kg], specific entropy s [J/(kgK)]:

wspV3HS(h, s)

where:

Function is based upon "Supplementary Release on Backward Equations for p(h, s) for Region 3, Equations as a Function of h and s for the Region Boundaries, and a Equation Tsat(h, s) for Wet Steam of the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam" (September 2004). Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

273.                     Specific volume in IF-97 region 3 [m3/kg] as function of pressure p [Pa], specific enthalpy h [J/kg]:

wspV3PH(p, h)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

274.                     Specific volume in IF-97 region 3 [m3/kg] as function of pressure p [Pa], specific entropy s [J/(kgK)]:

wspV3PS(p, s)

where:

Function is based upon IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam. The function uses the Newton method to determine the roots. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

275.                     Specific volume in IF-97 region 3 [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspV3PT(p, t)

where:

Function is based upon function wspR3PT. Since at some stages the function use iterations, to speed up calculations you may disable precision mode (functions wspGETTOLERANCEMODE and wspSETTOLERANCEMODE) or vary relative precision for internal iterations (functions wspGETTOLERANCE and wspSETTOLERANCE).

276.                     Specific volume in IF-97 region 5 [m3/kg] as function of pressure p [Pa], temperature t [K]:

wspV5PT(p, t)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam.

277.                     Properties calculation result in IF-97 region 1 as function of pressure p [Pa], temperature t [K]:

wspVUSHCVWDERPT1PT(p, t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam and returns the set of properties that speed up calculation considerably.

Note: In Mathcad the function have only parameters (p, t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (p, t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

278.                     Properties calculation result in IF-97 region 2 as function of pressure p [Pa], temperature t [K]:

wspVUSHCVWDERPT2PT(p, t, *v, *u, *s, *h, *Cv, *w, *DVDPt, *DUDPt, *DSDPt, *DHDPt, *DVDTp, *DUDTp, *DSDTp, *DHDTp)

where:

Function is based upon the IAPWS Industrial Formulation 1997 for Thermodynamic Properties of Water and Steam and returns the set of properties that speed up calculation considerably.

Note: In Mathcad the function have only parameters (p, t) and return an array containing output (return) parameters in SI units system.

Note: In Excel the function have only parameters (p, t) and return an array which contains output (return) parameters in SI units system. By default you will receive only first output parameter (only one element from array). To retrieve all values in array do next: 1) Enter formula with this function to any cell (as example B2). 2) Select cells: first cell with inputted formula and some cells on right hand (as example 2 cells from B2 to C2). 3) Press F2. 4) Press Ctrl+Shift+Enter. If you want to retrive values in vertical way (as example for cells from B2 to B3) please use the Excel built-in "TRANSPOSE" function.

279.                     Properties calculation result in IF-97 region 3 as function of pressure p [Pa], temperature t [K]: