Angular momentum preserving CFD on general grids
|
|
- Ethan Flowers
- 5 years ago
- Views:
Transcription
1 B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration Emmanuel Labourasse (CEA) Angular momentum preserving CFD on general grids p. 1 / 29
2 Section 1 Angular momentum preserving CFD on general grids p. 2 / 29
3 Question for conservation laws and Riemann addicts Compressible CFD is based on conservation laws tρ + (ρu) = 0, t(ρu) + (ρu 2 + pi) = 0, t(ρe) + (pu) = 0. 1D Riemann do at the same time excellently and poorly for multid flows. What about using angular momentum numerical conservation law which is a physically sound conservation law? Issue addressed already in astrophysics, but not that often in CFD. See however CFD for helicopters (topic I know almost nothing). Angular momentum preserving CFD on general grids p. 2 / 29
4 Angular momentum in nature, engineering and sport (internet provided) Angular momentum preserving CFD on general grids p. 3 / 29
5 Angular momentum The angular momentum w = u x satisfies an additional conservation laws t(ρw) + (ρu u x) + curl(px) = 0. The vorticity ω = u can be computed from the angular momentum. Let x 0 be a given point in the domain and consider w 0(x) = w(x) u(x) x 0 = u(x) (x x 0). Then w 0(x) = u(x) (x x 0) 2 u(x). Therefore one has the identity ω(x 0) = u(x 0) = 1 2 w0(x) Angular momentum preserving CFD on general grids p. 4 / 29
6 Some references (vorticity/am) : Roe and Morton, Preserving vorticity in finite-volume schemes, Finite volumes for complex applications II, , Hermes Sci. Publ., Paris, $ Caramana and Loubère, Curl-q : A vorticity damping artificial viscosity for essentially irrotational Lagrangian hydrodynamics calculations, $ Dukowicz-Meltz, Vorticity errors in multidimensional lagrangian codes Journal of Computational Physics, Volume 99, Issue 1, March Käppeli, Mishra, Structure preserving schemes. In Seminar for Applied Mathematics, ETH Zürich, Zhang-Zhang-C.W. Shu, Multistage interaction of a shock wave and a strong vortex, Phys. Fluids, Rault-Chiavassa-Donat, Shock-Vortex Interactions at High Mach Numbers, Oort, Angular in the Atmosphere-Ocean-Solid Earth System, Bulletin of the American Meteorological Society, 70 (10), 1, , Angular momentum preserving CFD on general grids p. 5 / 29
7 More references MultiD Riemann - Bourgeade-Le Floch-Raviart : An asymptotic expansion for the solution of the generalized Riemann problem. Part 2 : application to the equations of gas dynamics, Ben-Artzi-Falcovitz, Generalized Riemann Problems in Computational Fluid Dynamic, Brio-Zakharian-Webb, Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics, Balsara, Multidimensional Riemann Problem with Self-Similar Internal Structure-Part I. Application to Hyperbolic Conservation Laws on Structured Meshes, JCP 277 (2014) Low Mach flows - Guillard-Voizat, On the behaviour of upwind schemes in the low Mach number limit, Compters and Fluids, 28, Dellacherie,... Angular momentum preserving CFD on general grids p. 6 / 29
8 Meshes and computations - Multimaterial computation is courtesy from Delpino (JCP 2011). - 3D is from Chertkov-Levedev on internet. Angular momentum preserving CFD on general grids p. 7 / 29
9 First method We assume the thermodynamics is discretized cell-centered-fv-p 0 ρ j, T j, ε j = C v T j, p j = (γ 1)ρ j ε j R. Probably the velocity will be discretized cell-centered-fv-p 0 : u j R d. What can we do for the angular momentum? First idea : Angular momentum satisfies an independent conservation law, so one can use additional but not fully independent unknowns. Velocity locally in the cell V j made of all solid body displacement v j (x) = a j + b j x Define u j = 1 v j (x) V j V j = a j + b j x j, w j = 1 v j (x) x, V j V j = (a j + b j x j ) x j + 1 (b j (x x j )) (x x j ) V j V j Angular momentum preserving CFD on general grids p. 8 / 29
10 Inverse formulas Define the positive symmetric matrix 0 < H j = H t j R d d b R d, (H j b, b) := 1 V j One can prove H j is non-singular, even if H j = O(h 2 ). V j b (x x j ) 2. Therefore there is a one to one correspondence between the physical unknowns u j, w j and the couple a j, b j { b j = H 1 j (u j x j w j ), a j = u j b j x j. (1) Angular momentum preserving CFD on general grids p. 9 / 29
11 Second method : link with DG Define P n = u = 0 j 1 + +j d n x j x j d d a j1,...,j d a j1,...,j d R d are vectors Rd. The dimension of P n increases with n : dim(p 0 ) = d, dim(p 1 ) = d(d + 1) = d 2 + d and so on. Define a new space P 0 Q = Span {a + Bx} P 1 where a P 0 is arbitrary and the vectorial product is noted with the multiplication by an antisymmetric matrix B t = B. Therefore dim(q) = d + d(d 1) 2 = 1 2 dim(p1 ). The number of extra degrees of freedom with respect to P 0 is N extra = dim(q) dim(p 0 ) = d(d 1). 2 Angular momentum preserving CFD on general grids p. 10 / 29
12 DG procedure Consider the impulse equation DG yields That is d t t(ρu) + M = 0, M = M t = ρu u + pi d. Θ d t Θ ρu ũ + ( M) ũ = 0, ũ Q. Θ ρu ũ M : ũ = (Mũ, n) dσ, ũ Q. Θ Θ Proposition The interior integral vanishes : M : ũ = 0, ũ Q. Θ proof. Indeed ũ = B is an antisymmetric matrix. Since M is symmetric, the contraction vanishes M : B = 0. Angular momentum preserving CFD on general grids p. 11 / 29
13 Equations Take ũ = a. One obtains the conservative inertial momentum equation d t ρu ã = (Mã, n) dσ, ã R d. Θ Θ Take ũ = b x. One obtains the conservative angular momentum equation d t ρu b x = d t ρu x b ( ) = M b x, n dσ, b R d. Θ Θ Θ Both approaches yield the same result. Angular momentum preserving CFD on general grids p. 12 / 29
14 Section 3 Angular momentum preserving CFD on general grids p. 13 / 29
15 Impact is only on Muscl reconstruction Before : most production codes are 2nd order, and based on Muscl techniques. Usually the velocity is reconstructed in the cell V j New method : we decompose v j (x) = u j + u Muscl j (x x j ). u New j = D j + B j, D j = D t j = 1 2 ( ( u Muscl j + u Muscl j ) t ) and B j = B t j provided by the angular momentum degrees of freedom. This is the only place where the AM equation interacts with the usual scheme. Angular momentum preserving CFD on general grids p. 13 / 29
16 Our solver We start from Lagrange CFD ρd tτ u = 0, ρd tu + p = 0, ρd te + (pu) = 0. Our scheme GLACE is a cell-centered Lagrangian scheme, with corner-based Riemann. C jr = xr V j. m j d tτ j (t) = C jr u r, r N (j) m j d tu j (t) = C jr p jr, r N (j) m j d te j (t) = r N (j) C jr u r p jr. where m j, u j and e j are the Lagrangian mass, the velocity and the total energy, and d tϕ ( t + u ) ϕ is the Lagrangian derivative. The nodal unknowns u r (the velocity of the r-th node), and p jr (the nodal pressure) are computed thanks to a nodal Riemann solver for which many different versions exist nowadays. (2) Angular momentum preserving CFD on general grids p. 14 / 29
17 The corner-based Riemann solver It writes where n jr = C jr C jr p jr p j + ρ j c j (u r u j ) n jr = 0, C jr p jr = 0, j C(r) is the normalized corner direction. The first order nodal solver can be rewritten as { Ar u r = j A jr u j + C jr p j, F jr = C jr p j + A jr (u j u r ). In this system, the Glace and Eucclhyd schemes differ ultimately only by the definition of the R d d matrices A jr and of the vector F jr Angular momentum preserving CFD on general grids p. 15 / 29
18 Workflow 1 From the mesh, compute the geometrical features C jr, x j, x r, H j, 2 Compute b j = H, 1 j (u j x j w j ) from which B j is deduced. 3 Limit the reconstructed velocity, and deduce B # j and eventually D # j. 4 Compute v # j (x r ) = u j + u # j (x r x j ) 5 Construct 2nd order nodal pressure p # j (x r ) = p j + p # j (x r x j ). 6 Solve the system with the reconstructed quantities : { A r u # r = j A jr v # j (x r ) + ( C jr p # j (x r ), ) = C jr p # j (x r ) + A jr v # j (x r ) u # r. F # jr 7 Update u # j, w # j and e # j u # j = u j t F # jr m, j r w # j = w j t F # jr m x r. j r e # j = ej t F # jr u# r. m j r (3) 8... Angular momentum preserving CFD on general grids p. 16 / 29
19 Stability/entropy Consider the new AMC scheme with D j = 0 and no limitation on B j. The entropy balance of the Lagrangian scheme writes with (similar to first order scheme) Q j = r m j T j d ts j = Q j + R j (u r v j (x r ) A jr (u r v j (x r )) 0 and R j = m j b j u j d tx j + m j b j (d t(h j )b j /2). Since R j make take any sign, it is in practice controlled by a limiter technique. Angular momentum preserving CFD on general grids p. 17 / 29
20 Section 4 Angular momentum preserving CFD on general grids p. 18 / 29
21 Rotating ring Note : for central forces, the AMC is a local d t(ru θ ) = 0. Angular momentum preserving CFD on general grids p. 18 / 29
22 Rotating ring : standard 2nd Muscl scheme Angular momentum preserving CFD on general grids p. 19 / 29
23 Rotating ring : final mesh with the new method Angular momentum preserving CFD on general grids p. 20 / 29
24 Taylor vortex The mesh is fixed Eulerian, and the scheme is run Lagrange+remapp ρ(x) = (1 ) 1 (γ 1)β2 e 1 r 2 γ 1, p(x) = ρ(x) γ, u(x) = β 8γπ 2 2π e 1 r2 2 x x where β = 5 corresponds to the strength of the vortex, and r = x is the distance to the center. The analytic solution is stationnary 1D plot of the density (left) and the velocity (right) for the Taylor vortex. Angular momentum preserving CFD on general grids p. 21 / 29
25 Convergence table for the Taylor Vortex with amc no amc number of cells L 2 error order L 2 error order The order of our Muscl procedure is not degraded by the introduction of the angular momentum variable. The are even slightly better. Muscl shows more smearing of b = β e 1 r 2 2π 2 than new AMC technique. Angular momentum preserving CFD on general grids p. 22 / 29
26 Unexpected advantage of angular Claim : Angular momentum conservation lowers the mesh imprint on implosion problems. Example of a initial skewed mesh. The center of implosion is at M. Relevant for computational inertial fusion : NIF, ICF,.... Angular momentum preserving CFD on general grids p. 23 / 29
27 Principle Consider I = Ω ρu and W = Ω ρu x Assume for simplicity the boundary conditions are such that I and W are constant. Assume M be the center of implosion : t f u (x M) = 0. Consequently W = I M. such that x, one has This is a linear system with unknown M which cannot be fully recovered. Additional symmetry considerations yield M. The previous example corresponds to W = 0 0 w, I = 0 i 0, M = x M y M 0 = w = ix M = x M = w i. The symmetry assumption is y M = 0 since the center is on the main axis. Angular momentum preserving CFD on general grids p. 24 / 29
28 Results no AMC with AMC O1 O2 Our best result is O2-AMC. Angular momentum preserving CFD on general grids p. 25 / 29
29 Shock with rotation Set up Angular momentum preserving CFD on general grids p. 26 / 29
30 Reference solution AMC yields d dt (ru θ) = 0, d dt = t + u Angular momentum preserving CFD on general grids p. 27 / 29
31 Results no AMC with AMC O1 O2 Our best result is O2-AMC. Angular momentum preserving CFD on general grids p. 28 / 29
32 Conclusions Angular momentum is a natural conservation law. It can be easily incorporated in any standard 2nd order Muscl-FV code. It increases the accuracy. An issue is stability : the decomposition of the velocity gradient does the job in our case. My feeling is it could be also incorporated in most high order Eulerian CFD codes. D. Labourasse : see JCP online Angular momentum preserving CFD on general grids p. 29 / 29
On Fluid Maxwell Equations
On Fluid Maxwell Equations Tsutomu Kambe, (Former Professor, Physics), University of Tokyo, Abstract Fluid mechanics is a field theory of Newtonian mechanics of Galilean symmetry, concerned with fluid
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
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 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 informationNumerical resolution of a two-component compressible fluid model with interfaces
Numerical resolution of a two-component compressible fluid model with interfaces Bruno Després and Frédéric Lagoutière February, 25 Abstract We study a totally conservative algorithm for moving interfaces
More informationTwo-Dimensional Riemann Solver for Euler Equations of Gas Dynamics
Journal of Computational Physics 167, 177 195 (2001) doi:10.1006/jcph.2000.6666, available online at http://www.idealibrary.com on Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics M.
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationWell-balanced and Asymptotic Schemes for Friedrichs systems
B. Després (LJLL- University Paris VI) C. Buet (CEA) Well-balanced and Asymptotic Schemes for Friedrichs systems thanks to T. Leroy (PhD CEA), X. Valentin (PhD CEA), C. Enaux (CEA), P. Lafitte (Centrale
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More informationHigh-order uid-structure coupling for 2D Finite Volume Lagrange-Remap schemes
High-order uid-structure coupling for 2D Finite Volume Lagrange-Remap schemes Gautier Dakin 1 Bruno Després 2 Stéphane Jaouen 1 1 CEA, DAM, DIF, F-91297 Arpajon, France 2 Université Pierre et Marie Curie,
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationWhat is a flux? The Things We Does Know
What is a flux? Finite Volume methods (and others) (are based on ensuring conservation by computing the flux through the surfaces of a polyhedral box. Either the normal component of the flux is evaluated
More informationPositivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders
Positivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders François Vilar a,, Chi-Wang Shu a, Pierre-Henri Maire b a Division of Applied
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 informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
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 informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationSimple waves and a characteristic decomposition of the two dimensional compressible Euler equations
Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations Jiequan Li 1 Department of Mathematics, Capital Normal University, Beijing, 100037 Tong Zhang Institute
More informationEntropy Viscosity Method for the Single Material Euler Equations in Lagrangian Frame
Entropy Viscosity Method for the Single Material Euler Equations in Lagrangian Frame Jean-Luc Guermond a, Bojan Popov a, Vladimir Tomov b, a Department of Mathematics, Texas A&M University 3368 TAMU, College
More informationOn weak solution approach to problems in fluid dynamics
On weak solution approach to problems in fluid dynamics Eduard Feireisl based on joint work with J.Březina (Tokio), C.Klingenberg, and S.Markfelder (Wuerzburg), O.Kreml (Praha), M. Lukáčová (Mainz), H.Mizerová
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationCurl-q : A vorticity damping artificial viscosity for essentially irrotational Lagrangian hydrodynamics calculations
Journal of Computational Physics xxx (2005) xxx xxx Short note Curl-q : A vorticity damping artificial viscosity for essentially irrotational Lagrangian hydrodynamics calculations E.J. Caramana *, R. Loubère
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationA total Lagrangian discontinuous Galerkin discretization of the two-dimensional gas dynamics equations
A total Lagrangian discontinuous Galerkin discretization of the two-dimensional gas dynamics equations François Vilar, Pierre-Henri Maire, Remi Abgrall To cite this version: François Vilar, Pierre-Henri
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
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 informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationAbstract. We develop a new cell-centered control volume Lagrangian scheme for solving Euler
A cell-centered Lagrangian scheme with the preservation of symmetry and conservation properties for compressible fluid flows in two-dimensional cylindrical geometry Juan Cheng and Chi-Wang Shu Abstract
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationOn a simple model of isothermal phase transition
On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows
More informationFinite volumes schemes preserving the low Mach number limit for the Euler system
Finite volumes schemes preserving the low Mach number limit for the Euler system M.-H. Vignal Low Velocity Flows, Paris, Nov. 205 Giacomo Dimarco, Univ. de Ferrara, Italie Raphael Loubere, IMT, CNRS, France
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationNumerical methods for FCI Part IV Multi-temperature fluid models Hele-Shaw models
B. Després+ X. Blanc LJLL-Paris VI+CEA Thanks to same collegues as before plus C. Buet, H. Egly and R. Sentis methods for FCI Part IV Multi-temperature fluid s s B. Després+ X. Blanc LJLL-Paris VI+CEA
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
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 information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationOn the leapfrogging phenomenon in fluid mechanics
On the leapfrogging phenomenon in fluid mechanics Didier Smets Université Pierre et Marie Curie - Paris Based on works with Robert L. Jerrard U. of Toronto) CIRM, Luminy, June 27th 2016 1 / 22 Single vortex
More informationFrontiers in Mathematics. Bruno Després. Numerical Methods. for Eulerian and Lagrangian. Conservation. Laws
Frontiers in Mathematics Bruno Després Numerical Methods for Eulerian and Lagrangian Conservation Laws Frontiers in Mathematics Advisory Editorial Board Leonid Bunimovich (Georgia Institute of Technology,
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
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 information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationFluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers
Fluid-Structure Interaction Problems using SU2 and External Finite-Element Solvers R. Sanchez 1, D. Thomas 2, R. Palacios 1, V. Terrapon 2 1 Department of Aeronautics, Imperial College London 2 Department
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
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 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 informationJ OURNAL OF TURBULENCE
JOT J OURNAL OF TURBULENCE http://jot.iop.org/ Extension of the gridless vortex method into the compressible flow regime Monika Nitsche 1 and James H Strickland 2 1 Department of Mathematics and Statistics,
More informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
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 informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More information1/18/2011. Conservation of Momentum Conservation of Mass Conservation of Energy Scaling Analysis ESS227 Prof. Jin-Yi Yu
Lecture 2: Basic Conservation Laws Conservation Law of Momentum Newton s 2 nd Law of Momentum = absolute velocity viewed in an inertial system = rate of change of Ua following the motion in an inertial
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
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 informationPart 1: Numerical Modeling for Compressible Plasma Flows
Part 1: Numerical Modeling for Compressible Plasma Flows Dongwook Lee Applied Mathematics & Statistics University of California, Santa Cruz AMS 280C Seminar October 17, 2014 MIRA, BG/Q, Argonne National
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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 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 informationAcoustic Energy Estimates in Inhomogeneous Moving Media
Subject: FORUM ACUSTICUM 1999 Abstract Acoustic Energy Estimates in Inhomogeneous Moving Media, NASA Langley Research Center, Hampton, Virginia Mark Farris, Midwestern State University, Wichita Falls,
More informationarxiv: v2 [physics.class-ph] 4 Apr 2009
arxiv:0903.4949v2 [physics.class-ph] 4 Apr 2009 Geometric evolution of the Reynolds stress tensor in three-dimensional turbulence Sergey Gavrilyuk and Henri Gouin Abstract The dynamics of the Reynolds
More informationReceived: 5 September 2017; Accepted: 2 October 2017; Published: 11 October 2017
fluids Article A Finite Element Method for Incompressible Fluid Flow in a Noninertial Frame of Reference and a Numerical Study of Circular Couette Flow with Varying Angular Speed C. S. Jog * and Nilesh
More informationQUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES
1 QUELQUES APPLICATIONS D UN SCHEMA SCHEMA MIXTE-ELEMENT-VOLUME A LA LES, A L ACOUSTIQUE, AUX INTERFACES O. ALLAIN(*), A. DERVIEUX(**), I. ABALAKIN(***), S. CAMARRI(****), H. GUILLARD(**), B. KOOBUS(*****),
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More information