Periodic flow boiling in a non-uniformly heated microchannel heat sink

Hydrodynamic and thermal characteristics of flow boiling in a non-uniformly heated microchannel were studied. Experiments were performed with a single microchannel and a series of microheaters to study the microscale boiling of water under axially non-uniform heat input conditions. A simultaneous real time visualization of the flow pattern was performed with the measurement of experimental parameters. Tests were performed over a mass flux of 309.8 kg/m² s, and heat flux of 200–600 kW/m². Test results showed different fluctuations of heated wall temperature, pressure drop, and mass flux with variations of the heat input along the flow direction. The unique periodic flow boiling in a single microchannel was observed at all heat flux conditions except for the increasing heat input distribution case which is the nearly uniform effective heat input distribution condition. The instability is correlated with flow pattern transition. For the nearly uniform effective heating condition, no fluctuation of the wall temperature, pressure drop, or mass flux was observed. We can relieve the instability by increasing total heat input along the flow direction and predict the instability using the transition criteria and flow pattern map.     

Numerical simulation of compressible gas flow and heat transfer in a microchannel surrounded by solid media

In this study, a numerical model is developed to investigate the coupled compressible gas flow and heat transfer in a microchannel surrounded by solid media. To accommodate the varying flow cross-section, the compressible gas flow model is established in a non-orthogonal curvilinear coordinate system. An iterative numerical procedure is employed to solve the coupled heat transfer and gas flow equations. The computer code for the compressible gas flow is first validated against two test problems, and then extended by including the heat conduction in the solid media. The effect of the inlet Mach number on the Nusselt number is examined. It is found that the pressure difference from the pyrolysis front to the heated surface is induced essentially by the gas addition from the channel wall, instead from the pyrolysis front. The necessity of accounting for the gas compressibility is clearly demonstrated when severe heating is applied. The pressure distribution obtained along the channel axial direction is useful for further structural analysis of composite materials.     

A separated flow model for predicting two-phase pressure drop and evaporative heat transfer for vertical annular flow

A separated flow model has been developed that is applicable to vertical annular two-phase flow in the purely convective heat transfer regime. Conservation of mass, momentum, and energy are used to solve for the liquid film thickness, pressure drop, and heat transfer coefficient. Closure relationships are specified for the interfacial friction factor, liquid film eddy-viscosity, turbulent Prandtl number, and entrainment rate. Although separated flow models have been reported previously, their use has been limited, because they were tested over a limited range of flow and thermal conditions. The unique feature of this model is that it has been tested and calibrated against a vast array of two-phase pressure drop and heat transfer data, which include upflow, downflow, and microgravity flow conditions. The agreements between the measured and predicted pressure drops and heat transfer coefficients are, on average, better or comparable to the most reliable empirical correlations. This separated flow model is demonstrated to be a reliable and practical predictive tool for computing two-phase pressure drop and heat transfer rates. All of the datasets have been obtained from the open literature.     

Adiabatic and diabatic two-phase venting flow in a microchannel

The growth and advection of the vapor phase in two-phase microchannel heat exchangers increase the system pressure and cause flow instabilities. One solution is to locally vent the vapor formed by capping the microchannels with a porous, hydrophobic membrane. In this paper we visualize this venting process in a single 124 μm by 98 μm copper microchannel with a 65 μm thick, 220 nm pore diameter hydrophobic Teflon membrane wall to determine the impact of varying flow conditions on the flow structures and venting process during adiabatic and diabatic operation. We characterize liquid velocities of 0.14, 0.36 and 0.65 m/s with superficial air velocities varying from 0.3 to 8 m/s. Wavy-stratified and stratified flow dominated low liquid velocities while annular type flows dominated at the higher velocities. Gas/vapor venting can be improved by increasing the venting area, increasing the trans-membrane pressure or using thinner, high permeability membranes. Diabatic experiments with mass flux velocities of 140 and 340 kg/s/m² and exit qualities up to 20% found that stratified type flows dominate at lower mass fluxes while churn-annular flow became more prevalent at the higher mass-flux and quality. The diabatic flow regimes are believed to significantly influence the pressure-drop and heat transfer coefficient in vapor venting heat exchangers.     

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.     

Heat transfer to air–water annular flow in a horizontal pipe

Two-phase air–water flow and heat transfer in a 25 mm internal diameter horizontal pipe were investigated experimentally. The water superficial velocity varied from 24.2 m/s to 41.5 m/s and the air superficial velocity varied from 0.02 m/s to 0.09 m/s. The aim of the study was to determine the heat transfer coefficient and its connection to flow pattern and liquid film thickness. The flow patterns were visualized using a high speed video camera, and the film thickness was measured by the conductive tomography technique. The heat transfer coefficient was calculated from the temperature measurements using the infrared thermography method. It was found that the heat transfer coefficient at the bottom of the pipe is up to three times higher than that at the top, and becomes more uniform around the pipe for higher air flow-rates. Correlations on local and average Nusselt number were obtained and compared to results reported in the literature. The behavior of local heat transfer coefficient was analyzed and the role of film thickness and flow pattern was clarified.     

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included ground-based tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.     

