Flow pattern transition instability during flow boiling in a single microchannel

We studied the unique characteristics of flow boiling in a single microchannel, including the periodic pressure drop, mass flow rate, and temperature fluctuations, in terms of a long time period. Experiments were conducted using a single horizontal microchannel and deionized water to study boiling instabilities at very small mass and heat flow rate conditions. A Polydimethylsiloxane (PDMS) rectangular single microchannel had a hydraulic diameter of 103.5 μm and a length of 40 mm. A series of piecewise serpentine platinum microheaters were fabricated on the inner bottom wall of the rectangular microchannel to supply thermal energy to the test fluid. Real-time flow visualizations of the flow pattern inside the microchannel were performed simultaneously with measurements of the experimental parameters. Tests were performed for mass fluxes of 170 and 360 kg/m² s and heat fluxes of 200–530 kW/m². The test results showed that the heated wall temperature, pressure drop, and mass flux all fluctuated with a long period and large amplitude. These periodic fluctuations exactly matched the transition of two alternating flow patterns inside the microchannel: a bubbly/slug flow and an elongated slug/semi-annular flow. Therefore, the flow pattern transition instability in the single microchannel caused a cyclic behavior of the wall temperature, pressure drop, and mass flux, and this behavior had a very long period (100–200 s) and large amplitude.     

Numerical study of bubble growth and wall heat transfer during flow boiling in a microchannel

A numerical study has been performed to analyze the wall heat transfer mechanisms during growth of a vapor bubble inside a microchannel. The microchannel is of 200 μm square cross section and a vapor bubble begins to grow at one of the walls, with liquid coming in through the channel inlet. The complete Navier–Stokes equations along with continuity and energy equations are solved using the SIMPLER method. The liquid vapor interface is captured using the level set technique. Experiments have been conducted to validate the numerical model. The bubble growth rate and shapes show good agreement between numerical and experimental results. The numerical results show that the wall heat transfer increases with wall superheat but stays almost unaffected by the liquid flow rate. The liquid vapor surface tension value has little influence on bubble growth and wall heat transfer. However, the bubble with the lowest contact angle resulted in the highest wall heat transfer.     

Numerical study of condensation flow regime transition in wavy microchannels

A numerical research on flow regime transition in wavy microchannels was conducted. The model was based on the volume of fluid approach and user-defined routines including interfacial mass transfer and latent heat. The observed droplet flow, annular–wavy flow, injection flow, and slug–bubbly flow were qualitatively compared against experimental data and transition lines were established. The effects of inlet vapor velocity, wall heat flux, and microchannel geometry characteristics on the annular length, occurrence frequency of injection flow, initial slug volume, and bubble detachment frequency were investigated.     

Numerical simulation of thermofluid characteristics of two-phase slug flow in microchannels

A fundamental study of heat transfer characteristics of two-phase slug flow in microchannels is carried out employing the Volume-of-Fluid (VOF) method. Despite of the fact that numerical simulations of two-phase flows in microchannels have been attempted by many investigators, most efforts seem to have failed in correctly capturing the flow physics, especially those pertaining to the slug flow regime characteristics. The presence of a thin liquid film in the order of 10 μm around the bubble is a contributing factor to the above difficulty. Typically, liquid films have a significant effect on the flow field and heat transfer characteristics. In the simulations reported in this paper, the film is successfully captured and a very high local convective heat transfer coefficient is observed in the film region. A strong coupling between the conductive heat transfer in the solid wall and the convective heat transfer in the flow field is observed and characterized. Results showed that unsteady heat transfer through the solid wall in the axial direction is comparable to that in the radial direction. Results also showed that a fully developed condition could be achieved fairly quickly compared to single-phase flows. The fully developed condition is defined based on the Peclet number (Pe) and a dimensionless length of the liquid slug. Local and time-averaged Nusselt numbers for slug flows are reported for the first time. It was found that significant improvements in the heat transfer coefficient could be achieved by short slugs where the Nusselt number was found to be 610% higher than in single-phase flows. The study revealed new findings related to slug flow heat transfer in microchannels with constant wall heat flux.     

Local heat transfer distribution and effect of instabilities during flow boiling in a silicon microchannel heat sink

