Prandtl and capillary effects on heat transfer performance within laminar liquid–gas slug flows

This paper investigates a two-phase non-boiling slug flow regime for the purposes of enhancing heat transfer in microchannel heat sinks or compact heat exchangers. The primary focus is upon understanding the mechanisms leading to enhanced heat transfer and the effect of using different Prandtl number fluids, leading to variations in Capillary number. Experimental work was conducted using Infrared thermography and results are presented in the form of Graetz solution, spanning both the thermal entrance and fully developed flow regions. Nusselt numbers enhancements were observed throughout when data was reduced to account for void fraction. However, the gaseous void was also noted to demonstrate an artificial increase with greater thicknesses of the liquid film, due to higher Capillary numbers. Up to 600% enhancement in heat transfer rates were observed over conventional Poiseuille flow. This was verified through Nusselt number measurements over inverse Graetz number ranges from 10^?4 to 1 and slug length to channel diameter ratios from 0.88 to 32. Varying Prandtl and Capillary numbers caused notable effects in the transition region between entrance and fully developed flows. Significant Nu oscillations were observed for low Pr fluids due to internal circulation within the slug. However, these oscillations are observed to be damped out when higher Prandtl number fluids are employed. The thickness of the liquid film surrounding the gas bubbles is shown to have a significant influence on heat transfer performance. Overall, this study provides a greater understanding of the mechanisms leading to significant enhancements in heat exchange devices employing two-phase gas–liquid flows without boiling.     

Flow regime-based modeling of heat transfer and pressure drop in microchannel flow boiling

Local heat transfer coefficients and pressure drops during boiling of the dielectric liquid fluorinert FC-77 in parallel microchannels were experimentally investigated in recent work by the authors. Detailed visualizations of the corresponding two-phase flow regimes were performed as a function of a wide range of operational and geometric parameters. A new transition criterion was developed for the delineation of a regime where microscale effects become important to the boiling process and a conventional, macroscale treatment becomes inadequate. A comprehensive flow regime map was developed for a wide range of channel dimensions and experimental conditions, and consisted of four distinct regions – bubbly, slug, confined annular, and alternating churn/annular/wispy-annular flow regimes. In the present work, physics-based analyses of local heat transfer in each of the four regimes of the comprehensive map are formulated. Flow regime-based models for prediction of heat transfer coefficient in slug flow and annular/wispy-annular flow are developed and compared to the experimental data. Also, a regime-based prediction of pressure drop in microchannels is presented by computing the pressure drop during each flow regime that occurs along the microchannel length. The results of this study reveal the promise of flow regime-based modeling efforts for predicting heat transfer and pressure drop in microchannel boiling.     

Experimental analysis of the unit cell approach for two-phase flow dynamics in curved flow channels

Flow behavior of gas–liquid mixtures in thin channels has become increasingly important as a result of miniaturization of fluid and thermal systems. The present empirical study investigates the use of the unit cell or periodic boundary approach commonly used in two-phase flows. This work examines the flow patterns formed in small tube diameter (<3 mm) and curved geometry flow systems for air–water mixtures at standard conditions. Liquid and gas superficial velocities were varied from 0.1 to 7.0 (～±0.01) m/s and 0.03 to 14 (～±0.2) m/s for air and water respectively to determine the flow pattern formed in three geometries and dispersed bubble, plug, slug and annular flow patterns are reported using high-frame rate videography. Flow patterns formed were plotted on the generalized two-phase flow pattern map to interpret the effect of channel size and curvature on the flow regime boundaries. Relative to a straight a channel, it is shown that a ‘C shaped’ channel that causes a directional change in the flow induces chaotic advection and increases phase interaction to enhance gas bubble or liquid slug break-up thus altering the boundaries between the dispersed bubble and plug/slug flow regimes as well as between the annular and plug/slug flow regimes.     

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.     

Characteristics of Two-Phase Flows in a Rectangular Microchannel with a T-Junction Type Gas-Liquid Mixer

In this study, gas–liquid two-phase flows in a horizontal rectangular microchannel have been investigated. The rectangular microchannel has a hydraulic diameter of 0.235 mm, and a width and depth of 0.24 mm and 0.23 mm, respectively. A T-junction-type gas–liquid mixer was used to introduce gas and liquid in the channel. In order to know the effects of liquid properties, distilled water, ethanol, and HFE7200 were used as the test liquids, with nitrogen gas was used as the test gas. The flow pattern, the bubble length, the liquid slug length, and the bubble velocity in two-phase flow were measured with a high-speed video camera, and the void fraction was determined from the bubble velocity data and the superficial gas velocity data. In addition, the pressure drop was also measured with a calibrated differential pressure transducer. The bubble length data were compared with the calculation by the scaling law proposed by Garstecki et al. [7]. The bubble velocity data and/or the void fraction data were well correlated with the well-known drift flux model [12] with a new distribution parameter correlation developed in this study. The frictional pressure drop data were also well correlated with the Lockhart-Martinelli method with a correlation of the two-phase friction multiplier.     

