Effects of turbulence and secondary flows on subcooled flow boiling

Experiments are conducted on the influence of turbulence and longitudinal vortices on subcooled flow boiling in a vertical, rectangular channel. Different flow inserts are used to create turbulence and vortices in the channel. Studied boiling regimes range from the onset of nucleate boiling over the critical heat flux up to fully developed film boiling. A wide range of measuring techniques is applied: time averaged particle image velocimetry (PIV) is used in cold flows for the evaluation of the effects the inserts have on the flow, high speed PIV and photography are used to determine the effects on the fluid and vapor movement in boiling experiments. Digital Holographic Interferometry is used for the evaluation of temperature distributions in the boiling flow. Furthermore, optical microprobes are used to obtain pointwise measurements in areas inaccessible to the imaging techniques. The experiments show that the flow inserts can have considerable impact on the heat fluxes and the distribution of vapor and temperature along the channel. All used inserts lead to an increase in critical heat flux, which is more pronounced for stronger turbulence and higher flow rates and fluid subcoolings. The measuring techniques reveal both a better transport of vapor from the heater surface as well as an increase in mixing in the liquid phase with flow inserts.     

Recent achievements in the thermal hydraulics of high heat flux components in fusion reactors

Some components of fusion thermonuclear reactors, such as divertors, plasma limiters, or first-wall armor, are believed to be subjected to operating conditions characterized by extremely high thermal loads. It is therefore necessary to remove from the surface of these components very high heat fluxes, ranging from 2 to 60 MW/m². Water subcooled flow boiling, under conditions of high mass flux, high liquid subcooling, and small to intermediate channel diameter, can accomodate these very high heat fluxes. Further enhancement of the upper limit of cooling, the critical heat flux (CHF), can be obtained by making use of turbulence promoters such as twisted tapes and coiled wires even if coupled with a relevant increase in pressure drop. An overview is presented of recent achievements obtained in water subcooled flow boiling CHF under operating conditions of interest to the thermal hydraulic design of fusion reactors. Observed basic parametric trends—CHF as a function of mass flux, pressure, subcooling, and channel geometry—are outlined, together with findings on the use of CHF enhancement techniques. From experiments it was seen that water subcooled flow boiling allows CHF conditions as high as 228 MW/m² to be achieved under extreme geometric and thermal hydraulic conditions. On the other hand, design and engineering boundary conditions limit variation in these conditions, and a suitable compromise has not yet been reached. Predictive tools are presented for the evaluation of subcooled flow boiling CHF both in straight tubes and with twisted tapes, and are assessed with reference to recent available experimental data.

Although several indications for practical applications can be found in recent achievements, a full understanding of the basic mechanisms of heat transfer and CHF in subcooled flow boiling has not yet been achieved. Future research to overcome the present lack of knowledge in this field is suggested.     

PTV experiments of subcooled boiling flow through a vertical rectangular channel

Time resolved Particle Tracking Velocimetry (PTV) experiments were carried out to investigate turbulent, subcooled boiling flow of refrigerant HFE-301 through a vertical rectangular channel with one heated wall. Measurements were performed with liquid Reynolds numbers (based on the hydraulic diameter) of Re = 3309, 9929 and 16,549 over a wall heat flux range of 0.0–64.0 kW/m². Turbulence statistics are inferred from PTV full-field velocity measurements. Quantities such as: instantaneous 2D velocity fields, time-averaged axial and normal velocities, axial and normal turbulence intensities, and Reynolds stresses are obtained. The present results agree well with previous studies and provides new information due to the full-field nature of the technique. This work is an attempt to provide turbulent subcooled boiling flow data for validation and improvement of two-phase flow computational models.     

Thermal hydraulic modeling of a natural circulation loop

 The experiment was carried out on the test loop HRTL-5, which simulates the geometry and system design of a 5 MW nuclear heating reactor. The analysis was based on a one-dimensional two-phase flow drift model with conservation equations for mass, steam, energy and momentum. Clausius–Clapeyron equation was used for the calculation of flashing front in the riser. A set of ordinary equations, which describes the behavior of two-phase flow in the natural circulation system, was derived through integration of the above conservation equations for the subcooled boiling region, bulk boiling region in the heated section and for the riser. The method of time-domain was used for the calculation. Both static and dynamic results are presented. System pressure, inlet subcooling and heat flux are varied as input parameters. The results show that subcooled boiling in the heated section and void flashing in the riser have significant influence on the distribution of the void fraction, mass flow rate and flow instability of the system, especially at low pressure. The response of mass flow rate, after a small disturbance in the heat flux is shown, and based on it the instability map of the system is given through experiment and calculation. There exists three regions in the instability map of the investigated natural circulation system, namely, the stable two-phase flow region, the unstable bulk and subcooled boiling flow region and the stable subcooled boiling and single phase flow region. The mechanism of two-phase flow oscillation is interpreted. Received on 24 January 2000     

