Numerical simulation of particle focusing in a symmetrical serpentine microchannel

Di Jiang a, Wenlai Tanga, Nan Xiang*ab and Zhonghua Ni*a
aSchool of Mechanical Engineering, Jiangsu Key Laboratory for Design and Manufacture of Micro-Nano Biomedical Instruments, Southeast University, Nanjing, 211189, China. E-mail: nan.xiang@seu.edu.cn; nzh2003@seu.edu.cn
bState Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou, 310027, China

Received 1st April 2016 , Accepted 3rd June 2016

First published on 6th June 2016


Abstract

This paper elucidates the particle focusing mechanisms in a symmetrical serpentine microchannel. The particle trajectories and the fluid fields in serpentine microchannels are explored via the numerical simulation based on a lattice Boltzmann method (LBM)-immersed boundary method (IBM) model. Experiments are also performed to verify the obtained simulation results. The investigation results help us to clearly understand the important role of the Dean flow for particle focusing in symmetrical serpentine microchannels. The balance of the inertial lift force and the Dean flow drag force at lower flow intensities leads the particles to flow near the sidewalls. As the flow intensity becomes stronger, the alternation of the Dean flow direction has special hydrodynamic effects to focus or separate particles of different sizes. Because of the volume effects, large particles are more prone to rotate with the Dean vortex and can be separated from the small particles at specific flow intensities. We envision that the simulation results of the focusing behaviors of particles would be of great significance for efficiently designing and optimizing inertial microfluidics.


Introduction

The advent of inertial microfluidics1 has freed particle/cell manipulation from sheath flow controls2 and external force generators,3 which enables device miniaturization and operation simplification to be possible. Up to now, inertial microfluidics has been widely employed for particle/cell sheathless focusing,4 controllable trapping,5 size-dependent separation6 and other areas. The manipulation objects have also been diversified, which range from synthesized particles with various shapes7 to rare cells (e.g., circulating tumor cells (CTCs))8 in human blood. In addition to the classic straight channels, various channel geometrics (e.g., arc, spiral, contraction–expansion array and serpentine) have been successfully designed to induce an additional inertial effect (i.e., secondary Dean flow) for more flexible particle manipulation. Among these channel geometrics, the serpentine microchannels have attracted increasing research interests for precise particle/cell focusing because they have reduced equilibrium positions and can be parallel integrated more conveniently so that the processing throughput can be multiplied. Compared with the straight or spiral microchannels, flow patterns in serpentine microchannels appear in a more complex manner. The most distinctive difference between the serpentine channel flow and the spiral channel flow is that the direction of Dean flow in the former keeps alternating. This direction alternating leads the particle to swing in the channel cross-section, which promises to be a powerful method for hydrodynamic manipulation of particles/cells.

Di Carlo9 firstly proposed the serpentine microchannel for particle focusing, and found that in an asymmetrical serpentine channel, the four equilibrium positions observed in a common straight channel with a square cross-section could be reduced to one certain equilibrium position. Gossett et al.10 further investigated the variation of particle focusing positions from the top view under various flow intensities. However, the separation effect of particles by size in the asymmetrical serpentine microchannel is not obvious.11 Instead, the symmetrical serpentine channel may be a good candidate for the separation of differently sized particles. It has been confirmed in Zhang's experiments12 that the square-wave shaped serpentine channel could be employed to separate the particles of different sizes according to their different lateral focusing positions. Although a great success has been achieved in the application of serpentine channels, the particle focusing or separation mechanisms in these channel geometries remain unclear. Specifically, the above-mentioned experimental studies are incapable of showing the particle focusing behaviors in the cross-section and are difficult to study the interaction between the particle and the flows. These limitations can be easily overcome by using the numerical simulation.

The simulation research in this paper is implemented by using the lattice Boltzmann method coupled with immersed boundary method (LBM-IBM). The mesoscopic kinetic LBM has been paid increasing attentions nowadays for its complicated boundary adaptation, parallelizability13 and low memory consuming.14 Ladd15,16 firstly employed the LBM for the simulation of particle migration behaviors in straight channel flows and calculated the interaction forces through the momentum transfer between flow and solid nodes. Since the uniform mesh is employed to both flow field and solid particles in this method, the particle motion stability or simulation effectiveness is sacrificed. To overcome these shortcomings, Feng and Michaelides17 coupled the IBM18 with the LBM. In this method, the particle boundary is represented by the Lagrange nodes and the deformation between instantaneous and reference boundary define the force imposed on the particle, which guarantees the stability of the particle motion simulation.

