Direct numerical simulations of breaking waves and adaptive dispersive wave models
|
|
- Justin Cameron Ramsey
- 5 years ago
- Views:
Transcription
1 Direct numerical simulations of breaking waves and adaptive dispersive wave models Stéphane Popinet, Luc Deike, Ken Melville Institut Alembert, CNRS/UPMC, Paris Scripps Institution of Oceanography, UC San Diego February 8, 2016
2 Multiscale fluid mechanics: from interfacial flows... Edgerton, 1952 Breaking wave Savart Plateau Rayleigh instability, Jameson, 1971
3 Introduction Breaking waves Dispersive wave models... to geophysical fluid dynamics Dynamics and chemistry of black smokers, Walter 2012 Tohoku tsunami, 2011 Avalanches of granular material Idealised climate models
4 Outline 1 Introduction: equations, adaptivity, Gerris 2 A numerical study of breaking waves 3 Dispersive wave model: the Serre Green Naghdi equations
5 Outline 1 Introduction: equations, adaptivity, Gerris 2 A numerical study of breaking waves 3 Dispersive wave model: the Serre Green Naghdi equations
6 The Navier Stokes equations Conservation of mass and momentum Equations of state, for example t ρ + (ρu) = 0 ρ( t u + u u) = τ + S τ = p(ρ,...)i + η(ρ, u,...)( u + u T ) Incompressibility u = 0 diagnostic Poisson equation 1 ρ p =... Source terms S = gravity, Coriolis, surface tension, electric stresses,...
7 Unstructured statically refined mesh Adaptive in space Multigrid difficult
8 Regular Cartesian grid Not adaptive Multigrid easy
9 Dynamic refinement using quadtrees Adaptive in space and time Multigrid easy (require storage on non-leaf levels)
10 A natural multi-scale/frequency representation Efficient multigrid solvers for linear/nonlinear systems A large collection of other efficient divide-and-conquer algorithms: spatial indexing, compression etc... Formally linked to wavelets/multifractals ( multiresolution analysis )
11 A natural multi-scale/frequency representation Efficient multigrid solvers for linear/nonlinear systems A large collection of other efficient divide-and-conquer algorithms: spatial indexing, compression etc... Formally linked to wavelets/multifractals ( multiresolution analysis )
12 A natural multi-scale/frequency representation Efficient multigrid solvers for linear/nonlinear systems A large collection of other efficient divide-and-conquer algorithms: spatial indexing, compression etc... Formally linked to wavelets/multifractals ( multiresolution analysis )
13 Why adaptivity? Scaling of solution cost Number of degrees of freedom scales like C 4 (4 = 3 spatial dimensions + time) Moore s law Computing power doubles every two years combined with the above scaling gives Spatial resolution of climate models doubles every eight years
14 Does this work? 1e+09 1e+08 # degrees of freedom 1e+07 1e ECMWF MetOffice resolution doubles every 8 years resolution doubles every 10 years Year
15 Gerris A framework to solve partial differential equations on adaptive quad/octree meshes First released as Free Software in 2001 Free Software ideals and practice fit well with Open Science Reproductible results from publications Release early, release often no dark secrets General Public License ensures perennity of contributions A simple but powerful interface for research in fluid mechanics Fast setup and exploration of numerical experiments Multiphysics extensions A large community of users and developers share expertise on fluid mechanics in general Some numbers: > 400 subscribers to the mailing lists 10 active developers > 50 papers published in 2015
16 Outline 1 Introduction: equations, adaptivity, Gerris 2 A numerical study of breaking waves 3 Dispersive wave model: the Serre Green Naghdi equations
17 Geometric Volume-Of-Fluid on quad/octrees t c + u c = 0 C t t C a b C a C b Popinet, JCP 228, 2009
18 Spurious currents (circa 1995) Numerical imbalance between pressure gradient and surface tension
19 Balanced-force Continuum-Surface-Force model 0 = p + σκnδ S Renardy & Renardy, JCP, 2001; Francois et al, JCP, 2006 Continuum-Surface-Force (CSF) model: σκnδ s σκ C Equilibrium if σκ d C = d p Can κ be computed accurately from the discretisation of C?
20 Progress Amplitude of spurious currents as a function of spatial resolution
21 Context: ocean atmosphere interactions Global scale fluxes
22 Context: ocean atmosphere interactions Start with smaller waves (λ O(10) cm)
23 Independent parameters for gas liquid Stokes waves 1 Initial conditions: third-order Stokes waves of varying steepness ak (with a the amplitude and k the wave number) 2 Density ratio: set to air water value (1.2/1000) 3 Dynamic viscosity ratio: set to air water value (1/51) 4 Reynolds number Re = gλ 3 /ν w : set to (λ 5 cm) 5 Bond number Bo = ρga 2 /γ: varying from 1 to
24 Gravity capillary wave Bo = 5, ak = 0.3
25 Spilling-overturning breaker Bo = 15, ak = 0.45
26 Plunging breaker Bo = 1000, ak = 0.55
27 Systematic parametric study Criteria for the different regimes Spilling breaker θ 90 Overturning breaker θ 180 Re = , air water density and viscosity ratios Boundaries found iteratively using dichotomy on ak or Bo
28 State diagram
29 Breaking criteria Observation: the spilling and breaking boundaries scale like ɛ B = (ak) B (1 + Bo) 1/3
30 Breaking criteria Observation: the spilling and breaking boundaries scale like ɛ B = (ak) B (1 + Bo) 1/3 Relationship between vorticity, viscous stress tensor and pressure (taking the divergence of Navier Stokes) Γ 2 D 2 = 2 p ρ, at High Re Γ2 2 p ρ
31 Breaking criteria Observation: the spilling and breaking boundaries scale like ɛ B = (ak) B (1 + Bo) 1/3 Relationship between vorticity, viscous stress tensor and pressure (taking the divergence of Navier Stokes) Γ 2 D 2 = 2 p ρ, at High Re Γ2 2 p ρ The pressure due to capillary effects is 2 p γk 4 a
32 Breaking criteria Observation: the spilling and breaking boundaries scale like ɛ B = (ak) B (1 + Bo) 1/3 Relationship between vorticity, viscous stress tensor and pressure (taking the divergence of Navier Stokes) Γ 2 D 2 = 2 p ρ, at High Re Γ2 2 p ρ The pressure due to capillary effects is 2 p γk 4 a The vorticity is (Longuet-Higgins, 1992) with ω 2 = gk + γ/ρk 3 Γ = 2(ak) 2 ω
33 Breaking criteria Putting everything together gives ) (ak) 4 gk (1 + γk2 ρg γk4 a ρ
34 Breaking criteria Putting everything together gives ) (ak) 4 gk (1 + γk2 ρg γk4 a ρ The Bond number is Bo = ρg/(γk 2 ), so that (ak) (1 + Bo) 1/3
35 Energy dissipation: different regimes Theoretical 0.9 Stable ak=0.25 Total energy Spilling ak=0.35 Spilling-overturning ak= Plunging ak= Time Bo = 10
36 Statistical modelling of breaking wave dissipation Dissipation rate per unit length of breaking crest ɛ l = bρc 5 /g with b the non-dimensional breaking parameter (Duncan, 1981, Phillips, 1985)
37 Statistical modelling of breaking wave dissipation Dissipation rate per unit length of breaking crest ɛ l = bρc 5 /g with b the non-dimensional breaking parameter (Duncan, 1981, Phillips, 1985) Semi-empirical relationship for b (Drazen et al, 2008) b = 0.4(S 0.08) 5/2 with S = ak the slope of the wave.
38 Statistical modelling of breaking wave dissipation Dissipation rate per unit length of breaking crest ɛ l = bρc 5 /g with b the non-dimensional breaking parameter (Duncan, 1981, Phillips, 1985) Semi-empirical relationship for b (Drazen et al, 2008) b = 0.4(S 0.08) 5/2 with S = ak the slope of the wave. We fit the energy decay as E(t) = E 0 exp( ζt) which gives with c 2 = g/k + (γ/rho)k b = ɛ l g/ρc 5 = E 0 ζg/ρc 5
39 Breaking parameter as a function of wave slope Red curve b = 0.4(ɛ 0.08) 5/2
40 3D plunging breaker Bo = , 0.35 < ak < 0.45
41 Adaptive mesh refinement Mass conservation relative error is less than 1/25 000
42 Bubble formation and evolution
43 Bubble formation and evolution
44 Bubble formation and evolution
45 Bubble formation and evolution
46 Bubble formation and evolution
47 Bubble formation and evolution
48 Bubble size distribution Garett et al (2000): steady-state model N (r) Qɛ 1/3 r 10/3 with N (r) the size distribution per unit volume, ɛ the time-averaged turbulence dissipation rate and Q the constant air flow rate per unit volume of water. Experimental measurements
49 Evolution of bubble size distribution (from simulations) Evolution of N(r, t) Decay rate and bubble rise velocity
50 Time-averaged bubble size distributions N(r, t) r 10/3 exp [ 2 ( r ) ] 2 r m (t)
51 Conclusions (on breaking waves) Simple breaking criterion for gravity capillary Stokes waves Classical wave dissipation empirical models work also in the capillary limit Dissipative capillary effects are important for the global energy budget of air-water waves (at all scales!) We can capture (just!) the inertial bubble breakup cascade for 3D breaking waves Deike, Popinet & Melville, Capillary effects on wave breaking, 769: , JFM, Deike, Melville & Popinet, Air entrainment and bubble statistics in three-dimensional breaking waves, under review, JFM, 2016 Basilisk: a new, simpler, more efficient, framework for quadtree-adaptive numerical methods (basilisk.fr)
52 Outline 1 Introduction: equations, adaptivity, Gerris 2 A numerical study of breaking waves 3 Dispersive wave model: the Serre Green Naghdi equations
53 The Saint-Venant equations Conservative form t q = h hu hv Ω qdω =, f(q) = f(q) nd Ω hg z Ω hu hv hu gh2 huv huv hv gh2 System of conservation laws (with source terms) Ω Analogous to the 2D compressible Euler equations (with γ = 2) Hyperbolic characteristic solutions Godunov-type (colocated) 2nd-order finite-volume, shock-capturing Wetting/drying, positivity, lake-at-rest balance: scheme of Audusse et al (2004)
54 The Serre Green Naghdi equations with t η + (h V ) = 0 t V + ɛ V V + η = D (I + µt )(D) = µq( V ) T (D) = (h 3 D) +... and Q( V ) a (complicated) function of the first- and second-derivatives of V.
55 The Serre Green Naghdi equations with t η + (h V ) = 0 t V + ɛ V V + η = D (I + µt )(D) = µq( V ) T (D) = (h 3 D) +... and Q( V ) a (complicated) function of the first- and second-derivatives of V. No source term in the mass equation
56 The Serre Green Naghdi equations with t η + (h V ) = 0 t V + ɛ V V + η = D (I + µt )(D) = µq( V ) T (D) = (h 3 D) +... and Q( V ) a (complicated) function of the first- and second-derivatives of V. No source term in the mass equation Requires the inversion of a (spatially-coupled), time-dependent, 2nd-order linear system for (vector) D
57 The Serre Green Naghdi equations with t η + (h V ) = 0 t V + ɛ V V + η = D (I + µt )(D) = µq( V ) T (D) = (h 3 D) +... and Q( V ) a (complicated) function of the first- and second-derivatives of V. No source term in the mass equation Requires the inversion of a (spatially-coupled), time-dependent, 2nd-order linear system for (vector) D Can be recast into two scalar tridiagonal systems (on regular grids). Lannes and Marche, 2015, JCP
58 Geometric multigrid Fedorenko (1961), Brandt (1977) Convergence acceleration technique for iterative solvers e.g. Gauss Seidel converges in O(λ/ ) iterations wavelength decomposition of the problem on different grids 1 Given an initial guess u 2 Compute residual on fine grid: R = β L(u ) 3 Restrict residual to coarser grid: R R 2 4 Solve on coarse grid: L(δu 2 ) = R 2 5 Prolongate the correction onto fine grid: δu 2 δu 6 Smooth the correction (using e.g. Gauss Seidel iterations) 7 Correct the initial guess: u = u + δu Full multigrid has optimal computational cost O(N) Similar to Fourier (frequency domain) and closely-related to wavelet decomposition of the signal
59 Basilisk: a new quadtree-adaptive framework Free Software (GPL): basilisk.fr Principal objectives: Precision Simplicity Performance Generalised Cartesian grids : Cartesian schemes are turned seamlessly into quadtree-adaptive schemes Basic abstraction: operations only on local stencils a[ 1,1] a[0,1] a[1,1] a[ 1,0] a[0,0] a[1,0] a[ 1, 1] a[0, 1] a[1, 1] Code example: b = 2 a f o r e a c h ( ) b [ 0, 0 ] = ( a [ 0, 1 ] + a [ 1, 0 ] + a [0, 1] + a [ 1,0] 4. a [ 0, 0 ] ) / sq ( D e l t a ) ;
60 Restriction v o i d r e s t r i c t i o n ( s c a l a r v ) { v [ ] = ( f i n e ( v, 0, 0 ) + f i n e ( v, 1, 0 ) + f i n e ( v, 0, 1 ) + f i n e ( v, 1, 1 ) ) / 4. ; }
61 Prolongation v o i d p r o l o n g a t i o n ( s c a l a r v ) { / b i l i n e a r i n t e r p o l a t i o n from p a r e n t / v [ ] = ( 9. c o a r s e ( v, 0, 0 ) + 3. ( c o a r s e ( v, c h i l d. x, 0 ) + c o a r s e ( v, 0, c h i l d. y ) ) + c o a r s e ( v, c h i l d. x, c h i l d. y ) ) / 1 6. ; }
62 Boundary conditions Guarantees stencil consistency independently of neighborhood resolution active points restriction prolongation v o i d boundary ( s c a l a r v, i n t l e v e l ) { f o r ( i n t l = l e v e l 1 ; l <= 0 ; l ) f o r e a c h l e v e l ( l ) r e s t r i c t i o n ( v ) ; f o r ( i n t l = 0 ; l <= l e v e l ; l ++) f o r e a c h h a l o l e v e l ( l ) p r o l o n g a t i o n ( v ) ; }
63 Generic multigrid implementation in Basilisk v o i d m g c y c l e ( s c a l a r a, s c a l a r r e s, s c a l a r dp, v o i d ( r e l a x ) ( s c a l a r dp, s c a l a r r e s, i n t depth ), i n t n r e l a x, i n t m i n l e v e l ) { / r e s t r i c t r e s i d u a l / f o r ( i n t l = depth ( ) 1 ; l <= m i n l e v e l ; l ) f o r e a c h l e v e l ( l ) r e s t r i c t i o n ( p o i n t, r e s ) ; / m u l t i g r i d t r a v e r s a l / f o r ( i n t l = m i n l e v e l ; l <= depth ( ) ; l ++) { i f ( l == m i n l e v e l ) / i n i t i a l g u e s s on c o a r s e s t l e v e l / f o r e a c h l e v e l ( l ) dp [ ] = 0. ; e l s e / p r o l o n g a t i o n from c o a r s e r l e v e l / f o r e a c h l e v e l ( l ) p r o l o n g a t i o n ( dp ) ; boundary ( dp, l ) ; / r e l a x a t i o n / f o r ( i n t i = 0 ; i < n r e l a x ; i ++) { r e l a x ( dp, r e s, l ) ; boundary ( dp, l ) ; } } / c o r r e c t i o n / f o r e a c h ( ) a [ ] += dp [ ] ; }
64 Application to Poisson equation 2 a = b The relaxation operator is simply v o i d r e l a x ( s c a l a r a, s c a l a r b, i n t l ) { f o r e a c h l e v e l ( l ) a [ ] = ( a [ 1, 0 ] + a [ 1,0] + a [ 0, 1 ] + a [0, 1] sq ( D e l t a ) b [ ] ) / 4. ; } The corresponding residual is v o i d r e s i d u a l ( s c a l a r a, s c a l a r b, s c a l a r r e s ) { f o r e a c h ( ) r e s [ ] = b [ ] ( a [ 1, 0 ] + a [ 1,0] + a [ 0, 1 ] + a [0, 1] 4. a [ ] ) / sq ( D e l t a ) ; }
65 The Serre Green Naghdi residual [( h α d h 2 2 xy z b + xη y z b ) h 3 xd x + h [α d ( xη xz b + h2 ) 2x z b α ( d 3 x ) D y + h2 y z b xd y h2 3 2 xy D y h y D y ] + 1 D x+ ( xh + 1 )] 2 xz b = b x v o i d r e s i d u a l ( v e c t o r D, v e c t o r b, v e c t o r r e s ) { f o r e a c h ( ) f o r e a c h d i m e n s i o n ( ) { double hc = h [ ], dxh = dx ( h ), dxzb = dx ( zb ), d x e t a = dxzb + dxh ; double h l 3 = ( hc + h [ 1, 0 ] ) / 2. ; h l 3 = cube ( h l 3 ) ; double hr3 = ( hc + h [ 1, 0 ] ) / 2. ; hr3 = cube ( hr3 ) ; } } r e s. x [ ] = b. x [ ] ( a l p h a d / 3. ( hr3 D. x [ 1, 0 ] + h l 3 D. x [ 1,0] ( hr3 + h l 3 ) D. x [ ] ) / sq ( D e l t a ) + hc ( a l p h a d ( d x e t a dxzb + hc /2. d2x ( zb ) ) + 1. ) D. x [ ] + a l p h a d hc (( hc /2. d2xy ( zb ) + d x e t a dy ( zb )) D. y [ ] + hc /2. dy ( zb ) dx (D. y ) sq ( hc ) / 3. d2xy (D. y ) hc dy (D. y ) ( dxh + dxzb / 2. ) ) ) ;
66 Are dispersive terms important? (for large-scale tsunamis) Saint-Venant Serre Green Naghdi 30 minutes after the event Colorscale ±2 metres
67 Are dispersive terms important? (for large-scale tsunamis) Saint-Venant Serre Green Naghdi 60 minutes after the event Colorscale ±2 metres
68 Are dispersive terms important? (for large-scale tsunamis) Saint-Venant Serre Green Naghdi 90 minutes after the event Colorscale ±2 metres
69 Are dispersive terms important? (for large-scale tsunamis) Saint-Venant Serre Green Naghdi 120 minutes after the event Colorscale ±2 metres
70 Are dispersive terms important? (for large-scale tsunamis) Saint-Venant Serre Green Naghdi 150 minutes after the event Colorscale ±2 metres
71 Are dispersive terms important? Saint-Venant Serre Green Naghdi 30 minutes after the event
72 Are dispersive terms important? Saint-Venant Serre Green Naghdi 60 minutes after the event
73 Are dispersive terms important? Saint-Venant Serre Green Naghdi 90 minutes after the event
74 Are dispersive terms important? Saint-Venant Serre Green Naghdi 120 minutes after the event
75 Are dispersive terms important? Saint-Venant Serre Green Naghdi 150 minutes after the event
76 Conclusions The Serre Green Naghdi dispersive model can be implemented as a momentum source added to an existing Saint-Venant model Preserves the well-balancing, positivity of water depth (wetting/drying) of the orginal scheme Multigrid is simple and efficient for inverting the SGN operator on adaptive quadtree grids. Popinet, 2015, JCP The effective number of degrees of freedom of the physical problem scales like C d with d the effective (or fractal or information ) dimension. This leads to large gains in computational cost for a given error threshold for a wide range of problems (including wave dynamics). Current developments 4th-order quadtree-adaptive discretisations MPI parallelism (dynamic load-balancing etc...) GPUs
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
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 informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationTransport equation cavitation models in an unstructured flow solver. Kilian Claramunt, Charles Hirsch
Transport equation cavitation models in an unstructured flow solver Kilian Claramunt, Charles Hirsch SHF Conference on hydraulic machines and cavitation / air in water pipes June 5-6, 2013, Grenoble, France
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationNumerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
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 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 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 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 informationSimulation of free surface fluids in incompressible dynamique
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
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 informationNumerical dispersion and Linearized Saint-Venant Equations
Numerical dispersion and Linearized Saint-Venant Equations M. Ersoy Basque Center for Applied Mathematics 11 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 The Saint-Venant equations
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
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 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 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 informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
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 informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
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 information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
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 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 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 information1. INTRODUCTION TO CFD SPRING 2018
1. INTRODUCTION TO CFD SPRING 018 1.1 What is computational fluid dynamics? 1. Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
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 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 informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
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 informationA numerical study of breaking waves
PHYSICS OF FLUIDS VOLUME 16, NUMBER 7 JULY 2004 Chiyoon Song and Ana I. Sirviente a) Department of Naval Architecture and Marine Engineering, The University of Michigan, Ann Arbor, Michigan 48109-2145
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationModified Serre Green Naghdi equations with improved or without dispersion
Modified Serre Green Naghdi equations with improved or without dispersion DIDIER CLAMOND Université Côte d Azur Laboratoire J. A. Dieudonné Parc Valrose, 06108 Nice cedex 2, France didier.clamond@gmail.com
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 informationAn Experimental Investigation of the Effects of Surfactants on Spilling Breakers
An Experimental Investigation of the Effects of Surfactants on Spilling Breakers X. Liu and J. H. Duncan University of Maryland J. Kelly and G. M. Korenowski RPI Supported by the National Science Foundation
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 informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationThomas Pierro, Donald Slinn, Kraig Winters
Thomas Pierro, Donald Slinn, Kraig Winters Department of Ocean Engineering, Florida Atlantic University, Boca Raton, Florida Applied Physics Laboratory, University of Washington, Seattle, Washington Supported
More informationMA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies
MA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies 13th February 2015 Jorge Lindley email: J.V.M.Lindley@warwick.ac.uk 1 2D Flows - Shear flows Example 1. Flow over an inclined plane A
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 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 informationLecture: Wave-induced Momentum Fluxes: Radiation Stresses
Chapter 4 Lecture: Wave-induced Momentum Fluxes: Radiation Stresses Here we derive the wave-induced depth-integrated momentum fluxes, otherwise known as the radiation stress tensor S. These are the 2nd-order
More informationAlongshore Momentum Balance: Currents
Chapter 16 Alongshore Momentum Balance: Currents Two assumptions are necessary to get a simple equation for v. The first is that the flow is steady so that time derivatives can be neglected. Second, assume
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More information1. INTRODUCTION TO CFD SPRING 2019
1. INTRODUCTION TO CFD SPRING 2019 1.1 What is computational fluid dynamics? 1.2 Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationFluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow
OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail
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 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 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 informationIncompressible MHD simulations
Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline
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 informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationChapter 5. Methods for Solving Elliptic Equations
Chapter 5. Methods for Solving Elliptic Equations References: Tannehill et al Section 4.3. Fulton et al (1986 MWR). Recommended reading: Chapter 7, Numerical Methods for Engineering Application. J. H.
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
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 informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationPrediction of Coating Thickness
Prediction of Coating Tickness Jon D. Wind Surface Penomena CE 385M 4 May 1 Introduction Tis project involves te modeling of te coating of metal plates wit a viscous liquid by pulling te plate vertically
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationNumerical simulations of drop impacts
Numerical simulations of drop impacts Christophe Josserand Institut D Alembert, CNRS-UPMC L. Duchemin, Z. Jian, P. Ray and S. Zaleski Numerical simulations of drop impacts Christophe Josserand Institut
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
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 informationTopics in Fluid Dynamics: Classical physics and recent mathematics
Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:
More informationA Coupled VOF-Eulerian Multiphase CFD Model To Simulate Breaking Wave Impacts On Offshore Structures
A Coupled VOF-Eulerian Multiphase CFD Model To Simulate Breaking Wave Impacts On Offshore Structures Pietro Danilo Tomaselli Ph.d. student Section for Fluid Mechanics, Coastal and Maritime Engineering
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
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 informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationCCC Annual Report. UIUC, August 19, Argon Bubble Behavior in EMBr Field. Kai Jin. Department of Mechanical Science & Engineering
CCC Annual Report UIUC, August 19, 2015 Argon Bubble Behavior in EMBr Field Kai Jin Department of Mechanical Science & Engineering University of Illinois at Urbana-Champaign Introduction Argon bubbles
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
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 informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
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 information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More information