### 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 N_{x} 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 L_{2} error also shows super-convergence for certain values of N_{x}. This possibly motivates doing two calculations with different N_{x} 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.

Original language | English |
---|---|

Pages (from-to) | 261-282 |

Number of pages | 22 |

Journal | Applied Numerical Mathematics |

Volume | 145 |

Early online date | 21 Jun 2019 |

DOIs | |

Publication status | E-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.

Research output: Contribution to journal › Article

*Applied Numerical Mathematics*, vol. 145, pp. 261-282. https://doi.org/10.1016/j.apnum.2019.06.012

}

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

JF - Applied Numerical Mathematics

SN - 0168-9274

ER -