Yong Chen a,b,c,qiwang c,d, and Biao Li c,d

Transcription

1 Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation to Construct New Doubly-periodic Solutions of Nonlinear Evolution Equations Yong Chen abc QiWang cd and Biao Li cd a Department of Mathematics Ningbo University Ningbo China b Department of Physics Shanghai Jiao Tong University Shanghai China c Department of Applied Mathematics Dalian University of Technology Dalian China d MM Key Lab Chinese Academy of Sciences Beijing China Reprint requests to Y. Ch.; chenyong18@sjtu.edu.cn. Z. Naturforsch. 59a ; received May A new Jacobi elliptic function rational expansion method is presented by means of a new general ansatz and is very powerful with aid of symbolic computation to uniformly construct more new exact doubly-periodic solutions in terms of rational form Jacobi elliptic function of nonlinear evolution equations NLEEs. We choose a 2+1-dimensional dispersive long wave equation to illustrate the method. As a result we obtain the solutions found by most existing Jacobi elliptic function expansion methods and find other new and more general solutions at the same time. When the modulus of the Jacobi elliptic functions m 1 or 0 the corresponding solitary wave solutions and trigonometric function singly periodic solutions are also found. Key words: 2+1-dimensional Dispersive Long Wave Equation; Jacobi Elliptic Functions; Travelling Wave Solution; Soliton Solution; Periodic Solution. 1. Introduction Many activities aim at finding methods for exact solutions of nonlinear evolutions equations such as Bäcklund transformation Darboux transformation variable separation approach various tanh methods Painleve method generalized hyperbolic-function method homogeneous balance method similarity reduction method and so on [1 10]. For that there are some relations between the elliptic functions e.g. Jacobi elliptic functions and Weierstrass elliptic functions and some nonlinear evolution equations NLEEs. Many authors have paid attention to this subject to know whether elliptic functions can be used to express the exact solutions of more nonlinear equations. Recently Liu et al. [11] presented a Jacobi elliptic function expansion method that used three Jacobi elliptic functions to express exact solutions of some NLEEs. Fan [12] extended the Jacobi elliptic function method to some NLEEs and in particular special-type nonlinear equations for constructing their doubly periodic wave solutions. Such equations cannot be directly dealt with by this method and require some kinds of pre-possessing techniques. Yan [13] further developed an extended Jacobi elliptic function ex- pansion method by using 12 Jacobi elliptic functions. Based on the above idea a new Jacobi elliptic function rational expansion method is presented by means of a new general ansatz and is more powerful than the above Jacobi elliptic function method [12 13] to uniformly construct more new exact doubly-periodic solutions in terms of rational formal Jacobi elliptic functions of NLEEs. For illustration we apply the proposed method to the 2+1-dimensional dispersive long wave equationdlwe which reads u yt + v xx +uu x y = v t + u x +uv x + u xxy = The 2+1-dimensional DLWE 1 was first derived by Boiti et al. [14] as a compatibility for a weak Lax pair. Recently considerable effort has been devoted to the study of this system. In [15] Paquin and Winternitz showed that the symmetry algebra of the 2+1-dimensional DLWE 1 is infinite-dimensional and possesses a Kac-Moody-Virasoro structure. Some special similarity solutions are also given in [15] by using symmetry algebra and classical theoretical analysis. The more general symmetry algebra w symmetry algebra is given in [16]. Lou [17] has given / 04 / $ c 2004 Verlag der Zeitschrift für Naturforschung Tübingen

