Flow In Microchannel Research Articles

Mixing Performance Analysis and Optimal Design of a Novel Passive Baffle Micromixer.

Micromixers, as crucial components of microfluidic devices, find widespread applications in the field of biochemistry. Due to the laminar flow in microchannels, mixing is challenging, and it significantly impacts the efficiency of rapid reactions. In this study, numerical simulations of four baffle micromixer structures were carried out at different Reynolds numbers (Re = 0.1, Re = 1, Re = 10, and Re = 100) in order to investigate the flow characteristics and mixing mechanism under different structures and optimize the micromixer by varying the vertical displacement of the baffle, the rotation angle, the horizontal spacing, and the number of baffle, and by taking into account the mixing intensity and pressure drop. The results indicated that the optimal mixing efficiency was achieved when the baffle's vertical displacement was 90 μm, the baffle angle was 60°, the horizontal spacing was 130 μm, and there were 20 sets of baffles. At Re = 0.1, the mixing efficiency reached 99.4%, and, as Re increased, the mixing efficiency showed a trend of, first, decreasing and then increasing. At Re = 100, the mixing efficiency was 97.2%. Through simulation analysis of the mixing process, the structure of the baffle-type micromixer was effectively improved, contributing to enhanced fluid mixing efficiency and reaction speed.

Flow pattern-based heat transfer analysis of microchannel flow boiling: An experimental investigation

The present study investigates the flow boiling heat transfer and flow pattern of R245fa in a smooth, horizontal microchannel with an inner diameter of 1 mm. The experiments cover a range of mass flux from 300 to 700 kg/m2s and heat fluxes varying between 10 and 94 kW/m2, with corresponding saturation temperature of 23–54 °C. Heat transfer experiment results indicate that the local heat transfer coefficient (HTC) presents two distinct variation trends under different heat and mass flux conditions. In addition, the heat flux is positively correlated with the local HTC at each measure point. Visualization experiment results show that annular flow forms at low vapor quality and significantly influences heat transfer as the primary flow pattern. The flow boiling heat transfer of R245fa in the 1 mm microchannel can be characterized by three regimes: bubble/slug regime, liquid film evaporation regime, and intermittent dryout regime. Thin liquid film evaporation and intermittent dryout during annular flow predominately contribute to heat transfer. Finally, the experimental HTC is compared with five different predictive methods, revealing that the three-zone model proposed by Thome et al. [46] holds immense potential for estimating the local HTC in the microchannel.

Influence of randomly distributed wall roughness on the flow behaviors of liquid-liquid slug flow in a co-flowing microchannel

In present study, the influence of randomly distributed wall roughness on the flow behaviors of liquid-liquid slug flow in a co-flowing microchannel has been numerically investigated, as the manufacturing tolerance of microfluidic device in engineering application is unavoidable. A circular axisymmetric microchannel model is adopted to study the effect of roughness height and correlation length on the wall roughness generation and internal laminar flow characteristics in the microchannel. The wall roughness is modeled by well-known Gaussian distribution method. The prediction of slug velocity and liquid film thickness from numerical simulations has been well validated with the empirical and analytical correlations and experimental data in the literature. The obtained results show that the rough wall condition significantly affected the internal flow behaviors, such as slug moving velocity and slug breaking-off. According to the quantitative comparison of normalized surface areas between high (larger than 0.4 m/s) and low (smaller than 0.69 m/s) velocity region, the surface area of low velocity region under rough wall condition (wall roughness: D/l=3.333 and σ/D = 0.018) is 46 % larger than smooth wall condition, which means effective flow area of microchannel is reduced. The surface area of high velocity region under rough wall condition increases around 66 % compared with smooth wall condition. This study could be used as a promising reference and provide more understanding for roughness control in industrial microfluidic device.

Magnetic Microbead-Based Herringbone Chip for Sensitive Detection of Human Immunodeficiency Virus.

