Premium
Application of Richardson extrapolation to the numerical solution of partial differential equations
Author(s) -
Burg Clarence,
Erwin Taylor
Publication year - 2009
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.20375
Subject(s) - richardson extrapolation , mathematics , extrapolation , discretization , order of accuracy , partial differential equation , numerical partial differential equations , numerical analysis , exponential integrator , method of lines , differential equation , mathematical analysis , ordinary differential equation , differential algebraic equation
Richardson extrapolation is a methodology for improving the order of accuracy of numerical solutions that involve the use of a discretization size h . By combining the results from numerical solutions using a sequence of related discretization sizes, the leading order error terms can be methodically removed, resulting in higher order accurate results. Richardson extrapolation is commonly used within the numerical approximation of partial differential equations to improve certain predictive quantities such as the drag or lift of an airfoil, once these quantities are calculated on a sequence of meshes, but it is not widely used to determine the numerical solution of partial differential equations. Within this article, Richardson extrapolation is applied directly to the solution algorithm used within existing numerical solvers of partial differential equations to increase the order of accuracy of the numerical result without referring to the details of the methodology or its implementation within the numerical code. Only the order of accuracy of the existing solver and certain interpolations required to pass information between the mesh levels are needed to improve the order of accuracy and the overall solution accuracy. Using the proposed methodology, Richardson extrapolation is used to increase the order of accuracy of numerical solutions of the linear heat and wave equations and of the nonlinear St. Venant equations in one‐dimension. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009