AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
|
|
- Mark Miles
- 5 years ago
- Views:
Transcription
1 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 International Conference on Advances in Scientific Computing in Memory of Professor David Gottlieb December 6-8, 2009, Brown University, Providence, RI
2 1 Introduction High Order Methods SD Applications SD Stability and Motivation 2 DG Stability Proof DG and Nodal DG Formulation Nodal DG Energy Stability 3 SD Stability Proof SD Formulation SD Reformulation Energy Stability of SD Scheme 4 Conclusion 5 Closing Remarks Dedication Acknowledgement
3 INTRODUCTION
4 High Order Methods Overview of Traditional Methods Traditional Methods 1 Structured Finite Difference and Finite Volume Methods simple formulation not geometrically flexible 2 Unstructured finite difference and finite volume methods inefficient but geometrically flexible limited to 2nd order accuracy
5 High Order Methods Overview of High Order Methods Modern High Order Methods 1 Discontinuous Galerkin Method geometrically flexible and locally conservative finite element weak formulations with coupled equations 2 Spectral Volume Method geometrically flexible and locally conservative finite volume direct formulations with uncoupled equations 3 Spectral Difference Method geometrically flexible and locally conservative finite difference direct formulations with uncoupled equations
6 SD Applications Various SD Applications Turbulence Transition Shock Capture Moving Deforming Mesh Adaptive Mesh Refinement
7 SD Applications Shock Capturing Use of artificial bulk viscosity with a dilatation sensor Artificial viscosity is switched on only in regions of strong negative dilatation (shocks) Smooth variation of artificial viscosity Use of adaptive mesh refinement in combination with artificial viscosity
8 SD Applications Turbulence Transition Test case: Transitional Flow over SD7003 airfoil at 4 degrees AOA and Re=60000 Code: 3D, parallel, unstructured solver for Navier-Stokes equations, NO sub-grid models Results: Good agreement with computational and experimental results Data Set Freestream Separation Transition Reattachment Turbulence x sep x tr /c x r /c TU-BS (Experiment) 0.08% HFWT (Experiment) 0.1% Visbal (ILES) Uranga (ILES, DG) Present ILES
9 Introduction DG Stability Proof SD Stability Proof Conclusion Closing Remarks SD Applications Euler Vortex Propagation and Deformable Mesh Blending vor: E E E E E E E-02 vor: E E E E E E E E E-02 High Order Temporal and Spatial Accuracy Demonstrated Mesh Blending Technique with High Order Polynomial Numerical Results On Deformable Meshes Agree Well with Experiments Further Extension to FSI Problem Underway Plunging Airfoil and FSI on Deforming Mesh
10 SD Applications Adaptive Mesh Refinement Test Case: NACA0012 Airfoil, M=0.5, Steady Inviscid Flow Error Indicator 1: Unweighted Residual η k = R(Q H h ) Error Indicator 2: Entropy Adjoint η k = (V H h )T R(Q H h ) Two Mesh Adaptation Examples Illustrated
11 SD Stability and Motivation Spectral Difference Stability SD Stability Unstable for higher order triangular elements Sometimes weakly unstable in 1D depending on choices of flux collocation points
12 SD Stability and Motivation Motivation The aim is to prove: SD Stability in 1D The SD method for the 1D linear advection equation is stable for all orders of accuracy in an energy norm of Sobolev type. The norm is: Sobolev Norm u = (u 2 + cu (p)2 )dx where the coefficient c must be determined to cause a cancellation of terms (shown later).
13 SD Stability and Motivation STABILITY OF THE NODAL DISCONTINUOUS GALERKIN METHOD
14 Scalar Conservation Law Linear Advection Equation u t + u f (u) = x t + a u x = 0 we assume a > 0forrightrunningwave.
15 DG and Nodal DG Formulation Element Basis Function DG with Modal Basis NDG with Nodal Basis u h = nx u j Φ j j=1 u h = nx u j l j (x) j=1 where Φ is the modal basis where l is the Lagrange polynomial of degree p = n 1, and u j is the solution value at x j.
16 DG and Nodal DG Formulation Orthogonal Test Function The Galerkin formulations choose the same basis function to be the test function. The residual R h = u h t + x f (u h) is required to be orthogonal to the test function: Z xr x L R h Φ j dx = 0 or Z xr x L R h l j dx = 0
17 DG and Nodal DG Formulation Weak and Strong Forms Applying the technique of Integration-By-Parts to obtain the weak and strong forms of DG formulations: Weak Form Z xr x L uk h t Φ j dx Z xr x L f (uh k ) Φ j x dx + ˆf x R Φ j = 0 x L is required to be orthogonal to the test function: Strong Form Z xr x L uk h t Φ j dx + Φ j x L x f (uk h )dx +(ˆf f (uh k ))Φ x R j = 0 x L Z xr and similarly for Nodal DG formulation.
18 DG and Nodal DG Formulation Matrix Representation for DG Local Vectors u T =[u 1,...u n ] f = au Φ T =[Φ 1,...Φ n ] Mass and Stiffness Matrices M ij = S ij = Z xr x L Z xr x L Φ i Φ j dx Φ i Φ j dx Matrix Representation of the Weak Form M du dt S T f + ˆf x R Φ = 0 x L Matrix Representation of the Strong Form M du dt x R + Sf +(ˆf f )Φ = 0 x L
19 DG and Nodal DG Formulation Matrix Representations for Nodal DG Scheme Local Vectors u T =[u 1,...u n ] f = au l T =[l 1,...l n ] Mass and Stiffness Matrices M ij = S ij = Z 1 1 Z 1 1 l i l j dx l i l j dx Matrix Representation of the Weak Form M du dt S T f + ˆf l 1 1 = 0 Matrix Representation of the Strong Form M du dt + Sf +(ˆf f )l 1 1 = 0
20 DG and Nodal DG Formulation Matrix Representations for Nodal DG Scheme Multiplying by M 1,thestrongformofNodalDGformulation can be expressed as Matrix Representation of the Strong Form du dt + Df + M 1 l(ˆf f ) 1 1 = 0 where Differentiation Matrix D = M 1 S
21 Nodal DG Energy Stability Linear Advection Energy Estimate Multiplying the linear advection equation by u and integrating over x, Z b a u u t dx = a Z b a u u Z b dx = a x a ( u2 2 ) x dx Thus it satisfies the energy estimate d dt Z b a u 2 2 dx = 1 2 a(u2 a u 2 b )
22 Nodal DG Energy Stability Energy Estimate for the Discrete Solution Multiply the strong form by the local solution u T to obtain, u T M du dt + u T Sf + u T x R l(ˆf f ) = 0 x L Since M and S are pre-integrated exactly, this is equivalent to d dt Z xr x L Z u2 xr h 2 dx + a x L u u h h x dx + u x R h(ˆf au h ) = 0 x L Then integrate the middle term and combine it with the last term, d dt Z xr x L u2 h 2 dx = (u hˆf a u2 h 2 ) x R x L
23 Nodal DG Energy Stability Numerical Flux Element Interface Let u L and u R be values of u h on the left and right sides of a cell interface. For the numerical flux we now take ˆf = 1 2 a(u R + u L ) 1 2 α a (u R u L ), 0 α 1 where if α = 0wehaveacentralflux,andifα = 1wehavetheupwindflux.
24 Nodal DG Energy Stability Energy Estimate for the Discrete Solution On summing over all the elements: the L.H.S. of the energy estimate for the discrete DG solution derived previously yields: d dt Z xr x L u2 h 2 dx = d dt Z b a u 2 h 2 dx the R.H.S u 2 (u hˆf a h 2 ) has negative contribution at every element boundary, except the inflow boundary, which is strictly less than the boundary condition in the true solution, as shown in the next slide. x R x L
25 Nodal DG Energy Stability Energy Estimate for the Discrete Solution In the interior interface collecting the contributions elements on the left and right sides, there is a total negative contribution of u Rˆf a u 2 R 2 (u Lˆf a u2 L 2 ) = 1 2 a(u2 R u2 L ) 1 2 α a (u R u L ) a(u2 R u2 L ) = 1 2 α a (u R u L ) 2 = negative contribution at every interior interface At the two end boundaries, Extrapolated upwind flux f b = au h at outflow boundary = negative contribution Use the true flux f a = au a at inflow boundary = positive but less than the true inflow boundary contribution au a u h 1 2 au2 h = 1 2 au2 a 1 2 a(u a u h ) 2 < 1 2 au2 a
26 Nodal DG Energy Stability Energy Stability of DG Energy Stability of DG This completes the proof that the DG scheme is energy stable for the linear advection equation.
27 Nodal DG Energy Stability STABILITY OF THE SPECTRAL DIFFERENCE METHOD
28 SD Formulation SD Formulation - Solution Represent the SD scheme in terms of a reference element covering [-1,1]. Solution Points Distribution The discrete solution is locally represented by Lagrange interpolation on the solution collocation points x j as Solution Reconstruction u h (x) = nx u j l j (x) j=1 where for polynomials of degree p, n = p + 1.
29 SD Formulation SD Formulation - Flux At interior flux points f j = f (u h (ˆx j )); u h (ˆx j ) is interpolated from u h (x). Single valued numerical flux ˆf is used at boundaries f ( 1) and f (1). Flux Points Distribution Flux is represented by a polynomial of degree p + 1, Flux Reconstruction f h (x) = n+1 X j=1 f jˆlj (x) where ˆl j (x) are the Lagrange polynomials defined by the n + 1fluxcollocationpoints ˆx j,whichincludetheelementboundaries.
30 SD Formulation SD Formulation - Residual Flux Divergence Differentiate the flux polynomial at the solution collocation points to obtain Compute Residual du i dt = n+1 X j=1 f jˆl j (x i )
31 SD Reformulation Recast SD Formulation with Explicit Boundary Flux Interior Flux and Boundary Flux Correction f = au ˆf ( 1) =auh ( 1)+f CL, ˆf (1) =au h (1)+f CR f CL = ˆf ( 1) au h ( 1), f CR = ˆf (1) au h (1) Upon substitution, and noting u h (x) is a polynomial of degree p hence exactly represented by the sum, the flux becomes, Flux Reconstruction (Follow the Procedure Proposed by Huynh) f h (x) = n+1 X j=1 f jˆlj (x) = f h (x) =f CLˆl1 (x)+f CRˆln+1 (x)+a n+1 X j=1 = f h (x) =f CLˆl1 (x)+f CRˆln+1 (x)+au h (x) u h (ˆx j )ˆl j (x)
32 SD Reformulation Recast SD Scheme Rewrite the SD scheme as u h t = f h x = u h t = a u h x f CLˆl 1 f CRˆl n+1 Evaluating this at the solution points du i dt = a = a nx D ij u j f CLˆl 1 (x i ) f CRˆl n+1 (x i ) j=1 nx D ij u j f CLˆl 1 (x i ) for upwind numerical flux j=1 where D = M 1 S is the differentiation matrix associated with the solution collocation points, and is uniquely determined by the points location and thepolynomialdegreep.
33 SD Reformulation Matrix Form of SD Scheme The SD scheme can be converted to a form which resembles the nodal DG method by multiplying it by the mass matrix to produce SD Scheme in NDG Form X j M ij du j dt + a X j X S ij u j = f CL M ijˆl 1 (x j ) j Now since ˆl 1 ( 1) =1andˆl 1 (1) =0, Boundary Flux nx M ijˆl 1 (x j )= j=1 Z 1 1 = ˆl 1 l i 1 l i (x) 1 nx ˆl 1 (x j )l j (x)dx = j=1 Z 1 1 Z 1 1 l i (x)ˆl 1 (x)dx Z 1 l i (x)ˆl 1 (x)dx = l i ( 1) l i (x)ˆl 1 (x)dx 1
34 SD Reformulation SD Vs Nodal DG Upon substitution the resulting SD scheme becomes, SD Scheme X j M ij du j dt + a X j S ij u j = f CL l i ( 1)+ Z 1 1 «l i (x)ˆl 1 (x)dx which differs from the corresponding nodal DG equation, NDG Scheme X j M ij du j dt + a X j «S ij u j = f CL l i ( 1) only in the last term.
35 SD Reformulation SD Reformulation In order to make the last term in the SD scheme equal to the nodal DG scheme, we consider the following strategy by replacing the mass matrix M by a matrix Q > 0suchthatitwill eventually cause a cancellation of the desired terms. Choice of New Mass Matrix Q = M + C QD =(M + C)D = MD + CD = S = CD = 0 Thus each row of C must be orthogonal to every column of D.
36 SD Reformulation Choice of Matrix Q In order to find a row vector which is orthogonal to every column of D, consider p th difference operator d T nx j=1 d j R p (x j )=R (p) p for any polynomial of degree p. Then further operation by the differentiation matrix D gives X X d i i j D ij R p (x j )=R (p+1) p = 0 Thus we find the desired matrix C that is orthogonal to D where c is an arbitrary parameter. Q = M + C = M + cdd T, with C = cdd T
37 SD Reformulation SD Scheme using Mass Matrix Q SD Scheme using Mass Matrix Q X j Q ij du j dt + a X j X X S ij u j = f CL M ijˆl 1 (x j ) f CL C ijˆl 1 (x j ) j j Extra Term X X f CL C ijˆl 1 (x j )= c f CL d i j j (p+1) d jˆl 1 (x j )= c f CL ˆl l (p) 1 i SD Scheme using Mass Matrix Q X j Q ij du j dt + a X j S ij u j = f CL l i ( 1)+ Z 1 1 «l i (x)ˆl (p+1) 1 (x)dx c ˆl 1 l (p) i
38 SD Reformulation Choice of Cancellation Constant Now if we can choose c so that the last two terms on the right cancel, we can attain an energy estimate with the norm u T Qu replacing u T Mu in each element. Required Cancellation f CL Z 1 1 l i (x)ˆl 1 (x)dx = f CL c ˆl (p+1) 1 l (p) i
39 SD Reformulation Choice of Cancellation Constant Choose the interior flux collocation points as the zeros of the Legendre polynomial L p (x) of degree p. Then L.H.S Z 1 1 ˆl1 (x)l i (x)dx =( 1)p Z 1 1 xl p (x)l i (x)dx, where ˆl 1 (x) =( 1) p 1 2 (1 x)l p(x) Only the leading term in xl i (x) contributes to the integral. Let L.H.S L p (x) =c p x p +... = (2p 1) x p +... p! and xl i (x) =pa i x p +... where a i is the leading coefficient in l i (x).
40 SD Reformulation Choice of Cancellation Constant The left hand side of the extra term becomes: L.H.S Z 1 1 xl p (x)l i (x)dx = 2p 2p + 1 a i c p The right hand side of the extra term can be expressed as R.H.S and since l i (x) =a i x p +... ˆl (p+1) =( 1) p (p + 1)!c p l i (p) = p!a i
41 SD Reformulation Choice of Cancellation Constant Thus the desired cancellation is obtained by setting Cancellation Constant c = 2p 2p c 2 p 1 p!(p + 1)! > 0 In the case that the interface flux is not fully upwind, a similar calculation shows that the convection from the right boundary is correspondingly reduced, so that finally SD Scheme Cast in Nodal DG Form X j Q ij du j dt + a X j S ij u j = f CL l i ( 1) f CR l i (1)
42 Energy Stability of SD Scheme Energy Stability of the SD Scheme Since u h is a polynomial of degree p, i d iu i = u (p) h,andin each element X X u i Q ij u j = i j Z xr xl (uh 2 + cu(p)2 )dx h Now the same argument that was used to prove the energy stability of the nodal DG scheme establishes the energy stability of the SD scheme with the norm Z b a (uh 2 + cu(p)2 )dx h for the case of solution polynomials of degree p, provided that the interior flux collocation points are the zeros of L p (x).
43 Energy Stability of SD Scheme Behavior of the Cancellation Constant It can be seen that c decreases very rapidly with increasing p, as is illustrated by the following table of values of c. cvsp p c
44 Energy Stability of SD Scheme CONCLUSION
45 Conclusion The result is consistent with the conclusion of Van den Abeele, Lacor and Wang that the stability of the spectral difference method depends only on the location of the flux collocation points. While it establishes the stability of the SD scheme when the interior flux collocation points are the zeros of the Legendre polynomial L p (x), itdoesnotprecludethestabilityofthesd scheme with other choices of the flux collocation points, possibly in a different norm. However, extensive calculations for the second, third and fourth order cases (not included here) have indicated that the conditions for the cancellation of the last two terms of the boundary correction can only be satisfied by choosing the interior flux collocation points as the zeros of L p (x).
46 Conclusion As the order of accuracy is increased the norm in which the SD method has been shown to be stable asymptotically approaches the usual energy norm. In light of the theorem that all norms are equivalent in a finite dimensional space, the proof also establishes the stability of the SD schemes in any norm. Whether a similar proof of stability for all orders of accuracy can be established for the multidimensional case with either tensor product or simplex elements remains a subject for future research.
47 ACKNOWLEDGEMENT AND DEDICATION
48 Dedication Dedication Dedication This paper is dedicated to the memory of David Gottlieb, a man of extraordinary honor, integrity and generosity, whose seminal contributions throughout his career have laid the foundation for the emergence of practically useful higher order methods in scientific computing.
49 Acknowledgement Acknowledgement Acknowledgement The author is indebted to both the National Science Foundation and the Air Force Office of Scientific Research for their support of the research that led to this paper through grants and FA
High-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 informationOn the Non-linear Stability of Flux Reconstruction Schemes
DOI 10.1007/s10915-011-9490-6 TECHNICA NOTE On the Non-linear Stability of Flux econstruction Schemes A. Jameson P.E. Vincent P. Castonguay eceived: 9 December 010 / evised: 17 March 011 / Accepted: 14
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 New Class of High-Order Energy Stable Flux Reconstruction Schemes
DOI 10.1007/s10915-010-940-z A New Class of High-Order Energy Stable Flux Reconstruction Schemes P.E. Vincent P. Castonguay A. Jameson Received: May 010 / Revised: September 010 / Accepted: 5 September
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 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 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 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 informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
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 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 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 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 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 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 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 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 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 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 informationDiscontinuous Galerkin methods Lecture 2
y y RMMC 2008 Discontinuous Galerkin methods Lecture 2 1 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu y 1 0.75 0.5 0.25 0-0.25-0.5-0.75 y 0.75-0.0028-0.0072-0.0117 0.5-0.0162-0.0207-0.0252
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 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 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 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 informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
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 informationOn the Stability of the Flux Reconstruction Schemes on Quadrilateral Elements for the Linear Advection Equation
OI 10.1007/s10915-015-0102-8 On the Stability of the Flux Reconstruction Schemes on Quadrilateral Elements for the Linear Advection Equation Abhishe Sheshadri 1 Antony Jameson 1 Received: 25 March 2015
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 informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
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 informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
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 informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
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 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 informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
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 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 informationSPECTRAL METHODS ON ARBITRARY GRIDS. Mark H. Carpenter. Research Scientist. Aerodynamic and Acoustic Methods Branch. NASA Langley Research Center
SPECTRAL METHODS ON ARBITRARY GRIDS Mark H. Carpenter Research Scientist Aerodynamic and Acoustic Methods Branch NASA Langley Research Center Hampton, VA 368- David Gottlieb Division of Applied Mathematics
More informationOpportunities for efficient high-order methods based on the summation-by-parts property
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5- Opportunities for efficient high-order methods based on the summation-by-parts property Jason E. Hicken, a David
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 informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
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 informationHow Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm
AIAA 2 2673 How Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm Antony Jameson Stanford University Stanford, CA 9435 D. A. Caughey
More informationHp-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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
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 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 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 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 informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
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 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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving
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 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 informationFurther Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow
29th AIAA Applied Aerodynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3509 Further Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow Antony Jameson, John C. Vassberg,
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 informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
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 informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
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 informationContinuous adjoint based error estimation and r-refinement for the active-flux method
Continuous adjoint based error estimation and r-refinement for the active-flux method Kaihua Ding, Krzysztof J. Fidkowski and Philip L. Roe Department of Aerospace Engineering, University of Michigan,
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 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 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 information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
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 informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
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 informationJournal of Computational Physics
Journal of Computational Physics 230 (2011) 8134 8154 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Insights from von Neumann
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 informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 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 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 informationLecture - Implementation and numerical aspects of DG-FEM in 1D
Lecture - Implementation and numerical aspects of DG-FEM in 1D Ph.D. Course: An Introduction to DG-FEM for solving partial differential equations Allan P. Engsig-Karup Scientific Computing Section DTU
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 informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
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 informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
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 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 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 informationSTABLE SPECTRAL METHODS ON TETRAHEDRAL ELEMENTS
SIAM J. SCI. COMPUT. Vol. 21, No. 6, pp. 2352 2380 c 2000 Society for Industrial and Applied Mathematics STABLE SPECTRAL METHODS ON TETRAHEDRAL ELEMENTS J. S. HESTHAVEN AND C. H. TENG Abstract. A framework
More informationA discontinuous Galerkin spectral/hp method on hybrid grids
Applied Numerical Mathematics 33 (2000) 393 405 A discontinuous Galerkin spectral/hp method on hybrid grids R.M. Kirby, T.C. Warburton, I. Lomtev, G.E. Karniadakis Brown University, Center for Fluid Mechanics,
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
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 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 information