A fractal analysis of subcooled flow boiling heat transfer

A fractal model for the subcooled flow boiling heat transfer is proposed in this paper. The analytical expressions for the subcooled flow boiling heat transfer are derived based on the fractal distribution of nucleation sites on boiling surfaces. The proposed fractal model for the subcooled flow boiling heat transfer is found to be a function of wall superheat, liquid subcooling, bulk velocity of fluid (or Reynolds number), fractal dimension, the minimum and maximum active cavity size, the contact angle and physical properties of fluid. No additional/new empirical constant is introduced, and the proposed model contains less empirical constants than the conventional models. The proposed model takes into account all the possible mechanisms for subcooled flow boiling heat transfer. The model predictions are compared with the existing experimental data, and fair agreement between the model predictions and experimental data is found for different bulk flow rates.     

Experimental study of pool temperature effects on nucleate pool boiling

Nucleate pool boiling experiments with constant wall temperature were performed using pure R113 for subcooled, saturated, and superheated pool conditions. A microscale heater array and Wheatstone bridge circuits were used to maintain the constant wall temperature and to measure the instantaneous heat flow rate accurately with high temporal and spatial resolutions. Images of bubble growth were taken at 5000 frames per second using a high-speed CCD camera synchronized with the heat flow rate measurements. The bubble geometry was obtained from the captured bubble images. The effect of the pool conditions on the bubble growth behavior was analyzed using dimensionless parameters for the initial and thermal growth regions. The effect of the pool conditions on the heat flow rate behavior was also examined. The bubble growth behaviors during subcooled, saturated, and superheated pool boiling were analyzed using a modified Jakob number that we newly defined. Dimensionless time and bubble radius parameters with the modified Jakob number characterized the bubble growth behavior well. These phenomena require further analysis for various pool temperature conditions, and this study will provide good experimental data with precise constant wall temperature boundary condition for such works.     

Numerical investigation of static flow instability in a low-pressure subcooled boiling channel

A three-dimensional two-fluid model to predict subcooled boiling flow at low pressure is presented. The model is adopted to investigate the two-phase flow and heat transfer characteristics in a heated channel. The presence of bubbles as a consequence of heating flow through a vertical rectangular channel has a significant effect on the overall pressure drop along the channel. Numerical results were compared against a series experimental data performed at various conditions – mass flux, heat flux, inlet temperature and exit pressure. Good agreement on the overall pressure drop was achieved. The onset of flow instability velocity was also accurately determined when compared against measurements. Predicted results of void fraction provided useful information towards a more fundamental understanding of the occurrence of onset of nucleate boiling, onset of significant voiding and onset of flow instability. The phenomenon of boiling onset oscillations was also predicted through the use of the two-fluid model.     

Experimental and numerical study on nucleate bubble deformation in subcooled flow boiling

Experimental and numerical study was conducted to investigate the bubble behaviors in subcooled flow nucleate boiling. The bubble behaviors in subcooled flow boiling in an upward annular channel were investigated in the range of subcooling degree 5–30 K by visualization with high spatial and temporal resolutions using a high speed video camera and Cassegrain tele-microscope. Obvious deformation on the upstream side surface of the bubble during its growth process was frequently observed. This deformation phenomenon was caused by the condensation occurring at the upstream side bottom of the bubble, which results from the Marangoni flow along the bubble surface from the bubble bottom to the top. Since the Marangoni flow cannot be directly observed by the current experiments because it occurs in a very thin interface along the bubble surface, the numerical simulations of bubble growth and departure behaviors in subcooled flow boiling were carried out. As a result, it was confirmed that the bubble deformation was caused by the Marangoni flow along the bubble surface. Moreover, the phenomenon of wave propagation on the bubble surface during the condensation process was observed, and it can enhance the heat transfer between the bubble and the surrounding subcooled liquid.     

An experimental investigation of flow boiling in an asymmetrically heated rectangular microchannel

By using unique experimental techniques and carefully constructed experimental apparatus, the characteristics of flow boiling of water in microscale were investigated using a single horizontal rectangular microchannel. A polydimethylsiloxane rectangular microchannel (D_h = 103.5 and 133 μm) was fabricated by using the replica molding technique, a kind of soft lithography. A piecewise serpentine platinum microheater array on a Pyrex substrate was fabricated with the surface micromachining MEMS technique. Real time flow visualization of the phase change phenomena inside the microchannel was performed using a high speed CCD camera with microscope. The experimental local boiling heat transfer coefficients were studied, and single bubble inception, growth, and departure, as well as elongated bubble behavior were analyzed to elucidate the microscale heat transfer mechanisms. Tests were performed for mass fluxes of 77.5, 154.9, and 309.8 kg/m² s and heat fluxes of 180–500 kW/m². The effects of mass flux, heat flux, and vapor qualities on flow boiling heat transfer in a microchannel were studied.     

