Research Papers: Flows in Complex Systems

An Enhanced One-Layer Passive Microfluidic Mixer With an Optimized Lateral Structure With the Dean Effect

[+] Author and Article Information
Teng Zhou

School of Mechanical Engineering,
Yeungnam University,
Gyeongsan 712-749, South Korea

Yifan Xu, Zhenyu Liu

Changchun Institute of Optics
Fine Mechanics and Physics (CIOMP),
Chinese Academy of Science,
Changchun, Jilin 130033, China

Sang Woo Joo

School of Mechanical Engineering,
Yeungnam University,
Gyeongsan 712-749, South Korea
e-mail: swjoo@yu.ac.kr

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received January 14, 2015; final manuscript received March 29, 2015; published online May 19, 2015. Assoc. Editor: Shizhi Qian.

J. Fluids Eng 137(9), 091102 (May 19, 2015) (7 pages) Paper No: FE-15-1026; doi: 10.1115/1.4030288 History: Received January 14, 2015

Topology optimization method is applied to a contraction–expansion structure, based on which a simplified lateral flow structure is generated using the Boolean operation. A new one-layer mixer is then designed by sequentially connecting this lateral structure and bent channels. The mixing efficiency is further optimized via iterations on key geometric parameters associated with the one-layer mixer designed. Numerical results indicate that the optimized mixer has better mixing efficiency than the conventional contraction–expansion mixer for a wide range of the Reynolds number.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Lee, C.-Y. , Chang, C.-L. , Wang, Y.-N. , and Fu, L.-M. , 2011, “Microfluidic Mixing: A Review,” Int. J. Mol. Sci., 12(5), pp. 3263–3287. [CrossRef] [PubMed]
Cartier, C. A. , Drews, A. M. , and Bishop, K. J. M. , 2014, “Microfluidic Mixing of Nonpolar Liquids by Contact Charge Electrophoresis,” Lab Chip, 14(21), pp. 4230–4236. [CrossRef] [PubMed]
Afzal, A. , and Kim, K.-Y. , 2012, “Passive Split and Recombination Micromixer With Convergentcdivergent Walls,” Chem. Eng. J., 203, pp. 182–192. [CrossRef]
Ammar, H. , Ould el Moctar, A. , Garnier, B. , and Peerhossaini, H. , 2014, “Flow Pulsation and Geometry Effects on Mixing of Two Miscible Fluids in Microchannels,” ASME J. Fluids Eng., 136(12), p. 121101. [CrossRef]
Karami, M. , Shirani, E. , Jarrahi, M. , and Peerhossaini, H. , 2014, “Mixing by Time-Dependent Orbits in Spatiotemporal Chaotic Advection,” ASME J. Fluids Eng., 137(1), p. 011201. [CrossRef]
Maki, A.-J. , Hemmila, S. , Hirvonen, J. , Girish, N. N. , Kreutzer, J. , Hyttinen, J. , and Kallio, P. , 2014, “Modeling and Experimental Characterization of Pressure Drop in Gravity-Driven Microfluidic Systems,” ASME J. Fluids Eng., 137(2), p. 021105. [CrossRef]
Mohammadi, M. , and Sharp, K. V. , 2013, “Experimental Techniques for Bubble Dynamics Analysis in Microchannels: A Review,” ASME J. Fluids Eng., 135(2), p. 021202. [CrossRef]
Sen, A. K. , and Bhardwaj, P. , 2012, “Microfluidic System for Rapid Enumeration and Detection of Microparticles,” ASME J. Fluids Eng., 134(11), p. 111401. [CrossRef]
Solovitz, S. A. , Zhao, J. , Xue, W. , and Xu, J. , 2013, “Uniform Flow Control for a Multipassage Microfluidic Sensor,” ASME J. Fluids Eng., 135(2), p. 021101. [CrossRef]
Qian, S. , and Bau, H. H. , 2002, “A Chaotic Electroosmotic Stirrer,” Anal. Chem., 74(15), pp. 3616–3625. [CrossRef] [PubMed]
Qian, S. , and Bau, H. H. , 2005, “Magneto-Hydrodynamic Stirrer for Stationary and Moving Fluids,” Sens. Actuators B, 106(2), pp. 859–870. [CrossRef]
Liang, L. , and Xuan, X. , 2012, “Diamagnetic Particle Focusing Using Ferromicrofluidics With a Single Magnet,” Microfluid. Nanofluid., 13(4), pp. 637–643. [CrossRef]
Wen, C.-Y. , Liang, K.-P. , Chen, H. , and Fu, L.-M. , 2011, “Numerical Analysis of a Rapid Magnetic Microfluidic Mixer,” Electrophoresis, 32(22), pp. 3268–3276. [CrossRef] [PubMed]
Cardoso, V. F. , Knoll, T. , Velten, T. , Rebouta, L. , Mendes, P. M. , Lanceros-Mendez, S. , and Minas, G. , 2014, “Polymer-Based Acoustic Streaming for Improving Mixing and Reaction Times in Microfluidic Applications,” RSC Adv., 4(9), pp. 4292–4300. [CrossRef]
SadAbadi, H. , Packirisamy, M. , and Wuthrich, R. , 2013, “High Performance Cascaded pdms Micromixer Based on Split-and-Recombination Flows for Lab-On-A-Chip Applications,” RSC Adv., 3(20), pp. 7296–7305. [CrossRef]
Amini, H. , Lee, W. , and Di Carlo, D. , 2014, “Inertial Microfluidic Physics,” Lab Chip, 14(15), pp. 2739–2761. [CrossRef] [PubMed]
Zhang, J. , Li, W. , Li, M. , Alici, G. , and Nguyen, N.-T. , 2014, “Particle Inertial Focusing and Its Mechanism in a Serpentine Microchannel,” Microfluid. Nanofluid., 17(2), pp. 305–316. [CrossRef]
Amini, H. , Sollier, E. , Masaeli, M. , Xie, Y. , Ganapathysubramanian, B. , Stone, H. A. , and Di Carlo, D. , 2013, “Engineering Fluid Flow Using Sequenced Microstructures,” Nat. Commun., 4, p. 1826. [CrossRef] [PubMed]
Lee, M. G. , Choi, S. , and Park, J.-K. , 2009, “Rapid Laminating Mixer Using a Contraction–Expansion Array Microchannel,” Appl. Phys. Lett., 95(5), p. 051902. [CrossRef]
Lee, M. G. , Choi, S. , and Park, J.-K. , 2010, “Rapid Multivortex Mixing in an Alternately Formed Contraction-Expansion Array Microchannel,” Biomed. Microdevices, 12(6), pp. 1019–1026. [CrossRef] [PubMed]
Lee, M. G. , Choi, S. , and Park, J. K. , 2011, “Inertial Separation in a Contraction–Expansion Array Microchannel,” J. Chromatogr. A, 1218(27), pp. 4138–4143. [CrossRef] [PubMed]
Lee, M. G. , Shin, J. H. , Bae, C. Y. , Choi, S. , and Park, J.-K. , 2013, “Label-Free Cancer Cell Separation From Human Whole Blood Using Inertial Microfluidics at Low Shear Stress,” Anal. Chem., 85(13), pp. 6213–6218. [CrossRef] [PubMed]
Liu, Z. , Deng, Y. , Lin, S. , and Xuan, M. , 2012, “Optimization of Micro Venturi Diode in Steady Flow at Low Reynolds Number,” Eng. Optim., 44(11), pp. 1389–1404. [CrossRef]
Deng, Y. , Liu, Z. , Zhang, P. , Liu, Y. , Gao, Q. , and Wu, Y. , 2012, “A Flexible Layout Design Method for Passive Micromixers,” Biomed. Microdevices, 14(5), pp. 929–945. [CrossRef] [PubMed]
Deng, Y. , Zhang, P. , Liu, Y. , Wu, Y. , and Liu, Z. , 2013, “Optimization of Unsteady Incompressible Navier–Stokes Flows Using Variational Level Set Method,” Int. J. Numer. Methods Fluids, 71(12), pp. 1475–1493. [CrossRef]
Borrvall, T. , and Petersson, J. , 2003, “Topology Optimization of Fluids in Stokes Flow,” Int. J. Numer. Methods Fluids, 41(1), pp. 77–107. [CrossRef]
Gersborg-Hansen, A. , Sigmund, O. , and Haber, R. B. , 2005, “Topology Optimization of Channel Flow Problems,” Struct. Multidiscip. Optim., 30(3), pp. 181–192. [CrossRef]
Olesen, L. H. , Okkels, F. , and Bruus, H. , 2006, “A High-Level Programming-Language Implementation of Topology Optimization Applied to Steady-State Navier–Stokes Flow,” Int. J. Numer. Methods Eng., 65(7), pp. 975–1001. [CrossRef]
Donea, J. , and Huerta, A. , 2003, Finite Element Methods for Flow Problems, Wiley, Hoboken, NJ.
Deng, Y. , Liu, Z. , Zhang, P. , Liu, Y. , and Wu, Y. , 2011, “Topology Optimization of Unsteady Incompressible Navier–Stokes Flows,” J. Comput. Phys., 230(17), pp. 6688–6708. [CrossRef]
Liu, Z. , Gao, Q. , Zhang, P. , Xuan, M. , and Wu, Y. , 2011, “Topology Optimization of Fluid Channels With Flow Rate Equality Constraints,” Struct. Multidiscip. Optim., 44(1), pp. 31–37. [CrossRef]


