DOCSLIB.ORG

Artificial Viscosity Model to Mitigate Numerical Artefacts at Fluid Interfaces with Surface Tension

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Artificial Viscosity Model to Mitigate Numerical Artefacts at Fluid Interfaces with Surface Tension Computers and Fluids 143 (2017) 59–72 Contents lists available at ScienceDirect Computers and Fluids journal homepage: www.elsevier.com/locate/compfluid Artificial viscosity model to mitigate numerical artefacts at fluid interfaces with surface tension ∗ Fabian Denner a, , Fabien Evrard a, Ricardo Serfaty b, Berend G.M. van Wachem a a Thermofluids Division, Department of Mechanical Engineering, Imperial College London, Exhibition Road, London, SW7 2AZ, United Kingdom b PetroBras, CENPES, Cidade Universitária, Avenida 1, Quadra 7, Sala 2118, Ilha do Fundão, Rio de Janeiro, Brazil a r t i c l e i n f o a b s t r a c t Article history: The numerical onset of parasitic and spurious artefacts in the vicinity of fluid interfaces with surface ten- Received 21 April 2016 sion is an important and well-recognised problem with respect to the accuracy and numerical stability of Revised 21 September 2016 interfacial flow simulations. Issues of particular interest are spurious capillary waves, which are spatially Accepted 14 November 2016 underresolved by the computational mesh yet impose very restrictive time-step requirements, as well as Available online 15 November 2016 parasitic currents, typically the result of a numerically unbalanced curvature evaluation. We present an Keywords: artificial viscosity model to mitigate numerical artefacts at surface-tension-dominated interfaces without Capillary waves adversely affecting the accuracy of the physical solution. The proposed methodology computes an addi- Parasitic currents tional interfacial shear stress term, including an interface viscosity, based on the local flow data and fluid Surface tension properties that reduces the impact of numerical artefacts and dissipates underresolved small scale inter- Curvature face movements. Furthermore, the presented methodology can be readily applied to model surface shear Interfacial flows viscosity, for instance to simulate the dissipative effect of surface-active substances adsorbed at the inter- face. The presented analysis of numerical test cases demonstrates the efficacy of the proposed methodol- ogy in diminishing the adverse impact of parasitic and spurious interfacial artefacts on the convergence and stability of the numerical solution algorithm as well as on the overall accuracy of the simulation results. ©2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license. ( http://creativecommons.org/licenses/by/4.0/ ) 1. Introduction currence in the modelling of surface-tension-dominated flows. Pre- vious studies [3–6] have identified two distinct origins of para- The dynamics of an interface separating two immiscible fluids sitic currents: a) a discrete imbalance between pressure “jump” is governed by the behaviour of individual molecules, which typi- and surface tension and b) a numerically unbalanced evaluation cally have a size of less than one nanometer. Computational fluid of the interface curvature. The imbalance of the pressure gradi- dynamics (CFD), however, is based on continuum mechanics and ent and surface tension can be eliminated by employing so-called treats fluids as continuous media. Molecular scales cannot be re- balanced-force discretisation methods that assure a discrete bal- solved using continuum mechanics, which leads to a variety of ance between surface tension and pressure gradient, which has theoretical and numerical difficulties with respect to the represen- previously been proposed and successfully demonstrated by Re- tation of fluid interfaces using CFD. These difficulties manifest as nardy and Renardy [14] and Francois et al. [4] for sharp surface parasitic (of unphysical origin) or spurious (of physical origin but representations such as the ghost-fluid method [15,16] , and by numerically misrepresented) flow features in the vicinity of the in- Francois et al. [4] , Mencinger and Žun [5] and Denner and van terface. Numerical artefacts of particular practical interest in inter- Wachem [6] for methods using a continuum surface force formu- facial flow modelling that have attracted significant attention and lation [1,17] . A numerically unbalanced evaluation of the interface research efforts are parasitic currents [1–7] and spurious capillary curvature, which is typically associated with spatial aliasing errors waves [1,8–13] . resulting from the computation of the second derivative of a dis- The numerical onset of unphysical flow currents in the vicin- crete indicator function or reconstruction [3,6] , leads to an unphys- ity of the interface, so-called parasitic currents , are a common oc- ical contribution to the momentum equations (via surface tension) and, consequently, causes an unphysical acceleration of the flow in the vicinity of the interface [4,7] . Using a balanced-force discretisa- ∗ Corresponding author: tion of the surface tension, parasitic currents are solely dependent E-mail address: [email protected] (F. Denner). http://dx.doi.org/10.1016/j.compfluid.2016.11.006 0045-7930/© 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license. ( http://creativecommons.org/licenses/by/4.0/ ) 60 F. Denner et al. / Computers and Fluids 143 (2017) 59–72 on the evaluation of the interface curvature and can, for instance, methodology for finite volume methods and reported that this be eliminated ( i.e. reduced to machine precision) for certain cases semi-implicit formulation of surface tension allows to exceed the by imposing the geometrically exact curvature [4–7] . Although pre- capillary time-step constraint by up to factor five without desta- vious studies [7,18] have shown that, with a careful discretisation bilising the solution of the presented two-dimensional test cases. of the interface curvature, parasitic currents can be reduced to ma- Schroeder et al. [11] proposed a two-dimensional method with a chine precision in some cases once the interface has reached a nu- semi-implicit implementation of surface tension on a Lagrangian merical equilibrium, parasitic currents are still an issue of signif- interface mesh ( i.e. using an explicit representation of the inter- icant practical relevance for applications with evolving interfaces, face) coupled to a Eulerian mesh for the flow, presenting sta- as evident by the considerable body of literature on this subject ble results for time-steps up to t = 3 t σ . In a similar fashion, published over the past five years alone ( e.g. [7,19–22] ). Zheng et al. [12] recently proposed a fully-implicit coupling of a Another important numerical artefact in interfacial flows are Lagrangian interface mesh to a MAC grid and showed that the spurious capillary waves . In a numerical framework, the shortest method can yield stable results for time-steps up to t ≈ 10 3 t σ . wavelength unambiguously resolved by the computational mesh is However, with respect to methods that rely on an implicit interface λ = , min 2 x with x representing the mesh spacing. Because an representation, such as Volume-of-Fluid methods [24] or Level-Set adequate spatial resolution of waves requires at least 6 − 10 cells methods [25,26] , and the CSF method of Brackbill et al. [1] , Denner λ ≤ per wavelength [13], capillary waves with a wavelength of min and van Wachem [13] demonstrated that the temporal resolution λ ࣠ λ 3 min are not part of the physical solution and can, therefore, requirements associated with the propagation of capillary waves be considered to be spurious capillary waves, meaning that these is a result of the spatiotemporal sampling of capillary waves and waves are a response of the discretised governing equations to a is independent of whether surface tension is implemented explicit perturbation of the interface but are numerically misrepresented. or implicit. Simulating the thermocapillary migration of a spheri- Therefore, spurious capillary waves are, contrary to parasitic cur- cal drop, Denner and van Wachem [13] demonstrated that with- rents, not the result of numerical errors but the result of the lim- out external perturbations acting at the interface (such as para- itations associated with the finite resolution of the computational sitic currents), the capillary time-step constraint can be exceeded mesh. The origin of spurious capillary waves can be physical per- by several orders of magnitude without destabilising the solution, turbations of the interface, for instance due to the collision of the presenting stable results for t = 10 4 t σ using an explicit imple- interface with an obstacle, as well as numerical perturbations, for mentation of surface tension. Thus, in order to mitigate or lift the instance the finite accuracy of numerical algorithms, discretisation capillary time-step constraint for numerical methods that rely on errors or parasitic currents [13] . An adequate temporal resolution an implicit interface representation, the shortest capillary waves of the propagation of all spatially resolved capillary waves is essen- spatially resolved by the computational mesh have to be either fil- tial for a stable numerical solution [1,13] . The dispersion relation of tered out or damped with an appropriate method to mitigate their capillary waves in inviscid fluids is given as [23] impact. Issues regarding numerical artefacts are, however, not limited σ 3 2 k ω σ = , (1) to interfacial flows. For instance, numerical oscillations induced by ρ + ρ a b a high-order discretisation of advection terms is a longstanding is- from which the phase velocity follows as c σ = ω σ /k, where ω σ is sue in CFD as well as numerical
Load more

Recommended publications
  • Capillary Surface Interfaces, Volume 46, Number 7
    fea-finn.qxp 7/1/99 9:56 AM Page 770 Capillary Surface Interfaces Robert Finn nyone who has seen or felt a raindrop, under physical conditions, and failure of unique- or who has written with a pen, observed ness under conditions for which solutions exist. a spiderweb, dined by candlelight, or in- The predicted behavior is in some cases in strik- teracted in any of myriad other ways ing variance with predictions that come from lin- with the surrounding world, has en- earizations and formal expansions, sufficiently so Acountered capillarity phenomena. Most such oc- that it led to initial doubts as to the physical va- currences are so familiar as to escape special no- lidity of the theory. In part for that reason, exper- tice; others, such as the rise of liquid in a narrow iments were devised to determine what actually oc- tube, have dramatic impact and became scientific curs. Some of the experiments required challenges. Recorded observations of liquid rise in microgravity conditions and were conducted on thin tubes can be traced at least to medieval times; NASA Space Shuttle flights and in the Russian Mir the phenomenon initially defied explanation and Space Station. In what follows we outline the his- came to be described by the Latin word capillus, tory of the problems and describe some of the meaning hair. current theory and relevant experimental results. It became clearly understood during recent cen- The original attempts to explain liquid rise in a turies that many phenomena share a unifying fea- capillary tube were based on the notion that the ture of being something that happens whenever portion of the tube above the liquid was exerting two materials are situated adjacent to each other a pull on the liquid surface.
    [Show full text]
  • Hydrostatic Shapes
    7 Hydrostatic shapes It is primarily the interplay between gravity and contact forces that shapes the macroscopic world around us. The seas, the air, planets and stars all owe their shape to gravity, and even our own bodies bear witness to the strength of gravity at the surface of our massive planet. What physics principles determine the shape of the surface of the sea? The sea is obviously horizontal at short distances, but bends below the horizon at larger distances following the planet’s curvature. The Earth as a whole is spherical and so is the sea, but that is only the first approximation. The Moon’s gravity tugs at the water in the seas and raises tides, and even the massive Earth itself is flattened by the centrifugal forces of its own rotation. Disregarding surface tension, the simple answer is that in hydrostatic equi- librium with gravity, an interface between two fluids of different densities, for example the sea and the atmosphere, must coincide with a surface of constant ¡A potential, an equipotential surface. Otherwise, if an interface crosses an equipo- ¡ A ¡ A tential surface, there will arise a tangential component of gravity which can only ¡ ©¼AAA be balanced by shear contact forces that a fluid at rest is unable to supply. An ¡ AAA ¡ g AUA iceberg rising out of the sea does not obey this principle because it is solid, not ¡ ? A fluid. But if you try to build a little local “waterberg”, it quickly subsides back into the sea again, conforming to an equipotential surface. Hydrostatic balance in a gravitational field also implies that surfaces of con- A triangular “waterberg” in the sea.
    [Show full text]
  • The Shape of a Liquid Surface in a Uniformly Rotating Cylinder in the Presence of Surface Tension
    Acta Mech DOI 10.1007/s00707-013-0813-6 Vlado A. Lubarda The shape of a liquid surface in a uniformly rotating cylinder in the presence of surface tension Received: 25 September 2012 / Revised: 24 December 2012 © Springer-Verlag Wien 2013 Abstract The free-surface shape of a liquid in a uniformly rotating cylinder in the presence of surface tension is determined before and after the onset of dewetting at the bottom of the cylinder. Two scenarios of liquid withdrawal from the bottom are considered, with and without deposition of thin film behind the liquid. The governing non-linear differential equations for the axisymmetric liquid shapes are solved numerically by an iterative procedure similar to that used to determine the equilibrium shape of a liquid drop deposited on a solid substrate. The numerical results presented are for cylinders with comparable radii to the capillary length of liquid in the gravitational or reduced gravitational fields. The capillary effects are particularly pronounced for hydrophobic surfaces, which oppose the rotation-induced lifting of the liquid and intensify dewetting at the bottom surface of the cylinder. The free-surface shape is then analyzed under zero gravity conditions. A closed-form solution is obtained in the rotation range before the onset of dewetting, while an iterative scheme is applied to determine the liquid shape after the onset of dewetting. A variety of shapes, corresponding to different contact angles and speeds of rotation, are calculated and discussed. 1 Introduction The mechanics of rotating fluids is an important part of the analysis of numerous scientific and engineering problems.
    [Show full text]
  • Energy Dissipation Due to Viscosity During Deformation of a Capillary Surface Subject to Contact Angle Hysteresis
    Physica B 435 (2014) 28–30 Contents lists available at ScienceDirect Physica B journal homepage: www.elsevier.com/locate/physb Energy dissipation due to viscosity during deformation of a capillary surface subject to contact angle hysteresis Bhagya Athukorallage n, Ram Iyer Department of Mathematics and Statistics, Texas Tech University, Lubbock, TX 79409, USA article info abstract Available online 25 October 2013 A capillary surface is the boundary between two immiscible fluids. When the two fluids are in contact Keywords: with a solid surface, there is a contact line. The physical phenomena that cause dissipation of energy Capillary surfaces during a motion of the contact line are hysteresis in the contact angle dynamics, and viscosity of the Contact angle hysteresis fluids involved. Viscous dissipation In this paper, we consider a simplified problem where a liquid and a gas are bounded between two Calculus of variations parallel plane surfaces with a capillary surface between the liquid–gas interface. The liquid–plane Two-point boundary value problem interface is considered to be non-ideal, which implies that the contact angle of the capillary surface at – Navier Stokes equation the interface is set-valued, and change in the contact angle exhibits hysteresis. We analyze a two-point boundary value problem for the fluid flow described by the Navier–Stokes and continuity equations, wherein a capillary surface with one contact angle is deformed to another with a different contact angle. The main contribution of this paper is that we show the existence of non-unique classical solutions to this problem, and numerically compute the dissipation.
    [Show full text]
  • Arxiv:1808.01401V1 [Math.NA] 4 Aug 2018 Failure in Manufacturing Microelectromechanical Systems Devices, Is Caused by an Interfacial Tension [53]
    A CONTINUATION METHOD FOR COMPUTING CONSTANT MEAN CURVATURE SURFACES WITH BOUNDARY N. D. BRUBAKER∗ Abstract. Defined mathematically as critical points of surface area subject to a volume constraint, con- stant mean curvatures (CMC) surfaces are idealizations of interfaces occurring between two immiscible fluids. Their behavior elucidates phenomena seen in many microscale systems of applied science and engineering; however, explicitly computing the shapes of CMC surfaces is often impossible, especially when the boundary of the interface is fixed and parameters, such as the volume enclosed by the surface, vary. In this work, we propose a novel method for computing discrete versions of CMC surfaces based on solving a quasilinear, elliptic partial differential equation that is derived from writing the unknown surface as a normal graph over another known CMC surface. The partial differential equation is then solved using an arc-length continuation algorithm, and the resulting algorithm produces a continuous family of CMC surfaces for varying volume whose physical stability is known. In addition to providing details of the algorithm, various test examples are presented to highlight the efficacy, accuracy and robustness of the proposed approach. Keywords. constant mean curvature, interface, capillary surface, symmetry-breaking bifurcation, arc- length continuation AMS subject classifications. 49K20, 53A05, 53A10, 76B45 1. Introduction Determining the behavior of an interface between nonmixing phases is crucial for under- standing the onset of phenomena in many microscale systems. For example, interfaces induce capillary action in tubules [23], produce beading in microfluidics [24], and change the wetting properties of patterned substrates [34]. Additionally, stiction, the leading cause of arXiv:1808.01401v1 [math.NA] 4 Aug 2018 failure in manufacturing microelectromechanical systems devices, is caused by an interfacial tension [53].
    [Show full text]
  • Statics and Dynamics of Capillary Bridges
    Colloids and Surfaces A: Physicochem. Eng. Aspects 460 (2014) 18–27 Contents lists available at ScienceDirect Colloids and Surfaces A: Physicochemical and Engineering Aspects j ournal homepage: www.elsevier.com/locate/colsurfa Statics and dynamics of capillary bridges a,∗ b Plamen V. Petkov , Boryan P. Radoev a Sofia University “St. Kliment Ohridski”, Faculty of Chemistry and Pharmacy, Department of Physical Chemistry, 1 James Bourchier Boulevard, 1164 Sofia, Bulgaria b Sofia University “St. Kliment Ohridski”, Faculty of Chemistry and Pharmacy, Department of Chemical Engineering, 1 James Bourchier Boulevard, 1164 Sofia, Bulgaria h i g h l i g h t s g r a p h i c a l a b s t r a c t • The study pertains both static and dynamic CB. • The analysis of static CB emphasis on the ‘definition domain’. • Capillary attraction velocity of CB flattening (thinning) is measured. • The thinning is governed by capillary and viscous forces. a r t i c l e i n f o a b s t r a c t Article history: The present theoretical and experimental investigations concern static and dynamic properties of cap- Received 16 January 2014 illary bridges (CB) without gravity deformations. Central to their theoretical treatment is the capillary Received in revised form 6 March 2014 bridge definition domain, i.e. the determination of the permitted limits of the bridge parameters. Concave Accepted 10 March 2014 and convex bridges exhibit significant differences in these limits. The numerical calculations, presented Available online 22 March 2014 as isogones (lines connecting points, characterizing constant contact angle) reveal some unexpected features in the behavior of the bridges.
    [Show full text]
  • Georgia Institute of Technology in Partial Fulfillment Of
    M INVESTIGATION OF CAPILLARY SPREADING POWER OF EISJLSI01IS J A THESIS Presentad to the Faculty of the Division of Graduate Student: Georgia Institute of Technology In Partial Fulfillment of the Requirements for the Degree Master of Science in Chemical Engineering fey Boiling Gay Braviley September 19/49 0752fl rt 9 AN INVESTIGATION OF CAPILLARY SPREADING POWER OP EMULSIONS Approved: ^i S'-"—y<c\-— /. ./ J Date Approved by Chairman &, \j^e^^cZt^u^ <C~*Si. ^g*i #q t AC XBQRLgDGSWTS The author wishes to express his appreciation to the personnel of Southern Sizing Company for making this study possible. To Dr. J» M« DallaValle, my advisor, I am most deeply indebted for his suggestion of the study and his valuable guidance and interest in the problem. SABLE OF CONTENTS CHAPTER PAGE I. INTRODUCTION AMD SCOPE OF THE PROBLEM 1 II. REVIEW OF THE LITERATURE 3 Capillary Flow 3 Surface Tension 5 Wotting, Spreading, and Penetration 6 III. EQUIPMT 7 Selection 7 Description and Operation 8 IV. SELECTION, STABILITY, AND PHYSICAL PROPERTIES OF EMULSIONS . , 11 V. CAPILLARY SPREADING POWER 19 Development of Equations* •••••••• • 19 Effect of Temperature •• ••••• 26 Effect of Emulsion Concentration • 26 VI. SUT«RY £SB CONCLUSIONS 29 Summary 29 Conclusions •••••••••.••••••••••••• 30 BIBLIOGRAPHY 3J APPMDIX „ 36 * LIST OF TABLES TABLES PAGE I. Physical Properties of Test Emulsion 13 II. Time of Capillary Plow vs. Viscosity-Surface Tension Ratio •••• • • 21 III. Time of Capillary Plow vs. Concentration and Temperature #••.••••»••••••••••••*• ho IV. Variation of Physical Properties of Test Emulsion with Temperature ••••• •• hi V. Variation of Physical Properties of Test Emulsion with Concentration •••••••••••• •• i\Z ' LIST OP FIGURES FIGURES PAGE 1.
    [Show full text]
  • Sharp-Interface Limit for the Navier–Stokes–Korteweg Equations
    Sharp-Interface Limit for the Navier–Stokes–Korteweg Equations Dissertation zur Erlangung des Doktorgrades vorgelegt von Johannes Daube an der Fakult¨at f¨ur Mathematik und Physik der Albert-Ludwigs-Universit¨at Freiburg Februar 2017 Dekan: Prof. Dr. Gregor Herten 1. Gutachter: Prof. Dr. Dietmar Kr¨oner 2. Gutachter: Prof. Dr. Helmut Abels Datum der m¨undlichen Pr¨ufung: 09.11.2016 Contents Abstract v Acknowledgements vii List of Symbols ix 1 Introduction 1 1.1 PhaseTransitions................................ 1 1.2 CapillaryEffects ................................. 2 1.3 Sharp- and Diffuse-Interface Models and the Sharp-Interface Limit . 3 1.4 The Navier–Stokes–Korteweg Model . .... 4 1.5 TheStaticCase.................................. 10 1.6 ExistingResults ................................. 13 1.7 NewContributions ................................ 16 1.8 Outline ...................................... 16 2 Mathematical Background 19 2.1 Notation...................................... 19 2.2 Measures ..................................... 26 2.3 Functions of Bounded Variation . 30 3 The Diffuse-Interface Model 35 3.1 The Double-Well Potential . 35 3.2 TheNotionofWeakSolutions. 37 3.3 APrioriEstimates ................................ 43 3.4 CompactnessofWeakSolutions. 56 3.5 LimitingInterfaces .............................. 60 3.6 Remarks...................................... 63 4 The Sharp-Interface Model 67 4.1 Two-Phase Incompressible Navier–Stokes Equations with Surface Tension . 67 4.2 Hypersurfaces................................... 70
    [Show full text]
  • A Capillary Surface with No Radial Limits
    A CAPILLARY SURFACE WITH NO RADIAL LIMITS A Dissertation by Colm Patric Mitchell Master of Science, Wichita State University, 2009 Bachelor of Science, Wichita State University, 2005 Submitted to the Department of Mathematics and the faculty of the Graduate School of Wichita State University in partial fulfillment of the requirements for the degree of Doctor of Philosophy May 2017 ⃝c Copyright 2017 by Colm Patric Mitchell All Rights Reserved A CAPILLARY SURFACE WITH NO RADIAL LIMITS The following faculty members have examined the final copy of this dissertation for form and content, and recommend that it be accepted in partial fulfillment of the requirement for the degree of Doctor of Philosophy with a major in Applied Mathematics. Thomas DeLillo, Committee Chair James Steck, Committee Member Elizabeth Behrman, Committee Member Ziqi Sun, Committee Member Jason Ferguson, Committee Member Accepted for the College of Liberal Arts and Sciences Ron Matson, Dean Accepted for the Graduate School Dennis Livesay, Dean iii DEDICATION To my wife Jenny, without her gentle shove I may never have gone to college; and to my boys Caleb and Connor, who had many years of dad time taken up in my studies. iv ACKNOWLEDGEMENTS I wish to acknowledge the tireless efforts of Dr. Kirk Lancaster over the last several years. Without his guidance I may never have finished my dissertation. It is unfortunate that he could not be present to chair the committee at the defense. I wish to thank Dr. Thomas DeLillo for stepping in for Kirk to chair the committee at the last minute. v ABSTRACT We begin by discussing the circumstances of capillary surfaces in regions with a corner, both concave and convex.
    [Show full text]
  • Equilibria of Pendant Droplets: Spatial Variation and Anisotropy of Surface Tension
    Equilibria of Pendant Droplets: Spatial Variation and Anisotropy of Surface Tension by Dale G. Karr Department of Naval Architecture and Marine Engineering University of Michigan Ann Arbor, MI 48109-2145 USA Email: [email protected] 1 Abstract An example of capillary phenomena commonly seen and often studied is a droplet of water hanging in air from a horizontal surface. A thin capillary surface interface between the liquid and gas develops tangential surface tension, which provides a balance of the internal and external pressures. The Young-Laplace equation has been historically used to establish the equilibrium geometry of the droplet, relating the pressure difference across the surface to the mean curvature of the surface and the surface tension, which is presumed constant and isotropic. The surface energy per unit area is often referred to as simply surface energy and is commonly considered equal to the surface tension. The relation between the surface energy and the surface tension can be established for axisymmetric droplets in a gravitational field by the application of the calculus of variations, minimizing the total potential energy. Here it is shown analytically and experimentally that, for conditions of constant volume of the droplet, equilibrium states exist with surface tensions less than the surface energy of the water-air interface. The surface tensions of the interface membrane vary with position and are anisotropic. Keywords: capillary forces, droplets, surface energy, surface tension 2 Background For droplets of water in air with negligible inertia forces, the water mass is influenced by gravity and surface tension. These forces dominate at length scales of approximately 2mm, the characteristic capillary length (1-3).
    [Show full text]
  • Viscous Fluids 9
    Viscous Fluids 9 In the previous chapter on fluids, we introduced the basic ideas of pressure, fluid flow, the application of conservation of mass and of energy in the form of the continuity equa- tion and of Bernoulli’s equation, respectively, as well as hydrostatics. Throughout those discussions we restricted ourselves to ideal fluids, those that do not exhibit any frictional properties. Often these can be neglected and the results of the previous chapter applied without any modifications whatsoever. Clearly mass is conserved even in the presence of viscous frictional forces and so the continuity equation is a very general result. Real fluids, however, do not conserve mechanical energy, but over time will lose some of this well-ordered energy to heat through frictional losses. In this chapter we consider such behavior, known as viscosity, first in the case of simple fluids such as water. We study the effects of viscosity on the motion of simple fluids and on the motion of suspended bodies, such as macromolecules, in these fluids, with special attention to flow in a cylinder, the most important geometry of flow in biology. The complex nature of blood as a fluid is studied next leading into a description and physics perspective of the human circulatory system. We conclude the chapter with a discussion of surface tension and capillarity, two important surface phenomena in fluids. In Chapter 13 we return to the general notion of the loss of well-ordered energy to heat in the context of thermodynamics. 1. VISCOSITY OF SIMPLE FLUIDS Real fluids are viscous, having internal attractive forces between the molecules so that any relative motion of molecules results in frictional, or drag, forces.
    [Show full text]
  • A Numerical Study for Liquid Bridge Based Microgripping And
    © 2007 SANTANU CHANDRA ALL RIGHTS RESERVED A NUMERICAL STUDY FOR LIQUID BRIDGE BASED MICROGRIPPING AND CONTACT ANGLE MANIPULATION BY ELECTROWETTING METHOD A Dissertation Presented to The Graduate Faculty of the University of Akron In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy Santanu Chandra December, 2007 A NUMERICAL STUDY FOR LIQUID BRIDGE BASED MICROGRIPPING AND CONTACT ANGLE MANIPULATION BY ELECTROWETTING METHOD Santanu Chandra Dissertation Approved: Accepted: _________________________________ _________________________________ Advisor Department Chair Dr. Celal Batur Dr. Celal Batur _________________________________ _________________________________ Committee Member Dean of the College of Engineering Dr. Minel Braun Dr. George K. Haritos _________________________________ _________________________________ Committee Member Dean of the Graduate School Dr. Jiang “John” Zhe Dr. George R. Newkome _________________________________ _________________________________ Committee Member Date Dr. George G. Chase _________________________________ Committee Member Dr. Gerald Young ii ABSTRACT The last decade had witnessed some very outstanding research on Micro Electro Mechanical Systems (MEMS) that has vast impact on future technologies. However building a complete microsystem requires proper microassembly methods. But microassembly research is facing stiff resistance due to the presence of many dominant forces which appears due to the scaling law. Overcoming these forces has been found to be a major drawback in microgripper research. So the primary the challenge today researches face is the lack of proper manipulation schemes. Understanding the physical forces associated with micro-scale as well as devising techniques to control them is needed in order to design a proper micromanipulation scheme. The motivation of this study is to contribute to the microassembly research by studying the forces in microscale and propose a micromanipulation scheme to control these forces.
    [Show full text]