Efficient energy based modeling and experimental validation of liquid filling in planar micro-fluidic components and networks.

This paper presents a model that describes how liquid flow fills micro-fluidic components and networks. As an alternative to computational fluid dynamic (CFD) simulations, we use a constrained energy minimization approach. This approach is based on two assumptions that hold in many micro-fluidic devices: (i) The length scales are small, and we consider slow filling rates, hence fluid momentum and viscous terms are small compared to surface tension forces, consequently the liquid/gas interfaces can be viewed as a succession of quasi-steady equilibrium configurations. (ii) Any equilibrium configuration corresponds to a surface tension energy minima which is constrained by the device shape and the volume of liquid in the device. The model is developed for planar micro-fluidic devices, is based on a fundamental physical principle, and shows accurate agreement with experimental data. It takes us only a few minutes to evaluate the model for a planar component of any shape using the Surface Evolver software, and this is significantly less then the computer run time required for CFD simulations. Moreover, once a library of component models has been created (which takes less than an hour of computer time) it then takes only seconds to simulate different network architectures with thousands of components. This fast "reconfigure the network and simulate in seconds" capability is essential for the design of truly complex networks that will enable the next generation of passive, micro-fluidic, lab-on-a-chip systems.

[1]  Marc Madou,et al.  Experimental investigation and numerical simulation of injection molding with micro-features , 2002 .

[2]  J. Troutman.  Variational Calculus with Elementary Convexity , 1983 .

[3]  J. Monaghan.  Particle methods for hydrodynamics , 1985 .

[4]  Constantinos Theodoropoulos,et al.  Order Reduction of Nonlinear Dynamic Models for Distributed Reacting Systems , 1998 .

[5]  Adhemar Bultheel,et al.  Pade´ techniques for model reduction in linear system theory: a survey , 1986 .

[6]  Hyoung J. Cho,et al.  A New Fixed-Volume Metering Microdispenser Module based on sPROMs Technology , 2001 .

[7]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[8]  Peter M. Martin,et al.  Laminated plastic microfluidic components for biological and chemical systems , 1999 .

[9]  Rui Qiao,et al.  A compact model for electroosmotic flows in microfluidic devices , 2002 .

[10]  Constantinos Theodoropoulos,et al.  Order reduction for nonlinear dynamic models of distributed reacting systems , 2000 .

[11]  P. C. Hiemenz,et al.  Principles of colloid and surface chemistry , 1977 .

[12]  R. Leighton,et al.  Feynman Lectures on Physics , 1971 .

[13]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[14]  Benjamin Shapiro,et al.  Equilibrium and dynamic behavior of micro flows under electrically induced surface tension actuation forces , 2003, Proceedings International Conference on MEMS, NANO and Smart Systems.

[15]  J. Peraire,et al.  Balanced Model Reduction via the Proper Orthogonal Decomposition , 2002 .

[16]  Wendy D. Bennett,et al.  Microfluidic Devices on Polymer Substrates for Bioanalytical Applications , 1999 .

[17]  Jin-Woo Choi,et al.  Disposable smart lab on a chip for point-of-care clinical diagnostics , 2004, Proceedings of the IEEE.

[18]  R. Feynman,et al.  The Feynman Lectures on Physics Addison-Wesley Reading , 1963 .

[19]  R. Probstein.  Physicochemical Hydrodynamics: An Introduction , 1989 .

[20]  M. M. Athavale,et al.  Computational design of membrane pumps with active/passive valves for microfluidic MEMS , 1999, Design, Test, Integration, and Packaging of MEMS/MOEMS.

[21]  Jerzy M. Floryan,et al.  Numerical Methods for Viscous Flows With Moving Boundaries , 1989 .

[22]  S. Herminghaus,et al.  Wetting: Statics and dynamics , 1997 .

[23]  Diederich Hinrichsen,et al.  Modellreduktion mit Hilfe balancierter Realisierungen / Model reduction using balanced realizations , 1990 .

[24]  Thilo Penzl.  Algorithms for model reduction of large dynamical systems , 2006 .

