Shock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
|
|
- Elmer Summers
- 5 years ago
- Views:
Transcription
1 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 combines the good properties of the Discontinuous Galerkin method in smooth parts of the flow with the perfect properties of a total variation diminishing finite volume method for resolving shocks without spurious oscillations Due to the sub-cell approach the interior resolution of the Discontinuous Galerkin grid cell is nearly preserved and the number of degrees of freedom remains the same To fully benefit from the high order of the DG method we prefer rather large grid cells This is in contrast to the requirement of resolving small geometrical shapes or other local flow features With the possibilty to locally refine grids including non-conforming interfaces we close the gap between those two contrary requirements The overall implementation is designed to run very efficiently on massively parallel computers
2 Shock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells Matthias Sonntag, Claus-Dieter Munz Institute of Aerodynamics and Gas Dynamics University of Stuttgart SHARK-FV May 23, 2016
3 Motivation high order methods handle shocks, discontinuities, instabilities very efficiently on massive parallel systems Density x Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
4 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
5 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
6 Navier Stokes Equation Conservation Law u t + F (u, u) = 0 in Ω with the vector of conservative variables ρ u = ρv ρe and the convective and viscous fluxes F (u, u) = F c (u) + F v (u, u) Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
7 Discretization & Weak Formulation hexahedral meshes mapping to reference element J(ξ)u t + ξ F = 0 test function + integration Ju t, φ E + F, φ E = 0 partial integration weak formulation Ju t, φ E = f (u, u, n), φ E + F, φ E where f (u, u, n) is a numerical flux function (Riemann solver) Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
8 Discontinuous Galerkin Spectral Element Method (DGSEM) u i,1 ξ 2 1 ansatz functions u = N i,j,k=0 û ijk ψ ijk (ξ), with ψ ijk = l i (ξ 1 )l j (ξ 2 )l k (ξ 3 ) u 1,j where l i (ξ 1 ),... are the 1D Lagrange interpolation polynomials nodes of Gauss quadrature: {ξ i } N i=0 collocation: points of interpolation = points of integration Lagrange property l i (ξ j ) = δ ij i, j test functions = ansatz functions same approximation for fluxes: F = N i,j,k=0 ˆF ijk ψ ijk (ξ) û ij 0 1 ξ 1 Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
9 Time Derivative Integral Ju t, φ E = 1 1 ( 1 N ) J(ξ) û mno ψ mno (ξ) ψ ijk (ξ)dξ 1 dξ 2 dξ 3 dt m,n,o=0 inserting Gauss quadrature + ansatz function: Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
10 Time Derivative Integral Ju t, φ E = 1 1 ( 1 N ) J(ξ) û mno ψ mno (ξ) ψ ijk (ξ)dξ 1 dξ 2 dξ 3 dt m,n,o=0 inserting Gauss quadrature + ansatz function: = N N J(ξ λµν ) û mno l m (ξ dt }{{ λ) l n (ξ } µ) l o (ξν) ψ ijk (ξ λµν )ω λ ω µ ω ν }{{}}{{} λ,µ,ν=0 m,n,o=0 δ mλ δ nµ δ oν Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
11 Time Derivative Integral Ju t, φ E = 1 1 ( 1 N ) J(ξ) û mno ψ mno (ξ) ψ ijk (ξ)dξ 1 dξ 2 dξ 3 dt m,n,o=0 inserting Gauss quadrature + ansatz function: = N N J(ξ λµν ) û mno l m (ξ dt }{{ λ) l n (ξ } µ) l o (ξν) ψ ijk (ξ λµν )ω λ ω µ ω ν }{{}}{{} λ,µ,ν=0 m,n,o=0 δ mλ δ nµ δ oν = N J(ξ λµν )û λµν l i (ξ dt }{{ λ) l j (ξ } µ) l k (ξν) ω λ ω µ ω ν }{{}}{{} λ,µ,ν=0 δ iλ δ jµ δ kν Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
12 Time Derivative Integral Ju t, φ E = 1 1 ( 1 N ) J(ξ) û mno ψ mno (ξ) ψ ijk (ξ)dξ 1 dξ 2 dξ 3 dt m,n,o=0 inserting Gauss quadrature + ansatz function: = N N J(ξ λµν ) û mno l m (ξ dt }{{ λ) l n (ξ } µ) l o (ξν) ψ ijk (ξ λµν )ω λ ω µ ω ν }{{}}{{} λ,µ,ν=0 m,n,o=0 δ mλ δ nµ δ oν = N J(ξ λµν )û λµν l i (ξ dt }{{ λ) l j (ξ } µ) l k (ξν) ω λ ω µ ω ν }{{}}{{} λ,µ,ν=0 δ iλ δ jµ δ kν = J(ξ ijk ) û ijk dt ω iω j ω k Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
13 Semidiscrete DG Formulation weak formulation Ju t, φ E = f (u, u, n), φ E + F, φ E J ijk (û ijk ) t = + + [ [ [ f,ξ1 jk f,ξ2 ik f,ξ3 ij ˆl i ( 1) f,ξ1 jk ˆl j ( 1) f,ξ2 ik ˆl k ( 1) f,ξ3 ij ˆl i (1) + ˆl j (1) + ˆl k (1) + only 1D operations due to tensor product structure Gauss integration weigths are hidden in ˆ -terms explicit Runge Kutta time integration ] N ˆD iλ ˆF λjk 1 λ=0 N λ=0 N λ=0 ˆD jλ ˆF 2 iλk ] ˆD kλ ˆF 3 ijλ ] Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
14 Algorithm u Extrapolate u u 1. MPI 2. other Send S M : u explicit Runge Kutta time integration 1. MPI Lifting-Flux 2. other Lifting Volume Operator Lifting Surface Integral Extrapolate u u Volume Operator DG Fluxes 1. MPI 2. other 1. MPI 2. other Receive M S : u + Send M S : Lifting-Flux Receive S M : Lifting-Flux Send S M : u Receive M S : u + Send M S : Flux Surface integral 1. other 2. MPI Receive S M : Flux u t Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
15 Order Of Convergence (DG) poly. deg. cells L 2 error L 2 order L error L order theor. order e e e e N= e e e e e e e e e e N= e e e e e e e e e e N= e e e e e e e e e e N= e e e e e e Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
16 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
17 Finite Volume Sub-cells DG u i,1 ξ ξ 1 u 1,j û ij = interface flux points (Riemann solver) = interpolation/integration points Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
18 Finite Volume Sub-cells DG u i,1 ξ 2 1 FV sub-cells 0 1 ξ 1 u 1,j û ij = interface flux points (Riemann solver) = interpolation/integration points DG reference element is splitted into (N + 1) d sub-cells = FV sub-cell borders = inner flux points (Riemann solver) Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
19 Finite Volume Method integration over sub-cell Ju t dξ + κ ijk κ ijk F dξ = 0 i, j, k = 0,..., N divergence theorem + numerical flux function Ju t dξ + f (u, u, n) ds ξ = 0 κ ijk κ ijk choose nodal value û ijk as constant value of ijk-th FV sub-cell J ijk (û ijk ) t = ω [ jω k ω i ω j ω k f,ξ1 i 1 2,j,k f,ξ1 i+ 1 2,j,k ] + ω [ ] iω k f,ξ2,ξ2 f ω i ω j ω i,j 1 k 2,k i,j+ 1 2,k + ω [ ] iω j f,ξ3 f,ξ3 ω i ω j ω i,j,k 1 k 2 i,j,k+ 1 2 Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
20 Coupling of DG and FV P interface fluxes computed on FV points solution is interpolated to FV fluxes projected back to DG conservative / free-stream preserving even on curved meshes coupling of DG and FV via numerical fluxes Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
21 2nd Order Reconstruction first order FV is very dissipative improved by 2nd order reconstruction total variation diminishing (TVD) limiters: minmod, van Leer, Sweby,... computation of slopes: Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
22 Algorithm u Extrapolate u u 1. MPI 2. other Send S M : u, Type Receive M S : u +, Type explicit Runge Kutta time integration Lifting-Flux 1. MPI 2. other Lifting Vol. Op. / Inner Reconstr. Lifting Surf Int. / Surf. Reconstr. Extrapolate u u 1. MPI 2. other Vol. Op. DG / FV inner faces Fluxes 1. MPI 2. other Send M S : Lifting-Flux Receive S M : Lifting-Flux Send S M : u Receive M S : u + Send M S : Flux Surface integral 1. other 2. MPI Receive S M : Flux u t Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
23 Indicators detect shocks Persson indicator Jameson indicator (3D sub-cell modification) Ducros sensor (with Jameson to avoid turbulence detection) upper and lower threshold indicator value upper threshold lower threshold time DG FV DG Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
24 Timestep Restriction CFL condtion DG timestep: where: CFL = t λ x t = α RK (N) x DG 2N + 1 λ αrk - constant of explicit time integration λ - max propagation velocity FV timestep of sub-cell: with 2 N+1 = size of FV sub-cell 2 t = α RK (0) N + 1 x DG λ, Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
25 Order Of Convergence (FV sub-cells) no reconstruction minmod central poly. deg. cells L 2 error L 2 order L 2 error L 2 order L 2 error L 2 order e e e e e e N= e e e e e e e e e e e e e e e N= e e e e e e e e e e e e e e e N= e e e e e e e e e Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
26 Remarks same DOFs, same data structure straightforward implementation comparable performance (DG/FV sub-cells per element) conservative interpolation/projection between DG and FV sub-cells 2nd order reconstruction for FV sub-cells slope limiters, TVD Indicators decide where to use FV sub-cells, where DG Persson indicator JST scheme like indicator Ducros sensort... no additional timestep restriction Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
27 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
28 Mesh refinement with Mortar interfaces complex geometries mesh refinement especially for hex-only-meshes non-conforming interfaces Mortar approach 3 different types in 3D Type 2 Type 1 Type 3 Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
29 Mesh refinement with Mortar interfaces Type 1 is combination of types 2 and 3 Interpolation-operators M 0 1 and M 0 2: M 0 1 : big left half Projection-operators M 1 0 and M 2 0: M 0 2 : big right half M 1 0 : left half big M 2 0 : right half big operators reduce to matrix-vector-product use this operators in η or ξ direction additional operators M * * FV for FV sub-cell shock capturing Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
30 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
31 Parallel performance (strong scaling) speedup 6 elements per core 10 4 pure DG pure FV 87 95% half/half checker board 10 3 ideal cores polynomial degree N = DOFs Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
32 Performance index (PID) PID = wall-clock-time#cores #DOF#timesteps 24 to cores 3072 to 6 elements per core 8 7 PID (µs/dof) Load (DOF/core) pure DG pure FV half/half checker board Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
33 Parallel performance (strong scaling) speedup elements per core PID (µs/dof) cores Load (DOF/core) realistic test case (N=8, cells) strong scaling down to about 6 elements per core cache effects best performance Hazel Hen of HLRS Stuttgart (Cray XC40) Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
34 Forward Facing Step Mach 3 flow over step in wind tunnel Mortar interfaces, refinement at 1. step 6543 grid cells 2. additionally: top shear wave (red region) 8343 grid cells polynomial degree N = DOF computation by M. Hoffmann Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
35 Forward Facing Step Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
36 Forward Facing Step Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
37 Shockless explosion combustion test case validation test case only grid cells N = DOFs curved mesh mortar interfaces computation by N. Krais Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
38 Shockless explosion combustion test case Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
39 Towards Air-Breathing Hypersonic Vehicles (M. Atak) How effcient are DG methods for such class of problems? This will be the first DNS of a compressible supersonic turbulent boundary layer with DG Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
40 Shock Boundary Layer Interaction (M. Atak) DNS of a shock wave interacting with a flat plate turbulent boundary layer at M = 2.67 shock intensity p 2 /p 1 = 1.3 yields a steady state structured hexahedra mesh with 1,125,000 grid cells N = DOF indicator: Ducros with Jameson up to 56,250 processors more than 10 6 CPU hours no mortar interfaces computation & visualization by M. Atak Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
41 Shock Boundary Layer Interaction (M. Atak) Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
42 Content 1 Discontinuous Galerkin Spectral Element Method 2 Finite Volume Sub-cell Shock Capturing 3 Mesh refinement with Mortar interfaces 4 Numerical Results 5 Conclusion & Outlook Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
43 Conclusion & Outlook Conclusion Discontinuous Galkerkin Spectral Element Method Shock Capturing with Finite Volume Sub-cells Mortar interfaces for mesh refinement Numerical examples Outlook adaptiv mesh refinement (AMR) dynamic load balancing References M. Sonntag, C.-D. Munz: Efficient Parallelization of a Shock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells, Journal of Scientific Computing, in Review. Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
44 Thank you for your attention! Visit flexi-project.org nrg.iag.uni-stuttgart.de
45 Sod Shock Tube Initial conditions: Mesh: 10 DG elements N = 11 Jameson indicator I Jameson = left right ρ v 0 0 p N i,j,k=0 p min,ijk 2p ijk + p max,ijk ω i ω j ω k p min,ijk + 2p ijk + p max,ijk Vol(Element), where p min,ijk is the minimal pressure of the surrounding nodes t end = 0.2 Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
46 Sod Shock Tube 1 Exact solution DG part of solution FV part of solution Density x Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
47 Strong Shock-Vortex Interaction Mach Shock: 7.0 Mach Vortex: 1.7 Grid: h = 1/100 N = 5 FV distribution: equidistant Indicator: Jameson (with Persson on FV DG) Riemann: LF (HLLE at BC) Limiter: MinMod CFL: 0.5 Initial condition: numerical converged stationary 1D-shock Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
48 Strong Shock-Vortex Interaction Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
49 Strong Shock-Vortex Interaction Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
50 Double Mach Reflection send shock wave diagonally into a wall shock speed: Mach 10 Grid: h = 1/120 (480 x 120 cells) polynomial degree N = DOF Indicator: Jameson (with Persson on FV DG) Riemann: HLLE Limiter: MinMod no mortar interfaces Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
51 Double Mach Reflection Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
52 Double Mach Reflection Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
53 Double Mach Reflection 0.59 ρ Matthias Sonntag (IAG) Shock Capturing for DG using FV Sub-cells SHARK-FV May 23, / 39
A 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 informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
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 informationHybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics
Hybrid DG/FV schemes for magnetohydrodynamics and relativistic hydrodynamics Jonatan Núñez-De La Rosa, Claus-Dieter Munz To cite this version: Jonatan Núñez-De La Rosa, Claus-Dieter Munz. Hybrid DG/FV
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 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 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 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 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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
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 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 informationarxiv: v1 [math.na] 23 Jun 2017
On the use of kinetic energy preserving DG-schemes for large eddy simulation David Flad 1 arxiv:1706.07601v1 [math.na] 23 Jun 2017 Institute for Aerodynamics and Gas Dynamics, University of Stuttgart,
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
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 informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 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 Dept. of
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More 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 informationA Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method. Joshua Favors
A Comparison of Artificial Viscosity Sensors for the Discontinuous Galerkin Method by Joshua Favors A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
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 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 informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationarxiv: v1 [physics.comp-ph] 21 Nov 2018
Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows Ray Vandenhoeck a,, Andrea Lani b a Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300, 3001, Leuven,
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
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 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 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 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 informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
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 informationZFS - Flow Solver AIA NVIDIA Workshop
ZFS - Flow Solver AIA NVIDIA Workshop Andreas Lintermann, Matthias Meinke, Wolfgang Schröder Institute of Aerodynamics and Chair of Fluid Mechanics RWTH Aachen University Outline Softwareproject ZFS (applications
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationRecent Developments in the Flux Reconstruction Method
Recent Developments in the Flux Reconstruction Method Joshua Romero, PhD candidate Kartikey Asthana, PhD candidate Jonathan Bull, Postdoctoral Scholar Antony Jameson, Professor Aerospace Computing Laboratory,
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 informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationHigh Order Discontinuous Galerkin Methods for Aerodynamics
High Order Discontinuous Galerkin Methods for Aerodynamics Per-Olof Persson Massachusetts Institute of Technology Collaborators: J. Peraire, J. Bonet, N. C. Nguyen, A. Mohnot Symposium on Recent Developments
More informationPreconditioning for modal discontinuous Galerkin methods for unsteady 3D Navier-Stokes equations
Preconditioning for modal discontinuous Galerkin methods for unsteady 3D Navier-Stokes equations Philipp Birken a Gregor Gassner b Mark Haas b Claus-Dieter Munz b a University of Kassel, Department of
More informationHigh-Order Methods for Diffusion Equation with Energy Stable Flux Reconstruction Scheme
49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposi 4-7 January 11, Orlando, Florida AIAA 11-46 High-Order Methods for Diffusion Equation with Energy Stable Flux
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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationIncompressible Flow Solver by Means of Pseudo-Compressibility Method
Incompressible Flow Solver by Means of Pseudo-Compressibility Method Árpád Veress Department of Aircraft and Ships Budapest University of Technology and Economics H Budapest, Sztoczek u. 6. J. ép., Hungary
More informationInterior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations
Interior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations Will Pazner Brown University, 8 George St., Providence, RI, 9, U.S.A. Per-Olof Persson University of
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
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 informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
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 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 informationUtilising high-order direct numerical simulation for transient aeronautics problems
Utilising high-order direct numerical simulation for transient aeronautics problems D. Moxey, J.-E. Lombard, J. Peiró, S. Sherwin Department of Aeronautics, Imperial College London! WCCM 2014, Barcelona,
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 informationSpectral Difference Method for Unstructured Grids II: Extension to the Euler Equations
Journal of Scientific Computing, Vol. 32, No. 1, July 2007 ( 2006) DOI: 10.1007/s10915-006-9113-9 Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations Z. J. Wang, 1,4
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationarxiv: v3 [math.na] 10 Mar 2015
A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws Michael Dumbser a, Olindo Zanotti a, Raphaël Loubère b, Steven Diot c a Department of
More informationImplicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference Method
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4442 Implicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference
More informationThe Discontinuous Galerkin Method on Cartesian Grids with Embedded Geometries: Spectrum Analysis and Implementation for Euler Equations
The Discontinuous Galerkin Method on Cartesian Grids with Embedded Geometries: Spectrum Analysis and Implementation for Euler Equations by Ruibin Qin A thesis presented to the University of Waterloo in
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 High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 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 informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
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 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 informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationTAU Solver Improvement [Implicit methods]
TAU Solver Improvement [Implicit methods] Richard Dwight Megadesign 23-24 May 2007 Folie 1 > Vortrag > Autor Outline Motivation (convergence acceleration to steady state, fast unsteady) Implicit methods
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 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 informationProspects for High-Speed Flow Simulations
Prospects for High-Speed Flow Simulations Graham V. Candler Aerospace Engineering & Mechanics University of Minnesota Support from AFOSR and ASDR&E Future Directions in CFD Research: A Modeling & Simulation
More informationA novel discontinuous Galerkin method using the principle of discrete least squares
A novel discontinuous Galerkin method using the principle of discrete least squares Jan Glaubitz and Philipp Öffner October 4, 2017 In this work, a novel discontinuous Galerkin DG) method is introduced
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
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 informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationControl Theory Approach to Aero Shape Optimization
Control Theory Approach to Aero Shape Optimization. Overview. Finite Volume Method 3. Formulation of Optimal Design Problem 4. Continuous Adjoint Approach 5. Discrete Adjoint Approach Overview of Adjoint
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
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 informationImproving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on Matthew R. Norman Scientific Computing Group National Center for Computational Sciences Oak Ridge
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
More informationEmerging Numerical Methods for Atmospheric Modeling
Chapter 9 Emerging Numerical Methods for Atmospheric Modeling Ramachandran D. Nair, Michael N. Levy and Peter H. Lauritzen Abstract This chapter discusses the development of discontinuous Galerkin (DG)
More information