Fully discrete energy-dissipative and conservative discrete gradient particle methods for a class of continuity equations

TL;DR

This paper develops fully discrete particle methods for certain continuity equations that preserve energy dissipation and conservation laws, ensuring accurate and stable numerical simulations of complex physical systems.

Contribution

It introduces a novel formulation of particle methods using discrete gradient integrators that simultaneously dissipate energy and conserve mass, momentum, and kinetic energy at the fully discrete level.

Findings

Methods successfully dissipate energy and conserve mass in simulations.

Approach conserves momentum and kinetic energy for Landau equation.

Numerical examples demonstrate stability and accuracy of the methods.

Abstract

Structure-preserving particle methods have recently been proposed for a class of nonlinear continuity equations, including aggregation-diffusion equation in [J. Carrillo, K. Craig, F. Patacchini, Calc. Var., 58 (2019), pp. 53] and the Landau equation in [J. Carrillo, J. Hu., L. Wang, J. Wu, J. Comput. Phys. X, 7 (2020), pp. 100066]. One common feature to these equations is that they both admit some variational formulation, which upon proper regularization, leads to particle approximations dissipating the energy and conserving some quantities simultaneously at the semi-discrete level. In this paper, we formulate continuity equations with a density dependent bilinear form associated with the variational derivative of the energy functional and prove that appropriate particle methods satisfy a compatibility condition with its regularized energy. This enables us to utilize discrete gradient time integrators and show that the energy can be dissipated and the mass conserved simultaneously at the fully discrete level. In the case of the Landau equation, we prove that our approach also conserves the momentum and kinetic energy at the fully discrete level. Several numerical examples are presented to demonstrate the dissipative and conservative properties of our proposed method.