Critical heat flux of water in vertical round tubes at low pressure and low flow conditions

An experimental study on critical heat flux (CHF) has been performed for water flow in vertical round tubes under low pressure and low flow (LPLF) conditions to provide a systematic data base and to investigate parametric trends. Totally 513 experimental data have been obtained with Inconel-625 tube test sections in the following conditions: diameter of 6, 8, 10 and 12 mm; heated length of 0.31.77 m; pressure of 106951 kPa; mass flux of 20277 kg m⁻² s⁻¹; and inlet subcooling of 50654 kJ kg⁻¹, thermodynamic equilibrium critical quality of 0.3231.251 and CHF of 1081598 kW m⁻². Flow regime analysis based on Mishima & Ishii’s flow regime map indicates that most of the CHF occurred due to liquid film dryout in annular-mist and annular flow regimes. Parametric trends are examined from two different points of view: fixed inlet conditions and fixed exit conditions. The parametric trends are generally consistent with previous understandings except for the complex effects of system pressure and tube diameter. Finally, several prediction models are assessed with the measured data; the typical mechanistic liquid film dryout model and empirical correlations of (Shah, M.M., 1987. Heat Fluid Flow 8 (4), 326–335; Baek, W.P., Kim, H.G., Chang, S.H., 1997. KAIST critical heat flux correlation for water flow in vertical round tubes, NUTHOS-5, Paper No. AA5) show good predictions. The measured CHF data are listed in Appendix B for future reference.     

Effect of pressure on critical heat flux in uniformly heated vertical annulus under low flow conditions

Critical heat flux (CHF) experiments have been carried out in a wide range of pressure for an internally heated vertical annulus. The experimental conditions covered a range of pressure from 0.57 to 15.01 MPa, mass fluxes of 0 kg m⁻² s⁻¹ and from 200 to 650 kg m⁻² s⁻¹, and inlet subcoolings from 85 to 413 kJ kg⁻¹. Most of the CHFs were identified to the dryout of the liquid film in the annular-mist flow. For the mass fluxes of 550 and 650 kg m⁻² s⁻¹, the CHFs had a maximum value at a pressure of 2–3 MPa, and the pressure at the maximum CHF values had a trend moving toward the pressure at the peak value of pool boiling CHF as the mass flux decreased. The CHF data under a zero mass flux condition indicate that both the effects of pressure and inlet subcooling on the CHF were smaller, compared with those for the CHF with a net water upflow. The Doerffer correlation using the 1995 CHF look-up table and the Bowring correlation show a good prediction capability for the present CHF data.     

Critical heat flux under flow oscillation of water at low-pressure, low-flow conditions

In this work, the effect of flow oscillations on critical heat flux (CHF) is investigated for water flow in vertical round tubes at low-pressure, low-flow (LPLF) conditions. An experimental study has been conducted to investigate the difference in CHF between forced and natural circulations, and between stable and oscillating flow conditions with three vertical round tube test sections (5.0 mm ID×0.6 m in length, 6.6 mm ID×0.5 m in length, and 9.8 mm ID×0.6 m in length) for mass fluxes below 400 kg m⁻² s⁻¹ under near atmospheric pressure. It is found that flow oscillations can drastically reduce the CHF, in particular for natural-circulation conditions. In addition to the experiments, CHF correction factors for flow oscillation effects are developed for forced and natural circulations, respectively, based on the experimental data of the present work and others.     

Critical heat flux under choking flow conditions: Part II – Maximum values of flow parameters attained under choking flow conditions

Under conditions of forced convective boiling at low pressures and high mass fluxes, beyond a certain quality, choking flow may occur at the exit of a heated channel. An experimental investigation carried out by Olekhnovitch et al. (Olekhnovitch, A., Teyssedou, A., Tye, P., Champagne, P., 2000. Critical heat flux under choking flow conditions. Part I — Outlet pressure fluctuations. Nucl. Eng. Des., this issue) has shown that the occurrence of choking flow does not radically influence the values of the critical heat flux (CHF). However, once the choking flow conditions have occurred, for a given mass flux and quality, the outlet pressure cannot be lowered below a certain value that is fixed by the flow itself. A model that allows this pressure to be determined and which must be used in conjuction with correlations for the prediction of the CHF is presented.     

Critical heat flux under zero flow conditions in vertical annulus with uniformly and non-uniformly heated sections

