Flow boiling behaviors in hydrophilic and hydrophobic microchannels

Surface wettability is a critical parameter in small scale phenomena, especially two-phase flow, since the surface force becomes dominant as size decreases. In present study, experiments of water flow boiling in hydrophilic and hydrophobic rectangular microchannels were conducted to investigate the wettability effect on flow boiling in rectangular microchannels. The rectangular microchannels were fabricated with a photosensitive glass to visualize flow pattern. The hydrophilic bare photosensitive glass microchannel was chemically treated to obtain a hydrophobic microchannel. And, visualization of flow patterns was carried out. And boiling heat transfer and two-phase pressure drop was analyzed with visualization results. The boiling heat transfer coefficient in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with nucleation site density and liquid film motion. And the pressure drop in the hydrophobic rectangular microchannel was higher than that in the hydrophilic rectangular microchannel, which was highly related with unstable motions of bubble and liquid film. Finally, we find out the wettability is important parameter on the flow pattern, which were highly related with two-phase heat and mass transfer.     

Surface wettability effect on flow pattern and pressure drop in adiabatic two-phase flows in rectangular microchannels with T-junction mixer

Wettability is an important parameter in micro-scale flow patterns. Previous research has usually been conducted in conventional microtubes due to limitations of visualizing flow patterns and fabricating microchannels. However, most microchannels in practical applications have rectangular shape. Furthermore, pressure drop is closely related with flow pattern. Hence, we studied water liquid and nitrogen gas flows in rectangular microchannels with different wettabilities. The rectangular glass microchannels were fabricated from photosensitive glass, whose surface is hydrophilic. The surface of one was silanized using octadecyl-trichloro-silane (OTS) to prepare a hydrophobic microchannel. The two-phase flow pattern was visualized with a high-speed camera and a long distance microscope. The frictional pressure drop in the microchannel was measured directly with embedded pressure ports. The flow pattern and pressure drop in the hydrophobic microchannel were totally different from those in the hydrophilic microchannel. Finally, the two-phase frictional pressure drop was analyzed based on the flow patterns of different wettabilities.     

Experimental investigation of flow and heat transfer for the ethanol-water solution and FC-72 in rectangular microchannels

Rapid development of super scale integration circuit (IC) provides unprecedented challenge to thermal control for aviation electronic equipments. To solve the problem of cooling electronic chips and devices for aircraft avionics, this paper experimentally investigated the characteristics of single-phase forced convection heat transfer and flow resistance in rectangular microchannels with two liquid coolants. One was 30% of ethanol–water solution, the most commonly used coolant in aviation. The other was FC-72, the latest coolant for electronic equipments. Based on the experimental data collected and those available in the open literature, comparisons and analyses were carried out to evaluate the influences of liquid velocity, supercooling temperature, microchannel structures and wall temperature etc. on the heat transfer behaviors. And the correlations of flow resistance and heat transfer characteristics were provided for the ethanol–water solution and FC-72 respectively. The results indicate transition from laminar to turbulent flow occurs at the Reynolds number of 750–1,250 for FC-72, and the behaviors of flow and heat transfer in rectangular microchannels strongly depend on the kind of coolant and geometric configuration of microchannels.     

Experimental validation of a two equation RANS transitional turbulence model for compressible microflows

