arxiv: v1 [math.na] 21 Nov 2017
|
|
- Clifton Miles
- 5 years ago
- Views:
Transcription
1 High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete approximation to a partial differential equation PDE is constructed a discretization of the spatial operator with a truncation error τ is derived. This discrete operator should be semi-bounded for the scheme to be stable. Under these conditions the Lax Ricchtmyer equivalence theorem assures that the scheme converges and that the error will be, at most, of the order of τ. In most cases the error is in indeed of the order of τ. We demonstrate that for the Heat equation stable schemes can be constructed, whose truncation errors are τ, however, the actual errors are much smaller. This gives more degrees of freedom in the design of schemes which can make them more efficient more accurate or compact than standard schemes. In some cases the accuracy of the schemes can be further enhanced using post-processing procedures. Introduction Consider the differential problem: u ut = 0 = fx. = P u, x Ω R d,t 0 x where u=u,...,u m T and P / x is a linear differential operator with appropriate boundary conditions. It is assumed that this problem is well posed, i.e. Kt< such that ut Kt f, where typically Kt=Ke αt. Let Q be the discretization of P / x where we assume: Assumption : Assumption : Assumption 3: The discrete operator Q is based on grid points{x j }, j =,...,N. Q is semi bounded in some equivalent scalar product, H =,H, i.e w,qw H αw,w H = α w H. The local truncation error vector of Q is T e which is defined, at each entry j by T e j = Pwx j Qw j, 3 where wx is a smooth function and w is the projection of wx onto the grid. It is assumed that T e N 0. Consider the semi discrete approximation: A. Ditkowski School of Mathematical Sciences, Tel Aviv University, Tel Aviv 69978, Israel adid@post.tau.ac.il. The final publication is available at link.springer.com via 3. In later works these schemes are called Error Inhibiting Schemes EIS.
2 v = Qv, t 0 4 vt = 0 = f. Proposition: Under Assumptions 3 The semi discrete approximation converges. Proof: Let u is the projection of ux, t onto the grid. Then, from assumption 3, u = Pu = Qu+T e. 5 Let E = u v then by subtracting 4 from 5 one obtains the equation for E, namely E = QE+T e. 6 Next, by taking the H scalar product with E, using assumption and the Schwartz inequality the following estimate can be derived E, E = E,E H = E H E H = E,QE H +E,T e H αe,e H + E H T e H. Thus Therefore: H E H α E H + T e H. 7 E H t E H t = 0e αt + eαt α max T e H 0 τ t N 0. 8 Here we assumed that E H t = 0 is either 0, or at least of the order of machine accuracy. Equation 8 establishes the fact that if the scheme is stable and consistent, the numerical solution v converges to the projection of the exact solution onto the grid, u. Furthermore, it assures that the error will be at most in the truncation error T e H. This is one part of the landmark Lax Richtmyer equivalence theorem for semi discrete approximation. See e.g. [6]. Due to this and similar results the common way for constructing finite difference schemes is to derive a semi bounded Q with proper truncation error. Typically the error E H is indeed of the order of T e H. It should be noted, however, that 6 is the exact equation for the error dynamics, while 8 is an estimate. In this paper we present finite difference schemes in which the errors are smaller than their truncation errors. It is well known that boundary conditions can be of one order lower accuracy without destroying the convergence rate expected from the approximation at inner points, see e.g. [4], [5], [] and [0]. In [5] and [0] it was shown that for parabolic, incompletely parabolic and nd-order hyperbolic equations the boundary conditions can be of two order less. Here, however, we consider low order truncation errors in most or all of the grid points. This paper is constructed as follows; in Section we present a preliminary example and illustrate the mechanism which reduces the error. In Section 3 we present a two point block scheme which has a first order truncation error but has a second or third order error. In Section 4 two three point block schemes are presented. Discussion and remarks are presented in Section 5. Preliminary Example Consider the heat equation u = u + Fx,t, x [0,π,t 0 x ut = 0 = fx 9
3 with periodic boundary conditions. Let the scheme be: v j = D + D v j + j cv j + F j t ; x j = jh, h=π/n, N is even v j t = 0 = f j 0 where F j t is the projection of Fx,t onto the spatial grid. The truncation error is T e j = h u j xxxx + Oh j cv j = O. Formally, this scheme is not consistent. The scheme 0 was run with the initial condition fx = cosx, F = 0 and N = 3,64,...,04 with forward Euler time propagator. The plots of log 0 E vs. log 0 N, at t final = π, for c=0,0.5, are presented in Figure a. It can be seen that this is a second order scheme. a c=0, slope:.0069 c=0.5, slope:.005 c=.0, slope:.009 b c=0, slope:.0040 c=/6, slope: c= /6, slope:.895 c= /4, slope: Fig. Convergence plots, log 0 E vs. log 0 N, for different values of c. a: Scheme 0. b: Scheme 6. In order to understand this phenomenon let us consider one high frequency mode of the error 3. Denote T c j = j cv j. This term can also be written as: T c j = j ch α = ch α e inx j/ For the scheme 0 α = 0. The equation for the error term caused by T c is E c = D + D E c + T c 3 SinceE c j = Ê c e inx j/ π, Ê c C, then the equation for Ê c is 3 This scheme was presented for demonstrating the phenomenon that the error, due to high frequency modes, is lower than the truncation error. As this is not a practical scheme, full analysis of the error is not presented, only a demonstration of the dynamics of high frequency error modes is presented. Full analysis is given for the scheme presented in the next section 3
4 Ê c = N Êc + c h α 4 Therefore, E c t = Êc t = Êc 0e N t + Ê c 0e N t + O h α+ e N t c h α 5 N Note that the actual error, E c t, is two orders lower that the truncation error, T c. In the next sections we present practical schemes which utilize this idea. 3 Two point Block, 3rd Order Scheme Let the grid be: x j = j h, h=π/n+ and x j+/ = x j + h/, j= 0,...N. Altogether there are N+ points with spacing of h/. For simplicity, we assume that N is even. Consider the approximation: The truncation errors are: d dx u j h/ [ u j / u j + u j+/ + c u j / + 3u j 3u j+/ + u j+ d dx u j+/ h/ [ u j u j+/ + u j+ + c u j / 3u j + 3u j+/ u j+ T e j = T e j+ = h h u j xxxx + c [ h u j+ xxxx + c [ u j xxx + ] h u j xxxx + Oh 3 =Oh h u j+ xxx + h u j+ xxxx ] + Oh 3 =Oh 6 7 The motivation leading to this scheme is that the highly oscillating Oh error terms will be dissipated, as in the previous example, while the Oh terms will be canceled, for the proper value of c. 3. Analysis Let ω { N/,...,N/} and ω N+ ω > 0 ν = 8 ω+n+ ω 0 Then e iωx j = e iνx j and e iωx j+/ = e iνx j+/. 9 We look for eigenvectors in the form of: 4
5 ψ k ω= α k π. e iωx j e iωx j+/. + β k π. e iνx j e iνx j+/. 0 where, for normalization, it is require that α k + β k =, k =,. The expressions for α k, β k and the eigenvalues symbols ˆQ k are: α = [ + c cos4h/ω+4c cosh/ω+4c7c 8+cosh/ω+35c 3c+8 c sinh/ω+sinh/ω ] i8c 4cosh/ω+ β = csinh/ω+sinh/ωα [ c β = + sinh/ω+sinh/ω c cos4h/ω+4 c cosh/ω+4c7c 8+cosh/ω+35c 3c+8 ] 3 where and the α = icsinh/ω+sinh/ω 4 8ccosh/ω+ β = c cos4h/ω+38c + 8c 3c cosh/ω 3c ˆQ, ω= 4+ccosh/ω+3± h/. 6 It can be shown that the eigenvalues are real and non positive for c < /. It should be noted that by construction ψ k ω ψ k ν, ω ν. However, for all ω ψ ω and ψ ω are not orthogonal. Nevertheless it can be shown that they are almost perpendicular in the sense that Therefore the scheme is stable. For ωh the eigenvalues and eigenvectors are: cosθω = ψ ω,ψ ω ψ ω ψ ω > 0.9 ˆQ ω= ω + +4cω4 4c h + Oh 4 7 and c α = 3 c ωh 6 + O h 7, β = ic 4 8c ωh 3 + Oh 5 8 ˆQ ω= 4 8c h/ + 4cω + Oh 9 5
6 If the initial condition is α = ic c ωh + Oh 3, β = +Oh 30 v j 0=e iωx j,v j+ 0=e iωx j+ ; ω h 3 then v j t = e ω t +4cω t 4c ωh + Oh 4 e iωx j + ic 4 8c ωh 3 + Oh 5 e iνx j 3 The same expression hold for x j+. Therefore the scheme is, in general, nd order and it is 3rd order if c= /4. Note that by naive analysis of the truncation error terms, 7, one would expect to get 3rd order with c= /6. We used the approximation 6 for solving the heat equation 9, where Fx, t and the initial condition were chosen such that the exact solution is ux,t = expcosx t. The scheme was run with N = 3,64,...,04. 4th order Runge Kutta scheme was used for time integration. The plots of log 0 E vs. log 0 N for c=0,/6, /6, /4 are presented in Figure b. As can be seen, the results are as predicted by the analysis. 4 Three point Block In this section we briefly present two schemes which are based on three point block. Here we use the grid, x j = j h, h=π/n+ with the internal block nodes x j+/3 = x j + h/3 and x j+/3 = x j + h/3, j= 0,...N. Altogether there are 3N+ points with spacing of h/3. three point block, 3rd order scheme Consider the approximation: d dx u j = d dx u j+/3 = d dx u j+/3 = 4h/3 [ 4u j /3 8u j + 4u j+/3 + c u j /3 + 3u j 3u j+/3 + u j+/3 + Oh 4h/3 [ 4u j 8u j+/3 + 4u j+/3 + Oh 33 4h/3 [ 4u j+/3 8u j+/3 + 4u j+ + c u j 3u j+/3 + 3u j+/3 u j+ + Oh This scheme was run under the same conditions as the example in Section 3. As in the previous example the truncation error is of Oh, however this is a nd order scheme and 3rd order for c =.340. The convergence results are presented in Figure a. 6
7 a b c=0, slope: c=.0, slope: c=.0, slope: c= 0.385, slope: c=0, slope:.0009 c=.0, slope:.007 c=.0, slope:.597 c=.340, slope: Fig. Convergence plots, log 0 E vs. log 0 N, for different values of c. a: Scheme 33. b: Scheme 34. three point block, 5th order scheme d dx u j = d dx u j+/3 = d dx u j+/3 = h/3 [ u j /3 + 6u j /3 30u j + 6u j+/3 u j+/3 + c u j /3 5u j /3 + 0u j 0u j+/3 + 5u j+/3 u j + Oh 3 h/3 [ u j /3 + 6u j 30u j+/3 + 6u j+/3 u j + Oh 4 34 h/3 [ u j + 6u j+/3 30u j+/3 + 6u j+ u j+4/3 + c u j /3 + 5u j 0u j+/3 + 0u j+/3 5u j+ + u j+4/3 + Oh 3 This scheme was run under the same conditions as the previous examples, with the exception that now a 6th order Runge Kutta scheme was used for time integration. Here the truncation error is of Oh 3, however this is a 4th order scheme and 5th order for c= The convergence results are presented in Figure b. It should be noted that by taking c = the coefficients of u j /3 and u j+4/3 are 0. Therefore the scheme is more compact than standard explicit 4th order scheme in the sense that the scheme depends only on one term, on each side, out of the three point block u j,u j+/3,u j+/3. Potentially, this thinner stencil helps in the derivation of boundary schemes for initial boundary value problems. 5 Summary In this paper we presented a few block finite difference schemes in which the actual errors are much smaller than their truncation errors. This reduction of error was achieved by constructing the truncation errors to be oscillatory and using the dissipative property of the scheme. 7
8 A comparison between standard and block finite difference schemes in terms of the number of points out of the cell and operation count is presented in Table. As can be seen, in terms of accuracy and computational cost the 3rd and 5th order schemes are between the standards nd and 4th order and 4th and 6th order schemes respectively. Scheme Points out the block Number of operations Scheme Points out the block Number of operations at each side + at each side + standard nd order 3 point block 3rd order 3 4 standard 4th order point block 3rd order standard 6th order point block 5th order point block 4th order 4 5 compact Table Comparison between standard finite difference schemes and the ones presented on Sections 3,4. If c = /4 in the two point block, 3rd order scheme 6, the leading term in the error is highly oscillatory, see 3. It was suggested by Jennifer K. Ryan [8] that this term can be filtered by post processing. In this technique the high frequency error terms are filtered using convolution with a proper kernel. This method was successfully applied to the discontinuous Galerkin method. As this filtering is done only once, after the final time step, the cost is minimal, see e.g. [9], [], [3]. Here we used a global spectral filter and the local filter suggested in [3]. As can be seen in figure 3, the filtered scheme is 4th order accurate. The difference between both kernels is that the global filter is more computationally expensive, ON logn for proper values of N, but is accurate for small values of N while the local filter requires only ON operation but is accurate only for large values of N. The investigation of which is the optimal filter for these schemes is a topic for future research spectral filter, slope: * local filter, slope: no filter, slope: Fig. 3 Convergence plots, log 0 E vs. log 0 N, for scheme 6, c /4. with no post processing, with spectral filter and with the filter suggested in [3]. It should be noted that that finite difference representations of discontinuous Galerkin schemes for the heat equation have a similar form to the schemes presented above. They may also present similar enhancing of accuracy. See []. This manuscript was the inspiration to the current work. As was also pointed out in [] the increasing of accuracy may be related to a phenomenon called Supra Convergence [7]. Further research of the properties and implementations of these schemes as well as the existence of these kinds of schemes for hyperbolic problems are topics for future research. 8
9 Acknowledgments The author would like to thank Jennifer K. Ryan, Chi Wang Shu and Sigal Gottlieb for the fruitful discussions and their help. The author would also like to thank the anonymous reviewers for their useful remarks. References. Saul Abarbanel, Adi Ditkowski, and Bertil Gustafsson, On error bounds of finite difference approximations to partial differential equationstemporal behavior and rate of convergence, Journal of Scientific Computing 5 000, no., Bernardo Cockburn, Mitchell Luskin, Chi-Wang Shu, and Endre Suli, Post-processing of galerkin methods for hyperbolic problems, Springer, , Enhanced accuracy by post-processing for finite element methods for hyperbolic equations, Mathematics of Computation 7 003, no. 4, Bertil Gustafsson, The convergence rate for difference approximations to mixed initial boundary value problems, Mathematics of Computation 9 975, no. 30, , The convergence rate for difference approximations to general mixed initial-boundary value problems, SIAM Journal on Numerical Analysis 8 98, no., Bertil Gustafsson, Heinz-Otto Kreiss, and Joseph Oliger, Time-dependent problems and difference methods, John Wiley & Sons, H-O Kreiss, Thomas A Manteuffel, B Swartz, B Wendroff, and AB White, Supra-convergent schemes on irregular grids, Mathematics of Computation , no. 76, Jennifer Ryan, Privet communucation. 9. Jennifer Ryan, Chi-Wang Shu, and Harold Atkins, Extension of a post processing technique for the discontinuous galerkin method for hyperbolic equations with application to an aeroacoustic problem, SIAM Journal on Scientific Computing 6 005, no. 3, Magnus Svärd and Jan Nordström, On the order of accuracy for difference approximations of initial-boundary value problems, Journal of Computational Physics 8 006, no., Mengping Zhang and Chi-Wang Shu, An analysis of three different formulations of the discontinuous galerkin method for diffusion equations, Mathematical Models and Methods in Applied Sciences 3 003, no. 03,
Numerical Solution of Initial Boundary Value Problems. Jan Nordström Division of Computational Mathematics Department of Mathematics
Numerical Solution of Initial Boundary Value Problems Jan Nordström Division of Computational Mathematics Department of Mathematics Overview Material: Notes + GUS + GKO + HANDOUTS Schedule: 5 lectures
More informationSaul Abarbanel; Half a century of scientific work. Bertil Gustafsson, Uppsala University
Saul Abarbanel; Half a century of scientific work Bertil Gustafsson, Uppsala University Grew up in Tel Aviv Served in Israeli Army during the War of Independence 1948 1950 MIT 1952 1959 Ph.D 1959, Theoretical
More informationNumerical Methods for Partial Differential Equations CAAM 452. Spring 2005
Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
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 informationOn Error Bounds of Finite Difference Approximations to Partial Differential EquationsTemporal Behavior and Rate of Convergence
Journal of Scientific Computing, Vol. 15, No. 1, 2000 On Error Bounds of Finite Difference Approximations to Partial Differential EquationsTemporal Behavior and Rate of Convergence Saul Abarbanel, 1 Adi
More informationREVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS
REVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS MAGNUS SVÄRD AND JAN NORDSTRÖM Abstract. High-order finite difference methods are efficient, easy to program, scale well in multiple
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 informationNumerical resolution of discontinuous Galerkin methods for time dependent. wave equations 1. Abstract
Numerical resolution of discontinuous Galerkin methods for time dependent wave equations Xinghui Zhong 2 and Chi-Wang Shu Abstract The discontinuous Galerkin DG method is known to provide good wave resolution
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
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 informationDevelopment and stability analysis of the inverse Lax-Wendroff boundary. treatment for central compact schemes 1
Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes François Vilar 2 and Chi-Wang Shu 3 Division of Applied Mathematics, Brown University, Providence,
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
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 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 informationStrict Stability of High-Order Compact Implicit Finite-Difference Schemes: The Role of Boundary Conditions for Hyperbolic PDEs, I
Journal of Computational Physics 160, 42 66 (2000) doi:10.1006/jcph.2000.6420, available online at http://www.idealibrary.com on Strict Stability of High-Order Compact Implicit Finite-Difference Schemes:
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
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 informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationNumerical Smoothness and PDE Error Analysis
Numerical Smoothness and PDE Error Analysis GLSIAM 2013 Central Michigan University Tong Sun Bowling Green State University April 20, 2013 Table of contents Why consider Numerical Smoothness? What is Numerical
More informationOn the Convergence Rates of Energy- Stable Finite-Difference Schemes
Department of Mathematics On the Convergence Rates of Energy- Stable Finite-Difference Schemes Magnus Svärd and Jan Nordström LiTH-MAT-R--217/14--SE Department of Mathematics Linköping University S-581
More informationCharacterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients
1 Characterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients by Guang Wei Yu Supervisor: D. W. Zingg April 13, 2013 Abstract This paper presents the
More informationWell-posedness, Stability and Conservation for a Discontinuous Interface Problem
Department of Mathematics Well-posedness, Stability and Conservation for a Discontinuous Interface Problem Cristina La Cognata and Jan Nordström LiTH-MAT-R--214/16--SE Department of Mathematics Linköping
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
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 informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
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 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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
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 informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
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 informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationDiscontinuous Galerkin method for hyperbolic equations with singularities
Discontinuous Galerkin method for hyperbolic equations with singularities Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Qiang Zhang, Yang Yang and Dongming Wei Outline Introduction
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
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 informationIntroduction In this paper we review existing and develop new local discontinuous Galerkin methods for solving time dependent partial differential equ
Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives Jue Yan and Chi-Wang Shu 3 Division of Applied Mathematics Brown University Providence, Rhode Island
More informationNUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY
NUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY Stefan Johansson Department of Information Technology Scientific Computing
More informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I 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 informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
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 informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
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 information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
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 informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More informationHigh Order Difference Approximations for the Linearized Euler Equations
High Order Difference Approximations for the Linearized Euler Equations Stefan Johansson Abstract The computers of today make it possible to do direct simulation of aeroacoustics, which is very computationally
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
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 informationA stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
International Journal of Computer Mathematics Vol. 87, No. 11, September 2010, 2588 2600 A stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
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 informationEfficient wave propagation on complex domains
Center for Turbulence Research Annual Research Briefs 2006 223 Efficient wave propagation on complex domains By K. Mattsson, F. Ham AND G. Iaccarino 1. Motivation and objectives In many applications, such
More informationKrylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations
mathematics Article Krylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations Michael Machen and Yong-Tao Zhang * Department of Applied and Computational Mathematics and Statistics,
More informationOptimal diagonal-norm SBP operators
Optimal diagonal-norm SBP operators Ken Mattsson 1, Martin Almquist 1 and Mark H. Carpenter 2 1 Department of Information Technology, Uppsala University 2 Computational Aerosciences Branch, NASA Langley
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 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 informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationMethod of Lines. Received April 20, 2009; accepted July 9, 2009
Method of Lines Samir Hamdi, William E. Schiesser and Graham W. Griffiths * Ecole Polytechnique, France; Lehigh University, USA; City University,UK. Received April 20, 2009; accepted July 9, 2009 The method
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More 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 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 informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationSpectral analysis of the incompressible Navier-Stokes equations with different boundary conditions
Spectral analysis of the incompressible Navier-Stokes equations with different boundary conditions Cristina La Cognata, Jan Nordström Department of Mathematics, Computational Mathematics, Linköping University,
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 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 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 informationarxiv: v1 [math.na] 20 Nov 2013
REVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS MAGNUS SVÄRD AND JAN NORDSTRÖM arxiv:1311.4984v1 [math.na] 2 Nov 213 Abstract. High-order finite difference methods are efficient,
More informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
More informationA Wavelet-based Method for Overcoming the Gibbs Phenomenon
AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 A Wavelet-based Method for Overcoming the Gibbs Phenomenon NATANIEL GREENE Department of Mathematics Computer
More informationError formulas for divided difference expansions and numerical differentiation
Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
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 informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
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 informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationThe Heat Equation John K. Hunter February 15, The heat equation on a circle
The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate
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 informationDiagonal-norm upwind SBP operators
Diagonal-norm upwind SBP operators Ken Mattsson June 8, 16 Abstract High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
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 informationStability of the fourth order Runge-Kutta method for time-dependent partial. differential equations 1. Abstract
Stability of the fourth order Runge-Kutta method for time-dependent partial differential equations 1 Zheng Sun 2 and Chi-Wang Shu 3 Abstract In this paper, we analyze the stability of the fourth order
More information