Boiling flow characteristics in microchannels with very hydrophobic surface to super-hydrophilic surface

Boiling flow process plays a very important role to affect the heat transfer in a microchannel. Different boiling flow modes have been found in the past which leads to different oscillations in temperatures and pressures. However, a very important issue, i.e. the surface wettability effects on the boiling flow modes, has never been discussed. The current experiments fabricated three different microchannels with identical sizes at 105 × 1000 × 30000 μm but at different wettability. The microchannels were made by plasma etching a trench on a silicon wafer. The surface made by the plasma etch process is hydrophilic and has a contact angle of 36° when measured by dipping a water droplet on the surface. The surface can be made hydrophobic by coating a thin layer of low surface energy material and has a contact angle of 103° after the coating. In addition, a vapor–liquid–solid growth process was adopted to grow nanowire arrays on the wafer so that the surface becomes super-hydrophilic with a contact angle close to 0°. Different boiling flow patterns on a surface with different wettability were found, which leads to large difference in temperature oscillations. Periodic oscillation in temperatures was not found in both the hydrophobic and the super-hydrophilic surface. During the experiments, the heat flux imposed on the wall varies from 230 to 354.9 kW/m² and the flow of mass flux into the channel from 50 to 583 kg/m²s. Detailed flow regimes in terms of heat flux versus mass flux are also obtained.     

Effect of copper surface wettability on the evaporation performance: Tests in a flat-plate heat pipe with visualization

The effects of copper surface wettability on the evaporation performance of a copper mesh wick were experimentally studied in an operating flat-plate heat pipe. Different degrees of wettability were obtained by varying the exposure times in air after the wicked plates were taken out of the sintering furnace. Three different working fluids: water, methanol and acetone, which possess different figures of merit, were investigated at the same volumetric liquid charge. The surface wettability was quantified by the static contact angle of sessile water drops on a flat copper surface. While the static contact angles of water drops varied from 10° to 40° for different degrees of wettability, the methanol and acetone drops still fully wetted the copper surface. A two-layer 100 + 200 mesh copper wick, 0.26 mm in thickness, was sintered on a 3 mm-thick copper base plate. A glass plate was adopted as the top wall of the heat pipe for visualization. Uniform heating was applied to the base plate near one end, and a cooling water jacket was connected at the other end. With increasing heat load, the evaporative resistance decreased with liquid film recession until a critical heat load showing the minimum evaporative resistance. Afterwards, partial dryout began from the front end of the evaporator. With decreasing wettability, the evaporating water film receded faster with increasing heat load and the critical heat loads were significantly reduced. In contrast, the critical heat loads for methanol and acetone seemed hardly affected by different wettability conditions. The minimum evaporative resistances, however, remained unaffected by surface wettability for all the three working fluids.     

Performance of a proposed complete wetting surface counter flow channel type liquid desiccant air dehumidifier

An idea that improves the wettability over the surfaces of a cylindrical dehumidifier channel was proposed and experimentally proved. Fibrous sheets were attached to the inner surfaces of the channel. The capillary effect of fibers sustains the complete wetting of the heat and mass transfer surfaces. The air to be dehumidified and cooled flows upward in the annulus space between the two layers of fibrous sheets, which are saturated with the downward flowing desiccant solution. The permeability of the fibrous sheet was determined experimentally. It was 2.43 × 10^?10 m². The measured solution flow rate due to the capillary suction of the sheets was Γ_in,min = 1.12 kg/h m. The liquid desiccant tested was H₂O/CaCl₂ with salt concentration ratios ranging from 35 to 40%. The measured distribution of the solution flow rate along the circumference of the sheets at the outlet showed 5% deviation from the average flow rate. This is a good indication for the good wettability of walls inside the dehumidifier.Feeding the solution by this mechanism has many advantages over spray feeding. Beside sustaining complete surface wetting, it also prevents channel blockage with solution, which is a main factor in increasing the air pressure drop. About 95% of the air pressure drop is saved in this study by avoiding these problems. A simple theoretical model for the heat and mass transfer processes inside the dehumidifier was developed and experimentally validated. In general, there is good agreement between the predicted and measured data. The developed model was utilized to study the effect of the different parameters on the dehumidifier performance. For a 1 m height dehumidifier with an inlet specific humidity and air temperature of 0.0234 kg_v/kg_a, and 35 °C, respectively, the predicted outlet air specific humidity was 0.0102 kg_v/kg_k and the corresponding outlet air temperature was 27.4 °C. The inlet solution temperature and salt concentration were 25 °C and 40%, respectively.     

