A Stable Spectral Difference Method for Triangles
|
|
- Gladys Boyd
- 5 years ago
- Views:
Transcription
1 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 Technical University, Austria AIAA Aerospace Sciences Meeting January 4, 2011 Orlando, Florida Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 1 / 33
2 Outline 1 Background and Motivation 2 Spectral Difference(SD) Method 3 SD Method with Raviart-Thomas Elements 4 Stability Analysis 5 Results Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 2 / 33
3 Background and Motivation Spectral Difference (SD) high-order method for hyperbolic PDEs A quadrature free (pre-integrated) nodal Discontinuous Galerkin scheme Simple in formulation and implementation Found linearly unstable for triangles for order of accuracy > 2 Found stable with flux interpolation on Raviart-Thomas elements Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 3 / 33
4 SD Method for Triangles Hyperbolic conservation equation u( x,t) t + f (u) = 0 Transformation from reference element Φ : ξ x with J = x/ ξ (0,1) x1 x2 Φ :(ξ, η) (x,y) η (0,0) ξ (1,0) y x x3 Hyperbolic equation in reference domain ( u(ξ,t) t + 1 J ξ J J 1 f ) (u) Define Solution collocation nodes - ˆξ j Flux collocation nodes - ξ j = 0 Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 4 / 33
5 SD Method for Triangles Approximation of solution u h (ξ) = N m j=1 u jl j (ξ) l j P m l j (ˆξ k ) = δ jk u j = u h (ˆξj ) no. of degrees of freedom = N m = (m+1)(m+2) 2 Approximation of flux f h (ξ) = N m+1 j=1 f jˆlj (ξ) ˆlj P m+1 ˆlj ( ξ k ) = δ jk f j = f ) h ( ξj f j = { J J 1 f ( ξj ) ξj ˆT f num ξj ˆT f num n = h standard numerical flux no. of degrees of freedom = 2N m+1 Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 5 / 33
6 SD Method for Triangles Final form of the Spectral Difference scheme du (i) j dt + 1 N m+1 J (i) k=1 ) ξˆlk (i) (ˆξj f k = 0, j = 1,..., N m Linearly unstable for m 2 for triangles [Van den Abeele et al., 2008 ] Note - Each of the flux vectors need not be in P m+1 for the div to be in P m Raviart Thomas space Smallest space having div in P m Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 6 / 33
7 SD Method with Raviart-Thomas Elements Define Solution points - ˆξ j Flux points - ξ j, Directions - s j Approximation of solution u h (ξ) = N m j=1 u jl j (ξ) l j (ˆξ k ) = δ jk u j = u h (ˆξj ) Approximation of flux function in Raviart-Thomas (RT ) space f h (ξ) = N RT m j=1 f jψ j (ξ) ψ j ( ξ k ) s k = δ jk f j = f ) h ( ξj s j { ) J J 1 f ( ξj s j ξj f j = ˆT h ξj ˆT h standard numerical flux Nm RT = (m + 1) (m + 3) Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 7 / 33
8 SD Method with Raviart-Thomas Elements For a degree m, the RT space is defined as RT m = [P m ] 2 + (x, y) T P m. For m = 1, the monomials which form a basis in the RT space ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 x y x 2 xy,,,,,,, x y yx y 2 Less number of flux degrees of freedom compared to standard SD 2N m+1 N RT m = m + 3 Flux nodes distribution : m + 1 nodes on each edge and N RT m 3(m + 1) in the interior Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 8 / 33
9 SD Method with Raviart-Thomas Elements Final form of the new Spectral Difference scheme du (i) j dt + 1 Nm rt J (i) k=1 f (i) k ( ξ ψ ) ) k (ˆξj = 0, j = 1,..., N m Linearly stable for m = 1, 2, 3 in a simplified stability analysis Numerical experiments prove the viability of the scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 9 / 33
10 Linear Stability Analysis Linear advection equation u( x, t) t + f (u) = 0, f(u) = (u c cosθ, u c sinθ), θ [0, π 2 ] Consider Cartesian mesh with each element formed by fusing two triangles j i SD formulation, using upwind fluxes t U ( (i,j) = ν AU (i,j) + BU (i 1,j) + CU (i,j 1)), Linear stability analysis (LSA) Fourier transformation : u û Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 10 / 33
11 Linear Stability Analysis SD discretization of the Fourier mode ûe i(kxx+kyy) dû dt = ν t Zû Z = ( A + Be iσ + Ce iκ) (σ, κ) = (k x h, k y h) Full stability = Stability of spacial discretization + time discretization Stability of spacial discretization eigensystem of Z Stable flux nodes Re(λ(Z)) 0 Optimal flux nodes M ax( λ(z) ) (Spectral Radius) is minimum Stability is independent of the position of solution nodes Stability is independent of the position of flux nodes on the edges Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 11 / 33
12 LSA - Spatial Discretization RT 1 1 interior flux point at centroid - stable RT 2 3 interior points each with two ortho directions form 6 flux nodes 3 interior points are varied as ξ i = ξ c + α(ξi e ξ c ), i = 1, 2, 3 α [0, 1] Stable 0.5 α 0.521, considering θ [0, π 2 ] Stable and optimal α = 0.5 [higher order quadrature nodes] Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 12 / 33
13 LSA - Spatial Discretization RT 3 6 interior points each with two ortho directions form 12 flux nodes 6 interior points are varied as ξ i = ξ c + α(ξ e i ξ c ), i = 1, 2, 3 α [0, 1] ξ i = ξ c + β(ξ e i ξ c ), i = 4, 5, 6 β [0, 1] Stable and optimal α = nodes] β = [higher order quadrature Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 13 / 33
14 LSA - Spatial Discretization Stability and optimality for RT " Stable Region ! Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 14 / 33
15 LSA - Spatial Discretization y y y Stable and optimal flux points x RT x RT x RT 3 Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 15 / 33
16 LSA - Full Discretization Full discretization û n+1 = Gû n L 2 stability ρ(g) 1 Get allowable CFL number (= c t h ) Max CFL number for Shu-RK3 time discretization RT 1 RT 2 RT 3 CFL number ! Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 16 / 33
17 Results - 2D Linear Advection Equation Linear advection equation u( x, t) t + f (u) = 0, f(u) = (cx u, c y u) T Convergence study 1 2 RT1 RT2 RT3 3 Log(L Error) Log(N) Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 17 / 33
18 Results - 2D Euler Equations Euler equations in 2D f(u) x + g(u) = 0 y ρ ρu u = ρu ρv f = ρu 2 + p ρuv g = E u (E + p) NACA mesh elements 2.5 ρv ρuv ρv 2 + p v (E + p), Y X Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 18 / 33
19 Results - 2D Euler Equations Relaxation Technique Backward Euler / Damped Newton ( I t dr(u n ) ) U n = tr(u n ), du t Newton iteration (Quadratic convergence) Preconditioning - Incomplete LU factorization Linear system solver - Restarted GMRES algorithm Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 19 / 33
20 Results - 2D Euler Equations Test case 1 - Free stream Mach number - 0.3, Angle of attack - 0 degree Figure: Mach contours - RT 2 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 20 / 33
21 Results - 2D Euler Equations Test case 1 - Convergence of residual Log(Res) RT 1 RT 2 RT Iterations Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 21 / 33
22 Results - 2D Euler Equations Test case 1 -Comparison of 4th order DG and RT 3 schemes Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 22 / 33
23 Results - 2D Euler Equations Test case 1 - Free stream Mach number - 0.3, Angle of attack - 0 degree Figure: Mach contours - RT 1 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 23 / 33
24 Results - 2D Euler Equations Test case 1 - Free stream Mach number - 0.3, Angle of attack - 0 degree Figure: Mach contours - RT 2 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 24 / 33
25 Results - 2D Euler Equations Test case 1 - Free stream Mach number - 0.3, Angle of attack - 0 degree Figure: Mach contours - RT 3 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 25 / 33
26 Results - 2D Euler Equations Test case 2 - Free stream Mach number - 0.4, Angle of attack - 5 degree Figure: Mach contours - RT 2 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 26 / 33
27 Results - 2D Euler Equations Test case 2 - Convergence of Residual Log(Res) RT 1 RT 2 RT Iterations Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 27 / 33
28 Results - 2D Euler Equations Test case 2 -Comparison of 4th order DG and RT 3 Schemes Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 28 / 33
29 Results - 2D Euler Equations Test case 2 - Free stream Mach number - 0.4, Angle of attack - 5 degree Figure: Mach contours - RT 1 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 29 / 33
30 Results - 2D Euler Equations Test case 2 - Free stream Mach number - 0.4, Angle of attack - 5 degree Figure: Mach contours - RT 2 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 30 / 33
31 Results - 2D Euler Equations Test case 2 - Free stream Mach number - 0.4, Angle of attack - 5 degree Figure: Mach contours - RT 3 scheme Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 31 / 33
32 Summary and Outlook Difference between SD and new SD New SD formulation is found to be linearly stable Numerical results show the viability Needs to be extended to solve NS equations, also to simulate transonic flows Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 32 / 33
33 Acknowledgement Financial support from the Deutsche Forschungsgemeinschaft (German Research Association) through grant GSC 111, and by the Air Force Office of Scientific Research, Air Force Materiel Command, USAF, under grant number FA , is gratefully acknowledged Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 33 / 33
34 LSA - Spatial Discretization Stability and optimality for RT 2 max ",#,$ Re(%) ! max ",# Re($) 1 x ! max, (Z) Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 34 / 33
35 Runge-Kutta Time stepping The solution at (n + 1)-th iteration, U n+1, is obtained from U n as For stability analysis w (0) = U n, k 1 w (k) = α kl w (l) + tβ kl R (l) k = 1,..., p, l=0 w (0) = Û n, U n+1 = w (p), k 1 w (k) = α kl w (l) + νβ kl Zw (l) k = 1,..., p, (1) l=0 Û n+1 = w (p). If G (k) is the amplification matrix in the k-th intermediate step, then one obtains k 1 G (0) = I, G (k) = (α kl I + νβ kl Z) G (l) k = 1,..., p. l=0 Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 35 / 33
36 Runge-Kutta Time Stepping Shu-RK3 1 α = , β = stage 4th order SSP α = , β = Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 36 / 33
37 Backward-Euler Time Stepping Backward Euler Taylor series for R U n+1 U n t = R(U n+1 ), Rewrite R(U n+1 ) = R(U n ) + dr(u n ) t +... dt dr(u n ) dt t = dr(u n ) du du dt t dr(u n ) du (U n+1 U n ). If (U n+1 U n ) is denoted as U n, the implicit scheme is ( I t dr(u n ) ) U n = tr(u n ), du Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 37 / 33
38 Stability Analysis Numerical schemes should posess non-linear stability properties Eg. Total Variation Diminishing (TVD) N T T V (ū n ) = ū n i+1 ū n i, i=1 T V ( ū n+1) T V (ū n ). TVD property convergence A conservative numerical scheme can be made to satisfy TVD by using limiters Linearly unstable limiter will act on smooth regions affect order of accuracy Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 38 / 33
39 Finding Transfer Matrix Solution interpolation Solution at flux nodes Dubiner basis u h (ξ) = χ (ξ) = Dubiner basis at flux nodes N m j=1 u j l j (ξ) u h (ˆξj ) = u j u h ( ξk ) = N m j=1 χ j l j (ξ) ) χ ( ξk = N m j=1 N m j=1 u j l j ( ξk ) χ j l j ( ξk ) ) χ j = χ (ˆξj Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 39 / 33
40 Finding Differentiation Matrix The monomials in the RT space at solution nodes φ n (ˆξj ) = Nm RT k=1 a n,k ψk (ˆξj ) a n,k = φ n ( ξk ) s k Its divergence ξ φ n (ˆξj ) = Nm RT k=1 a n,k ( ξ ψ k ) (ˆξj ) Aravind Balan (AICES, RWTH Aachen) New Spectral Difference Jan 4, Orlando 40 / 33
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
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 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 Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes
A Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes N. Chalmers and L. Krivodonova Abstract When the discontinuous Galerkin (DG) method is applied to hyperbolic problems in
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 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 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 informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
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 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 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 numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
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 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 informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
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 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 New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements
DOI 10.1007/s10915-011-9505-3 A New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements P. Castonguay P.E. Vincent A. Jameson Received: December 010 / Revised: 7 April
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 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 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 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 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 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 informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
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 informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationFinite 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 informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
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 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 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 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 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 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 informationRunge-Kutta Residual Distribution Schemes
Runge-Kutta Residual Distribution Schemes Andrzej Warzyński, Matthew E. Hubbard, Mario Ricchiuto RESEARCH REPORT N 837 September 23 Project-Teams BACCHUS ISSN 249-6399 ISRN INRIA/RR--837--FR+ENG Runge-Kutta
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
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 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 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 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 informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More informationActive Flutter Control using an Adjoint Method
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an
More informationFirst order BSSN formulation of Einstein s field equations
David Brown 1 Peter Diener 2 3 Jan Hesthaven 4 Frank Herrmann 3 Abdul Mroué 5 Olivier Sarbach 6 Erik Schnetter 7 Manuel Tiglio 3 Michael Wagman 4 1 North Carolina State University 2 Louisiana State University
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 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 informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
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 informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationAccuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization
Preprint accepted in Journal of Computational Physics. https://doi.org/0.06/j.jcp.07.04.075 Accuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization Hiroaki Nishikawa and Yi
More information1 PART1: Bratu problem
D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationFinite Element methods for hyperbolic systems
Finite Element methods for hyperbolic systems Eric Sonnendrücker Max-Planck-Institut für Plasmaphysik und Zentrum Mathematik, TU München Lecture notes Wintersemester 14-15 January 19, 15 Contents 1 Introduction
More informationPerformance tuning of Newton-GMRES methods for discontinuous Galerkin discretizations of the Navier-Stokes equations
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2685 Performance tuning of Newton-GMRES methods for discontinuous Galerkin
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
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 informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
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 informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
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 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 informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
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 informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
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 informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
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 informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
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 informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationarxiv: v1 [math.na] 25 Jan 2017
Stage-parallel fully implicit Runge-Kutta solvers for discontinuous Galerkin fluid simulations arxiv:1701.07181v1 [math.na] 25 Jan 2017 Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics,
More informationADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES
ADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES Niles A. Pierce and Michael B. Giles Applied & Computational Mathematics, California Institute of Technology Computing Laboratory,
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
More informationApproximate tensor-product preconditioners for very high order discontinuous Galerkin methods
Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics, Brown University, Providence, RI, 02912
More informationarxiv: v1 [math.na] 22 Nov 2018
Asymptotic preserving Deferred Correction Residual Distribution schemes Rémi Abgrall and Davide Torlo arxiv:1811.09284v1 [math.na] 22 Nov 2018 Abstract This work aims to extend the residual distribution
More informationModel reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature. Masayuki Yano
Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony Patera;
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationContents of lecture 2b. Lectures 2a & 2b. Physical vs. computational coordinates [2] Physical vs. computational coordinates [1]
Contents of lecture b Lectures a & b P. A. Zegeling Mathematical Institute Utrecht University The Netherlands Parameter-free non-singular moving grids in D: Theory & properties Application to resistive
More informationCurved grid generation and DG computation for the DLR-F11 high lift configuration
ECCOMAS Congress 2016, Crete Island, Greece, 5-10 June 2016 Curved grid generation and DG computation for the DLR-F11 high lift configuration Ralf Hartmann 1, Harlan McMorris 2, Tobias Leicht 1 1 DLR (German
More information