Discretization of flux-limited gradient flows: $\Gamma$-convergence and numerical schemes
HTML articles powered by AMS MathViewer
- by Daniel Matthes and Benjamin Söllner;
- Math. Comp. 89 (2020), 1027-1057
- DOI: https://doi.org/10.1090/mcom/3492
- Published electronically: November 19, 2019
- HTML | PDF | Request permission
Abstract:
We study a discretization in space and time for a class of nonlinear diffusion equations with flux limitation. That class contains the so-called relativistic heat equation, as well as other gradient flows of Renyi entropies with respect to transportation metrics with finite maximal velocity. Discretization in time is performed with the JKO method, thus preserving the variational structure of the gradient flow. This is combined with an entropic regularization of the transport distance, which allows for an efficient numerical calculation of the JKO minimizers. Solutions to the fully discrete equations are entropy dissipating, mass conserving, and respect the finite speed of propagation of support.
First, we give a proof of $\Gamma$-convergence of the infinite chain of JKO steps in the joint limit of infinitely refined spatial discretization and vanishing entropic regularization. The singularity of the cost function makes the construction of the recovery sequence significantly more difficult than in the $L^p$-Wasserstein case. Second, we define a practical numerical method by combining the JKO time discretization with a “light speed” solver for the spatially discrete minimization problem using Dykstra’s algorithm, and demonstrate its efficiency in a series of experiments.
References
- Martial Agueh, Existence of solutions to degenerate parabolic equations via the Monge-Kantorovich theory, Adv. Differential Equations 10 (2005), no. 3, 309–360. MR 2123134
- Luigi Ambrosio, Nicola Gigli, and Giuseppe Savaré, Gradient flows in metric spaces and in the space of probability measures, 2nd ed., Lectures in Mathematics ETH Zürich, Birkhäuser Verlag, Basel, 2008. MR 2401600
- F. Andreu, V. Caselles, and J. M. Mazón, Some regularity results on the ‘relativistic’ heat equation, J. Differential Equations 245 (2008), no. 12, 3639–3663. MR 2462697, DOI 10.1016/j.jde.2008.06.024
- Fuensanta Andreu, Vicent Caselles, Jose M. Mazón, and Salvador Moll, The Dirichlet problem associated to the relativistic heat equation, Math. Ann. 347 (2010), no. 1, 135–199. MR 2593288, DOI 10.1007/s00208-009-0428-3
- Jean-David Benamou, Guillaume Carlier, Quentin Mérigot, and Édouard Oudet, Discretization of functionals involving the Monge-Ampère operator, Numer. Math. 134 (2016), no. 3, 611–636. MR 3555350, DOI 10.1007/s00211-015-0781-y
- Jérôme Bertrand and Marjolaine Puel, The optimal mass transport problem for relativistic costs, Calc. Var. Partial Differential Equations 46 (2013), no. 1-2, 353–374. MR 3016512, DOI 10.1007/s00526-011-0485-9
- Adrien Blanchet, Vincent Calvez, and José A. Carrillo, Convergence of the mass-transport steepest descent scheme for the subcritical Patlak-Keller-Segel model, SIAM J. Numer. Anal. 46 (2008), no. 2, 691–721. MR 2383208, DOI 10.1137/070683337
- Yann Brenier, Extended Monge-Kantorovich theory, Optimal transportation and applications (Martina Franca, 2001) Lecture Notes in Math., vol. 1813, Springer, Berlin, 2003, pp. 91–121. MR 2006306, DOI 10.1007/978-3-540-44857-0_{4}
- Guillaume Carlier, Vincent Duval, Gabriel Peyré, and Bernhard Schmitzer, Convergence of entropic schemes for optimal transport and gradient flows, SIAM J. Math. Anal. 49 (2017), no. 2, 1385–1418. MR 3635459, DOI 10.1137/15M1050264
- J. A. Carrillo, V. Caselles, and S. Moll, On the relativistic heat equation in one space dimension, Proc. Lond. Math. Soc. (3) 107 (2013), no. 6, 1395–1423. MR 3149851, DOI 10.1112/plms/pdt015
- José Antonio Carrillo, Katy Craig, and Francesco S. Patacchini, A blob method for diffusion, Calc. Var. Partial Differential Equations 58 (2019), no. 2, Paper No. 53, 53. MR 3913840, DOI 10.1007/s00526-019-1486-3
- José A. Carrillo, Bertram Düring, Daniel Matthes, and David S. McCormick, A Lagrangian scheme for the solution of nonlinear diffusion equations using moving simplex meshes, J. Sci. Comput. 75 (2018), no. 3, 1463–1499. MR 3798109, DOI 10.1007/s10915-017-0594-5
- 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/10), no. 6, 4305–4329. MR 2566595, DOI 10.1137/080739574
- V. Caselles, Convergence of the ‘relativistic’ heat equation to the heat equation as $c\to \infty$, Publ. Mat. 51 (2007), no. 1, 121–142. MR 2307149, DOI 10.5565/PUBLMAT_{5}1107_{0}6
- Lénaïc Chizat, Gabriel Peyré, Bernhard Schmitzer, and François-Xavier Vialard, Scaling algorithms for unbalanced optimal transport problems, Math. Comp. 87 (2018), no. 314, 2563–2609. MR 3834678, DOI 10.1090/mcom/3303
- M. Cuturi, Sinkhorn distances: Lightspeed computation of optimal transport. In Advances in neural information processing systems, pages 2292–2300, 2013.
- Richard Jordan, David Kinderlehrer, and Felix Otto, The variational formulation of the Fokker-Planck equation, SIAM J. Math. Anal. 29 (1998), no. 1, 1–17. MR 1617171, DOI 10.1137/S0036141096303359
- Oliver Junge, Daniel Matthes, and Horst Osberger, A fully discrete variational scheme for solving nonlinear Fokker-Planck equations in multiple space dimensions, SIAM J. Numer. Anal. 55 (2017), no. 1, 419–443. MR 3614682, DOI 10.1137/16M1056560
- O. Junge and B. Söllner, A convergent Lagrangian discretization of $p$-Laplace equations in one space dimension, Archive preprint arxiv:1906.01321, 2019.
- Philippe Laurençot and Bogdan-Vasile Matioc, A gradient flow approach to a thin film approximation of the Muskat problem, Calc. Var. Partial Differential Equations 47 (2013), no. 1-2, 319–341. MR 3044141, DOI 10.1007/s00526-012-0520-5
- Jan Maas and Daniel Matthes, Long-time behavior of a finite volume discretization for a fourth order diffusion equation, Nonlinearity 29 (2016), no. 7, 1992–2023. MR 3521637, DOI 10.1088/0951-7715/29/7/1992
- Daniel Matthes and Horst Osberger, Convergence of a variational Lagrangian scheme for a nonlinear drift diffusion equation, ESAIM Math. Model. Numer. Anal. 48 (2014), no. 3, 697–726. MR 3177862, DOI 10.1051/m2an/2013126
- Daniel Matthes and Horst Osberger, A convergent Lagrangian discretization for a nonlinear fourth-order equation, Found. Comput. Math. 17 (2017), no. 1, 73–126. MR 3600850, DOI 10.1007/s10208-015-9284-6
- Robert J. McCann and Marjolaine Puel, Constructing a relativistic heat flow by transport time steps, Ann. Inst. H. Poincaré C Anal. Non Linéaire 26 (2009), no. 6, 2539–2580 (English, with English and French summaries). MR 2569908, DOI 10.1016/j.anihpc.2009.06.006
- Felix Otto, Lubrication approximation with prescribed nonzero contact angle, Comm. Partial Differential Equations 23 (1998), no. 11-12, 2077–2164. MR 1662172, DOI 10.1080/03605309808821411
- Felix Otto, The geometry of dissipative evolution equations: the porous medium equation, Comm. Partial Differential Equations 26 (2001), no. 1-2, 101–174. MR 1842429, DOI 10.1081/PDE-100002243
- Gabriel Peyré, Entropic approximation of Wasserstein gradient flows, SIAM J. Imaging Sci. 8 (2015), no. 4, 2323–2351. MR 3413589, DOI 10.1137/15M1010087
- G. Peyré and M. Cuturi, Computational optimal transport, to appear in Foundations and Trends in Machine Learning, 2018.
- P. Rosenau, Tempered diffusion: A transport process with propagating fronts and inertial delay, Physical Review A, 46:R7371, (1992), no. 12.
Bibliographic Information
- Daniel Matthes
- Affiliation: Technische Universität München, Zentrum Mathematik Boltzmannstr. 3 D-85748 Garching, Germany
- MR Author ID: 722279
- Email: matthes@ma.tum.de
- Benjamin Söllner
- Affiliation: Technische Universität München, Zentrum Mathematik Boltzmannstr. 3 D-85748 Garching, Germany
- Email: soellneb@ma.tum.de
- Received by editor(s): January 31, 2019
- Received by editor(s) in revised form: August 7, 2019
- Published electronically: November 19, 2019
- Additional Notes: This research was supported by the DFG Collaborative Research Center TRR 109 “Discretization in Geometry and Dynamics”
- © Copyright 2019 American Mathematical Society
- Journal: Math. Comp. 89 (2020), 1027-1057
- MSC (2010): Primary 35K20, 35K65, 49M29, 65M12
- DOI: https://doi.org/10.1090/mcom/3492
- MathSciNet review: 4063311