This paper intends to study the particle focusing behaviors in a symmetrical serpentine microchannel via the numerical simulation and the experiments. In our numerical simulation, the particle structure is expressed by the finite element method (FEM) and the flow field is calculated by the lattice Boltzmann method (LBM). Then, the fluid part and the particle part are coupled using the immersed boundary method (IBM). Through applying the developed IB-LBM codes, the inertial focusing behaviors of 10 μm and 5 μm particles in a serpentine microchannel are analyzed under various flow intensities, and the interactions between particle and surrounding fluid field are also explored. Then, the experiments are performed to verify the findings obtained from the simulation.

IB-LBM model

The flow field in the microchannel is calculated by the standard D3Q19 single relaxation time LBGK model19 with the external force term,20 as follows:
 
image file: c6ra08374a-t1.tif(1)
Here, i represents the 19 discrete directions in the D3Q19 model, fi(x, t) is the velocity distribution function which defines the possibility of a fluid particle at position x and time t moving in the ith direction. ci is the lattice velocity and its 19 velocity directions are
 
image file: c6ra08374a-t2.tif(2)
where c = Δxtx is the space step and Δt is the time step). τ is the dimensionless relaxation time. feqi is the equilibrium distribution function, which is discretized from the Maxwell–Boltzmann equilibrium and has a form of
 
image file: c6ra08374a-t3.tif(3)
where u is the flow velocity and the weight coefficients ω0 = 1/3, ω1–6 = 1/18, ω7–18 = 1/36.

Fi(x, t) is the external force term and can be defined as

 
image file: c6ra08374a-t4.tif(4)
where F is the force acting on the fluid. The fluid density ρ and momentum ρu are respectively defined as
 
image file: c6ra08374a-t5.tif(5)

The kinematic viscosity can be calculated via the relaxation time τ as follows:

 
image file: c6ra08374a-t6.tif(6)

In order to learn the effects of forces on particles, a 3D finite element membrane model with the energy method21 is applied. Specifically, the particle surface is divided into triangular faces and the vertices of each triangular faces are the boundary nodes. External work W is calculated as W = ES + EB + EA + EV22,23 where the ES is the strain energy, EB is the bending energy, EA is the area energy and EV is the volume energy. The strain energy EST of a single deformed triangular face can be calculated as24

 
image file: c6ra08374a-t7.tif(7)
where I1 and I2 are respectively the first and second strain invariant and can be estimated as I1 = λ12 + λ22 − 2, I2 = (λ1 + λ2)2 − 1 (λ1, λ2 are principle stains). κS is the surface elastic shear modulus and κα is the area dilation modulus. The total strain energy ES is the sum of strain energies EST of each triangular face as
 
image file: c6ra08374a-t8.tif(8)

The bending energy EB describes the curvature contribution and can be written as23

 
image file: c6ra08374a-t9.tif(9)
Here κB is binding modulus. θ is the angle between the normals of a triangular face and its neighbor and θ0 represents the angle before the deformation of particles. The EB is the sum of angle change energy at each side of particle surfaces.

Similarly, EA and EV describe the surface and volume energy stored in the deformed particles and can respectively be described as23,25

 
image file: c6ra08374a-t10.tif(10)
 
image file: c6ra08374a-t11.tif(11)
Here A, A0, V, V0 are respectively the particle surface area and volume after and before deformation. κA is the surface area modulus and κV is the volume modulus. The elasticity moduli in each energy equations determine the energy scale and describe the stiffness of a particle. Elastic particle boundary is represented by Lagrange nodes n, which can be represented as Euler coordinate xn. The restoring force exerted on a boundary node is calculated as23,26
 
f(xn) = −∂W/∂xn. (12)

In IBM method,18 the force of a fluid node xf transferred from neighbor particle boundary nodes xn is represented as

 
image file: c6ra08374a-t12.tif(13)
 
Dirac delta function D(xnxf) = δ(xnxf)δ(ynyf)δ(znzf), (14)
 