Impinging jet study of the deposition of colloidal particles on modified polycarbonate and poly(ethylene terephthalate) surfaces

Main focus of this study was on characterization of surface properties of virgin and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) phospholipid layer coated poly(ethylene terephthalate) (PET) and polycarbonate (PC) planar articles. Surface properties were followed measuring static contact angles of wetting by means of sessile drop method and deposition of negative polystyrene (PS) colloidal particles followed by impinging jet method at defined flow regimes. It was found that phospholipid coating of both studied samples (PET, PC) let to the vigorous increase of the surface free energy. For coated samples major part of the surface free energy was dominated by polar component in contrast to the non-treated materials. Here the dispersive component was dominating. Results of the deposition experiments of polystyrene particles of 3 μm diameter correspond with trends obtained by contact angle measurements, i.e. the surface treated materials exhibited higher surface activity reflected in increased particle deposition rates. Simultaneously there was confirmed the fact, that with increasing magnitude of the Reynolds number of the dispersion flux the higher deposition rates were observed.     

Theoretical model for annular flow condensation in rectangular micro-channels

This study examines the pressure drop and heat transfer characteristics of annular condensation in rectangular micro-channels with three-sided cooling walls. A theoretical control-volume-based model is proposed based on the assumptions of smooth interface between the annular liquid film and vapor core, and uniform film thickness around the channel’s circumference. Mass and momentum conservation are applied to control volumes encompassing the liquid film and the vapor core separately. The model accounts for interfacial suppression of turbulent eddies due to surface tension with the aid of a new eddy diffusivity model specifically tailored to shear-driven turbulent films. The model predictions are compared with experimental pressure drop and heat transfer data for annular condensation of FC-72 along 1 × 1 mm² parallel channels. The condensation is achieved by rejecting heat to a counterflow of water. The data span FC-72 mass velocities of 248–367 kg/m² s, saturation temperatures of 57.8–62.3 °C, qualities of 0.23–1.0, and water mass flow rates of 3–6 g/s. The data are also compared to predictions of previous separated flow mini/micro-channel and macro-channel correlations. While some of the previous correlations do provide good predictions of the average heat transfer coefficient, they fail to capture axial variation of the local heat transfer coefficient along the channel. The new model accurately captures the pressure drop and heat transfer coefficient data in both magnitude and trend, evidenced by mean absolute error values of 3.6% and 9.3%, respectively.     

Adiabatic two-phase flow in rectangular microchannels with different aspect ratios: Part I – Flow pattern,pressure drop and void fraction

An aspect ratio is an important parameter for two-phase flow in a rectangular microchannel. To study the aspect ratio effect on the flow pattern, pressure drop and void fraction, experiments of adiabatic liquid water and nitrogen gas two-phase flow in rectangular microchannels were conducted. The widths and heights of rectangular microchannels are 510 μm × 470 μm, 608 μm × 410 μm, 501 μm × 237 μm and 503 μm × 85 μm. Therefore, the aspect ratios of the rectangular microchannels are 0.92, 0.67, 0.47 and 0.16; and the hydraulic diameters of the rectangular microchannels were 490, 490, 322 and 143 μm, respectively. Experimental ranges were liquid superficial velocities of 0.06–1.0 m/s and gas superficial velocities of 0.06–71 m/s. Visible rectangular microchannels were fabricated using a photosensitive glass. And pressure drop in microchannels was directly measured through embedded ports. The visualization of the flow pattern was carried out with a high-speed camera and a long distance microscope. Typical flow patterns in the rectangular microchannels observed in this study were bubble flow, transitional flow (multiple flow) and liquid ring flow. As the aspect ratio decreased, the bubble flow regime became dominant due to the confinement effect and the thickness of liquid film in corner was decreased. A void fraction in the rectangular microchannels has a linear relation with the volumetric quality. And the two-phase flow becomes homogeneous with decreasing aspect ratio owing to the reduction of the liquid film thickness. Like Zhang et al.’s [19] correlation, as the confinement number increased, the C-value in Lockhart and Martinelli correlation decreased. And a frictional pressure drop in the rectangular microchannels was highly related with the flow pattern.     

