Numerical Simulation of Bubbles Deformation,Flow, and Coalescence in a Microchannel Under Pseudo-Nucleation Conditions

This paper reports the results of numerical study on bubbles deformation, flow, and coalescence under pseudo-nucleate boiling conditions in horizontal mini-/microchannels. The numerical simulation, which is based on the multiphase model of volume of fluid method, aims to study the corresponding flow behaviors of nucleate bubbles generated from the tube walls in mini-/microchannels so as to understand the effect of confined surfaces/walls on nucleate bubbles and heat transfer. Under the pseudo- or quasi-nucleate boiling condition, superheated small vapor bubbles are injected at the wall to ensure that the bubbles generation is under a similar condition of real nucleation. The numerical study examined the fluid mechanics of bubble motion with heat transfer, but the mass transfer across the bubble–liquid interface is not simulated in the present work.     

Contribution of thin-film evaporation during flow boiling inside microchannels

Flow boiling through microchannels is characterized by nucleation and growth of vapor bubbles that fill the entire channel cross-sectional area. As the bubbles nucleate and grow inside the microchannel, a thin film of liquid or a microlayer gets trapped between the bubbles and the channel walls. The heat transfer mechanism present at the channel walls during flow boiling is studied numerically. It is then compared to the heat transfer mechanisms present during nucleate pool boiling and in a moving evaporating meniscus. Increasing contact angle improved wall heat transfer in case of nucleate boiling and moving evaporating meniscus but not in the case of flow boiling inside a microchannel. It is shown that the thermal and the flow fields present inside the microchannel around a bubble are fundamentally different as compared to nucleate pool boiling or in a moving evaporating meniscus. It is explained why thin-film evaporation is the dominant heat transfer mechanism and is responsible for creating an apparent nucleate boiling effect inside a microchannel.     

Experimental and theoretical studies on subcooled flow boiling of pure liquids and multicomponent mixtures

To improve the design of modern industrial reboilers, accurate knowledge of boiling heat transfer coefficients is essential. In this study flow boiling heat transfer coefficients for binary and ternary mixtures of acetone, isopropanol and water were measured over a wide range of heat flux, subcooling, flow velocity and composition. The measurements cover the regimes of convective heat transfer, transitional boiling and fully developed subcooled flow boiling. Two models are presented for the prediction of flow boiling heat transfer coefficients. The first model is the combination of the Chen model with the Gorenflo correlation and the Schlünder model for single and multicomponent boiling, respectively. This model predicts flow boiling heat transfer coefficients with acceptable accuracy, but fails to predict the nucleate boiling fraction NBF reasonably well. The second model is based on the asymptotic addition of forced convective and nucleate boiling heat transfer coefficients. The benefit of this model is a further improvement in the accuracy of flow boiling heat transfer coefficient over the Chen type model, simplicity and the more realistic prediction of the nucleate boiling fraction NBF.     

The effect of heated wall thickness and materials on nucleate boiling at high heat flux

The present work is to numerically investigate the effect of heater side factors on the nucleate boiling at high heat flux, which is characterized by the existence of macrolayer. Two-region equations are proposed to study both thermo-capillary driven flow in the liquid layer and heat conduction in the solid wall. The numerical results indicate that the thermo-capillary driven flow in the macrolayer and evaporation at the vapor-liquid interface constitute a very efficient heat transfer mechanism to explain the high heat transfer coefficient of nucleate boiling heat transfer near CHF. For a very thin wall and/or wall with a poor thermal conductivity (heat side factors) are found to have significant effect on flow pattern in the liquid layer and the temperature distribution in the heated wall.     

Single-phase and two-phase heat transfer characteristics of low temperature hybrid micro-channel/micro-jet impingement cooling module