Stereo-PIV measurements of vapor-induced flow modifications in confined jet impingement boiling

A single subcooled jet of water which undergoes boiling upon impingement on a discrete heat source is studied experimentally using time-resolved stereo particle image velocimetry (PIV). The impinging jet issues from a 3.75 mm diameter sharp-edged orifice in a confining orifice plate positioned 4 orifice diameters from the target surface. The behavior at jet Reynolds numbers of 5,000 and 15,000 is compared for a constant jet inlet subcooling of 10 °C. Fluorescent illumination allows for simultaneous imaging of both the flow tracers and the vapor bubbles in the flow. Flow structure, time-averaged velocities, and turbulence statistics are reported for the liquid regions within the confinement gap for a range of heat inputs at both Reynolds numbers, and the effect of the vapor generation on the flow is discussed. Vapor generation from boiling is found to modify the liquid velocities and turbulence fluctuations in the confinement gap. Flow in the confinement gap is dominated by vapor flow, and the vapor bubbles disrupt both the vertical impinging jet and horizontal wall jet flow. Moreover, vapor bubbles are a significant source of turbulence kinetic energy and dissipation, with the bubbly regions above the heated surface experiencing the most intense turbulence modification. Spectral analysis indicates that a Strouhal number of 0.023 is characteristic of the interaction between bubbles and turbulent liquid jets.     

Heat transfer at the phase interface of bubbles collapsing in subcooled liquids and during subcooled boiling

Heat transfer phenomena with bubble condensation in subcooled liquids and during subcooled boiling are discussed, which were studied by using holographic interferometry. This inertialess and non-intrusive measuring technique allows interesting insights into fluid dynamics and thermodynamics of phase change processes. The resistance for the heat transport through the phase boundaries is mainly on the liquid side and therefore the heat transfer coefficient can expressed by simple correlations following the procedure in forced convective flow. For the sake of simple boundary conditions first bubble condensation in subcooled liquids with constant and homogenous temperature are studied and than the more complex situation with large temperature gradients in the wall near boundary during subcooled boiling are discussed.Dedicated to Prof. Dieter Mewes on the occasion of his 65th birthday.     

Interface oscillation of subcooled flow boiling in locally heated microchannels

An investigation was conducted to understand flow boiling of subcooled de-ionized water in locally heated parallel microchannels. High-speed visualization technology was employed to visually observe the transient phase change process in an individual microchannel. Signal analysis method was employed in studying the interface movement and phase change process. The phase change at locally heated condition was different from those at entirely heated condition where elongated bubble(s) stayed quasi-stable for a long time without venting out. Diversified and intensive interface oscillation was observed occurring on both of the upstream and downstream bubble caps. Evaporation and condensation modes were characterized with distinguished oscillation frequencies. The film-driven oscillations of both evaporating and condensing interfaces generally operated at higher frequencies than the oscillations driven by nucleation or dropwise condensation.     

A mechanistic critical heat flux model for subcooled flow boiling based on local bulk flow conditions

The critical heat flux (CHF) mechanisms for subcooled flow boiling are reviewed. Based on experimental observations reported by previous investigators, the authors have developed a new mechanistic CHF model for vertical subcooled flow at high pressure and high mass velocity. This model is based on the dryout of a thin liquid layer (sublayer) beneath an intermittent vapor blanket due to a Helmholtz instability at the sublayer-vapor interface. The parametric trends of CHF have been explored qualitatively and quantitatively with respect to variations in pressure, mass velocity, subcooling and tube diameter. Comparisons of the model predictions with experimental data for water show good agreement in the simulation of subcooled flow conditions of pressurized water reactors (PWRs).     

Measured turbulent entropy production with large eddy particle image velocimetry

In this paper, statistical post-processing of measured velocity, dissipation rate and turbulence data is performed to establish whole-field distributions of entropy production within a channel. Thermal irreversibilities arising from temperature variations were not included in the study, as the experiments were conducted between unheated plexiglass plates in an essentially isothermal water tunnel. Unlike velocity or temperature, the measurement of entropy cannot be performed directly, so entropy production is measured indirectly through spatial differencing of measured velocities in large eddy PIV. In contrast to single-point methods of anemometry, large eddy PIV enables whole-field, time-varying measurements of the velocity field, which can be post-processed to yield entire spatial variations of the entropy production rate. An uncertainty analysis is performed to estimate measurement uncertainties with the new experimental technique. The uncertainties are decomposed into systematic and random components, including a propagated uncertainty, due to spatial differencing of the velocity field. Close comparisons between measured results of turbulence dissipation and direct numerical simulations provide useful verification of the formulation, before post-processed results of dissipation rates are used to determine entropy production within a channel.     