Experimental and analytical study of the effect of contact angle on liquid convective heat transfer in microchannels

In this paper we examine the effect of contact angle (or surface wettability) on the convective heat transfer coefficient in microchannels. Slip flow, where the fluid velocity at the wall is non-zero, is most likely to occur in microchannels due to its dependence on shear rate or wall shear stress. We show analytically that for a constant pressure drop, the presence of slip increases the Nusselt number. In a microchannel heat exchanger we modified the surface wettability from a contact angle of 20°–120° using thin film coating technology. Apparent slip flow is implied from pressure and flow rate measurements with a departure from classical laminar friction coefficients above a critical shear rate of approximately 10,000 s⁻¹. The magnitude of this departure is dependant on the contact angle with higher contact angles surfaces exhibiting larger pressure drop decreases. Similarly, the non-dimensional heat flux is found to decrease relative to laminar non-slip theory, and this decrease is also a function of the contact angle. Depending on the contact angle and the wall shear rate, variations in the heat transfer rate exceeding 10% can be expected. Thus the contact angle is an important consideration in the design of micro, and even more so, nano heat exchangers.     

Analysis of flow patterns emerging during evaporation in parallel microchannels

The evaporation processes of 2-propanol and water in cyclo olefin polymer (COP) and silicon microchannels of square cross-section are studied with a high-speed camera. The COP channels with a cross-section of 50 μm × 50 μm are rather smooth, whereas the 30 μm × 30 μm silicon channels have comparatively rough surfaces. For the COP channels, two different evaporation modes are identified, both with oscillating liquid–vapor menisci. One of these modes is characterized by an extremely rapid evaporation and a corresponding discontinuous shift of the meniscus. In the silicon channels four different evaporation modes are observed. Oscillatory motion of the liquid fronts also dominates here, and depending on the total mass flow and the wall temperature the oscillations in different channels are synchronized or desynchronized. Besides the flow patterns also the velocity trajectories of the evaporating liquid fronts are analyzed in detail and show a rather good reproducibility over different channels and different cycles. Compared to most other studies reported in this field, bubble nucleation is found to be of secondary importance for the evaporation processes.     

The importance of turbulence during condensation in a horizontal circular minichannel

Three-dimensional simulations of condensation of refrigerant R134a in a horizontal minichannel are presented. Mass fluxes ranging from 50 kg m^?2 s^?1 up to 1000 kg m^?2 s^?1 are considered in a circular minichannel of 1 mm diameter, and uniform wall and vapour–liquid interface temperatures are imposed as boundary conditions. The Volume of Fluid (VOF) method is used to track the vapour–liquid interface; the effects of interfacial shear stress, gravity and surface tension are taken into account. The influence of turbulence in the condensate film is analysed and compared against the assumption of laminar condensate flow by employing different computational approaches and validating the results against experimental data. Under the assumption of laminar condensate flow, experimental heat transfer coefficient values at low mass fluxes can be predicted, but the computed heat transfer coefficient is found to be almost independent of mass flux and vapour quality. Only when turbulence in the condensate film is taken into account does the numerical model capture the influence of mass flux that is observed in the experimental measurements.     

Theoretical and experimental study of the effects of spray inclination on two-phase spray cooling and critical heat flux