This study examines the single-phase and two-phase cooling performance of a hybrid micro-channel/micro-jet impingement cooling scheme using HFE 7100 as working fluid. This scheme consists of supplying coolant from a series of jets that deposit liquid into the micro-channels. A single-phase numerical scheme that utilizes the k–ε turbulent model and a method for determining the extent of the laminarized wall layer shows very good predictions of measured wall temperatures. It is shown jet velocity has a profound influence on single-phase cooling performance. High jet velocities enable jet fluid to penetrate the axial micro-channel flow and produce a strong impingement effect at the wall. On the other hand, the influence of jets at low jet velocities is greatly compromised compared to the micro-channel flow. During nucleate boiling, vapor layer development along the micro-channel in the hybrid module is fundamentally different from that encountered in conventional micro-channels. Here, subcooled jet fluid produces repeated regions of bubble growth followed by bubble collapse, rather than the continuous growth common to conventional micro-channel flow. By reducing void fraction along the micro-channel, the hybrid scheme contributes greater wall temperature uniformity. Increasing subcooling and/or flow rate delay the onset of boiling to higher heat fluxes and higher wall temperatures, and also increase critical heat flux considerably. A nucleate boiling heat transfer coefficient correlation is developed that fits the present data with a mean absolute error of 6.10%.     

Modeling and numerical investigation of hydrogen nucleate flow boiling heat transfer

Liquid hydrogen flow boiling heat transfer in tubes is of great importance in the hydrogen applications such as superconductor cooling, hydrogen fueling. In the present study, a numerical model for hydrogen nucleate flow boiling based on the wall partition heat flux model is established. The key parameters in the model such as active nucleation site density, bubble departure diameter and frequency are carefully discussed and determined to facilitate the modeling and simulation of hydrogen flow boiling. Simulation results of the numerical model show reasonably well agreement with experimental data from different research groups in a wide operation condition range with the means absolute error (MAE) of 10.6% for saturated and 5.3% for subcooled flow boiling. Based on the model, wall heat flux components and void fraction distribution of hydrogen flow boiling are studied. Effects of mass flow rate and wall heat flux on the flow boiling heat transfer performance are investigated. It is found that in the hydrogen nucleate flow boiling, the predominated factor is the Boiling number, rather than the vapor quality. A new simple correlation is proposed for predicting hydrogen saturated nucleate flow boiling Nusselt number. The MAE between the correlation predicted and experimentally measured Nusselt number is 13.6% for circular tubes and 12.5% for rectangular tubes. The new correlation is applicable in the range of channel diameter 4–6.35 mm, Reynolds number 64000–660,000, saturation temperature 22–29 K, Boiling number 8.37 × 10^?5–2.33 × 10^?3.     

Numerical Simulation of Single Bubble Dynamics During Flow Boiling Conditions on a Horizontal Surface

Complete three-dimensional numerical simulations of single bubble dynamics during flow boiling conditions are carried out using the computational fluid dynamics code FLOW3D based on the volume-of-fluid method. The analyses include a numerically robust kinetic phase-change model and transient wall heat conduction. The simulation approach is calibrated by comparison with available experimental and theoretical data. It is found that the observed hydrodynamics (i.e., bubble shape, departure, and deformation) are simulated very well. The comparison with high-resolution transient temperature measurements during a heating foil experiment indicates that the modeling of the spatiotemporal heat sink distribution during bubble growth requires major attention. The simulation tool is employed for single bubble dynamics during flow boiling on a horizontal heating wall, and the agreement is excellent with published experimental data. The numerical results indicate how bulk flow velocity and wall heat transfer influence the bubble dynamics and heat transfer characteristics.     

Saturated flow boiling heat transfer of refrigerant R-410A in a horizontal annular finned duct

An experiment is conducted here to investigate the saturated flow boiling heat transfer characteristics of ozone friendly refrigerant R-410A in a horizontal annular finned duct. Meanwhile the associated bubble characteristics in the duct are also inspected from the flow visualization. The experimental data are presented in terms of saturated flow boiling curves, boiling heat transfer coefficients and flow photos. In addition, empirical correlation equations for the saturated flow boiling heat transfer coefficient and mean bubble departure diameter are proposed. The saturated flow boiling curves show that boiling hysteresis is insignificant in the flow and the wall superheat needed for the onset of nucleate boiling is slightly affected by the refrigerant mass flux. Besides, the boiling curves are mainly affected by the imposed heat flux and refrigerant mass flux. Moreover, the measured saturated flow boiling heat transfer coefficient increases with the imposed heat flux and refrigerant mass flux. Furthermore, at a higher refrigerant mass flux the departing bubbles are smaller.     

