Abstract

We make an analytic study of the diffusive, dispersive and overall errors, which arise when using semi-implicit semi-Lagrangian (SISL) finite difference methods to approximate those travelling wave solutions of the one-dimensional Burgers' equation with small diffusion, which develop sharp fronts. For the case of a fixed uniform spatial mesh, with piecewise linear interpolation, a backward error analysis approach is used to construct a precise formal analytic description of the front profile of the numerical approximation to this solution. From this description it is possible to obtain precise estimates of the front width and the front speed in terms of the spatial and temporal step size and to express the overall solution error in terms of these. These formal estimates agree closely with numerical calculations, both qualitatively and quantitatively, and display a roughly periodic behaviour as the number Nx of mesh points increases, and the CFL number passes through integer values. In particular, they show that despite the otherwise poor resolution of the method, the front width is closely approximated when the CFL number is close to an integer, and the front speed is closely approximated when it is close to a half integer. The overall L2 error also shows super-convergence for certain values of Nx. This possibly motivates doing two calculations with different Nx when using the SISL method on such problems to separately minimise the diffusive and dispersive errors. Similar errors in the front width and speed are observed for a number of different interpolation schemes with and without flux limiters.

LanguageEnglish
Pages261-282
Number of pages22
JournalApplied Numerical Mathematics
Volume145
Early online date21 Jun 2019
DOIs
StatusE-pub ahead of print - 21 Jun 2019

Keywords

  • Burgers' equation
  • Error estimates
  • Modified equation
  • Numerical weather prediction
  • Semi-implicit
  • Semi-Lagrangian

ASJC Scopus subject areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

Cite this

Error estimates for semi-Lagrangian finite difference methods applied to Burgers' equation in one dimension. / Cook, Stephen; Budd, Chris; Hill, Adrian; Melvin, Thomas.

In: Applied Numerical Mathematics, Vol. 145, 01.11.2019, p. 261-282.

Research output: Contribution to journalArticle

@article{561e3fef7a56494d9f50410233c95d47,
title = "Error estimates for semi-Lagrangian finite difference methods applied to Burgers' equation in one dimension",
abstract = "We make an analytic study of the diffusive, dispersive and overall errors, which arise when using semi-implicit semi-Lagrangian (SISL) finite difference methods to approximate those travelling wave solutions of the one-dimensional Burgers' equation with small diffusion, which develop sharp fronts. For the case of a fixed uniform spatial mesh, with piecewise linear interpolation, a backward error analysis approach is used to construct a precise formal analytic description of the front profile of the numerical approximation to this solution. From this description it is possible to obtain precise estimates of the front width and the front speed in terms of the spatial and temporal step size and to express the overall solution error in terms of these. These formal estimates agree closely with numerical calculations, both qualitatively and quantitatively, and display a roughly periodic behaviour as the number Nx of mesh points increases, and the CFL number passes through integer values. In particular, they show that despite the otherwise poor resolution of the method, the front width is closely approximated when the CFL number is close to an integer, and the front speed is closely approximated when it is close to a half integer. The overall L2 error also shows super-convergence for certain values of Nx. This possibly motivates doing two calculations with different Nx when using the SISL method on such problems to separately minimise the diffusive and dispersive errors. Similar errors in the front width and speed are observed for a number of different interpolation schemes with and without flux limiters.",
keywords = "Burgers' equation, Error estimates, Modified equation, Numerical weather prediction, Semi-implicit, Semi-Lagrangian",
author = "Stephen Cook and Chris Budd and Adrian Hill and Thomas Melvin",
year = "2019",
month = "6",
day = "21",
doi = "10.1016/j.apnum.2019.06.012",
language = "English",
volume = "145",
pages = "261--282",
journal = "Applied Numerical Mathematics",
issn = "0168-9274",
publisher = "Elsevier",

}

TY - JOUR

T1 - Error estimates for semi-Lagrangian finite difference methods applied to Burgers' equation in one dimension

AU - Cook, Stephen

AU - Budd, Chris

AU - Hill, Adrian

AU - Melvin, Thomas

PY - 2019/6/21

Y1 - 2019/6/21

N2 - We make an analytic study of the diffusive, dispersive and overall errors, which arise when using semi-implicit semi-Lagrangian (SISL) finite difference methods to approximate those travelling wave solutions of the one-dimensional Burgers' equation with small diffusion, which develop sharp fronts. For the case of a fixed uniform spatial mesh, with piecewise linear interpolation, a backward error analysis approach is used to construct a precise formal analytic description of the front profile of the numerical approximation to this solution. From this description it is possible to obtain precise estimates of the front width and the front speed in terms of the spatial and temporal step size and to express the overall solution error in terms of these. These formal estimates agree closely with numerical calculations, both qualitatively and quantitatively, and display a roughly periodic behaviour as the number Nx of mesh points increases, and the CFL number passes through integer values. In particular, they show that despite the otherwise poor resolution of the method, the front width is closely approximated when the CFL number is close to an integer, and the front speed is closely approximated when it is close to a half integer. The overall L2 error also shows super-convergence for certain values of Nx. This possibly motivates doing two calculations with different Nx when using the SISL method on such problems to separately minimise the diffusive and dispersive errors. Similar errors in the front width and speed are observed for a number of different interpolation schemes with and without flux limiters.

AB - We make an analytic study of the diffusive, dispersive and overall errors, which arise when using semi-implicit semi-Lagrangian (SISL) finite difference methods to approximate those travelling wave solutions of the one-dimensional Burgers' equation with small diffusion, which develop sharp fronts. For the case of a fixed uniform spatial mesh, with piecewise linear interpolation, a backward error analysis approach is used to construct a precise formal analytic description of the front profile of the numerical approximation to this solution. From this description it is possible to obtain precise estimates of the front width and the front speed in terms of the spatial and temporal step size and to express the overall solution error in terms of these. These formal estimates agree closely with numerical calculations, both qualitatively and quantitatively, and display a roughly periodic behaviour as the number Nx of mesh points increases, and the CFL number passes through integer values. In particular, they show that despite the otherwise poor resolution of the method, the front width is closely approximated when the CFL number is close to an integer, and the front speed is closely approximated when it is close to a half integer. The overall L2 error also shows super-convergence for certain values of Nx. This possibly motivates doing two calculations with different Nx when using the SISL method on such problems to separately minimise the diffusive and dispersive errors. Similar errors in the front width and speed are observed for a number of different interpolation schemes with and without flux limiters.

KW - Burgers' equation

KW - Error estimates

KW - Modified equation

KW - Numerical weather prediction

KW - Semi-implicit

KW - Semi-Lagrangian

UR - http://www.scopus.com/inward/record.url?scp=85068035562&partnerID=8YFLogxK

U2 - 10.1016/j.apnum.2019.06.012

DO - 10.1016/j.apnum.2019.06.012

M3 - Article

VL - 145

SP - 261

EP - 282

JO - Applied Numerical Mathematics

T2 - Applied Numerical Mathematics

JF - Applied Numerical Mathematics

SN - 0168-9274

ER -