Flow boiling of the perfluorinated dielectric fluid FC-77 in a silicon microchannel heat sink is investigated. The heat sink contains 60 parallel microchannels each of 100 μm width and 389 μm depth. Twenty-five evenly distributed temperature sensors in the substrate yield local heat transfer coefficients. The pressure drop across the channels is also measured. Experiments are conducted at five flow rates through the heat sink in the range of 20–80 ml/min with the inlet subcooling held at 26 K in all the tests. At each flow rate, the uniform heat input to the substrate is increased in steps so that the fluid experiences flow regimes from single-phase liquid flow to the occurrence of critical heat flux (CHF). In the upstream region of the channels, the flow develops from single-phase liquid flow at low heat fluxes to pulsating two-phase flow at high heat fluxes during flow instability that commences at a threshold heat flux in the range of 30.5–62.3 W/cm² depending on the flow rate. In the downstream region, progressive flow patterns from bubbly flow, slug flow, elongated bubbles or annular flow, alternating wispy-annular and churn flow, and wall dryout at highest heat fluxes are observed. As a result, the heat transfer coefficients in the downstream region experience substantial variations over the entire heat flux range, based on which five distinct boiling regimes are identified. In contrast, the heat transfer coefficient midway along the channels remains relatively constant over the heat flux range tested. Due to changes in flow patterns during flow instability, the heat transfer is enhanced both in the downstream region (prior to extended wall dryout) and in the upstream region. A previous study by the authors found no effect of instabilities during flow boiling in a heat sink with larger microchannels (each 300 μm wide and 389 μm deep); it appears therefore that the effect of instabilities on heat transfer is amplified in smaller-sized channels. While CHF increases with increasing flow rate, the pressure drop across the channels has only a minimal dependence on flow rate once boiling is initiated in the microchannels, and varies almost linearly with increasing heat flux.     

Boiling heat transfer and critical heat flux of ethanol–water mixtures flowing through a diverging microchannel with artificial cavities

This paper presents an experimental study on the convective boiling heat transfer and the critical heat flux (CHF) of ethanol–water mixtures in a diverging microchannel with artificial cavities. The results show that the boiling heat transfer and the CHF are significantly influenced by the molar fraction (x_m) as well as the mass flux. For the single-phase convection region except for the region near the onset of nucleate boiling with temperature overshoot, the single-phase heat transfer coefficient is independent of the wall superheat and increases with a decrease in the molar fraction. After boiling incipience, the two-phase heat transfer coefficient is much higher than that of single-phase convection. The two-phase heat transfer coefficient shows a maximum in the region of bubbly-elongated slug flow and deceases with a further increase in the wall superheat until approaching a condition of CHF, indicating that the heat transfer is mainly dominated by convective boiling. A flow-pattern-based empirical correlation for the two-phase heat transfer coefficient of the flow boiling of ethanol–water mixtures is developed. The overall mean absolute error of the proposed correlation is 15.5%, and more than 82.5% of the experimental data were predicted within a ±25% error band. The CHF increases from x_m = 0–0.1, and then decreases rapidly from x_m = 0.1–1 at a given mass flux of 175 kg/m² s. The maximum CHF is reached at x_m = 0.1 due to the Marangoni effect, indicating that small additions of ethanol into water could significantly increase the CHF. On the other hand, the CHF increases with increasing the mass flux at a given molar fraction of 0.1. Moreover, the experimental CHF results are compared with existing CHF correlations of flow boiling of the mixtures in a microchannel.     

Numerical simulation of the pressure drop and heat transfer of two phase slug flows in microtubes using moving frame of reference technique

Enhancement in heat transfer using two phase slug flow in microtubes and microchannels has encouraged researchers to focus on this topic as one of the potential methods for miniaturizing heat sinks and exchangers. Numerical simulation of two phase slug flow is time consuming, so some researchers conduct their numerical studies using the moving frame of reference technique for a unit cell consisting of only one slug, i.e. a single phase study, in order to accelerate the simulation process. Both single phase and two phase simulation methods have been performed in the present study and results have been compared. This shows to what extent the moving frame of reference assumption is valid in the case of two phase flow in hydrophilic microtubes when a thin liquid film exists around moving gas bubbles. The present comparison has been conducted for pressure drop and heat transfer for two thermal boundary conditions i.e. constant wall temperature and constant wall heat flux. It has been shown that the moving frame of reference method overpredicts pressure drop and heat transfer and possible reasons have been discussed. This also shows that in a slug flow with no film around bubbles more heat transfer could be achieved.     