Simulation on nucleate boiling in micro-channel

Using the VOF multiphase flow model, numerical simulations are conducted to investigate the nucleate boiling of water in micro-channels. The Marangoni heat transfer through the bubble surface is analyzed, and is compared with the incipient heat flux at the onset of nucleate boiling in micro-channels. The bubble growth in the channel is divided into two stages. At the initial stage, bubble growth is controlled by surface tension, while at the second stage the incipient heat transfer dominated the boiling process. In the results, the full process of bubble generating, growing, departing, combining, and shrinking in the channel is displayed. The simulated results with similar condition are agreed well with some experimental results in references. The method and discussion in the paper are helpful to the investigation of the mechanism of micro-scale two-phase flow and heat transfer.     

High resolution measurements of wall temperature distribution underneath a single vapour bubble under low gravity conditions

Heat transfer performance in nucleate boiling crucially depends on a circular thin film area near the foot of a vapour bubble where high heat fluxes and thus high local evaporation rates occur. The corresponding wall temperature drop close to the tiny thin film area is computed using an existing nucleate boiling model. To verify the predicted temperature distribution, an experiment is designed with a thin electrically heated wall featuring two-dimensional, high resolution temperature measurement using unencapsulated thermochromic liquid crystals. By means of temperature measurements during a parabolic flight under low-g conditions, the validity of the model used to calculate the temperature distribution in the tiny thin film area could be confirmed.     

Fundamental Issues Related to Flow Boiling and Two-Phase Flow Patterns in Microchannels – Experimental Challenges and Opportunities

ABSTRACT

Flow boiling heat transfer in microchannels is used today in many diverse applications. The previous studies addressing the effect of channel size, heat flux, vapor quality, and mass flux on heat transfer during flow boiling are reviewed in the present paper. The relationship between flow characteristics and flow boiling heat transfer was studied experimentally for refrigerant R-C318 at moderate reduced pressures where the contribution of nucleate boiling is decisive. Flow boiling mechanisms were identified using an annular microchannel with transparent outer wall for successive visualization of boiling. The considerable suppression of nucleate boiling heat transfer was observed at transition to annular flow and explained by formation of a liquid flow with thin film and dry spots. A general equation for prediction of two-phase flow boiling heat transfer inside the circular, annular, and rectangular microchannels is proposed and verified using the experimental data. This equation accounts for the nucleate boiling suppression, forced convection, and thin film evaporative heat transfer in the form that allows to distinguish more clearly the contribution of each mechanism of heat transfer under the conditions, when it is predominant. A new approach for prediction of transition to the annular flow is proposed and verified, using the experimental data.     

Time periodic flow boiling heat transfer of R-134a and associated bubble characteristics in a narrow annular duct due to flow rate oscillation

An experiment is conducted here to investigate the effects of the imposed time periodic refrigerant flow rate oscillation in the form of nearly a triangular wave on refrigeriant R-134a flow boiling heat transfer and associated bubble characteristics in a horizontal narrow annular duct with the duct gap fixed at 2.0 mm. The results indicate that when the imposed heat flux is close to that for the onset of stable flow boiling, intermittent flow boiling appears in which nucleate boiling on the heated surface does not exist in an entire periodic cycle. At somewhat higher heat flux persistent boiling prevails. Besides, the refrigerant flow rate oscillation only slightly affects the time-average boiling curves and heat transfer coefficients. Moreover, the heated wall temperature, bubble departure diameter and frequency, and active nucleation site density are found to oscillate periodically in time as well and at the same frequency as the imposed mass flux oscillation. Furthermore, in the persistent boiling the resulting heated wall temperature oscillation is stronger for a longer period and a larger amplitude of the mass flux oscillation. And for a larger amplitude of the mass flux oscillation, stronger temporal oscillations in the bubble characteristics are noted. The effects of the mass flux oscillation on the size of the departing bubble and active nucleation site density dominate over the bubble departure frequency, causing the heated wall temperature to decrease and heat transfer coefficient to increase at reducing mass flux in the flow boiling, opposing to that in the single-phase flow. But they are only mildly affected by the period of the mass flux oscillation. However, a short time lag in the wall temperature oscillation is also noted. Finally, a flow regime map is provided to delineate the boundaries separating different boiling regimes for the R-134a flow boiling in the annular duct.     