2 530 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation nine types of the two-dimensional partial differential equation reductions and thirteen types of the ordinary differential equation reductions by means of the direct and nonclassical method. The systems 1 have no Painleve property though they are Lax or IST integrable [18]. More recently Tang et al. [19] by means of the variable separation approach the abundant localized coherent structures of the system 1 have given out. In [20] the possible chaotic and fractal localized structures are revealed for the system 1. This paper is organized as follows: In Sect. 2 we summarize the Jacobi elliptic function rational expansion method. In Sect. 3 we apply the generalized method to 2+1-dimensional dispersive long wave equations and bring out many solutions. Conclusions will be presented in Section Summary of the Jacobi Elliptic Function Rational Expansion Method In the following we outline the main steps of our method: Step 1. For a nonlinear evolution equation with physical fields u i xyt in three variables x y t F i u i u it u ix u iy u itt u ixt u iyt u ixx u iyy u ixy =0 2.1 we use the transformation u i xyt=u i ξ ξ = kx + ly λ t 2.2 where k l and λ are constants to be determined later. Then 2.1 is reduced to the nonlinear ordinary differential equation ODE G i u i u i u i = Step 2. We introduce a new ansatz in terms of a finite Jacobi elliptic function rational expansion: u i ξ =a i0 + m i a ij sn j ξ µsnξ +1 j sn j 1 ξ cnξ + b ij µsnξ +1 j. 2.4 Notice that mi du i dξ = dnξ a ij jsnξ j 1 cnξ b ij snξ j 2 b ij µsnξ j 1 µsnξ +1 j+1 m i dnξ b ij jsnξ + j 2 b ij jsnξ j µsnξ +1 j where snξ cnξdnξnsξcsξdsξ etc. are Jacobi elliptic functions which are double periodic and posses the following properties: 1. Properties of triangular functions: cn 2 ξ + sn 2 ξ = dn 2 ξ + m 2 sn 2 ξ = Derivatives of the Jacobi elliptic functions: sn ξ = cnξ dnξ cn ξ = snξ dnξ 2.7 dn ξ = m 2 snξ cnξ where m is a modulus. The Jacobi-Glaisher functions for elliptic functions can be found in [21 23]. Step 3. The underlying mechanism for a series of fundamental solutions such as polynomial exponential solitary wave rational triangular periodic Jacobi and Weierstrass doubly periodic solutions to occur is that different effects that act to change the wave forms in many nonlinear equations i. e. dispersion dissipation and nonlinearity either separately or in various combinations are able to balance out. We define the degree of u i ξ as D[u i ξ ] = n i which gives rise to the degrees of other expressions as D[u α i ]=n i + α D[u β i uα j s ]=n i β +α + n j s. 2.8 Therefore we can get the value of m i in 2.4. If n i is a nonnegative integer then we first make the transformation u i = ω n i.

