Prepare 
& Submit Calculation Form
The thermodynamic and transport properties of  heavy water can be calculated according to the property equations presented in: 
[1] Hill P.G., MacMillan R.D.C., Lee V., "Tables of Thermodynamic Properties of Heavy Water in S.I. Units", Department of Mechanical Engineering, The University of British Columbia, Atomic Energy of Canada Limited, AECL-7531, 1981.
[2] Hill P.G., MacMillan R.D.C., Lee V., "A Fundamental Equation of State for Heavy Water", Journal of Physical Chemistry Reference Data, 11(1):1, 1982; Errata: 12(4):1065, 1983.
[3] Matsunaga N., Nagashima A., "Transport Properties of Liquid and Gaseous Deuterium Oxide over a Wide Range of Temperature and Pressure", Journal of Physical Chemistry Reference Data, 12(4):933, 1983.
[4] Straub J., Rosner N. und Grigull U., "Oberflaechenspannung von Leichtem und Schwerem Wasser", Waerme- und Stoffuebertragung, 13:241, 1980.
Based mainly on the above references property equations have been implemented in the HEAVY_WASP (HEAVY Water And Steam Properties) FORTRAN 77 code in the course of the following PhD work (in Greek):
Petropoulos N.P., "THERMOPHYSICAL PROPERTIES OF NUCLEAR POWER REACTOR COOLANTS: COMPARATIVE STUDY AND CALCULATION CODES", Nuclear Engineering Section, Mechanical Engineering School, National Technical University of Athens, 2000 (in progress).
The set of equations used provide reliable thermodynamic and transport properties results for the industrial use of heavy water. The following validity range is covered:
Pressure: 0.006601 - 1000 bar
Temperature: 3.8 - 500 
oC
 
A small region at the vicinity of the critical 
point is excluded:
 
|Ô - Ô*| £ 10 Ê, |ñ/ñ* - 1| £ 0.3
where T*= 643.89 K and ñ*= 358.0 kg/m3
Calculable Properties:
| Specific volume v in m3/kg | Density ñ in kg/m3 | 
| Internal energy u in kJ/kg | Specific enthalpy h in kJ/kg | 
| Entropy s in kJ/(kgK) | Isobaric heat capacity cp in kJ/(kgK) | 
| Isochoric heat capacity cv in kJ/(kgK) | Viscosity ç in kg/(ms) | 
| Thermal conductivity ë in W/(mK) | Surface tension ó in N/m. | 
Code accepts the following input pair:
(P,T) in bars and degrees Centigrade respectively.
Addendum 
This WWW page 
is maintained regularly to ensure best performance and error-free accurate 
calculations for the thermodynamic and transport properties of light 
(ordivary)  water. All efforts have been made to ensure that the 
present  code  is as accurate implementation as possible of the 
property equations included at the above literature. Authors welcome your 
comments and suggestions for further improvement along with any bug reports and 
they will acknowledge any reported bug in this very page under the link reported 
bugs regardless whether they are fixed or not.  
Please, contact:
Nick P. Petropoulos (npetro@nuclear.ntua.gr)
Research Assistant
Nuclear Engineering Section
Mechanical Engineering Dept.
National Technical University of Athens
15780 Athens, GREECE
tel: +302107722939
fax: +302107722914
URL: http://arcas3.nuclear.mech.ntua.gr/
However, authors will not try to fix all reported bugs 
promptly, and they will certainly not be ivolved in an effort to improve 
performance of the property equations in pressure - temperature regions, where 
their results are ambiguous or not to be trusted.   Furthermore, 
authors cannot accept ANY responsibility if some or all calculations based on 
this code may be proven wrong or inadequate.Visitors of this page, who would 
like to receive free copies of the source code for their own private use or for 
evaluating purposes are encouraged to place a request at npetro@nuclear.ntua.gr, stating also 
their intented implementation platform (i.e. FORTRAN make & version, UNIX or 
Windows based application etc). Not all requests are to be granted. This 
version of the software is not commercially available. 
Back to the 
Previous Page