Computer Programs

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHODS, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES, AUXILIARIES, STATUS, REFERENCES, HARDWARE REQUIREMENTS, LANGUAGE, SOFTWARE REQUIREMENTS, OTHER RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES

[ top ]

[ top ]

To submit a request, click below on the link of the version you wish to order. Rules for end-users are
available here.

Program name | Package id | Status | Status date |
---|---|---|---|

D2O V.1.0 | NEA-1535/01 | Tested | 30-OCT-2000 |

Machines used:

Package ID | Orig. computer | Test computer |
---|---|---|

NEA-1535/01 | PC-80486 | PC Pentium III 500,Linux-based PC,IBM RISC6000 WS |

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

A computer program for the fast computation of the thermodynamic and transport properties of heavy water (D2O) at saturation, in subcooled liquid and superheated vapor states. Specific volume (or density), specific enthalpy, specific entropy, constant-pressure specific heat and temperature at saturation are calculated by a number of piecewise continuous approximation functions of (and their derivatives are calculated with respect to) pressure whereas pressure at saturation is calculated by a piecewise continuous approximation function of temperature for heavy water. Density in subcooled liquid state, specific volume in super-heated vapor state, specific enthalpy, specific entropy and constant-pressure specific heat in both of these states are calculated by some piecewise continuous approximation functions of pressure and temperature for heavy water.

The correlations used in the calculation of these thermodynamic properties of heavy water were derived by fitting some appropriate curves to the data given in the steam tables by Hill et al (1981). The whole set of correlations and the approximation method used in their derivation are presented by Durmayaz (1997).

Dynamic viscosity and thermal conductivity for heavy water are calculated as functions of temperature and density with the correlations given by Hill et al (1981), by Matsunaga and Nagashima (1983) and by Kestin et al (1984).

Surface tension for heavy water is calculated as a function of temperature with the correlation given by Crabtree and Siman-Tov (1993).

A computer program for the fast computation of the thermodynamic and transport properties of heavy water (D2O) at saturation, in subcooled liquid and superheated vapor states. Specific volume (or density), specific enthalpy, specific entropy, constant-pressure specific heat and temperature at saturation are calculated by a number of piecewise continuous approximation functions of (and their derivatives are calculated with respect to) pressure whereas pressure at saturation is calculated by a piecewise continuous approximation function of temperature for heavy water. Density in subcooled liquid state, specific volume in super-heated vapor state, specific enthalpy, specific entropy and constant-pressure specific heat in both of these states are calculated by some piecewise continuous approximation functions of pressure and temperature for heavy water.

The correlations used in the calculation of these thermodynamic properties of heavy water were derived by fitting some appropriate curves to the data given in the steam tables by Hill et al (1981). The whole set of correlations and the approximation method used in their derivation are presented by Durmayaz (1997).

Dynamic viscosity and thermal conductivity for heavy water are calculated as functions of temperature and density with the correlations given by Hill et al (1981), by Matsunaga and Nagashima (1983) and by Kestin et al (1984).

Surface tension for heavy water is calculated as a function of temperature with the correlation given by Crabtree and Siman-Tov (1993).

[ top ]

4. METHODS

A group of pressure-enthalpy (P-h) pairs can be given in an input data file or assigned in the main program without knowing the state in which fluid takes place. In this case, first, the enthalpies at saturation corresponding to the given pressure are computed. Second, the state is determined by comparing the given enthalpy to the saturation enthalpies. Then, the properties are computed.

Program D2O can also be linked as a subroutine of another main program. In this case, any thermodynamic property at saturation can be computed directly as a function of saturation pressure as well as saturation pressure is determined directly as a function of saturation temperature. Similarly, if a P-h pair is known in subcooled liquid or superheated vapor state, any thermodynamic or transport property can be computed directly. If a ressure- temperature (P-T) pair is known in subcooled liquid or superheated vapor state, enthalpy can also be computed directly.

A group of pressure-enthalpy (P-h) pairs can be given in an input data file or assigned in the main program without knowing the state in which fluid takes place. In this case, first, the enthalpies at saturation corresponding to the given pressure are computed. Second, the state is determined by comparing the given enthalpy to the saturation enthalpies. Then, the properties are computed.

Program D2O can also be linked as a subroutine of another main program. In this case, any thermodynamic property at saturation can be computed directly as a function of saturation pressure as well as saturation pressure is determined directly as a function of saturation temperature. Similarly, if a P-h pair is known in subcooled liquid or superheated vapor state, any thermodynamic or transport property can be computed directly. If a ressure- temperature (P-T) pair is known in subcooled liquid or superheated vapor state, enthalpy can also be computed directly.

[ top ]

[ top ]

6. TYPICAL RUNNING TIME

Running time depends on the number of the input P-h data pairs and the type of the FORTRAN compiler and computer used. For example, it is about 100 seconds on an IBM/PC-AT compatible 486DX2-66 and about 13 seconds on an IBM/PC-AT compatible Pentium 166 for 7487 input P-h data pairs (records). This sample run includes the reading of the 7487 P-h data pairs from D2OINP.DAT, the fourth test input data file, which also includes the determination of the states, the computation of the thermodynamic and transport properties, the writing of the outputs (14722 records) for superheated vapor into D2OVOU.DAT file and the writing of the outputs (409 records) for saturated vapor into D2OSOU.DAT file.

