A theoretica anaysis of the resoution due to diffusion and size-dispersion of partices in deterministic atera dispacement devices arxiv:7.347v [physics.fu-dyn] 2 Nov 27 Martin Heer and Henrik Bruus MIC Department of Micro and Nanotechnoogy, Technica University of Denmark DTU Buiding 345 East, DK-28 Kongens Lyngby, Denmark (Dated: 2 November 27) We present a mode incuding diffusion and partice-size dispersion for separation of partices in deterministic atera dispacement devices aso known as bumper arrays. We determine the upper critica diameter for diffusion-dominated motion and the ower critica diameter for pure convection-induced dispacement. Our mode expains the systematic deviation, observed for sma partices in severa experiments, from the critica diameter for separation given by simpe aminar fow considerations. I. INTRODUCTION In 24 Huang et a. [] deveoped the eegant method of partice separation by deterministic atera dispacement in so-caed microfuidic bumper arrays. The method, which reies on the aminar fow properties characteristic of microfuidcs, shows a great potentia for fast and accurate separation of partices on the micrometer scae [, 2, 3, 4]. Among the key assets of the deterministic atera dispacement separation principe are that cogging can be avoided because partices much smaer than the feature size of the devices can be separated, that the devices are passive, i.e. the partices bump into soid obstaces or bumpers, and that the separation process is continuous. More precisey, partice transport in microfuidic bumper arrays is primariy governed by convection due to the fuid fow and dispacement due to interaction with the bumpers in the array []. These processes are deterministic and the critica diameter for separation of reativey arge partices in these devices is we understood in terms of the width of fow anes bifurcating around the bumpers in the periodic arrays [3]. However, if bumper arrays and partices are scaed down, diffusion wi infuence the separation process and affect the critica partice size significanty. Previousy reported data on separation of partices in bumper arrays a show a bias towards arger critica partice size than that given by the width of the fow anes nearest to the bumpers of the array [, 2, 3]. In this work we extend existing modes by adding diffusion and taking partice-diameter dispersion into account, and thereby expain the observed discrepancy. In bumper arrays partices are convected by the fuid fow through an array of bumpers paced in coumns separated by the distance λ in the fow direction, see Fig. (a). For a given integer N, The array is made N-periodic in the fow direction by dispacing the bumpers in a given coumn a distance λ/n perpendicuar to the fow direc- (a) III g II Fow 2 3 2 3 2 3 2 3 2 3 2 IV III λ λ I N II λ 2 3 4 c (b) w (3) w (2) w () w () FIG. : (a) An array of bumpers (back dots) with the definition of the ane number (sma arabic numbers), the coumn number c (arge arabic numbers), and the gap number g (roman numbers). The shaded region iustrates how the shift in the position of gap II foows the geometry of the array. (b) Cose-up of a singe gap between two bumpers (disks) in the array. Each of the four fow anes carries the same fow rate. Due to the paraboic fow profie in the gap region, the width w () of fow-ane depends on the position in the gap. tion with respect to the bumper positions in the previous coumn. Due to this periodicity of the array and the aminarity of the fow, the stream can naturay be divided into N anes, each carrying the same amount of fuid fux, and each having a specific path through the device, see Ref. []. For a given steady pressure drop, the fuid in the device moves with an average veocity u. Assuming a paraboic veocity profie u(x) in the gap of width w g between two neighboring bumpers, see Fig. (b), ) x u(x) = 6u ( xwg, () w g the tota fow rate Q tot is given by Q tot = wg u(x)dx = w g u. (2) URL: http://www.mic.dtu.dk/microfuidics For an N-periodic array, the N fow anes in a given gap carry the same fow rate Q tot /N. The width w () of ane
2 is found by soving Q x () +w () tot N = u(x)dx, (3) x ().8 Mode based on data from Larsen et a. Huang et a. Ingis et a. Width of first ane where x () = j= w(j) is the starting position of ane. In the simpe bifurcating fow-ane mode [, 3] the critica diameter d c is given as d c /2 = w (). A sma partice with d < d c wi never eave its initia fow ane and wi thus be convected in the genera fow direction foowing a so-caed zigzag path. Large partices with d > d c wi quicky bump against a bumper and from then on be forced by consecutive bumping to foow the skew direction of the array geometry, the so-caed dispacement path. When a partice gets bumped by a bumper in the array it wi be dispaced perpendicuar to the fow direction unti its center is ocated one partice radius d/2 from the surface of the bumper. This corresponds to n anes of dispacement, n = N w g u d/2 u(x)dx = N d2 4w 2 g ( 3 d ). (4) w In the buk fuid, where the anes are assumed to have equa width λ/n, see Fig. (a), the dispaced distance disp is therefore disp (d) = n λ N = λ d2 4w 2 g (3 dwg ). (5) In this work we extend the simpe bifurcating fowane mode by incuding diffusion and partice-diameter dispersion. II. MODEL INCLUDING DIFFUSION During the average time τ = λ/u it takes a partice to move by convection from one coumn to the next, it aso diffuses. We assume that the diffusion process perpendicuar to the fow direction is normay distributed with mean vaue zero and variance σ 2 = 2Dτ, (6) where the diffusivity D is given by the Stokes Einstein expression D = k BT 3πηd. (7) Here k B is Botzmann s constant, T is the temperature and η is the viscosity of the fuid. A. Diffusion mode In order to escape bumping, a partice must in the time interva τ diffuse more than the difference disp λ/n d/wg.6.4.2 Experimenta data points Bump, Larsen et a. Mixed, Larsen et a. Zig-zag, Larsen et a. Bump, Huang et a. Zig-zag, Huang et a. Bump, Ingis et a. Zig-zag, Ingis et a...2.3.4.5 /N FIG. 2: Our mode appied to the experimenta data of Ingis et a. [3], Huang et a. [] and Larsen [5]. Partice diameter d over gap width w g is potted as a function of the inverse period /N. The back ines show the theoreticay predicted critica partice size for the bumper arrays used by Ingis et a. [3], Huang et a. [] and Larsen [5], respectivey. The experimenta data points are representing partices foowing bumper paths (open symbos), zigzag paths (soid back symbos), and neither of these paths (soid gray symbos). between the buk dispacement and the shift in position of the next bumper. The probabiity p esc for this to happen is given in terms of the error function as p esc (d) = 2 ( 2 erf disp (d) N λ ), (8) 2σ(d) where we have introduced the d-dependence expicity. When a partice is transported through a bumper array it must bump at every bumper within one period of the array in order to be dispaced one gap at the outet. Thus, if the partice escapes at east one time in N attempts, it wi not be dispaced. We define the critica partice size d c to be the size for which haf of the partices escape bumping as they pass one period of the array. Thus d c can be found by soving N k= N p k esc (d c ) k[ ] N k p esc (d c ) = 2. (9) In Fig. 2 we have potted the critica partice size as a function of the bumper period for parameter vaues corresponding to the bumper arrays used by Huang et a. [] (dashed ine), Ingis et a. [3] (dotted ine), and Larsen [5] (fu ine).
3 B. Comparison to experiments The observation that the critica partice size in a bumper device is arger than the width of the first fow ane is aso supported by the experimenta data in Fig. 2 of Ref. [3]. Our mode suggests that the deviation of the critica partice size from the width of the first fow ane can be expained by diffusion of the partices. In Fig. 2 above, it is seen how we the ines are predicting the transition between zigzag paths and dispacement paths: the thick ine is dividing open and cosed trianges, the dashed ine is dividing open and cosed squares, and, to a esser extent, the dotted ine is dividing open and cosed circes. Using parameter vaues corresponding to the bumper device presented by Huang et a. [] our mode predicts a critica partice diameter of.45 times the width of the gap for the partices traveing through the device with an average veocity of 4µm/s. This is in good agreement with Fig. 2(a) in Ref. []. III. A DISCRETE MODEL INCLUDING DIFFUSION AND DISPERSION The partices used in the experiments on separation of partices are not mono-disperse. Their average diameters are distributed around a certain mean vaue with a reative standard deviation d/d, which typicay is 2 %, % and 5 % for partices with d = 25nm, d = nm and d = 5 nm, respectivey. In the foowing we introduce a discrete mode of the transport of partices with different diameters d in an N- periodic bumper array taking convection, diffusion and size-dispersion into account. The mode aows us to study the reative infuence of a three phenomena on the separation efficiency in a fast and simpe manner. In particuar our resuts suggest that the critica size for separation or dispacement, studied above, must be suppemented by a smaer critica size beow which pure diffusion governs the motion of the partices in the bumper array. This prediction has not yet been tested experimentay. A. Definition of the discrete mode At any instant, a partice is assumed to be positioned in the center of a fow ane situated in a gap g in some coumn c of the array. For simpicity we further assume that the size distribution of any given set of partices is a norma distribution with a mean vaue given by the size quoted by the manufacturer and a reative standard deviation of %. By convection any given partice moves from one coumn to the next. If it ends up in the ast ane ( = N ) in one gap, it wi be shifted to the first ane ( = ) in the subsequent gap. Otherwise it wi stay in the current gap and move up one ane. In our mode pure convection is therefore described by the discrete map { (c +, g +, ), if = N, (c, g, ) () (c +, g, + ), otherwise. Because of the finite diameter d of the partice there is a minimum and a maximum ane number that it can occupy. The minimum ane number is the smaest integer min that satisfies min w () > d 2. (a) = Simiary, the maximum ane number max is the argest integer that satisfies N = max w () > d 2. (b) Consequenty, the simpe convection mapping from Eq. () needs to be modified to account for the finite size of the partices (c +, g +, min ), if = N, (c, g, ) (c +, g, + ), if < max, (2) (c +, g, max ), otherwise. The above convection scheme accounts for the separation of partices in deterministic atera dispacement devices according to size. The critica partice radius predicted by this mode is d c = 2w () (3) in accordance with the geometric arguments of Ref. [3]. To characterize the separation quantitativey, we define the reative partice numbers r, r and r d as r = reative number of partices foowing the zigzag path. r = reative number of partices foowing the dispacement path. r d = reative number of a other partices. (4a) (4b) (4c) With these definitions the sum r +r +r d is aways unity. If r = a partices foow the zigzag path and if r = a partices foow the dispacement path. If r d some of the partices end up at positions not expained by the deterministic anaysis of the separation process. In Fig. 3 we have potted the reative partice numbers r, r, and r d as a function of average partice diameter d. B. Pure mono-disperse convection For mono-disperse and non-diffusing partices, the mode resuts, as expected, in two modes: the zigzag
4 r r r d.5 diff on, disp on diff off, disp on diff off, disp off diff on, disp off 2 3 4 5 6 7 d/d c.5 diff on, disp on diff off, disp on diff off, disp off diff on, disp off 2 3 4 5 6 7 d/d c.5 diff on, disp on diff off, disp on diff off, disp off diff on, disp off 2 3 4 5 6 7 d/d c FIG. 3: Reative number r, r and r d of partices foowing the zigzag path, the dispacement path, and neither of these two paths, respectivey, potted as a function of the normaized, average partice diameter d/d c. The parameters of the bumper array is taken form Ref. [5]: N =, λ = 8 µm, w g = µm, L = 2Nλ = 6 mm, and u = 25 µm/s. The buffer iquid is water at room temperature. Negecting diffusion (soid symbos), the partices foow the zigzag path if d < d c = 8 nm. If d > d c they foow the dispacement path. Incuding diffusion (open symbos), the sma partices are not infuenced by the bumpers. For d > 2. dc the infuence of the bumpers set in, and for d > 3.4 dc the partices foow the dispacement path. The thick back vertica ines indicate the partice size when sma partices stop behaving purey diffusive (eft-most ines) and when arge partices begin a purey deterministic dispacement (right-most ines). mode and the dispacement mode, see the cosed circes in Fig. 3. For d < d c we have r =, and for d d c we have r =, whie we aways have r d =. The reative partice numbers can therefore be written as by the array geometry r = dc (s d) 2 exp 2π( d) 2 2( d) 2 ds. (6a) (r, r d, r ) = { (,, ) for d < dc (,, ) for d d c. (5) C. Infuence of size dispersion If we assume that the partices are not mono-disperse, but are distributed around a mean size d with standard deviation d, the shift as a function of d from the zigzag mode to the dispacement mode happens graduay instead of abrupty at a certain critica size d c (Fig. 3, cosed squares). The reative number of partices foowing the zigzag path r can be found by integrating over a partice sizes smaer than the critica diameter given Simiary, the reative number of partices foowing the dispacement path r can be found by integrating over a partice sizes arger than d c r = (s d) 2 exp 2π( d) 2 2( d) 2 ds. d c (6b) The system is sti a bi-moda system because r d = for a partice sizes. D. Infuence of diffusion In D during the time step τ a partice diffuses the distance, the average of which is the size-dependent
5 diffusion ength σ(d) given by σ(d) = = 2Dτ = 2 k BT λ. (7) 3πηd u In our mode we discretize the transverse diffusion as the propery rounded number of fow anes crossed by the partice during diffusion, n jump = N λ. (8) The addition of diffusion smears out the dispacement of the partices and causes the critica radius to be arger than in the diffusion-ess case (Fig. 3, open symbos).. Bumping criterion for sma partices Very sma partices are competey dominated by diffusion, and the partice distribution at the end of the array is simpy given by the transverse diffusion of the partices during the time L/u it takes for the partice to be convected a the way L through the array. For sma partices we therefore have r = 3 λ 2 N 3 λ 2 N x 2 exp 2πσ 2 2σ 2 dx, (9) where σ 2 = 2DL/u. As the partice diameter d is increased, the bumpers begin to become important as the diffusion ength σ(d), Eq. (7), is decreased and becomes equa to the dispacement ength disp, Eq. (5). Using the criterion σ(d ) = disp, with the parameter vaues used in Fig. 3, we find that partices stop behaving as sma diffusion-dominated partices and start interacting with the bumpers when d = 2. d c. This cross-over vaue is indicated by the eft-most vertica ines in Fig. 3, and it fits we with the simuation data. 2. Bumping criterion for arge partices Large partices wi interact with the bumpers at every row in the array and their position is thus reset at every bump to be disp. Diffusion can therefore be negected for such partices, they a foow the dispacement path, and r =. As the partice diameter is owered, the probabiity p esc that a partice escapes the dispacement path can be estimated by the probabiity of diffusing from the dispaced position to the ast fow ane in the gap, i.e. the distance disp λ N. This probabiity p esc is given by Eq. (8). In order to foow the dispacement path, a partice must bump at each row in the array. If we consider a N- periodic array with m fu periods, the partices wi have mn bumping opportunities as they pass through the entire array. If a partice evades bumping at a bumper, it wi be convected by the fow through a fu period of the array before bumping is possibe again. Because of the escape, it wi miss N bumping opportunities and end up one gap from the dispacement path. Consequenty, if a partice escapes one time, it wi ony have (m )N bumping opportunities and has therefore escaped bumping with a probabiity of /[(m )N]. The upper critica partice size d 2 for convection-induced dispacement is defined using Eq. (8) as p esc (d 2 ) = (m )N. (2) For the device used in the experiments by Larsen [5] we find d 2 = 3.4d c. The predicted upper imit for diffusion dominated motion and the ower imit for convectioninduced dispacement compares we with the experimenta observation by Larsen [5] that some partices end up between the dispacement path and the zigzag path (Fig. 2, gray trianges). IV. CONCLUSION Experimenta data on separation of partices in bumper arrays a show a systematic deviation from predictions made from the bifurcating fow ane mode [, 3, 5]. By adding diffusion to the bifurcating fow ane mode, we expain the discrepancy. In addition, we have suggested a simpe discrete mode for simuating partice separation in bumper arrays. The presented mode takes partice diffusion and size dispersion into account and has been vaidated against experimenta data for a bumper device with period N =. The transition from zigzag path to dispacement path happens at a partice size 2. to 3.4 times arger than the critica partice size predicted from geometrica arguments. This is in correspondence with the experimenta observations from Larsen [5]. Our discrete mode and the estimates presented in this paper suggest that partices of twice the size of the geometrica critica size of the N = bumper device behave diffusivey and are not affected by the bumpers because the sma diffusive partices rarey come into contact with the bumpers due to random Browninan motion. V. ACKNOWLEDGEMENT We thank Asger Vig Larsen, Anders Kristensen, Jason Beech, and Jonas Tegenfedt for inspiring discussions on experimenta issues. Martin Heer was supported by the Danish Research Counci for Technoogy and Production Sciences grant no. 26-4-74.
6 [] L. R. Huang, E. C. Cox, R. H. Austin, and J. C. Sturm, Science 34, 987 (24). [2] S. Zheng, R. Yung, Y.-C. Tai, and H. Kasdan, in 8th IEEE Internationa Conference on Micro Eectro Mechanica Systems (Miami, USA, 25), pp. 85 854. [3] D. W. Ingis, J. A. Davis, R. H. Austin, and J. C. Sturm, Lab Chip 6, 655 (26). [4] J. Davis, D. Ingis, K. Morton, D. Lawrence, L. Huang, S. Chou, J. Sturm, R. Austin, PNAS 3, 4779 (26). [5] A. V. Larsen, MSc thesis, MIC Department of Micro and Nanotechnoogy, Technica University of Denmark (26), URL www.mic.dtu.dk/nse-optofuidics.