Forces on Particles and Bubbles

Total Page:16

File Type:pdf, Size:1020Kb

Forces on Particles and Bubbles Title Forces on Particles and Bubbles M. Sommerfeld Mechanische Verfahrenstechnik Zentrum für Ingenieurwissenschaften Martin-Luther-Universität Halle-Wittenberg D-06099 Halle (Saale), Germany www-mvt.iw.uni-halle.de Martin-Luther-Universität Halle-Wittenberg Content of the Lecture BBO equation and particle tracking Forces acting on particles moving in fluids Drag force Pressure, virtual mass and Basset Transverse lift forces Electrostatic force Thermophoretic and Brownian force Importance of the different forces Particle response time and Stokes number Behaviour of bubbles and forces Particle response to oscillatory flow filed Martin-Luther-Universität Halle-Wittenberg Equation of Motion 1 The equation of motion for particles in a quiescent fluid was first derived by Basset (1888), Bousinesque (1885), and Oseen (1927) BBO-equation. A rigorous derivation of the equation of motion for non-uniform Stokes flow was performed by Maxey and Riley (1983). The BBO-equation without the Faxen terms (due to curvature of the velocity field) is given by: d u 18 µ m Du d u P = F − − P − ∇ + ∇τ + F − P mP 2 mP (uF uP ) ( p ) 0.5 mF d t ρP DP ρP Dt d t Importance of d uF d uP the different t * − * ρ µ m d t d t * (u − u ) forces ??? + 9 F F P dt + F P 0 + m g ∫ * 1 2 P π ρP DP 0 (t − t ) t Accounts for d Derivative along D Substantial : : initial condition dt particle path Dt derivative 1: drag force 2: pressure term 3: added mass 4: Basset force (with initial condition) 5: gravity force Martin-Luther-Universität Halle-Wittenberg Equation of Motion 2 The calculation of particle trajectories requires the solution of several partial differential equations: particle location particle velocity particle angular velocity d x p d u d ω = u p p = p mp = ∑ Fi Ip T dt dt dt The consideration of heat and mass transfer requires the solution of two additional partial differential equations for droplet diameter and droplet temperature. Analytical solutions for the different forces acting on the particle are only available for small particle Reynolds numbers (Stokes regime). For higher particle Reynolds numbers empirical correlations are needed to calculate the different forces. In the flowing the different forces for describing the particle motion are introduced and discussed (see also Sommerfeld 2000, 2008 and 2010). Martin-Luther-Universität Halle-Wittenberg Drag Force 1 The drag force is the most important one and includes the friction and form drag or pressure drag. For highly viscous flow (low particle Reynolds number) an analytic solution can be obtained for the drag force (Stokes 1851): F = 3 π µ D (u − u ) D P F P At higher particle Reynolds numbers however, empirical correlations are needed which are expressed as a drag coefficient: F C = D D ρ v2 A 2 rel P This results in the following expression for the drag force: ρ π F = F D2 C (u − u ) u − u D 2 4 P D F P F P Martin-Luther-Universität Halle-Wittenberg Drag Force 2 Fitting of drag coefficient data for spherical solid particles: Stokes regime Standard drag correlation (Stokes 1851): (Schiller and Naumann, 1933): 24 C = 24 24 D = + 0.687 = ReP CD (1 0.15 ReP ) fD ReP ReP 2 measurements 10 Stokes regime correlation Newton regime: 101 Newton regime [-] CD ≈ 0.44 D C 100 Critical Reynolds number -1 10 -1 0 1 2 3 4 5 6 ρ D (u − u ) 10 10 10 10 10 10 10 10 F P F P ReP = Re P [-] µF Martin-Luther-Universität Halle-Wittenberg Drag Force 3 Summary of the drag coefficient obtained in different experimental studies (Crowe et al. 2012): Martin-Luther-Universität Halle-Wittenberg Drag Force 4 The average drag coefficient for non-spherical particles (for a certain stable orientation) is determined from empirical correlations fitted to experimental data (Haider and Levenspiel 1989; Thompson and Clark 1991). 1000 24 B C CD = (1+ A ReP )+ ReP D 1+ 100 φ = 0.025 ReP φ = 0.05 [ - ] φ = 0.1 D φ = 0.2 c 10 φ = 0.3 φ = 0.5 φ = 0.7 1 φ = 0.9 A = exp (2.3288 − 6.4581 φ + 2.4486 φ2 ) φ = 1.0 -1 0 1 2 3 4 5 B = 0.0964 + 0.5565 φ 10 10 10 10 10 10 10 Re P [ - ] C = exp (4.905 −13.8944 φ +18.4222 φ2 −10.2599 φ3 ) D = exp (1.4681+12.2584 φ − 20.7322 φ2 +15.8855 φ3 ) sphericity S Martin-Luther-Universität φ = V Halle-Wittenberg SPartikel Drag Force 5 RarefactionThe importance effects: of rarefaction effects may be estimated on the basis of the Knudsen number (Crowe 2006): λ π γ Ma (u − u ) Kn = = P = F P P Ma P Molecules DP 2 ReP a Mean free path of the gas molecules: µF 8 p λ = c = Particle Mol π ρ 0.499 cMol ρF F Classification of the different regimes: Continuum: KnP < 0.015 Slip flow: 0.015 < KnP < 0.15 Transitional: 0.15 < KnP < 4.5 Free molecule: KnP > 4.5 Martin-Luther-Universität Halle-Wittenberg Drag Force 6 Cunningham correction for rarefaction effects (Davies 1945): C 1 1 D = = < CD,Stokes Cu 0.55 ReP 0.25 1+ Kn 2.514 + 0.8 exp− P 0.1< Kn <1000 Kn P P 1 [ - ] 0.1 D, Stokes 0.01 / C D C 1E-3 1E-4 1E-3 0.01 0.1 1 10 100 1000 Kn [ - ] Martin-Luther-Universität Halle-Wittenberg Drag Force 7 Other effects on the drag force: Turbulence of the surrounding fluid reduces the critical Reynolds number to about 1000. Surface roughness also causes a reduction of the critical Reynolds number. The porosity of particles results in a reduction of the drag coefficient. With increasing particle concentration the drag is considerably increased (hydrodynamic interaction). Particle motion in the vicinity of walls (normal and parallel). Martin-Luther-Universität Halle-Wittenberg Pressure Force Force on the particle due to a local pressure gradient and the shear stress in the flow: mP FP = (− ∇ p + ∇ τ) ρP With the Navier-Stokes equation one finds that: D u − ∇ + ∇ τ = ρ F − p F g D t This yields the total pressure force as: ρ D u F = m F F − g The second term is p P the buoyancy force ρP D t Martin-Luther-Universität Halle-Wittenberg Virtual Mass and Basset Force 1 For higher particle Reynolds numbers the virtual (added) mass and the Basset (history) force may be expressed as: mP d FA = 0.5 CA ρF (u F − u P ) ρP d t d The second term in the (u − u ) ρ µ m t d t F P (u − u ) Basset force accounts = F F P τ + F P 0 FB 9 CB 1 2 d for the initial slip (Reeks π ρ D ∫ (t − τ) t P P 0 & McKee 1984) Coefficients given by Odar and Hamilton (1964): 0.132 Acceleration number C = 2.1− A A2 + 0.12 C 2 ReP < 60 u F − u P 0.52 AC = CB = 0.48 + d u − u (A +1)3 D F P C P d t Martin-Luther-Universität Halle-Wittenberg Virtual Mass and Basset Force 2 Recent studies of Michaelides and Roig (2011) showed that the added mass coefficient is constant (CA =1) and the Basset coefficient should be expressed in dependence of the Strouhal number: 0.82 2.5 CB = 2.0 −1.0533 [1− exp(− 0.14 ReP Sr ) ] The Strouhal number describes the behaviour of oscillatory flows. In this context the fluid time scale is the reciprocal value of the characteristic fluid oscillations. In turbulent flows this time scale corresponds to the integral time scale of turbulence. 1 T Sr = = I St τP Consequently the Strouhal number is the reciprocal value of the particle turbulent Stokes number. Martin-Luther-Universität Halle-Wittenberg Other Forces Other important forces (field or external forces) are: Gravity force: = Fg mp g Centrifugal and Coriolis force (appears only if the equation of motion is written in polar co-ordinates): 2 vP,ϕ vP,r vP,ϕ FZ = mP FC = mP rP rP Electrostatic force: Fe = qP E Thermophoretic force Magnetic force Martin-Luther-Universität Halle-Wittenberg Electrostatic Forces 1 In an electrostatic precipitator, for example, the particles are charged by an ion-bombardment created by a negative corona discharge in the vicinity of a charging wire (Löffler 1988). The charging of the particles is caused by two mechanisms: t Field charging occurs 2 qp (t) = n ⋅e = π ε0 p DP E1 due to the convective t + τq motion of the ions and is 2 relevant for particles larger qmax = n e = π ε0 p DP E1 than about 0.5 to 1 µm. D − 1 (conducting particles: p = 3, D ∞) p = 1+ 2 D + 2 (non-conducting particles p = 1.5 ÷ 2) D: relative dielectric constant of particle ε0: absolute dielectric constant in vacuum Diffusion charging is the e: elementary charge result of the thermal motion E1: field strength in the charging region of the ions and is relevant for particles with a diameter k ⋅T D D c N e2 t = π ε P + P 0 smaller than about 0.2 µm. qP (t) 4 0 ln 1 e 2 8 ε0 k T Martin-Luther-Universität Halle-Wittenberg Electrostatic Forces 2 Calculated flow field, electrostatic potential and particle trajectories in an electrostatic precipitator (Böttner and Sommerfeld 2003) Uav = 0.5 m/s VH = 60 kV Particle charge: 0.66 ⋅ qmax Channel: 300 mm x 600 mm 0.15 Separation wires: 150 mm 0.10 Y [ m ] 0.05 Particles 12 – 20 µm 0.00 0.0 0.2 0.4 0.6 0.8 1.0 0.15 X [ m ] 0.10 Particles 0.9 – 1.5 µm Y [ m ] 0.05 0.00 0.0 0.2 0.4 0.6 Martin0.8 -Luther1.0 -Universität X [ m ] Halle-Wittenberg Thermophoretic Force Thermophoretic force due to temperature gradient (LBM): κF 6 π 1.17 + 2.18 Kn κP 2 1 FT = dPart ρF νF ∇TF κ κ κ T F, P: thermal (1+ 3.42 Kn) 1+ 2 F + 4.36 Kn F conductivity of κ P fluid and particles dT/dh=1K/mm hot cold Molecules dT/dh=10 K/mm Particle DP = 100 nm Martin-Luther-Universität Halle-Wittenberg Brownian Motion 1 Brownian particle motion (LBM): 216⋅k ν T = ζ Boltz F F FB mP 2 5 ρ ρF dP ∆t ζ Random number (0 ≤ ζ ≤ 1) π P i CCun ρF Only drag force drag force and Brownian motion Martin-Luther-Universität Halle-Wittenberg
Recommended publications
  • Rectilinear Oscillations of a Sphere Immersed in a Bounded Viscous Fluid Kenneth George Mcconnell Iowa State University
    Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1963 Rectilinear oscillations of a sphere immersed in a bounded viscous fluid Kenneth George McConnell Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Applied Mechanics Commons Recommended Citation McConnell, Kenneth George, "Rectilinear oscillations of a sphere immersed in a bounded viscous fluid " (1963). Retrospective Theses and Dissertations. 2547. https://lib.dr.iastate.edu/rtd/2547 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. This dissertation has been 64—3883 microfilmed exactly as received McCONNELL, Kenneth George, 1934— RECTILINEAR OSCILLATIONS OF A SPHERE IMMERSED IN A BOUNDED VISCOUS FLUID. Iowa State University of Science and Technology Ph.D., 1963 Engineering Mechanics University Microfilms, Inc., Ann Arbor, Michigan KiSCxïLIiNJiÀK. USUlLLâXlOaS OF A ornhiKti MSieRSêv IN A BOUNDED VISCOUS FLUID by Kenneth George McConnell A Dissertation Submitted to the Graduate Faculty in Partial Fulfillment of The Requirements for the Degree of DOCTOR OF PHILOSOPHY Major Subject: Theoretical and Applied Mechanics Approved: Signature was redacted for privacy. Signature was redacted for privacy. Head of Major Department Signature was redacted for privacy. Iowa State University Of Science and Technology Ames, Iowa 1963 ii TABLE OF CONTENTS Page I. INTRODUCTION 1 A. The Phenomenon 1 B. Survey of Literature 6 II.
    [Show full text]
  • Introduction to Added Mass
    An Internet Book on Fluid Dynamics Introduction to Added Mass Whenever acceleration is imposed on a fluid flow either by acceleration of a body in the fluid or by acceleration externally imposed on the fluid, additional fluid forces will act on the surfaces in contact with the fluid. These fluid inertial forces can be of considerable importance in many practical situations. In this and the sections which follow we will review the state of knowledge of these forces and, in particular, identify the added mass matrices that can be used to characterize them. Perhaps the most fundamental view of the phenomenon of added mass is that it defines the necessary work that is needed to change the kinetic energy associated with a fluid motion. Any fluid motion such as that which occurs when a body moves through the fluid implies a certain, positive, non-zero kinetic energy associated with the fluid motion. We will confine attention to an incompressible fluid of density, ρ, in which case the total kinetic energy, T ,isgivenby ρ 2 2 2 ρ T = (u1 + u2 + u3)dv = ujujdv (Bmba1) 2 V 2 V where uj, j =1, 2, 3 are the Cartesian components of the fluid velocity and V is entire volume or domain of fluid. If the motion of the body is one of steady rectilinear translation at a velocity U through a fluid otherwise at rest then clearly the total kinetic energy is finite and constant; it must in fact be equal to the work that had to be done on the body to get it up to that velocity after starting form rest with all velocities equal to zero.
    [Show full text]
  • Experimental Study of the Fluid Drag on a Torus at Low Reynolds Number. Dharmaratne Amarakoon Louisiana State University and Agricultural & Mechanical College
    Louisiana State University LSU Digital Commons LSU Historical Dissertations and Theses Graduate School 1982 Experimental Study of the Fluid Drag on a Torus at Low Reynolds Number. Dharmaratne Amarakoon Louisiana State University and Agricultural & Mechanical College Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_disstheses Recommended Citation Amarakoon, Dharmaratne, "Experimental Study of the Fluid Drag on a Torus at Low Reynolds Number." (1982). LSU Historical Dissertations and Theses. 3743. https://digitalcommons.lsu.edu/gradschool_disstheses/3743 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Historical Dissertations and Theses by an authorized administrator of LSU Digital Commons. For more information, please contact [email protected]. INFORMATION TO USERS This reproduction was made from a copy of a document sent to us for microfilming. While the most advanced technology has been used to photograph and reproduce this document, the quality o f the reproduction is heavily dependent upon the quality of the material submitted. The following explanation of techniques is provided to help clarify markings or notations which may appear on this reproduction. t.The sign or “target” for pages apparently lacking from the document photographed is "Missing Page(s)” . If it was possible to obtain the missing page(s) or section, they are spliced into the film along with adjacent pages. This may have necessitated cutting through an image and duplicating adjacent pages to assure complete continuity. 2. When an image on the film is obliterated with a round black mark, it is an indication of either blurred copy because of movement during exposure, duplicate copy, or copyrighted materials that should not have been filmed.
    [Show full text]
  • Sedimentation of Finite-Size Spheres in Quiescent and Turbulent Environments
    Under consideration for publication in J. Fluid Mech. 1 Sedimentation of finite-size spheres in quiescent and turbulent environments Walter Fornari1y, Francesco Picano2 and Luca Brandt1 1Linn´eFlow Centre and Swedish e-Science Research Centre (SeRC), KTH Mechanics, SE-10044 Stockholm, Sweden 2Department of Industrial Engineering, University of Padova, Via Venezia 1, 35131 Padua, Italy (Received ?; revised ?; accepted ?. - To be entered by editorial office) Sedimentation of a dispersed solid phase is widely encountered in applications and envi- ronmental flows, yet little is known about the behavior of finite-size particles in homo- geneous isotropic turbulence. To fill this gap, we perform Direct Numerical Simulations of sedimentation in quiescent and turbulent environments using an Immersed Boundary Method to account for the dispersed rigid spherical particles. The solid volume fractions considered are φ = 0:5 − 1%, while the solid to fluid density ratio ρp/ρf = 1:02. The par- ticle radius is chosen to be approximately 6 Komlogorov lengthscales. The results show that the mean settling velocity is lower in an already turbulent flow than in a quiescent fluid. The reduction with respect to a single particle in quiescent fluid is about 12% and 14% for the two volume fractions investigated. The probability density function of the particle velocity is almost Gaussian in a turbulent flow, whereas it displays large positive tails in quiescent fluid. These tails are associated to the intermittent fast sedimentation of particle pairs in drafting-kissing-tumbling motions. The particle lateral dispersion is higher in a turbulent flow, whereas the vertical one is, surprisingly, of comparable mag- nitude as a consequence of the highly intermittent behavior observed in the quiescent fluid.
    [Show full text]
  • Added Mass and Aeroelastic Stability of a Flexible Plate Interacting With
    Added Mass and Aeroelastic Stability of a Flexible Plate Rajeev K. Jaiman1 Interacting With Mean Flow Assistant Professor Department of Mechanical Engineering, in a Confined Channel National University of Singapore, 117576 Singapore e-mail: [email protected] This work presents a review and theoretical study of the added-mass and aeroelastic instability exhibited by a linear elastic plate immersed in a mean flow. We first present a Manoj K. Parmar combined added-mass result for the model problem with a mean incompressible and com- pressible flow interacting with an elastic plate. Using the Euler–Bernoulli model for the Research Assistant plate and a 2D viscous potential flow model, a generalized closed-form expression of Scientist University of Florida, added-mass force has been derived for a flexible plate oscillating in fluid. A new com- Gainesville, FL 32611 pressibility correction factor is introduced in the incompressible added-mass force to account for the compressibility effects. We present a formulation for predicting the criti- Pardha S. Gurugubelli cal velocity for the onset of flapping instability. Our proposed new formulation considers Graduate Research Assistant tension effects explicitly due to viscous shear stress along the fluid-structure interface. In National University of Singapore, general, the tension effects are stabilizing in nature and become critical in problems 117576 Singapore involving low mass ratios. We further study the effects of the mass ratio and channel height on the aeroelastic instability using the linear stability analysis. It is observed that the proximity of the wall parallel to the plate affects the growth rate of the instability, however, these effects are less significant in comparison to the mass ratio or the tension effects in defining the instability.
    [Show full text]
  • Dynamic Tidal Power (Dtp): a Review of a Promising Technique for Harvesting Sustainable Energy at Sea
    DYNAMIC TIDAL POWER (DTP): A REVIEW OF A PROMISING TECHNIQUE FOR HARVESTING SUSTAINABLE ENERGY AT SEA From : Harmen Talstra, Tom Pak (Svašek Hydraulics) To : ir. W.L. Walraven (Stichting DTP Netherlands) Date : 29 July 2020 Reference : 2037/U20232/A/HTAL Checked by : A.J. Bliek Status : Draft on behalf of whitepaper 1 INTRODUCTION This memorandum describes the technique of Dynamic Tidal Power (DTP), a conceptually new way of harvesting large-scale tidal energy at open sea; in particular we focus upon the hydrodynamic aspects of it. At present, most existing installations exploiting tidal energy encompass a structure at the mouth of a river, estuary or tidal basin. This rather small scale limits the amount of sustainable energy that can be extracted, whereas these type of exploitations may cause conflicts with other (e.g. economical or ecological) functions of vulnerable coastal or estuarine waters. The concept of DTP includes a truly large-scale sustainable energy production by utilizing the tidal wave propagation at open sea. This can be done by creating a water head difference over a (very) long dike, roughly perpendicular to the local tidal flow direction, taking advantage of the oscillatory dynamic behaviour of tidal waves. These dikes can be considered as long sequences of pre-fab solid dike modules, containing a large concentration of energy turbines. Dikes can be either attached to an existing coast line, or be constructed at a detached “stand-alone” location at open sea (for instance in combination with an offshore wind farm). The presence of such long dikes at open sea can possibly be utilized for additional economical and environmental functionalities as well.
    [Show full text]
  • 2.016 Hydrodynamics Fluid Forces on Bodies
    2.016 Hydrodynamics Reading #5 2.016 Hydrodynamics Prof. A.H. Techet Fluid Forces on Bodies 1. Steady Flow In order to design offshore structures, surface vessels and underwater vehicles, an understanding of the basic fluid forces acting on a body is needed. In the case of steady viscous flow, these forces are straightforward. Lift force, perpendicular to the velocity, and Drag force, inline with the flow, can be calculated based on the fluid velocity, U , force coefficients, CD and C L , the object’s dimensions or area, A , and fluid density, ρ . For viscous flows the drag and lift on a body are defined as follows 1 F = ρU2 AC (5.1) Drag 2 D 1 F = ρU2 AC (5.2) Lift 2 L These equations can also be used in a quiescent (stationary) fluid for a steady translating body, where U is the body velocity instead of the fluid velocity, since U is still the relative velocity of the fluid with respect to the body. The drag force arises due to viscous rubbing of the fluid. The fluid may be thought of as comprised of several “layers” which move relative to one another. The layer at the surface of the body “sticks” to the surface due to the no-slip condition. The next layer of fluid away from the surface rubs against the layer below, and this rubbing requires a certain amount of force because of viscosity. One would expect that in the absence of viscosity, the force would go to zero. Jean Le Rond d'Alembert (1717-1783) performed a series of experiments to measure the drag on a sphere in a flowing fluid, and on the basis of the potential flow analysis he expected that the force would approach zero as the viscosity of the fluid approached zero.
    [Show full text]
  • (English Edition) a Review on the Flow Instability of Nanofluids
    Appl. Math. Mech. -Engl. Ed., 40(9), 1227–1238 (2019) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-019-2521-9 A review on the flow instability of nanofluids∗ Jianzhong LIN†, Hailin YANG Department of Mechanics, State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China (Received Mar. 22, 2019 / Revised May 6, 2019) Abstract Nanofluid flow occurs in extensive applications, and hence has received widespread attention. The transition of nanofluids from laminar to turbulent flow is an important issue because of the differences in pressure drop and heat transfer between laminar and turbulent flow. Nanofluids will become unstable when they depart from the thermal equilibrium or dynamic equilibrium state. This paper conducts a brief review of research on the flow instability of nanofluids, including hydrodynamic instability and thermal instability. Some open questions on the subject are also identified. Key words nanofluid, thermal instability, hydrodynamic instability, review Chinese Library Classification O358, O359 2010 Mathematics Subject Classification 76E09, 76T20 1 Introduction Nanofluids have aroused significant interest over the past few decades for their wide applica- tions in energy, machinery, transportation, and healthcare. For example, low concentration of particles causes viscosity changes[1–2], decreases viscosity with increasing shear rate[3], changes friction factors and pressure drops[4–6], and improves heat transfer[7–9]. Moreover, the unique flow properties of nanofluids are determined by the flow pattern. Heat transfer[10] and pressure drop[11] are also much lower and higher, respectively, in laminar flow than in turbulent flow. By adding particles to the fluid, the thermal entropy generation and friction are of the same order of magnitude as in turbulent flow, while the effect of heat transfer entropy generation strongly outweighs that of the friction entropy generation in laminar flow.
    [Show full text]
  • Suspensions of Finite-Size Rigid Particles in Laminar and Turbulent
    Suspensions of finite-size rigid particles in laminar and turbulentflows by Walter Fornari November 2017 Technical Reports Royal Institute of Technology Department of Mechanics SE-100 44 Stockholm, Sweden Akademisk avhandling som med tillst˚andav Kungliga Tekniska H¨ogskolan i Stockholm framl¨agges till offentlig granskning f¨or avl¨aggande av teknologie doctorsexamenfredagen den 15 December 2017 kl 10:15 i sal D3, Kungliga Tekniska H¨ogskolan, Lindstedtsv¨agen 5, Stockholm. TRITA-MEK Technical report 2017:14 ISSN 0348-467X ISRN KTH/MEK/TR-17/14-SE ISBN 978-91-7729-607-2 Cover: Suspension offinite-size rigid spheres in homogeneous isotropic turbulence. c Walter Fornari 2017 � Universitetsservice US–AB, Stockholm 2017 “Considerate la vostra semenza: fatti non foste a viver come bruti, ma per seguir virtute e canoscenza.” Dante Alighieri, Divina Commedia, Inferno, Canto XXVI Suspensions offinite-size rigid particles in laminar and tur- bulentflows Walter Fornari Linn´eFLOW Centre, KTH Mechanics, Royal Institute of Technology SE-100 44 Stockholm, Sweden Abstract Dispersed multiphaseflows occur in many biological, engineering and geophysical applications such asfluidized beds, soot particle dispersion and pyroclastic flows. Understanding the behavior of suspensions is a very difficult task. Indeed particles may differ in size, shape, density and stiffness, their concentration varies from one case to another, and the carrierfluid may be quiescent or turbulent. When turbulentflows are considered, the problem is further complicated by the interactions between particles and eddies of different size, ranging from the smallest dissipative scales up to the largest integral scales. Most of the investigations on this topic have dealt with heavy small particles (typically smaller than the dissipative scale) and in the dilute regime.
    [Show full text]
  • Advection Diffusion Model for Particles Deposition in Rayleigh-Bénard Turbulent Flows
    International Conference Nuclear Energy for New Europe 2005 Bled, Slovenia, September 5-8, 2005 Advection Diffusion Model for Particles Deposition in Rayleigh-Bénard Turbulent Flows Paolo Oresta, Antonio Lippolis DIASS – Politecnico di Bari Viale del Turismo 8, 74100 Taranto, Italia [email protected], [email protected] Roberto Verzicco DIMeG and CEMeC – Politecnico di Bari Via Re David, 200, 70125 Bari, Italia [email protected] Alfredo Soldati DEM and CIFI – Università degli Studi di Udine Via delle Scienze, 208, 33100 Udine, Italia [email protected] ABSTRACT In this paper, Direct Numerical Simulation (DNS) and Lagrangian Particle Tracking are used to precisely investigate the turbulent thermally driven flow and particles dispersion in a closed, slender cylindrical domain. The numerical simulations are carried out for Rayleigh (Ra) and Prandtl numbers (Pr) equal to Ra = 2·108 and Pr = 0.7, considering three sets of particles with Stokes numbers, based on Kolmogorov scale, equal to Stk = 1.3, Stk = 0.65 and Stk = 0.13. This data are used to calculate a priori the drift velocity and the turbulent diffusion coefficient for the Advection Diffusion model. These quantities are function of the Stokes, Froude, Rayleigh and Prandtl numbers only. One dimensional, time dependent, Advection- Diffusion Equation (ADE) is presented to predict particles deposition in Rayleigh-Bénard flow in the cylindrical domain. This archetype configuration models flow and aerosol dynamics, produced in case of accident in the passive containment cooling system (PCCS) of a nuclear reactor. ADE results show a good agreement with DNS data for all the sets of particles investigated. 1 INTRODUCTION The passive containment cooling system (PCCS) consists of an inner steel shell and an outer concrete shell [1].
    [Show full text]
  • On Dimensionless Numbers
    chemical engineering research and design 8 6 (2008) 835–868 Contents lists available at ScienceDirect Chemical Engineering Research and Design journal homepage: www.elsevier.com/locate/cherd Review On dimensionless numbers M.C. Ruzicka ∗ Department of Multiphase Reactors, Institute of Chemical Process Fundamentals, Czech Academy of Sciences, Rozvojova 135, 16502 Prague, Czech Republic This contribution is dedicated to Kamil Admiral´ Wichterle, a professor of chemical engineering, who admitted to feel a bit lost in the jungle of the dimensionless numbers, in our seminar at “Za Plıhalovic´ ohradou” abstract The goal is to provide a little review on dimensionless numbers, commonly encountered in chemical engineering. Both their sources are considered: dimensional analysis and scaling of governing equations with boundary con- ditions. The numbers produced by scaling of equation are presented for transport of momentum, heat and mass. Momentum transport is considered in both single-phase and multi-phase flows. The numbers obtained are assigned the physical meaning, and their mutual relations are highlighted. Certain drawbacks of building correlations based on dimensionless numbers are pointed out. © 2008 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved. Keywords: Dimensionless numbers; Dimensional analysis; Scaling of equations; Scaling of boundary conditions; Single-phase flow; Multi-phase flow; Correlations Contents 1. Introduction .................................................................................................................
    [Show full text]
  • Dynamics of Heavy and Buoyant Underwater Pendulums
    This draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics 1 Dynamics of heavy and buoyant underwater pendulums Varghese Mathai1,2y, Laura A. W. M. Loeffen2y, Timothy T. K. Chan2,3, and Sander Wildeman4,2z 1School of Engineering, Brown University, Providence, RI 02912, USA. 2Physics of Fluids Group and Max Planck Center for Complex Fluids, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands. 3Department of Physics, The Chinese University of Hong Kong, Shatin, Hong Kong. 4Institut Langevin, ESPCI, CNRS, PSL Research University, 1 rue Jussieu, 75005 Paris, France. (Received xx; revised xx; accepted xx) The humble pendulum is often invoked as the archetype of a simple, gravity driven, oscillator. Under ideal circumstances, the oscillation frequency of the pendulum is inde- pendent of its mass and swing amplitude. However, in most real-world situations, the dynamics of pendulums is not quite so simple, particularly with additional interactions between the pendulum and a surrounding fluid. Here we extend the realm of pendulum studies to include large amplitude oscillations of heavy and buoyant pendulums in a fluid. We performed experiments with massive and hollow cylindrical pendulums in water, and constructed a simple model that takes the buoyancy, added mass, fluid (nonlinear) drag, and bearing friction into account. To first order, the model predicts the oscillation frequencies, peak decelerations and damping rate well. An interesting effect of the nonlinear drag captured well by the model is that for heavy pendulums, the damping time shows a non-monotonic dependence on pendulum mass, reaching a minimum when the pendulum mass density is nearly twice that of the fluid.
    [Show full text]