NNIN/LNF Microfluidic Workshop

Total Page:16

File Type:pdf, Size:1020Kb

NNIN/LNF Microfluidic Workshop NNIN/LNF Microfluidic Workshop March 6, 2014 Pilar Herrera-Fierro Behrouz Shiari What do we mean by Micro and Nanofluidics? Microfluidics The science and engineering of systems in which fluid behavior differs from conventional flow theory primarily due to the small length scale of the system. Nanofluidics The science of building microminiaturized devices with chambers and tunnels for the containment and flow of fluids measured at the nanometer level. Low flow (nanoliter, picoliter and femtoliter) fluid delivery systems. Flow in Micro/Nanochannels • Continuum model fails when the gradients of macroscopic variables become so steep that the length scale is of the order of average distance traveled by the molecules between collision. • Knudsen number ( Kn / L ) is typical parameter used to classify the length scale and flow regimes: Kn < 0.01: Continuum approach with traditional Navier-Stokes and no-slip boundary conditions are valid. 0.01<Kn<0.1: Slip flow regime and Navier-Stokes with slip boundary conditions are applicable 0.1<Kn<10: Transition regime – Continuum approach completely breaks – Molecular Dynamic Simulation Kn > 10 : Free molecular regime – The collisionless Boltzman equation is applicable. Continuum Slip Transition Free-molecular 0 0.01 0.1 10 Kn Microfluidic Channels Microfluidics is fluid mechanics in extremely confined space. What do we gain? • Minimize physical size and space • Low power and low production cost per device • Efficient use of reagents and reactants H. G. Craighead, Science, 290, 1532 (2000). • Fast response time • Precise volumetric control • Utilize microfluidic phenomena Micro/Nanofluids Modeling Family Tree Simulation Models Molecular Models Continuum Models Deterministic Statistical Navier-Stoks MD DSMC LBM Stokes Navier–Stokes Equations, Coupled with other Equations ANSYS CFD (ANSYS Fluent and ANSYS CFX), , COMSOL, CFD- ACE+, Coventor are Structure some of commercially available software packages for simulating some of multiphysics models. Model Geometry Inlets (p=p0) Modeling Process To characterize the fluid flow’s turbulent behavior scientists generally use the Reynolds number where ρ is the fluid density, u is flow velocity, L is a characteristic length, and η is the fluid’s dynamic viscosity. Turbulent flow takes place when the Reynolds number is high, typically when Re > 2000. In the device microchannels , for water-like substances Re is less than unity. The fluid flow is thus clearly laminar, so effective mixing of fluids in the devices requires other means. Modeling: Coupled Physics We build the model with two interfaces: Laminar Flow interface solves for the fluid flow (incompressible Navier-Stokes equations) : where ρ is fluid density, u = (u, v, w) is the flow-velocity field, p is fluid pressure, I is the unit diagonal matrix, η is the fluid’s dynamic viscosity, and F = (fx, fy, fz) is a volume force affecting the fluid. F = 0 because there are no volume forces. Transport of Diluted Species interface solves for the concentration field in the fluid (convection-diffusion equation): where c is the concentration, D is the diffusion coefficient, and R=0 is the reaction rate. Modeling The selection for the Transport of Diluted Species interface ensures that no artificial accumulation of concentration arises when solving the model. Because the concentration does not affect the fluid flow, it is not necessary to solve the interfaces simultaneously. By solving them sequentially, starting with the Navier-Stokes equations, we improve the solution’s convergence and reduce the solution time. Using this technique it takes only a few minutes to solve this model. On the microchannel and mixing-chamber walls, the no-slip boundary condition applies. However, in the vertical direction, due to the geometry, we can use a symmetry boundary condition. Modeling Steps Fluid Velocity m/s Fluid Streamline Vorticity Magnitude .
Recommended publications
  • Delft University of Technology DSMC Investigation of Rarefied Gas Flow
    Delft University of Technology DSMC investigation of rarefied gas flow through diverging micro- and nanochannels Ebrahimi, Amin; Roohi, E. DOI 10.1007/s10404-017-1855-1 Publication date 2017 Document Version Accepted author manuscript Published in Microfluidics and Nanofluidics Citation (APA) Ebrahimi, A., & Roohi, E. (2017). DSMC investigation of rarefied gas flow through diverging micro- and nanochannels. Microfluidics and Nanofluidics, 21(2), [18]. https://doi.org/10.1007/s10404-017-1855-1 Important note To cite this publication, please use the final published version (if applicable). Please check the document version above. Copyright Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons. Takedown policy Please contact us and provide details if you believe this document breaches copyrights. We will remove access to the work immediately and investigate your claim. This work is downloaded from Delft University of Technology. For technical reasons the number of authors shown on this cover page is limited to a maximum of 10. This is an Accepted Author Manuscript of an article published by Springer in the journal Microfluidics and Nanofluidics, available online: http://dx.doi.org/10.1007/s10404-017-1855-1 DSMC investigation of rarefied gas flow through diverging micro- and nanochannels Amin Ebrahimi1, Ehsan Roohi2 High-Performance Computing (HPC) Laboratory, Department of Mechanical Engineering, Faculty of Engineering, Ferdowsi University of Mashhad, Mashhad, P.O. Box 91775-1111, Iran.
    [Show full text]
  • Effect of Viscous Dissipation on Slip Boundary Layer Flow of Non-Newtonian Fluid Over a Flat Plate with Convective Thermal Boundary Condition
    Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 7 (2017), pp. 3403-3432 © Research India Publications http://www.ripublication.com Effect of Viscous Dissipation on Slip Boundary Layer Flow of Non-Newtonian Fluid over a Flat Plate with Convective Thermal Boundary Condition Shashidar Reddy Borra Department of Mathematics and Humanities, Mahatma Gandhi Institute of Technology, Gandipet, R.R.District-500075, Telangana, India. Abstract The purpose of this paper is to investigate the magnetic effects of a steady, two dimensional boundary layer flow of an incompressible non-Newtonian power-law fluid over a flat plate with convective thermal and slip boundary conditions by considering the viscous dissipation. The resulting governing non-linear partial differential equations are transformed into non linear ordinary differential equations by using similarity transformation. The momentum equation is first linearized by using Quasi-linearization technique. The set of ordinary differential equations are solved numerically by using implicit finite difference scheme along with the Thomas algorithm. The solution is found to be dependent on six governing parameters including Knudsen number Knx, heat transfer parameter , magnetic field parameter M, power-law fluid index n, Eckert number Ec and Prandtl number Pr. The effects of these parameters on the velocity and temperature profiles are discussed. The special interest are the effects of the Knudsen number Knx, heat transfer parameter and Eckert number Ec on the skin friction f )0( , temperature at the wall )0( and the rate of heat transfer )0( . The numerical results are tabulated for , and )0( . Keywords: Magnetic field parameter, Knudsen number, Heat transfer parameter, Prandtl number, Viscous Dissipation and Non-Newtonian fluid.
    [Show full text]
  • Mhd Slip Flow Between Parallel Plates Heated with a Constant Heat Flux
    Isı Bilimi ve Tekniği Dergisi, 33,1, 11-20, 2013 J. of Thermal Science and Technology ©2013 TIBTD Printed in Turkey ISSN 1300-3615 MHD SLIP FLOW BETWEEN PARALLEL PLATES HEATED WITH A CONSTANT HEAT FLUX Ayşegül ÖZTÜRK Trakya University, Department of Mechanical Engineering, 22180 Edirne, Turkey, [email protected] (Geliş Tarihi: 12. 04. 2011, Kabul Tarihi: 06. 07. 2011) Abstract: The steady fully developed laminar flow and heat transfer of an electrically conducting viscous fluid between two parallel plates heated with a constant heat flux is investigated analytically in the presence of a transverse magnetic field. The momentum and energy equations are solved under the first order boundary conditions of slip velocity and temperature jump. The effects of the Knudsen number, Brinkman number and Hartmann number on the velocity and temperature distribution and heat transfer characteristics are discussed. Keywords: Parallel plate, MHD, slip-flow, viscous dissipation SABİT ISI AKISINA MARUZ PARALEL PLAKALAR ARASINDAN MHD KAYMA AKIŞI Özet: Sabit ısı akısına maruz iki paralel plaka arasından, elektrik iletkenliği olan bir viskoz akışkanın, daimi, tam gelişmiş laminer akışı ve ısı transferi akışa dik bir manyetik alanın varlığında analitik olarak incelenmiştir. Momentum ve enerji denklemleri birinci mertebe kayma hızı ve sıcaklık sıçraması sınır şartları altında çözülmüştür. Knudsen sayısı, Brinkman sayısı ve Hartmann sayısının hız, sıcaklık dağılımı ve ısı transferi karakteristikleri üzerindeki etkisi tartışılmıştır. Anahtar Kelimeler: Paralel
    [Show full text]
  • Dimensionless Numbers
    Dimensionless Numbers 3.185 November 24, 2003 Note: you are not responsible for knowing the different names of the mass transfer dimensionless numbers, just call them, e.g., “mass transfer Prandtl number”, as many people do. Those names are given here because some people use them, and you’ll probably hear them at some point in your career. Heat Transfer Mass Transfer Biot Number M.T. Biot Number Ratio of conductive to convective H.T. resistance Ratio of diffusive to reactive or conv MT resistance Determines uniformity of temperature in solid Determines uniformity of concentration in solid Bi = hL/ksolid Bi = kL/Dsolid or hD L/Dsolid Fourier Number Ratio of current time to time to reach steady­state Dimensionless time in temperature curves, used in explicit finite difference stability criterion 2 2 Fo = αt/L Fo = Dt/L Knudsen Number: Ratio of gas molecule mean free path to process lengthscale Indicates validity of line­of­sight (> 1) or continuum (< 0.01) gas models λ kT Kn = = √ L 2πσ2 P L Reynolds Number: Ratio of convective to viscous momentum transport Determines transition to turbulence, dynamic pressure vs. viscous drag Re = ρU L/µ = U L/ν Prandtl Number M.T. Prandtl (Schmidt) Number Ratio of momentum to thermal diffusivity Ratio of momentum to species diffusivity Determines ratio of fluid/HT BL thickness Determines ratio of fluid/MT BL thickness Pr = ν/α = µcp/k Pr = ν/D = µ/ρD Nusselt Number M.T. Nusselt (Sherwood) Number Ratio of lengthscale to thermal BL thickness Ratio of lengthscale to diffusion BL thickness Used to calculate heat transfer
    [Show full text]
  • Lecture 1: on Fluids, Molecules and Probabilities
    Lecture 1: On fluids, molecules and probabilities September 2, 2015 1 1 Introduction Computational fluid dynamics (CFD) is commonly associated to continuum mechanics, which surely offers an economic and powerful representation of the physics of fluids. However, the picture is deeper and broader than that; in fact modern science is increasingly concerned with flowing matter at micro and nanoscales, where the atomistic/molecular nature of the materials is fully exposed. This is why modern CFD must be informed on all three levels above. In this lecture we provide an introduction to the three main descriptions of matter: the Macroscopic (Ma) level of continuum fields, the microscopic (Mi) level of atoms and molecules and the mesoscopic (Me) level of probabil- ity distribution functions. For completeness, we also mention the quantum level, which lies underneath the atomistic ones, although this course shall not discuss moving matter below the classical atomistic scales. Which one of three levels is most apt to the solution of a given problem depends on the nature of the problem itself; in general continuum mechanics holds at macroscopic scales, say micrometers upwards, but no general rule exists; sometimes the continuum picture can be brought down to nearly nanometer scales. Conversely, under conditions of strong non-equilibrium, say large gradients in the flow configuration, the continuum picture can break down at micrometric scales, well above the atomistic ones. This is where the mesocsale description becomes essential. So the main take-home message is: there's more to fluid dynamcis than continuum mechanics. 2 The four level hierarchy The quantitative description of matter relies on two main pillars: continuum and atomistic mechanics.
    [Show full text]
  • DSMC Investigation of Rarefied Gas Flow Through Diverging Micro- and Nanochannels
    DSMC investigation of rarefied gas flow through diverging micro- and nanochannels Amin Ebrahimi1, Ehsan Roohi2 High-Performance Computing (HPC) Laboratory, Department of Mechanical Engineering, Faculty of Engineering, Ferdowsi University of Mashhad, Mashhad, P.O. Box 91775-1111, Iran. Abstract Direct simulation Monte Carlo (DSMC) method with simplified Bernoulli-trials (SBT) collision scheme has been used to study the rarefied pressure-driven nitrogen flow through diverging microchannels. The fluid behaviours flowing between two plates with different divergence angles ranging between 0o to 17o are described at different pressure ratios (1.5≤Π≤2.5) and Knudsen numbers (0.03≤Kn≤12.7). The primary flow field properties, including pressure, velocity, and temperature, are presented for divergent microchannels and are compared with those of a microchannel with a uniform cross-section. The variations of the flow field properties in divergent microchannels, which are influenced by the area change, the channel pressure ratio and the rarefication are discussed. The results show no flow separation in divergent microchannels for all the range of simulation parameters studied in the present work. It has been found that a divergent channel can carry higher amounts of mass in comparison with an equivalent straight channel geometry. A correlation between the mass flow rate through microchannels, the divergence angle, the pressure ratio, and the Knudsen number has been suggested. The present numerical findings prove the occurrence of Knudsen minimum phenomenon in micro- and Nano- channels with non-uniform cross-sections. Keywords: Divergent micro/Nanochannel; Rarefied gas flow; DSMC; simplified Bernoulli-trials; Knudsen minimum. 1 - Current affiliation: Department of Materials Science & Engineering, Delft University of Technology, Mekelweg 2, 2628 CD Delft, The Netherlands.
    [Show full text]
  • Kinetic Theory of Aerosols
    NCAR Manuscript No. 471 KINETIC THEORY OF AEROSOLS G. M. Hidy National Center for Atmospheric Research Boulder, Colorado October 1967 NATIONAL CENTER FOR ATMOSPHERIC RESEARCH Boulder, Colorado These preprint copies exclude the photographs (Plates I A, B, C). TABLE OF CONTENTS Page CAPTION LIST .......... iii I. INTRODUCTION . A. Some Phenomenological Observations . 1. Optical Effects .. ........ 2. Electrical Charging . ..... 5 3. Physicochemical Phenomena .. 5 4. Dynamical Properties . .. .... 8 5. Mechanics of Assemblies of Particles . 1 Extent . 13 6. Physical Chemistry of Aerosols and . .. their Kinetic Theory . 2. 4 II. AN IDEALIZED DYNAMIC MODEL OF AEROSOLS .. ... 13 Parameters and Characteristic Ratios. 14 A. Basic . 11 1. Mechanical Restrictions... .. 15 . 12 2. Single Particles in a Gas of Infinite Extent . 16 3. Assemblies of Particles . &0. 022. 22 . B. Single Spheres in a Continuum . .. ... 24 1. Steady Motion and Stokes Law . ... 25 2. Accelerated Motion of Particles. 27 3. Heat and Mass Transfer . 29 4. Quasi-stationary Processes and Evaporation S . 31 C. Single Spheres in Non-continuous Media . a. 33 1. The Free Molecule Regime .. ..... 33 2. The Zone of Transition in Knudsen Number .. 39 i D. Rate Processes in Particulate Clouds .. ........ 49 1. Brownian Motion .. .............. .51 2. Conservation Equations ............. 54 3. Homogeneous Nucleation . .... ... .... 61 4. Particle Interaction -- Coagulation . ...... 64 III. SOME POTENTIALLY SIGNIFICANT FEATURES OF THE PHYSICAL CHEMISTRY OF AEROSOLS . ... 71 A. Uncertainties in Dynamic Behavior . .. ..... 73 1. Accommodation and Evaporation Coefficients . 73 2. More Influences on Growth Rates by Vaporization . 75 3. Collision and Agglomeration .. .... 76 B. Size Dependence of Thermodynamic Properties . ... 77 C. Surface and Structural Heterogenieties ........ 82 IV. SUMMARY .
    [Show full text]
  • Differences Between Liquid and Gas Transport at the Microscale
    BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 53, No. 4, 2005 Differences between liquid and gas transport at the microscale M. GAD-EL-HAK¤ Virginia Commonwealth University, Richmond, Virginia 23284-3015, USA Abstract. Traditional fluid mechanics edifies the indifference between liquid and gas flows as long as certain similarity parameters – most prominently the Reynolds number – are matched. This may or may not be the case for flows in nano- or microdevices. The customary continuum, Navier-Stokes modelling is ordinarily applicable for both air and water flowing in macrodevices. Even for common fluids such as air or water, such modelling bound to fail at sufficiently small scales, but the onset for such failure is different for the two forms of matter. Moreover, when the no-slip, quasi-equilibrium Navier – Stokes system is no longer applicable, the alternative modelling schemes are different for gases and liquids. For dilute gases, statistical methods are applied and the Boltzmann equation is the cornerstone of such approaches. For liquid flows, the dense nature of the matter precludes the use of the kinetic theory of gases, and numerically intensive molecular dynamics simulations are the only alternative rooted in first principles. The present article discusses the above issues, emphasizing the differences between liquid and gas transport at the microscale and the physical phenomena unique to liquid flows in minute devices. Key words: liquid and gas transport, microscale. 1. Introduction of 45 nanometers, and approaches 10 nm in research labora- tories (source: International Technology Roadmap for Semi- Almost three centuries apart, the imaginative novelists quoted conductors, http://public.itrs.net).
    [Show full text]
  • Review of Rarefied Gas Effects in Hypersonic Applications
    Review of Rarefied Gas Effects in Hypersonic Applications Eswar Josyula and Jonathan Burt U.S. Air Force Research Laboratory Wright-Patterson Air Force Base, Ohio, USA [email protected], [email protected] ABSTRACT Rarefied gas phenomena are found in a wide variety of hypersonic flow applications, and accurate characterization of such phenomena is often desired for engineering design and analysis. In atmospheric flows involving high Mach numbers and/or low densities, underlying assumptions in a continuum fluid dynamics description tend to break down, and distributions of both molecule velocity and internal energy modes can depart significantly from equilibrium. In this paper, characteristics of rarefied and nonequilibrium gas flows are discussed, and several numerical methods for rarefied flow simulation are outlined. A variety of hypersonic flow applications are presented, and additional discussion is provided for directions of current and future research. 1.0 INTRODUCTION The field of rarefied gas dynamics is concerned with gas flow problems that do not conform to a continuum description; in particular, flows for which the Knudsen number – defined as the ratio of the mean free path to some characteristic length – is greater than that required to approximate the flow as a continuous medium. The mean free path is in turn defined as the average distance, typically in a reference frame which moves with the flow, that a molecule travels before colliding with another molecule. Traditional continuum gas dynamics models tend to break down due to gradient transport assumptions for diffusion terms in continuum equations of mass, momentum, and energy conservation. More generally, in a rarefied flow, the effect of intermolecular collisions to drive the distribution of molecular velocities toward the equilibrium limit (called a Maxwellian velocity distribution) is insufficient to maintain a near- equilibrium distribution.
    [Show full text]
  • Role of Induced Magnetic Field on MHD Natural Convection Flow In
    Alexandria Engineering Journal (2016) xxx, xxx–xxx HOSTED BY Alexandria University Alexandria Engineering Journal www.elsevier.com/locate/aej www.sciencedirect.com ORIGINAL ARTICLE Role of induced magnetic field on MHD natural convection flow in vertical microchannel formed by two electrically non-conducting infinite vertical parallel plates Basant K. Jha, Babatunde Aina * Department of Mathematics, Ahmadu Bello University, Zaria, Nigeria Received 2 March 2016; revised 25 May 2016; accepted 25 June 2016 KEYWORDS Abstract The present work consists of theoretical investigation of MHD natural convection flow Microchannel; in vertical microchannel formed by two electrically non-conducting infinite vertical parallel plates. Transverse magnetic field; The influence of induced magnetic field arising due to the motion of an electrically conducting fluid Induced magnetic field; is taken into consideration. The governing equations of the motion are a set of simultaneous ordi- Velocity slip and temperature nary differential equations and their exact solutions in dimensionless form have been obtained for jump the velocity field, the induced magnetic field and the temperature field. The expressions for the induced current density and skin friction have also been obtained. The effects of various non- dimensional parameters such as rarefaction, fluid wall interaction, Hartmann number and the mag- netic Prandtl number on the velocity, the induced magnetic field, the induced current density, and skin friction have been presented in graphical form. It is found that the effect of Hartmann number and magnetic Prandtl number on the induced current density is found to have a decreasing nature at the central region of the microchannel. Ó 2016 Faculty of Engineering, Alexandria University.
    [Show full text]
  • CONTENTS Page ABSTRACT 1
    CONTENTS Page ABSTRACT 1 NOMENCLATURE 1 INTRODUCTION 2 SUMMARY 2 Conclusions 4 Recommendations 4 APPLICATION AREAS 4 RAREFIED GAS EFFECTS 5 Calculation of Knudsen Number 5 Application to Fluid Temperature Sensor 6 MATERIAL EXPLORATIONS 13 Carbon/Carbon Composites 13 Oxidation of Unprotected Carbon/Carbon and Barrier Considerations 13 Coefficient of Thermal Expansion Considerations 14 Application of Carbon/Carbon to the Temperature Sensor 15 Ceramic Matrix Composites 17 Toughness Considerations 17 Creep Resistance 18 Fatigue 19 Oxidation 19 Strength 19 Material Candidates 20 Long Term Material 20 Materials For Short Term Use 21 Refractory Metals 23 PRESSURE TRANSDUCERS 24 Piezoelectric Transducers 24 Potential Suppliers and Products 24 Solid State Transducers 25 i Page ANALYTICAL TIME RESPONSE STUDIES 25 Time Response 25 Dissociation Effects 29 FABRICATION CAPABILITIES 34 Lockheed Missiles and Space Co. Inc. 34 Rosemount Inc., Aerospace Division 34 Schlumberger Industries, Aerospace Division 35 TSI Incorporated 35 Volna Engineering 35 CONCLUSIONS 36 Pressure Transducers 36 Materials 36 Rarefied Gas Effects 36 Fabrication Capability 37 APPENDIX A - Fluid Temperature Sensor: A Brief History A-1 APPENDIX B - References B-1 APPENDIX C - Sources of Carbon/Carbon Composites C-1 APPENDIX D - Pressure Transducer Manufacturers D-1 APPENDIX E - Bibliography E-1 ii ILLUSTRATIONS Figure Page 1 Knudsen Number vs Altitude and Mach Number 12 2 3-D Carbon/Carbon Orthogonal Weave 16 3 4-D Carbon/Carbon Hexagonal Weave 16 4 5-D Carbon/Carbon Weave 16 5 Ceramic
    [Show full text]
  • Numerical Simulations of High Knudsen Number Gas Flows and Microchannel Electrokinetic Liquid Flows Fang Yan Bakhtier Farouk, Ph
    Numerical Simulations of High Knudsen Number Gas Flows and Microchannel Electrokinetic Liquid Flows A Thesis Submitted to the Faculty of Drexel University by Fang Yan in partial fulfillment of the requirements for the degree of Doctor of Philosophy May 2003 ii ACKNOWLEDGMENTS I express my gratitude to my advisor, Prof. B. Farouk at Drexel University for his valuable guidance, support, encouragement and patience during my studies at Drexel University. I am also grateful to Prof. N. Cernansky, Prof. Y. Gogotsi, Prof. M. Vallieres, Prof. D. Wootton and Prof. K. Choi from Drexel University for their Valuable advice on this thesis as committee members. I acknowledge the financial support from Department of Mechanical Engineering and Mechanics, Drexel University that provide me the opportunity to pursue my doctoral study. I appreciate the help from the fellow laboratory members Dr. K. Bera, Dr. J. Yi, and Mr. M. Aktas. I wish to express my special gratitude to my parents and wife for their understanding and encouragement during my stay in the U. S. A. iii TABLE OF CONTENTS LIST OF TABLES......................................................................................................................... v LIST OF FIGURES...................................................................................................................... vi ABSTRACT .................................................................................................................................. x 1. INTRODUCTION ...............................................................................................................
    [Show full text]