Geometrically exact time‐integration mesh‐free schemes for advection‐diffusion problems derived from optimal transportation theory and their connection with particle methods

We develop an optimal transportation mesh-free particle method for advection-diffusion in which the concentration or density of the diffusive species is approximated by Dirac measures. We resort to an incremental variational principle for purposes of time discretization of the diffusive step. This principle characterizes the evolution of the density as a competition between the Wasserstein distance between two consecutive densities and entropy. Exploiting the structure of the Euler–Lagrange equations, we approximate the density as a collection of Diracs. The interpolation of the incremental transport map is effected through mesh-free max-ent interpolation. Remarkably, the resulting update is geometrically exact with respect to advection and volume. We present three-dimensional examples of application that illustrate the scope and robustness of the method.

[1]  D. Sulsky,et al.  A particle method for history-dependent materials , 1993 .

[2]  J. K. Chen,et al.  An improvement for tensile instability in smoothed particle hydrodynamics , 1999 .

[3]  Bernd Schmidt,et al.  Convergence Analysis of Meshfree Approximation Schemes , 2011, SIAM J. Numer. Anal..

[4]  David J. Benson,et al.  Momentum advection on a staggered mesh , 1992 .

[5]  Sara Daneri,et al.  Lecture notes on gradient flows and optimal transport , 2010, Optimal Transport.

[6]  J. Marsden,et al.  Asynchronous Variational Integrators , 2003 .

[7]  J. Marsden,et al.  Variational time integrators , 2004 .

[8]  J. Marsden,et al.  Symplectic-energy-momentum preserving variational integrators , 1999 .

[9]  J. Marsden,et al.  Discrete mechanics and variational integrators , 2001, Acta Numerica.

[10]  Magdalena Ortiz,et al.  Local maximum‐entropy approximation schemes: a seamless bridge between finite elements and meshfree methods , 2006 .

[11]  M. Knott,et al.  On the optimal mapping of distributions , 1984 .

[12]  D. Kinderlehrer,et al.  Dynamics of the fokker-planck equation , 1999 .

[13]  Y. Brenier.  Polar Factorization and Monotone Rearrangement of Vector-Valued Functions , 1991 .

[14]  D. Kinderlehrer,et al.  Free energy and the Fokker-Planck equation , 1997 .

[15]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[16]  A. Peirce.  Computer Methods in Applied Mechanics and Engineering , 2010 .

[17]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[18]  C. Villani.  Topics in Optimal Transportation , 2003 .

[19]  Michael Ortiz,et al.  Convergent meshfree approximation schemes of arbitrary order and smoothness , 2012 .

[20]  Wing Kam Liu,et al.  Reproducing kernel particle methods , 1995 .

[21]  Romesh C. Batra,et al.  Modified smoothed particle hydrodynamics method and its application to transient problems , 2004 .

[22]  Bo Li,et al.  Optimal transportation meshfree approximation schemes for fluid and plastic flows , 2010 .

[23]  Daniel Matthes,et al.  Cahn–Hilliard and thin film equations with nonlinear mobility as gradient flows in weighted-Wasserstein metrics , 2012, 1201.2367.

[24]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[25]  D. Kinderlehrer,et al.  THE VARIATIONAL FORMULATION OF THE FOKKER-PLANCK EQUATION , 1996 .

[26]  A Wasserstein approach to the numerical solution of the one-dimensional Cahn-Hilliard equation , 2010 .