![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
The homotopy analysis method is applied to a b-Novikov equation in order to obtain analytic approximations of its peakon solution for various values of b. The results demonstrate that there is an excellent agreement between the approximate solution and the known exact peakon solution of the equation. Moreover, the amplitude of the peakon is approximated and a conservation property of the obtained solution is validated.
1. Introduction
A considerable amount of work has been devoted to finding and studying travelling wave solutions of nonlinear partial differential equations (PDEs), including solitons and peakons. In this paper, analytic approximations of the peakon solution of the b-Novikov equation
(1)
(1) will be determined for various integer values of
. Equation (Equation1
(1)
(1) ) generalizes the well-known Novikov equation
(2)
(2) which is deduced from (Equation1
(1)
(1) ) for b = 3. At the same time, Equation (Equation1
(1)
(1) ) is a special case (for k = 2) of the generalized Camassa-Holm (g-kbCH) equation
(3)
(3) which was studied in [Citation1], as well as (for p = 2, C = 1, A = b + 1) of the four parameter equation
(4)
(4) studied in [Citation2].
Equation (Equation2(2)
(2) ) was derived in [Citation3] and has been extensively studied mostly from a theoretical point of view. For a brief review of the existing results for (Equation2
(2)
(2) ) the reader is referred to [Citation4–6] and the references therein.
In [Citation1], it was proved among other things, that every member of the family of Equation (Equation3(3)
(3) ) has peakon travelling wave solutions, which on the line, have the form
(5)
(5) where
is the velocity of the wave. Obviously, (Equation1
(1)
(1) ) admits the peakon solution
(6)
(6) Thus, peakon solutions are considered as weak solutions to the PDE under consideration. However, they are orbitally stable (see [Citation7] for the Novikov equation). Moreover, the peakon (Equation5
(5)
(5) ) belongs to any Sobolev space
with
, [Citation1]. In addition, peakon solutions can be considered as “reminiscent of Stokes waves” [Citation8] and thus are extremely important on the study of water wave phenomena.
Regarding Equation (Equation1(1)
(1) ) there exist very few papers focusing on theoretical results. In [Citation9], the local well posedness for the strong solutions to the Cauchy problem of (Equation1
(1)
(1) ) was studied. In addition, it was proved that (Equation6
(6)
(6) ) is a global weak solution to (Equation1
(1)
(1) ). The local well posedness to the Cauchy problem of (Equation1
(1)
(1) ) in critical Besov spaces, was also studied in [Citation10]. Most recently in [Citation11], the instability and nonuniqueness for the initial value problem of (Equation1
(1)
(1) ), when the initial data belong to Sobolev spaces
with
, was studied both on the line and the circle, proving that for b>2, (Equation1
(1)
(1) ) is ill-posed for
. This result was accomplished by constructing appropriate two-peakon solutions with arbitrary small initial size data that collide in arbitrary small time. Moreover, a complete group classification of (Equation1
(1)
(1) ) was carried out in [Citation12], which in turn led to the construction of solutions to (Equation1
(1)
(1) ) and the establishment of some conservation laws. Studying Equation (Equation4
(4)
(4) ) was the subject of [Citation2]. More precisely, low-order conservation laws were classified, the admittance of peakon and multi-peakon solutions to (Equation4
(4)
(4) ) was rigorously proved and their dynamics was investigated.
Another equation strongly connected with (Equation1(1)
(1) ) is the b-equation
(7)
(7) which for b = 2 reduces to the well-known Camassa Holm equation and can be deduced from (Equation3
(3)
(3) ) for k = 1. For a brief review of the results concerning (Equation7
(7)
(7) ), the reader is referred to [Citation9] and the references therein. Equation (Equation7
(7)
(7) ) was also numerically studied in [Citation13,Citation14]. However, to the best of the authors' knowledge, there aren't any numerical results available for (Equation1
(1)
(1) ). In addition, it is worth mentioning that there are very few papers where peakons are numerically computed, compared to papers concerned with smooth solitons, and in these few papers more elaborate techniques are utilized. The reader may refer to [Citation15–17].
Since Equation (Equation1(1)
(1) ) is highly nonlinear, finding its solutions analytically is quite difficult. However, analytic approximations of it may be obtained using approximation techniques. One analytical technique which has been proved extremely successful for the evaluation of peakon solutions, is the Homotopy Analysis Method (HAM) which is thoroughly described in books [Citation18,Citation19]. Indeed, the HAM was applied in [Citation20,Citation21] for determining peakon solutions to the Degasperis-Procesi and the Camassa-Holm equation, respectively. In both cases the amplitude of the peakon was also approximated. More recently in [Citation4], the HAM was implemented in order to deduce the peakon and antipeakon solutions to equation
(8)
(8) which is the Novikov equation including the term
.
Applying HAM to problems of nonlinear water waves is one of the five challenging problems included in [Citation22, p. 1–33]. In the present paper, the work of [Citation20,Citation21] is continued for Equation (Equation1(1)
(1) ), which has a cubic nonlinearity, whereas both the Camassa-Holm and the Degasperis-Processi equation have a quadratic nonlinearity.
The rest of the paper is organized as follows: In Section 2, the mathematical formulation of the problem, as well as a thorough description of HAM applied to Equation (Equation1(1)
(1) ) are given. In addition, the peakon solution in the special case b = 2 is analytically obtained. In Section 3, a discussion of the obtained results is given, including comparisons with the exact solution (Equation6
(6)
(6) ), as well as estimates of the amplitude of the peakon and the value of the first derivative at the crest. In addition, the results indicate that the quantity
, which is conserved by solutions of the Novikov equation, is also conserved by solutions of (Equation1
(1)
(1) ), which was already stated in [Citation2,Citation12]. Finally, in Section 4 conclusions are presented.
2. Mathematical formulation and implementation of HAM
2.1. Mathematical formulation
It is well-known that solitons are travelling wave solutions to a PDE of the form
(9)
(9) where
is the velocity of the wave, satisfying the boundary conditions
(10)
(10) In order to find the peakon solution to Equation (Equation1
(1)
(1) ) and on the same time incorporate the amplitude a of the wave, instead of (Equation9
(9)
(9) ) the transformation
(11)
(11) will be used. Thus, peakons of (Equation1
(1)
(1) ) are obtained as solutions to the ordinary differential equation (ODE)
(12)
(12) Although
, it suffices to study (Equation12
(12)
(12) ) only for
, since the simple change of variables
(13)
(13) leaves (Equation12
(12)
(12) ) unchanged. Moreover, it is easy to see that if g is a solution to (Equation12
(12)
(12) ), so is
, something already depicted in [Citation9] where it was proved that not only (Equation6
(6)
(6) ), but also
is a global weak solution to (Equation1
(1)
(1) ).
Due to the form of the known solution (Equation6(6)
(6) ), it is natural to change the dependent variable of (Equation12
(12)
(12) ) by setting
(14)
(14) and study instead the ODE
(15)
(15) from which c is eliminated. Of course, the
denotes differentiation with respect to η.
Thus, in order to obtain peakon solutions of (Equation1(1)
(1) ) above the horizontal axis, it suffices to calculate the solution to (Equation15
(15)
(15) ) subject to the conditions:
(16)
(16) As it will be shown in the next paragraph, the other two boundary conditions
will be automatically satisfied. It should be noted that the first derivative of the peakon at the crest is not continuous and thus an initial condition on
cannot be utilized in this framework. This makes the estimation of
especially important.
Before describing the implementation of HAM in the next paragraphs, it is worth mentioning that in the special case b = 2, Equation (Equation15(15)
(15) ) is analytically solvable. Indeed, for b = 2, (Equation15
(15)
(15) ) becomes:
which after one integration with respect to η and taking into account that w,
,
, as
, gives
(17)
(17) Since solitons are non constant solutions to a PDE, (Equation17
(17)
(17) ) reduces to
the general solution of which is
(18)
(18) where
,
are arbitrary constants. Due to the boundary conditions at
, (Equation18
(18)
(18) ) finally takes the form
where
is an arbitrary constant, which is of the form (Equation6
(6)
(6) ), as expected. Thus, in the rest of the paper only
will be considered. It should be noted that a rigorous proof of the existence of peakon (and multi-peakon) solutions to (Equation4
(4)
(4) ) is given in [Citation2].
2.2. Short description of HAM
The Homotopy Analysis Method proposed by Liao (see [Citation18,Citation19]) which will be implemented for solving the boundary value problem (Equation15(15)
(15) )–(Equation16
(16)
(16) ), is based on the notion of homotopy. If
is an initial guess of the solution
of (Equation15
(15)
(15) ), satisfying conditions (Equation16
(16)
(16) ) and
is a parameter called embedding parameter, the HAM is based on a kind of continuous mappings
,
, such that, as q increases from 0 to 1,
varies from
to the exact solution
and
, for q = 1 gives the amplitude a.
In order to ensure this, the following suitable homotopy is constructed
(19)
(19) where
is an auxiliary linear operator with the property
(20)
(20) ℏ is a nonzero auxiliary parameter called convergence control parameter,
is a nonzero auxiliary function and
is a nonlinear operator. Homotopy (Equation19
(19)
(19) ) has the following two important properties:
If
then
(21)
(21)
If
then
(22)
(22)
Thus, if the nonlinear operator is suitably chosen in such a way so that the solution of (Equation22
(22)
(22) ) is
, then the solution
of the family of equations
(23)
(23) will vary from
(for q = 0) to
(for q = 1). Equation (Equation23
(23)
(23) ) is usually called zero-order deformation equation and due to conditions (Equation16
(16)
(16) ) it should be accompanied by conditions
(24)
(24) By expanding
and
in Taylor series with respect to q around 0 and using the notations
(25)
(25)
(26)
(26) these series take the form
(27)
(27)
(28)
(28) Differentiating (Equation23
(23)
(23) ) and (Equation24
(24)
(24) ) m times with respect to q, setting q = 0 and finally dividing with
, the linear equations, with respect to
, are obtained
(29)
(29) which are called mth-order deformation equations. In (Equation29
(29)
(29) ),
is defined by
(30)
(30) and
(31)
(31) Equation (Equation29
(29)
(29) ) are accompanied by suitable conditions which will be given in the next subsection (see relation (Equation39
(39)
(39) ) and the discussion therein). By solving the higher-order deformation equations, the
are found. In order to find the
, specific algebraic equations are solved which are connected with the conditions accompanying (Equation29
(29)
(29) ) and will also be described in the next subsection.
Once and
are determined, the solution
and the amplitude a are also determined, provided that the series (Equation27
(27)
(27) ) and (Equation28
(28)
(28) ) converge at q = 1, by the relations
(32)
(32) and
(33)
(33)
2.3. Implementation of HAM
In order to implement HAM for the solution to (Equation15(15)
(15) )–(Equation16
(16)
(16) ), suitable choices should be made for
,
,
,
and ℏ. Moreover, an appropriate set, called the set of base functions, describing the solution to the problem should also be chosen. For problems regarding peakon solutions to PDEs, the set
(34)
(34) has been proven very useful (see [Citation20,Citation21]). According to HAM, if
can be expressed as a linear combination of members of
and
,
and
are chosen in such a way that
are also expressed as linear combinations of members of
, then the solution to (Equation15
(15)
(15) )–(Equation16
(16)
(16) ) can be represented in a series of the form
(35)
(35) where
are coefficients to be determined. This is called the rule of solution expression.
For the problem discussed in the present paper, the linear operator is chosen to be
(36)
(36) whereas the initial approximation
is taken to be
(37)
(37) where ε is a parameter to be determined later (see Section 3). Obviously (Equation37
(37)
(37) ) is of the form (Equation35
(35)
(35) ) and satisfies the conditions (Equation16
(16)
(16) ). It should be emphasized that (Equation37
(37)
(37) ) was also used in [Citation20,Citation21] as an initial approximation. However, operator (Equation36
(36)
(36) ) differs from the corresponding operators used in the aforementioned papers.
The auxiliary function is usually taken identically equal to 1 and this also stands here. Finally,
is chosen to be
(38)
(38) Using (Equation38
(38)
(38) ) and definition (Equation30
(30)
(30) ), the term
is found to be
Since
defined by (Equation37
(37)
(37) ) satisfies conditions (Equation16
(16)
(16) ) it suffices to consider the boundary conditions
(39)
(39) to accompany (Equation29
(29)
(29) ). In this way, the obtained solution given in the form (Equation32
(32)
(32) ), also satisfies the required conditions (Equation16
(16)
(16) ).
Thus, beginning with and solving the boundary value problems (Equation29
(29)
(29) ), (Equation39
(39)
(39) ) several
can be found. The general solution of (Equation29
(29)
(29) ) with
given by (Equation36
(36)
(36) ), has the form
(40)
(40) where
are arbitrary constants and
is a particular solution to (Equation29
(29)
(29) ). According to the boundary condition (Equation39
(39)
(39) ) at infinity and the rule of solution expression, the constants
must be zero. The unknown
is obtained by solving the linear algebraic equation
(41)
(41) which is deduced from the first condition of (Equation39
(39)
(39) ).
Regarding the convergence of the method, the following proposition holds.
Theorem 2.1
If the series (Equation32(32)
(32) ), (Equation33
(33)
(33) ) converge, then (Equation32
(32)
(32) ) is a solution to (Equation15
(15)
(15) )–(Equation16
(16)
(16) ).
The proof of this theorem is omitted since it is similar to the proof of the general convergence theorem given in [Citation18]. In practice, M terms of the homotopy series are calculated, in order to obtain an adequate approximation of and α. The corresponding series
(42)
(42)
(43)
(43) are called approximate solutions of Mth order, where
.
The convergence of the series at q = 1 can be guaranteed if ℏ is appropriately chosen. One way to choose ℏ, is to use the so called ℏ-curves, which are graphs of specific, usually important, quantities associated with the solution to the problem, against ℏ (see [Citation18]). For the problem discussed in this paper, the ℏ-curves will be the graphs of the amplitude α. If (Equation33(33)
(33) ) is convergent, a horizontal line segment will appear in the ℏ-curves and if a value of ℏ is chosen from this region, it is quite certain that (Equation32
(32)
(32) ) converges.
Another way to choose ℏ, is to find the optimum ℏ for which the ‘discrete squared residual’ error defined by
(44)
(44) is minimized [Citation19], where
,
and N = 20. This approach will be used in the present paper.
3. Results and discussion
In this section, analytic approximations of the solution to (Equation15(15)
(15) )–(Equation16
(16)
(16) ) will be determined via HAM, for some values of b. The implementation of HAM has been carried out via Mathematica, following the outline given in Sections 2.2 and 2.3.
To begin with, the case b = 3 will be considered, since then Equation (Equation1(1)
(1) ) becomes the Novikov Equation (Equation2
(2)
(2) ). As it is obvious from the discussion in sections 2.2 and 2.3, both
and
,
depend on the convergence control parameter ℏ and the parameter ε. By appropriately choosing these two parameters, the series (Equation32
(32)
(32) ) and (Equation33
(33)
(33) ) converge. For a given ℏ, the influence of ε on the convergence of the amplitude a can be investigated by plotting the curve of a versus ε. This is depicted in Figure , where the fourth order approximation of the solution was utilized for
,
and
. Clearly the same can be done for various values of ℏ.
It is obvious from Figure that which agrees with the expected value and a valid region for ε is roughly between
and 2/5. For the rest of the paper and for all the values of b that have been taken into consideration, the value
will be used, since several simulations pointed out this value as reasonable.
After having fixed ε, the next step is to determine an appropriate value for ℏ. This was done by minimizing the errors defined by (Equation44(44)
(44) ). The second, third and fifth column of Table , show the optimum values for ℏ, the corresponding errors and the values of a in the case where b = 3, for several values of the order M from 8 up to 25. As the order increases, the error reduces as expected and a converges to the rounded value 0.974.
Table 1. and corresponding errors for
and the optimum ℏ.
As pointed out in Section 2.1, the first derivative of the peakon at the crest is not continuous and thus an initial condition on cannot be utilized in the framework of HAM. This makes the estimation of
especially important. The Mth order approximation of
is given by the value of
depicted at the fourth column of Table , showing a very good agreement between the calculated by HAM value and the expected one.
Another important quantity is the integral , i.e. the
norm of the solution of (Equation1
(1)
(1) ), which is a conserved quantity as shown in [Citation2,Citation12]. This conserved quantity can be considered as conservation of energy for the peakon solution of (Equation1
(1)
(1) ). Since the peakon decays exponentially to zero, it suffices to calculate the integral
instead of I. The Mth order approximation of
is given by the value of
depicted at the sixth column of Table , showing again a very good agreement between the calculated by HAM values and the expected ones.
Moreover, in Figure the ℏ-curve for a is given for the 25th order approximation of the solution. It is evident that the values for the optimum ℏ are in accordance with this ℏ-curve. The graph of the 25th order approximation of the solution via HAM for the corresponding optimum ℏ is given in Figure , against the graph of the exact solution (Equation6(6)
(6) ) for c = 1. It is observed that there is an excellent agreement between the two.
Figure 3. Comparison of the HAM solution in the case b = 3 for M = 25 with the exact solution for c = 1.
![Figure 3. Comparison of the HAM solution in the case b = 3 for M = 25 with the exact solution for c = 1.](/cms/asset/d4d73ac2-223b-44ca-a2ef-ff551ba291d3/gipe_a_2099850_f0003_ob.jpg)
The same analysis has been carried out by the authors for various values of b, with results analogous to the obtained results for b = 3. Indicatively, tables similar to Table are given for the values b = 4, b = 5 and b = 7 (see Tables ). The results indicate that although HAM is successful in all four cases, the errors increase as b increases albeit they still remain small. This could possibly be overcome by changing ε or even .
Table 2. and corresponding errors for
and the optimum ℏ.
Table 3. and corresponding errors for
and the optimum ℏ.
Table 4. and corresponding errors for
and the optimum ℏ.
4. Conclusions
In this paper, the HAM was applied in order to analytically obtain several approximations of the peakon solutions to the b-Novikov Equation (Equation1(1)
(1) ). Four values of b were taken into consideration, including b = 3 which corresponds to the well known Novikov equation. The obtained results are in excellent agreement with the exact ones.
This work complements earlier works of [Citation20,Citation21], where peakon solutions were studied in a similar way for the Degasperis-Processi and Camassa-Holm equations, which include quadratic nonlinearities. Since (Equation1(1)
(1) ) is a special case of the more general Equation (Equation3
(3)
(3) ), HAM provides an excellent tool for further investigation of peakon solutions to not only (Equation3
(3)
(3) ), but also to other nonlinear PDEs.
Acknowledgments
The authors would like to thank the reviewers for their valuable comments which improved the presentation of the paper and also strengthened the reported findings.
Disclosure statement
The authors report that there are no competing interests to declare.
References
- Grayshan K, Himonas AA. Equations with peakon traveling wave solutions. Adv Dyn Syst Appl. 2013;8:217–232.
- Anco SC, da Silva PL, Freire IL. A family of wave-breaking equations generalizing the Camassa-Holm and Novikov equations. J Math Phys. 2015;56(9):Article ID 091506. 21.
- Novikov V. Generalizations of the Camassa-Holm equation. J Phys A. 2009;42(34):Article ID 342002.
- Efstathiou AG, Petropoulou EN. Peakons of Novikov equation via the homotopy analysis method. Symmetry. 2021;13(5):738.
- Himonas AA, Holliman C, Kenig C. Construction of 2-peakon solutions and ill-posedness for the Novikov equation. SIAM J Math Anal. 2018;50(3):2968–3006.
- Lundmark H, Shuaib B. Ghostpeakons and characteristic curves for the Camassa-Holm, Degasperis-Procesi and Novikov equations. SIGMA Symmetry Integrab Geom Methods Appl. 2019;15:017.
- Liu X, Liu Y, Qu C. Stability of peakons for the Novikov equation. J Math Pures Appl. 2014;101(2):172–187.
- Hone ANW, Lundmark H, Szmigielski J. Explicit multipeakon solutions of Novikov's cubically nonlinear integrable Camassa-Holm type equation. Dyn Partial Differ Equ. 2009;6(3):253–289.
- Mi Y, Mu C. On the Cauchy problem for the modified Novikov equation with peakon solutions. J Differ Equ. 2013;254(3):961–982.
- Zhou S, Chen R. A few remarks on the generalized Novikov equation. J Inequal Appl. 2013;2013(1):560.
- Himonas AA, Holliman C. Instability and nonuniqueness for the b-Novikov equation. J Nonlinear Sci. 2022;32(4), Paper 46 (29 pp). doi:10.1007/s00332-022-09798-6.
- da Silva PL, Freire IL. On the group analysis of a modified Novikov equation. In Cojocaru MG, Kotsireas IS, Makarov RN, et al. editors, Interdisciplinary topics in applied mathematics, modeling and computational science, Springer Proceedings in Mathematics & Statistics, Vol. 117. 2015. p. 161–166.
- Holm DD, Staley MF. Wave structure and nonlinear balances in a family of evolutionary PDEs. SIAM J Appl Dyn Syst. 2003;2(3):323–380.
- Holm DD, Staley MF. Nonlinear balance and exchange of stability in dynamics of solitons, peakons, rambs/cliffs and leftons in a 1−−1 nonlinear evolutionary PDE. Phys Lett A. 2003;308(5–6):437–444.
- Artebrand R, Schroll HJ. Numerical simulation of Camassa–Holm peakons by adaptive upwinding. Appl Numer Math. 2006;56(5):695–711.
- Kalisch H, Lenells J. Numerical study of traveling-wave solutions for the Camassa–Holm equation. Chaos Solitons Fract. 2005;25(2):287–298.
- Xu Y, Shu C-W. Local discontinuous Galerkin method for the Camassa–Holm equation. SIAM J Numer Anal. 2006;46(4):1998–2021.
- Liao SJ. Beyond perturbation: introduction to the homotopy analysis method. Boca Raton: Chapman & Hall/CRC, CRC Press LLC; 2004.
- Liao SJ. Homotopy analysis method in nonlinear differential equations. Heidelberg: Springer & Higher Education Press; 2012.
- Abbasbandy S, Parkes EJ. Solitary-wave solutions of the Degasperis–Procesi equation by means of the homotopy analysis method. Int J Comput Math. 2010;87(10):2303–2313.
- Wu W, Liao SJ. Solving solitary waves with discontinuity by means of the homotopy analysis method. Chaos Solitons Fract. 2005;26(1):177–185.
- Liao SJ. Chance and challenge: a brief review of homotopy analysis method. In Liao S, editor. Advances in the homotopy analysis method. Hackensack, NJ: World Sci. Publ.; 2014.