Grahic Jump Location
Fig. 1

Schematic of CAE in (a) front view and (b) Isometric drawing. Blue shade shows the optimization domain. Ha: width of expansion domain; La: length of expansion domain; Lp: distance between the periods of the mixer; Hb: width of contraction domain; and D: channel depth.

Grahic Jump Location
Fig. 2

The structure of the design domain for the simulation (a) before and (b) after optimization. Initially rigid domain, shown in solid black, is partially converted to fluid channel, shown in white, through optimization process.

Grahic Jump Location
Fig. 3

Lateral structure. R and r are the radii of the semicircles centered at O1 and O2, respectively. Lo1 is the position of center of semicircles in x direction, and Lo1 = R = 1/2 La. Lb and Lc are widths of two rectangles, and Lb + Lc = 2R. Ho1 depicts the position of O2 and the height of the small rectangle. Hb is the width for the contraction channel.

Grahic Jump Location
Fig. 4

(a) The streamline. The color denotes the velocity. The cross plane A-A is 20 μm from the inlet to lateral structure, B-B cross plane is about 20 μm to the outlet of lateral structure. The velocity field and the (v, w) arrow plot of the cross planes (b) A-A and (c) B-B.

Grahic Jump Location
Fig. 5

Streamlines in (a) y direction with distance from the wall of 10 μm, 1/5 of channel width (b) z direction with distance from the bottom of 50 μm, 4/5 of channel depth, the color denotes the distance from the wall in z direction above