Unsteady steam condensation flow patterns inside a miniature tube

Unsteady steam condensation inside a single miniature tube has been studied. The visualization of different instantaneous and periodically two-phase flow is conducted for different experimental conditions. The two-phase flow characterization is obtained using the image processing. Annular, slug bubbly, spherical bubbly, and wavy flows are observed by varying the steam inlet pressure and cooling heat transfer. The cycle of the periodically flows are compared. It is shown that increasing the cooling heat flow rate reduces the number of the instabilities and the injected bubbles. The axial vapor velocity decreases during the waves growth. The local distribution of the condensate film thickness is analyzed. It is shown that the liquid film becomes thinner near the meniscus-like interface because of the surface tension effect. The reverse annular flow is observed at the end of each periodic flow when the bubbles leave the channel. It can be concluded from experimental results that the stratification effect is not significant during the condensation inside the miniature tube. The capillary pressure evolution is measured. The maximum values are obtained in the waves locations and near the meniscus of the annular flow.     

Experimental investigations of flow boiling heat transfer and pressure drop in straight and expanding microchannels – A comparative study

Flow boiling experiments were conducted in straight and expanding microchannels with similar dimensions and operating conditions. Deionized water was used as the coolant. The test vehicles were made from copper with a footprint area of 25 mm × 25 mm. Microchannels having nominal width of 300 μm and a nominal aspect ratio of 4 were formed by wire cut Electro Discharge Machining process. The measured surface roughness (Ra) was about 2.0 μm. To facilitate easier comparison with the straight microchannels and also to simplify the method of fabrication, the expanding channels were formed with the removal of fins at selected location from the straight microchannel design, instead of using a diverging channel. Tests were performed on both the microchannels over a range of mass fluxes, heat fluxes and an inlet temperature of 90 °C. It was observed that the two-phase pressure drop across the expanding microchannel heat sink was significantly lower as compared to its straight counterpart. The pressure drop and wall temperature fluctuations were seen reduced in the expanding microchannel heat sink. It was also noted that the expanding microchannel heat sink had a better heat transfer performance than the straight microchannel heat sink, under similar operating conditions. This phenomenon in expanding microchannel heat sink, which was observed in spite of it having a lower convective heat transfer area, is explained based on its improved flow boiling stability that reduces the pressure drop oscillations, temperature oscillations and hence partial dry out.     

Optimum thermal design of microchannel heat sinks by the simulated annealing method

By adopting the simulated annealing method, a three-dimensional numerical simulation is executed to minimize the thermal resistance of the microchannel heat sink corresponding to the optimum specification under the fixed flow power. The depths of the microchannel heat sink in this study are fixed at either 1 cm or 2 cm. Based on the theory of the fully developed flow, the pressure drop between the inlet and exit in each single channel can be analytically derived if the flow power and the associated specification of the microchannel heat sink are fixed in advance. Then, this pressure drop will be used as the input condition to calculate the temperature distribution of the microchannel heat sink. For the first part of the optimum analysis, the fin width, and channel width are chosen as the design variables to find their optimum sizes. As to the second part of the present analysis, three design variables including channel height, fin width and channel width are individually prescribed as a suitable range to search for their optimum geometric configuration when the other specifications of the microchannel heat sink are fixed as 24 different cases.     

Numerical Analysis of Flow and Thermal Performance of Liquid-Cooling Microchannel Heat Sinks with Bifurcation

Single-phase liquid-cooling microchannels have received great attention to remove the gradually increased heat loads of heat sinks. Proper changes of the flow path and/or heat transfer surface can result in much better thermal performance of microchannel heat sinks. In this study, a kind of rectangular straight microchannel heat sink with bifurcation flow arrangement has been designed, and the corresponding laminar flow and heat transfer have been investigated numerically. Four different configurations are considered. The effects of the bifurcation ratio (the initial channel number over the bifurcating channel number) and length ratio (the channel length before bifurcation over the bifurcation channel length) on laminar heat transfer, pressure drop, and thermal resistance are considered and compared with those of the traditional straight microchannel heat sink without bifurcation flow. The overall thermal resistances subjected to inlet Reynolds number and pumping power are compared for the five microchannel heat sinks. Results show that the thermal performance of the microchannel heat sink with bifurcation flow is better than that of the corresponding straight microchannel heat sink. The heat sinks with larger bifurcation ratio and length ratio provide much better thermal performance. It is suggested to employ bifurcation flow path in the liquid-cooling microchannel heat sinks to improve the overall thermal performance by proper design of the bifurcation position and number of channels.     

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.     