image file: c6ra08374a-t13.tif(15)
which limits the fluid area under the interaction force computing consideration.

Using the no slip boundary condition, the velocity of the particle membrane node can be obtained from its neighbor fluid nodes. Fluid velocity uf(xf) is calculated from LBM and the solid node position xn can then be updated by ∂xn/∂t = up(xn). The particle nodes velocity is then determined as image file: c6ra08374a-t14.tif.

Symmetrical serpentine channel model

Fig. 1(a) sketches the geography of the symmetrical serpentine microchannel employed in this paper. The simulation area, as illustrated in Fig. 1(b), is a basic structure unit of the whole serpentine microchannel. The unit length is 280 μm. The channel width W is fixed at 80 μm and the height H is fixed at 40 μm, which makes the aspect ratio (AR) 0.5. The Reynolds number Re (Re = Udh/υ, where U is the characteristic velocity of the flow and dh = 2WH/(W + H) is the hydraulic diameter) varies in the simulation to study the particle focusing behaviors under different flow intensities in this channel unit. Simulation parameters were chosen as space step Δx = 1.25 × 10−6 m, time step Δt = 1.5 × 10−8 s and υ = 10−6 m2 s−1. Particle moduli were set as κS = 3.2 × 10−1 N m−1, κα = 3.2 × 10−1 N m−1, κB = 3.2 × 10−13 N m, κA = 3.2 × 10−2 N m−1 and κV = 3.2 × 104 N m−2, which are large enough to make the particles be considered as rigid ones. At the beginning of the simulation, six 5 μm particles and six 10 μm particles were released separately in the fluid field. Actuating force F drives the fluid field and the acceleration a from the actuating force imposed on the fluid node is obtained from the Poiseuille equation and can be simplified as a = 32υ2Re/dh3. Periodic boundary was applied to the inlet and outlet of the curved unit to simulate the particles flowing continuously in a serpentine microchannel.
image file: c6ra08374a-f1.tif
Fig. 1 (a) Sketch of the symmetrical serpentine microchannel and (b) curved microchannel unit for simulation.

Results and discussion

Simulation results

To study the particle focusing behaviors in the serpentine microchannel, particle motions and flow fields in the curved microchannel unit were simulated under Re = 25, 50, 75, and 100, separately. Fig. 2 depicts the time history of the vertical position of particles, 5 μm particles are marked in blue and 10 μm particles are marked in black. Particle flowing through a curved microchannel unit shows a wavy trajectory in the vertical direction. After passing through several curved microchannel units, the particle migrations trend to be stable. Vertical equilibrium positions of small particles oscillate at the heights of 10 μm and −10 μm, while larger particles flow closer to the vertical center slightly. However, there are a few particles staying in the vertical center plane on account of the vertical symmetry of the microchannel. The appearance of these trajectories is determined by the initial position of the particles and the particles initially located near the vertical center plane are prone to stay there, which may be unlikely events for large particles and minor events for small particles under a real condition. The trajectories of these particles will not be considered in the following simulation analyses. We chose the last period of simulation time to draw the top view of the stacked particle trajectories, as illustrated in Fig. 3. Time interval between each pair of neighbor particles is 9 × 10−6 s and the sequence of particles gives an idea on how fast the particle goes through the curved unit. It is found that both the large and small particles are focused in a perfect single line when Re = 75. Either higher or lower flow velocity makes the differently sized particles flow in their own trajectories. When Re increases from 25 to 50, the separated trajectories of small particles along the sidewalls begin to migrate towards the midline of the channel, and finally merge to form a single curved line at Re = 75. The migration of large particles to the channel center is much faster as compared with that of small particles. At Re = 100, the trajectories of some small particles deflect from the focused trajectories of large particles and separate into symmetrical lines again.
image file: c6ra08374a-f2.tif
Fig. 2 Time history of vertical position of the particles under (a) Re = 25, (b) Re = 50, (c) Re = 75 and (d) Re = 100 in a symmetrical serpentine microchannel. The number of the released 5 μm particles at Re = 100 was doubled due to the complicated focusing pattern of small particles at high flow intensities.

