Simulation of flow boiling in micro-channels: Effects of inlet flow rate and hot-spots

This paper presents the findings of a numerical study on the flow boiling in a micro-channel heat sink. The Navier-Stokes equations, energy equation, and the continuity equation are solved in a finite-volume framework using the front-tracking method. The numerical method is validated by comparison with the experimental results for a slug bubble growth, and vertical flow boiling. The numerical method is then used to study the effect of changing the inflow mass-velocity on the heat transfer coefficient, bubble size distribution, and the bubble nucleation frequency for a constant heat flux. The mean heat transfer coefficient of all the cases is found to be nearly twice that of the single-phase heat transfer coefficient. The bubble nucleation frequency is found to increase monotonically with the inflow mass-velocity. The bubble size distribution along the channel is found to become flatter as the mass-velocity is increased. We identify three distinct phases of the bubble evolution, namely the initial rapid growth phase, the boiling dominant phase, and finally the condensation dominant phase. Subsequently, the numerical method is used to study the effect of having a hot-spot near the bubble nucleation site on the heat transfer characteristics. It is found that the bubble nucleation frequency increases and the bubbles’ maximum volume decreases as the intensity of the hot-spot is increased for a fixed inlet flow rate. It is also observed that the average heat transfer coefficient does not change significantly with changing the intensity of the hot-spot, and that the bubble size distribution along the channel becomes flatter as the intensity of the hot-spot is increased.     

二维水平通道内流动沸腾换热的格子Boltzmann模拟

利用格子Boltzmann方法模拟二维水平通道内水的流动沸腾过程,获得不同壁面过热度下流型特点和不同因素对换热过程的影响规律。结果表明,随着壁面过热度升高,流道内流型依次经历从泡状流、弹状流到反环流的转变,平均热流密度和平均换热系数先增大后减小。入口流速降低会使流道内出现受限气泡流,核态沸腾受到抑制。提高入口流速能够有效促进气泡脱离,壁面平均换热系数随入口流速增大而增大,但增长速率有所减小。减小通道宽度有利于汽化现象发生,核态沸腾得到强化,壁面平均换热系数有所提高。     

Study on heat transfer performance of spray cooling: model and analysis

In consideration of droplet–film impaction, film formation, film motion, bubble boiling (both wall nucleation bubbles and secondary nucleation bubbles), droplet–bubble interaction, bulk air convection and radiation, a model to predict the heat and mass transfer in spray cooling was presented in this paper. The droplet–film impaction was modeled based on an empirical correlation related with droplet Weber number. The film formation, film motion, bubble growth, and bubble motion were modeled based on dynamics fundamentals. The model was validated by the experimental results provided in this paper, and a favorable comparison was demonstrated with a deviation below 10%. The film thickness, film velocity, and non-uniform surface temperature distribution were obtained numerically, and then analyzed. A parameters sensitivity analysis was made to obtain the influence of spray angle, surface heat flux density, and spray flow rate on the surface temperature distribution, respectively. It can be concluded that the heat transfer induced by droplet–film impaction and film-surface convection is dominant in spray cooling under conditions that the heated surface is not superheated. However, the effect of boiling bubbles increases rapidly while the heated surface becomes superheated.     

Studies on heat transfer and flow characteristics in subcooled flow boiling—Part 1. Boiling characteristics

Measurement of wall temperature profile and photographic observation are performed for R-113 subcooled boiling flow in a channel with heat fluxes up to the CHF. The incipient boiling superheats measured are little affected by mass velocity and liquid subcooling. Hysteresis in boiling observed by increasing and decreasing heat flux seems to be ascribed to variation in size of active nucleation cavities on the wall. Increasing heat flux up to the CHF, the bubble density on the heated surface increases and remarkably large coalescent bubbles appear periodically near the heating section outlet.     

Bubble nucleation characteristics in pool boiling of a wetting liquid on smooth and rough surfaces