Experiments were conducted with PF-5052 liquid sprays impacting a 1.0 × 1.0 cm² heated test surface at different inclination angles, flow rates, and subcoolings. Inclination angle had no noticeable effect on the single-phase or two-phase regions of the boiling curve. Maximum CHF was always achieved with the spray impinging normal to the test surface; increasing angle of inclination away from the normal decreased CHF appreciably. Video analysis showed inclined sprays produced lateral liquid film flow towards the farthest downstream region of the test surface. The film liquid provided partial resistance to dryout despite the weak volumetric spray flux in the downstream region. A new theoretical model of the spray’s impact area and volumetric flux proves this decrease is the result of a sharp reduction in the fraction of the test surface area that is directly impacted by the spray. Combining the model and video results with a previous point-based CHF correlation for normal sprays is shown to accurately predict the effects of orientation angle on CHF for different nozzles and operating conditions.     

Determination of annular condensation heat transfer coefficient of steam in microchannels with trapezoidal cross sections

Experiments have been carried out to determine annular condensation heat transfer coefficient of steam in two silicon microchannels having trapezoidal cross sections with the same aspect ratio of 3.15 at 54 < G < 559 kg/m² s under 3-side cooling conditions. A semi-analytical method, based on turbulent flow boundary layer theory of liquid film with correlations of pressure drop and void fraction valid for microchannels, is used to derive the annular local condensation heat transfer coefficients. The predicted values based on the semi-analytical model are found within ±20% of 423 data points. It is shown that the annular condensation heat transfer coefficient in a microchannel increases with mass flux and quality and decreases with the hydraulic diameter.     

A novel time strip flow visualisation technique for investigation of intermittent dewetting and dryout in elongated bubble flow in a microchannel evaporator

A flow visualisation study of flow boiling of R245fa in silicon multi-microchannels at low mass flux and moderate heat flux has been carried out with a high speed digital camera. The micro-evaporator had 67 channels of length 20 mm, width 223 μm, and height 680 μm while the fin width between adjacent channels was 80 μm. The base heat flux ranged from 2 to 26 W cm^?2 for a mass velocity of 100 kg s^?1 m^?2, resulting in exit vapour qualities ranging from 10% to 70%. In particular, a novel time strip technique was developed to analyse the recorded image sequences and significantly highlight the various phenomena occurring along given channels. Notably, this technique was able to reveal profound details regarding the intermittent dryout mechanism of liquid films trapped between the elongated bubbles and the heated channel walls. The results show that the intermittent dryout of the evaporating liquid film is comprised of four stages with distinct time scales and dynamics: (i) the growth of liquid film thinning perturbations to a critical amplitude causing the rupture of the metastable liquid film, (ii) a dewetting stage involving expanding dry spots leading to a rivulet flow regime, (iii) evaporation of the rivulets leading to full dryout, and (iv) a rewetting stage. This intermittent dryout mechanism appears to explain the many seemingly contradictory heat transfer coefficient trends observed with changes in vapour quality in microchannels, thus resolving an important heat transfer dilemma. Furthermore, since dryout is an undesirable event during the practical application of a microchannel evaporator, it is important to delay or even suppress the initial rupture of the liquid film that leads to dryout. This can be achieved by manufacturing or treating the channel surfaces to be highly wettable with the chosen refrigerant.     

Hydrodynamics of quenching with impinging free-surface jet

The hydrodynamics of jet impingement quenching of a stainless steel specimen has been studied experimentally. The specimen is heated to an initial temperature of about 900 °C and then quenched by a subcooled free-surface water jet. High-speed imaging shows that the free-surface of the water film in the wetted region is smooth. The water film outside the wetted region is deflected away from the surface and then breaks into droplets due to surface tension and shear forces. The splashed droplet velocity is found to be low up to a wetting front radius of 6 mm (r/d_J ≈ 2), beyond which it increases rapidly before reaching a constant value at a wetting front radius of about 8 to 10 mm (2.67 ? r/d_J ? 3.34). The water film velocity at the wetting front is calculated using the single-phase boundary layer model suggested by Watson [2]. At moderate subcooling, the splashed droplet velocity up to a wetting front radius of 10 mm (r/d_J ≈ 3.34) is found to be much lower than the estimated single-phase film velocity. The study shows that although the wetted region may appear devoid of any bubbles, strong two-phase flow occurs within this region.     