Unsteady-state fluctuations analysis during bubble growth in a “rectangular” microchannel

Boiling in microchannels shows great potential for cooling systems and compact heat removal applications. However for confidence in this cooling technique, it is essential that any excursions from typical flow boiling are understood and predicted. Confined bubble growth can cause pressure fluctuations which interfere with bubble nucleation and growth and can also lead to flow reversal and instances of temperature excursions. Boiling experiments are performed in a single rectangular microchannel of hydraulic diameter 771 μm, using n-Pentane as the working fluid. A heating technique was incorporated on the exterior walls of the microchannel; a transparent, metallic, conductive deposit, which allows simultaneous uniform heating and visualisation to be achieved. In conjunction with obtaining high-speed imaging, an infrared camera is used to record the temperature profile at the microchannel wall, and sensitive pressure sensors are used to record the pressure drop across the microchannel over time. During flow boiling in the microchannel periodic and non-periodic fluctuations in both the channel pressure drop and channel temperature profile over time are apparent. In this paper we provide a full analysis of the temperature measurements and pressure data obtained during the growth of a vapour bubble in the microchannel. An augmentation of the heat transfer coefficient of over 216% has been achieved during periodic two-phase flow boiling in the microchannel. However overpressure (over 410% increase) in the microchannel occurs at corresponding instances to the heat transfer enhancement. The two time steps during the periodic bubble dynamics, namely the bubble expansion time period and the waiting time period in-between the bubble expansion fluctuations, are also investigated and modelled. It was determined that both the bubble dynamics and the channel wall heating time period are responsible for the pressure and temperature fluctuation time periods observed.     

ADIABATIC GAS–LIQUID FLOW IN MICROCHANNELS

This article presents a review of adiabatic two-phase flow in minichannels and microchannels. Differences between them are identified and explained based on this review and our own research. Several channels of decreasing diameter were used in our experiments to determine the effect of the channel size on the two-phase flow of nitrogen gas and water. The effect of channel geometry was examined by characterizing the two-phase flow in a circular and square microchannel of similar size. Only slug flow was observed in the microchannels. Four new sub-classes of slug flow were subsequently defined. A new correlation was developed for the time-averaged void fraction data in the microchannels. The two-phase pressure drop in microchannels was predicted by treating the two phases as being separate with a large velocity difference. Regarding the effect of microchannel geometry, the transition boundaries on the two-phase flow regime maps were shifted for the slug flow subcategories.     

Flow Pattern and Heat Transfer Behavior of Boiling Two—Phase flow in Inclined Pipes

Movable Electrical Conducting Probe (MECP), a kind of simple and reliable measuring transducer, used for predicting full-flow-path flow pattern in a boiling vapor/liquid two-phase flow is introduced in this paper. When the test pipe is set at different inclination angles, several kinds of flow patterns, such as bubble, slug, churn, intermittent, and annular flows, may be observed in accordance with the locations of MECP. By means of flow pattern analysis, flow field numerical calculations have been carried out, and heat transfer coefficient correlations along full-flow-path derived. The results show that heat transfer performance of boiling two-phase flow could be significantly augmented as expected in some flow pattern zones.The results of the investigation, measuring techniques and conclusions contained in this paper would be a useful reference in foundational research for prediction of flow pattern and heat transfer behavior in boiling two-phase flow, as well as for turbine vane liquid-cooling design.     

Flow pattern transition instability in a microchannel with CO2 bubbles produced by chemical reactions

The present study investigates experimentally the two-phase flow in a rectangular microchannel with CO₂ bubbles generated by chemical reactions of sulfuric acid (H₂SO₄) and sodium bicarbonate (NaHCO₃). The microchannel with a hydraulic diameter of 132.7 μm is prepared using bulk micromachining and anodic bonding process. Evolution of two-phase flow patterns in the microchannel was observed using a high speed video camera and the corresponding pressure drop was investigated. It is found that the inlet concentration and flow rate of reactants have a significant effect on the evolution of two-phase flow characteristics and slug flow is the dominant flow pattern. The flow pattern transition instability between bubbly-slug and slug flow takes place for the cases with highest inlet concentration, i.e., C = 0.8 mol/L, and low flow rates of this study. The oscillation frequency is from 0.024 to 0.041 Hz and the magnitude of oscillation in pressure drop is from 10 to 15 kPa. A mechanism based on the flow circulation in the liquid slug may reasonably explain the flow pattern transition instability. Small amplitude, high frequency oscillations with a frequency of about 45 Hz are superimposed on the low frequency flow pattern transition as well as prevail for other cases without the flow pattern transition instability. The two-phase flow pressure drop increases with increase in both flow rate and inlet concentration.     

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.     

