Simulating the dynamic behavior of immiscible binary fluids in three-dimensional chemically patterned microchannels

Physical Review E - Statistical, Nonlinear, and Soft Matter Physics

Volumn 68, Issue 1 1, 2003, Pages 115021-115028

  Kuksenok, Olga ^a   Balazs, Anna C ^a  

^a University of Pittsburgh   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

CHANNEL FLOW; COMPUTER SIMULATION; INTERFACES (MATERIALS); REYNOLDS NUMBER; SOLUBILITY; THERMODYNAMICS;

MICROFLUIDIC DEVICES;

BINARY MIXTURES;

EID: 13144288089     PISSN: 1063651X     EISSN: None     Source Type: Journal    
DOI: None     Document Type: Article

References (29)

1. D.R. Reyes, D. Iosifidis, A. Manz, Anal. Chem. 74, 2623 (2002).
2. Deirdre R. Meldrum and Mark R. Hall, Science 297, 1197 (2002).
3. G.M. Whitesides and A.D. Stroock, Phys. Today 54(6), 42 (2001).
4. P.B. Umbanhowar, V. Prasad, and D.A. Weitz, Langmuir 16, 347 (2000), and references therein.
5. B. Zhao, J. Moore, and D.J. Beebe, Science 291, 1023 (2001).
6. D.E. Kataoka and S. Troian, Nature (London) 402, 794 (1999).
7. S. Puri and S.L. Frisch, J. Phys.: Condens. Matter 9, 2109 (1997).
8. S. Puri and K. Binder, Phys. Rev. A 46, R4487 (1992); Phys. Rev. E 49, 5359 (1994).
9. S. Puri and K. Binder, Phys. Rev. A 46, R4487 (1992); Phys. Rev. E 49, 5359 (1994).
10. H. Tanaka, Phys. Rev. Lett. 70, 53 (1993); 70, 2770 (1993).
11. H. Tanaka, Phys. Rev. Lett. 70, 53 (1993); 70, 2770 (1993).
12. S. Dietrich, J. Phys.: Condens. Matter 8, 9127 (1996).
13. A. Karim et al., Phys. Rev. E 57, R6273 (1998).
14. C. Bauer and S. Dietrich, Phys. Rev. E 60, 6919 (1999).
15. M.K. Chaudhury and G.M. Whitesides, Science 256, 1539 (1992).
16. O. Kuksenok, J.M. Yeomans, and A. Balazs, Langmuir 17, 7186 (2001).
17. O. Kuksenok, J.M. Yeomans, and A. Balazs, Phys. Rev. E 65, 031502 (2002).
18. A.D. Stroock et al., Science 295, 647 (2002).
19. A.J. Bray, Adv. Phys. 43, 357 (1994).
20. A.C. Balazs et al., J. Phys. Chem. B 104, 3411 (2000).
21. Other possible choices for the fluid-substrate interaction involve a short-range interaction that is introduced only through the appropriate boundary condition (for example, as in [8]), or long-range interactions due to van der Waals forces or double layer forces [20]. One can also assume a long-ranged Lennard-Jones type of interaction between the fluid and substrate elements (for example, as in Ref. [12]). For small dimensions, for example if the channel height is 1 μm or lower, even a short-ranged interaction already extends significantly into the bulk. In Eq. (2), we chose the simplest form of the interaction potential. This potential basically implies a preferential value of the order parameter on the substrate and decays away from the substrate as a function of distance. It therefore gives a reasonable phenomenological description without specifying the exact type of interaction or atomic structure of the substrate. We note that using a power law decay instead of an exponential decay in Eq. (2) does not change the results qualitatively, if one chooses parameters to maintain the same effective strength and range of interaction.
22. P.G. de Gennes, Rev. Mod. Phys. 57, 827 (1985).
23. D. Jasnow and J. Vinals, Phys. Fluids 8, 660 (1996).
24. Y. Oono and S. Puri, Phys. Rev. A 38, 434 (1988).
25. D. Gottlieb and S. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications (Arrowsmith, Bristol, 1977).
26. i = 0 at the boundaries. (The name "tau method" originates from the fact that the resulting approximation [Eq. (8)] is the exact solution to a modified problem, which differs from the original one by a small (tau) term [23]).
27. -3). In the present paper, we study the general behavior of the system; the investigation of the regime where the instability and droplet formation occur is the subject of a separate study [O. Kusenok, D. Jasnow, J. Yeomans, and A. Balazs (unpublished)].
28. These peaks are higher for higher values of C; in fact, all distortions of the original Poiseuille velocity profile increase with increases in C. (Nonetheless, we observed that the overall behavior of the system is qualitatively similar for the higher values of C studied here.)
29. T. Ohta, D. Jasnow, and K. Kawasaki, Phys. Rev. Lett. 49, 1223 (1982).