Volume 46, Number 3, May-June 2012Special volume in honor of Professor David Gottlieb
|Page(s)||515 - 534|
|Published online||11 January 2012|
Department of Mathematics, North Carolina State
2 Control and Dynamical Systems, California Institute of Technology, Pasadena, 91125 CA, USA
3 Control and Dynamical Systems, California Institute of Technology, Pasadena, 91125 CA, USA
The purpose of this paper is to apply particle methods to the numerical solution of the EPDiff equation. The weak solutions of EPDiff are contact discontinuities that carry momentum so that wavefront interactions represent collisions in which momentum is exchanged. This behavior allows for the description of many rich physical applications, but also introduces difficult numerical challenges. We present a particle method for the EPDiff equation that is well-suited for this class of solutions and for simulating collisions between wavefronts. Discretization by means of the particle method is shown to preserve the basic Hamiltonian, the weak and variational structure of the original problem, and to respect the conservation laws associated with symmetry under the Euclidean group. Numerical results illustrate that the particle method has superior features in both one and two dimensions, and can also be effectively implemented when the initial data of interest lies on a submanifold.
Mathematics Subject Classification: 35708 / 37K10 / 65M25 / 74J35 / 76B15
Key words: Solitons / peakons / integrable Hamiltonian systems / particle methods / weak solutions / variational principle / momentum maps / shallow water and internal waves
© EDP Sciences, SMAI, 2012
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.