Analysis on the heat transfer in a square cavity with a semitransparent wall: Effect of the roof materials

A theoretical study on conjugated heat transfer (natural convection, radiation and conduction) in a square cavity with turbulent flow is presented. The cavity is a representation of a room, where the left wall is isothermal, the right wall is semitransparent (glass), the lower wall is considered as insulated and on the upper opaque wall heat conduction is present. Both conductive walls (opaque and semitransparent) interact with the ambient. The semitransparent wall is subject to a constant heat flux (G₂ = 736 W/m²) whereas on the opaque wall a constant heat flux (G₁ = 875 W/m²) falls perpendicularly. The sizes of the cavity under study were 5.0, 4.0, 3.0 and 2.0 m. The upper opaque wall was considered as a mixture of concrete and a composite material (concrete–expanded polystyrene) with different thicknesses and diverse types of water-repellent coatings on top of it. From the results, it was found that the white coating on top of the opaque wall significantly reduces the amount of energy towards the inside of the cavity. It was also determined that the opaque wall with a 20 cm thickness shows the best thermal performance and it is the most adequate to reduce thermal gains inside the cavity. Correlations for the total heat transfer as a function of the cavity size, the type of coating and material of the opaque upper wall are proposed.     

Heat transfer characteristics of flow boiling in a single horizontal microchannel

Heat transfer characteristics of subcooled flow boiling of FC-72 in a single horizontal circular cross-section microchannel (480 μm i.d., 800 μm o.d., 102 mm long) are presented. Different flow patterns, both in the stable and unstable flow boiling regimes, have been captured using high speed video camera. Data in small, medium, high and very high heat flux cases under small, medium and high mass flux has been presented. Convective heat transfer coefficients in each flow boiling situation have been calculated and presented. Stable flow boiling with alternating bubbly/slug flow, slug/annular flow and annular/mist flow have been observed for heat flux of 150 kW/m² or higher and mass flux of 1500 kg/m² s or higher. Back and forth oscillations with flow instabilities have been observed in cases of lower heat and mass fluxes. However, no complete reverse flow in upstream direction has been observed.     

Subcooled flow boiling and microbubble emission boiling phenomena in a partially heated microchannel

A simultaneous visualization and measurement study has been carried out to investigate subcooled flow boiling and microbubble emission boiling (MEB) phenomena of deionized water in a partially heated Pyrex glass microchannel, having a hydraulic diameter of 155 μm, which was integrated with a Platinum microheater. Effects of mass flux, inlet water subcooling and surface condition of the microheater on subcooled flow boiling in microchannels are investigated. It is found that MEB occurred at high inlet subcoolings and at high heat fluxes, where vapor bubbles collapsed into microbubbles after contacting with the surrounding highly subcooled liquid. In the fully-developed MEB regime where the entire microheater was covered by MEB, the mass flux, the inlet water subcooling and the heater surface condition have only small effects on the boiling curves. The occurrence of MEB in microchannel can remove a large amount of heat flux, as high as 14.41 MW/m² at a mass flux of 883.8 kg/m² s, with only a moderate rise in wall temperature. Therefore, MEB is a very promising method for cooling of microelectronic chips. Heat transfer in the fully-developed MEB in the microchannel is presented, which is compared with existing subcooled flow boiling heat transfer correlations for macrochannels.     

The effect of axial conduction on heat transfer in a liquid microchannel flow

Analysis is presented for conjugate heat transfer in a parallel-plate microchannel. Axial conduction in the fluid and in the adjacent wall are included. The fluid is a constant property liquid with a fully-developed velocity distribution. The microchannel is heated by a uniform heat flux applied to the outside of the channel wall. The analytic solution is given in the form of integrals by the method of Green’s functions. Quadrature is used to obtain numerical results for the local and average Nusselt number for various flow velocities, heating lengths, wall thicknesses, and wall conductivities. These results have application in the optimal design of small-scale heat transfer devices in areas such as biomedical devices, electronic cooling, and advanced fuel cells.     

Hydrodynamics and heat transfer during flow boiling instabilities in a single microchannel