Running time depends on the number of the input P-h data pairs and the type of the FORTRAN compiler and computer used. For example, it is about 100 seconds on an IBM/PC-AT compatible 486DX2-66 and about 13 seconds on an IBM/PC-AT compatible Pentium 166 for 7487 input P-h data pairs (records). This sample run includes the reading of the 7487 P-h data pairs from D2OINP.DAT, the fourth test input data file, which also includes the determination of the states, the computation of the thermodynamic and transport properties, the writing of the outputs (14722 records) for superheated vapor into D2OVOU.DAT file and the writing of the outputs (409 records) for saturated vapor into D2OSOU.DAT file.

[ top ]

7. UNUSUAL FEATURES

D2O can be used in the two-phase thermalhydraulic system analysis of the nuclear reactors that are cooled and/or moderated by heavy water because the correlations used for the thermodynamic properties in this program are functions of pressure and temperature. Temperature can also be determined by using a root finding algorithm in case pressure and enthalpy are the independent variables.

D2O can be used in the two-phase thermalhydraulic system analysis of the nuclear reactors that are cooled and/or moderated by heavy water because the correlations used for the thermodynamic properties in this program are functions of pressure and temperature. Temperature can also be determined by using a root finding algorithm in case pressure and enthalpy are the independent variables.

[ top ]

[ top ]

10. REFERENCES

- Durmayaz, A.

Approximate functions for the fast computation of the thermo- dynamic properties of heavy water.

Nuclear Engineering and Design, Volume 178 (1997), pp. 309-329.

- Durmayaz, A.

A Loss-of Coolant Accident Analysis for Near Stagnation Flow in the Horizontal Fuel Channel of CANDU Reactor.

PhD Thesis, ITU Institute for Nuclear Energy (May 4 1995)

Istanbul

- Hill, P.G., MacMillan, R.D. and Lee, V.

Tables of thermodynamic properties of heavy water in S.E Units.

AECL-7531 (Dec. 1981).

- Kestin, J., Sengers, J.V., Kamgar-Parsi, B. Levelt-Sengers, J.M.H.

Thermophysical properties of fluid D2O.

Journal of Physical Chemistry Reference Data, Volume 13 (1984), No.

2, pp. 601-609.

- Matsunaga, N., Nagashima, A.

Transport properties of liquid and gases D20 over a wide range of

temperature and pressure.

Journal of Physical Chemistry Reference Data, Volume 12 (1983),

No.4, pp. 933-965.

- Crabtree, A., Siman-Tov, M.

Thermophysical properties of saturated light and heavy water for

advances neutron source application.

ORNL/TM-12322, May 1989

- Durmayaz, A.

Approximate functions for the fast computation of the thermo- dynamic properties of heavy water.

Nuclear Engineering and Design, Volume 178 (1997), pp. 309-329.

- Durmayaz, A.

A Loss-of Coolant Accident Analysis for Near Stagnation Flow in the Horizontal Fuel Channel of CANDU Reactor.

PhD Thesis, ITU Institute for Nuclear Energy (May 4 1995)

Istanbul

- Hill, P.G., MacMillan, R.D. and Lee, V.

Tables of thermodynamic properties of heavy water in S.E Units.

AECL-7531 (Dec. 1981).

- Kestin, J., Sengers, J.V., Kamgar-Parsi, B. Levelt-Sengers, J.M.H.

Thermophysical properties of fluid D2O.

Journal of Physical Chemistry Reference Data, Volume 13 (1984), No.

2, pp. 601-609.

- Matsunaga, N., Nagashima, A.

Transport properties of liquid and gases D20 over a wide range of

temperature and pressure.

Journal of Physical Chemistry Reference Data, Volume 12 (1983),

No.4, pp. 933-965.

- Crabtree, A., Siman-Tov, M.

Thermophysical properties of saturated light and heavy water for

advances neutron source application.

ORNL/TM-12322, May 1989

[ top ]

[ top ]

[ top ]

[ top ]

[ top ]

NEA-1535/01

OUTTP: D20VOU.DAT Sample output data file for propert. in superheated vapour stateSRCTP: COM.BL Common blocks that are included in the source program

SRCTP: D20FOR Source file

DATTP: D20INP.DAT Sample input data file for properties at saturation

OUTTP: D20LOU.DAT Sample output data file for properties in subcooled liquid state

OUTTP: D20SOU.DAT Sample output data file for properties in superheated vapour stat

MISTP: ABSTRACT.DOC Abstract file

MISTP: READ.ME Text file for instructions to users

MISTP: ReadMe.doc File for instructions to users in MS Word Version 7.0

DATTP: D20INP.DAT Third test input data file for prop. in subcooled liquid state

DATTP: D20INP.DAT First input data file for properties in liquid plase at saturatio

DATTP: D20INP.DAT Second input data file for properties in gas phase at saturation

OUTTP: D20SOU.DAT First test outp. data file for prop. in liquid phase at saturatio

OUTTP: D20SOU.DAT Second test outp. data file for prop. in gas phase at saturation

OUTTP: D20LOU.DAT 3rd test output data file for prop. in subcooled liquid state

OUTTP: D20SOU.DAT 3rd test outp. data files for prop. in liqu. phase at saturation.

OUTTP: D20SOU.DAT Fourth test data output file for prop. at gas phase at saturation

OUTTP: D20VOU.DAT Fourth test output data file for prop. in superheated vapor state

DATTP: D20INP.DAT Fourth test input data file for prop. in superheated vapor stat

Keywords: heavy water, thermodynamic properties.