On the transition from partial to fully developed subcooled flow boiling

The subject of the present study is to relate the boiling heat transfer process with experimentally observed bubble behaviour during subcooled flow boiling of water in a vertical heated annulus. It presents an attempt to explain the transition from partial to fully developed flow boiling with regard to bubble growth rates and to the time that individual bubbles spend attached to the heater surface.Within the partial nucleate boiling region bubbles barely change in size and shape while sliding a long distance on the heater surface. Such behaviour indicates an important contribution of the microlayer evaporation mechanism in the overall heat transfer rate. With increasing heat flux, or reducing flow rate at constant heat flux, bubble growth rates increase significantly. Bubbles grow while sliding, detach from the heater, and subsequently collapse in the bulk fluid within a distance of 1-2 diameters parallel to the heater surface. This confirms that bubble agitation becomes a leading heat transfer mode with increasing heat flux. There is however, a sharp transition between the two observed bubble behaviours that can be taken as the transition from partial to fully developed boiling. Hence, this information is used to develop a new model for the transition from partial to fully developed subcooled flow boiling.     

A Study of Bubble Behavior and Boiling Heat Transfer Enhancement under Electric Field

The effect of a D.C. electric field on nucleate boiling heat transfer for refrigerants, R11, R113, and FC72, was investigated experimentally in a single-tube shell/tube heat exchanger by using the temperature control method of wall superheat. Also the behavior of bubble under nonuniform electric field produced by wire electrodes was studied by numerical calculation. For R11, the electrohydrodynamic (EHD) enhancement for boiling heat transfer was observed for all ranges of wall superheat tested. However, the enhancement in boiling heat transfer disappeared if the wall superheat exceeded 13°C for R113, and no electric field effect on the boiling heat transfer was observed for FC72. An application of approximately 5 kV was enough to eliminate the boiling hysteresis for R11 and R113. Numerical study of the electric field in a single medium has hinted that the bubbles are forced away from the heating surface and toward the electrostatic stagnation point by the dielectrophoretic force. Such modified bubble motion turns out to promote the boiling heat transfer if one uses proper electrode configuration.     

Effects of heater size and working fluids on nucleate boiling heat transfer

The present study is an experimental investigation of nucleate boiling heat transfer mechanism in pool boiling from wire heaters immersed in saturated FC-72 coolant and water. The vapor volume flow rate departing from a wire during nucleate boiling was determined by measuring the volume of bubbles from the wire utilizing the consecutive-photo method. The effects of the wire size on heat transfer mechanism during a nucleate boiling were investigated, varying 25 μm, 75 μm, and 390 μm, by measuring vapor volume flow rate and the frequency of bubbles departing from a wire immersed in saturated FC-72. One wire diameter of 390 μm was selected and tested in saturated water to investigate the fluid effect on the nucleate boiling heat transfer mechanism. Results of the study showed that an increase in nucleate boiling heat transfer coefficients with reductions in wire diameter was related to the decreased latent heat contribution. The latent heat contribution of boiling heat transfer for the water test was found to be higher than that of FC-72. The frequency of departing bubbles was correlated as a function of bubble diameters.     

Nucleate Pool Boiling of Pure Liquids and Binary Mixtures：PartI－Analytical Model for Boiling Heat Transfer of Pure Liquids on Smooth Tubes