Boiling in microchannels is widely considered as one of the front runners in process intensification heat removal. Flow boiling heat transfer in microchannel geometry and the associated flow instabilities are not well understood, further research is necessary into the flow instabilities adverse effect on heat transfer.Boiling is induced in microchannel geometry (hydraulic diameter 727 μm) to investigate several flow instabilities. A transparent, metallic, conductive deposit has been developed on the exterior of rectangular microchannels, allowing simultaneous heating and visualisation.Presented in this paper is data for a particular case with a uniform heat flux of 4.26 kW/m² applied to the microchannel and inlet liquid mass flowrate, held constant at 1.13 × 10^?5 kg/s. In conjunction with obtaining high-speed images, a sensitive infrared camera is used to record the temperature profiles on the exterior wall of the microchannel, and a data acquisition system is used to record the pressure fluctuations over time. Various phenomena are apparent during the flow instabilities; these can be characterised into timescales occurring at 100’s seconds, 10’s seconds, several seconds and finally milliseconds. Correlation of pressure oscillations with temperature fluctuations as a function of the heat flux applied to the microchannel is possible.From analysis of our results, images and video sequences with the corresponding physical data obtained, it is possible to follow simultaneously particular flow, pressure and temperature conditions leading to nucleate boiling, flow instabilities and transition regimes during flow boiling in a microchannel. The investigation allowed us to quantify and characterise the timescales of various observed instabilities during flow boiling in a microchannel. High speed imaging revealed some of the controlling physical mechanisms responsible for the observed instabilities.     

A study on the simulation and experiment of a microchannel counter-flow heat exchanger

Numerical simulations and experimental tests were carried out to study the fluid flow and heat transfer characteristics for a rectangular-shaped microchannel heat exchanger. Moreover, influences of gravity to heat transfer and pressure drop behaviors of the microchannel heat exchanger were presented by variation of the physical inclinations of the microchannel heat exchanger system used for experiments. For experimental results, a heat flux of 17.4 W/cm² was achieved for the heat exchanger. Besides, the results obtained for the actual effectiveness and for the effectiveness (the so-called effectiveness-NTU method) were determined. In this study, the pressure drop decreases as the water temperature rises. As the pressure drop increases from 880 to 4400 Pa, the mass flow rate increases from 0.1812 to 0.8540 g/s. In addition, the results obtained from numerical analyses were in good agreement with those obtained from experiments, with discrepancies of the heat transfer coefficient estimated to be less than 9%.     

Flow Boiling Heat Transfer of Refrigerant R-134a in Copper Microchannel Heat Sink

In this paper we present experimental data on heat transfer and pressure drop characteristics at flow boiling of refrigerant R-134a in a horizontal microchannel heat sink. The primary objective of this study was to experimentally establish how the local heat transfer coefficient and pressure drop correlate with the heat flux, mass flux, and vapor quality. The copper microchannel heat sink contains 21 microchannels with 335 × 930 μm² cross section. The microchannel plate and heating block were divided by the partition wall for the local heat flux measurements. Distribution of local heat transfer coefficients along the length and width of the microchannel plate was measured in the range of external heat fluxes from 50 to 500 kW/m²; the mass flux varied within 200–600 kg/m²-s, and pressure varied within 6–16 bar. The obvious impact of heat flux on the magnitude of heat transfer coefficient was observed. It showed that nucleate boiling is the dominant mechanism for heat transfer. A new model of flow boiling heat transfer, considering nucleate boiling suppression and liquid film evaporation, was proposed and verified experimentally in this paper.     

A theoretical study of evaporative heat transfer in high‐velocity two‐phase flow of air–water in a small vertical tube

A theoretical study was performed to investigate the evaporative heat transfer of high‐velocity two‐phase flow of air–water in a small vertical tube under both heating conditions of constant wall temperature and constant heat flux. A simplified two‐phase flow boundary layer model was used to evaluate the evaporative heat transfer characteristics of the annular two‐phase flow. The analytical results show that the gravitational force, the gas–liquid surface tension force, and the inertial force are much smaller than the frictional force and hence can be neglected for a small tube. The evaporative heat transfer characteristics of the small tube with constant wall temperature are quite close to those of the small tube with constant heat flux. The mechanism of the heat transfer enhancement is the forced convective evaporation on the surface of the thin liquid film. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(5): 430–444, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10110     

Effects of superheat and temperature-dependent thermophysical properties on evaporating thin liquid films in microchannels