3 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation 531 Step 4. Substitute 2.4 into 2.3 along with 2.5 and 2.7 and then set all coefficients of sn i ξ cn j ξ i = 12 ; j = 01 to zero to get an overdetermined system of algebraic equations with respect to λ l µ a i0 a ij and b ij i = 12 ; j = 12 m i. Step 5. Solving the over-determined system of nonlinear algebraic equations by use of Maple we would end up with explicit expressions for λ l µ k a i0 a ij and b ij i = 12 ; j = 12 m i. In this way we can get double periodic solutions with Jacobi elliptic functions. Since: lim snξ = tanhξ lim m 1 cnξ = sechξ m 1 lim dnξ = sechξ m 1 lim snξ = sinξ lim m 0 cnξ = cosξ m 0 lim dnξ = m 0 u 1 and u 2 degenerate respectively as the following form: 1. solitary wave solutions: u i ξ =a i0 + m i a ij tanh j ξ µ tanhξ +1 j tanh j 1 ξ sechξ + b ij µ tanhξ +1 j. 2. triangular function formal solution: u i ξ =a i0 + m i a ij sin j ξ µ sinξ +1 j sin j 1 ξ cosξ + b ij µ sinξ +1 j So the Jacobi elliptic function rational expansion method is more powerful than the method by Liu et al. [11] the method by Fan [12] and the method extended by Yan [13]. The solutions which contain solitary wave solutions singular solitary solutions and triangular function formal solutions can be obtained by the extended method. Remark1: If we replace the Jacobi elliptic functions snξ and cnξ in the ansatz 2.4 with other pairs of Jacobi elliptic functions such as snξ and dnξ ; nsξ and csξ ; nsξ and dsξ ; scξ and ncξ ; dcξ and ncξ ; sdξ and ndξ ; cdξ and ndξ [13 14] it is necessary to point out that the above combinations only require solving the recurrent coefficient relation or derivative relation for the terms of polynomial to close the computation. Therefore other new double periodic wave solutions solitary wave solutions and triangular functional solutions can be obtained for some systems. For simplicity we omit them here. 3. Exact Solutions of the 2+1-dimensional Dispersive Long Wave Equation DLWE Let us consider the 2+1-dimensional dispersive long wave equation DLWE i. e. { uyt + v xx +uu x y = v t + u x +uv x + u xxy = 0. According to the above method to seek travelling wave solutions of 3.1 we make the transformation uxyt=φξ vxyt=σξ ξ = x + ly λt 3.2 where l and λ are constants to be determined later and thus 3.1 become { λ lφ + σ + lφ 2 + lφφ = 0 λσ + φ +φσ + lφ 3.3 = 0. According to Step 1 in Section 2 by balancing φ ξ and σξ φξ in 3.3 and by balancing σ ξ and φξ φ ξ in 3.3 we suppose that 3.3 has the following formal solutions snξ φξ =a 0 + a 1 µsnξ +1 + b 1cnξ µsnξ +1 snξ σξ =A 0 + A 1 µsnξ +1 + B cnξ µsnξ +1 sn 2 ξ + A 2 µsnξ B snξ cnξ 2 µsnξ +1 2 where a 0 a 1 b 1 A 0 A 1 A 2 B 1 B 2 are constants to be determined later. With the aid of Maple substituting 3.4 along with 2.6 and 2.7 into 3.3 yields a set of algebraic equations for sn i ξ cn j ξ i = 01 ; j = 01. Setting the coefficients of these terms sn i ξ cn j ξ

4 532 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation to zero yields a set of over-determined algebraic equations with respect to a 0 a 1 b 1 A 0 A 1 B 1 A 2 B 2 l and λ. By use of the Maple soft package Charsets by Dongming Wang which is based on the Wuelimination method [24] solving the over-determined algebraic equations we get the following results. Case 3. Case 1. Case 2. A 1 = B 1 = µ = 0 a 0 = λ a 1 = ±m b 1 = ±im A 0 = l 1 A 2 = im 2 B 2 = ±ilm b 1 = A 1 = B 1 = B 2 = µ = 0 a 0 = λ a 1 = ±2m A 0 = l 1 + lm 2 A 2 = 2im Case 4. a 1 = A 1 = B 1 = B 2 = µ = 0 a 0 = λ b 1 = ±2im A 0 = l 1 A 2 = 2im Case 5. Case 6. a 0 = 8λ m2 + 2λ ± i 20m m m 2 a 1 = ± m 2 b 1 = ± 4 2 A 0 = 24lm2 16m 2 + l m 2 B 1 = ±il 1 4m 2 20m m 2 A 1 = i2m2 + 3l A 2 = lm l B 2 = ± l 1 4m 2 20m µ = ± i. a 0 = 8m6 λ + 20m 4 λ + 10λ m 2 2λ ± 24m m m 2 m m m m 4 + 5m m a 1 = ± m m 2 m 4m 2 b 1 = ± m m m 2 A 0 = 16m6 + 32lm 6 40m m 4 l + 4lm 2 20m l 44m m 4 + 5m 2 µ = ± 1 A 1 = ± 2m4 + 6m 2 + 1l m m 2 1 4m 4 + 2m m m m 2 A 2 = 3 l8m6 + 16m m 2 16m 4 + 2m B 1 = ± l m m m 2 24m m m 2 m m m m m 4 + 4m 2 1 B 2 = ± l m m m 2 24m m m 2 8m 4 + 2m a 0 = ±µ3 µ+ m 2 + µ 4 µ 2 µ 2 m 2 λ m 2 + µ 4 µ 2 µ 2 m 2 a 1 = ± m 2 + µ 4 µ 2 µ 2 m 2 b 1 = ± m 2 + µ 2 A 0 = lµ4 + 2lµ 2 m 2 µ 2 + m 2 lm 2 µ 2 + m 2 A 1 = lm 2 µ +2lµ 3 lµ A 2 = lm 2 lµ 4 +lµ 2 +lµ 2 m 2 B 1 = ± l m 2 + µ 2 µ 1 + µ 2 m 2 + µ 4 µ 2 µ 2 m 2 B 2 = ±l m 2 + µ 4 µ 2 µ 2 m 2 m 2 + µ

5 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation 533 Case 7. m a 0 = 2 + µ 4 µ 2 µ 2 m 2 µ + µm 2 m 2 + µ 4 µ 2 µ 2 m 2 2 m 2 + µ 4 µ 2 µ 2 m 2 µ 3 m 2 µ 4 + µ 2 + µ 2 m 2 λ m 2 µ 4 λ + µ 2 λ + µ 2 m 2 λ m 2 µ 4 + µ 2 + µ 2 m 2 a 1 = 2 m 2 + µ 4 µ 2 µ 2 m 2 b 1 = B 1 = B 2 = 0A 0 = µ2 m 2 µ 4 + µ 2 m 2 6lµ 2 m 2 + 3lµ 4 + m 4 l 2lµ 6 + 3lm 2 µ 4 + lm 2 m 2 µ 4 + µ 2 + µ 2 m 2 A 1 = 2lµ 2lm 2 µ + 4lµ 3 A 2 = 2lm 2 + 2lµ 2 m 2 2lµ 4 + 2lµ From and Case 1 7 we obtain the following solutions for 3.1. Family 1. From 3.5 we obtain the following rational formal doubly periodic solutions for the DLWE: u 1 xyt=λ ± msnξ ± imcnξ v 1 xyt=l 1 lm 2 sn 2 ξ ± ilm 2 snξ cnξ 3.12 where ξ = x + ly λt l and λ are arbitrary constants. Family 2. From 3.6 we obtain the following rational formal doubly periodic solutions for the DLWE: u 2 xyt=λ ± 2msnξ v 2 xyt=l 1 + lm 2 2lm 2 sn 2 ξ 3.13 where ξ = x + ly λt l and λ are arbitrary constants. Family 3. From 3.7 we obtain the following rational formal doubly periodic solutions for the DLWE: u 3 xyt=λ ± 2imcnξ v 3 xyt=l 1 2lm 2 sn 2 ξ 3.14 where ξ = x + ly λt l and λ are arbitrary constants. Family 4. From 3.8 we obtain the following rational formal doubly periodic solutions for the DLWE: u 4 xyt= 8λ m2 + 2λ ± i 20m m 2 ± v 4 xyt= 24lm2 16m 2 + l m m 2 + 5snξ ±2isnξ +4 ± 1 4m 2 cnξ ±isnξ +2 + i2m2 + 3lsnξ ±2isnξ +4 ± ilcnξ ±2isnξ lm l sn 2 ξ ± 1 2 ± l 1 4m 2 20m 2 + 5snξ cnξ 2 isnξ +1 8 ± isnξ where ξ = x + ly λt l and λ are arbitrary constants. Family 5. From 3.9 we obtain the following rational formal doubly periodic solutions for the DLWE: 24m u 5 xyt=a 0 ± m m 2 snξ m ± m m 2 cnξ m 2 m + 1 ± m 2 1snξ m 1+m 2 m ± 2 +12m 2 1snξ m 2

6 534 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation v 5 xyt=a 0 ± 2m4 + 6m 2 + 1l m m 2 1snξ m 4 + 2m 2 m + 1 ± m 2 1snξ m 2 3l8m m m 2 sn 2 ξ 16m 4 + 2m 2 m + 1 ± 2 +12m 2 1snξ m 2 ± l m m m 2 24m m m 2 m m 2 1cnξ 4m m m 4 + 4m 2 m 1 ± m 2 1snξ m 2 ± l m m m 2 24m m m 2 snξ cnξ m 4 + 2m 2 m + 1 ± 2 +12m 2 1snξ m 2 where ξ = x + ly λt a 0 = 8m6 λ +20m 4 λ +10λ m 2 2λ ± 24m 6 +48m m 2 m 2 +12m m 6 +10m 4 +5m 2 1 A 0 = 16m6 +32lm 6 40m 4 +44m 4 l+4lm 2 20m 2 +4+l l and λ are arbitrary constants. 44m 6 +10m 4 +5m 2 1 Family 6. From 3.10 we obtain the following rational formal doubly periodic solutions for the DLWE: m u 6 xyt=± 2 + µ 4 µ 2 µ 2 m 2 snξ m ± 2 + µ 2 cnξ µsnξ +1 µsnξ +1 v 6 xyt= lµ4 + 2lµ 2 m 2 µ 2 + m 2 lm 2 µ 2 + m 2 + lm2 µ + 2lµ 3 lµsnξ µsnξ ±µ3 µ+ m 2 + µ 4 µ 2 µ 2 m 2 λ m 2 + µ 4 µ 2 µ 2 m 2 ± l m 2 + µ 2 µ 1 + µ 2 cnξ m 2 + µ 4 µ 2 µ 2 m 2 µsnξ +1 + lm2 lµ 4 + lµ 2 + lµ 2 m 2 sn 2 ξ µsnξ ± l m 2 + µ 4 µ 2 µ 2 m 2 m 2 + µ 2 snξ cnξ µsnξ +1 2 where ξ = x + ly λt l and λ are arbitrary constants. Family 7. From 3.11 we obtain the following rational formal doubly periodic solutions for the DLWE: m u 7 xyt=a µ 4 µ 2 µ 2 m 2 snξ µsnξ +1 v 7 xyt=a 0 + 2lµ 2lm2 µ + 4lµ 3 snξ µsnξ lm2 + 2lµ 2 m 2 2lµ 4 + 2lµ 2 sn 2 ξ µsnξ m where ξ = x + ly λt a 0 = 2 +µ 4 µ 2 µ 2 m 2 µ+µm 2 m 2 +µ 4 µ 2 µ 2 m 2 2 m 2 +µ 4 µ 2 µ 2 m 2 µ 3 m 2 µ 4 +µ 2 +µ 2 m 2 A 0 = µ2 m 2 µ 4 +µ 2 m 2 6lµ 2 m 2 +3lµ 4 +m 4 l 2lµ 6 +3lm 2 µ 4 +lm 2 l and λ are arbitrary constants. m 2 µ 4 +µ 2 +µ 2 m 2 Remark 2: 1 The solutions 3.12 reproduce the solution 15 in [22] when A 0 = 2C 1+1m m 2 λ 2 2 The other solutions obtained here to our knowledge are all new families of doubly periodic solutions of the DLWE. 3 Three figures are drawn to illustrate the properties for some solutions.

7 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation 535 Fig. 1. The doubly periodic solutions u 2 and v 2 where l = 1 λ = 1 m = 1 2 t = 0. Fig. 2. The doubly periodic solutions u 6 and v 6 where l = 1 λ = 1 m = 1 2 t = 0 µ = 2. Fig. 3. The doubly periodic solutions u 7 and v 7 where l = 1 λ = 1 m = 1 2 t = 0 µ = 2.

8 536 Y. Chen et al. Jacobi Elliptic Function Rational Expansion Method with Symbolic Computation 4. Summary and Conclusions In this paper we have presented the new Jacobi elliptic function rational expansion method. The method is more powerful than the method proposed recently by Liu [11] Fan [12] and Yan [13]. The 2+1- dimensional dispersive long wave equation DLWE is chosen to illustrate the method such that seven families of new Jacobi elliptic function solutions are obtained. When the modulus m 1 some of these obtained solutions degenerate as solitary wave solutions. The algorithm can also be applied to many nonlinear evolution equations in mathematical physics. Further work about various extensions and improvement of the Jacobi function method is needed to find more general ansätze or the more general subequation. Acknowledgements The work is supported by the Outstanding Youth Foundation No [1] M. J. Ablowitz and P. A. Clarkson Soliton Nonlinear Evolution Equations and Inverse Scattering Cambridge University Press New York [2] V. B. Matveev and M. A. Salle Darbooux Transformation and Soliton Spinger-Verlag Berlin [3] S. Y. Lou and J. Z. Lu J. Phys. A ; S. Y. Lou Phys. Lett. A ; S. Y. Lou and H. Y. Ruan J. Phys. A [4] R. Hirta and J. Satsuma Phys. Lett. A [5] H. W. Tam W. X. Ma X. B. Hu and D. L. Wang J. Phys. Soc. Japan [6] Y. Chen B. Li and H. Q. Zhang Chaos Solitons and Fractals ; Y. Chen Z. Y. Yan and H. Q. Zhang Theor. Math. Phys ; Y. Chen and Y. Zheng IJMPC ; B. Li Y. Chen and H. Q. Zhang Phys. Lett. A [7] M. L. Wang Phys Lett. A [8] E. J. Parkes and B. R. Duffy Comp. Phys. Commun ; E. J. Parkes and B. R. Duffy Phys. Lett. A [9] Y. T. Gao and B. Tian Comput. Phys. Commun [10] E. Fan Computer Physics Communications ; E. Fan Chaos Solitons and Fractals ; E. Fan Physics Letters A [11] S. K. Liu et al. Phys. Lett. A ; S. K. Liu et al. Phys. Lett. A [12] E. Fan Engui Hon and Y. C. Benny Phys. Lett. A ; E. Fan Engui H. Q. Zhang and Jian Phys. Lett. A [13] Z. Y. Yan Chaos Solitons and Fractals ; Z. Y. Yan Chaos Solitons and Fractals ; Z. Y. Yan Computer Physics Communications [14] P. Boiti J. J. P. Leon M. Manna and F. Pempinelli Inverse Problem ; ibid [15] G. Paquin and P. Winternitz Physica D [16] S.-Y. Lou J. Phys. A: Math. Gen [17] S.-Y. Lou Math. Methods Appl. Sci [18] S.-Y. Lou Phys. Lett. A [19] X.-Y. Tang C.-L. Chen and S.-Y. Lou J. Phys. A: Math. Gen. 35 L [20] X.-Y. Tang and S.-Y. Lou Chaos Solitons and Fractals [21] K. Chandrasekharan Elliptic Function Springer- Verlag Berlin [22] D. V. Patrick Elliptic Function and Elliptic Curves At the University Press Cambridge [23] Z. X. Wang and X. J. Xia Special Functions World Scientific Singapore [24] W. Wu Algorithms and Computation Springer-Verlag Berlin 1994 p.1.