Author:
Carrillo José A.,Craig Katy,Wang Li,Wei Chaozhen
Abstract
AbstractCombining the classical theory of optimal transport with modern operator splitting techniques, we develop a new numerical method for nonlinear, nonlocal partial differential equations, arising in models of porous media, materials science, and biological swarming. Our method proceeds as follows: first, we discretize in time, either via the classical JKO scheme or via a novel Crank–Nicolson-type method we introduce. Next, we use the Benamou–Brenier dynamical characterization of the Wasserstein distance to reduce computing the solution of the discrete time equations to solving fully discrete minimization problems, with strictly convex objective functions and linear constraints. Third, we compute the minimizers by applying a recently introduced, provably convergent primal dual splitting scheme for three operators (Yan in J Sci Comput 1–20, 2018). By leveraging the PDEs’ underlying variational structure, our method overcomes stability issues present in previous numerical work built on explicit time discretizations, which suffer due to the equations’ strong nonlinearities and degeneracies. Our method is also naturally positivity and mass preserving and, in the case of the JKO scheme, energy decreasing. We prove that minimizers of the fully discrete problem converge to minimizers of the spatially continuous, discrete time problem as the spatial discretization is refined. We conclude with simulations of nonlinear PDEs and Wasserstein geodesics in one and two dimensions that illustrate the key properties of our approach, including higher-order convergence our novel Crank–Nicolson-type method, when compared to the classical JKO method.
Publisher
Springer Science and Business Media LLC
Subject
Applied Mathematics,Computational Theory and Mathematics,Computational Mathematics,Analysis
Reference106 articles.
1. L. Ambrosio, N. Fusco, and D. Pallara, Functions of bounded variation and free discontinuity problems, vol. 254, Clarendon Press Oxford, 2000.
2. L. Ambrosio, N. Gigli, and G. Savaré, Gradient flows: in metric spaces and in the space of probability measures, Springer Science & Business Media, 2008.
3. R. Bailo, J. A. Carrillo, and J. Hu, Fully discrete positivity-preserving and energy-dissipating schemes for aggregation–diffusion equations with a gradient-flow structure, Commun. Math. Sci., 18 (2020), pp. 1259–1303.
4. D. Balagué, J. A. Carrillo, T. Laurent, and G. Raoul, Dimensionality of local minimizers of the interaction energy, Arch. Ration. Mech. Anal., 209 (2013), pp. 1055–1088.
5. D. Balagué, J. A. Carrillo, T. Laurent, and G. Raoul, Nonlocal interactions by repulsive–attractive potentials: radial ins/stability, Phys. D, 260 (2013), pp. 5–25.
Cited by
26 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献