Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

TL;DR

This paper develops and analyzes discontinuous Galerkin methods for one-dimensional nonlocal diffusion problems, ensuring convergence to classical schemes as the nonlocal interaction radius vanishes, with rigorous stability and error proofs.

Contribution

It introduces DG methods with penalty terms for nonlocal diffusion that are asymptotically compatible and provides theoretical analysis and numerical validation.

Findings

Methods are stable and accurate for various horizons.

Proposed schemes converge to classical DG methods as horizon vanishes.

Numerical experiments confirm theoretical error estimates.

Abstract

In this paper, we consider a class of discontinuous Galerkin (DG) methods for one-dimensional nonlocal diffusion (ND) problems. The nonlocal models, which are integral equations, are widely used in describing many physical phenomena with long-range interactions. The ND problem is the nonlocal analog of the classic diffusion problem, and as the interaction radius (horizon) vanishes, then the nonlocality disappears and the ND problem converges to the classic diffusion problem. Under certain conditions, the exact solution to the ND problem may exhibit discontinuities, setting it apart from the classic diffusion problem. Since the DG method shows its great advantages in resolving problems with discontinuities in computational fluid dynamics over the past several decades, it is natural to adopt the DG method to compute the ND problems. Based on [Du-Ju-Lu-Tian-CAMC2020], we develop the DG methods with different penalty terms, ensuring that the proposed DG methods have local counterparts as the horizon vanishes. This indicates the proposed methods will converge to the existing DG schemes as the horizon vanishes, which is crucial for achieving asymptotic compatibility. Rigorous proofs are provided to demonstrate the stability, error estimates, and asymptotic compatibility of the proposed DG schemes. To observe the effect of the nonlocal diffusion, we also consider the time-dependent convection-diffusion problems with nonlocal diffusion. We conduct several numerical experiments, including accuracy tests and Burgers' equation with nonlocal diffusion, and various horizons are taken to show the good performance of the proposed algorithm and validate the theoretical findings.