Simulation of free surface fluids in incompressible dynamique
|
|
- Dylan Newman
- 6 years ago
- Views:
Transcription
1 Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor Interfaces in Fluid Flows 12 th December 2016 Institut Henri Poincaré 1 / 24
2 Free-surface incompressible flows 2 / 24
3 Incompressible fluids Incompressible Navier Stokes equations : ρ ( u + (u )u) µ u + p = ρf t div u = 0 u(0, x) = u 0(x). Parameters : ρ density, µ dynamic viscosity Unknowns : u velocity, p pressure Data : u 0 initial velocity, f density of external force (NS) 3 / 24
4 Incompressible fluids Incompressible Navier Stokes equations : ρ ( u + (u )u) µ u + p = ρf t div u = 0 u(0, x) = u 0(x). Parameters : ρ density, µ dynamic viscosity Unknowns : u velocity, p pressure Data : u 0 initial velocity, f density of external force Difficulties (NS) 1 Mathematical complexity of the Navier Stokes equations 2 Time evolution of the fluid domain 3 / 24
5 Incompressible fluids Incompressible Navier Stokes equations : ρ ( u + (u )u) µ u + p = ρf t div u = 0 u(0, x) = u 0(x). Parameters : ρ density, µ dynamic viscosity Unknowns : u velocity, p pressure Data : u 0 initial velocity, f density of external force Difficulties (NS) 1 Mathematical complexity of the Navier Stokes equations 2 Time evolution of the fluid domain Approaches 1 Study of shallow water approximations of incompressible Navier Stokes/Euler (Advantages : simplicity, dimension reduction) 2 Study directly incompressible Navier Stokes/Euler (Advantages : precision, general regime) 3 / 24
6 First approach : Shallow water Shallow water models : Saint-Venant (1871), Green Naghdi (1967) used mainly in Oceanography 1 The water surface is presented by an unknown function : No extra numerical difficulty to treat the free surface 2 Reduction of dimension 3 Valid in shallow water regime, sometimes for irrotational fluids 4 The water surface is presented by a spatial function (Pb : Breaking waves) 4 / 24
7 Second approach Free surface incompressibles fluid equations (Euler or Navier Stokes) : (Similar techniques as for bi-fluide Navier Stokes incompressible flows) It covers more regimes, free surface is not a function The treatment of the free surface is a numerical problem : 1 MAC methods : fixed grid, finite difference methods, free surface represented by markers transported by the flow velocity 2 ALE : Particles are transported by flow velocity (along the characteristic curves), moving mesh to deal with topology changes and sever displacements. 3 VF : Surface represented by a discontinuous characteristic function (weak formulation of advection equation), good to deal with topology changes, tedious to implement on unstructured mesh 4 LS : Surface represented by a continuous characteristic function : good to deal with topology changes, easy to implement on unstructured meshes. 5 / 24
8 Plan 1 Mathematical model 2 Numerical algorithm General idea Numerical tools Global scheme 3 Numerical results 6 / 24
9 Mathematical model 7 / 24
10 Mathematical model Fluid equations : ρ ( u + (u )u) µ u + p = ρf in Ω t t, div u = 0 in Ω t, u(0, x) = u 0(x) in Ω 0. (1) 7 / 24
11 Mathematical model Fluid equations : ρ ( u + (u )u) µ u + p = ρf in Ω t t, div u = 0 in Ω t, u(0, x) = u 0(x) in Ω 0. Boundary conditions : u n = 0 on Ω t \ Γ t, [ ( αu + µ u + t u ) n ] = 0 on Ωt \ Γt, tan ( ( µ u + t u ) p ) n = (γκ + p a ) n in Γ t, (1) (2) 7 / 24
12 Mathematical model Fluid equations : ρ ( u + (u )u) µ u + p = ρf in Ω t t, div u = 0 in Ω t, u(0, x) = u 0(x) in Ω 0. Boundary conditions : u n = 0 on Ω t \ Γ t, [ ( αu + µ u + t u ) n ] = 0 on Ωt \ Γt, tan ( ( µ u + t u ) p ) n = (γκ + p a ) n in Γ t, (1) (2) Kinematic boundary condition : { (1, u) is tangential to the free surface (t, Γ t), t 0. (3) 7 / 24
13 Level set formulation to characterize the fluid domain Ω t is characterized by a continuous function φ such that 8 / 24
14 Level set formulation to characterize the fluid domain Ω t is characterized by a continuous function φ such that The kinematic boundary condition on the free surface is now replaced by { tφ + ũ φ = 0 (x, t) D [0, T ], φ(0, x) = φ 0(x) on D, (4) where ũ is a smooth field defined on D such that ũ = u on Γ t, and φ 0 is a level set function for the initial fluid domain. 8 / 24
15 Level set formulation to characterize the fluid domain Ω t is characterized by a continuous function φ such that The kinematic boundary condition on the free surface is now replaced by { tφ + ũ φ = 0 (x, t) D [0, T ], φ(0, x) = φ 0(x) on D, (4) where ũ is a smooth field defined on D such that ũ = u on Γ t, and φ 0 is a level set function for the initial fluid domain. Question : How to construct the velocity field ũ defined on D from u defined on Ω t? 8 / 24
16 Velocity extension Question : How to construct the velocity field ũ defined on D from u defined on Ω t? a ~u + ~u = 0 in D \ Ω t ~u n = 0 on D (5) ~u = u on Ω t 9 / 24
17 Velocity extension Question : How to construct the velocity field ũ defined on D from u defined on Ω t? a ~u + ~u = 0 in D \ Ω t ~u n = 0 on D (5) ~u = u on Ω t a > 0 is large enough such that the extended velocity is smooth enough. a > 0 is small enough such that the values of the flux on the free surface do not interfere with its values outside or inside the fluid. 9 / 24
18 Velocity extension Question : How to construct the velocity field ũ defined on D from u defined on Ω t? a ~u + ~u = 0 in D \ Ω t ~u n = 0 on D (5) ~u = u on Ω t a > 0 is large enough such that the extended velocity is smooth enough. a > 0 is small enough such that the values of the flux on the free surface do not interfere with its values outside or inside the fluid. Numerical resolution : Lagrange finite elements method 9 / 24
19 Plan 1 Mathematical model 2 Numerical algorithm General idea Numerical tools Global scheme 3 Numerical results 10 / 24
20 General idea The problem is time dependent and we intend to solve it on [0, T ]. 11 / 24
21 General idea The problem is time dependent and we intend to solve it on [0, T ]. We divide the problem time interval, [0, T ] to N sub-intervals [t n, t n+1 ] n {0,..,N 1}. 11 / 24
22 General idea The problem is time dependent and we intend to solve it on [0, T ]. We divide the problem time interval, [0, T ] to N sub-intervals [t n, t n+1 ] n {0,..,N 1}. Iteration n Resolution of the problem on [t n, t n+1 ]. 11 / 24
23 General idea The problem is time dependent and we intend to solve it on [0, T ]. We divide the problem time interval, [0, T ] to N sub-intervals [t n, t n+1 ] n {0,..,N 1}. Iteration n Resolution of the problem on [t n, t n+1 ]. At each iteration : 1 We move the fluid domain following the fluid velocity using the 1 velocity extension (construction of velocity field of the level set advection equation) 2 distancing (construction of the convenient initial condition for the level set advection equation) 3 resolution of the advection equation (to find the level set corresponding to the advected domain). 2 We solve the time-discretized fluid equation on the advected domain using a convenient mesh adaptation. 11 / 24
24 Distancing The level-set advection equation on the time sub-interval [t n, t n+1 ] writes { tφ + ũ φ = 0 (x, t) D [t n, t n+1 ], (Adv) φ(t n, x) = φ t n(x) on D, where φ t n is a level-set function for the fluid domain Ω t n at time t n. 12 / 24
25 Distancing The level-set advection equation on the time sub-interval [t n, t n+1 ] writes { tφ + ũ φ = 0 (x, t) D [t n, t n+1 ], (Adv) φ(t n, x) = φ t n(x) on D, where φ t n is a level-set function for the fluid domain Ω t n at time t n. For stability reasons, φ t n must satisfy in the vicinity of the free surface. φ t n = 1, (6) 12 / 24
26 Distancing The level-set advection equation on the time sub-interval [t n, t n+1 ] writes { tφ + ũ φ = 0 (x, t) D [t n, t n+1 ], (Adv) φ(t n, x) = φ t n(x) on D, where φ t n is a level-set function for the fluid domain Ω t n at time t n. For stability reasons, φ t n must satisfy φ t n = 1, (6) in the vicinity of the free surface. Distancing is the construction of a level-set function φ t n(x) of Ω t n that (6) is satisfied [Dapogny, Frey 2012 ]. such 12 / 24
27 Distancing Problem : Loose of φ = 1 property through the advection equation Solution : Redistancing process before the resolution of the advection equation at each iteration of the algorithm : φ (x, t) + sgn(φ0)( φ 1) = 0 t φ(x, 0) = φ 0 (x) x D with φ(x) if x T Γ φ 0 (x) = + if x T h \ T Γ if x (D h \ T h ) \ T Γ where T Γ is a vicinity of the free surface. [Dapogny,Frey 2012] 13 / 24
28 Method of characteristics for advection equation Under regularity assumption for ~u and φ 0, the unique solution φ to (Adv) is constant along characteristic curves i.e. φ(x, t) = φ 0(X(x, t; 0)), (7) 14 / 24
29 Method of characteristics for advection equation Under regularity assumption for ~u and φ 0, the unique solution φ to (Adv) is constant along characteristic curves i.e. φ(x, t) = φ 0(X(x, t; 0)), (7) where s X(x, t; s) the characteristic curve of ~u passing through x at time t : dx(x, t; s) = ~u(x(x, t; s), t) for s R dt (8) X(x, t; t) = x 14 / 24
30 Method of characteristics for advection equation Under regularity assumption for ~u and φ 0, the unique solution φ to (Adv) is constant along characteristic curves i.e. φ(x, t) = φ 0(X(x, t; 0)), (7) where s X(x, t; s) the characteristic curve of ~u passing through x at time t : dx(x, t; s) = ~u(x(x, t; s), t) for s R dt (8) X(x, t; t) = x The characteristic equation is solved using a fourth order Runge-Kutta scheme. 14 / 24
31 Time discretization of the Navier Stokes equations we use the fact that ( ) u du(x(x, t; s), s) + (u )u = s=t, t ds where X(x, t; s) is the characteristic curve (parametrized by s) which cross x at time t. 15 / 24
32 Time discretization of the Navier Stokes equations we use the fact that ( ) u du(x(x, t; s), s) + (u )u = s=t, t ds where X(x, t; s) is the characteristic curve (parametrized by s) which cross x at time t. Discretized NS equation at iteration n is then ρ un (x) u n 1 X n 1 (x) µ u n (x) + p n (x) = ρf n in Ω t n t divu n (x) = 0 in Ω t n The time-discretized equation is solved by a standard finite elements method using the Uzawa algorithm. Remark : Slip boundary conditions are taken into account using a penalization technique. 15 / 24
33 Mesh adaptation The fluid free surface is explicitly discretized. The mesh is adapted such that the approximation, geometric and interpolation errors are controlled. Mesh elements outside the fluid (far from the free surface) are large since no value is of interest in this region. 16 / 24
34 Global scheme Notations : D n represents a convenient meshing of the computational domain D at time t n. T n is a subset of D n and represents a convenient meshing of the fluid domain Ω t n at time t n 17 / 24
35 Global scheme Notations : Scheme D n represents a convenient meshing of the computational domain D at time t n. T n is a subset of D n and represents a convenient meshing of the fluid domain Ω t n at time t n 1 At time t 0 = 0 : T 0 D 0 characterizing Ω t 0, u 0. 2 For n = 0,..., N 1, mesh input output Velocity extension (D n \ T n ) Γ n u n ~u n Distancing D n (D n, T n ) φ n Advection D n (~u n, φ n ) φ n+1 Mesh adaptation D n (φ n, D n ) D n+1 Interpolation D n+1 ~u n ~u n NS resolution T n+1 ~u n (u n+1, p n+1 ) 3 We return u N, p N, φ N, D N. Related paper : Pascal Frey, Dena Kazerani, Thi Thanh Mai Ta. An adaptative numerical scheme for solving incompressible two-phase and free-surface flows. Submitted. 17 / 24
36 Plan 1 Mathematical model 2 Numerical algorithm General idea Numerical tools Global scheme 3 Numerical results 18 / 24
37 Viscous fluid in the unit circular domain Parameters : ρ = 1, µ = 1, γ = 7.2, α = 10, p a = 0. Data : f = (0, 100) and the initial solution is given by : Space and time step : t 10 3, h min = h max = 0.3 and Final time : T = / 24
38 Dam break test cases The dam break test cases considered here are the same as in Marcela A. Cruchaga, Diego J. Celentano, Tayfun E. Tezduyar. Collapse of a liquid column : numerical simulation and experimental validation. Published in Comput. Mech., / 24
39 Dam break test cases The dam break test cases considered here are the same as in Marcela A. Cruchaga, Diego J. Celentano, Tayfun E. Tezduyar. Collapse of a liquid column : numerical simulation and experimental validation. Published in Comput. Mech., The test cases are done for shampoo water 20 / 24
40 Dam break test cases The dam break test cases considered here are the same as in Marcela A. Cruchaga, Diego J. Celentano, Tayfun E. Tezduyar. Collapse of a liquid column : numerical simulation and experimental validation. Published in Comput. Mech., The test cases are done for shampoo water Our results are compared with the physical experiments the numerical results presented in this reference. 20 / 24
41 Dam break test cases The dam break test cases considered here are the same as in Marcela A. Cruchaga, Diego J. Celentano, Tayfun E. Tezduyar. Collapse of a liquid column : numerical simulation and experimental validation. Published in Comput. Mech., The test cases are done for shampoo water Our results are compared with the physical experiments the numerical results presented in this reference. Differences Our framework Reference s framework Mathematical problem bi-fluid problem free-surface fluid Fluid equation computational domain fluid domain Meshing homogenous fixed meshing mobile unstructured meshing Advection equation ETILT method of characteristics Time step t = /10 times larger We use the same friction law as well as the same turbulent model (for the dam break with water problem) as in this reference. 20 / 24
42 Dam break with shampoo Parameters : ρ = 1024, µ = 8, γ = 0.07, p a = 0 and α = 10 2 (computed following the law α = ρ U C f 2 with C f = 190 and U = 0.4). Data : f = (0, 9.8). Domain dimensions : Computational domain is 0.42m 0.44m, fluid domain is 0.114m 0.114m. Time step : t = 0.01 or the first ten iterations and t = 0.02 for other iterations. Final time : T = 0.5. Space step : h min = , h max = 1 and hgrad = 2.5. t = 0.1 t = 0.2 t = 0.3 t = 0.4 t = / 24
43 Dam break with water Parameters : ρ = 1000, µ = 0.001, α = 10 6, p a = 0 and γ = Data : f = (0, 0.98). Specificity : Reynolds number is high turbulent effects energy dissipation we correct the model by modifying the viscosity as ) µ mod = min (µ + lmix 2 ρ D(u) : D(u)/2; µ max, where l mix = C t h UGN such that C t is a modelling parameter, h UGN is a characteristic element size and µ max is a cut-off value. We set here C t = 3.57, h UGN = hmin = , µ max = / 24
44 Dam break with water Parameters : ρ = 1000, µ = 0.001, α = 10 6, p a = 0 and γ = Data : f = (0, 0.98). Specificity : Reynolds number is high turbulent effects energy dissipation we correct the model by modifying the viscosity as ) µ mod = min (µ + lmix 2 ρ D(u) : D(u)/2; µ max, where l mix = C t h UGN such that C t is a modelling parameter, h UGN is a characteristic element size and µ max is a cut-off value. We set here C t = 3.57, h UGN = hmin = , µ max = 1.5. In practice : we take µ = for the first iteration of the algorithm. Then, we take the cut-off value µ max for viscosity i.e / 24
45 Dam break with water Parameters : ρ = 1000, µ = 0.001, α = 10 6, p a = 0 and γ = Data : f = (0, 0.98). Specificity : Reynolds number is high turbulent effects energy dissipation we correct the model by modifying the viscosity as ) µ mod = min (µ + lmix 2 ρ D(u) : D(u)/2; µ max, where l mix = C t h UGN such that C t is a modelling parameter, h UGN is a characteristic element size and µ max is a cut-off value. We set here C t = 3.57, h UGN = hmin = , µ max = 1.5. In practice : we take µ = for the first iteration of the algorithm. Then, we take the cut-off value µ max for viscosity i.e Time step : t = for the first iteration of the algorithm and t = for others. Final time : T= 1.5. Space step : hmax = 1 and hgrad = / 24
46 Dam break with water t = t = t = / 24
47 Thank you for your attention! 24 / 24
.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
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 informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationA Finite Element Approximation of a Coupled Model for Two Turbulent Fluids
A Finite Element Approximation of a Coupled Model for Two Turbulent Fluids Driss Yakoubi: GIREF, Université Laval, Québec. Canada joint works with Tomas Chacon Rebollo: Laboratoire Jacques-Louis Lions,
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 informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationPast, present and space-time
Past, present and space-time Arnold Reusken Chair for Numerical Mathematics RWTH Aachen Utrecht, 12.11.2015 Reusken (RWTH Aachen) Past, present and space-time Utrecht, 12.11.2015 1 / 20 Outline Past. Past
More informationThe Effect of Undulated Wall on Natural Convection Heat Transfer in Vertical Channels
Proceedings of the 3 rd International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT 16) Ottawa, Canada May 2 3, 2016 Paper No. 123 The Effect of Undulated Wall on Natural Convection Heat Transfer
More informationPhysically Based Simulations (on the GPU)
Physically Based Simulations (on the GPU) (some material from slides of Mark Harris) CS535 Fall 2014 Daniel G. Aliaga Department of Computer Science Purdue University Simulating the world Floating point
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 informationα 2 )(k x ũ(t, η) + k z W (t, η)
1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationA MOTION OF FREELY OSCILLATING DROPLET OF A YIELD STRESS FLUID: ANALYSIS AND NUMERICAL STUDIES
A MOTION OF FREELY OSCILLATING DROPLET OF A YIELD STRESS FLUID: ANALYSIS AND NUMERICAL STUDIES A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial
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 informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationNumerical Simulations of the Fluid-Structure Coupling in Physiological Vessels Mini-Workshop I
Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels INRIA Rocquencourt, projet REO miguel.fernandez@inria.fr CRM Spring School, Montréal May 16, 2005 Motivations It is generally
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 informationLevel Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
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 informationNumerical simulations and mesh adaptation : exercise sessions Charles Dapogny 1, Pascal Frey 2, Yannick Privat 2
Numerical simulations and mesh adaptation : exercise sessions Charles Dapogny Pascal Frey 2 Yannick Privat 2 The purpose of this exercise session is to give an overview of several aspects of scientific
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationApplications of parabolized stability equation for predicting transition position in boundary layers
Appl. Math. Mech. -Engl. Ed., 33(6), 679 686 (2012) DOI 10.1007/s10483-012-1579-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2012 Applied Mathematics and Mechanics (English Edition) Applications
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationEnergy dissipation caused by boundary layer instability at vanishing viscosity
Energy dissipation caused by boundary layer instability at vanishing viscosity Marie Farge, Ecole Normale Supérieure, Paris Kai Schneider, Université d Aix-Marseille in collaboration with Romain Nguyen-Nouch
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
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 informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationA Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles
A Fluctuating Immersed Boundary Method for Brownian Suspensions of Rigid Particles Aleksandar Donev Courant Institute, New York University APS DFD Meeting San Francisco, CA Nov 23rd 2014 A. Donev (CIMS)
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
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 informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
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 informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationOptimal rotary control of the cylinder wake using POD reduced order model p. 1/30
Optimal rotary control of the cylinder wake using POD reduced order model Michel Bergmann, Laurent Cordier & Jean-Pierre Brancher Michel.Bergmann@ensem.inpl-nancy.fr Laboratoire d Énergétique et de Mécanique
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationSECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES
SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES Yosuke Hasegawa Institute of Industrial Science The University of Tokyo Komaba 4-6-1, Meguro-ku, Tokyo 153-8505, Japan ysk@iis.u-tokyo.ac.jp
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
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 informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationSupplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer
Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer Thomas P. Forbes and Jason G. Kralj National Institute of Standards and Technology,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationOn the local well-posedness of compressible viscous flows with bounded density
On the local well-posedness of compressible viscous flows with bounded density Marius Paicu University of Bordeaux joint work with Raphaël Danchin and Francesco Fanelli Mathflows 2018, Porquerolles September
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationFLOWS IN LIQUID FOAMS
FLOWS IN LIQUID FOAMS A finite element approach A. SAUGEY, E. JANIAUD, W. DRENCKHAN, S. HUTZLER, D. WEAIRE Physics Department,, TRINITY COLLEGE DUBLIN Facing a problem of fluid flow in a physical system,
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More information5.1 Fluid momentum equation Hydrostatics Archimedes theorem The vorticity equation... 42
Chapter 5 Euler s equation Contents 5.1 Fluid momentum equation........................ 39 5. Hydrostatics................................ 40 5.3 Archimedes theorem........................... 41 5.4 The
More informationFinite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationM E 320 Supplementary Material Pralav Shetty
M E 320 Supplementary Material Pralav Shetty Note: In order to view the demonstrations below, you must first download CDF player to your PC/Mac/Linux. Link for CDF player http://www.wolfram.com/cdf-player/
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationNavier-Stokes equations
1 Navier-Stokes equations Introduction to spectral methods for the CSC Lunchbytes Seminar Series. Incompressible, hydrodynamic turbulence is described completely by the Navier-Stokes equations where t
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
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 informationDNS of the Taylor-Green vortex at Re=1600
DNS of the Taylor-Green vortex at Re=1600 Koen Hillewaert, Cenaero Corentin Carton de Wiart, NASA Ames koen.hillewaert@cenaero.be, corentin.carton@cenaero.be Introduction This problem is aimed at testing
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 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 informationRelativistic Hydrodynamics L3&4/SS14/ZAH
Conservation form: Remember: [ q] 0 conservative div Flux t f non-conservative 1. Euler equations: are the hydrodynamical equations describing the time-evolution of ideal fluids/plasmas, i.e., frictionless
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
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 information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
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 informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More information