Quantitative measurements are obtained from high-speed visualizations of pool boiling at atmospheric pressure from smooth and roughened surfaces, using a perfluorinated hydrocarbon (FC-77) as the working fluid. The boiling surfaces are fabricated from aluminum and prepared by mechanical polishing in the case of the smooth surface, and by electrical discharge machining (EDM) in the case of the roughened surface. The roughness values (R_a) are 0.03 and 5.89 μm for the polished and roughened surfaces, respectively. The bubble diameter at departure, bubble departure frequency, active nucleation site density, and bubble terminal velocity are measured from the monochrome movies, which have been recorded at 8000 frames per second with a digital CCD camera and magnifying lens. Results are compared to predictions from existing models of bubble nucleation behavior in the literature. Wall superheat, heat flux, and heat transfer coefficient are also reported.     

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

Enhancement of pool boiling heat transfer in a horizontal water layer through surface roughness and screen coverage

An experimental investigation was carried out to study the augmentation of heat transfer in saturated pool boiling of a liquid water layer on a heated horizontal stainless steel plate by roughing the surface and/or covering it with a single layer of stainless steel screen. The results were presented in terms of the boiling curves. Effects of various parameters – the surface roughness, liquid level and size of the stainless steel mesh on the boiling heat transfer were examined in detail. The measured data clearly indicated that a lowering of the liquid level from 60 to 5?mm in water depth causes heat transfer reduction. Roughing the surface was found to sig- nificantly enhance the heat transfer. Use a layer of metal screen to cover the heated surface was shown to substantially augment the heat transfer especially for a shallow water layer if the mesh size is comparable with the bubble departure diameter. Covering the rough surface with the metal mesh, however, reduced the heat transfer.     

Flow boiling behaviors in hydrophilic and hydrophobic microchannels

Surface wettability is a critical parameter in small scale phenomena, especially two-phase flow, since the surface force becomes dominant as size decreases. In present study, experiments of water flow boiling in hydrophilic and hydrophobic rectangular microchannels were conducted to investigate the wettability effect on flow boiling in rectangular microchannels. The rectangular microchannels were fabricated with a photosensitive glass to visualize flow pattern. The hydrophilic bare photosensitive glass microchannel was chemically treated to obtain a hydrophobic microchannel. And, visualization of flow patterns was carried out. And boiling heat transfer and two-phase pressure drop was analyzed with visualization results. The boiling heat transfer coefficient in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with nucleation site density and liquid film motion. And the pressure drop in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with unstable motions of bubble and liquid film. Finally, we find out the wettability is important parameter on the flow pattern, which were highly related with two-phase heat and mass transfer.     

Effects of bubble size on heat transfer enhancement by sub-millimeter bubbles for laminar natural convection along a vertical plate

Injection of sub-millimeter bubbles is considered a promising technique for enhancing natural convection heat transfer for liquids. So far, we have experimentally investigated heat transfer characteristics of laminar natural convection flows with sub-millimeter bubbles. However, the effects of the bubble size on the heat transfer have not yet been understood. The purpose of this study is to clarify the effects of the bubble size on the heat transfer enhancement for the laminar natural convection of water along a vertical heated plate with uniform heat flux. Temperature and velocity measurements, in which thermocouples and a particle tracking velocimetry technique are, respectively used, are conducted to investigate heat transfer and flow characteristics for different bubble sizes. Moreover, two-dimensional numerical simulations are performed to comprehensively understand the effects of bubble injection on the flow near the heated plate. The result shows that the ratio of the heat transfer coefficient with sub-millimeter-bubble injection to that without injection ranges from 1.3 to 2.2. The result also shows that for a constant bubble flow rate, the heat transfer coefficient ratio increases with a decrease in the mean bubble diameter. It is expected from our estimation based on both experimental data and simulation results that this increase results from an increase in the advection effect due to bubbles.     

Artificial boiling heat transfer in the free convection to carbonic acid solution

Free convection phenomenon has been experimentally investigated around a horizontal rod heater in carbonic acid solution. Because of the tendency of the solution to desorb carbon dioxide gas when temperature is increased, bubbles appear when cylinder surface is heated. The bubbles consists mainly carbon dioxide and also a negligible amount of water vapor. The results present that dissolved carbon dioxide in water significantly enhances the heat transfer coefficient in compare to pure free convection regime. This is mainly due to the microscale mixing on the heat transfer surface, which is induced by bubble formation. In this investigation, experiments are performed at different bulk temperatures between 288 and 333 K and heat fluxes up to 400 kW m⁻² at atmospheric pressure. Bubble departure diameter, nucleation site density and heat transfer coefficient have been experimentally measured. A model has been proposed to predict the heat transfer coefficient.     