The microfluidic chip-based nucleic acid detection method significantly improves the sensitivity since it precisely controls the microfluidic flow in microchannels. Nonetheless, significant challenges still exist in improving the detection efficiency to meet the demand for rapid detection of trace substances. This work provides a novel magnetic herringbone (M-HB) structure in a microfluidic chip, and its advantage in rapid and sensitive detection is verified by taking complementary DNA (cDNA) sequences of human immunodeficiency virus (HIV) detection as an example. The M-HB structure is designed based on controlling the magnetic field distribution in the micrometer scale and is formed by accumulation of magnetic microbeads (MMBs). Hence, M-HB is similar to a nanopore microstructure, which has a higher contact area and probe density. All of the above is conducive to improving sensitivity in microfluidic chips. The M-HB chip is stable and easy to form, which can linearly detect cDNA sequences of HIV quantitatively ranging from 1 to 20 nM with a detection limit of 0.073 nM. Compared to the traditional herringbone structure, this structure is easier to form and release by controlling the magnetic field, which is flexible and helps in further study. Results show that this chip can sensitively detect the cDNA sequences of HIV in blood samples, demonstrating that it is a powerful platform to rapidly and sensitively detect multiple nucleic acid-related viruses of infectious diseases.

Understanding Inconsistencies in Thermohydraulic Characteristics Between Experimental and Numerical Data for DI Water Flow Through a Rectangular Microchannel

Abstract Facing discrepancies between numerical simulation, experimental measurement, and theory is common in studies of fluid flow and heat transfer in microchannels. The cause of these discrepancies is often linked to the transition from the macroscale to the microscale, where the flow dynamics might be expected to deviate due to possible changes in dominant forces. In this work, an attempt is made to achieve agreement between experiment, numerical simulation, and theoretical description within the usual framework of laminar flow theory. For this purpose, the pressure drop, friction factor, and Poiseuille number under isothermal conditions and the temperature profile, heat transfer coefficient, Nusselt number, and thermal performance index under diabatic conditions (heating power of 10 W) in a heat sink with a stainless steel microchannel with a hydraulic diameter of 850 μm were investigated numerically and experimentally for mass flow rates between 1 and 68 gmin−1. The source of inconsistencies in pressure drop characteristics is found to be linked to the geometrical details of the utilized microchannel, for example, the design of inlet/outlet manifolds, the artifacts of manufacturing technique, and other features of the experimental test rig. For the heat transfer characteristics, it is identified that an appropriate estimation of the outer boundary condition for the numerical simulation remains the crucial challenge to obtain a reasonable agreement. The paper provides a detailed overview of how to account for these details to mitigate the discrepancies and to establish a handshake between experiments, numerical simulations, and theory.

Improved flow boiling performance and temperature uniformity in counter-flow interconnected microchannel heat sink

Microchannel flow boiling is one of the most effective solutions to the problem of heat dissipation caused by high power electronic devices. However, inherent drawbacks such as low critical heat flux due to dry-out of the liquid film near the channel exit, relatively high two-phase flow pressure drop, severe temperature nonuniformity and boiling instability at high heat flux still prevent widespread commercial application. In this study, based on the idea of uniform vapor quality along the flow direction, a copper Counter-Flow Interconnected Microchannel (CFIM) was proposed to address the main issues currently faced by two-phase microchannel heat sink. A comprehensive comparative study of CFIM, Co-Current Microchannel (CCM) and Counter-Flow Microchannel (CFM) was carried out from the aspects of boiling stability, boiling heat transfer characteristics and temperature uniformity. Three slots of 0.2 mm, 0.5 mm and 0.8 mm in width, denoted as IM0.2, IM0.5 and IM0.8, respectively, are considered. Experiments are conducted at mass fluxes of 291–618 kg/m2·s, saturation temperature of 40 °C and effective heat fluxes of 3–293 W/cm2, using the dielectric fluid R1233zd(E) as the working fluid. The results show that the CHF and average two-phase heat transfer coefficient (HTC) of CFIM are increased by 40.6 ∼ 100.5 % and 70.4–83.6 %, respectively, while the two-phase pressure drop is decreased by 17–32.6 % compared to CCM. The design of the counter-flow interconnected structure can manipulate the void fraction and two-phase flow pattern, resulting in a near uniform void fraction and flow pattern. The structure of these connecting slots promotes nucleate boiling and fluid mixing of neighboring channels. More importantly, the flow boiling instability and temperature non-uniformity in CFIM are well suppressed. In addition, IM0.2 has higher heat dissipation than IM0.5 and IM0.8 due to enhanced bubble nucleation, reducing superheat requirements and periodic rewetting to the center of the channel. The results of this research reveal that utilizing a microchannel heat sink with a counter-flow interconnected configuration can effectively improve heat transfer performance during flow boiling. Furthermore, it helps to mitigate issues such as flow boiling instability and temperature non-uniformity.

