Advective and Diffusive Turbulent Mixing
|
|
- Cuthbert McKinney
- 5 years ago
- Views:
Transcription
1 Advective and Diffusive Turulent Mixing J. J. Monaghan School of Mathematical Sciences Monash University, Clayton Melourne, Australia J. B. Kajtar School of Mathematical Sciences Monash University, Clayton Melourne, Australia Astract Recent studies of a turulence model for SPH [], [2] show that, for two dimensional turulence in a square container with no slip walls, it gives results that are in good agreement with those from a high accuracy spectral method. A key feature of the model is that the fluid is advected with a velocity v otained y smoothing the velocity v. The equations of motion are otained from a Lagrangian in which the kinetic energy per unit mass is 2 v v, and the movement of an SPH particle is given y dr/ = v. The linear and angular momentum and a discrete form of the circulation are conserved. The simulation of turulent flow in a square container shows that the velocity correlation function for low resolution is close to that calculated using a resolution twice as fine. The reason for the improvement is that the noise in the velocity field is greatly reduced y the smoothing. However, the smoothing of the velocity field affects the mixing ecause the velocity variations on a short length scale are removed, and there are clear differences in the flow. Some of the difference can e attriuted to the fact that the turulent flow is chaotic and the smoothing makes initial small changes to the advection and these changes are amplified. However, part of the difference is due to the effect of smoothing on the small length scale motion. In this paper we analyse the time variation of the mixing at different length scales and we show the results otained y including a diffusion equation which descries the diffusion of material due to small velocity fluctuations. This equation takes the same form as that given y Monaghan [3] in connection with the freezing of inary solutions. A key feature of the formulation is that the difference etween v and v determines the velocity fluctuation, and this determines the coefficient of diffusion. I. INTRODUCTION When a laminar flow ecomes turulent mixing occurs. In particular the linear and angular momentum are mixed and this changes the stresses in the fluid. Any scalar quantity the fluid carries is also mixed. The mixing occurs on all scales ut it varies with the length scale with the small scales eing less easily mixed than large scales. If turulence is simulated using high resolution, the turulent flow is found to eventually involve thin sheets and rions that decay y viscosity. To simulate these structures correctly it is necessary to resolve down to the Reynolds length l R which is given in 2D y l R = L/R /2 where L is the macroscopic scale and the Reynolds numer is R = V L/ν where V is the large scale velocity and ν the kinematic viscosity coefficient. These thin structures are seen in the direct SPH simulations of 2D turulence [], [2], [4]. Unfortunately, in many simulations, the resolution required to descrie the turulence correctly cannot e afforded, even with large parallel clusters. In the case of finite difference calculations the way of escape is to use a turulence model where the original fluid equations are smoothed in space. The method is called Large Eddy Simulation LES and the grid scale is larger than grid spacing of l R. The use of a coarse grid is partly compensated y including a su-grid stress that must e guessed. It has een shown that a recently proposed turulence model for SPH has a numer of desirale properties []. These include conservation of linear and angular momentum and circulation approximately together with the aility to recover the correlation functions of direct numerical simulations while using a much coarser resolution. The equations of motion are derived in a similar way to that used in a continuum Lagrangian turulence model called the Lagrangian Averaged Navier-Stokes alpha model LANS-α due to Holm and his colleagues [5], [6], and for further references see Lunasin et al. [7]. LANS was originally developed from a lengthy consideration of turulent fluctuations. However, the end result is simple. A smoothed, or regularised, velocity v is calculated y a linear operation on the un-smoothed velocity v, and then the Eulerian equations of motion are determined from a Lagrangian with kinetic energy per unit mass 2v v. The particles are moved according to dr/ = v, and the comination of v and v guarantees that circulation is conserved in the asence of dissipation. The overall result is that, in driven turulence in the asence of viscosity, the energy in small length scales is redistriuted to larger scales. Comined with a viscous term, the equations give a very satisfactory description of turulence in periodic domains [6], gyres relevant to oceanography [8], and mixing [9]. II. GOVERNING EQUATIONS The smoothed velocity v a and the un-smoothed velocity v a for an SPH liquid particle a are related y []: v a = v a + ɛ m M v v K r a r, l. Here m is the mass of particle and M is a mass closely related to the typical mass of a particle and 0 ɛ <. The function K is a smoothing function with typical length scale l. It is similar to a Gaussian. The integral of K over the space of the simulation is equal to l d where d is the numer of dimensions. The reader familiar with SPH will recognise this smoothing as the XSPH variant of SPH [0]. It can e shown that, if the velocities are expanded in a Fourier series, the
2 7 th international SPHERIC workshop Prato, Italy, May 29-3, 202 values of the coefficients for wavenumers greater than /l are reduced y a factor ɛ y the smoothing []. Typical calculations use ɛ = 0.8 in which case the coefficients of high wave numer modes are reduced y a factor 5 every time step. The Lagrangian leading to the Euler equations is the particle equivalent of Eckart s Lagrangian []. Taking into account the smoothing, it has the form L = m 2 v v uρ, s, 2 where ρ is the density of particle, u is the thermal energy/mass and s assumed constant is its entropy. Particle moves with the velocity v. By writing L in terms of v and r for each particle it is straightforward, though lengthy, to show that the Lagrange equations for particle c take the form dv c = m Pc ρ 2 c + ɛ 2 + P ρ 2 c W c m M v c v 2 c K c, 3 where P is the pressure and the equation is the SPH equivalent of the Euler equations []. Taking the limit where the numer of particles goes to infinity the SPH equation ecomes v t + v v = ρ P + ɛ ρ vr vr 2 Kr r, ldr, 4 2M where dr denotes a volume element. The last term on the right hand side is the extra stress that appears in the acceleration equation ecause of the smoothing. It is the equivalent of the Reynolds stress. If the integral involving K is expanded in a Taylor series it can e shown [] that it reduces to the form appearing in the LANS model [5]. It is straightforward to show that this equation conserves energy E in the form E = m 2 v v + uρ, s, 5 and the equations conserve linear and angular momentum when they should. A simple analysis, ased on the invariance of the system to the exchange of neighouring particles around a loop [], [3], shows that a discrete form of the circulation is conserved with an accuracy that depends on the scale of the velocity field relative to the loop over which the circulation is calculated. To model the Navier-Stokes equations an SPH viscous term is added. Because the strength of the smoothing and the new stress terms depend on the parameter ɛ we call this the SPH-ɛ turulence model. This SPH turulence model has een applied to twodimensional turulence which, though different to turulence in three dimensions, has the same features of eing non-linear and disordered []. The results of the SPH turulence model for the decaying energy and enstrophy are in good agreement with those from the experiments and other numerical simulations [], [2], [3]. In this paper we focus on the mixing that occurs in turulent flow. We will explore whether the diffusion of a material associated with different parts of the fluid in a low resolution calculation can reproduce the mixing in a high resolution calculation. We refer to this scalar quantity as the concentration associated with an element of fluid. An analogy would e to consider the concentration as a mass fraction of salt in a ody of fresh water, ut in the present case the concentration does not influence the dynamics. We egin y considering a general case, where there are more than one sustances in the fluid. The concentration of sustance j is denoted C j, where the mass of the sustance in a mass M j of fluid is CM j. The rate of change of the concentration takes a form similar to a heat conduction equation [3] such that dc j = ρ Dj C j, 6 where D j is the diffusion coefficient for sustance j. The SPH form of this equation for the rate of change of concentration of sustance j for particle a, C a j, is given y dc a j = m D a j + D j C a j C j F a, 7 ρ a ρ 2 where r a F a = a W a. This SPH form of the equation ensures that the total mass of the sustances in question is conserved. In this paper we consider only the case where there is one sustance in the fluid, and hence the superscript j can e dropped. The diffusion coefficient for particle a, D a, must have dimensions ML T. We chose a coefficient ased on the differences etween the smoothed and unsmoothed velocities of a so that D a = βρ a ṽh, 8 where β is dimensionless parameter chosen y experimentation, and we typically use β = 5. For these calculations we have used a speed ṽ = v a v a. 9 Other choices are possile, for example, a speed ased on the smoothed differences such as m ṽ = v v ρ 2 W a. 0 III. TIME STEPPING We use the same time stepping as that descried y [] ut we use a different technique to calculate the smoothed velocity field. The direct use of requires an iteration since it depends on the particle coordinates which in turn depend on v. In almost all cases that we have examined the iteration converges quickly 3 or 5 iterations. However in some cases it fails to converge. As an alternative, which is faster, and stale
3 7 th international SPHERIC workshop Prato, Italy, May 29-3, 202 at least for 0 ɛ 0.85 we convert into a differential equation. The time derivative of is d v a = f a + ɛ m M [f ak a + v a v a a K a ], where, f a = dv a / and, for any vector A, the notation A a = A a A has een used. This differential equation is integrated together with the other differential equations. IV. STUDY OF DECAYING TURBULENCE Here we consider a square container of fluid with no-slip walls. The domain is similar to that considered y Monaghan [], in which it was shown that the turulence model for SPH gives results that are good agreement with a high accuracy spectral method. The velocity field is initialised y a 4 4 set of Gaussian vortices. The vortices were initially centred on a grid of squares with separation 0.2, ut then randomly shifted in oth x and y y 0.02η, where η is a uniformly distriuted random numer in the range < η <. A given initial seed value generates a set of random numers such that the simulations are exactly reproducile. The sign of the rotation was positive or negative, and followed a checkeroard pattern. The velocity at position r due to a vortex with centre at R is vr = Ωe z r R, 2 where e z is a unit vector normal to the fluid, and Ω is given y d Ω = 2π r R 2 e r R 2 /d 2, 3 and d = The no-slip condition requires that the velocity field goes to zero at the oundary. This is achieved y smoothing the velocity field calculated from 2 []. The Reynolds numer was R = V L/ν = 000, where the characteristic velocity is V = 0., which is maximum speed at the instant the vortices are placed, and the characteristic length is L = 0.5, or half the wih of the ox. We ran the calculations at two resolutions. For the high resolution case, there were fluid particles in the domain. It has een shown [] that at this resolution, an SPH simulation using a standard algorithm reproduces the results otained y other researchers for a similar prolem. For the low resolution there were fluid particles. The turulence model plus concentration diffusion were used for the low resolution simulation. The particles in a square region initially are coloured light lue, and the outer particles are dark lue. The colouring is used simply to indicate how the advection would mix a contaminant. Fig. shows the fluid in similar states for the two resolutions at a short period of time after the vortices are applied. There is a difference in the elapsed time etween the two ecause the kinetic energy decays more rapidly for the low resolution case. Fig. a is the high resolution case without the turulence model, and Fig. is the low resolution case with the turulence model. Although there are some differences, it can e seen that the low resolution case captures many of the same features. Fig. 2 shows the low resolution case with concentration diffusion. The fluid particles in Fig. 2 and Fig. are in identical positions. The central square region initially had concentration C = dark lue, and the outer region had C = 0 yellow. As the fluid moves, the concentration diffuses and the intermediate colours are generated. The diffused concentration traces the features of the high resolution calculation more accurately than the particles alone in the low resolution calculation. As an example, consider the plume in the top left corner that curls around in an anti-clockwise direction. In Fig. a, the plume is represented y small ut steady stream of light lue particles. For the equivalent structure in Fig., it is represented y significantly few particles. However, in Fig. 2, it can e seen that the same structure is represented more accurately y the light green stream. Ultimately we would like to understand whether a lower resolution calculation with the turulence model and concentration diffusion can reproduce the same mixing properties as a direct numerical simulation. As an attempt investigate this, we divided the domain into 6 6 square cells and computed the mean concentration in each cell. For example, if a cell has ten particles of which nine particles have concentration C = 0 and one particle has C =, then the mean concentration for that cell is C = 0.. Initially, 25% of the cells have mean concentration C = the central square, and the remaining 75% have C = 0. We ran each resolution with four different sets of initial random shifts for the centres of the Gaussian vortices. This amounts to simply choosing four different initial random numer seeds. With the four sets, we then produced an ensemle average of the mean concentration level against the percentage of cells with that concentration. Fig. 3 shows the histogram as a result of this procedure. The curves are given for the high resolution case, and for the low resolution, oth with and without the consideration of the concentration diffusion. The graph shows, for example, that approximately 40% of the cells contain mean concentration 0 C < 0. in each case. If the entire domain were perfectly mixed, then 00% of the cells would contain mean concentration C = 0.25 since that was the initial distriution on the domain scale. Since all three curves are very similar, it is difficult to determine whether diffusion scheme improves the mixing measure for lower resolution. However it is interesting to note that the mixing for the lower resolution cases is as good as the high resolution case. Other measures for mixing will e explored in future, such as the more general measure considered y Roinson et al. [4]. V. CONCLUSION We have presented a scheme of turulence modelling and concentration diffusion that attempts to capture the mixing properties of a turulent flow. The turulence model for SPH is the same as that presented y Monaghan [], ut here the
4 7th international SPHERIC workshop Prato, Italy, May 29-3, 202 a Fig. 2. Concentration map for the particle case, with the turulence model and concentration diffusion. Dark lue corresponds to concentration C = and yellow corresponds to C = Percentage of cells Mean concentration Fig.. Particle positions 5s after the Gaussian vortices are applied. The particles in a central square region initially are coloured light lue, and the outer particles are dark lue. a is with fluid particles and no turulence model, is with particles and the turulence model. smoothed velocity is computed y integrating a differential equation, rather than iterating. As a result, the computation time is reduced y 25%. A concentration diffusion equation was introduced as a way of modelling the small-scale mixing features of the turulent flow. It was shown that the diffusion effectively captures the structures seen in a high resolution calculation despite the significantly smaller numer of particles in those regions with coarser resolution. The level of mixing on the scale we chose was the same for all three cases that were considered: low resolution oth with and without diffusion and high resolution. R EFERENCES [] J.J. Monaghan, A turulence model for Smoothed Particle Hydrodynamics. Eur. J. Mechanics B/Fluids. 30, , 20. Fig. 3. Histogram of mean concentration for a given cell against the percentage of cells in the domain with that concentration after t 30. Blue denotes the high resolution case, green denotes the low resolution case with diffusion, and red denotes the low resolution case without diffusion. [2] A. Valizadeh, J.J. Monaghan, Smoothed particle hydrodynamics simulations of turulence in fixed and rotating oxes in two dimensions with no-slip oundaries. Phys. Fluids, 24, 03507, 202. [3] J.J. Monaghan, Smoothed particle hydrodynamics. Rep. Prog. Phys., 68, , [4] M. Roinson, J.J. Monaghan, Direct numerical simulation of decaying two-dimensional turulence in a no-slip square ox using Smoothed Particle Hydrodynamics. Int. J. Numer. Meth. Eng., DOI: 0.002/fld.2677, 20. [5] D.D. Holm, Fluctuation effects on 3D Lagrangian mean and Eulerian mean fluid motion. Physica, 33, , 999. [6] S. Chen, D.D. Holm, L.G. Margolin, R. Zhang, Direct numerical simulations of the Navier Stokes alpha model. Physica D., 33, 66-83, 999. [7] E. Lunasin, S. Kurien, M.A. Taylor, E.S. Titi, A study of the Navier Stokes model for two dimensional turulence. J. Tur., 830, -2, [8] D.D. Holm, B.T. Nagida, Modeling mesocale turulence in the arotropic domain. J. Phys. Oceangr., 333, [9] B.J. Geurts, D.D. Holm, Leray and LANS-α modelling of turulent mixing. J. Tur., 70, -33, [0] J.J. Monaghan, On the prolem of penetration in particle methods. J. Computat. Phys., 82, -5, 989.
5 7 th international SPHERIC workshop Prato, Italy, May 29-3, 202 [] C. Eckart, Variation principles of hydrodynamics. Phys. Fluids, 3, [2] H.J.H. Clercx, S.R. Maassen, G.J.F. van Heijst, Decaying twodimensional turulence in square containers with no-slip or stress-free oundaries. Phys. of Fluids, 3, 6-626, 999. [3] G.J.F. van Heijst, H.J.H. Clercx, and D. Molenaar, The effects of solid oundaries on confined two-dimensional turulence. J. Fluid Mech., 554, 4-43, [4] M. Roinson, P. Cleary, J.J. Monaghan, Analysis of mixing in a Twin- Cam mixer using Smoothed Particle Hydrodynamics. AIChe Journal, 548, 2008.
On Decaying Two-Dimensional Turbulence in a Circular Container
Frontiers of Computational Sciences Y. Kaneda, H. Kawamura and M. Sasai (Eds.) Springer, 2007, pp. 89-95 On Decaying Two-Dimensional Turbulence in a Circular Container Kai Schneider and Marie Farge Univesité
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationProject Topic. Simulation of turbulent flow laden with finite-size particles using LBM. Leila Jahanshaloo
Project Topic Simulation of turbulent flow laden with finite-size particles using LBM Leila Jahanshaloo Project Details Turbulent flow modeling Lattice Boltzmann Method All I know about my project Solid-liquid
More informationDecaying 2D Turbulence in Bounded Domains: Influence of the Geometry
Computational Physics and New Perspectives in Turbulence Y. Kaneda (Ed.) Springer, 2007, pp. 241-246 Decaying 2D Turbulence in Bounded Domains: Influence of the Geometry Kai Schneider 1 and Marie Farge
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationValidation of an Entropy-Viscosity Model for Large Eddy Simulation
Validation of an Entropy-Viscosity Model for Large Eddy Simulation J.-L. Guermond, A. Larios and T. Thompson 1 Introduction A primary mainstay of difficulty when working with problems of very high Reynolds
More informationEnergy dissipating structures generated by dipole-wall collisions at high Reynolds number
Energy dissipating structures generated by dipole-wall collisions at high Reynolds number Duncan Sutherland 1 Charlie Macaskill 1 David Dritschel 2 1. School of Mathematics and Statistics University of
More informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationES265 Order of Magnitude Phys & Chem Convection
ES265 Order of Magnitude Phys & Chem Convection Convection deals with moving fluids in which there are spatial variations in temperature or chemical concentration. In forced convection, these variations
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationNonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto
Nonequilibrium Dynamics in Astrophysics and Material Science 2011-11-02 @ YITP, Kyoto Multi-scale coherent structures and their role in the Richardson cascade of turbulence Susumu Goto (Okayama Univ.)
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationAC & DC Magnetic Levitation and Semi-Levitation Modelling
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, March 24-26, 2003 AC & DC Magnetic Levitation and Semi-Levitation Modelling V. Bojarevics, K. Pericleous Abstract
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationUNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationModeling of turbulence in stirred vessels using large eddy simulation
Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,
More informationLecture 4: The Navier-Stokes Equations: Turbulence
Lecture 4: The Navier-Stokes Equations: Turbulence September 23, 2015 1 Goal In this Lecture, we shall present the main ideas behind the simulation of fluid turbulence. We firts discuss the case of the
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationMODELLING OF FLOW IN POROUS MEDIA AND RESIN TRANSFER MOULDING USING SMOOTHED PARTICLE HYDRODYNAMICS
Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melourne, Australia 6-8 Decemer 1999 MODELLING OF FLOW IN POROUS MEDIA AND RESIN TRANSFER MOULDING USING SMOOTHED PARTICLE
More informationCVS filtering to study turbulent mixing
CVS filtering to study turbulent mixing Marie Farge, LMD-CNRS, ENS, Paris Kai Schneider, CMI, Université de Provence, Marseille Carsten Beta, LMD-CNRS, ENS, Paris Jori Ruppert-Felsot, LMD-CNRS, ENS, Paris
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationEFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW
th International Symposium on Turulence and Shear Flo Phenomena (TSFP, Chicago, USA, July, 7 EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW Mitsuhiro Nagata Mechanical
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationOPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW
OPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW Shingo Motoki, Genta Kawahara and Masaki Shimizu Graduate School of Engineering Science Osaka University 1-3 Machikaneyama, Toyonaka, Osaka
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationStatistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders
Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Toshiki Morita 1, Takashi Watanabe 2 and Yorinobu Toya 3 1. Graduate School of Information
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationTURBULENCE IN STRATIFIED ROTATING FLUIDS Joel Sommeria, Coriolis-LEGI Grenoble
TURBULENCE IN STRATIFIED ROTATING FLUIDS Joel Sommeria, Coriolis-LEGI Grenoble Collaborations: Olivier Praud, Toulouse P.H Chavanis, Toulouse F. Bouchet, INLN Nice A. Venaille, PhD student LEGI OVERVIEW
More information4.2 Concepts of the Boundary Layer Theory
Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationEffect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
More informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationSmoothed Dissipative Particle Dynamics: theory and applications to complex fluids
2015 DPD Workshop September 21-23, 2015, Shanghai University Smoothed Dissipative Particle Dynamics: Dynamics theory and applications to complex fluids Marco Ellero Zienkiewicz Centre for Computational
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationApplication of Viscous Vortex Domains Method for Solving Flow-Structure Problems
Application of Viscous Vortex Domains Method for Solving Flow-Structure Problems Yaroslav Dynnikov 1, Galina Dynnikova 1 1 Institute of Mechanics of Lomonosov Moscow State University, Michurinskiy pr.
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationModeling of Water Flows around a Circular Cylinder with the SPH Method
Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 1 2, pp. 39 60 DOI: 10.1515/heem-2015-0003 IBW PAN, ISSN 1231 3726 Modeling of Water Flows around a Circular Cylinder with
More informationAn initiation to SPH
An initiation to SPH Lucas Braune, Thomas Lewiner Department of Mathematics, PUC Rio Rio de Janeiro, Brazil Figure 1. An SPH simulation of a dam reak. Astract The recent expansion of particle-ased methods
More informationO. A Survey of Critical Experiments
O. A Survey of Critical Experiments 1 (A) Visualizations of Turbulent Flow Figure 1: Van Dyke, Album of Fluid Motion #152. Generation of turbulence by a grid. Smoke wires show a uniform laminar stream
More informationSpectrally condensed turbulence in two dimensions
Spectrally condensed turbulence in two dimensions Hua Xia 1, Michael Shats 1, Gregory Falovich 1 The Australian National University, Canberra, Australia Weizmann Institute of Science, Rehovot, Israel Acnowledgements:
More informationThe Evolution of SPH. J. J. Monaghan Monash University Australia
The Evolution of SPH J. J. Monaghan Monash University Australia planetary disks magnetic fields cosmology radiation star formation lava phase change dam break multiphase waves bodies in water fracture
More informationChapter 1: Basic Concepts
What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information6. Basic basic equations I ( )
6. Basic basic equations I (4.2-4.4) Steady and uniform flows, streamline, streamtube One-, two-, and three-dimensional flow Laminar and turbulent flow Reynolds number System and control volume Continuity
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS
ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS Karsten Lindegård Jensen 1, B. Mutlu Sumer 1, Giovanna Vittori 2 and Paolo Blondeaux 2 The pressure field in an oscillatory boundary layer
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationVortex statistics for turbulence in a container with rigid boundaries Clercx, H.J.H.; Nielsen, A.H.
Vortex statistics for turbulence in a container with rigid boundaries Clercx, H.J.H.; Nielsen, A.H. Published in: Physical Review Letters DOI: 0.03/PhysRevLett.85.752 Published: 0/0/2000 Document Version
More informationThe Reynolds experiment
Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).
More informationComputational Fluid Dynamics Modelling of Natural Convection in Copper Electrorefining
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Abstract Computational Fluid Dynamics Modelling of Natural Convection in Copper Electrorefining A computational
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
More informationMax Planck Institut für Plasmaphysik
ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More informationLocal flow structure and Reynolds number dependence of Lagrangian statistics in DNS of homogeneous turbulence. P. K. Yeung
Local flow structure and Reynolds number dependence of Lagrangian statistics in DNS of homogeneous turbulence P. K. Yeung Georgia Tech, USA; E-mail: pk.yeung@ae.gatech.edu B.L. Sawford (Monash, Australia);
More informationCHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM
CHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM Suheni and Syamsuri Department of Mechanical Engineering, Adhi Tama Institute of Technology Surabaya, Indonesia E-Mail:
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More information15. Physics of Sediment Transport William Wilcock
15. Physics of Sediment Transport William Wilcock (based in part on lectures by Jeff Parsons) OCEAN/ESS 410 Lecture/Lab Learning Goals Know how sediments are characteried (sie and shape) Know the definitions
More information