Experiment and prediction of a cavity type separated flow

A detail study involving flow visualization, Laser Doppler Velocimeter (LDV) measurements and numerical prediction is presented. The visualization experiments revealed striking results of a pulsatile motion in the separated flow region associated with the formation and passage of large eddy structures. Measurements of mean velocities and turbulence intensity profiles across the separated flow field, provided information about the separated shear layer development and the recirculating flow pattern. The numerical predictions, obtained with a two-layer turbulence model in conjunction with the SIMPLE algorithm, failed to reproduce the coherent eddies and the pulsatile motion, but the mean velocities are reasonably reproduced.     

Whole-field density measurements by ‘synthetic schlieren’

 This paper outlines novel techniques for producing qualitative visualisations of density fluctuations and for obtaining quantitative whole-field density measurements in two-dimensional density-stratified flows. These techniques, which utilise image processing technology, are much simpler to set up than the classical schlieren and interferometry methods, and provide useful information in situations where shadowgraph is of little or no value. Moreover, they may be set-up to analyse much larger domains than is feasible with the classical approaches, and do not require high quality optical windows in the experimental apparatus. Ultimately the greatest strength of these techniques is the ability to extract accurate, quantitative measurements of the density field. Application of these techniques is illustrated by an internal wave field produced by an oscillating cylinder. Recent theoretical advances for this classical problem make it the ideal test bed. Results are presented for both a circular and a square cylinder oscillating vertically in a linear stratification. Further aspects of the techniques are illustrated by considering thermal convection from a hand and flow over an obstacle towed through a density stratified fluid. Received: 27 October 1998 /Accepted: 24 May 1999     

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state,and the flow is frequently unsteady and turbulent.To assess the state-of-the-art of computational capabilities for unsteady cavitating flows,different cavitation and turbulence model combinations are conducted.The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics.The kε turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively.We have also blended these alternative cavitation and turbulence treatments,to illustrate that the eddy viscosity near the closure region can significantly influence the capture of detached cavity.From the experimental validations regarding the force analysis,frequency,and the cavity visualization,no single model combination performs best in all aspects.Furthermore,the implications of parameters contained in different cavitation models are investigated.The phase change process is more pronounced around the detached cavity,which is better illus-trated by the interfacial dynamics model.Our study provides insight to aid further modeling development.     

Turbulent kinetic energy balance measurements in the wake of a low-pressure turbine blade

The turbulent kinetic energy budget in the wake generated by a high lift, low-pressure two-dimensional blade cascade of the T106 profile was investigated experimentally using hot-wire anemometry. The purpose of this study is to examine the transport mechanism of the turbulent kinetic energy and provide validation data for turbulence modeling. Point measurements were conducted on a high spatial resolution, two-dimensional grid that allowed precise derivative calculations. Positioning of the probe was achieved using a high accuracy traversing mechanism. The turbulent kinetic energy (TKE) convection, production, viscous diffusion and turbulent diffusion were all obtained directly from experimental measurements. Dissipation and pressure diffusion were calculated indirectly using techniques presented and validated by previous investigators. Results for all terms of the turbulent kinetic energy budget are presented and discussed in detail in the present work.     

A correlation for heat transfer during subcooled boiling on a single tube with forced crossflow

Many heat exchangers, such as shell and tube heat exchangers and kettle reboilers, involve boiling with flow across tubes. For rational design of such heat exchangers, it is desirable to be able to predict heat transfer on a single tube. The dimensionless correlation presented here agrees well with available data for subcooled boiling during crossflow on a single tube. The correlating parameters are the same as those used for boiling inside tubes¹⁶. The data correlated include three fluids, four tube materials, tube diameters from 1.2 to 25.4 mm, subcooling from 0 to 80°C, and velocities from 0.02 to 7.8 m/s. The mean deviation of 334 data points is 9.5%. Hence the new correlation appears to be usable over a wide range of parameters.     

沸腾换热的分形分析

研究沸腾换热过程是安全、高效地利用能源的基础.简要评述了沸腾换热(池内沸腾、流动沸腾、临界热流密度和纳米流体沸腾换热)的研究进展;详细论述了采用分形理论和方法研究沸腾换热分析解的理论和方法;指出了采用分形理论和方法有可能解决其它尚未解决的有关沸腾换热的若干课题和方向.