Non-Newtonian Characteristics of Water Flow in Microchannels

Non-Newtonian Characteristics of Water Flow in Microchannels

A Comparative Study of Two Numerical Methods Applied for 3D Liquid-Liquid Taylor Flow in a Microchannel

The laminar nature of flow in mini and microchannels has pushed researchers to develop novel solutions to overcome reaction rate reduction and heat/mass transfer issues. In this regard, Taylor flow is one of the possible solutions that could be used to enhance mixing inside mini and microchannels with reasonable pressure drop. The hydrodynamics of Taylor liquid-liquid flow is numerically studied in this work by employing two different droplet generation methods, specifically T-junction and patching methods. To this end, a three-dimensional model of rectangular microchannel flow is considered. The computational domain was designed and meshed by ICEM CFD and then simulated with commercial software ANSYS Fluent. The interface between the two phases was captured using the Volume of Fluid (VOF) method. The generating and development process of water droplets dispersed in an ethylene/propylene glycol carrier phase for both methods is discussed in detail. According to the results, both methods show satisfactory performance regarding liquid film and droplet shape, with only a slight difference. However, the patching method was found to be more economical in terms of computational time. This study would improve the state of knowledge on two-phase flow simulation in microchannels and thus contribute to the understanding of Taylor flow hydrodynamics.

MODELING OF TWO-PHASE GAS-LIQUID SLUG FLOWS IN MICROCHANNELS

Even though numerical simulation of two-phase flows in microchannels has been attempted by many investigators, most efforts seem to have failed in correctly capturing the flow physics, especially the slug-flow regime characteristics. The presence of a thin liquid film in the order of 10 μm around the bubble (sometimes called gas pocket or gas slug) may be a contributing factor to the above difficulty. Typically, liquid films have a significant effect on the flow field. Thus, there is a strong motivation to employ numerical simulation methods to avoid some of the experimental difficulties. In this paper, the characteristics of two-phase slug flows in microchannels are calculated with the help of the volume-of-fluid (VOF) method. Formation of the slugs for different superficial velocities, capillary numbers, and gas volume fractions are investigated. The minimum mesh resolution required to capture the liquid film surrounding the gas bubble is reported employing a dynamic mesh adaptation methodology with interface tracking. Results are shown to be in good agreement with experimental data and empirical correlations.

Cross-streamline migration and near-wall depletion of elastic fibers in micro-channel flows.

The complex dynamics of elastic fibers in viscous fluids are central to many biological and industrial systems. Fluid-structure interactions underlying these dynamics govern the shape and transport of flexible fibers, and understanding these interactions can help tune flow properties in applications such as microfluidic separation, printing and clogging. In this work, we use slender-body theory to study micromechanical dynamics that arise from the coupling between the elastic backbone of a fiber and the local straining flow that contributes to filament flipping and cross-streamline migration. The resulting transverse drift is unbiased in either direction in simple shear flow. However, a non-uniform shear rate results in bias towards regions of high shear, which we connect to the shape transitions during flips. We discover a depletion layer that forms near the boundaries of pressure-driven channel flow due to the competition between such a cross-streamline drift and steric exclusion from the walls. Finally, we develop scaling laws for the curvature of filaments during flip events, demonstrating the origin of the drift bias in non-uniform flows, and confirm this behavior from our simulations. Put together, these results shed light on the role of a local and dominant coupling between elasticity and viscous resistance in dictating long-term dynamics and transport of elastic fibers in confined flows.

Upstream mobility and swarming of light activated micromotors.

Micromotors have been proposed for applications such as targeted drug delivery, thrombolysis, or sensing. However, single micrormotors are limited in the amount of payload they can deliver or force they can exert. Swarms of micromotors can overcome many of these challenges, however creating and controlling such swarms presents many challenges of its own. In particular, utilizing swarms in fluid flows is of significant importance for biomedical or lab-on-chip applications. Here, the upstream mobility and swarm formation of light driven micromotors in microchannel flows is demonstrated with maximum speeds around 0.1 mm s-1. Additionally, the light actuated microrobots operate in fairly low concentrations of hydrogen peroxide of approximately 1%. The micromotors form swarms at the boundary of the locally applied light pattern and the swarms can be moved by translating the light up or downstream.

