Flow In Microchannel Research Articles

Magnetohydrodynamic unsteady natural convection slip flow in a vertical parallel plate microchannel heated with constant heat flux.

This research presents a numerical study over the unsteady natural convection of an electrically conducting fluid in an open-ended vertical parallel plate microchannel under uniform and asymmetric heat flux subjected to a uniform lateral magnetic field. Slip velocity, as well as temperature jump at channel walls, are modeled using a first-order model. The effects of Knudsen number)Kn(, heat flux ratio)rq(, Grashof number)Gr(, and Hartmann number)M(on mass flow rate, the maximum temperature of the wall, and average Nusselt (Nu) as a function of time are discussed. In this research, the flow in the limit of 0≥ M≥2; 0≥Kn≥0.1; 0.1≥ rq≥1, and 0.0001≥ L≥0.01 were simulated numerically. Furthermore, and study is conducted on steady-state velocity and temperature profiles. The research showed that increasing the Ha value, causes a reduction in the mass flow and an increase in the maximum temperature of the wall with respect to time. As M increases, the average Nusselt number as a function of time diminishes. Also, these variables present at least one overshoot or undershoot in time. The overshoot and undershoot decrease with increasing Ha. It is also revealed that the magnetic force effect decreases the fluid flow and thermal behavior with an increase in Kn. In addition, in higher Gr, magnetic force influence on profiles of velocity and temperature decreases for all rq and Kn values. Moreover, the magnetic field has no significant effect on the time at which steady state condition is attained for all rq and Gr number values.

Shapes of surfactant-laden Taylor bubbles in a square microchannel

Experiments on contaminated Taylor flows in a square microchannel were carried out to investigate the effects of surfactant on the bubble shape in the nose and tail regions for different surfactant properties. The nose curvature was found to be proportional to the bubble length at low surfactant concentrations, while it was independent of the concentration at high concentrations. The rate of increase in the nose curvature at the former concentrations can be expressed in terms of the surface coverage ratio. The bubble velocity decreased with increasing the nose curvature, whereas the surface tension reduced by surfactant adsorption worked better to correlate the velocity data. The curvature of the bubble tail increased steeply at low concentrations as a consequence of the early coverage due to interfacial advection. The tail curvature also had a strong correlation with the surface coverage ratio.

Analysis of the flow and heat transfer performance of nonlinear variable cross-section microchannels based on analytical method

Analysis of the flow and heat transfer performance of nonlinear variable cross-section microchannels based on analytical method

Dual effects of pressure difference for long-wave stratified flow in horizontal microchannels

In this work, the linear stability analysis combining with energy balance theory is performed to clarify the suppressing and boosting long-wave instabilities driven by pressure difference in horizontal microchannels as well as its corresponding physical meanings. It is found that, for a specific fluid combination, there exists a certain relative liquid film thickness (hN,c) at which the disturbance work done by interfacial shear stress over unit wavelength (denoted by ISS) is exactly equal to the total viscous dissipation rate over unit wavelength within both phase bulk flows (denoted by DIStot). When the liquid film thickness (hN) exceeds hN,c, ISS is smaller than DIStot, indicating that the initial disturbance energy cannot be replenished from the primary flow until its completely disappear due to the viscosity dissipation. Conversely, ISS is larger than DIStot, indicating that the initial disturbance can acquire energy from the primary flow and grow up. The strength of suppressing and boosting effect can be influenced by liquid film thickness, flow rate, and channel height. However, the demarcation of the dual effects remains unchanged. An analytical expression of hN,c with respect to the viscosity ratio is established by means of long-wave approximation method, which matches well with the numerical results. The energy balance mechanisms revealed in the present study provide a new insight into the wave mode in the confined space and help the design of microfluidic systems.

Cavitation characteristics in rectangular flow restriction microchannel seals

The shaft seal of the primary pump relies on microchannel flow to prevent coolant leakage, but cavitation increases resistance, reduces efficiency, and damages equipment, affecting safety. Therefore, studying microchannel cavitation is crucial for optimizing seal performance and equipment reliability. This study explores the cavitation flow characteristics within a 0.1 mm rectangular flow restriction microchannel seal through experimental methods and numerical simulations. First, using a microfluidic experimental platform and visualization techniques, we directly observed the cavitation phenomena in the microchannel and quantified the evolution of cavitation bubbles under different conditions. Second, by establishing a numerical model of the microchannel, we analyzed the formation mechanisms of cavitation and its impact on flow characteristics under various conditions. The results show that lower cavitation numbers significantly increase the length of cavitation jets and the distribution range of cavitation bubbles, intensifying the cavitation phenomena within the microchannel. Additionally, microchannels with closely spaced double flow restriction structures are more effective at inhibiting cavitation than channels with a single flow restriction structure. These findings offer valuable insights into the mechanisms of microchannel cavitation and provide scientific guidance for the design and optimization of shaft seals.

Linear temporal instabilities and transient energy growth in rotating curved microchannel flow

Linear temporal instabilities and transient energy growth in rotating curved microchannel flow

A slip length analysis of microchannel flow with smooth and structured surfaces

One of the most significant challenges in microchannel design is the reduction of pressure drops within miniature fluidic pathways. A promising approach involves the manipulation of boundary conditions, mainly through the integration of structured hydrophobic surfaces. The mathematical characterization of these hydrophobic surfaces is achieved through the application of the Navier boundary condition, with the slip length identified as the critical parameter of interest. This paper delves into an in-depth exploration of how varying slip lengths impact the liquid flow dynamics within microchannels with a rectangular cross section. We consider a diverse range of microchannel hydraulic diameters 5–200 μm and aspect ratios from 1:1 to 1:20 for Reynolds numbers in the laminar range 10–1000. Three-dimensional calculations are performed on both conventional smooth microchannels and those equipped with air-filled groove-structured surfaces. The results are then compared with the analytical solution and experimental results. The application of a hydrophobic structure to a single microchannel surface is resulted in the friction factor reduction by over 30%, while the application of the hydrophobic structure to two opposing walls led to a reduction by over 60%. The maximum throughput is shown to be achieved for a bubble protrusion angle of approximately 0° for a microchannel with a single hydrophobic wall. The greatest reduction in the friction factor was achieved when the bubbles were positioned in a staggered configuration at bubble protrusion angles of approximately –25° for microchannel with two opposing hydrofobic walls.

Enhanced flow boiling by manipulating two-phase flow in Tesla channel heat sink using HFE-7100

Flow boiling of dielectric fluids in copper microchannel heat sinks is highly desirable for cooling large-sized insulated gate bipolar transistor (IGBT) power electronic modules. However, dielectric fluids present challenges in flow boiling because of their unfavorable thermophysical properties. These factors make it difficult to enhance critical heat flux (CHF) without precooling. To address this, we developed a large copper heat sink (10 cm × 5 cm) with Tesla microchannels designed to suppress vapor backflow and promote intense fluid mixing. The microchannels have a high length to hydrodynamic diameter ratio of approximately 220, significantly higher than those in previous studies. Flow boiling experiments using HFE-7100 were conducted for both Tesla and plain-wall microchannels. Tesla channels demonstrated a 26.2 % increase in CHF and a 120 % improvement in heat transfer coefficient (HTC). These enhancements are attributed to the vapor backflow suppression and improved fluid mixing. Moreover, the standard deviation of wall temperature in plain-wall microchannels was 10 times higher than in Tesla channels, highlighting the effectiveness of the periodic Tesla valves in reducing two-phase flow instabilities. Flow pattern visualization was conducted to further understand the mechanism behind vapor regulation, clarifying the role of Tesla valves in controlling vapor backflow. This study demonstrates the potential of dielectric fluids in Tesla microchannels for flow boiling applications, offering a promising solution for cooling large electronics.

Femtosecond laser ablation of nitrocellulose for spatio-temporal flow control in µPADs

Microfluidic paper-based analytical devices (μPADs) are gaining popularity due to their low cost and ease of use, but controlling fluid flow for more complex biochemical assays within these devices remains challenging. This study investigates femtosecond laser ablation of nitrocellulose (NC), a preferred material for µPADs, to create mechanically switchable barriers and flow controllers. We investigated NC ablation using single laser pulses and spatially overlapping pulses that generate lines. Single pulse ablation thresholds were determined for wavelengths of 1,030 nm, 515 nm and 343 nm. Line ablation characteristics were investigated as a function of the temporal and spatial pulse separation and laser wavelength. High aspect ratio grooves (up to 4.26) were achieved under specific conditions. These grooves can be used to define the spatial separation of the flow in separated microchannels or to form a barrier line perpendicular to the microchannel that can modulate the temporal behavior of the fluid flow. This barrier introduced an additional high flow resistance slowing down the flow or, if it was designed to cut through nitrocellulose at the entire depth, completely stopped the liquid flow. It was further shown that a barrier formed in this way could be switched by mechanically bending the µPAD at the barrier position. The femtosecond laser patterning method presented here provides precise spatio-temporal control not only for flow branching and multiplexing, but also for controlling flow speed and switching flow on and off within the same manufacturing process. Our results open up new possibilities for complex, multi-step assays on µPADs.

Characterization of the Dynamic Flow Response in Microfluidic Devices.

The purpose of this study is to characterize the dynamic response of fluid flow in microchannels, which can show significant delay times before reaching steady flow conditions. Two main sources of these delays are numerically and experimentally investigated, the hydraulic compliance which originates from the flexibility of the system components (microchannel, tubing, syringe, etc.), and the compressibility of the liquid dead volume in the setup, also known as the "bottleneck effect". A fluid-structure interaction model is presented for the compliance of rectangular PDMS microchannels that is used to form a numerically based relation for the compliance as a function of the pressure and geometry. This relation is successfully able to predict the dynamics of the flow inside PDMS microchannels in stop-flow experiments. The time delays associated with the bottleneck effect is also shown when using different syringe volumes, microchannel resistances, and liquid types. In these tests, the bottleneck effect has a much larger effect compared to the compliance of the PDMS microchannels. This is true even when using softer PDMS by increasing the monomer-to-curing agent mixing ratio. The characterization that is presented here allows for a simple analysis of microfluidic networks using the hydraulic-circuit approach.

Interaction between vapor bubbles during flow boiling heat transfer in microchannels

Microchannel flow boiling is an efficient cooling solution for high-power-density miniaturized systems. Many studies on microchannel flow boiling focused on the dynamics of single vapor bubbles, while neglecting the interaction between bubbles, which is important in relevant applications. Here, numerical simulations are carried out to study the interaction between multiple vapor bubbles in microchannel flow boiling. The results show that for different numbers of bubbles in the microchannels with the same initial size and position of leading bubbles, the bubble size in a single-bubble microchannel is larger compared to the leading bubble of multiple-bubble cases because of heat absorption by the vaporization at the rear bubbles. As the initial volume ratio between the leading bubble and the rear bubble decreases, the leading bubble size in the downstream becomes smaller because of the reduced contact with the superheated thermal boundary layer. With increasing the Reynolds number, both the leading and the trailing bubbles increase slightly in size in the upstream of the heated region, because the bubbles at higher Reynolds number move faster and firstly get in contact with the superheated fluid. The increase in the bottom wall thickness increases the growth rate of the multiple bubble sizes with earlier bubble coalescence because of the higher upstream wall temperature by heat conduction in the solid wall.

Flow boiling of HFE-7100 for cooling Multi-Chip modules using manifold microchannels

Two-phase manifold microchannel heat sinks have gained significant interest for their efficient use of the coolant’s latent heat. Despite this, manifold microchannel flow boiling heat transfer has not yet been applied to wide-bandgap semiconductor power modules with multiple chips, primarily due to the complexity of the process and challenges related to electrical insulation. This study introduces a novel embedded manifold microchannel design that ensures uniform mass flow distribution, tailored for the thermal management of multi-chip power modules. The microchannels were laser-etched onto direct bonded copper (DBC) to reduce thermal resistance while maintaining electrical insulation. Thermal test vehicles (TTVs) for multi-chip power modules, featuring two different microchannel widths, were assembled using silver sintering. Experimental tests were then conducted using HFE-7100 as the coolant to evaluate two-phase heat dissipation performance. Results indicate that during flow boiling heat transfer, the temperature difference between chips in a power module strongly correlates with the exit quality. A higher coolant mass flow rate significantly reduces temperature variation between chips, particularly under high chip heat flux. At a coolant mass flow rate of 9 g/s, with a chip heat flux of 357 W/cm2 and total heat dissipation of 536 W, the minimum thermal resistance reached 0.15 cm2∙K/W, yielding a COP of 1391. With a slight sacrifice in thermal resistance, 0.17 cm2∙K/W and 0.20 cm2∙K/W were achieved at mass flow rates of 6 g/s and 3 g/s, respectively. Correspondingly, the COPs reached 2179 and 6749. This research offers valuable insights for applying flow boiling heat transfer with dielectric coolants to cool multi-chip power modules.

Thermal and hydrodynamic characteristics of microencapsulated phase change materials slurry flow in wavy microchannels

Microencapsulated phase change material slurry (MPCS) is a novel cooling fluid with promising performance. This study numerically investigates heat transfer and flow characteristics of MPCS within wavy microchannels. MPCS is modeled as a homogeneous, Newtonian fluid. Five wavy geometries were examined across a Reynolds number range of 50 to 250 (laminar flow regime), varying in amplitude and wavelength. The model results show that with increase in Reynolds number and decrease in the amplitude and radius of curvature of wavy microchannels, the pressure drop and Nusselt number also increase. Furthermore, the results reveal that Dean vortices intensify with increasing wave amplitude and Reynolds number. Conversely, these vortices weaken as channel wavelength and radius of curvature increase. The formation of Dean vortices enhances fluid mixing and consequently improves the thermal performance of the slurry. The study concludes that in wavy microchannels employing MPCS, increasing the Reynolds number, increasing the channel amplitude with respect to wavelength, and decreasing the radius of curvature improves the overall performance. Also, the results reveal that in higher Reynolds numbers, the radius of curvature is the most effective parameter on the overall performance of wavy microchannels.

A microfluidic approach for controllable synthesis of TiO2 submicrospheres

TiO2 submicrospheres are used in dye-sensitized solar cells and Li-ion batteries to enhance energy conversion efficiency and volumetric performance. In this study, a microfluidic chip was designed for the realization of the continuous controllable synthesis of TiO2 submicrospheres. The chip utilizes two T-channels to discrete dispersed phases into droplets, and a Y-channel to fuse the two dispersed phase droplets. The droplets generation process, the influencing factors of droplet length and the fusion process of droplets were investigated. The effects of the flow rate ratio of titanium oxysulfate and ammonia solutions, microchannel geometrical profiles including the Y-channel entrance angle on the size of the synthesized TiO2 were studied. The length of the generated droplets was found to decrease with increasing continuous-phase flow rate, and increase with increasing dispersed-phase flow rate and microchannel cross-section area. When the flow rate ratio of titanium oxysulfate solution and ammonia solution was 1:1, the average diameter of TiO2 spheres synthesized was the smallest, around 300 nm. Increasing the microchannel cross-section and Y-channel entrance angle were beneficial in obtaining larger TiO2 submicrospheres. This proposed microfluidic strategy has the potential for applications required continuous synthesis of TiO2 submicrospheres with controllable sizes.

An analytical gradient method to model interfacial heat and mass transfer in multi-field CFD codes

Phase change at the interface between a liquid and a gas, commonly referred to as bulk condensation/boiling or interface condensation/boiling, is often overshadowed by wall condensation/boiling. However, interfacial phase change plays a critical role in various industrial applications, including boiling flows in microchannels and boiling liquid metal flows. In this article, we introduce a novel analytical model based on the gradient method. This model accounts for interfacial phase change within multi-field codes. It is validated against both analytical and experimental cases involving bulk boiling, demonstrating excellent agreement. Notably, the mesh convergence aligns with the methods employed in single-fluid codes. Additionally, we propose a hybrid approach that combines the high accuracy of this model with the computational efficiency of dispersed methods.

Development of a sensor for liquid film thickness measurements during annular flow in microchannels

Development of a sensor for liquid film thickness measurements during annular flow in microchannels

Numerical Simulations of Flow and Heat Exchange in Zigzag-Shaped Microchannels

Cooling of electronic components is of great importance currently. One of the most important methods of integrated circuit cooling is the use of microchannel heat sinks. The aim of this work is to analyze the cooling capacity of zigzag shaped microchannels. The heat exchange for four different microchannels was tested depending on the flow rate of the cooling liquid (water) and the temperature of the cooled element. The system of differential equations that describe fluid flow and heat transport is presented in the paper. The equations were solved using the finite volume method. The work showed that both the increase in the number of zigzag sections and the increase in the flow velocity cause an increase in the flow resistance. In contrast, an increase in the temperature of the cooled element causes a decrease in fluid viscosity, which results in a decrease in flow resistance. From the point of heat exchange, both the increase in the number of segments, the flow rate, and the temperature of the cooled element increase the cooling capacity of the microchannel. Research shows that zigzag shaped microchannels are excellent cooling elements, especially for geometrically small heat sources such as electronic components.

Unsteady numerical calculation of flow boiling heat transfer in microchannels with Kagome structures

A three-dimensional microchannel flow boiling model was developed based on the VOF model and the Lee model, and the flow and heat transfer in microchannels with Kagome structures were numerically simulated. The characteristics of bubble flow patterns, gas phase distribution, and heat transfer were analyzed under the heat fluxes of 5–8 MW/m2 and flow velocities of 0.2–0.5 m/s. The results show that at a flow velocity of 0.5 m/s, the disturbances caused by fluid impact on the Kagome structure result in bubble break-up, effectively suppressing the formation of large bubbles. However, at a flow velocity of 0.2 m/s, bubbles tend to be adhered on the surface of the Kagome structure, causing local blockage and deteriorating heat transfer. Compared to the flow velocity of 0.5 m/s, the average heat transfer coefficient in the microchannel decreases by up to 75.4 % and the pressure loss decreases by 45.38 % at the flow velocity of 0.2 m/s.

Predictive Model for Cell Positioning during Periodic Lateral Migration in Spiral Microchannels.

The periodic lateral migration of submicrometer cells is the primary factor leading to low precision in a spiral microchannel during cell isolation. In this study, a mathematical predictive model (PM) is derived for the lateral position of cells during the periodic lateral migration process. We analyze the relationship of migration period, migration width, and starting point of lateral migration with microchannel structure and flow conditions and determine the empirical coefficients in PM. Results indicate that the aspect ratio of the microchannel and the Reynolds number (Re) are key factors that influence the periodicity of the cell lateral migration. The lateral migration width is jointly affected by Re, the cell blockage ratio, and the microchannel curvature radius. The inlet structure of the microchannel and the ratio of the cell sample to the sheath flow rate are critical parameters for regulating the initial position. Moreover, the structure of the pressure field at the inlet constrains the distribution range of the starting point of the lateral migration. Regardless of whether the particles/cells undergo 0.5, 1, or multiple lateral migration cycles, the lateral positions predicted by PM align well with the experimental observations, thus verifying the accuracy of PM. This research helps to elucidate the characteristics of periodic lateral migration of cells in spiral microchannels and can provide practical guidance for the development and optimization of miniature spiral microchannel chips for precise cell isolation.

Heat and mass transfer analysis of s-PTT nanofluid in microchannels under combined electroosmotic and pressure-driven flows with wall slip using the homotopy perturbation method

The heat and mass transfer of the electroosmotic flow in microchannel transporting viscoelastic nanofluid is investigated considering Brownian motion of nanoparticles and slip boundary conditions. The simplified Phan-Thien-Tanner model is employed to describe the rheological behavior of fluid and the nonlinear Navier model with non-zero slip critical shear stress is considered at walls. The governing nonlinear momentum, mass, and heat transfer equations are solved using the Homotopy Perturbation Method. The study reveals that increasing the fluid elasticity, nanoparticle concentration, and size significantly enhances the flow rate, heat and mass transfer. Additionally, elasticity and Reynolds number decrease the friction factor. Reducing the double-layer thickness and increasing the Reynolds number lead to higher flow rates and fluid velocities. Notably, the findings emphasize the critical role of the slip conditions on the Sherwood and Nusselt numbers.