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

We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multidimensional 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² - 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, we are able to prove consistency in the sense that if the discrete solutions are regular independently of the mesh size, then they converge to a classical solution. Finally, we present results from numerical experiments in space dimension d = 2.