Impact of aspect ratio on flow boiling of water in rectangular microchannels

In this paper we focus on the impact of varying the aspect ratio of rectangular microchannels, on the overall pressure drop involving water boiling. An integrated system comprising micro-heaters, sensors and microchannels has been realized on (1 1 0) silicon wafers, following CMOS compatible process steps. Rectangular microchannels were fabricated with varying aspect ratios (width [W] to depth [H]) but constant hydraulic diameter of 142 ± 2 μm and length of 20 mm. The invariant nature of the hydraulic diameter is confirmed through two independent means: physical measurements using profilometer and by measuring the pressure drop in single-phase fluid flow. The experimental results show that the pressure drop for two-phase flow in rectangular microchannels experiences minima at an aspect ratio of about 1.6. The minimum is possibly due to opposing trends of frictional and acceleration pressure drops, with respect to aspect ratio. In a certain heat flux and mass flux range, it is observed that the two-phase pressure drop is lower than the corresponding single-phase value. This is the first study to investigate the effect of aspect ratio in two-phase flow in microchannels, to the best of our knowledge. The results are in qualitative agreement with annular flow model predictions. These results improve the possibility of designing effective heat-sinks based on two-phase fluid flow in microchannels.     

Investigations of phase inversion and frictional pressure gradients in upward and downward oil–water flow in vertical pipes

The present study has attempted to investigate phase inversion and frictional pressure gradients during simultaneous vertical flow of oil and water two-phase through upward and downward pipes. The liquids selected were white oil (44 mPa s viscosity and 860 kg/m³ density) and water. The measurements were made for phase velocities varying from 0 to 1.24 m/s for water and from 0 to 1.87 m/s for oil, respectively. Experiments were carried either by keeping the mixture velocity constant and increasing the dispersed phase fraction or by keeping the continuous phase superficial velocity constant and increasing the dispersed phase superficial velocity. From the experimental results, it is shown that the frictional pressure gradient reaches to its lower value at the phase inversion point in this work. The points of phase inversion are always close to an input oil fraction of 0.8 for upward flow and of 0.75 for downward flow, respectively. A few models published in the literature are used to predict the phase inversion point and to compare the results with available experimental data. Suitable methods are suggested to predict the critical oil holdup at phase inversion based on the different viscosity ratio ranges. Furthermore, the frictional pressure gradient is analyzed with several suitable theoretical models according to the existing flow patterns. The analysis reveals that both the theoretical curves and the experimental data exhibit the same trend and the overall agreement of predicted values with experimental data is good, especially for a high oil fraction.     

Heat transfer and fluid dynamics of air–water two-phase flow in micro-channels

Heat transfer, pressure drop, and void fraction were simultaneously measured for upward heated air–water non-boiling two-phase flow in 0.51 mm ID tube to investigate thermo–hydro dynamic characteristics of two-phase flow in micro-channels. At low liquid superficial velocity j_l frictional pressure drop agreed with Mishima–Hibiki’s correlation, whereas agreed with Chisholm–Laird’s correlation at relatively high j_l. Void fraction was lower than the homogeneous model and conventional empirical correlations. To interpret the decrease of void fraction with decrease of tube diameter, a relation among the void fraction, pressure gradient and tube diameter was derived. Heat transfer coefficient fairly agreed with the data for 1.03 and 2.01 mm ID tubes when j_l was relatively high. But it became lower than that for larger diameter tubes when j_l was low. Analogy between heat transfer and frictional pressure drop was proved to hold roughly for the two-phase flow in micro-channel. But satisfactory relation was not obtained under the condition of low liquid superficial velocity.     

Positive frictional pressure gradient in vertical gas-high viscosity oil slug flow

The present study describes the wall shear stress and the falling liquid film behavior in upward vertical slug flow of air and high viscosity oil. The frictional pressure gradient is directly related to the wall shear stress, and it is usually negative (opposite to the overall flow direction). However, in vertical slug flow, the average total wall shear stress of a slug unit may be negative (in the same direction of the overall flow), resulting in a positive frictional pressure gradient. However, this does not mean, by any way, generation of additional energy or violation of the second law of thermodynamics.The positive frictional pressure gradient phenomenon, reasons and required conditions were explained in this paper. A simplified model was developed and validated against recent experimental data of air-high viscosity oil slug flow in a 50.8 mm ID vertical pipe. The oil viscosity was in the range of 127 mPa s to 580 mPa s. Positive frictional pressure gradient appears when the liquid film wall shear stress supersede the wall shear stress in the slug body. The rate of increase of both wall shear stresses (with respect to the mixture Reynolds number) depend, not only, on the mixture Reynolds number but also, highly, on the liquid viscosity.     