image file: c6ra08374a-f3.tif
Fig. 3 Particle trajectories from the top view under (a) Re = 25, (b) Re = 50, (c) Re = 75 and (d) Re = 100.

The cross-section views of the stacked particle trajectories were also drawn to find the equilibrium positions as illustrated in Fig. 4. Specifically, at Re = 25, small particles migrate close to the sidewalls with circle trajectories in the cross-section while large particles focus along the sidewalls with their equilibrium positions closer to the channel center, as compared with the small ones. There are four clear equilibrium rings in the cross-section for small and large particles due to the vertical and horizontal symmetry of the fluid field in the serpentine microchannel at low flow intensities. When Re = 50, particles focus close to the horizontal center of the channel. At Re = 75, the observed four equilibrium positions further reduce to two and the particle trajectories in the cross-section appear as “V” shapes. The symmetrical “V” leads the differently sized particles to be aligned from the top view. However, when the Re increases to 100, trajectories of large particles are stretched in the cross-section and small particles are defocused into four equilibrium positions once more. Compared with the large particles, small particles are more easily dragged by the strong effects of Dean flow with two counter-rotating vortices, and the direction alternating Dean flow also focus the small particles at two symmetrical trajectories from the top view. The simulation results of Reynolds number-focusing position agree well with the bifurcation focusing pattern in symmetrical serpentine channel from Gossett and Di Carlo's experiments.10


image file: c6ra08374a-f4.tif
Fig. 4 Particle trajectories in the cross-section under (a) Re = 25, (b) Re = 50, (c) Re = 75 and (d) Re = 100. The directions of small particle trajectories in (a) are shown by orange arrows.

In order to explain the above particle focusing phenomena, we quote the magnitude ratio of inertial lift force FL to Dean drag force FD,9,27

 
image file: c6ra08374a-t15.tif(16)
to show the relationship between two forces. Curvature ratio δ = dh/2r and r is the channel radius. Particle diameter d is in proportion to the FL. With increasing Re, the FD will become the dominant force in a curved channel.

The competition of FL and FD decides the focusing pattern of various sized particles. We drew the trajectories of a pair of particles in symmetric positions and the fluid field around the lower particle (see Fig. 5) to study the basic rule of particle focusing. The serpentine channel is divided into four sections to analyze particle migration. As cross-referenced in Fig. 4, at Re = 25, small particles flowing through Section 1 are mainly driven by the Dean flow and move close to the sidewall. When the small particles enter into Section 2, the Dean flow direction turns opposite, but the high flow intensity here strengthen the FL, which slows down the rate of particle going towards the vertical center plane and keeps the particle flowing at the corner of the channel. In Sections 3 and 4, FD plays the dominant role again and the particles flow towards the vertical center plane first and flow in the opposite direction. When Re increases to 50 and 75, small particles migrate to the center of the channel under FD, which swings the particles horizontally in the cross-section. FL here has relatively slight influence on the particle migration in vertical direction. At Re = 100, the strong Dean flow drags the small particles to the sidewall again, and the FL in the high flow intensity area has obvious effects on the small particles.


image file: c6ra08374a-f5.tif
Fig. 5 The particle trajectories and rotation in the serpentine channel. (a) 5 μm particle trajectories in the serpentine channel and time interval between each pair of neighbor small particles is 9 × 10−6 s. (b) The rotation of small particles under various channel Re numbers. (c) 10 μm particle trajectories in the serpentine channel with double time interval in (a). (d) The rotation of large particles under various channel Re numbers. In (b) and (d), the blue lines are the particle center trajectories and the red lines are the trajectories of a fixed node on the particle surface. Small particles in (b) are stretched twice their size for better differentiation.