Heat transfer and two-phase flow characteristics of refrigerants flowing under varied heat flux in a double-pipe evaporator

A mathematical model based on the annular flow pattern is developed to simulate the evaporation of refrigerants flowing under varied heat flux in a double tube evaporator. The finite difference form of governing equations of this present model is derived from the conservation of mass, energy and momentum. The experimental set-up is designed and constructed to provide the experimental data for verifying the simulation results. The test section is a 2.5 m long counterflow double tube heat exchanger with a refrigerant flowing in the inner tube and heating water flowing in the annulus. The inner tube is made from smooth horizontal copper tubing of 9.53 mm outer diameter and 7.1 mm inner diameter. The agreement of the model with the experimental data is satisfactory. The present model can be used to investigate the axial distributions of the temperature, heat transfer coefficient and pressure drop of various refrigerants. Moreover, the evaporation rate or the other relevant parameters that is difficult to measure in the experiment are predicted and presented here. The results from the present mathematical model show that the saturation pressure and temperature of refrigerant decrease along the tube due to the tube wall friction and the flow acceleration of refrigerant. The liquid heat transfer coefficient increases with the axial length due to reducing the thickness of the liquid refrigerant film. Due to increase of the liquid heat transfer coefficient, increasing wall heat flux is obtained.Finally, the evaporation rate of refrigerant increases with increasing wall heat flux.     

Carbon dioxide flow boiling in a single microchannel - Part II: Heat transfer

Flow boiling heat transfer coefficients of CO₂ have been measured in a single microchannel. Experiments were carried out in a horizontal stainless steel tube of 0.529 mm inner diameter, for three temperatures (−10, −5 and 0 °C), with the mass flux ranging from 200 to 1200 kg/m² s and the heat flux varying from 10 to 30 kW/m². The investigation covered qualities from zero to the dryout inception, i.e. pre-dryout conditions. Compared to larger microchannels and positive temperatures, a higher contribution of convective boiling was found, with a larger heat transfer coefficient than for pure nucleate boiling. Mainly two heat transfer regimes were found, depending on the boiling number (Bo). For Bo > 1.1 × 10⁻⁴, the heat transfer coefficient was highly dependent on the heat flux and moderately influenced by the quality and the mass flux. For Bo < 1.1 × 10⁻⁴, the heat transfer coefficient was hardly affected by the heat flux but strongly influenced by the quality and the mass flux. In addition, dryout results were reported. The effect of the mass flux on the dryout inception quality was found to be highly dependent on the heat flux and the saturation temperature.     

The critical role of channel cross-sectional area in microchannel flow boiling heat transfer

Experiments are conducted with a perfluorinated dielectric fluid, Fluorinert FC-77, to identify the critical geometric parameters that affect flow boiling heat transfer and flow patterns in microchannels. In recent work by the authors (Harirchian and Garimella, 2009), seven different silicon test pieces containing parallel microchannels of widths ranging from 100 to 5850 μm, all with a depth of 400 μm were tested and it was shown that for a fixed channel depth, the heat transfer coefficient was independent of channel width for microchannels of widths 400 μm and larger, with the flow regimes in these microchannels being similar; nucleate boiling was also found to be dominant over a wide range of heat fluxes. In the present study, experiments are performed with five additional microchannel test pieces with channel depths of 100 and 250 μm and widths ranging from 100 to 1000 μm. Flow visualizations are performed using a high-speed digital video camera to determine the flow regimes, with simultaneous local measurements of the heat transfer coefficient and pressure drop. The aim of the present study is to investigate as independent parameters the channel width and depth as well as the aspect ratio and cross-sectional area on boiling heat transfer in microchannels, based on an expanded database of experimental results. The flow visualizations and heat transfer results show that the channel cross-sectional area is the important governing parameter determining boiling mechanisms and heat transfer in microchannels. For channels with cross-sectional area exceeding a specific value, nucleate boiling is the dominant mechanism and the boiling heat transfer coefficient is independent of channel dimensions; below this threshold value of cross-sectional area, vapor confinement is observed in all channels at all heat fluxes, and the heat transfer rate increases as the microchannel cross-sectional area decreases before premature dryout occurs due to channel confinement.     

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

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

Slip-flow heat transfer in a microchannel with viscous dissipation

Forced convection flow in a microchannel with constant wall temperature is studied, including viscous dissipation effect. The slip-flow regime is considered by incorporating both the velocity-slip and the temperature-jump conditions at the surface. The energy equation is solved for the developing temperature field using finite integral transform. To increase β_v Kn is to increase the slip velocity at the wall surface, and hence to decrease the friction factor. Effects of the parameters β_v Kn, β, and Br on the heat transfer results are illustrated and discussed in detail. For a fixed Br, the Nusselt number may be either higher or lower than those of the continuum regime, depending on the competition between the effects of β_v Kn and β. At a given β_v Kn the variation of local Nusselt number becomes more even when β becomes larger, accompanied by a shorter thermal entrance length. The fully developed Nusselt number decreases with increasing β irrelevant to β_v Kn. The increase in Nusselt number due to viscous heating is found to be more pronounced at small β_v Kn.     