Two-phase pressure drop and flow visualization of FC-72 in a silicon microchannel heat sink

The rapid development of two-phase microfluidic devices has triggered the demand for a detailed understanding of the flow characteristics inside microchannel heat sinks to advance the cooling process of micro-electronics. The present study focuses on the experimental investigation of pressure drop characteristics and flow visualization of a two-phase flow in a silicon microchannel heat sink. The microchannel heat sink consists of a rectangular silicon chip in which 45 rectangular microchannels were chemically etched with a depth of 276 μm, width of 225 μm, and a length of 16 mm. Experiments are carried out for mass fluxes ranging from 341 to 531 kg/m² s and heat fluxes from 60.4 to 130.6 kW/m² using FC-72 as the working fluid. Bubble growth and flow regimes are observed using high speed visualization. Three major flow regimes are identified: bubbly, slug, and annular. The frictional two-phase pressure drop increases with exit quality for a constant mass flux. An assessment of various pressure drop correlations reported in the literature is conducted for validation. A new general correlation is developed to predict the two-phase pressure drop in microchannel heat sinks for five different refrigerants. The experimental pressure drops for laminar-liquid laminar-vapor and laminar-liquid turbulent-vapor flow conditions are predicted by the new correlation with mean absolute errors of 10.4% and 14.5%, respectively.     

Frictional pressure drop during vapour–liquid flow in minichannels: Modelling and experimental evaluation

Condensation in minichannels is widely used in air-cooled condensers for the automotive and air-conditioning industry, in heat pipes and other applications for system thermal control. The knowledge of pressure drops in such small channels is important in order to optimize heat transfer surfaces. This paper presents a model for calculation of the frictional pressure gradient during condensation or adiabatic liquid–gas flow inside minichannels with different surface roughness. In order to account for the effects of surface roughness, new experimental frictional pressure gradient data associated to single-phase flow and adiabatic two-phase flow of R134a inside a single horizontal mini tube with rough wall has been used in the modelling. It is a Friedel (1979) [Friedel, L., 1979. Improved friction pressure drop correlations for horizontal and vertical two-phase pipe flow. In: Proceedings of the European Two-Phase Flow Group Meeting, Ispra, Paper E2] based model and it takes into account mass velocity, vapor quality, fluid properties, reduced pressure, tube diameter, entrainment ratio and surface roughness. With respect to the flow pattern prediction capability, it has been built for shear dominated flow regimes inside pipes, thus, annular, annular-mist and mist flow are here predicted. However, the suggested procedure is extended to the intermittent flow in minichannels and it is also applied with success to horizontal macro tubes.     

Flow pattern based correlations of two-phase pressure drop in rectangular microchannels

Numerous pressure drop correlations for microchannels have been proposed; most of them can be classified as either a homogeneous flow model (HFM) or a separated flow model (SFM). However, the predictions of these correlations have not been compared directly because they were developed in experiments conducted under a range of conditions, including channel shape, the number of channels, channel material and the working fluid. In this study, single rectangular microchannels with different aspect ratios and hydraulic diameters were fabricated in a photosensitive glass. Adiabatic water-liquid and Nitrogen-gas two-phase flow experiments were conducted using liquid superficial velocities of 0.06–1.0 m/s, gas superficial velocities of 0.06–72 m/s and hydraulic diameters of 141, 143, 304, 322 and 490 μm. A pressure drop in microchannels was directly measured through embedded ports. The flow pattern was visualized using a high-speed camera and a long-distance microscope. A two-phase pressure drop in the microchannel was highly related to the flow pattern. Data were used to assess seven different HFM viscosity models and ten SFM correlations, and new correlations based on flow patterns were proposed for both HFMs and SFMs.     

Flow patterns and pressure drop of ionic liquid–water two-phase flows in microchannels