The experimental study of water CHF (critical heat flux) under zero flow conditions has been carried out in an annulus flow channel with uniformly and non-uniformly heated sections over a pressure range of 0.52–14.96 MPa. In the present boiling system, the CHFs occur in the upper region of the heated section, in contrast to the results in the experiments for boiling tubes conducted by several investigators. The general trend of the CHF with pressure is that the CHF increases up to a medium pressure of about 6–8 MPa and decreases as the pressure is further increased. A comparison of the present data with the existing flooding CHF correlations shows that the correlations depend greatly on the effect of the heat flux distribution. When the correction terms with the density ratio and the effect of the heat flux distribution proposed in the present work are used with the CHF correlation based on the Wallis flooding correlation, it predicts the measured flooding CHF within an RMS error of 9.0%.     

Critical heat flux prediction in rod bundles under upward low mass flux densities

The analysis of experimental data and results of calculations for heat transfer crisis in heated channels under low upward coolant mass flux densities is presented. This analysis allows the determination of the basic features of the boiling crisis phenomenon. It is shown that the methods currently used for critical heat flux (CHF) prediction have insufficient accuracy in the given range of parameters. A new relationship for the CHF calculation is presented. It should be used for the water–water energy reactor (WWER) and uran–graphite channel reactor—Chernobyl-type (RBMK) rod bundles, and is verified by the test data. The comparison of results obtained by a new CHF correlation and the relationship used in RELAP5/MOD3.1 Code is presented. It is shown that the latter overpredicts the CHF values at atmospheric pressure and for x_cr>0.4 and does not provide conservative estimations for the RBMK fuel bundles.     

An experimental and analytical investigation of uncovered core heat transfer under high pressure, low heat flux conditions

Experimental and analytical results are reported from two series of high pressure core uncovering experiments. It was determined that the uncovered core is cooled primarily by convection and radiation to dry steam and that droplets are confined to the immediate vicinity of the mixture level. Spacer grids substantially increased heat transfer at and downstream of the grid. A simple heat transfer model is presented which accurately predicts uncovered core heat transfer at modified wall Reynolds numbers greater than 2000. Results are expected to be use in modelling small break loss of coolant accidents.     

Critical heat flux and flow pattern for water flow in annular geometry

An experimental study on critical heat flux (CHF) and two-phase flow visualization has been performed for water flow in internally-heated, vertical, concentric annuli under near atmospheric pressure. Tests have been done under stable forced-circulation, upward and downward flow conditions with three test sections of relatively large gap widths (heated LENGTH = 0.6 m, inner DIAMETER = 19 mm, outer DIAMETER = 29, 35 and 51 mm). The outer wall of the test section was made up of the transparent Pyrex tube to allow the observation of flow patterns near the CHF occurrence. The CHF mechanism was changed in the order of flooding, churn-to-annular flow transition and local dryout under a large bubble in churn flow as the flow rate was increased from zero to higher values. Observed parametric trends are consistent with the previous understanding except that the CHF for downward flow is considerably lower than that for the upward flow. In addition to the experiment, selected CHF correlations for annuli are assessed based on 1156 experimental data from various sources. The Doerffer et al. (1994); Barnett (1966); Jannsen and Kervinen (1963); Levitan and Lantsman (1977) correlations show reasonable predictions for wide parameter ranges, among which the Doerffer et al. (1994) correlation shows the widest parameter ranges and a possibility of further improvement. However, there is no correlation predicting the low-pressure, low-flow CHF satisfactorily.     

Two-phase flow regimes and mechanisms of critical heat flux under subcooled flow boiling conditions

A literature review of critical heat flux (CHF) experimental visualizations under subcooled flow boiling conditions was performed and systematically analyzed. Three major types of CHF flow regimes were identified (bubbly, vapor clot and slug flow regime) and a CHF flow regime map was developed, based on a dimensional analysis of the phenomena and available experimental information. It was found that for similar geometric characteristics and pressure, a Weber number (We)/thermodynamic quality (x) map can be used to predict the CHF flow regime.Based on the experimental observations and the review of the available CHF mechanistic models under subcooled flow boiling conditions, hypothetical CHF mechanisms were selected for each CHF flow regime, all based on a concept of wall dry spot overheating, rewetting prevention and subsequent dry spot spreading. Even though the selected concept has not received much attention (in term or theoretical developments and applications) as compared to other more popular DNB models, its basis have often been cited by experimental investigators and is considered by the authors as the “most-likely” mechanism based on the literature review and analysis performed in this work. The selected modeling concept has the potential to span the CHF conditions from highly subcooled bubbly flow to early stage of annular flow and has been numerically implemented and validated in bubbly flow and coupled with one- and three-dimensional (CFD) two-phase flow codes, in a companion paper. [Le Corre, J.M., Yao, S.C., Amon, C.H., in this issue. A mechanistic model of critical heat flux under subcooled flow boiling conditions for application to one and three-dimensional computer codes. Nucl. Eng. Des.].     