[25]  B. Shapiro,et al.  Equilibrium behavior of sessile drops under surface tension, applied external fields, and material variations , 2003 .

[26]  Paul W. Tuinenga,et al.  SPICE (3rd ed.): a guide to circuit simulation and analysis using PSpice , 1995 .

[27]  Jijun Xiong,et al.  Passive valves based on hydrophobic microfluidics , 2003 .

[28]  A Manz,et al.  Modular approach to fabrication of three-dimensional microchannel systems in PDMS-application to sheath flow microchips. , 2001, Lab on a chip.

[29]  Wr Graham,et al.  OPTIMAL CONTROL OF VORTEX SHEDDING USING LOW-ORDER MODELS. PART I-OPEN-LOOP MODEL DEVELOPMENT , 1999 .

[30]  Stanley Osher,et al.  Level Set Methods , 2003 .

[31]  E. B. Dussan,et al.  LIQUIDS ON SOLID SURFACES: STATIC AND DYNAMIC CONTACT LINES , 1979 .

[32]  R. A. Brown,et al.  The shape and stability of rotating liquid drops , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[33]  D. Quéré,et al.  Self-propelling slugs , 2002, Journal of Fluid Mechanics.

[34]  C Gärtner,et al.  Polymer microfabrication methods for microfluidic analytical applications , 2000, Electrophoresis.

[35]  D. J. Harrison,et al.  Electroosmotic pumping and electrophoretic separations for miniaturized chemical analysis systems , 1994 .

[36]  N F de Rooij,et al.  Multi-layer microfluidic glass chips for microanalytical applications , 2001, Fresenius' journal of analytical chemistry.

[37]  Naresh K. Sinha,et al.  Reduced-order models for linear systems , 1992, [Proceedings] 1992 IEEE International Conference on Systems, Man, and Cybernetics.

[38]  Paul W. Tuinenga,et al.  SPICE: A Guide to Circuit Simulation and Analysis Using PSpice , 1988 .

[39]  Kwang-Cheol Lee,et al.  Micro-channel filling flow considering surface tension effect , 2002 .

[40]  M.-A. Grétillat,et al.  A New Fabrication Method for Borosilicate Glass Capillary Tubes with Lateral Inlets and Outlets , 1997 .

[41]  R. Fedkiw,et al.  A numerical method for two-phase flow consisting of separate compressible and incompressible regions , 2000 .

[42]  D. Sorensen,et al.  A Survey of Model Reduction Methods for Large-Scale Systems , 2000 .

[43]  B. Shapiro,et al.  Modeling and Passive Control of Channel Filling for Micro Fluidic Networks with Thousands of Components , 2003 .

[44]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[45]  L. Lin,et al.  A water-powered micro drug delivery system , 2004, Journal of Microelectromechanical Systems.

[46]  Stephen J. Wright,et al.  Numerical Optimization (Springer Series in Operations Research and Financial Engineering) , 2000 .

[47]  P. Gennes.  Wetting: statics and dynamics , 1985 .

[48]  Kenneth A. Brakke,et al.  The Surface Evolver , 1992, Exp. Math..

[49]  A J Przekwas,et al.  Computational Modeling of Microfluid Devices with Free Surface Liquid Handling , 1998 .

[50]  R. G. Christensen,et al.  Fabrication of plastic microfluid channels by imprinting methods. , 1997, Analytical chemistry.

[51]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[52]  Naresh K. Sinha,et al.  Reduced Order Models for Complex Systems—A Critical Survey , 1990 .

[53]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[54]  Frédéric Reymond,et al.  Polymer microchips bonded by O2‐plasma activation , 2002, Electrophoresis.

[55]  B. Moore.  Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[56]  Fan-Gang Tseng,et al.  Fluid filling into micro-fabricated reservoirs , 2002 .

[57]  Eric James Grimme,et al.  Krylov Projection Methods for Model Reduction , 1997 .

[58]  James M. Hyman,et al.  Numerical methods for tracking interfaces , 1984 .