The two-phase flow of a hydrophobic ionic liquid and water was studied in capillaries made of three different materials (two types of Teflon, FEP and Tefzel, and glass) with sizes between 200 μm and 270 μm. The ionic liquid was 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide, with density and viscosity of 1420 kg m⁻³ and 0.041 kg m⁻¹ s⁻¹, respectively. Flow patterns and pressure drop were measured for two inlet configurations (T- and Y-junction), for total flow rates of 0.065–214.9 cm³ h⁻¹ and ionic liquid volume fractions from 0.05 to 0.8. The continuous phase in the glass capillary depended on the fluid that initially filled the channel. When water was introduced first, it became the continuous phase with the ionic liquid forming plugs or a mixture of plugs and drops within it. In the Teflon microchannels, the order that fluids were introduced did not affect the results and the ionic liquid was always the continuous phase. The main patterns observed were annular, plug, and drop flow. Pressure drop in the Teflon microchannels at a constant ionic liquid flow rate, was found to increase as the ionic liquid volume fraction decreased, and was always higher than the single phase ionic liquid value at the same flow rate as in the two-phase mixture. However, in the glass microchannel during plug flow with water as the continuous phase, pressure drop for a constant ionic liquid flow rate was always lower than the single phase ionic liquid value. A modified plug flow pressure drop model using a correlation for film thickness derived for the current fluids pair showed very good agreement with the experimental data.     

Effect of gravity on the mechanical properties of lunar regolith tested using a low gravity simulation device

Simulations of the bearing capacity and shear strength of regolith under Earth’s gravity produce different results from those under low gravity. A low-gravity simulation device was developed in this study, and an internal stress model of regolith simulant was established to correct the errors. The model revealed additional force on both shear plane in the shear test and the press plate area in the pressure–sinkage test. The sinkage and shear test results showed that low gravity decreased the deformable index n, frictional modulus k_φ and cohesion c, whereas there were no obvious changes to the cohesive modulus k_c and internal friction angle φ. The sinkage generally increased as the gravity decreased under a consistent normal load larger than 50 N, but when the wheel load was lower than 50 N, the sinkage of the TYII-1 simulant was larger under 1 G than 1/6 G. Gravity had little effect on the shear strength of the regolith. However, the tractive thrust of the TYII-1 simulant was lower under 1/6 G than 1 G. The smaller difference was due to differences in the way the soils responded to changes in the gravity level for the TYII-2 simulant.     

Pressure drop prediction in annular two-phase flow in macroscale tubes and channels

A new prediction method for the frictional pressure drop in annular two-phase flow is presented. This new prediction method focuses on the aerodynamic interaction between the liquid film and the gas core in annular flows, and explicitly takes into account the asymmetric liquid film distribution in the tube cross section induced by the action of gravity in horizontal tubes operated at low mass fluxes. The underlying experimental database contains 6291 data points from the literature with 13 fluid combinations (both single-component saturated fluids such as water, carbon dioxide and refrigerants R12, R22, R134a, R245fa, R410a, R1234ze, and two-component fluids such as water-argon, water-nitrogen, alcohol-argon, water plus alcohol-argon and water-air), vertical and horizontal tubes and annuli with diameters from 3 mm to 25 mm, and both adiabatic and evaporating flow conditions. The new prediction method is very simple to implement and use, is physically based and outperforms existing pressure drop correlations (mean absolute error of 12.9%, and 7 points out of 10 captured to within ±15%).     

Gas–liquid pressure drop in vertical internally wavy 90° bend

Experiments of air water two-phase flow pressure drop in vertical internally wavy 90° bend have been carried out. The tested bends are flexible and made of stainless steel with inner diameter of 50 mm and various curvature radiuses of 200, 300, 400 and 500 mm. The experiments were performed under the following conditions of two-phase parameters; mass flux from 350 to 750 kg/m² s. Gas quality from 1% to 50% and system pressure from 4 to 7.5 bar. The results demonstrate that the effect of the above-mentioned parameters is very significant at high ranges of mass flow quality. Due to the increasing of two-phase flow resistance, energy dissipations, friction losses and interaction of the two-phases in the vertical internally wavy 90° bend the total pressure drops are perceptible about 2–5 times grater than that in smooth bends. Based on the mass and energy balance as well as the presented experimental results, new empirical correlation has been developed to calculate the two-phase pressure drop and hence the two-phase friction factor of the tested bends. The correlation includes the relevant primary parameter, fit the data well, and is sufficiency accurate for engineering purposes.     

Low-GWP refrigerants flow boiling heat transfer in a 5 PPI copper foam