Effect of parameters-dependent viscosity and thermal conductivity on the thermo-hydraulic performance of Al2O3-based nanofluids in a rectangular microchannel

The multi-parameter dependence of nanofluid viscosity and thermal conductivity on the flow and heat transmission properties of Al2O3-H2O nanofluidin microchannels is investigated in this paper using a numerical simulation approach, exploring the effects of various parameters on flow and heat dissipation characteristics, such as flow velocity, pressure drop, temperature, and heat transfer coefficient, under different Reynolds numbers, etc. These parameters include the nanoparticles phericity (0.5~1.0), volume fraction (0.6%~6.0%), and temperature (290K~370K) of the nanofluids, as well as the boundary parameter Reynolds number (200~1000). The results show that the multi-parameter dependence of nanofluid viscosity and thermal conductivity setting has a significant impact on the flow and heat dissipation characteristics of nanofluids, comparing analysis with the nanofluid viscosity and thermal conductivity only related to the volume fraction under the same condition. Under the combined action of sphericity, volume fraction, and temperature, increasing the sphericity increases the pressure drop and decreases the heat transfer coefficient. For example, at a Reynolds number of 1000, the maximum rates of change for pressure drop and heat transfer coefficient are 0.85 and 4.26, respectively. The nanofluid thermo-hydraulic performance is sensitive to temperature, volume fraction and sphericity in turn. Setting up viscosity and thermal conductivity equations with multiple parameter dependencies can provide more accurate results for there search of nanofluids, further deepening the application research of nanofluids.

The adoption of novel cooling liquids for enhancing heat transfer performance of three-dimensional integrated circuits with embedded micro-channel

The nanofluids (including MWCNT based nanofluid and SWCNT based nanofluid) and liquid metal Ga68In20Sn12 are proposed to replace the conventional water as cooling liquid of micro-channel for enhancing heat transfer performance of three-dimensional integrated circuits (3-D ICs) in this paper. An equivalent thermal model of 3-D ICs with integrated micro-channel is established to investigate the heat transfer performances for using different cooling liquids. The results show that the steady-state temperature for MWCNT based nanofluid, SWCNT based nanofluid and Ga68In20Sn12 as cooling liquids can be reduced over 25.698%, 28.771% and 35.735% than the conventional water scheme in a four-layers stacked chip, respectively. Besides, it is found that the steady-state temperature of all die layer in 3-D ICs can be further reduced by increasing the micro-channel size and flow velocity of cooling liquid. Therefore, the proposed novel materials (i.e., nanofluids and Ga68In20Sn12) as cooling liquids have excellent application prospect in solving thermal problems of 3-D ICs.

Laser-induced microbubble as an in vivo valve for optofluidic manipulation in living Mice's microvessels.

Optofluidic regulation of blood microflow in vivo represents a significant method for investigating illnesses linked to abnormal changes in blood circulation. Currently, non-invasive strategies are limited to regulation within capillaries of approximately 10 μm in diameter because the adaption to blood pressure levels in the order of several hundred pascals poses a significant challenge in larger microvessels. In this study, using laser-induced microbubble formation within microvessels of the mouse auricle, we regulate blood microflow in small vessels with diameters in the tens of micrometers. By controlling the laser power, we can control the growth and stability of microbubbles in vivo. This controlled approach enables the achievement of prolonged ischemia and subsequent reperfusion of blood flow, and it can also regulate the microbubbles to function as micro-pumps for reverse blood pumping. Furthermore, by controlling the microbubble, narrow microflow channels can be formed between the microbubbles and microvessels for assessing the apparent viscosity of leukocytes, which is 76.9 ± 11.8 Pa·s in the in vivo blood environment. The proposed design of in vivo microbubble valves opens new avenues for constructing real-time blood regulation and exploring cellular mechanics within living organisms.

A normal-mode approach for high-speed rarefied plane Couette flow

