Convergence rates for upwind schemes with rough coefficients

Abstract

This paper is concerned with the numerical analysis of the explicit upwind finite volume scheme for numerically solving continuity equations. We are interested in the case where the advecting velocity field has spatial Sobolev regularity and initial data are merely integrable. We estimate the error between approximate solutions constructed by the upwind scheme and distributional solutions of the continuous problem in a Kantorovich-Rubinstein distance, which was recently used for stability estimates for the continuity equation by Seis [23]. Restricted to Cartesian meshes, our estimate shows that the rate of weak convergence is at least of order $1/2$ in the mesh size. The proof relies on a probabilistic interpretation of the upwind scheme Delarue and Lagoutière [9]. We complement the weak convergence result with an example that illustrates that for rough initial data no rates can be expected in strong norms. The same example suggests that the weak order $1/2$ rate is optimal.

Revised and extended version with optimality section containing lower bounds on the convergence rate

Cited by in corpus (6)

• Nonlocal-interaction equation on graphs: gradient flow structure and continuum limit
• Analysis of the implicit upwind finite volume scheme with rough coefficients
• Variational convergence of the Scharfetter-Gummel scheme to the aggregation-diffusion equation and vanishing diffusion limit
• Error estimates for a finite volume scheme for advection-diffusion equations with rough coefficients
• Second-order finite difference approximations of the upper-convected time derivative
• Convergent finite difference schemes for stochastic transport equations