This paper reports an experimental investigation of the heat transfer performance of the new low-GWP refrigerants, R1234yf and R1234ze(E), during flow boiling heat transfer inside a horizontal high porosity copper foam with 5 Pores Per Inch (PPI). Metal foams are a class of cellular structured materials consisting of a stochastic distribution of interconnected pores; these materials have been proposed as effective solutions for heat transfer enhancement during both single and two-phase heat transfer. R1234yf and R1234ze(E) refrigerants are appealing alternatives of the more traditional R134a by virtue of their negligible values of GWP and normal boiling temperatures close to that of R134a, which make them suitable solution in several different applications, such as: refrigeration and air conditioning and electronic thermal management. This work compares the two-phase heat transfer behaviour of these new HFO refrigerants, studying the boiling process inside a porous medium and permitting to understand their effective heat transfer capabilities. The experimental measurements were carried out by imposing three different heat fluxes: 50, 75, and 100 kW m⁻², at a constant saturation temperature of 30 °C; the refrigerant mass velocity was varied between 50 and 200 kg m⁻² s⁻¹, whilst the mean vapour quality varied from 0.2 to 0.95. The two-phase heat transfer and pressure drop performance of the two new HFO refrigerants is compared against that of the more traditional R134a.     

Experimental study on the effect of drag reducing polymer on flow patterns and drag reduction in a horizontal oil–water flow

In this study, a HMW anionic co-polymer of 40:60 wt/wt NaAMPS/acrylamide was used as a drag reducing polymer (DRP) for oil–water flow in a horizontal 25.4 mm ID acrylic pipe. The effect of polymer concentration in the master solution and after injection in the main water stream, oil and water velocities, and pipe length on drag reduction (DR) was investigated. The injected polymer had a noticeable effect on flow patterns and their transitions. Stratified and dual continuous flows extended to higher superficial oil velocities while annular flow changed to dual continuous flow. The results showed that as low as 2 ppm polymer concentration was sufficient to create a significant drag reduction across the pipe. DR was found to increase with polymer concentration increased and reached maximum plateau value at around 10 ppm. The results showed that the drag reduction effect tends to increase as superficial water velocity increased and eventually reached a plateau at U_sw of around 1.3 m/s. At U_sw > 1.0 m/s, the drag reduction decreased as U_so increased while at lower water velocities, drag reduction is fluctuating with respect to U_so. A maximum DR of about 60% was achieved at U_so = 0.14 m/s while only 45% was obtained at U_so = 0.52 m/s. The effectiveness of the DRP was found to be independent of the polymer concentration in the master solution and to some extent pipe length. The friction factor correlation proposed by Al-Sarkhi et al. (2011) for horizontal flow of oil–water using DRPs was found to underpredict the present experimental pressure gradient data.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Investigation of a continuous adjoint-based optimization procedure for aeroacoustic control of plane jets

A control optimization technique using the continuous adjoint of the compressible Navier–Stokes equations is implemented for aeroacoustic optimization of plane jet flows. The purpose of the adjoint equations is to provide sensitivity information, which is afterwards used in a gradient-based minimization of a prescribed cost functional, designed to describe the far-field sound pressure level (SPL). The objective of the present paper is to demonstrate the ability to reduce the sound in the near far-field of plane jets. Furthermore, as the continuous adjoint approach can become inaccurate, due to inconsistencies between the continuous and the discretized system, the accuracy of the continuous adjoint approach is investigated. The considered cases exhibit a nozzle exit Reynolds number of Re_jet = ρu_jetD/μ = 2000 and a Mach number of M_jet = 0.9, performed using two-dimensional direct numerical simulation and three-dimensional large-eddy simulation, respectively. A comparison of the obtained gradient via adjoint and finite differences is presented and it is shown, that in order to obtain reliable gradient directions, the length of the optimization time needs to be restricted. Furthermore, a receding horizon optimization for the two-dimensional plane jet simulation is used to obtain a sound reduction over much longer time intervals. The influence of different formulations of the viscosity in the adjoint equations is finally investigated.     

Evolution of turbulence characteristics from straight to curved pipes

Fully developed, statistically steady turbulent flow in straight and curved pipes at moderate Reynolds numbers is studied in detail using direct numerical simulations (DNS) based on a spectral element discretisation. After the validation of data and setup against existing DNS results, a comparative study of turbulent characteristics at different bulk Reynolds numbers Re_b = 5300 and 11,700, and various curvature parameters κ = 0, 0.01, 0.1 is presented. In particular, complete Reynolds-stress budgets are reported for the first time. Instantaneous visualisations reveal partial relaminarisation along the inner surface of the curved pipe at the highest curvature, whereas developed turbulence is always maintained at the outer side. The mean flow shows asymmetry in the axial velocity profile and distinct Dean vortices as secondary motions. For strong curvature a distinct bulge appears close to the pipe centre, which has previously been observed in laminar and transitional curved pipes at lower Re_b only. On the other hand, mild curvature allows the interesting observation of a friction factor which is lower than in a straight pipe for the same flow rate.All statistical data, including mean profile, fluctuations and the Reynolds-stress budgets, is available for development and validation of turbulence models in curved geometries.     

Transient heating effects in high pressure Diesel injector nozzles

The tendency of today’s fuel injection systems to reach injection pressures up to 3000 bar in order to meet forthcoming emission regulations may significantly increase liquid temperatures due to friction heating; this paper identifies numerically the importance of fuel pressurization, phase-change due to cavitation, wall heat transfer and needle valve motion on the fluid heating induced in high pressure Diesel fuel injectors. These parameters affect the nozzle discharge coefficient (C_d), fuel exit temperature, cavitation volume fraction and temperature distribution within the nozzle. Variable fuel properties, being a function of the local pressure and temperature are found necessary in order to simulate accurately the effects of depressurization and heating induced by friction forces. Comparison of CFD predictions against a 0-D thermodynamic model, indicates that although the mean exit temperature increase relative to the initial fuel temperature is proportional to (1 − C_d²) at fixed needle positions, it can significantly deviate from this value when the motion of the needle valve, controlling the opening and closing of the injection process, is taken into consideration. Increasing the inlet pressure from 2000 bar, which is the pressure utilized in today’s fuel systems to 3000 bar, results to significantly increased fluid temperatures above the boiling point of the Diesel fuel components and therefore regions of potential heterogeneous fuel boiling are identified.     

Coalescence of sessile droplets of varying viscosities for line printing

Coalescence of sessile droplets is studied experimentally with water–glycerin mixtures of different viscosities. Effects of viscosity on the dimensionless spreading length (Ψ) and the center-to-center distance (L) are investigated for two droplets; the first droplet (D_s) is stationary on a substrate and the second droplet (D₀) landing at a center-to-center distance L from the first droplet. For a low viscosity fluid, Ψ is maximum when L approaches zero (or λ → 1, where λ = 1 − L/D_s), which represents a head-on collision. For a high viscosity fluid, Ψ is minimum when λ → 0.6. The effect of λ on line printing for various viscosities is also examined by printing multiple droplets. We found that the larger the viscosity, the less the breakup between droplets; viscosities smaller than 60 wt% glycerin yielded line breakup. The overlap ratio of λ > 0.3 produced not a line, but a bigger droplet or puddle because of coalescence. Data obtained in this work can provide insights for the fabrication of conductive microtracks or microinterconnects in printed-electronics applications where a line breakup between droplets would lead to an electrical circuit short.     

One-dimensional interfacial area transport of vertical upward bubbly flow in narrow rectangular channel

The design and safety analysis for miniature heat exchangers, the cooling system of high performance microelectronics, research nuclear reactors, fusion reactors and the cooling system of the spallation neutron source targets requires the knowledge of the gas–liquid two-phase flow in a narrow rectangular channel. In this study, flow measurements of vertical upward air–water flows in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm were performed at seven axial locations by using the imaging processing technique. The local frictional pressure loss gradients were also measured by a differential pressure cell. In the experiment, the superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. The developing two-phase flow was characterized by the significant axial changes of the local flow parameters due to the bubble coalescence and breakup in the tested flow conditions. The existing two-phase frictional multiplier correlations such as Chisholm, 1967, Mishima et al., 1993 and Lee and Lee (2001) were verified to give a good prediction for the measured two-phase frictional multiplier. The predictions of the drift-flux model with the rectangular channel distribution parameter correlation of Ishii (1977) and several existing drift velocity correlations of Ishii, 1977, Hibiki and Ishii, 2003 and Jones and Zuber (1979) agreed well with the measured void fractions and gas velocities. The interfacial area concentration (IAC) model of Hibiki and Ishii (2002) was modified by taking the channel width as the system length scale and the modified IAC model could predict the IAC and Sauter mean diameter acceptably.