doi:10.1016/j.matcom.2007.01.016
Copyright © 2007 IMACS Published by Elsevier Ltd.
Padé numerical method for the Rosenau–Hyman compacton equation
aDepartamento de Lenguajes y Ciencias de la Computación, Universidad de Málaga, Campus de Teatinos, 29071 Málaga, Spain
Available online 21 January 2007.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
Three implicit finite difference methods based on Padé approximations in space are developed for the Rosenau–Hyman K(n,n) equation. The analytical solutions and their invariants are used to assess the accuracy of these methods. Shocks which develop after the interaction of compactons are shown to be independent of the numerical method and its parameters indicating that their origin may not be numerical. The accuracy in long-time integrations of high-order Padé methods is shown.
Keywords: Padé methods; K(n,n) equation; Compactons; Dispersive shocks
Fig. 2. Left plot, shocked zero-mass ripple after the collision of two compactons with c1=2 and c2=1; center plot, zoom in of the shock region in the left plot; and right plot, compactons and shocks generated from a “dilated” compacton with c=1 and dilation equal to 10 (the dashed line shows the initial data). The three graphics for the K(2,2) equation with Method 1,
,
, c0=1 and
=10−4.
Table 1.
Numerical errors of the invariants and the solution for the one-compacton solution of the K(2,2) equation calculated with Methods 1–3, Methods 1–3 and Methods 1–3 with
, L=90,
=0, and c=c0=1, where Err (Ii)=Ii,n(0)−Ii,n(150), for i=1,2, Err (Ii)=max[0,150]|Ii−Ii,n|, for i=3,4, Err (U)=
uc(150)−U(150)
∞, Err (Û)=
ûc(150)−U(150)
∞, uc is the analytical solution, U the numerical one, and ûc is the analytical solution relocated at the position of the numerical one