Laminar-to-turbulent flow transition in microchannels can be useful to enhance mixing and heat transfer in microsystems. Typically, the small characteristic dimensions of these devices hinder in attaining higher Reynolds numbers to limit the total pressure drop. This is true especially in the presence of a liquid as a working medium. On the contrary, due to lower density, Reynolds number larger than 2000 can be easily reached for gas microflows with an acceptable pressure drop. Since microchannels are used as elementary building blocks of micro heat exchangers and micro heat-sinks, it is essential to predict under which conditions, the laminar-to-turbulent flow transition inside such geometries can be expected. In this paper, experimental validation of a two equations transitional turbulence model, capable of predicting the laminar-to-turbulent flow transition for internal flows as proposed by Abraham etal. (2008), is presented for the first time for microchannels. This is done by employing microchannels in which Nitrogen gas is used as a working fluid. Two different cross-sections namely circular and rectangular are utilized for numerical and experimental investigations. The inlet mass flow rate of the gas is varied to cover all the flow regimes from laminar to fully turbulent flow. Pressure loss experiments are performed for both cross-sectional geometries and friction factor results from experiments and numerical simulations are compared. From the analysis of the friction factor as a function of the Reynolds number, the critical value of the Reynolds number linked to the laminar-to-turbulent transition has been determined. The experimental and numerical critical Reynolds number for all the tested microchannels showed a maximum deviation of less than 12%. These results demonstrate that the transitional turbulence model proposed by Abraham etal. (2008) for internal flows can be extended to microchannels and proficiently employed for the design of micro heat exchangers in presence of gas flows.     

Conjugate forced convection and heat conduction in circular microchannels

The effects of axial heat conduction in the solid walls of microchannels of circular cross-sections are analyzed here. A systematic approach is adopted, with the aim of pointing out the influence of geometrical parameters and of solid wall thermal conductivity on microchannel heat transfer. The reliability of a commonly adopted criterium, based on the so-called axial conduction number, to assess the relevance of axial heat conduction is also discussed. Numerical simulations concern the simultaneously developing laminar flow of a constant property fluid in microchannels of different length, wall thickness and wall material, heated with a uniform heat flux at the outer surface, for different values of the Reynolds number. Moreover, since often in experimental tests the two end sections of the microchannel wall are not perfectly insulated, the effects of heat losses through these sections are also considered. A hybrid finite element procedure, which implies the step-by-step solution of the parabolized momentum equations in the fluid domain, followed by the solution of the energy equation in the entire domain, corresponding to both the solid and the fluid parts, is used for the numerical simulations.     

Axial static pressure measurements of water flow in a rectangular microchannel

An experimental investigation of water flow through an aluminum rectangular microchannel with a hydraulic diameter of 169 μm was conducted over a Reynolds number (based upon mean velocity and hydraulic diameter) range from 230 to 4,740. Pressure measurements were simultaneously acquired at eight different axial locations within the channel along with pressure measurements in the inlet and outlet ports. The 27 μm pressure taps were more densely packed near the channel entrance in order to study the developing flow region. The average Poiseuille number for laminar flows was 86.4, which is in excellent agreement with the theoretical value of 86.9. The average critical Reynolds number was found to be 2,370. The limited turbulent friction factor data were in good agreement with the Haaland equation. The inlet to the channel was not well rounded and pressure distributions near the channel entrance show a region of pressure recovery. Entrance length and some minor loss coefficient data were not in agreement with theory, but the cause of these deviations were primarily a function of the inlet geometry and pressure recovery in the microchannel rather than a microscale effect.     

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.     

Rotating electroosmotic flow of an Eyring fluid

A perturbation analysis is presented in this paper for the electroosmotic(EO) flow of an Eyring fluid through a wide rectangular microchannel that rotates about an axis perpendicular to its own. Mildly shear-thinning rheology is assumed such that at the leading order the problem reduces to that of Newtonian EO flow in a rotating channel, while the shear thinning effect shows up in a higher-order problem.Using the relaxation time as the small ordering parameter,analytical solutions are deduced for the leading-as well as first-order problems in terms of the dimensionless Debye and rotation parameters. The velocity profiles of the Ekman–electric double layer(EDL) layer, which is the boundary layer that arises when the Ekman layer and the EDL are comparably thin, are also deduced for an Eyring fluid. It is shown that the present perturbation model can yield results that are close to the exact solutions even when the ordering parameter is as large as order unity. By this order of the relaxation time parameter, the enhancing effect on the rotating EO flow due to shear-thinning Eyring rheology can be significant.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

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.     

Free-Molecular Gas Flow in a Channel with Curving Boundary

