Geometric Numerical Integration for Peakon b-Family Equations

Geometric Numerical Integration for Peakon b-Family Equations

Year:    2016

Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 24–52

Abstract

In this paper, we study the Camassa-Holm equation and the Degasperis-Procesi equation. The two equations are in the family of integrable peakon equations, and both have very rich geometric properties. Based on these geometric structures, we construct the geometric numerical integrators for simulating their soliton solutions. The Camassa-Holm equation and the Degasperis-Procesi equation have many common properties, however, they also have the significant differences, for example, there exist the shock wave solutions for the Degasperis-Procesi equation. By using the symplectic Fourier pseudo-spectral integrator, we simulate the peakon solutions of the two equations. To illustrate the smooth solitons and shock wave solutions of the DP equation, we use the splitting technique and combine the composition methods. In the numerical experiments, comparisons of these two kinds of methods are presented in terms of accuracy, computational cost and invariants preservation.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.171114.140715a

Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 24–52

Published online:    2016-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    29

Keywords: