Introduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
|
|
- Milo Whitehead
- 5 years ago
- Views:
Transcription
1 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
2 Introduction to Finite Volume projection methods The finite-volume principle Upwinding and dissipation for convection-dominated flows Divergence constraints from the zero-mach number limit Projection and auxiliary variable schemes
3 Finite Volume Principle Mass Conservation for arbitrary Ω Ω Ω Ω σ n t v n v t (v n) Change of a volume s mass with time Figure 1.1: Change of a volume s mass with time M(t; Ω) = If, in addition, ϱv satisfies the conditions of Gauß integral theorem, (see appendix A), then Ω (ϱv) n dσ = Ω (ϱv) dv,and Ω ( ) ϱ t + (ϱv) dv =0 for arbitrary t2 Gauß domains Ω. (1.4) Ω M(t 2 ; Ω) M(t 1 ; Ω) = ρ(t, x) dv (ρv) n dσ dt This equation can hold, for continuously differentiablet 1 fields Ω ϱ, v and for arbitrary control volumes Ω only, if pointwise the following partial differential equation is satisfied:
4 Finite Volume Principle General Conservation Law for arbitrary Ω Ω Ω Ω σ n t v n v t (v n) Change of a volume s mass with time Figure 1.1: Change of a volume s mass with time Φ(t; Ω) = If, in addition, ϱv satisfies the conditions of Gauß integral theorem, (see appendix A), then Ω (ϱv) n dσ = Ω (ϱv) dv,and Ω ( ) ϱ t + (ϱv) dv =0 for arbitrary t2 Gauß domains Ω. (1.4) Ω Φ(t 2 ; Ω) Φ(t 1 ; Ω) = ρφ(t, x) dv F φ n dσ dt This equation can hold, for continuously differentiablet 1 fields Ωϱ, v and for arbitrary control volumes Ω only, if pointwise the following partial differential equation is satisfied:
5 Finite Volume Principle General Conservation Law for cells of a grid Change of a volume s mass with time Φ(t; Ω) = Ω Φ(t 2 ; Ω) Φ(t 1 ; Ω) = ρφ(t, x) dv t2 t 1 Ω F φ n dσ dt
6 Finite Volume Principle Tasks: Approximation of fluxes F φ Recovery of detailed φ-distributions from cell averages Robust & accurate time integration Key advantage: Based on fundamental conservation principle valid, in principle, even for under-resolved solutions example: shock capturing schemes for hyperbolic systems Disadvantages: Recovery difficult on general unstructured meshes (but: see, e.g., recent work by M. Dumbser et al.) Rigorous theory less advanced than that of FEM methods (but: see recent developments of Discountinuous Galerkin methods)
7 Introduction to Finite Volume projection methods The finite-volume principle Upwinding and dissipation for convection-dominated flows Divergence constraints from the zero-mach number limit Projection and auxiliary variable schemes
8 Upwinding and Dissipation Burgers Equation or u t + ( ) u 2 /2 = 0 x ( t + u x ) u = 0 smooth solutions: u constant along characteristics dx ch /dt = u. characteristic-based solution u(t, x) = u(0, x t u(t, x)) There cannot be any over/undershoots beyond the initial data!!
9 Upwinding and Dissipation Upwind scheme: F u = ( u 2 /2 ) upwind Solution at time t = 0.1 Solution at time t = 0.2 Solution at time t = q 0 q 0 q x Central scheme: F u = 1 2 Solution at time t = 0.1 [( u 2 /2 ) L + ( u 2 /2 ) R] x Solution at time t = x q q 0 your guess what happens next x x
10 Upwinding and Dissipation Conclusion: Upwinding introduces physically correct dissipation
11 Introduction to Finite Volume projection methods The finite-volume principle Upwinding and dissipation for convection-dominated flows Divergence constraints from the zero-mach number limit Projection and auxiliary variable schemes
12 Divergence Constraints Incompressible flow with variable density Leading order result: p = P + M 2 p Consequences: ρe = p γ 1 + M2ρv2 2 P γ 1 (M 0) (ρe) t + ([ρe + p]v) = 0 v (0) = 0 ρ t + (ρv) = 0 ρ (0) t + v (0) ρ (0) = 0 (ρv) t + (ρv v) + 1 M2 p = 0 (ρv)(0) t + (ρv v) (0) + p = 0 P and p must be maintained as primary variables in a low-mach numerical implementation
13 Introduction to Finite Volume projection methods The finite-volume principle Upwinding and dissipation for convection-dominated flows Divergence constraints from the zero-mach number limit Projection and auxiliary variable schemes
14 Projection and Aux-Variable schemes Compressible flow equations: ρ t + (ρv) = 0 (ρv) t + (ρv v) + ΓP π = ρgk P t + (P v) = 0 P = p 1 γ = ρθ, π = p/p, Γ = c p /R
15 Projection and Aux-Variable schemes Pseudo-incompressible model: (Durran (1988)) ρ t + (ρv) = 0 (ρv) t + (ρv v) + ΓP π = ρgk + (P v) = 0 P P (z), ρθ = P (z), θ = θ(z) + θ
16 Projection and Aux-Variable schemes Auxiliary variable projection scheme predictor step Solve over t by second-order (upwind FV) scheme: ρ t + (ρv) = 0 (ρv) t + (ρv v) + ΓP π n = ρgk Provides: P t + (P v) = 0 second-order advection of {θ, v} no divergence control imposed! second-order gravity term first-order pressure gradient effect P P n = O(( t) 2 ) incompr. case, 2nd-ord.: Schneider et al., JCP 99; 4th-ord.: Kadioglu, Minion, Klein, JCP 08
17 Projection and Aux-Variable schemes Auxiliary variable projection scheme predictor step Solve over t by second-order (upwind FV) scheme: ρ t + (ρv) = 0 (ρv) t + (ρv v) + ΓP π n = ρgk Provides: P t + (P v) = 0 second-order advection of {θ, v} no divergence control imposed! second-order gravity term first-order pressure gradient effect P P n = O(( t) 2 )
18 Projection and Aux-Variable schemes Auxiliary variable projection scheme predictor step Solve over t by second-order (upwind FV) scheme: ρ t + (ρv) = 0 (ρv) t + (ρv v) + ΓP π n = ρgk Provides: P t + (P v) = 0 second-order advection of {θ, v} no divergence control imposed! second-order gravity term first-order pressure gradient effect P P n = O(( t) 2 )
19 Projection and Aux-Variable schemes Auxiliary variable projection scheme MAC-projection Constraint (P v) n+1 2 = 0 δπ = π n+1 2 π n Advection velocity correction (at cell interfaces) Divergence control (P v) n+1 2 = (P v) t 2 P θ δπ 1 P n (P θ δπ) = 2 1 ( t) 2 P (P P n ) n Post-correction of all advective flux contributions
20 Projection and Aux-Variable schemes Auxiliary variable projection scheme MAC-projection MAC-projection controls advection fluxes across grid cell interfaces Advection velocities result from some interpolation from the cell-centers Danger: pressure-velocity decoupling Previous approaches: Rhie-Chow stabilization Exact projection Approximate projection compromises mass continuity non-compact stencils inexact divergence control
21 Projection and Aux-Variable schemes Alternative route : cell-centered exact projection Exact control of divergence on dual cells Piecewise linear, discontinuous ansatz functions for momentum on primary cells Bi-/trilinear pressure ansatz functions S. Vater, R.K., Num. Math., 113, (2009)
22 Projection and Aux-Variable schemes Stabilizing Second projection : inf-sup-stable Petrov-Galerkin scheme for generalized saddle-point problem Find (u, p) (X 2 M 1 ), such that a(u, v) + b 1 (v, p) = f, v v X 1 b 2 (u, q) = g, q q M 2 where a(u, v) = (u, v) b 1 (v, p) = (v, p) b 2 (u, q) = ( u, q) (result currently available for 2D shallow water / homentropic gasdynamics) S. Vater, R.K., Num. Math., 113, (2009); using Nicolaïdes (1982), Bernardi et al. (1988)
23 Projection and Aux-Variable schemes Summary of auxiliary-variable projection method Predictor auxiliary unconstrained system second-order upwind finite volume scheme well-balanced discretization of gravity term MAC-projection controls advective flux divergence brings P back to P (z) cell-centered projection stabilizes checkerboard p v modes generates second-order accuracy w.r.t. pressure Background nowhere invoked explicitly!!
24 Projection and Aux-Variable schemes Last remark and warning: Projection does not preserve the no-slip wall boundary condition! Therefore, the quality of maintaining this b.c. depends on the accuracy of the predictor step. See BCM (2001) for a detailed discussion. Brown, Cortez & Minion, JCP, 168, (2001)
Finite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
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 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 informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
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 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 informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationDiscretization Techniques for Surface Tension Forces in Fully Conservative Two-Phase Flow Finite Volume Methods
Discretization Techniques for Surface Tension Forces in Fully Conservative Two-Phase Flow Finite Volume Methods Matthias Waidmann 1, Stephan Gerber 2 Michael Oevermann 3, Rupert Klein 1 1 Freie Universität
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 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 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 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 informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
More informationExtension of finite volume compressible flow solvers to multi-dimensional, variable density zero Mach number flows
Extension of finite volume compressible flow solvers to multi-dimensional, variable density zero Mach number flows T. Schneider yz and N. Botta yz 2 and K. J. Geratz ffl and R. Klein yfiz y Freie Universität
More informationFinite element exterior calculus framework for geophysical fluid dynamics
Finite element exterior calculus framework for geophysical fluid dynamics Colin Cotter Department of Aeronautics Imperial College London Part of ongoing work on UK Gung-Ho Dynamical Core Project funded
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 informationComparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models
Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationFailsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model
Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model Matthias Möller 1 joint work with: Dmitri Kuzmin 2 1 Institute of Applied Mathematics
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
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 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 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 informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationDiscretisation of the pseudo-incompressible Equations with implicit LES for Gravity Wave Simulation DACH2010. Felix Rieper, Ulrich Achatz
Discretisation of the pseudo-incompressible Equations with implicit LES for Gravity Wave Simulation DACH2010, Ulrich Achatz Institut für Atmosphäre und Umwelt Aim of the project Analysis: Propagation and
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 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 informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
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 informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
More informationOn the regime of validity of sound-proof model equations for atmospheric flows
On the regime of validity of sound-proof model equations for atmospheric flows Rupert Klein Mathematik & Informatik, Freie Universität Berlin ECMWF, Non-hydrostatic Modelling Workshop Reading, November
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 informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
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 informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
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 informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
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 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 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 informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
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 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 informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
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 informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
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 informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
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 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 informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationSkew-Symmetric Schemes for compressible and incompressible flows
Skew-Symmetric Schemes for compressible and incompressible flows Julius Reiss Fachgebiet für Numerische Fluiddynamik TU Berlin CEMRACS J. Reiss (TU Berlin) Skew-Symmetric Schemes August 8 th, 202 / 46
More informationColumnar Clouds and Internal Waves
Columnar Clouds and Internal Waves Daniel Ruprecht 1, Andrew Majda 2, Rupert Klein 1 1 Mathematik & Informatik, Freie Universität Berlin 2 Courant Institute, NYU Tropical Meteorology Workshop, Banff April
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationThermodynamic Consistency of a Pseudoincompressible Approximation for General Equations of State
MARCH 202 K L E I N A N D P A U L U I S 96 Thermodynamic Consistency of a Pseudoincompressible Approximation for General Equations of State RUPERT KLEIN Mathematik und Informatik, Freie Universität Berlin,
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
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 informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
More informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
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 informationarxiv: v2 [math.na] 9 Apr 2018
oname manuscript o. will be inserted by the editor) The BR1 Scheme is Stable for the Compressible avier-stokes Equations Gregor J. Gassner Andrew R. Winters Florian J. Hindenlang David A. Kopriva arxiv:1704.03646v
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationAn efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes
An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes Christian Klingenberg, Frank Pörner, Yinhua Xia Abstract
More informationHigh-resolution finite element schemes for an idealized Z-pinch implosion model
High-resolution finite element schemes for an idealized Z-pinch implosion model Matthias Möller a, Dmitri Kuzmin b, John N. Shadid c a Institute of Applied Mathematics (LS III), Dortmund University of
More informationNumerik 2 Motivation
Numerik 2 Motivation P. Bastian Universität Heidelberg Interdisziplinäres Zentrum für Wissenschaftliches Rechnen Im Neuenheimer Feld 368, D-69120 Heidelberg email: Peter.Bastian@iwr.uni-heidelberg.de April
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 informationHIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 HIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
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 informationSolving the Navier-Stokes Equations
FMIA F. Moukalled L. Mangani M. Darwish An Advanced Introduction with OpenFOAM and Matlab This textbook explores both the theoretical oundation o the Finite Volume Method (FVM) and its applications in
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationUsing the sound-proof limit for balanced data initialization
Using the sound-proof limit for balanced data initialiation Rupert Klein, Tommaso Benacchio, Warren O Neill FB Mathematik & Informatik, Freie Universität Berlin Arnimallee 6, 14195 Berlin, Germany rupert.klein@math.fu-berlin.de
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
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 Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
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 information