Compared with the small particles, large particles in Fig. 5(c) have similar behaviors in the serpentine channel. However, along with the increase of Re, the side to center migration processes of large particles are much faster. When the Re is up to 100, large particles are more difficult to be separated into multi-lines than small particles. At low flow rates, particles flow along two sidewalls under the mixing effects of FL and FD28 and focus at channel midline by the strong direction alternating FD29 at high flow rates. However, the size difference also plays an important role for focusing pattern, and it is found that the large particles are more likely to focus at the channel horizontal center. From eqn (16), we found that the small particles are much easier to be influenced by FD, which explains why the small particles separate easily under a high Re. Considering the volume effects of large particles in a low AR channel, rotation of the particle in the Dean vortex is a possible explanation for the faster side to center migration. Large particles are more likely to be rotated by the counter-direction flow of the Dean vortex and sucked in it. When the flow intensity is weak, the rotation strengthens the Dean flow effects on the large particle and focus them closer to the channel center. The particle rotation here can be defined by the vector p. This vector p starts from the particle center and points to a fixed node on the particle surface. φ is defined as the angle between the projection of p on the channel cross-section plane and the z-axis. The rotations of the small and large particles are shown in Fig. 5(b) and (d). Blue lines are the trajectories of particle centers and the red lines are the trajectories of the fixed nodes. The rotation direction of the small particles has no obvious rules, but the large particles, meanwhile, rotate with the direction of Dean flow vortex as expected. The rotation of the large particles is also strengthened with the increase of Re. When the Re is higher, the rotation directions of large particles agree better with the Dean flow vortex.

To understand the force and rotation of particles deeper, particles and the fluid field in the channel cross-section A–A marked in Fig. 5(c) are extracted and illustrated in Fig. 6 to study the interaction between particles and the flow around. It can be concluded that at Re = 25, the 5 μm particle near the corner is balanced under FL and FD, which are in opposite directions. The FL acted on the small particle points to the channel corner, which holds the particle focusing close to the sidewall. In the case of 10 μm particles, there is more interaction between the particle and fluid, and the particles are too large to balance as the small particles close to the sidewall. The Dean flow easily rotates the large particles horizontally with its counter-direction flow near the Dean flow vortex and therefore the large particles focus closer to the channel center than small ones. From the eqn (16), the larger particle has heavier FL influence and that is why the large particles usually flow closer to the inner wall of a spiral channel.30 What is interesting is that smaller particles here flow closer to the sidewalls in the serpentine channel, which means that there should be more consideration when analyzing particle behaviors in complicated fluid field. Due to the volume effects, different sized particles may flow under different fluid field environment in the same channel.


image file: c6ra08374a-f6.tif
Fig. 6 The cross-section view of the fluid field in the Fig. 5 A–A section and tangent vectors around 5 μm and 10 μm diameter particles.

When Re increases up to 50, 75 and 100, the Dean flow in the serpentine microchannel is stronger and particles are driven left and right by the direction alternating Dean flow. Even if at a higher Re, the FL here is still in a considerable quantity. The FL leads the particles to flow towards the upper and bottom walls, which forms the bottom of the “V” trajectory in the cross-section at Re = 75. Large particles rotate in the Dean flow vortex center, which strengthens the particle focusing performance. Instead, the 5 μm particle is too small to rotate with the Dean flow, which makes the focusing of small particles in the serpentine channel more unstable.

The Dean flow is the main motivation for the particles focusing. To analyze the Dean flow vortex's effects on the particles migration in the whole serpentine channel, the vertical velocity fluid field in z-direction is drawn to show the Dean flow distribution from the top view. Fig. 7 shows typical upper Dean vortex flow fields in a serpentine channel. Positive value means that the flow rises up and negative value shows that the flow goes down. The vertical rotation flow around the small particle is also shown in the figure. The fluid velocity gratitude of the main flow rotates the particle and the rotation flow is stronger than the Dean flow when the flow intensity is weak. However, at Re = 25, small particles focus at the vortex center when passing the inner wall side in a curved channel. As the flow intensity in the serpentine channel is stronger, particles begin to focus near the Dean flow vortex center in the whole serpentine channel, which means that the Dean flow plays the dominant role for particle focusing in this situation.


image file: c6ra08374a-f7.tif
Fig. 7 The flow velocity field in vertical direction of upper Dean flow vortex in the serpentine channel.

From the analysis above, the effect of FD is to regularly swing the particles in the cross-section of the serpentine channel and the FL which is relatively strong at lower flow intensity is more like a trap to keep the particle flow near the sidewalls. Stronger flow intensity frees the particle from the FL trap and intensifies the swing effects of FD. Large particles are more likely to rotate in the cross-section with Dean flow vortex, which makes the focusing of large particles is more stable than small ones. In the focusing process, the equilibrium positions are dynamically changing with the fluid condition, not fixed at a point in the cross-section as the equilibrium positions in a straight microchannel.