Free-molecular gas flow through a microchannel with moving walls curved in accordance with the wave law is simulated numerically. It is shown that the probability of passage of the gas molecules through such a channel depends significantly on the dimensionless ratio of the channel wall wave velocity and the characteristic thermal velocity of the gas molecules. It is revealed that the probabilities of passage are also significantly different when the gas flows “along with” and “against” the direction of wave propagation on the boundary. Applications of this effect to both creating microseparating devices and designing micropumps are discussed. The effect of the problem parameters on the efficiency of these devices is investigated.     

Heat transfer in circular microchannels during volumetric heating with magnetic field

Convective heat transfer within circular microchannels in a rectangular solid substrate with heat generation due to imposed magnetic field was studied. A detailed parametric study was performed by varying Reynolds number, magnetic field strength, working fluid, and the diameter of the channel. It was found that the heat transfer coefficient decreases downstream along the channel. Nusselt number increased with Reynolds number. The tube diameter, properties of the working fluid, and magnetic field strength affected the temperature distribution and heat transfer rate at the solid-fluid interface.     

平板流道交变流动速度环状效应分析

为了分析管内交变流动速度环状效应的特性, 通过数值模拟对平板流道内可压缩交变流动的速度环状效应进行研究. 根据速度环状效应的产生机理, 分析得出在黏性流体管内交变流动中必然存在速度环状效应, 并且在流道截面中心速度为零相位下最容易观察到速度环状效应. 为了定量描述速度环状效应, 根据流道中心速度为零相位下的速度分布曲线的斜率情况, 提出定量评价指标速度环状效应系数, 并针对平行平板流道内层流交变流动,利用该评价指标定量分析了无量纲参数(包括瓦伦西数Va和最大雷诺数Re_max)速度环状效应的影响.      

Pulsatile magneto-biofluid flow and mass transfer in a non-Darcian porous medium channel

A numerical study of pulsatile flow and mass transfer of an electrically conducting Newtonian biofluid via a channel containing porous medium is considered. The conservation equations are transformed and solved under boundary conditions prescribed at both walls of the channel, using a finite element method with two-noded line elements. The influence of magnetic field on the flow is studied using the dimensionless hydromagnetic number, Nm, which defines the ratio of magnetic (Lorentz) retarding force to the viscous hydrodynamic force. A Darcian linear impedance for low Reynolds numbers is incorporated in the transformed momentum equation and a second order drag force term for inertial (Forchheimer) effects. Velocity and concentration profiles across the channel width are plotted for various values of the Reynolds number (Re), Darcy parameter (λ), Forchheimer parameter (Nf), hydro-magnetic number (Nm), Schmidt number (Sc) and also with dimensionless time (T). Profiles of velocity varying in space and time are also provided. The conduit considered is rigid with a pulsatile pressure applied via an appropriate pressure gradient term. Increasing the hydromagnetic number (Nm) from 1 to 15 considerably depresses biofluid velocity (U) indicating that a magnetic field can be used as a flow control mechanism in, for example, medical applications. A rise in Nf from 1 to 20 strongly retards the flow development and decreases the velocity, U, across the width of the channel. The effects of other parameters on the flowfield are also discussed at length. The flow model also has applications in the analysis of electrically conducting haemotological fluids flowing through filtration media, diffusion of drug species in pharmaceutical hydromechanics, and also in general fluid dynamics of pulsatile systems.     

Effects of microchannel geometry on pulsed flow mixing

Although the mixing of reagents is often crucial in many microfluidic devices, good mixing in these laminar, low Reynolds number, flows remains a challenge. It was shown in Refs. [Glasgow, I., Aubry, N., 2003. Lab on a Chip 3, p. 114; Glasgow, I., Batton, J., Aubry, N., 2004. Lab on a Chip 4, p. 558] that pulsing can induce mixing at the confluence of two inlet microchannels in an efficient manner. In this paper, we show that this mixing is affected by both the geometry of the confluence and the inclusion of features in the channels, which induce secondary flow. More specifically, we study mixing in 200 μm wide by 120 μm deep channels, at flow rates from 48 nl s⁻¹ to 4.8 μl s⁻¹, corresponding to Reynolds numbers of 0.3–30. For the parameter values studied, the pulsed flow technique is more effective at mixing than the secondary flow induced by the channel geometry features, and combining both methods leads to even better mixing. In addition, pulsing the reagents such that they pass multiple times through the spatial features, which induce secondary flow leads to mixing over shorter distances.     