Two-phase flow distribution of air–water annular flow in a parallel flow heat exchanger

The air and water flow distribution are experimentally studied for a round header – flat tube geometry simulating a parallel flow heat exchanger. The number of branch flat tube is 30. The effects of tube outlet direction, tube protrusion depth as well as mass flux, and quality are investigated. The flow at the header inlet is identified as annular. For the downward flow configuration, the water flow distribution is significantly affected by the tube protrusion depth. For flush-mounted configuration, most of the water flows through frontal part of the header. As the protrusion depth increases, more water is forced to the rear part of the header. The effect of mass flux or quality is qualitatively the same as that of the protrusion depth. Increase of the mass flux or quality forces the water to rear part of the header. For the upward flow configuration, however, most of the water flows through rear part of the header. The protrusion depth, mass flux, or quality does not significantly alter the flow pattern. Possible explanations are provided based on the flow visualization results. Negligible difference on the water flow distribution was observed between the parallel and the reverse flow configuration.     

Numerical study of unsteady flow and heat transfer of circular tangential direction jets flowing over the inner cylinder surface in the annular chamber

In order to understand the dynamics of vortices on heat transfer, the unsteady flow field of tangential direction jets flowing in the annular chamber is numerically investigated by scale-adaptive simulation (SAS). The jet Reynolds number is 332,000 based on the jet’s diameter and inflow velocity for a specific geometric model. The analogy theory is used to obtain the convective heat transfer coefficient distribution on the hub surface. Spectral analysis via fast Fourier transform (FFT) is used to analyze frequency information that flows inside the chamber. The proper orthogonal decomposition (POD) method is performed on the velocity field in the chamber and the convective heat transfer coefficient on the hub surface using a snapshot method. The fast Fourier transform helps find the dominant frequency of the unsteady flow in the chamber. The time sequence of velocity fields on the radial plane shows the presence of cyclic flapping of the jet. The proper orthogonal decomposition analysis indicates that the unsteady periodic flow phenomenon in the chamber and unsteady heat transfer on the hub surface are mainly related to the dynamics of the counter-rotating vortices caused by the jet.     

Algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow, focusing in particular on momentum and heat transfer (so-called ‘convective boiling’) through the annular liquid film. In contrast with single-phase wall-bounded flow theory, shear-driven annular liquid films are assumed here to behave as fluid-bounded flows, mostly interacting with the shearing gas-entrained droplets core flow. Besides providing velocity and temperature profiles through the liquid film, the turbulence model proposed here predicts key parameters such as the average liquid film thickness, the void fraction and the convective boiling heat transfer coefficient with accuracies comparable or better than those of leading design correlations. This turbulence model is part of a unified annular flow modeling suite that includes methods to predict the entrained liquid fraction and the axial frictional pressure gradient. The underlying heat transfer database covers nine fluids (water, two hydrocarbons and six refrigerants) for vertical and horizontal tubes of 1.03-14.4 mm diameter and pressures of 0.1-7.2 MPa. Importantly, this study shows that there appears to be no macro-to-microscale transition when it comes to annular flow. Simply better physical modeling is required to span this range.     

Heat transfer in an oscillating vertical annular liquid column open to atmosphere

The heat transfer from a surface heated with constant heat flux to an oscillating vertical annular liquid column having an interface with the atmosphere is investigated experimentally in the present paper. The analysis is carried out for the case of different oscillation frequencies while the displacement amplitude remains constant. Based on the experimental data a correlation equation is obtained for the cycle-averaged Nusselt number as a function of kinetic Reynolds number.     

Two-phase flow modelling: the closure issue for a two-layer flow

The fully-developed, laminar flow of two fluid layers in a horizontal channel is studied by means of the height-averaged balance equations. The closure issue is addressed and the closure relations for the wall and interfacial shear stresses are given for some particular cases.     

Studying disturbance waves in vertical annular flow with high-speed video

In many annular two-phase gas–liquid flows, large disturbance waves propagate liquid mass. These waves are important for modeling of gas-to-liquid momentum transfer and liquid film behavior. High-speed videos of vertical upflow have been analyzed to extract individual and average wave data. Two types of structures, coherent waves and piece waves, have been identified in these flows. Velocities, lengths, and temporal spacings of individual waves and average velocities, lengths, frequencies, and intermittencies have been studied as functions of both gas and liquid flow rates. Velocity and frequency increase with liquid and gas flow rates, length decreases with increasing gas flow and increases with increasing liquid flow, and intermittency is predominantly an increasing function of liquid flow.