Grahic Jump Location
Fig. 6

Arrangement of lateral structures (a) in two sides, G denotes the distance between the periodic structures. (b) The serpentine structures, Hb, denote the width of the turning, same as the width of main channel. Here, Ga = Gb, Ga + Gb + Hb = G.

Grahic Jump Location
Fig. 7

(a) The streamlines in the turning. The cross planes A-A and B-B to the turning is about 20 μm. The velocity field and (v, w) arrow plot for the cross planes A-A and B-B with (b) Re = 5 and (c) 150.

Grahic Jump Location
Fig. 8

(a) The arrangement of the lateral structure for the mixer. The benchmark is the original contraction–expansion mixer. T12-1 and T12-2 denote the chips with structures arranged in two sides for the period with one and two lateral structures. One side corresponds to the chip with all the lateral structures in the same side. S12-1 and S12-2 denote serpentine structure with a period of one and two lateral structures. (b) Mixing efficiency of chips with different arrangements versus Reynolds number, Re, compared against the contraction–expansion chip before optimization. (c) Mixing efficiency of chips with different depth versus Reynolds number, Re.

Grahic Jump Location
Fig. 9

Mixing efficiency for different number of periods for S12-1 mixer. (a), (b), (c), and (d) represent the cross sections of the first, second, third, and fourth contraction regions after the lateral structure, respectively, for Re = 25.

Grahic Jump Location
Fig. 10

Mixing efficiency for different number of periods for the original contraction–expansion mixer. (a)–(d) represent the cross sections of the first, second, third, and fourth contraction regions after the expansion structure, respectively, for Re = 25.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In