ＮｕｃｌｅａｔｅＰｏｏｌＢｏｉｌｉｎｇｏｆＰｕｒｅＬｉｑｕｉｄｓａｎｄＢｉｎａｒｙＭｉｘｔｕｒｅｓ：ＰａｒｔＩ－ＡｎａｌｙｔｉｃａｌＭｏｄｅｌｆｏｒＢｏｉｌｉｎｇＨｅａｔＴｒａｎｓｆｅｒｏｆＰｕｒｅＬｉｑｕｉｄｓｏｎＳｍｏｏｔｈＴｕｂｅｓＧｕｏｑｉｎ...     

Nucleate Pool Boiling of Pure Liquids and Binary Mixtures ：Part I—Analytical Model for Boiling Heat Transfer of Pure Liquids on Smooth Tubes

A mechanism is proposed for nucleate pool boiling heat transfer along with a general model for both pure liquids and binary mixtures. A combined physical model of bubble growth is also proposed along with a corresponding bubble growth model for pure liquids on smooth tubes. Using the general model and the bubble growth model for pure liquids, an analytical model for nucleate pool boiling heat transfer of pure liquids on smooth tubes is developed.     

Evaporation of refrigerants in a horizontal tube: an improved flow pattern dependent heat transfer model compared to ammonia data

Few experimental test data are available for evaporation of ammonia inside tubes and numerous new data have been measured and presented here. An improved approach to the prediction of flow boiling heat transfer in horizontal tubes has been proposed through the study of each flow pattern separately, incorporating a new criterion defining the onset of nucleate boiling as a function of the critical convective heat transfer coefficient representative of the location where nucleate boiling might occur. A new function, based on a pseudo-Biot number delineates two different mean heat fluxes on the perimeter of the tube in stratified types of flow, one in contact with the liquid and one in contact with the vapor. Considering pure convective heat transfer, or mixed convective and nucleate heat transfer, this division allows the use of a common criterion to be applied to each flow pattern. Even if the database showed that the flow conditions in the annular liquid film were close to, or in the turbulent to laminar flow transition, and even if the major part of the experimental points where purposely obtained close to the various flow pattern transitions, the new model showed very good agreement with the experimental database of refrigerants HFC-134a and ammonia. Due to the precision of the new flow pattern map and the effectiveness of the onset on nucleate boiling criterion, this new heat transfer model accurately predicts the heat transfer conditions during evaporation.     

Models for enhanced boiling heat transfer by unusual Marangoni effects under microgravity conditions

A new approach is suggested to enhance boiling heat transfer through introduction of unusual surface tension effects. The surface tension of aqueous solutions of alcohols with a chain length longer than four carbon atoms offers a positive gradient with temperature when the temperature exceeds a certain value. Moreover, the positive gradient near the boiling point has a very large value. This will generate a considerable driving force for bubble departure. As a result, in the nucleate boiling of these solutions, the Marangoni effect around the bubble surface will not impede the bubble detachment from the heater surface but rather will provide an additional driving force for the bubble departure. This effect combines with the buoyancy under normal gravity and acts as a main driving force of bubble departure in microgravity. Models for predictions of the bubble detachment diameters, the nucleate boiling heat transfer coefficient, and the critical heat flux are developed.     

Numerical Investigation of Saturated Flow Boiling in Microchannels by the Lattice Boltzmann Method

Bubble formation in saturated flow boiling in 2D microchannels, generated from a microheater under constant wall heat flux or constant wall temperature conditions, is studied numerically based on a newly developed lattice Boltzmann model for liquid-vapor phase change. Simulations are carried out to study effects of inlet velocity, contact angle, and heater size on saturated flow boiling of water under constant wall heat flux conditions. Important information, such as effects of static contact angle on nucleation time and nucleation temperature, which was unable to be obtained by other numerical simulation methods, is obtained. Furthermore, effects of inlet velocity, contact angle, and superheat on nucleate boiling heat transfer in steady flow boiling of water under constant wall temperature conditions are also presented. It is found that the nucleate boiling heat transfer at the microheater is higher if the heater surface is more hydrophilic, because the superheated vapor at the hydrophilic wall has a thinner thermal boundary layer and a larger thermal conductivity.