Based on gas kinetic theory, a linear stability analysis method for low-speed rarefied flows was developed by Zou et al. [“A new linear stability analysis approach for microchannel flow based on the Boltzmann Bhatnagar–Gross–Krook equation,” Phys. Fluids 34, 124114 (2022) and “A novel linear stability analysis method for plane Couette flow considering rarefaction effects,” J. Fluid Mech. 963, A33 (2023)]. In the present study, we extended the method to high-speed rarefied flows using the Bhatnagar–Gross–Krook model. The Chebyshev spectral method is employed to discretize physical space, and the Gauss–Hermite and fourth-order Newton–Cotes quadrature methods are used to discretize velocity space. The fourth-order Newton–Cotes quadrature method was found to have sufficient accuracy for the stability analysis, laying the foundation for future research on hypersonic flows. The stability analysis of compressible rarefied Couette flow showed that acoustic modes are reflected between the wall and the relative sonic line, and the variation in their phase speed and growth rate with the wavenumber is not affected by the Mach number (Ma) and the Knudsen number (Kn). Increasing Kn has a stabilizing effect on both the acoustic and viscous modes, but as Ma increases, the attenuation rate of each mode's growth rate gradually decreases. In subsonic and sonic flows, the least stable viscous mode dominates in the case of small numbers. As Kn increases, the viscous mode gradually dominates over all wavenumber ranges considered in subsonic flow. In sonic flow, mode 1 is dominant in the region beyond the range of small wavenumbers. In supersonic flow, mode 2 is the least stable in the large wavenumber ranges, while mode 1 is the least stable in other wavenumber ranges. At a fixed wavenumber, as Kn increases, the decay rate of the growth rate of mode 2 is the highest. Additionally, under different Knudsen numbers, the growth rates of mode 1, mode 2, and the least stable viscous mode monotonically increase with an increase in Ma, with mode 2 showing the most significant increase.

Dean vortex-enhanced blood plasma separation in self-driven spiral microchannel flow with cross-flow microfilters

Point-of-care (POC) diagnostic devices have been developing rapidly in recent years, but they are mainly using saliva instead of blood as a test sample. A highly efficient self-separation during the self-driven flow without power systems is desired for expanding the point-of-care diagnostic devices. Microfiltration stands out as a promising technique for blood plasma separation but faces limitations due to blood cell clogging, resulting in reduced separation speed and efficiency. These limitations are mainly caused by the high viscosity and hematocrit in the blood flow. A small increment in the hematocrit of the blood significantly increases the pressure needed for the blood plasma separation in the micro-filters and decreases the separation speed and efficiency. Addressing this challenge, this study explores the feasibility of diluting whole blood within a microfluidic device without external power systems. This study implemented a spiral microchannel utilizing the inertial focusing and Dean vortex effects to focus the red blood cells and extract the blood with lower hematocrit. The inertial migration of the particles during the capillary flow was first investigated experimentally; a maximum of 88% of the particles migrated to the bottom and top equilibrium positions in the optimized 350 × 60 μm (cross-sectional area, 5.8 aspect ratio) microchannel. With the optimized dimension of the microchannel, the whole blood samples within the physiological hematocrit range were tested in the experiments, and more than 10% of the hematocrit reduction was compared between the outer branch outlet and inner branch outlet in the 350 × 60 μm microchannel.

Dynamic behavior of DNA molecules in microchannels: exploring deflective, elliptical, and spin motions induced by Saffman and Magnus forces.

Precise manipulation of individual DNA molecules entering and leaving the channel ports, as well as their smooth passage across the channel, is essential for the detection and screening of DNA molecules using nano-/micro-fluidic technologies. In this paper, by combining single-molecule fluorescence imaging and numerical simulations, the motion states of DNA molecules translocating through a microfluidic channel under the action of the applied electric field are monitored and analyzed in detail. It is found that, under certain conditions of the applied electric field DNA molecules exhibit various motion states, including translation crossing, deflection outflow, reverse outflow, reciprocal movement, and elliptical movement. Simulations indicate that, under the action of Saffman force, DNA molecules can only undergo deflective motion when they experience a velocity gradient in the microchannel flow field; and they can only undergo elliptical motion when their deflective motion is accompanied by a spin motion. In this case, the Magnus force also plays an important role. The detailed study and elucidation of the movement states, dynamic characteristics and mechanisms of DNA molecules such as the deflective and elliptical motions under the actions of Saffman and Magnus forces have helpful implications for the development of related DNA/gene nano-/microfluidic chips, and for the separation, screening and detection of DNA molecules.

