An official website of the United States government
When fluids flow through straight channels sustained turbulence occurs only at high Reynolds numbers [typically Re∼O(1000)]. It is difficult to mix multiple fluids flowing through a straight channel in the low Reynolds number laminar regime [Re<O(100)] because in the absence of turbulence, mixing between the component fluids occurs primarily via the slow molecular diffusion process. This Letter reports a simple way to significantly enhance the low Reynolds number (in our case Re≤10) passive microfluidic flow mixing in a straight microchannel by introducing asymmetric wetting boundary conditions on the floor of the channel. We show experimentally and numerically that by creating carefully chosen two-dimensional hydrophobic slip patterns on the floor of the channels, we can introduce stretching, folding, and/or recirculation in the flowing fluid volume, the essential elements to achieve mixing in the absence of turbulence. We also show that there are two distinctive pathways to produce homogeneous mixing in microchannels induced by the inhomogeneity of the boundary conditions. It can be achieved either by (1) introducing stretching, folding and twisting of fluid volumes, i.e., via a horse-shoe type transformation map, or (2) by creating chaotic advection, achieved through manipulation of the hydrophobic boundary patterns on the floor of the channels. We have also shown that by superposing stretching and folding with chaotic advection, mixing can be optimized in terms of significantly reducing mixing length, thereby opening up new design opportunities for simple yet efficient passive microfluidic reactors.
more »
« less
Flow Modeling in Long Surface Patterned Micromixers Using Division in Multiple Geometrical Subunits
Optimization of mixing in microfluidic devices is a popular application of computational fluid dynamics software packages, such as COMSOL Multiphysics, with an increasing number of studies being published on the topic. On one hand, the laminar nature of the flow and lack of turbulence in this type of devices can enable very accurate numerical modeling of the fluid motion and reactant/particle distribution, even in complex channel geometries. On the other hand, the same laminar nature of the flow, makes mixing, which is fundamental to the functionality of any microfluidic reactor or assay system, hard to achieve, as it forces reliance on the slow molecular diffusion, rather than on turbulence. This in turn forces designers of microfluidic systems to develop a broad set of strategies to enable mixing on the microscale, targeted to the specific applications of interest. In this context, numerical modeling can enable efficient exploration of a large set of parameters affecting mixing, such as geometrical characteristics and flow rates, to identify optimal designs. However, it has to be noted that even very performant mixing topologies, such as the use of groove-ridge surface features, require multiple mixing units. This in turn requires very high resolution meshing, in particular when looking for solutions for the convection-diffusion equation governing the reactant or chemical species distribution. For the typical length of microfluidic mixing channels, analyzed using finite element analysis, this becomes computationally challenging due to the large number of elements that need to be handled. In this work we describe a methodology using the COMSOL Computational Fluid Dynamics and Chemical Reaction Engineering modules, in which large geometries are split in subunits. The Navier-Stokes and convection-diffusion equations, are then solved in each subunit separately, with the solutions obtained being transferred between them to map the flow field and concentration through the entire geometry of the channel. As validation, the model is tested against data from mixers using periodic systems of groove-ridge features in order to engineer transversal mixing flows, showing a high degree of correlation with the experimental results. It is also shown that the methodology can be extended to long mixing channels that lack periodicity and in which each geometrical mixing subunit is distinct.
more »
« less
- Award ID(s):
- 1659541
- PAR ID:
- 10084593
- Date Published:
- Journal Name:
- Comsol-Newsletter
- ISSN:
- 1869-5388
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.more » « less
-
Microfluidic cell sorters have shown great potential to revolutionize the current technique of enriching rare cells. In the past decades, different microfluidic cell sorters have been developed by researchers for separating circulating tumor cells, T-cells, and other biological markers from blood samples. However, it typically takes months or even years to design these microfluidic cell sorters by hand. Thus, researchers tend to use computer simulation (usually finite element analysis) to verify their designs before fabrication and experimental testing. Despite this, conducting precision finite element analysis of microfluidic devices is computationally expensive and labor-intensive. To address this issue, we recently presented a microfluidic simulation method that can simulate the behavior of fluids and particles in some typical microfluidic chips instantaneously. Our method decomposes the chip into channels and intersections. The behavior of fluid in each channel is determined by leveraging analogies with electronic circuits, and the behavior of fluid and particles in each intersection is determined by querying a database containing 92,934 pre-simulated channel intersections. While this approach successfully predicts the behavior of complex microfluidic chips in a fraction of the time required by existing techniques, we nonetheless identified three major limitations with this method: (1) the library of pre-simulated channel intersections is unnecessarily large (only 2,072 of 92,934 were used); (2) the library contains only cross-shaped intersections (and no other intersection geometries); and (3) the range of fluid flow rates in the library is limited to 0 to 2 cm/s. To address these deficiencies, in this work we present an improved method for instantaneously simulating the trajectories of particles in microfluidic chips. Firstly, inspired by dynamic programming, our new method optimizes the generation of pre-simulated intersection units and avoids generating unnecessary simulations. Secondly, we constructed a cloud database (http://cloud.microfluidics.cc) to share our pre-simulated results and to let users become contributors and upload their simulation results into the cloud database as a benefit to the whole microfluidic simulation community. Lastly, we investigated the impact of different channel angles and different fluid flow rates on predicting the trajectories of particles. We found a wide range of device geometries and flow rates over which our existing simulation results can be extended without having to perform additional simulations. Our method should accelerate the simulation of particles in microfluidic chips and enable researchers to design new microfluidic cell sorter chips more efficiently.more » « less
-
Our recently published papers reporting results of Direct Numerical Simulation (DNS) of forced convection flows in porous media suggest that in a porous medium the size of turbulent structures is restricted by the pore scale. Since the turbulent kinetic energy is predominantly contained within large eddies, this suggests that turbulent flow in a porous medium may carry less energy that its counterpart in a clear fluid domain. We use this insight to develop a practical model of turbulent flow in composite porous/fluid domains. In such domains, most of the flow is expected to occur in the clear fluid region; therefore, in most cases the flow in the porous region either remains laminar or starts its transition to turbulence even if the flow in the clear fluid region is fully turbulent. This conclusion is confirmed by comparing appropriate Reynolds numbers with their critical values. Therefore, for most cases, using the Forchheimer term in the momentum equation and the thermal dispersion term in the energy equation may result in a sufficiently good model for the porous region. However, what may really affect turbulent convection in composite domains is the roughness of the porous/fluid interface. If particles or fibers that constitute the porous medium (and the pores) are relatively large, the impact of the roughness on convection heat transfer in composite porous/fluid domains may be much more significant than the impact of possible turbulence in the porous region. We use the above considerations to develop a practical model of turbulent flow in a composite porous/fluid domain, concentrating on the effect of interface roughness on turbulence.more » « less
-
The phenomenon of ferrofluid-water mixing is investigated using a double-layer magnetic micromixer, in which a layer of micromagnet bars is placed immediately below the fluid layer. A wavy pattern of the ferrofluid–water interface is surprisingly observed at each micromagnet responsible for improved mixing. The mechanism causing the wavy mixing is discovered and analyzed through experimental measurements and numerical simulations, and the mixing efficiency under different flow conditions is discussed. For flows with Re ≪ 1, the resultant steep gradient of opposing magnetic forces by micromagnets in the ferrofluid region gives rise to a local pressure source that induces a transverse/spanwise pressure gradient and activates momentum transfer between fluids. The current finding enables effective localized mixing of ferrofluids with a small footprint and, thus, has great potential to achieve fast mixing for high-throughput flows with an integrated parallel system of multiple microfluidic channels and micromagnets.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- The DOI is not currently available.
