NEA-0682 H2OTP.

1. NAME OR DESIGNATION OF PROGRAM:  H2OTP.

3. DESCRIPTION OF PROBLEM OR FUNCTION

H2OTP is a FORTRAN subroutine for the calculation of a full set of light-water thermodynamic and transport properties after the 1967 IFC recommendations for industrial use are calculated.

The following specific parameters can be calculated for steam- or water or saturated-conditions:
     enthalpy, energy, volume, entropy
     Gibbs function, Helmholtz function
     relative isobaric volume expansion
     isobaric and isochoric heat capacity
     isentropic exponent, sonic velocity
     relative isothermal compressibility
     Joule and Joule-Kelvin coefficient
     heat of evaporation
     saturation temperature given by pressure
     saturation pressure given by temperature
     surface tension
     dynamic and kinematic viscosity
     thermal conductivity and Prandtl number

The independent variables are temperature and pressure, British or S.I. units.

4. METHOD OF SOLUTION

The data source is the 1967 IFC formulations for Industrial use. The nature of this formulation is of the form of a highly complex mathematical expression for the Gibbs (or Helmholz) function dependent on the temperature and pressure (spc. volume). The thermodynamic quantities are derived from this basic function as simple functions of the derivatives in temperature and pressure up to the second order of differentiation.
The transport properties are calculated from the recommended expressions in the 1967 IFC formulation. These expressions are numerical approximations to the experimental results. The set of transport properties calculations presented by the 1967 IFC formulation is not complete, so it has been found necessary to complement with some approximations taken from ref. Risoe-M-669 below, and with some approximations made especially for this code.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The temperature range is 0.01 deg. C to 800 deg. C and the pressure range is 610 N/m**2 to 500.0E5 (1000.0E5) N/m**2.
The use of double-precision words with a mantissa length of 72 bits, corresponding to approximately 22 decimal digits on the Burroughs B-7800 computer of Riso National Laboratory, is essential  and has made it possible to use numerical differentiation. Computers of other makes with a similar or higher accuracy can be applied.

6. TYPICAL RUNNING TIME

The calculations are slow and are meant as a reference base. The CPU-time per call of the subroutine H2OTP is approximately 0.025 s on the Burroughs B-7800.

NEA-DB executed the test case on VAX-11/780 in 1.1 sec. of CPU time.

7. UNUSUAL FEATURES OF THE PROGRAM:

8. RELATED AND AUXILIARY PROGRAMS:

10. REFERENCES

- H. Schmidt,
  Properties of Water and Steam in SI-Units,
  VDI, Verein Deutscher Ingenieure, Springer 1969.
- Thermodynamic and Transport Properties of Steam,
  ASME, American Society of Mechanical Engineers, 1967.
- M. Hjelm-Hansen,
  Polynomial Approximation of Thermodynamic Properties of Light
  Water,
  Risoe-M-669, December 1967.

- F.W. Cortzen:
  H2OTP, "A Subroutine in FORTRAN for a Full Set of Light-Water
  Thermodynamic and Transport Properties after the 1967  IFC
  Recommendations for Industrial Use"
  Risoe-M-1844 (March 1976)

11. MACHINE REQUIREMENTS

To execute the test case on VAX-11/780, 35,328 bytes of main storage are required.

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED

The test case included in the package was executed on VAX-11/780 under VMS-V3.4.

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

15. NAME AND ESTABLISHMENT OF AUTHOR

          F. W. Cortzen
          Risoe National Laboratory
          Dept. of Energy Technology
          DK-4000 Roskilde
          Denmark

File name	File description	Records
NEA0682_01.003	H2OTP INFORMATION FILE	41
NEA0682_01.004	H2OTP SOURCE	3411
NEA0682_01.005	H2OTP TEST CASE OUTPUT	144