Incompressible and compressible flows through rectangular microorifices entrenched in silicon microchannels

Incompressible and compressible flows through indispensable configurations such as rectangular microorifices entrenched in microchannels have been experimentally investigated. The current endeavor evaluates the effects of microorifice and microchannel size, estimates the discharge coefficients associated with both compressible and incompressible flows, examines the contraction coefficients, probes subsonic and supercritical gas flows, and explores the presence of any anomalous effects such as those reported for microchannels. The discharge coefficient in incompressible flow, using deionized (DI) water as the working fluid, rises and peaks at a critical Reynolds number, (200/spl les/Re/sub Crit//spl les/500). The reported range of the transitional Reynolds number compares favorably with the values observed in conventional scale studies and suggests the absence of any irregular scaling effects. Furthermore, nitrogen flows through various microorifices suggests that the constriction element rather than the microchannel area determines the flow rate. Additionally, the critical pressure ratio at choking is close to the isentropic value (0.47/spl les/(P/sub 2//P/sub 1/)/sub Crit//spl les/0.64) and no anomalous scale or slip effects have been observed. Unlike macroscale compressible flows through an orifice, the losses seem minimal and the discharge coefficients are close to unity. The geometry acts as a smooth converging-diverging nozzle and the mass flow rate trends appear similar to the data obtained in micronozzle flows.     

Two-Dimensional Analysis of Electrokinetically Driven Out-of-Plane Vortices in a Microchannel Liquid Flow Using Patterned Surface Charge

The technology developed for photolithographically patterning the electric surface charge to be negative, positive, or neutral enables the realization of complex liquid flows even in straight and uniform microchannels with extremely small Reynolds number. A theoretical model to analyze a steady incompressible electrokinetically driven two-dimensional liquid flow in a microchannel with an inhomogeneous surface charge under externally applied electric field is derived. The flow field is obtained analytically by solving the biharmonic equation with the Helmholtz-Smoluchowski slip boundary condition using the Fourier series expansion method. The model has been applied to study three basic out-of-plane vortical flow fields: single vortex and a train of corotating and a series of counterrotating vortex pairs. For model verification, the solution for the single vortex has been tested against numerical computations based on the full Navier-Stokes equations revealing the dominant control parameters. Two interesting phenomena have been observed in out-of-plane multivortex dynamics: merging of corotating vortices and splitting of counterrotating vortices. The criteria for the onset of both phenomena are discussed     

Microscale thermal fluid transport process in a microchannel integrated with arrays of temperature and pressure sensors

One of the most important components in a microfluidic system is the microchannel which involves complicated flow and transport process. This work presents microscale thermal fluid transport process inside a microchannel with a height of 37 μm. The channel can be heated on the bottom wall and is integrated with arrays of pressure and temperature sensors which can be used to measure and determine the local heat transfer and pressure drop. A more simplified model with modification of Young’s Modulus from the experimental test is used to design and fabricate the arrays of pressure sensors. Both the pressure sensors and the channel wall use polymer materials which greatly simplifies the fabrication process. In addition, the polymer materials have a very low thermal conductivity which significantly reduces the heat loss from the channel to the ambient that the local heat transfer can be accurately measured. The airflow in the microchannel can readily become compressible even at a very low Reynolds number condition. Therefore, simultaneous measurement of both the local pressure drop and the temperature on the heated wall are required to determine the local heat transfer. Comparison of the local heat transfer for a compressible airflow in microchannel is made with the theoretical prediction based on incompressible airflow in large scale channel. The comparison has clarified many of the conflicting results among different works.     

Acoustically generated flows in microchannel flexural plate wave sensors: Effects of compressibility

Acoustically generated flowfields in flexural plate wave sensors filled with a Newtonian liquid (water) are considered. A computational model based on compressible flow is developed for the sensor with a moving wall for pumping and mixing applications in microchannels. For the compressible flow formulation, an isothermal equation of state for water is employed. The velocity and pressure profiles for different parameters including flexural wall frequency, channel height, amplitude of the wave and wave length are investigated for four microchannel height/length geometries. It is found that the flowfield becomes pseudo-steady after sufficient number of flexural cycles. Both instantaneous and time averaged results show that an evanescent wave is generated in the microchannel. The predicted flows generated by the FPWs are compared with results available in the literature. The proposed device can be exploited to integrate micropumps with complex microfluidic chips improving the portability of micro-total-analysis systems.     

Uniformly Effective Numerical Methods for Hyperbolic Systems

A unified method to compute compressible and incompressible flows is presented. Accuracy and efficiency do not degrade as the Mach number tends to zero. A staggered scheme solved with a pressure correction method is used. The equation of state is arbitrary. A Riemann problem for the barotropic Euler equations with nonconvex equation of state is solved exactly and numericaly. A hydrodynamic flow with cavitation in which the Mach number varies between 10⁻³ and 20 is computed. Unified methods for compressible and incompressible flows are further discussed for the flow of a perfect gas. The staggered scheme with pressure correction is found to have Mach-uniform accuracy and efficiency, and for the fully compressible case the accuracy is comparable with that of established schemes for compressible flows. Received October 20, 1999; revised May 26, 2000     

Lattice Boltzmann modeling of microchannel flows in the transition flow regime

Owing to its kinetic nature and distinctive computational features, the lattice Boltzmann method for simulating rarefied gas flows has attracted significant research interest in recent years. In this article, a lattice Boltzmann (LB) model is presented to study microchannel flows in the transition flow regime, which have gained much attention because of fundamental scientific issues and technological applications in various micro-electro-mechanical system (MEMS) devices. In the model, a Bosanquet-type effective viscosity is used to account for the rarefaction effect on gas viscosity. To match the introduced effective viscosity and to gain an accurate simulation, a modified second-order slip boundary condition with a new set of slip coefficients is proposed. Numerical investigations demonstrate that the results, including the velocity profile, the non-linear pressure distribution along the channel, and the mass flow rate, are in good agreement with the solution of the linearized Boltzmann equation, the direct simulation Monte Carlo (DSMC) results, and the experimental results over a broad range of Knudsen numbers. It is shown that taking the rarefaction effect on gas viscosity into consideration and employing an appropriate slip boundary condition can lead to a significant improvement in the modeling of rarefied gas flows with moderate Knudsen numbers in the transition flow regime.     

Improved modified Reynolds equation for thin-film gas lubrication from an extended slip velocity boundary condition

An extended slip velocity boundary condition is derived from the regularized 13 moment equations firstly. Different from the existing slip velocity boundary condition, the slip coefficients of the extended one are not fixed, which will change with the wall accommodation coefficient and the Knudsen number of the gas flow. Using the extended slip velocity condition, an improved modified Reynolds equation for thin-film gas lubrication is established. From solving the improved modified Reynolds equation, the pressure distribution of the slider gas bearing is obtained and has a better agreement with that from the direct simulation Monte Carlo method under different pitch angles and wall velocities. It is found that the improved modified Reynolds equation can predict a more accurate pressure distribution of the slider gas bearing than the Fukui and Kaneko’s lubrication model from the linearized Boltzmann equation in the near transition regime.     

Fluid flows in microchannels with cavities

Pressure-driven gas and liquid flows through microchannels with cavities have been studied using both experimental measurements and numerical computations. Several microchannels with cavities varying in shape, number and dimensions have been fabricated. One set of microdevices, integrated with sensors on a silicon wafer, is used for flow rate and pressure distribution measurements in gas flows. Another set of microdevices, fabricated using glass-to-silicon wafer bonding, is utilized for visualization of liquid flow patterns. The cavity effect on the flow in the microchannel is found to be very small, with the mass flow rate increasing slightly with increasing number of cavities. The flow pattern in the cavity depends on two control parameters; it is fully attached only if both the reduced Reynolds number and the cavity number are small. A flow regime map has been constructed, where the critical values for the transition from attached to separated flow are determined. The numerical computations reveal another control parameter, the cavity aspect ratio. The flow in the cavity is similar only if all three control parameters are the same. Finally, the vorticity distribution and related circulation in the cavity are analyzed. [1546].     

Investigations on characteristics of micro/meso scale gas foil journal bearings for 100–200 W class micro power systems using first order slip velocity boundary conditions and the effective viscosity model

In this paper, a modified compressible Reynolds equation for micro/meso scale gas foil journal bearings considering first order slip and effective viscosity under rarefied flow conditions is presented. The influence of rarefaction effect on the load carrying capacity, attitude angle, speed and frequency dependent stiffness and damping coefficients, modal impedance, natural frequencies and unbalance response is studied. From numerical analysis, it has been found that there is significant change in all the static and dynamic characteristics predicted by the no-slip model and model with effective viscosity. There is also a considerable difference between the values predicted by a model with effective viscosity and a model without effective viscosity. For a given eccentricity ratio, the influence of effective viscosity on load carrying capacity and attitude angle is more significant for the typical operating speed range of micro/meso scale gas turbines. The influence of effective viscosity decreases with increase in compliance of bearing structure. Similarly, the influence of effective viscosity on frequency dependent stiffness and damping coefficients increases with excitation frequency ratio. Significant change in natural frequency, modal impedance and unbalance response for model with no slip and slip with effective viscosity is observed. The influence of effective viscosity is found to be significant with increase in Knudsen number.     

Liquid and gas flows in microchannels of varying cross section: a comparative analysis of the flow dynamics and design perspectives

This paper presents a comparative study of the flow of liquid and gases in microchannels of converging and diverging cross sections. Towards this, the static pressure across the microchannels is measured for different flow rates of the two fluids. The study includes both experimental and numerical investigations, thus providing several useful insights into the local information of flow parameters as well. Three different microchannels of varying angles of convergence/divergence (4°, 8° and 12°) are studied to understand the effect of the angle on flow properties such as pressure drop, Poiseuille number and diodicity. A comparison of the forces involved in liquid and gas flows shows their relative significance and effect on the flow structure. A diodic effect corresponds to a difference in the flow resistance in a microchannel of varying cross section, when the flow is subjected alternatively to converging and diverging orientations. In the present experiments, the diodic effect is observed for both liquid and gas as working fluids. The effect of governing parameters—Reynolds number and Knudsen number, on the diodicity is analysed. Based on these results, a comparison of design perspectives that may be useful in the design of converging/diverging microchannels for liquid and gas flows is provided.     

Comparative study between computational and experimental results for binary rarefied gas flows through long microchannels

A comparative study between computational and experimental results for pressure-driven binary gas flows through long microchannels is performed. The theoretical formulation is based on the McCormack kinetic model and the computational results are valid in the whole range of the Knudsen number. Diffusion effects are taken into consideration. The experimental work is based on the Constant Volume Method, and the results are in the slip and transition regime. Using both approaches, the molar flow rates of the He–Ar gas mixture flowing through a rectangular microchannel are estimated for a wide range of pressure drops between the upstream and downstream reservoirs and several mixture concentrations varying from pure He to pure Ar. In all cases, a very good agreement is found, within the margins of the introduced modeling and measurement uncertainties. In addition, computational results for the pressure and concentration distributions along the channel are provided. As far as the authors are aware of, this is the first detailed and complete comparative study between theory and experiment for gaseous flows through long microchannels in the case of binary mixtures.     

Analytical solution of gas lubricated slider microbearing

The analytical solution of a two-dimensional, isothermal, compressible gas flow in a slider microbearing is presented. A higher order accuracy of the solution is achieved by applying the boundary condition of Kn ² order for the velocity slip on the wall, together with the momentum equation of the same order (known as the Burnett equation). The analytical solution is obtained by the perturbation analysis. The order of all terms in continuum and momentum equations and in boundary conditions is evaluated by incorporating the exact relation between the Mach, Reynolds and Knudsen numbers in the modelling procedure. Low Mach number flows in microbearing with slowly varying cross-sections are considered, and it is shown that under these conditions the Burnett equation has the same form as the Navier–Stokes equation. Obtained analytical results for pressure distribution, load capacity and velocity field are compared with numerical solutions of the Boltzmann equation and some semi-analytical results, and excellent agreement is achieved. The model presented in this paper is a useful tool for the prediction of flow conditions in the microbearings. Also, its results are the benchmark test for the verifications of various numerical procedures.     

Efficient heat transfer enhancement by elastic turbulence with polymer solution in a curved microchannel

Heat and mass transfer in microscale flows are limited due to extremely low Reynolds number (Re). In a curved microchannel, however, complex flow behaviors, such as elastic instability and elastic turbulence, can be induced via viscoelastic fluid at vanishingly low-Re conditions, which is of great potential to enhance the heat transfer performance. The influence of elastic instabilities and turbulence on heat dissipation of exothermic components is experimentally investigated in this study. The heat transfer performance of both viscoelastic (polymer solutions) and Newtonian (sucrose solutions) fluid flows in a curved microchannel with a square cross section is experimentally characterized. Titanium–platinum (Ti–Pt) thin films embedded at the bottom wall of the polydimethylsiloxane (PDMS) microchannel serve as both microheater and temperature sensor. For viscoelastic fluids, the spectrum of outlet temperature fluctuation in broad frequency (f) region fits the power law of f ^?1.1. Heat transfer enhancement due to the elastic turbulence in a curved microchannel is thereby identified by the drastic growth of the Nusselt number (Nu, the ratio of convective to conductive heat transfer normal to the boundary) with the increase in the Weissenberg number (Wi, the ratio of elastic stress to viscous stress). The mechanism of heat transfer enhanced by the convection effect of elastic turbulence is also elucidated.     

Numerical simulation of electroosmotic effect in serpentine channels

Numerical simulation of flow through a three-dimensional serpentine microchannel, subjected to a voltage perpendicular to the flow direction is presented here. Commercial CFD software CFD-ACE+ is used for the numerical analysis. A parametric study is conducted to investigate the effect of radius of curvature, Reynolds number, zeta potential and Debye length on the pressure drop and friction factor. Each case is compared with flow without electroosmotic effect. It is found that electroosmosis induces secondary flow patterns in the straight portion of the channel in addition to secondary vortices at bends. This electroosmosis-induced secondary flow causes additional pressure drop as compared to flow without electroosmosis. For flow with and without electroosmosis, the pressure drop increases with Reynolds number and the nature of variation is qualitatively similar in both the cases. It is also found that the pressure drop increases as the Debye length is reduced.     

Second-order gaseous slip flow models in long circular and noncircular microchannels and nanochannels

This paper significantly extends previous studies to the transition regime by employing the second-order slip boundary conditions. A simple analytical model with second-order slip boundary conditions for a normalized Poiseuille number is proposed. The model can be applied to either rarefied gas flows or apparent liquid slip flows. The developed simple models can be used to predict the Poiseuille number, mass flow rate, tangential momentum accommodation coefficient, pressure distribution of gaseous flow in noncircular microchannels and nanochannels by the research community for the practical engineering design of microchannels and nanochannels. The developed second-order models are preferable since the difficulty and “investment” is negligible compared with the cost of alternative methods such as molecular simulations or solutions of Boltzmann equation. Navier–Stokes equations with second-order slip models can be used to predict quantities of engineering interest such as the Poiseuille number, tangential momentum accommodation coefficient, mass flow rate, pressure distribution, and pressure drop beyond its typically acknowledged limit of application. The appropriate or effective second-order slip coefficients include the contribution of the Knudsen layers in order to capture the complete solution of the Boltzmann equation for the Poiseuille number, mass flow rate, and pressure distribution. It could be reasonable that various researchers proposed different second-order slip coefficients because the values are naturally different in different Knudsen number regimes. It is analytically shown that the Knudsen’s minimum can be predicted with the second-order model and the Knudsen value of the occurrence of Knudsen’s minimum depends on inlet and outlet pressure ratio. The compressibility and rarefaction effects on mass flow rate and the curvature of the pressure distribution by employing first-order and second-order slip flow models are analyzed and compared. The condition of linear pressure distribution is given.     

Experiment and simulation of mixed flows in a trapezoidal microchannel

This paper presents experimental and numerical results of mixed electroosmotic and pressure driven flows in a trapezoidal shaped microchannel. A micro particle image velocimetry (μPIV) technique is utilized to acquire velocity profiles across the microchannel for pressure, electroosmotic and mixed electroosmotic-pressure driven flows. In mixed flow studies, both favorable and adverse pressure gradient cases are considered. Flow results obtained from the μPIV technique are compared with 3D numerical predictions, and an excellent agreement is obtained between them. In the numerical technique, the electric double layer is not resolved to avoid expensive computation, rather a slip velocity is assigned at the channel surface based on the electric field and electroosmotic mobility. This study shows that a trapezoidal microchannel provides a tapered-cosine velocity profile if there is any pressure gradient in the flow direction. This result is significantly different from that observed in rectangular microchannels. Our experimental results verify that velocity distribution in mixed flow can be decomposed into pressure and electroosmotic driven components.     

An immersed-boundary method for compressible viscous flows

This paper combines a state-of-the-art method for solving the preconditioned compressible Navier-Stokes equations accurately and efficiently for a wide range of the Mach number with an immersed-boundary approach which allows one to use Cartesian grids for arbitrarily complex geometries. The method is validated versus well documented test problems for a wide range of the Reynolds and Mach numbers. The numerical results demonstrate the efficiency and versatility of the proposed approach as well as its accuracy, from incompressible to supersonic flow conditions, for moderate values of the Reynolds number. Further improvements, obtained via local grid refinement or non-linear wall functions, can render the proposed approach a formidable tool for studying complex three-dimensional flows of industrial interest.     

Microchannel Pressure Measurements Using Molecular Sensors

Fluid mechanics on the microscale is an important subject for researchers who are interested in studying microdevices since physical phenomena change from macroscale to microscale. Channel flow is a fundamental topic for fluid mechanics. By using a molecular sensor known as pressure-sensitive paint (PSP), detailed pressure data can be obtained inside the microchannel and at the channel entrance. The achievable spatial resolution of the acquired pressure map can be as high as 5 mum. PSP measurements are obtained for various pressure ratios from 1.76 to 20, with Knudsen number (K n) varying from 0.003 to 0.4. Compressibility and rarefaction effects can be seen in the pressure data inside the microchannel and at the channel entrance.     

Demarcating wetting states in textured microchannels under flow conditions by Poiseuille number

In this work, flow friction in microchannels decorated with micropillars was investigated experimentally, with an interest to understand the wetting transition through two simple means: Poiseuille number and scaling laws. Different wetting states were demarcated by qualitatively assessing the behaviour of Poiseuille number (Po = f·Re, where f is friction factor and Re is Reynolds number), which are further corroborated by confocal microscopy-based measurements and numerical simulations. The wetting transition ensued smoothly with an increase in Re, independent of the gas fraction (a ratio of area covered by the liquid–gas interface to the total projected area), for moderate gas fractions, whereas an early breakdown of the Cassie–Baxter state occurred irrespective of Re at high gas fractions. Additionally, the scaling laws were found to correlate well with the underlying state of the flow. Our observations revealed that the liquid–gas interface exhibits a partial slip, contrary to the common notion that it is shear free. It is inferred that an increase in effective flow area leads to a reduction in flow friction in textured microchannels. The present work underlines three important outcomes. The first is the identification of wetting states in flow conditions shown by tracking the Poiseuille number. The second is that the liquid–gas interface is deduced to behave like a partial slip boundary. The third is that a textured microchannel can be worse than an enlarged dimension microchannel.     

Computational method for turbulent supersonic flows

The calculation features of the turbulent flows described by the Reynolds equations and the two-equation model of turbulence are examined for an explicit high-order accurate Godunov method. Under these features, a new version of a high order of Godunov’s method is developed for calculating the compressible turbulent flows. To illustrate the capability of the new method, some results of the calculation are shown for a supersonic turbulent jet with a complex shock-wave structure and for a separate flow in a plane nozzle.