Analysis of the critical heat flux in round vertical tubes under low pressure and flow oscillation conditions. Applications of artificial neural network

Artificial neural networks (ANNs) for predicting critical heat flux (CHF) under low pressure and oscillation conditions have been trained successfully for either natural circulation or forced circulation (FC) in the present study. The input parameters of the ANN are pressure, mean mass flow rate, relative amplitude, inlet subcooling, oscillation period and the ratio of the heated length to the diameter of the tube, L/D. The output is a nondimensionalized factor F, which expresses the relative CHF under oscillation conditions. Based on the trained ANN, the influences of principal parameters on F for FC were analyzed. The parametric trends of the CHF under oscillation obtained by the trained ANN are as follows: the effects of pressure below 500 kPa are complex due to the influence of other parameters. F will increase with increasing mean mass flow rate under any conditions, and will decrease generally with an increase in relative amplitude. F will decrease initially and then increase with increasing inlet subcooling. The influence curves of mean mass flow rate on F will be almost the same when the period is shorter than 5.0 s or longer than 15 s. The influence of L/D will be negligible if L/D>200. It is found that the minimum number of neurons in the hidden layer is a product of the number of neurons in the input layer and in the output layer.     

Liquid metal heat transfer in heat exchangers under low flow rate conditions

The present paper describes the liquid metal heat transfer in heat exchangers under low flow rate conditions. Measured data from some experiments indicate that heat transfer coefficients of liquid metals at very low Péclet number are much lower than what are predicted by the well-known empirical relations. The cause of this phenomenon was not fully understood for many years. In the present study, one countercurrent-type heat exchanger is analyzed using three, separated countercurrent heat exchanger models: one is a heat exchanger model in the tube bank region, while the upper and lower plena are modeled as two heat exchangers with a single heat transfer tube. In all three heat exchangers, the same empirical correlation is used in the heat transfer calculation on the tube and the shell sides. The Nusselt number, as a function of the Péclet number, calculated from measured temperature and flow rate data in a 50 MW experimental facility was correctly reproduced by the calculation result, when the calculated result is processed in the same way as the experiment. Finally, it is clarified that the deviation is a superficial phenomenon which is caused by the heat transfer in the plena of the heat exchanger.     

An utilization of liquid sublayer dryout mechanism in predicting critical heat flux under low pressure and low velocity conditions in round tubes

From a theoretical assessment of extensive critical heat flux (CHF) data under low pressure and low velocity (LPLV) conditions, it was found out that lots of CHF data would not be well predicted by a normal annular film dryout (AFD) mechanism, although their flow patterns were identified as annular–mist flow. To predict these CHF data, a liquid sublayer dryout (LSD) mechanism has been newly utilized in developing the mechanistic CHF model based on each identified CHF mechanism. This mechanism postulates that the CHF occurrence is caused by dryout of the thin liquid sublayer resulting from the annular film separation or breaking down due to nucleate boiling in annular film or hydrodynamic fluctuation. In principle, this mechanism well supports the experimental evidence of residual film flow rate at the CHF location, which can not be explained by the AFD mechanism. For a comparative assessment of each mechanism, the CHF model based on the LSD mechanism is developed together with that based on the AFD mechanism. The validation of these models is performed on the 1406 CHF data points ranging over P=0.1–2 MPa, G=4–499 kg m⁻² s⁻¹, L/D=4–402. This model validation shows that 1055 and 231 CHF data are predicted within ±30 error bound by the LSD mechanism and the AFD mechanism, respectively. However, some CHF data whose critical qualities are <0.4 or whose tube length-to-diameter ratios are <70 are considerably overestimated by the CHF model based on the LSD mechanism. These overestimations seem to be caused by an inadequate CHF mechanism classification and an insufficient consideration of the flow instability effect on CHF. Further studies for a new classification criterion screening the CHF data affected by flow instabilities as well as a new bubble detachment model for LPLV conditions, are needed to improve the model accuracy.     

Experimental investigation on critical heat flux in horizontal channel for low pressure low flow (LPLF) condition