Experimental validation

The 40 μm high and 80 μm wide serpentine channel used in our experiment has the same dimension with the model illustrated in Fig. 1 and the whole serpentine channel is constructed with 18 curved channel units. The symmetrical serpentine channel device was fabricated by the standard soft-lithography. Specifically, the photomask was printed on a transparency film using a high resolution laser photoplotter (RP525-SST, EIE) and aligned to the silicon wafer covered with a layer of photoresist (SU-8 2050, MicroChem Corp.). After exposing using an ultraviolet exposure machine (MJB4, SUSS) and careful development, a SU-8 master mold could be obtained. The channel structure was created by molding with polydimethylsiloxane (PDMS). Then, the PDMS blockage consisted with an open channel was bonded with a glass slide.

Polystyrene fluorescent particles with average diameters of 4.8 μm (G0500, Thermo Fisher Scientific Inc.) and 9.9 μm (G1000, Thermo Fisher Scientific Inc.) were analyzed in this experiment. The 4.8 μm particle sample of 1% solid content was diluted 200 times with phosphate-buffered saline (PBS) solution (Sigma-Aldrich) and the 9.9 μm particle sample of 1% solid content was diluted 40 times. Tween 20 (Sigma-Aldrich) was added to prevent the particle from aggregation. Mixed particle suspension with different particle sizes was the mixture of equal volume of two sized particle suspension. The ethyl alcohol with a low surface energy was loaded into the symmetrical serpentine channel firstly before the particle focusing experiments to clean the channel walls and prevent the bubbles from producing to disturb the particle flowing. In the particle focusing experiment, the flow in serpentine microchannel was driven by a micro-injection pump (KDS270, KD Scientific Inc.) and observed by a fluorescence inverse microscope (IX71, Olympus). 10× objective lens (UPLFLN, 10×/0.3) and fluorescence excitation module were applied to observe the particle behavior in light field and dark field. The equipped 14 bit high-speed CCD camera (Exi Blue, Qimaging) and software IMAGE-PRO EXPRESS were used to capture the particle behaviors as video sequences.

To avoid the random factor, multiple image frames under the same operating condition were stacked to obtain the statistical particle distribution during a certain period of time. The video files were loaded to free ImageJ software (version 1.48v, NIH) and the image sequences were stacked by “Sum Slices” method to obtain the composite fluorescence image illustrating the particle trajectories. As for the bright field images, the obtained sequences were stacked by “Min Intensity” method to obtain the statistical particle distribution across the channel width.

Particle focusing patterns were recorded under increasing flow intensities in the experiment and the results are shown in Fig. 8. When the volume flow rate Q is 80 μL min−1 (Re ≈ 22.2), both the 4.8 μm and 9.9 μm particles focus into two symmetrical curved lines along the microchannel sidewalls and the space between the trajectories of large particles is narrower than that of small particles. At Q = 180 μL min−1 (Re ≈ 50), 9.9 μm particles focus to the channel center, but the 4.8 μm particles obviously focus in two lines at the microchannel center and there are some small particles still remain in the close wall area. It can be seen that before all of the small particles focus into a line in the center of the channel, the particles in center already begin to be separated by the Dean flow. The small particles near the center are easily controlled by the Dean flow and separate into two symmetrical trajectories and others close to the walls are still under the influence of FL. Complicated flow condition in the serpentine channel makes different dominant forces act on the small particles at different positions. With the stronger Dean flow at increased Re, the small particles begin to focus in clear two curved lines (Q = 220 μL min−1, Re ≈ 61.1), and the two symmetric focused trajectories of large particles appear when the Re is higher (Q = 300 μL min−1, Re ≈ 83.3). The Dean flow has more obvious influence on the small particles. However, the focusing performance of large particles is much better, which perhaps can be ascribed to the rotation in the Dean vortex as predicted in the simulation. At a high Re, the rotation of large particles can absorb parts of the energy that drags the particle towards the sidewalls. Because of the gravitation to the Dean vortex center, the large particles have a more stable focusing performance and flow closer to the channel center than small particles. Some stubborn small particles here are still unfocused and form the background of the focused trajectories.


