Deformation of a viscoelastic droplet passing through a microfluidic contraction

A computational fluid dynamics algorithm that uses a volume of fluid method to model immiscible phase behaviour is adapted to simulate the flow of viscoelastic fluids. Two-dimensional Cartesian simulations of a viscoelastic droplet passing through a contraction are compared against experimental results generated using a three-dimensional planar contraction. Viscoelastic effects are simulated using the Oldroyd-B rheological model. The physical and numerical models capture the deformation process well, particularly the characteristic forked tail that develops on the viscoelastic droplet while it is within the contraction. Encapsulation of continuous phase fluid is observed in the experiments, but not predicted by the simulations. We suggest this is due to inadequate mesh refinement at the length scale of a fine thread entrainment process. The numerical results show the existence of an elastic stress within the droplet that could cause the entrainment and subsequent encapsulation. The process has application potential in the manufacture of micron-sized drug and chemical delivery capsules.

[1]  Yuriko Renardy,et al.  Development and implementation of VOF-PROST for 3D viscoelastic liquid–liquid simulations , 2006 .

[2]  Darwin R. Reyes,et al.  Micro total analysis systems. 2. Analytical standard operations and applications. , 2002, Analytical chemistry.

[3]  M. F. Webster,et al.  Numerical prediction of extensional flows in contraction geometries: hybrid finite volume/element method , 2002 .

[4]  J. Rothstein,et al.  Stable jets of viscoelastic fluids and self-assembled cylindrical capsules by hydrodynamic focusing. , 2005, Langmuir : the ACS journal of surfaces and colloids.

[5]  Daniel R. Lester,et al.  Drop formation dynamics of constant low-viscosity, elastic fluids , 2002 .

[6]  David F. Fletcher,et al.  A hydrodynamic and thermodynamic simulation of droplet impacts on hot surfaces, Part II: validation and applications , 2001 .

[7]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[8]  M. Davidson,et al.  A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets , 2006 .

[9]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[10]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

[11]  J. Azaiez,et al.  Numerical simulation of viscoelastic flows through a planar contraction , 1996 .

[12]  Gareth H. McKinley,et al.  Role of the elasticity number in the entry flow of dilute polymer solutions in micro-fabricated contraction geometries , 2007 .

[13]  M. Rudman.  A parametric study of droplet deformation through a microfluidic contraction D. J. E. Harvie M. R. Davidson † J. J. Cooper-White ‡ , 2005 .

[14]  S. Pillapakkam,et al.  A level-set method for computing solutions to viscoelastic two-phase flow , 2001 .

[15]  Daniel D. Joseph,et al.  Fluid Dynamics Of Viscoelastic Liquids , 1990 .

[16]  Victor M Ugaz,et al.  Microfabricated electrophoresis systems for DNA sequencing and genotyping applications: current technology and future directions , 2004, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[17]  S. White,et al.  Review of the entry flow problem: Experimental and numerical , 1987 .

[18]  D. V. Boger.  Viscoelastic Flows Through Contractions , 1987 .

[19]  J. L. Doricio,et al.  Solving viscoelastic free surface flows of a second‐order fluid using a marker‐and‐cell approach , 2007 .

[20]  Simulations of viscoelastic droplet deformation through a microfluidic contraction , 2006 .

[21]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[22]  S. Quake,et al.  Microfluidics: Fluid physics at the nanoliter scale , 2005 .

[23]  J. Li,et al.  Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method , 2000 .

[24]  M. Davidson,et al.  Simulations of pendant drop formation of a viscoelastic liquid , 2006 .

[25]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[26]  Javad Mostaghimi,et al.  Deposition of tin droplets on a steel plate: simulations and experiments , 1998 .

[27]  M. Rudman.  A Volume-Tracking Method for Incompressible Multifluid Flows with Large Density Variations , 1998 .

[28]  Roland Keunings,et al.  An algorithm for the simulation of transient viscoelastic flows with free surfaces , 1986 .

[29]  E. Verpoorte.  Microfluidic chips for clinical and forensic analysis , 2002, Electrophoresis.

[30]  A parametric study of droplet deformation through a microfluidic contraction: Shear thinning liquids , 2007 .

[31]  J. Cooper-White,et al.  Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration , 2006 .

[32]  Stéphane Popinet,et al.  A front-tracking algorithm for accurate representation of surface tension , 1999 .

[33]  J. M. Rallison,et al.  Creeping flow of dilute polymer solutions past cylinders and spheres , 1988 .

[34]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[35]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[36]  D. V. Boger,et al.  Experimental removal of the re‐entrant corner singularity in tubular entry flows , 1994 .

[37]  M. Davidson,et al.  An analysis of parasitic current generation in Volume of Fluid simulations , 2005 .

[38]  Howard A. Stone,et al.  ENGINEERING FLOWS IN SMALL DEVICES , 2004 .

[39]  F. Pinho,et al.  Visualizations of Boger fluid flows in a 4:1 square/square contraction , 2005 .

[40]  M. F. Tomé,et al.  Die-swell, splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows , 2007 .

[41]  Roger I. Tanner,et al.  Computing transient flows with high elasticity , 2002 .

[42]  L. G. Leal,et al.  The deformation of a viscoelastic drop subjected to steady uniaxial extensional flow of a Newtonian fluid , 1999 .

[43]  Helen Song,et al.  Millisecond kinetics on a microfluidic chip using nanoliters of reagents. , 2003, Journal of the American Chemical Society.