Studies reported in the past on critical heat flux (CHF) are mostly limited to vertical flow, large channel diameter, high pressure and high mass flux. Only few investigations are reported in the literature for horizontal flow CHF especially under low pressure and low flow conditions. Hence, predictive methods of CHF for horizontal flow are scarce. There is a need for understanding CHF in horizontal flow under low pressure and low flow conditions because they are commonly encountered in nuclear reactor fuel channels of pressurized heavy water reactor (PHWR) under loss of coolant accidental (LOCA) conditions. The present work investigates CHF of horizontal flow for low flow rates (mass flux of 100–400 kg/m² s) at nearly atmospheric pressure conditions. Parameters covered in this study are diameter (5.5 mm, 7.5 mm and 9.5 mm), length (0.45 m and 0.8 m) and a inlet temperature of 32 °C. The first occurrence of ‘red hot’ spot on the test section is considered as the onset of critical heat flux condition in the present work. Experimental results obtained are compared with Groeneveld et al. (2007) look up table data for vertical flow after applying correction factor given by Wong et al. (1990). The deviation of experimental CHF data from those predicted using Groeneveld et al. (2007) look up table and Wong et al. (1990) correction factor is more than 50%.     

Critical heat flux and heat transfer above mixture level under high-pressure boil-off conditions in PWR type and tight-lattice type fuel bundles

Heat transfer tests were conducted in PWR 17 × 17 type and tight-lattice type fuel bundles under high-pressure boil-off (very-low flow, mass fluxes lower than 100 kg/m²s) conditions. There is almost no significant difference in both critical heat flux (CHF) (or dryout point) data and convective heat transfer data above the mixture level between the PWR type and tight-lattice type bundles. The “complete vaporization equation” predicts well the CHF data, i.e. the dryout occurs nearly at the elevation where the thermal-equilibrium quality reaches 1.0. The Groeneveld CHF table used in the RELAP5/MOD3 code should be improved in the region of mass flux between 10 and 100 kg/m²s. The radiative heat transfer has an important contribution to total heat transfer above the mixture level. The Dittus-Boelter correlation, with use of the film temperature in evaluating steam properties, predicts well the convective heat transfer above the mixture level.     

Critical heat flux at high velocity channel flow with high subcooling

A quantitative analysis of critical heat flux (CHF) in heated channels under high mass flux with high subcooling was successfully carried out by applying a new flow model to the existing CHF model of a macro-water-sublayer on the heated wall and steam blankets over it. The CHF correlation proposed could correctly predict the existing experimental data for circular tubes of 0.33–4 mm in diameter with mass flux of 124–90 000 kg (m² s)⁻¹ and inlet water subcooling of 35–210 K at 0.1–7.1 MPa, resulting in CHF of 4.2–224 MW m⁻², and for rectangular channels of 3–20 mm gap with a mass flux of 940–27 000 kg (m² s)⁻¹ and inlet water subcooling of 13–166 K at 0.1–3.0 MPa, resulting in CHF of 2.0–62 MW m⁻². An error of the CHF correlation has also been estimated.     

Analysis of critical heat flux during subcooled boiling for finned fuel elements

Experimental and analytical studies were performed to determine the critical heat flux (CHF) during subcooled boiling on finned fuel elements. Tests were conducted in a vertical, concentric-annulus test section consisting of a glass tube containing a finned heater element with either six, eight, or ten longitudinal fins. The phenomena leading to CHF are described and the parametric trends are discussed.A two-dimensional finite-element heat transfer model using the Galerkin method was used to analyse the experimental data to obtain CHF values. A dimensionless correlation was derived to predict the CHF values during subcooled boiling. Over 90% of the predicted CHF values agreed with those obtained from the two-dimensional analysis within ±30%.     

Critical heat flux data in a vertical tube at low and medium pressures

