Integration of the EPDiff equation by particle methods

Alina Chertock; Philip Du Toit; Jerrold Eldon Marsden

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2012)

  • Volume: 46, Issue: 3, page 515-534
  • ISSN: 0764-583X

Abstract

top
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.

How to cite

top

Chertock, Alina, Toit, Philip Du, and Marsden, Jerrold Eldon. "Integration of the EPDiff equation by particle methods." ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique 46.3 (2012): 515-534. <http://eudml.org/doc/273089>.

@article{Chertock2012,
abstract = {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.},
author = {Chertock, Alina, Toit, Philip Du, Marsden, Jerrold Eldon},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {solitons; peakons; integrable hamiltonian systems; particle methods; weak solutions; variational principle; momentum maps; shallow water and internal waves; integrable Hamiltonian systems; Euler-Poincaré differential equation; numerical results},
language = {eng},
number = {3},
pages = {515-534},
publisher = {EDP-Sciences},
title = {Integration of the EPDiff equation by particle methods},
url = {http://eudml.org/doc/273089},
volume = {46},
year = {2012},
}

TY - JOUR
AU - Chertock, Alina
AU - Toit, Philip Du
AU - Marsden, Jerrold Eldon
TI - Integration of the EPDiff equation by particle methods
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2012
PB - EDP-Sciences
VL - 46
IS - 3
SP - 515
EP - 534
AB - 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.
LA - eng
KW - solitons; peakons; integrable hamiltonian systems; particle methods; weak solutions; variational principle; momentum maps; shallow water and internal waves; integrable Hamiltonian systems; Euler-Poincaré differential equation; numerical results
UR - http://eudml.org/doc/273089
ER -

References

top
  1. [1] M.S. Alber, R. Camassa, Y.N. Fedorov, D.D. Holm and J.E. Marsden, The complex geometry of weak piecewise smooth solutions of integrable nonlinear PDE’s of shallow water and Dym type. Commun. Math. Phys.221 (2001) 197–227. Zbl1001.37062MR1846907
  2. [2] R. Artebrant and H.J. Schroll, Numerical simulation of Camassa-Holm peakons by adaptive upwinding. Appl. Numer. Math.56 (2006) 695–711. Zbl1156.65313MR2211502
  3. [3] R. Beals, D.H. Sattinger and J. Szmigielski, Peakon-antipeakon interaction. J. Nonlin. Math. Phys. 8 (2001) 23–27; Nonlinear evolution equations and dynamical systems, Kolimbary (1999). Zbl0977.35106MR1821503
  4. [4] R. Camassa and D.D. Holm, An integrable shallow water equation with peaked solitons. Phys. Rev. Lett.71 (1993) 1661–1664. Zbl0972.35521MR1234453
  5. [5] R. Camassa, D.D. Holm and J.M. Hyman, A new integrable shallow water equation. Adv. Appl. Mech.31 (1994) 1–33. Zbl0808.76011
  6. [6] R. Camassa, J. Huang and L. Lee, On a completely integrable numerical scheme for a nonlinear shallow-water wave equation. J. Nonlin. Math. Phys.12 (2005) 146–162. MR2117177
  7. [7] R. Camassa, J. Huang and L. Lee, Integral and integrable algorithms for a nonlinear shallow-water wave equation. J. Comput. Phys.216 (2006) 547–572. Zbl1220.76016MR2235383
  8. [8] S. Chen, C. Foias, D.D. Holm, E. Olson, E.S. Titi and S. Wynne, Camassa-Holm equations as a closure model for turbulent channel and pipe flow. Phys. Rev. Lett.81 (1998) 5338–5341. Zbl1042.76525MR1745983
  9. [9] A. Chertock and A. Kurganov, On a practical implementation of particle methods. Appl. Numer. Math.56 (2006) 1418–1431. Zbl1103.65103MR2245465
  10. [10] A. Chertock and D. Levy, Particle methods for dispersive equations. J. Comput. Phys.171 (2001) 708–730. Zbl0991.65008MR1848732
  11. [11] A. Chertock and D. Levy, A particle method for the KdV equation. J. Sci. Comput.17 (2002) 491–499. Zbl1001.76079MR1910746
  12. [12] A.J. Chorin, Numerical study of slightly viscous flow. J. Fluid Mech.57 (1973) 785–796. MR395483
  13. [13] G.M. Coclite, K.H. Karlsen and N.H. Risebro, A convergent finite difference scheme for the Camassa-Holm equation with general H1 initial data. SIAM J. Numer. Anal.46 (2008) 1554–1579. Zbl1172.35310MR2391006
  14. [14] A. Cohen and B. Perthame, Optimal approximations of transport equations by particle and pseudoparticle methods. SIAM J. Math. Anal. 32 (2000) 616–636 (electronic). Zbl0972.65058MR1786160
  15. [15] G.-H. Cottet and P.D. Koumoutsakos, Vortex methods. Cambridge University Press, Cambridge (2000). Zbl0953.76001MR1755095
  16. [16] G.-H. Cottet and S. Mas-Gallic, A particle method to solve transport-diffusion equations, Part 1 : the linear case. Tech. Report 115, Ecole Polytechnique, Palaiseau, France (1983). Zbl0678.35077
  17. [17] G.-H. Cottet and S. Mas-Gallic, A particle method to solve the Navier-Stokes system. Numer. Math.57 (1990) 805–827. Zbl0707.76029MR1065526
  18. [18] P. Degond and S. Mas-Gallic, The weighted particle method for convection-diffusion equations. I. The case of an isotropic viscosity. Math. Comput. 53 (1989) 485–507. Zbl0676.65121MR983559
  19. [19] P. Degond and S. Mas-Gallic, The weighted particle method for convection-diffusion equations. II. The anisotropic case. Math. Comput. 53 (1989) 509–525. Zbl0676.65122MR983560
  20. [20] P. Degond and F.-J. Mustieles, A deterministic approximation of diffusion equations using particles. SIAM J. Sci. Statist. Comput.11 (1990) 293–310. Zbl0713.65090MR1037516
  21. [21] S. Gottlieb, C.-W. Shu and E. Tadmor, High order time discretization methods with the strong stability property. SIAM Rev.43 (2001) 89–112. Zbl0967.65098MR1854647
  22. [22] O.-H. Hald, Convergence of vortex methods, Vortex methods and vortex motion. SIAM, Philadelphia, PA (1991) 33–58. MR1095603
  23. [23] A.N. Hirani, J.E. Marsden and J. Arvo, Averaged Template Matching Equations, EMMCVPR, Lecture Notes in Computer Science 2134. Springer (2001) 528–543. Zbl1001.68646
  24. [24] H. Holden and X. Raynaud, Convergence of a finite difference scheme for the Camassa-Holm equation. SIAM J. Numer. Anal. 44 (2006) 1655–1680 (electronic). Zbl1122.76065MR2257121
  25. [25] H. Holden and X. Raynaud, A convergent numerical scheme for the Camassa-Holm equation based on multipeakons. Discrete Contin. Dyn. Syst.14 (2006) 505–523. Zbl1111.35061MR2171724
  26. [26] D.D. Holm and J.E. Marsden, Momentum maps and measure-valued solutions (peakons, filaments, and sheets) for the EPDiff equation, The breadth of symplectic and Poisson geometry, Progr. Math. 232. Birkhäuser Boston, Boston, MA (2005) 203–235. MR2103008
  27. [27] D.D. Holm and M.F. Staley, Wave structure and nonlinear balances in a family of evolutionary PDEs. SIAM J. Appl. Dyn. Syst. 2 (2003) 323–380 (electronic). Zbl1088.76531MR2031278
  28. [28] D.D. Holm and M.F. Staley, Interaction dynamics of singular wave fronts, under “Recent Papers” at http://cnls.lanl.gov/~staley/. 
  29. [29] D.D Holm, J.T. Ratnanather, A. Trouvé and L. Younes, Soliton dynamics in computational anatomy. NeuroImage 23 (2004) S170–S178. 
  30. [30] H.-P. Kruse, J. Scheurle and W. Du, A two-dimensional version of the Camassa-Holm equation, Symmetry and perturbation theory. World Sci. Publ., Cala Gonone, River Edge, NJ (2001) 120–127. Zbl1061.35114MR1875473
  31. [31] A.K. Liu, Y.S. Chang, M.-K. Hsu and N.K. Liang, Evolution of nonlinear internal waves in the east and south China Sea. J. Geophys. Res.103 (1998) 7995–8008. 
  32. [32] J.E. Marsden and T.S. Ratiu, Introduction to mechanics and symmetry, Texts in Applied Mathematics 17, 2nd edition. Springer-Verlag, New York (1999). Zbl0933.70003MR1723696
  33. [33] R.I. McLachlan and P. Atela, The accuracy of symplectic integrators. Nonlinearity5 (1992) 541–562. Zbl0747.58032MR1158385
  34. [34] R. McLachlan and S. Marsland, N-particle dynamics of the Euler equations for planar diffeomorfism. Dyn. Syst.22 (2007) 269–290. Zbl1149.37038MR2354966
  35. [35] P.-A. Raviart, An analysis of particle methods. Numerical methods in fluid dynamics (Como, 1983), Lecture Notes in Math. 1127. Springer, Berlin (1985) 243–324. Zbl0598.76003MR802214
  36. [36] G.D. Rocca, M.C. Lombardo, M. Sammartino and V. Sciacca, Singularity tracking for Camassa-Holm and Prandtl’s equations. Appl. Numer. Math.56 (2006) 1108–1122. Zbl1096.76036MR2234843
  37. [37] S.F. Singer, Symmetry in mechanics : A gentle, modern introduction. Birkhäuser Boston Inc., Boston, MA (2001). Zbl0970.70003MR1816059
  38. [38] Y. Xu and C.-W. Shu, A local discontinuous Galerkin method for the Camassa-Holm equation. SIAM J. Numer. Anal.46 (2008) 1998–2021. Zbl1173.65063MR2399405

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.