Convergence of a variational lagrangian scheme for a nonlinear drift diffusion equation
ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 48 (2014) no. 3, pp. 697-726.

We study a Lagrangian numerical scheme for solution of a nonlinear drift diffusion equation on an interval. The discretization is based on the equation's gradient flow structure with respect to the Wasserstein distance. The scheme inherits various properties from the continuous flow, like entropy monotonicity, mass preservation, metric contraction and minimum/ maximum principles. As the main result, we give a proof of convergence in the limit of vanishing mesh size under a CFL-type condition. We also present results from numerical experiments.

DOI : https://doi.org/10.1051/m2an/2013126
Classification : 65N12,  49Q20,  65N06
Mots clés : lagrangian discretization, nonlinear Fokker-Planck equation, gradient flow, Wasserstein metric
@article{M2AN_2014__48_3_697_0,
author = {Matthes, Daniel and Osberger, Horst},
title = {Convergence of a variational lagrangian scheme for a nonlinear drift diffusion equation},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
pages = {697--726},
publisher = {EDP-Sciences},
volume = {48},
number = {3},
year = {2014},
doi = {10.1051/m2an/2013126},
zbl = {1293.65119},
mrnumber = {3177862},
language = {en},
url = {http://www.numdam.org/item/M2AN_2014__48_3_697_0/}
}
Matthes, Daniel; Osberger, Horst. Convergence of a variational lagrangian scheme for a nonlinear drift diffusion equation. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 48 (2014) no. 3, pp. 697-726. doi : 10.1051/m2an/2013126. http://www.numdam.org/item/M2AN_2014__48_3_697_0/

[1] L. Ambrosio, N. Gigli and G. Savaré, Gradient flows in metric spaces and in the space of probability measures. Lectures in Mathematics. ETH Zürich, Birkhäuser Verlag, Basel (2005). | MR 2129498 | Zbl 1145.35001

[2] L. Ambrosio, S. Lisini and G. Savaré, Stability of flows associated to gradient vector fields and convergence of iterated transport maps. Manuscripta Math. 121 (2006) 1-50. | MR 2258529 | Zbl 1099.49027

[3] J.-D. Benamou and Y. Brenier, A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem. Numer. Math. 84 (2000) 375-393. | MR 1738163 | Zbl 0968.76069

[4] A. Blanchet, V. Calvez and J.A. Carrillo, Convergence of the mass-transport steepest descent scheme for the subcritical Patlak-Keller-Segel model. SIAM J. Numer. Anal. 46 (2008) 691-721. | MR 2383208 | Zbl 1205.65332

[5] C.J. Budd, G.J. Collins, W.Z. Huang and R.D. Russell, Self-similar numerical solutions of the porous-medium equation using moving mesh methods. R. Soc. Lond. Philos. Trans. Ser. A Math. Phys. Eng. Sci. 357 (1999) 1047-1077. | MR 1694702 | Zbl 0931.76079

[6] M. Burger, J.A. Carrillo and M.-T. Wolfram, A mixed finite element method for nonlinear diffusion equations. Kinet. Relat. Models 3 (2010) 59-83. | MR 2580954 | Zbl 1194.35026

[7] J.A. Carrillo and J.S. Moll, Numerical simulation of diffusive and aggregation phenomena in nonlinear continuity equations by evolving diffeomorphisms. SIAM J. Sci. Comput. 31 (2009/2010) 4305-4329. | MR 2566595 | Zbl 1205.65265

[8] F. Cavalli and G. Naldi, A Wasserstein approach to the numerical solution of the one-dimensional Cahn-Hilliard equation. Kinet. Relat. Models 3 (2010) 123-142. | MR 2580956 | Zbl 1191.80042

[9] B. Düring, D. Matthes and J.P. Milišić, A gradient flow scheme for nonlinear fourth order equations. Discrete Contin. Dyn. Syst. Ser. B 14 (2010) 935-959. | MR 2670179 | Zbl 1206.65221

[10] L.C. Evans, O. Savin and W. Gangbo, Diffeomorphisms and nonlinear heat flows. SIAM J. Math. Anal. 37 (2005) 737-751. | MR 2191774 | Zbl 1096.35061

[11] E. Giusti, Minimal surfaces and functions of bounded variation, vol. 80, Monographs in Mathematics. Birkhäuser Verlag, Basel (1984). | MR 775682 | Zbl 0545.49018

[12] L. Gosse and G. Toscani, Identification of asymptotic decay to self-similarity for one-dimensional filtration equations. SIAM J. Numer. Anal. 43 (2006) 2590-2606 (electronic). | MR 2206449 | Zbl 1145.76048

[13] R. Jordan, D. Kinderlehrer and F. Otto, The variational formulation of the Fokker-Planck equation. SIAM J. Math. Anal. 29 (1998) 1-17. | Zbl 0915.35120

[14] D. Kinderlehrer and N.J. Walkington, Approximation of parabolic equations using the Wasserstein metric. ESAIM: M2AN 33 (1999) 837-852. | Numdam | Zbl 0936.65121

[15] M. Leven, Gradientenfluß-basierte diskretisierung parabolischer gleichungen, diplomarbeit, Universität Bonn (2002).

[16] R.C. Maccamy and E. Socolovsky, A numerical procedure for the porous media equation. Hyperbolic partial differential equations, II. Comput. Math. Appl. 11 (1985) 315-319. | Zbl 0608.76083

[17] R.J. Mccann, A convexity principle for interacting gases. Adv. Math. 128 (1997) 153-179. | Zbl 0901.49012

[18] F. Otto, The geometry of dissipative evolution equations: the porous medium equation. Comm. Partial Differential Equations 26 (2001) 101-174. | Zbl 0984.35089

[19] T. Roessler, Discretizing the porous medium equation based on its gradient flow structure - a consistency paradox, Technical report 150, Sonderforschungsbereich 611, May 2004.Available online at http://sfb611.iam.uni-bonn.de/uploads/150-komplett.pdf.

[20] R. Rossi and G. Savaré, Tightness, integral equicontinuity and compactness for evolution problems in Banach spaces. Ann. Sc. Norm. Super. Pisa Cl. Sci. 2 (2003) 395-431. | Numdam | MR 2005609 | Zbl 1150.46014

[21] G. Russo, Deterministic diffusion of particles. Commun. Pure Appl. Math. 43 (1990) 697-733. | MR 1059326 | Zbl 0713.65089

[22] S. Serfaty, Gamma-convergence of gradient flows on Hilbert and metric spaces and applications, Discrete Contin. Dyn. Syst. 31 (2011) 1427-1451. | MR 2836361 | Zbl 1239.35015

[23] C. Villani, Topics in optimal transportation, in vol. 58 of Graduate Studies in Mathematics. American Mathematical Society, Providence, RI (2003). | MR 1964483 | Zbl 1106.90001

[24] M. Westdickenberg and J. Wilkening, Variational particle schemes for the porous medium equation and for the system of isentropic Euler equations. ESAIM: M2AN 44 (2010) 133-166. | Numdam | MR 2647756