A model based on the augmented Young–Laplace equation and the Clausius–Clapeyron equation was developed to describe the extended evaporating meniscus in a microchannel. The effects of the adsorbed film thickness, channel height and temperature-dependent thermophysical properties of the fluid are included in the model at wall superheats up to 50 K. The liquid flow is coupled with the vapor flow to obtain the mass transport across the liquid–vapor interface. The results show that the constant thermophysical property model greatly overestimates the liquid pressure difference and the total thin film heat transfer rate at higher superheats compared with the variable thermophysical property model. The adsorbed film thickness, which is controlled by the disjoining pressure limit, reaches a minimum near about 20 K superheat for water. The maximum film curvature and liquid pressure difference then decrease at superheats larger than 20 K. The effects of the capillary pressure limit produced by the channel height can be reduced by increasing the superheat.     

On the dynamics and heat transfer of bubble train in micro-channel flow boiling

The dynamics and heat transfer characteristics of flow boiling bubble train moving in a micro channel is studied numerically. The coupled level set and volume of fluid (CLSVOF) is utilized to track interface and a non-equilibrium phase change model is applied to calculate the interface temperature as well as heat flux jump. The working fluid is R134a and the wall material is aluminum. The fluid enters the channel with a constant mass flux (335 kg/m² 1 s), and the boundary wall is heated with constant heat flux (14 kW/m²). The growth of bubbles and the transition of flow regime are compared to an experimental visualization. Moreover, the bubble evaporation rate and wall heat transfer coefficient have been examined, respectively. Local heat transfer is significantly enhanced by evaporation occurring vicinity of interface of the bubbles. The local wall temperature is found to be dependent on the thickness of the liquid film between the bubble train and the wall.     

An experimental study of flow boiling instability in a single microchannel

A simultaneous visualization and measurement study has been carried out to investigate stable and unstable flow boiling phenomena of deionized water in a single microchannel having a hydraulic diameter of 155 µm with a bottom Pyrex glass wall. Fifteen platinum serpentine microheaters, bonded on the Pyrex glass wall, were used to measure local instantaneous wall temperatures. At low mass flux, a syringe pump was used to drive the subcooled water passing through the microchannel. Stable and unstable flow boiling modes in the single microchannel are identified, and flow pattern maps in terms of heat flux and mass flux as well as in term of exit vapor quality are presented respectively. It was found that unstable flow boiling occurred in the single microchannel if the exit vapor quality x_e > 0.013.     

Microchannel size effects on local flow boiling heat transfer to a dielectric fluid

Heat transfer with liquid–vapor phase change in microchannels can support very high heat fluxes for use in applications such as the thermal management of high-performance electronics. However, the effects of channel cross-sectional dimensions on the two-phase heat transfer coefficient and pressure drop have not been investigated extensively. In the present work, experiments are conducted to investigate the local flow boiling heat transfer of a dielectric fluid, Fluorinert FC-77, in microchannel heat sinks. Experiments are performed for mass fluxes ranging from 250 to 1600 kg/m² s. Seven different test pieces made from silicon and consisting of parallel microchannels with nominal widths ranging from 100 to 5850 μm, all with a nominal depth of 400 μm, are considered. An array of temperature sensors on the substrate allows for resolution of local temperatures and heat transfer coefficients. The results of this study show that for microchannels of width 400 μm and greater, the heat transfer coefficients corresponding to a fixed wall heat flux as well as the boiling curves are independent of channel size. Also, heat transfer coefficients and boiling curves are independent of mass flux in the nucleate boiling region for a fixed channel size, but are affected by mass flux as convective boiling dominates. A strong dependence of pressure drop on both channel size and mass flux is observed. The experimental results are compared to predictions from a number of existing correlations for both pool boiling and flow boiling heat transfer.     

Fluid flow and convective heat transfer in flat microchannels

This study investigates the design, construction and instrumentation of an experimental microchannel, with a rectangular cross-section and large aspect ratio, that allows characterization of the flow and convective heat transfer under well defined and precise conditions and makes it possible to vary the hydraulic diameter of the microchannel. The flow friction coefficient is estimated by direct pressure drop measurements inside the microchannel in a zone where the flow is fully developed. Since the wall thermal conditions inside the microchannel can not be measured directly, their estimation requires temperature measurements in the wall thickness and an inverse heat conduction method. The thermal and hydrodynamic results obtained by varying the hydraulic diameter between 1 mm and 100 μm do not deviate from the theory or empirical correlations for large-scale channels. These results let us confirm that for smooth walls the continuum mechanics laws for convection and fluid mechanics remain valid in microchannels of hydraulic diameter greater than or equal to 100 μm.