Abstract
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 L2 - 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.
Original language | English |
---|---|
Pages (from-to) | 419-423 |
Number of pages | 5 |
Journal | SIAM Journal on Numerical Analysis |
Volume | 55 |
Issue number | 1 |
DOIs | |
State | Published - 2017 |
Keywords
- Fokker-Planck equation
- Gradient flow
- Lagrangian discretization
- Structure preservation
- Wasserstein metric