AECL Research and École Polytechnique have been cooperating on the validation of the critical heat flux (CHF) look-up table (D.C. Groeneveld et al., Heat Transfer Eng. 7(1–2) (1986) 46–62). For low and medium pressures the values in the table have been obtained by extrapolation and curve fitting; therefore, errors could be expected. To reduce these possible extrapolation errors, CHF experiments are being carried out in water cooled 8 mm internal diameter (ID) tubes, at conditions where the data are scarce. This paper presents some of the experimental CHF data obtained for vertical up flow in an 8 mm ID test section, for a wide range of exit qualities (5–70%) and the exit pressure ranging from 5 to 30 bar. The experiments were carried out for heated lengths of 0.75, 1, 1.4 and 1.8 m. In general, the collected data show parametric trends similar to those described in the open literature. However, it was observed that for low pressure conditions CHF depends on the heated length; this dependence begins to disappear for exit pressure of about 30 bar. The CHF data have also been compared with predictions of well-known correlations (L. Biasi et al., Energia Nucl. 14(9) (1967) 530–536; R. Bowring, Br. Report AEEW-R789, Winfrith, UK, 1972; Y. Khatto and H. Ohno, Int. J. Heat Mass Transfer 27 (1984) 1641–1648) and those of the look-up table given by Groeneveld et al. For low pressures and low mass fluxes the look-up table seems to yield better predictions of the CHF than the correlations. However, for medium pressures and mass fluxes the correlations perform better than the look-up table; among those tested, Katto and Ohno's correlation gives the best results.     

Void fraction under high pressure, low flow conditions in rod bundle geometry

Void-fraction data are reported from a series of high pressure, low heat and mass flux experiments. Testing was performed in a heated vertical rod bundle with internal dimensions similar to a PWR fuel bundle. The results are of interest in analyses of small break loss of coolant accidents. The experiments showed that, at a given pressure, void-fraction data could be fitted to a drift-flux equation with a constant drift-velocity. The drift-velocity was observed to decrease with increasing pressure and was independent of void fraction; a characteristic normally associated with churn-turbulent flow. However, relevant drift-flux correlations found in the literature gave relatively poor predictions of void fraction. The best predictions were obtained from an empirical correlation based on dimensional analysis.     

A mechanistic model of critical heat flux under subcooled flow boiling conditions for application to one- and three-dimensional computer codes

Based on a review of visual observations at or near critical heat flux (CHF) under subcooled flow boiling conditions and consideration of CHF triggering mechanisms, presented in a companion paper [Le Corre, J.M., Yao, S.C., Amon, C.H., 2010. Two-phase flow regimes and mechanisms of critical heat flux under subcooled flow boiling conditions. Nucl. Eng. Des.], a model using a two-dimensional transient thermal analysis of the heater undergoing nucleation was developed to mechanistically predict CHF in the case of a bubbly flow regime. The model simulates the spatial and temporal heater temperature variations during nucleation at the wall, accounting for the stochastic nature of the boiling phenomena. It is postulated that a high local wall superheat occurring underneath a nucleating bubble at the time of bubble departure can prevent wall rewetting at CHF (Leidenfrost effect). The model has also the potential to evaluate the post-DNB heater temperature up to the point of heater melting.Validation of the proposed model was performed using detailed measured wall boiling parameters near CHF, thereby bypassing most needed constitutive relations. It was found that under limiting nucleation conditions; a peak wall temperature at the time of bubble departure can be reached at CHF preventing wall cooling by quenching. The simulations show that the resulting dry patch can survive the surrounding quenching events, preventing further nucleation and leading to a fast heater temperature increase. The model was applied at CHF conditions in simple geometry coupled with one-dimensional and three-dimensional (CFD) codes. It was found that, within the range where CHF occurs under bubbly flow conditions (as defined in Le Corre et al., 2010), the local wall superheat underneath nucleating bubbles is predicted to reach the Leidenfrost temperature. However, a better knowledge of statistical variations in wall boiling parameters would be necessary to correctly capture the CHF trends with mass flux (or Weber number).     

An experimental study on the critical heat flux for low flow of water in a non-uniformly heated vertical rod bundle over a wide range of pressure conditions

An experimental study of the critical heat flux (CHF) has been performed for a water flow in a non-uniformly heated vertical 3 × 3 rod bundle under low flow and a wide range of pressure conditions. The experiment was especially focused on the parametric trends of the CHF and the applicability of the conventional CHF correlations to a return-to-power conditions of a main steam line break accident whose conditions might be a low mass flux, intermediate pressure, and a high inlet subcooling. The effects of the mass flux and pressure on the CHF are relatively large and complicated in the low pressure conditions. At a high mass flux or a low critical quality, the local heat flux at the CHF location sharply decreases with an increasing local critical quality. However, at a low mass flux or a high critical quality, the local heat flux at the CHF location shows a nearly constant value regardless of the increase of the critical quality. The CHF data at the very low mass flux conditions are correlated well by the churn-to-annular flow transition criterion or the flow reversal phenomena. Several conventional CHF correlations predict the present return-to-power CHF data with reasonable accuracies. However, the prediction capabilities become worse in a very low mass flux of below about 100 kg/(m² s).