A convergent Lagrangian discretization for <inline-formula><tex-math id="M1">\begin{document}$ p $\end{document}</tex-math></inline-formula>-Wasserstein and flux-limited diffusion equations

Abstract

We study a Lagrangian numerical scheme for solving a nonlinear drift diffusion equations of the form $ \partial_t u = \partial_x(u \cdot ({\sf c}^*)^\prime[\partial_x \mathit{h}^\prime(u)+ \mathit{v}^\prime]) $, like Fokker-Plank and $ q $-Laplace equations, on an interval. This scheme will consist of a spatio-temporal discretization founded on the formulation of the equation in terms of inverse distribution functions. It is based on the gradient flow structure of the equation with respect to optimal transport distances for a family of costs that are in some sense $ p $-Wasserstein like. Additionally we will show that, under a regularity assumption on the initial data, this also includes a family of discontinuous, flux-limiting cost inducing equations like Rosenau's relativistic heat equation. We show that this discretization inherits various properties from the continuous flow, like entropy monotonicity, mass preservation, a minimum/maximum principle and flux-limitation in the case of the corresponding cost. Convergence in the limit of vanishing mesh size will be proven as the main result. Finally we will present numerical experiments including a numerical convergence analysis.