Effect of acoustic cavitation on boiling heat transfer

Boiling heat transfer on a horizontal circular copper tube in an acoustical field is investigated experimentally and the relation between the liquid cavitation, the boiling and the micro bubble radii are analyzed theoretically. The results show that cavitation bubbles have an important influence on the nucleation, growth and collapse of vapor embryo within cavities on the heat transfer surface and that the enhancement of boiling heat transfer by acoustic cavitation mainly depends on whether the vapor embryo is activated by the cavitation bubbles to initiate boiling.     

Experimental investigation of flow boiling characteristics in a cross-linked microchannel heat sink

This paper experimentally investigates flow boiling characteristics in a cross-linked microchannel heat sink at low mass fluxes and high heat fluxes. The heat sink consists of 45 straight microchannels each with a hydraulic diameter of 248 μm and heated length of 16 mm. Three cross-links, of width 500 μm, are introduced in the present microchannel heat sink to achieve better temperature uniformity and to avoid flow mal-distribution. Flow visualization, flow instability, two-phase pressure drop, and two-phase heat transfer measurements are conducted using the dielectric coolant FC-72 over a range of heat flux from 7.2 to 104.2 kW/m², mass flux from 99 to 290 kg/m² s, and exit quality from 0.01 to 0.71. Thermochromic liquid crystals are used in the present study as full-field surface temperature sensors to map the temperature distribution on the heat sink surface. Flow visualization studies indicate that the observed flow regime is primarily slug. Visual observations of flow patterns in the cross-links demonstrate that bubbles nucleate and grow rapidly on the surface of the cross-links and in the tangential direction at the microchannels’ entrance due to the effect of circulations generated in those regions. The two-phase pressure drop strongly increases with the exit quality, at x_e,o < 0.3, and the two-phase frictional pressure drop increases by a factor of 1.6–2 compared to the straight microchannel heat sink. The flow boiling heat transfer coefficient increases with increasing exit quality at a constant mass flux, which is caused by the dominance of the nucleation boiling mechanism in the cross-link region.     

Study of transient heat transfer and synchronized flow visualizations during sub-cooled flow boiling in a small aspect ratio microchannel

The challenges that microchannel flow boiling technology faces are the lack of understanding of underlying mechanisms of heat transfer during various flow boiling regimes and a dearth of analytical models that can predict heat transfer. This paper aims to understand flow boiling heat transfer mechanisms by analyzing results obtained by synchronously captured high-speed flow visualizations with local, transient temperature data. Using Inverse Heat Conduction Problem (IHCP) solution methodology, the transient wetted surface heat flux and temperature as well as heat transfer coefficient are calculated. These are then correlated with the visual data. Experiments are performed on a single microchannel embedded with fast response temperature sensors located (630 µm) below the wetted surface. The height, width and length of the microchannel are 0.42 mm, 2.54 mm and 25.4 mm respectively. De-ionized, de-gassed water is used as the working fluid. Two heat fluxes are tested at each of the mass fluxes of 182 kg/(m²s) and 380 kg/(m²s). Because of vapor confinement, slug flow is observed for the tested conditions. The present study provides detailed insights into the effect of various events such as passage of vapor slug, 3-phase contact line, partial-dry-out and liquid slug on transient heat transfer coefficient. Transient heat transfer coefficient peaks when thin film evaporation mechanism is prevalent. The peak value is influenced by the distance of bubble incipience as well as downstream events obstructing the flow. Heat transfer coefficient during the passage of liquid slug and 3-phase contact line were relatively lower for the tested experimental conditions.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

Bubble growth in saturated pool boiling in water and surfactant solution

The presence of surfactant additives in water was found to enhance the boiling heat transfer significantly. The objective of the present investigation is to compare the bubble growth in water to that of a surfactant solution with negligible environmental impact. The study was conducted at two values of heat fluxes to clarify the effect of the heat flux on the dynamics of bubble nucleation. The bubble growth under condition of pool boiling in water and non-ionic surfactant solution was studied using high speed video technique. The bubble generation was studied on a horizontal flat surface; and the natural roughness of the surface was used to produce the bubbles.     

Local heat transfer from a hot plate to a water jet