Heat transfer of gas–liquid mixture in micro-channel heat sink

The main objective of the present investigation is to study heat transfer in parallel micro-channels of 0.1 mm in size. Comparison of the results of this study to the ones obtained for two-phase flow in “conventional” size channels provides information on the complex phenomena associated with heat transfer in micro-channel heat sinks. Two-phase flow in parallel micro-channels, feeding from a common manifold shows that different flow patterns occur simultaneously in the different micro-channels: liquid alone (or single-phase flow), bubbly flow, slug flow, and annular flow (gas core with a thin liquid film, and a gas core with a thick liquid film). Although the gas core may occupy almost the entire cross-section of the triangular channel, making the side walls partially dry, the liquid phase always remained continuous due to the liquid, which is drawn into the triangular corners by surface tension. With increasing superficial gas velocity, a gas core with a thin liquid film is observed. The visual observation showed that as the air velocity increased, the liquid droplets entrained in the gas core disappeared such that the flow became annular. The probability of appearance of different flow patterns should be taken into account for developing flow pattern maps. The dependence of the Nusselt number, on liquid and gas Reynolds numbers, based on liquid and gas superficial velocity, respectively, was determined in the range of Re_LS = 4–56 and Re_GS = 4.7–270. It was shown that an increase in the superficial liquid velocity involves an increase in heat transfer (Nu_L). This effect is reduced with increasing superficial gas velocity, in contrast to the results reported on two-phase heat transfer in “conventional size” channels.     

Correlation of critical heat flux and two-phase friction factor for subcooled convective boiling in structured surface microchannels

Critical heat flux (CHF) and pressure drop of subcooled flow boiling are measured for a microchannel heat sink containing 75 parallel 100 μm × 200 μm structured surface channels. The heated surface is made of a Cu metal sheet with/without 2 μm thickness diamond film. Tests and measurements are conducted with de-ionized water, de-ionized water +1 vol.% MCNT additive solution, and FC-72 fluids over a mass velocity range of 820–1600 kg/m² s, with inlet temperatures of 15(8.6)°C, 25(13.6)°C, 44(24.6)°C, and 64(36.6)°C for DI water (FC-72), and heat fluxes up to 600 W/cm². The CHF of subcooled flow boiling of the test fluids in the microchannels is measured parametrically. The two-phase pressure drop is also measured. Both CHF and the two-phase friction factor correlation for one-side heating with two other side-structured surface microchannels are proposed and developed in terms of the relevant parameters.     

Wall heat flux partitioning during subcooled forced flow film boiling of water on a vertical surface

Subcooled flow film boiling experiments were conducted on a vertical flat plate, 30.5 cm in height, and 3.175 cm wide with forced convective upflow of subcooled water at atmospheric pressure. Data have been obtained for mass fluxes ranging from 0 to 700 kg/m²s, inlet subcoolings ranging from 0 to 25 °C and wall superheats ranging from 200 to 400 °C. Correlations for wall heat transfer coefficient and wall heat flux partitioning were developed as part of this work. These correlations derive their support from simultaneous measurements of the wall heat flux, fluid temperature profiles, liquid side heat flux and interfacial wave behavior during steady state flow film boiling. A new correlation for the film collapse temperature was also deduced by considering the limiting case of heat flux to the subcooled liquid being equal to the wall heat flux. The premise of this deduction is that film collapse under subcooled conditions occurs when there is no net vapor generation. These correlations have also been compared with the data and correlations available in the literature.     

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.     

Transient analysis for the cathode gas diffusion layer of PEM fuel cells