Interactions between a rectangular cylinder and a free-surface flow

The salient features of the interaction between a free-surface flow and a cylinder of rectangular cross-section are investigated and discussed. Laboratory-scale experiments are performed in a water channel under various flow conditions and elevations of the cylinder above the channel floor. The flow field is characterized on the basis of time-averaged and fluctuating local velocity measurements. Dynamic loadings on the cylinder are measured by two water-insulated dynamometers placed inside the cylinder structure. Starting from frequency and spectral analyses of the force signals, insights on the relationship between force dominant frequencies and the Strouhal number of the vortex shedding phenomenon are provided. Experimental results highlight the strong influence of the asymmetric configuration imposed by the two different boundary conditions (free surface and channel floor) on (i) the mean force coefficients and (ii) the vortex shedding frequencies. We provide an analysis of the nature of the dependence of average force coefficients on relevant dimensionless groups, i.e., the Reynolds number, normalized flow depth and cylinder submersion.     

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.     

Laminar flow and heat transfer in a periodically converging-diverging channel

A finite difference solution for laminar viscous flow through a sinusoidally curved converging-diverging channel is presented. The physical wavy domain is transformed into a rectangular computational domain in order to simplify the application of boundary conditions on the channel walls. The discretized conservation equations for mass, momentum and energy are derived on a control volume basis. The pseudo-diffusive terms that arise from the co-ordinate transformation are treated as source terms, and the resulting system of equations is solved by a semi-implicit procedure based on line relaxation. Results are obtained for both the developing and the fully developed flow for a Prandtl number of 0.72, channel maximum width-to-pitch ratio of 1.0, Reynolds number ranging from 100 to 500 and wall amplitude-to-pitch ratio varying from 0.1 to 0.25. Results are presented here for constant fluid properties and for a prescribed wall enthalpy only.     

Viscous dissipation effect on heat transfer characteristics of rectangular microchannels under slip flow regime and H1 boundary conditions

Viscous dissipation effect on heat transfer characteristics of a rectangular microchannel is studied. Flow is governed by the Navier–Stokes equations with the slip flow and temperature jump boundary conditions. Integral transform technique is applied to derive the temperature distribution and Nusselt number. The velocity distribution is taken from literature. The solution method is verified for the case where viscous dissipation is neglected. It is found that, the viscous dissipation is negligible for gas flows in microchannels, since the contribution of this effect on Nu number is about 1%. However, this effect should be taken into account for much more viscous flows, such as liquid flows. Neglecting this effect for a flat microchannel with an aspect ratio of 0.1 for Br=0.04 underestimates the Nu number about 5%.     

Theoretical and experimental investigations of the flow of a viscous fluid near the line of intersection of a cylindrical surface and a flat surface

Flow along a corner was investigated at large Reynolds numbers in, for example, [1–3]. The present author [4] considered flow in the neighborhood of a corner formed by the intersection of a plane and a concave cylindrical surface, the main attention being devoted to the formation of the three-dimensional boundary layer on the plane near the corner. It was shown that the curvature of one of the intersecting surfaces changes the flow pattern qualitatively. In the present paper, we report an investigation of the formation of the flow on a concave cylindrical surface near such a corner and consider how the flow is rearranged in the neighborhood of a corner in, for example, a channel of rectangular cross section that has an initial straight section and then a bend with a discontinuity of the curvature of the line of intersection of the concave and flat sides of the channel. The results are given of some experimental investigations of flow near the line of intersection of a flat wall and a curved (concave and convex) wall at a bend in a rectangular channel.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 64–68, January–February, 1983.I thank G. M. Bam-Zelikovich for constant interest in the work and A. I. Ruban for a number of extremely helpful comments.