image file: c6ra08374a-f8.tif
Fig. 8 (a) Microscope photographs of focusing patterns of 4.8 μm and 9.9 μm particles under different flow intensities and (b) the comparison of particle distributions, obtained by experiments and simulations. (c) Microscope photographs of mixed particles near the serpentine channel outlet under different flow intensities.

To compare the particle distribution of simulation and experiments, we draw the particle focusing positions in Fig. 8(b). The normalized fluorescence intensity data are extracted from Fig. 8(a) cross-section B–B. Since the relationship between driving force and Re was calculated form the Poiseuille's law, which is defined for a straight channel, the Re in the simulation is higher than the actual value. In the particle distribution comparison of simulation and experiments, the Re of the simulation has been divided by 1.5. For large particles, the simulation focusing positions match the experimental results well as a bifurcation pattern. The small particles in the experiment have more complicated patterns. Before small particles appear at the channel center, they focus well along the sidewalls. When majority of small particles focus at the channel center at higher Re, there is still a small part of them focus close to the sidewalls. These small particles may still be trapped by the FL when the flow intensity is relatively low (Q = 180 μL min−1). When the FD is the dominant force in the channel (Q = 220 μL min−1, Q = 300 μL min−1), these small particles may flow at the channel vertical center plane and are pushed by the Dean flow further to the sidewalls.

Dark and bright field photographs of the mixture of differently sized particles focused at serpentine microchannel outlet are shown in Fig. 8(c). As the prediction in the simulation, the large particles focus closer to the channel center when the flow intensity is slow (Q = 80 μL min−1). At Q = 140 μL min−1 (Re ≈ 38.9), the large particles focus in the channel center first and the small ones still in two curved lines aside. This time difference of different sized particles focusing is useful for the particle separation. Here FD is the relatively strong force for the large particles which are trapped in the Dean flow vortex and the small particles are still controlled by the balance of FL and FD. When Q = 180 μL min−1, mixed particles focus in a trajectory and some small particles still flow near the sidewalls. However, the separation of differently sized particles is not obvious under the stronger Dean flow (Q = 220 μL min−1) because their trajectories are close and the disturbance of the defocused small particles is difficult to avoid.

Conclusions

In this paper, particle focusing behaviors in a symmetrical serpentine channel are studied and explored through numerical simulation and experimental approaches. The competition of FL and FD acting on the particles are illuminated clearly. Lower Re strengthens the influence of FL and the balance of FL and FD focuses the particles close to the sidewalls. The small particles are found to focus closer to the sidewalls than the large particles, which means that the large particles are not always under the heavy influence of FL and they are more prone to rotate in the Dean vortex center in a low AR confined microchannel. A reasonable higher Re intensifies the FD, and the Dean flow swings particles out of the trap of FL and focuses particles in the horizontal center of the channel near the Dean flow vortex center. In addition, large particles are inclined to rotate with the Dean flow vortex and flow to the channel center before the small particles, which makes the serpentine channel a good candidate for the separation of differently sized particles. When the Re is high enough, the Dean flow plays the absolute dominant role and the particles focus in two symmetrical trajectories from the top view again. This hydrodynamic control of micro-particles promises to be an effective and efficient manipulation method for particle focusing and separation.

Acknowledgements

This research work is supported by the National Natural Science Foundation of China (51505082 and 51375089), the Natural Science Foundation of Jiangsu Province (BK20150606), the “333” Project of Jiangsu Province (BRA2015291) and the Open Foundation of the State Key Laboratory of Fluid Power and Mechatronic Systems (GZKF-201501).

