A fully discrete variational scheme for solving nonlinear Fokker-Planck equations in higher space dimensions

Abstract

We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multi-dimensional unit cube. This discretization is based on two structural properties of these equations: the first is the representation as a gradient flow of an entropy functional in the $L^2$-Wasserstein metric, the second is the Lagrangian nature, meaning that solutions can be written as the push forward transformation of the initial density under suitable flow maps. The resulting numerical scheme is entropy diminishing and mass conserving. Further, the scheme is weakly stable, which allows us to prove convergence under certain regularity assumptions. Finally, we present results from numerical experiments in space dimension $d=2$.

21 pages, 2 figures

Cited by in corpus (5)

• Numerical Study of a Particle Method for Gradient Flows
• Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms
• A Lagrangian scheme for the solution of nonlinear diffusion equations using moving simplex meshes
• Convergence of a particle method for diffusive gradient flows in one dimension
• Discrete approximation of the minimizing movement scheme for evolution equations of Wasserstein gradient flow type with nonlinear mobility