High order numerical schemes for one-dimension non-local conservation laws

TL;DR

This paper develops high-order numerical schemes, specifically Discontinuous Galerkin and FV-WENO, to accurately approximate solutions of one-dimensional non-local conservation laws, addressing oscillations and stability issues in applications like traffic flow and sedimentation.

Contribution

It introduces high-order DG and FV-WENO schemes tailored for non-local conservation laws, improving accuracy and stability over classical methods.

Findings

DG schemes yield the best numerical accuracy.

FV-WENO schemes allow larger time steps.

Quadratic reconstructions are essential for high-order convolution evaluation.

Abstract

This paper focuses on the numerical approximation of the solutions of non-local conservation laws in one space dimension. These equations are motivated by two distinct applications, namely a traffic flow model in which the mean velocity depends on a weighted mean of the downstream traffic density, and a sedimentation model where either the solid phase velocity or the solid-fluid relative velocity depends on the concentration in a neighborhood. In both models, the velocity is a function of a convolution product between the unknown and a kernel function with compact support. It turns out that the solutions of such equations may exhibit oscillations that are very difficult to approximate using classical first-order numerical schemes. We propose to design Discontinuous Galerkin (DG) schemes and Finite Volume WENO (FV-WENO) schemes to obtain high-order approximations. As we will see, the DG schemes give the best numerical results but their CFL condition is very restrictive. On the contrary, FV-WENO schemes can be used with larger time steps. We will see that the evaluation of the convolution terms necessitates the use of quadratic polynomials reconstructions in each cell in order to obtain the high-order accuracy with the FV-WENO approach. Simulations using DG and FV-WENO schemes are presented for both applications.