References

  1. D. Di Carlo, Lab Chip, 2009, 9, 3038–3046 RSC.
  2. G. B. Lee, C. C. Chang, S. B. Huang and R. J. Yang, J. Micromech. Microeng., 2006, 16, 1024 CrossRef.
  3. I. Doh and Y. H. Cho, Sens. Actuators, A, 2005, 121, 59–65 CrossRef CAS.
  4. S. C. Hur, H. T. K. Tse and D. Di Carlo, Lab Chip, 2010, 10, 274–280 RSC.
  5. E. Sollier, D. E. Go, J. Che, D. R. Gossett, S. O'Byrne, W. M. Weaver, N. Kummer, M. Rettig, J. Goldman, N. Nickols, S. McCloskey, R. P. Kulkarni and D. Di Carlo, Lab Chip, 2014, 14, 63–77 RSC.
  6. S. S. Kuntaegowdanahalli, A. A. S. Bhagat, G. Kumar and I. Papautsky, Lab Chip, 2009, 9, 2973–2980 RSC.
  7. S. C. Hur, S. E. Choi, S. Kwon and D. D. Carlo, Appl. Phys. Lett., 2011, 99, 044101 CrossRef.
  8. C. Jin, S. M. McFaul, S. P. Duffy, X. Deng, P. Tavassoli, P. C. Black and H. Ma, Lab Chip, 2014, 14, 32–44 RSC.
  9. D. Di Carlo, D. Irimia, R. G. Tompkins and M. Toner, Proc. Natl. Acad. Sci. U. S. A., 2007, 104, 18892–18897 CrossRef CAS PubMed.
  10. D. R. Gossett and D. D. Carlo, Anal. Chem., 2009, 81, 8459–8465 CrossRef CAS PubMed.
  11. D. Di Carlo, J. F. Edd, D. Irimia, R. G. Tompkins and M. Toner, Anal. Chem., 2008, 80, 2204–2211 CrossRef CAS PubMed.
  12. J. Zhang, S. Yan, R. Sluyter, W. Li, G. Alici and N. T. Nguyen, Sci. Rep., 2014, 4, 4527 Search PubMed.
  13. S. Chen and G. D. Doolen, Annu. Rev. Fluid Mech., 1998, 30, 329–364 CrossRef.
  14. D. Kandhai, D. Vidal, A. Hoekstra, H. Hoefsloot, P. Iedema and P. Sloot, Int. J. Numer. Methods Fluids, 1999, 31, 1019–1033 CrossRef CAS.
  15. A. J. Ladd, J. Fluid Mech., 1994, 211, 285–309 CrossRef.
  16. A. J. Ladd, J. Fluid Mech., 1994, 271, 311–339 CrossRef CAS.
  17. Z. G. Feng and E. E. Michaelides, J. Comput. Phys., 2004, 195, 602–628 CrossRef.
  18. C. S. Peskin, Acta Numerica, 2002, 11, 479–517 CrossRef.
  19. Y. Qian, D. d'Humières and P. Lallemand, Europhys. Lett., 1992, 17, 479 CrossRef.
  20. Z. Guo, C. Zheng and B. Shi, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2002, 65, 046308 CrossRef PubMed.
  21. Y. Sui, Y. T. Chew, P. Roy and H. T. Low, J. Comput. Phys., 2008, 227, 6351–6371 CrossRef.
  22. D. K. Sun and Z. Bo, Int. J. Heat Mass Transfer, 2015, 80, 139–149 CrossRef.
  23. T. Krüger, Computer simulation study of collective phenomena in dense suspensions of red blood cells under shear, Vieweg+Teubner Verlag, Wiesbaden, 2012 Search PubMed.
  24. R. Skalak, A. Tozeren, R. P. Zarda and S. Chien, Biophys. J., 1973, 13, 245–264 CrossRef CAS PubMed.
  25. K. Tsubota and S. Wada, Int. J. Mech. Sci., 2010, 52, 356–364 CrossRef.
  26. J. Charrier, S. Shrivastava and R. Wu, J. Strain Anal. Eng. Des., 1989, 24, 55–74 CrossRef.
  27. J. Matas, J. Morris and E. Guazzelli, Oil Gas Sci. Technol., 2004, 59, 59–70 Search PubMed.
  28. J. Zhang, S. Yan, W. Li, G. Alici and N.-T. Nguyen, RSC Adv., 2014, 4, 33149–33159 RSC.
  29. J. Zhang, W. Li, M. Li, G. Alici and N.-T. Nguyen, Microfluid. Nanofluid., 2014, 17, 305–316 CrossRef CAS.
  30. R. Aman, K. G. Amit, N. Sunitha, C. Dino Di, F. E. Jon and T. Mehmet, New J. Phys., 2009, 11, 075025 CrossRef PubMed.

Footnote

These authors contributed equally to this work.

This journal is © The Royal Society of Chemistry 2016
Click here to see how this site uses Cookies. View our privacy policy here.