Jet impingement boiling is very efficient in cooling of hot surfaces as a part of the impinging liquid evaporates. Because of its importance to many cooling procedures, investigations on basic mechanisms of jet impingement boiling heat transfer are needed. Until now, most of the experimental studies, carried out under steady-state conditions, used a heat flux controlled system and were limited by the critical heat flux (CHF). The present study focuses on steady-state experiments along the entire boiling curve for hot plate temperatures of up to 700°C. A test section has been built up simulating a hot plate. It is divided into 8 independently heated modules of 10 mm length to enable local heat transfer measurements. By means of temperature controlled systems for each module local steady-state experiments in the whole range between single phase heat transfer and film boiling are possible. By solving the two dimensional inverse heat conduction problem, the local heat flux and the corresponding wall temperature on the surface of each module can be computed. The measurements show important differences between boiling curves measured at the stagnation line and those obtained in the parallel flow region. At the stagnation line, the transition boiling regime is characterised by very high heat fluxes, extended to large wall superheats. Inversely, boiling curves in the parallel flow region are very near to classical ones obtained for forced convection boiling. The analysis of temperature fluctuations measured at a depth of 0.8 mm from the boiling surface enables some conclusions on the boiling mechanism in the different boiling regimes.     

Effects of channel dimension,heat flux,and mass flux on flow boiling regimes in microchannels

Experiments are conducted with a perfluorinated dielectric fluid, Fluorinert FC-77, to investigate the effects of channel size and mass flux (225–1420 kg/m²s) on microchannel flow boiling regimes by means of high-speed photography. Seven different silicon test pieces with parallel microchannels of widths ranging from 100 to 5850 μm, all with a depth of 400 μm, are considered. Flow visualizations are performed with a high-speed digital video camera while local measurements of the heat transfer coefficient are simultaneously obtained. The visualizations and the heat transfer data show that flow regimes in the microchannels of width 400 μm and larger are similar, with nucleate boiling being dominant in these channels over a wide range of heat flux. In contrast, flow regimes in the smaller microchannels are different and bubble nucleation at the walls is suppressed at a relatively low heat flux for these sizes. Two types of flow regime maps are developed and the effects of channel width on the flow regime transitions are discussed.     

Influence of surface topography in the boiling mechanisms

The present paper addresses the qualitative and quantitative analysis of the pool boiling heat transfer over micro-structured surfaces. The surfaces are made from silicon chips, in the context of pool boiling heat transfer enhancement of immersion liquid cooling schemes for electronic components. The first part of the analysis deals with the effect of the liquid properties. Then the effect of surface micro-structuring is discussed, covering different configurations, from cavities to pillars being the latter used to infer on the potential profit of a fin-like configuration. The use of rough surfaces to enhance pool boiling mainly stands on the arguments that the surface roughness will increase the liquid–solid contact area, thus enhancing the convection heat transfer coefficient and will promote the generation of nucleation sites. However, one should not disregard bubble dynamics. Indeed, the results show a strong effect of bubble dynamics and particularly of the interaction mechanisms in the overall cooling performance of the pair liquid–surface. The inaccurate control of these mechanisms leads to the formation of large bubbles and strong vertical and horizontal coalescence effects promote the very fast formation of a vapor blanket, which causes a steep decrease of the heat transfer coefficient. This effect can be strong enough to prevail over the benefit of increasing the contact area by roughening the surface. For the micro-patterns used in the present work, the results evidence that one can reasonably determine guiding pattern characteristics to evaluate the intensity of the interaction mechanisms and take out the most of the patterning to enhance pool boiling heat transfer, when using micro-cavities. Instead, it is far more difficult to control the appearance of active nucleation sites and the optimization of the patterns allowing a reasonable control of the interaction mechanisms and in particular of horizontal coalescence, when dealing with the patterns based on micro-pillars. Hence, providing an increase of the liquid contact area by an effective increase of the roughness ratio is not enough to assure a good performance of the micro-structured surface. Despite it was not possible to clearly evidence a pin–fin effect or of an additional cooling effect due to liquid circulation between the pillars, the results show a significant increase of the heat transfer coefficient of about 10 times for water and 8 times for the dielectric fluid, in comparison to the smooth surface, when the micro-patterning based on pillars is used.