Numerical model development for critical region flow and thermodynamic transitions of CO2 fluids in microchannel

Dramatic changes of thermophysical properties in the vicinity of critical point brings significant accelerating effect on the thermodynamic transitions of fluid. This study performs a confined geometric design with CO2 fluids passing through the critical zone in microchannel. Such a problem formulation combines the compressible Navier-Stokes equations with different property formulations to obtain critical anomaly effects at complex critical transport coefficients. Three equations of state (van der Waals (vdW), Redlich-Kwong (R-K) and Peng-Robinson (P-R) state equations) and their corresponding thermodynamic relations are derived in solving the governing equations of flow and heat transfer equations. Basic results are found to be consistent with theoretical asymptotic analysis, with linearly curves of pressure and a sudden shift of density (425 kg/m3–502 kg/m3) and velocity (0.82 m/s to 0.97 m/s) profiles when the microchannel flow across the critical region. The cascade effect between weakened radial thermal conductivity (λ ∼ (ρ – 1)-1/2, cp ∼ (ρ - 1)−2) and enhanced fluid expansion (βp ∼ (ρ - 1)−2) in the critical region is also shown to cause a sharp negative transition in temperature. In addition, the simulations demonstrate a qualitative comparison for three state equations, and the distinctions in quantitative terms (negative temperature value are −0.1 (vdW), −0.17 (R-K), and − 0.23(P-R)) are found to originate from stronger divergence characteristics (the specific heat cp and thermal expansion coefficients βp) in P-R gases than vdW and R-K gases.

Influence of convective heat transfer and wall thermal capacity on dynamic interactions between wall temperature and pressure drop oscillations during microchannel flow boiling

Flow boiling in microchannel heat sinks is attractive for use in industrial thermal management applications owing to highly effective heat dissipation and compactness. However, flow boiling is subject to two-phase flow instabilities, specifically pressure drop oscillations, that may deteriorate performance by causing transient wall temperature oscillations or initiating premature critical heat flux. While there has been extensive model-based investigation of the stability and flow dynamics associated with pressure drop oscillations in microchannels, prior work has not focused on the coupled interactions between these flow oscillations and the dynamic heat transfer processes, as is required to predict thermal performance degradation and ensure reliable operation. This study predicts the influence of wall thermal capacity and heat transfer coefficient on the dynamic interaction between pressure drop oscillations and wall temperature. First, we demonstrate and explain how added wall thermal capacity acts to reduce the amplitude of pressure drop oscillations due to coupling between the thermal and flow dynamics. Whereas previous work has attributed this effect to a thermal time constant of the wall, the amplitude of the pressure drop oscillations is shown to scale with the Biot number by scanning various combinations of wall thermal capacity and heat transfer coefficient. A Biot number threshold is thereby suggested as the criterion for judging whether the thermal wall capacity must be considered in dynamic models, based on the strength coupling between the wall temperature and pressure drop oscillations.

Heat Transfer of Air-Water Flow in Wavy Microchannels

Effective cooling of the electronic devices is of great importance to their performance. Of particular interest is when the size of the equipment has scaled down to micro level. Despite a continuous improvement over the years, there is only a handful information on gas-liquid flow in wavy microchannels. This paper presents an experimental study of the air-water two-phase flow in a wavy microchannels heatsink typically used for cooling of electronic devices. The heatsink was made from copper with 26 wavy channels. The width and the depth of each channel were 800 and 500 μm, respectively. Two cartridge heaters were used to generate a constant heat flux. Pressure and temperatures of the working fluid at the inlet and outlet of the test section were measured using a pressure transducer and thermocouples. To evaluate the benefit of using the wavy microchannels, experiment was performed using water as the working fluid. The results showed that the wavy microchannels could elevate the Nusselt number by up to 112% compared to that obtained from the straight microchannels. Experiments at various air and water flowrates in the wavy microchannels suggested that the advantages of using the air-water flow would be obvious at the liquid Reynolds number above 379.