Transient flow patterns and bubble slug lengths in parallel microchannels with oxygen gas bubbles produced by catalytic chemical reactions

Transient flow patterns and bubble slug lengths were investigated with oxygen gas (O₂) bubbles produced by catalytic chemical reactions using a high speed camera bonded with a microscope. The microreactor consists of an inlet liquid plenum, nine parallel rectangular microchannels followed by a micronozzle, using the MEMS fabrication technique. The etched surface was deposited by the thin platinum film, which is acted as the catalyst. Experiments were performed with the inlet mass concentration of the hydrogen peroxide from 50% to 90% and the pressure drop across the silicon chip from 2.5 to 20.0 kPa. The silicon chip is directly exposed in the environment thus the heat released via the catalytic chemical reactions is dissipated into the environment and the experiment was performed at the room temperature level. It is found that the two-phase flow with the catalytic chemical reactions display the cyclic behavior. A full cycle consists of a short fresh liquid refilling stage, a liquid decomposition stage followed by the bubble slug flow stage. At the beginning of the bubble slug flow stage, the liquid slug number reaches maximum, while at the end of the bubble slug flow stage the liquid slugs are quickly flushed out of the microchannels. Two or three large bubbles are observed in the inlet liquid plenum, affecting the two-phase distributions in microchannels. The bubble slug lengths, cycle periods as well as the mass flow rates are analyzed with different mass concentrations of hydrogen peroxide and pressure drops. The bubble slug length is helpful for the selection of the future microreactor length ensuring the complete hydrogen peroxide decomposition. Future studies on the temperature effect on the transient two-phase flow with chemical reactions are recommended.     

Adiabatic two-phase flow in rectangular microchannels with different aspect ratios: Part II – bubble behaviors and pressure drop in single bubble

One of the major flow patterns in a microchannel is an elongated bubble flow, which is similar to a long slug bubble. Behaviors and pressure drop for a single bubble in a rectangular microchannel were studied. Based on the experiments in Part I of this paper, data for liquid superficial velocities of 0.06–0.8 m/s, gas superficial velocities of 0.06–0.66 m/s and AR of 0.92, 0.67, 0.47 and 0.16 were analyzed. The velocity, length, number, and frequency of the single bubble in the rectangular microchannel were obtained from image processing based on a unit cell model. The bubble velocities were proportional to total superficial velocity. As the aspect ratio decreased, the portion of the bubble area increased due to the corner effect. New correlation of the bubble velocity for different aspect ratio was proposed. Also, bubble and liquid slug length, the number of the unit cell and bubble frequency were analyzed with different aspect ratios. The pressure drop for the single bubble in the rectangular microchannels was evaluated using the information of the bubble behavior. The pressure drop in the single elongated bubble was proportional to the bubble velocity. The pressure drop in the single elongated bubble in the rectangular microchannel increased as the aspect ratio decreased.     

Two-Phase Flow Patterns and Heat Transfer in Parallel Microchannels

MicroChannel heat sinks with two-phase flow can satisfy the increasing heat removal requirements of modern micro electronic devices. One of the important aspects associated with two- phase flows in microchannels is to study the bubble behavior. However, in the literature most of the reports present data of only a single channel. This does not account for flow mixing and hydrodynamic instability that occurs in parallel microchannels, connected by common inlet and outlet collectors. In the present study, experiments were performed for air- water and steam- water flow in parallel triangular microchannels with a base of 200-300μ m. The experimental study is based on systematic measurements of temperature and flow pattern by infrared radiometry and high-speed digital video imaging. In air-water flow, different flow patterns were observed simultaneously in the various microchannels at a fixed values of water and gas flow rates. In steam-water flow, instability in uniformly heated microchannels was observed.     

Two-Phase Flow Patterns in Microchannel Vaporization of CO2 at Near-Critical Pressure

Carbon dioxide (CO ₂ /R-744) is receiving renewed interest as a refrigerant, in many cases for systems with microchannel heat exchangers that have high pressure capability, efficient heat transfer, and compact design. A good understanding of two-phase flow of evaporating CO ₂ in microchannels is needed to analyze and predict heat transfer. A special test rig was built in order to observe two-phase flow patterns using a horizontal glass tube with ID 0.98 mm. Flow visualization experiments were conducted for temperatures 20°C and 0°C and for mass flux ranging from 100 to 580 kg m^?2 s^?1 . The observations showed a dominance of intermittent (slug) flow at low x and wavy annular flow with entrainment of droplets at higher x. The aggravated dryout problem reported from heat transfer experiments at high mass flux could be explained by increased entrainment. The flow pattern observations did not fit generalized maps or transition lines showed in the literature.