A one-dimensional, non-isothermal, two-phase transient model has been developed to study the transient behaviour of water transport in the cathode gas diffusion layer of PEM fuel cells. The effects of four parameters, namely the liquid water saturation at the interface of the gas diffusion layer and flow channels, the proportion of liquid water to all of the water at the interface of the cathode catalyst layer and the gas diffusion layer, the current density, and the contact or wetting angle, on the transient distribution of liquid water saturation in the cathode gas diffusion layer are investigated. Especially, the time needed for liquid water saturation to reach steady state and the liquid water saturation at the interface of the cathode catalyst layer and gas diffusion layer are plotted as functions of the above four parameters. The ranges of water vapour condensation and liquid water evaporation are identified across the thickness of the gas diffusion layer. In addition, the effects of the above four parameters on the steady state distributions of gas phase pressure, water vapour concentration, oxygen concentration and temperature are also presented. It is found that increasing any one of the first three parameters will increase the water saturation at the interface of the catalyst layer and gas diffusion layer, but decrease the time needed for the liquid water saturation to reach steady state. When the liquid water saturation at the interface of the gas diffusion layer and flow channels is high enough (≥0.1), the liquid water saturation at steady state is almost uniformly distributed across the thickness of the gas diffusion layer. It is also found that, under the given initial and boundary conditions in this paper, evaporation takes place within the gas diffusion layer close to the channel side and is the major process for water phase change at low current density (<2000 A m⁻²); condensation occurs close to the catalyst layer side within the gas diffusion layer and dominates the phase change at high current density (>5000 A m⁻²). For hydrophilic gas diffusion layers, both the time needed for liquid water saturation to reach steady state and the water saturation at the interface of the catalyst layer and gas diffusion layer will increase when the contact angle increases; but for hydrophobic gas diffusion layers, both of them decrease when the contact angle increases.     

Boiling two-phase flow and efficiency of co- and counter-current microchannel heat exchangers with gas heating

To learn how to utilize the exhaust heat from a high-temperature gas product of a methanol reformer, the present study experimentally investigates the boiling two-phase flow in co- and counter-current microchannel heat exchangers (MCHE) with gas heating. Boiling two-phase flow patterns, two-phase flow instability, and efficiency are explored. The working fluid on the hot and cold sides are helium and liquid methanol, respectively. The silicon-based MCHE, which has dimensions of 20 mm × 20 mm, is designed with 18 parallel microchannels on both sides and is prepared using microfabrication processes. Four types of two-phase flow patterns – bubbly-elongated slug flow, annular flow, annular flow with liquid film breakup, and dryout are identified in both types of MCHE that are studied. A flow pattern map is then constructed on the plane of the methanol mass flux versus heat flux for both types of MCHE. In the counter-current MCHE, the efficiency increases significantly with an increase in the mass flux in both the single- and two-phase flow regions, while the effect of mass flux is insignificant in the co-current MCHE. In the two-phase flow region, the efficiency of both types of MCHEs gradually increases with an increase in the hot-side thermal power until the CHF is approached. The highest efficiency obtained in the present study is about 0.85 and 0.90 for the co- and counter-current MCHEs, respectively.     

An experimental investigation of sessile droplets evaporation on hydrophilic and hydrophobic heating surface with constant heat flux

Droplet evaporation widely exists in the daily life and industrial production. In most of previous experimental studies, the evaporation of sessile droplets was conducted under a constant substrate temperature condition. However, drops often evaporating on a heating surface under a constant heat flux condition in many practical applications. In this paper, we have carried out an experiment on sessile 3 μl DI water droplets evaporated on hydrophilic and hydrophobic heating surfaces under constant heat flux in the range from 1153 W/m² to 6919 W/m². A high-speed camera was used to record the changing shapes of two sessile droplets on a hydrophilic and a hydrophobic heating surface placed side by side. The droplet height, dynamic contact angle, droplet contact diameter, evaporation mode and evaporation rate are presented.