Skip to content
Publicly Available Published by De Gruyter September 24, 2015

Exact Travelling Wave Solutions of the Nonlinear Evolution Equations by Auxiliary Equation Method

  • Melike Kaplan EMAIL logo , Arzu Akbulut and Ahmet Bekir

Abstract

The auxiliary equation method presents wide applicability to handling nonlinear wave equations. In this article, we establish new exact travelling wave solutions of the nonlinear Zoomeron equation, coupled Higgs equation, and equal width wave equation. The travelling wave solutions are expressed by the hyperbolic functions, trigonometric functions, and rational functions. It is shown that the proposed method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering. Throughout the article, all calculations are made with the aid of the Maple packet program.

1 Introduction

Nonlinear evolution equations (NLEEs) are encountered in various fields of mathematics, physics, chemistry, biology, engineering, and numerous applications. Exact solutions of NLEEs play a vital role in the proper understanding of qualitative features of many phenomena and processes in various areas of natural science. Thus, many effective and powerful methods have been established and improved. Among these are inverse scattering method [1], Hirota bilinear method [2], Jacobi elliptic function method [3, 4], Darboux and Backlünd transform [5], Wronskian technique [6], tanh method [7], extended tanh method [8], sine-cosine method [9], homogeneous balance method [10], F-expansion method [11], variational iteration method [12], exp-function method [13], (G′/G) -expansion method [14, 15], modified simple equation method [16], transformed rational function method [17], mapping method [18], first integral method [19, 20], trial equation method [21], and so on.

The adopted algorithm was analysed earlier in a nonlinear mechanics journal [22]. More generally, Frobenius decompositions [23] and the refined invariant subspaces method [24] present more diverse exact solutions to nonlinear partial differential equations (PDEs). Also, integrable equations possess Hirota bilinear forms, which take a crucial role in presenting exact solutions. Recently, the Hirota bilinear derivatives were generalised to much wider cases [25], and the corresponding nonlinear equations have resonant soliton solutions.

This article investigates the applicability and effectiveness of the auxiliary equation method on nonlinear evolution equations and nonlinear evolution equation systems. The article is organised as follows. In Section 2, the mentioned method is described. In Section 3, we illustrate this method in detail with the Zoomeron equation, coupled Higgs equation, and equal width wave (EW) equation. Finally, some conclusions are given.

2 Auxiliary Equation Method

The auxiliary equation method is a direct method to find more and new multiple travelling wave solutions for nonlinear PDEs. This method was defined by Sirendaoreji [26, 27] and lately Abdoulkary et al. obtained exact travelling wave solutions of the nonlinear Schrödinger equation [28].

The main steps of using this method are summarised as follows:

Step 1: We consider the following general nonlinear partial differential of the type

(1)P(u,ut,ux,utt,uxt,uxx,)=0, (1)

where u is an unknown function and P is a polynomial of u and its partial derivatives. The travelling wave wariable

u(x,t)=U(ξ),ξ=x±dy±wt

carries (1) to a nonlinear ordinary differential equation

(2)Q(U,U,U,)=0, (2)

where the prime denotes the derivation with respect to ξ.

Step 2: Suppose that the solution of the ordinary differential equation (ODE) (2) can be expressed as

(3)U(ξ)=i=02NaiFi(ξ) (3)

where ai (i=0, …, 2N) are constants to be determined later, and the positive integer N can be determined by using the homogeneous balance method between the highest-order derivatives and the nonlinear terms appearing in ODE (2). Here, F(ξ) satisfies the following variable separated ordinary differential equation

(4)F2(ξ)=aF2(ξ)+bF4(ξ)+cF6(ξ) (4)

where a, b, and c are parameters to be determined.

Step 3: Substituting ansatz (3) along with (4) into (2) and equating the coefficients of all powers of F(ξ) to zero yields a set of algebraic equations for unknowns a, b, c, d, ai (i=0, …, 2N), and w. We solve the set of algebraic equations with the aid of Maple and substitute the obtained solutions in this step back into (2) so as to obtain the exact travelling wave solutions for (1).

3 Applications

3.1 Zoomeron Equation

In this subsection, we implement the auxiliary equation method to solve the Zoomeron equation in the form

(5)(uxyu)tt(uxyu)xx+2(u2)xt=0, (5)

where u(x, y, t) is the amplitude of the relative wave mood. This equation is one of incognito evolution equation. The equation was introduced by Calogero and Degasperis [29]. Using the wave wariable

(6)u(x,y,t)=U(ξ),ξ=x+dywt, (6)

Equating (5) is carried to an ODE

(7)d(w21)U2wU3rU=0, (7)

where the prime denotes the derivation with respect to ξ and r is the integration constant. Balancing the highest-order derivative term U″ with the nonlinear term U3 of (7) yields N=1.

Thus, we have

(8)U(ξ)=a0+a1F+a2F2. (8)

Here, a0, a1, and a2 are constants to be determined, while F(ξ) is an unknown function to be determined. From (8), it is easy to see that

(9)U(ξ)=12(2aF+4bF3+6cF5)(a1+2a2F)+2a2(aF2+bF4+cF6) (9)

and

(10)U3(ξ)=(a0+a1F+a2F2)3. (10)

Then we substitute (8)–(10) into (7) and collect all the terms with the same power of Fi(i=0, 1, …, 6). Equating each coefficient to zero yields a set of the following algebraic equations:

(11)F0:2wa03ra0=0,F1:daa1+dw2aa16wa02a1ra1=0,F2:4daa2+4dw2aa26wa0a12ra26wa02a2=0,F3:2dba112wa0a1a2+2dw2ba12wa13=0,F4:6dba26wa0a22+6dw2ba26wa12a2=0,F5:3dca1+3dw2ca16wa1a22=0,F6:8dw2ca28dca22wa23=0. (11)

Solving the set of algebraic equations by Maple, we have the following results:

(12)a0=±r2w,a1=0,a2=±2rw(w21)bda=r2d(w21),c=(w21)b2d2r (12)

Substituting the values of a and c from (12) into (1), we obtain an ODE. Solution of this ODE for F(ξ) yields

(13)F(ξ)=rdw2bdbc1(coshaξ+sinhaξ)2. (13)

Finally we substitute (13) into (8) and we get the exact solution of (5):

(14)U(ξ)=r2w2rw(w21)bdrc1(coshaξ+sinhaξ)2dw2b+db (14)

where ξ=x+dywt,a=r2d(w21) and c1 is an arbitrary constant.

Note comparing our solutions with [30–32], it can be seen that by choosing suitable values for the parameters, similiar solutions can be verified.

3.2 Coupled Higgs Equation

We use the auxiliary method to obtain new and general exact solutions of the coupled Higgs equation

(15)uttuxx+|u|2u2uv=0,vtt+vxxx(|u|2)xx=0. (15)

Tajiri [33] obtained N-soliton solutions to the system (15). Zhao [34] and Bekir [35] constructed more general travelling wave solutions of the equation system (15).

Using the travelling wave transformation

(16)u(x,t)=eiθU(ξ),v(x,t)=V(ξ),θ=px+rt,ξ=x+wt, (16)

equation system (15) is reduced to the following ODE system:

(17)(w21)U+(p2r2)U2UV+U3=0,(w2+1)V2(U)22UU=0. (17)

Here the prime denotes derivation with respect to ξ. Integrating the second equation in the system (17) and neglecting the constant of integration we find

(18)(w2+1)V=U2, (18)

then substituting (18) into the first equation of the system and integrating we find

(19)(w41)U+(w2+1)(p2r2)U+(w21)U3=0 (19)

where the prime denotes the derivation with respect to ξ. Balancing the highest-order derivative term U″ and the nonlinear term U3 from (19) yields 3N=N+2, which gives N=1. Hence, U(ξ) takes the form (8). Substituting (8)–(10) into (19), and equating the coefficients of Fi(i=0, 1, …, 6) to zero, we separately obtain

(20)F0:p2a0r2a0w2r2a0+w2p2a0a03+w2a03=0,F1:w4aa1+3w2a02a1w2r2a1+w2p2a13a02a1+p2a1aa1r2a1=0,F2:3w2a02a23a02a2r2a2+3w2a0a12w2r2a23a0a12+p2a2+w2p2a2+4w4aa24aa2=0,F3:a13+6w2a0a1a2+w2a136a0a1a22ba1+2w4ba1=0,F4:3w2a12a2+6w4ba2+3w2a0a223a0a223a12a26ba2=0,F5:3ca1+3w2a1a22+3w4ca13a1a22=0,F6:8w4ca2+w2a23a238ca2=0. (20)

Solving the algebraic equation system by Maple, we get

(21)a0=(w2+1)(r2p2)w21,a1=0,a2=±2bw41r2p2 (21)

and

(22)a=p2r22(w21),c=(w21)b22(p2r2). (22)

Substituting (22) into (4) and solving the ODE, we have

(23)F(ξ)=p2r2c1(cosh2aξsinh2aξ)+b(1d2). (23)

Then we substitute (23) into (8) and we get the exact travelling wave solutions of equation system (15):

(24)U(ξ)=(w2+1)(r2p2)w21±2b(1w4)(p2r2)c1(cosh2aξsinh2aξ)+b(1d2) (24)

and from (18)

(25)V(ξ)=((w2+1)(r2p2)w212b(1w4)(p2r2)c1(cosh2aξsinh2aξ)+b(1d2))2w2+1, (25)

where a=p2r22(w21),ξ=x+wt and c1 is an arbitrary constant.

Note that our solutions are more extensive than the given ones [35–37]. To our knowledge, other solutions of (15) that we acquired in this article are new and are not trackable in the previous literature.

3.3 Equal Width Wave Equation

We finally consider the equal width wave (EW) equation in the form [38]

(26)ut+uuxuxxt=0. (26)

applying the transformation

(27)u(x,y,t)=U(ξ),ξ=xdt, (27)

to (27) yields

(28)dU+12U2+dU=0, (28)

where the prime denotes the derivative with respect to ξ. Balancing the highest-order derivative term U″ and the nonlinear term U2 of (29) yields N=2. Consequently, (3) takes the form

(29)U(ξ)=a0+a1F+a2F2+a3F3+a4F4. (29)

Here, ai (i=0, 1, …, 4) are constants to be determined, while F(ξ) is an unknown function to be determined. From (30) along with the auxiliary equation (4), it is easy to see that

(30)U(ξ)=(aF+2bF3+3cF5)(a1+2a2F+3a3F2+4a4F3)+(aF2+bF4+cF6)(2a2+6a3F+12a4F2), (30)

and

(31)U2(ξ)=(a0+a1F+a2F2+a3F3+a4F4)2. (31)

We substitute (29)–(31) and collect all the terms with the same power of Fi (i=0, 1, …, 8). Equating each coefficient to zero yields a set of the following algebraic equations:

(32)F8:24cda4+12a42=0F7:15cda3+a3a4=0F6:a2a4+20bda4+8cda2+12a32=0F5:12bda3+a1a4+a2a3+3cda1=0F4:16ada4da4+12a22+a0a4+6bda2+a1a3=0F3:9ada3+2bda1+a1a2+a0a3da3=0F2:da2+a0a2+12a12+4ada2=0F1:ada1+a0a1da1=0F0:12a02da0=0 (32)

Solving the equation system (32), we have the following results:

Case 1: When

(33)a0=a1=a3=0,a2=12bd,a4=48db2 (33)

and

(34)a=14,c=b2. (34)

Substituting (34) into (4), we obtain an ODE. Solution of this equation yields:

(35)F(ξ)=±1c1(coshξsinhξ)2b. (35)

Then, we substitute (35) into (30) and we find the exact solutions of (27) as follows:

U(ξ)=24bdc1(coshξsinhξ)+2b48db2(c1(coshξsinhξ)+2b)2,

where c1 is an arbitrary constant and ξ=x–dt.

Case 2: When

(36)a0=2d,a1=a3=0,a2=24bd,a4=48db2 (36)

and

(37)a=14,c=b2. (37)

Substituting (37) into (4), we obtain an ODE. Solution of this equation yields

(38)F(ξ)=1c1(cosξisinξ)+2b, (38)

and

(39)F(ξ)=1c1(cosξ2isinξ2)c1(cosξisinξ)1. (39)

Then, we substitute (38) into (30) and we find the exact solutions of (27) as follows:

U(ξ)=2d+24dbc1(cosξisinξ)2b+48db2(c1(cosξisinξ)+2b)2.

Also by substituting (39) into (30), we get the exact solutions of (27) as follows:

U(ξ)=2d24dbc1(cosξ2isinξ2)2(c1(cosξisinξ)1)+48db2c1(cosξ2isinξ2)4(c1(cosξisinξ)1)2

where c1 is an arbitrary constant and ξ=x–dt.

Case 3: When

(40)a0=2d,a1=a3=a4=0,a2=12bd (40)

and

(41)a=14,c=0. (41)

Substituting (41) into (4), we obtain an ODE. Solution of this equation yields

(42)F(ξ)=±12b(tan(ξ2+c1)2+1)btan(ξ2+c1) (42)

Then, we substitute (42) into (30) and we find the exact solutions of (27) as follows:

U(ξ)=2d3d(tan(ξ2+c1)2+1)tan((ξ2)+c1)2,

where c1 is an arbitrary constant and ξ=x–dt.

Note that our solutions are new and more extensive than the given ones in Raslan [39]. When the parameters are given special values, the solitary waves are derived from the travelling waves.

The transformed rational function method [17] paves a way of handling the solution process of nonlinear equations, unifying the tanh-function-type methods, exp-function method, mapping method, homogeneous balance method, and F-expansion-type methods. Also, its key point is to find rational solutions not only to variable-coefficient ODEs but also constant coefficient ones transformed from given nonlinear PDEs. The method given by Ma and Fuchssteiner [22] is a kind of combination of the direct method and the ansatz method. In this study, we obtain various travelling wave solutions including bright and dark solitons, periodic solutions, and exponential solutions.This method helps us find solutions to nonlinear equations through quadratures. In this article, for special cases of a, b, and c, we obtain exponential and trigonometric solutions of referred equations. For example, taking c=0 in (4), we get exponential solutions. The newly proposed technique has many advantages; it is straightforward and concise.

4 Conclusions

In this article, we used the auxiliary equation method to derive new exact solutions of travelling wave solutions of the nonlinear Zoomeron equation, coupled Higgs equation, and EW equation. All the solutions for this article have been checked by Maple software. We showed that the solutions we found are different from the solutions presented by other authors. We foresee that our results can be found potentially useful for applications in mathematical physics and engineering. Note that this method is more general and simpler; we end up very quickly with solutions using calculations with Maple. Also, exact solutions of many nonlinear differential equation and equation systems can be found with the aid of the mentioned method.


Corresponding author: Melike Kaplan, Art-Science Faculty, Department of Mathematics-Computer, Eskisehir Osmangazi University, Eskisehir, Turkey, Phone: +90 222 2393750, Fax: +90 222 2393578, E-mail:

References

[1] M. J. Ablowitz and P. A. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse Scattering Transform, Cambridge Univ. Press, Cambridge, UK 1990.10.1017/CBO9780511623998Search in Google Scholar

[2] R. Hirota, in: Backlünd Transformations, (Eds. R. Bullough, P. Caudrey), Springer, Berlin 1980.Search in Google Scholar

[3] G. T. Liu and T. Y. Fan. Phys. Lett. A 345, 161 (2005).10.1016/j.physleta.2005.07.034Search in Google Scholar

[4] M. Inc, E. Fendoglu, H. Triki, and A. Biswas. Nonlinear Anal. Model. 19, 209 (2014).10.15388/NA.2014.2.5Search in Google Scholar

[5] V. B. Matveev and M. A. Salle, Darboux Transformation and Solitons, Springer, Berlin 1991.10.1007/978-3-662-00922-2Search in Google Scholar

[6] N. C. Freeman and J. J. C. Nimmo. Proc. R. Soc. London A 389, 319 (1983).10.1098/rspa.1983.0112Search in Google Scholar

[7] A. M. Wazwaz. Appl. Math. Comput. 154, 713 (2004).10.1016/S0096-3003(03)00745-8Search in Google Scholar

[8] S. Shukri and K. Al-Khaled. Appl. Math. Comput. 217, 1997 (2010).10.1016/j.amc.2010.06.058Search in Google Scholar

[9] A. Bekir. Physica Scripta 77, 045008 (2008).10.1088/0031-8949/77/04/045008Search in Google Scholar

[10] E. Fan and H. Zhang. Phys. Lett. A 246, 403 (1998).10.1016/S0375-9601(98)00547-7Search in Google Scholar

[11] W. W. Li, Y. Tian, and Z. Zhang. Appl. Math. Comput. 219, 1135 (2012).Search in Google Scholar

[12] J. H. He and X. H. Wu. Chaos Soliton. Fract. 29, 108 (2006).10.1016/j.chaos.2005.10.100Search in Google Scholar

[13] J. H. He and M. A. Abdou. Chaos Soliton. Fract. 34, 1421 (2006).10.1016/j.chaos.2006.05.072Search in Google Scholar

[14] F. K. Chand and A, Malik, Int. J. Nonlinear Sci. 14, 416 (2012).Search in Google Scholar

[15] G. Ebadi, N. Y. Fard, H. Triki, and A. Biswas. Nonlinear Anal. Model. 17, 280 (2012).10.15388/NA.17.3.14056Search in Google Scholar

[16] M. Mirzazadeh. Inf. Sci. Lett. 3, 1 (2014).10.12785/isl/030101Search in Google Scholar

[17] W. X. Ma and J. H. Lee. Chaos Soliton. Fract. 42, 1356 (2009).10.1016/j.chaos.2009.03.043Search in Google Scholar

[18] Y. Z. Peng. Chin. J. Phys. 41, 103 (2003).10.1097/00000446-200301001-00008Search in Google Scholar

[19] N. Taghizadeh, M. Mirzazadeh, and F. Farahrooz. J. Math. Anal. Appl. 374, 549 (2011).10.1016/j.jmaa.2010.08.050Search in Google Scholar

[20] H. Jafari, A. Sooraki, Y. Talebi, and A. Biswas. Nonlinear Anal. Model. 17, 182 (2012).10.15388/NA.17.2.14067Search in Google Scholar

[21] Y. Gurefe, E. Misirli, A. Sonmezoglu, and M. Ekici. Appl. Math. Comput. 219, 5253 (2013).10.1016/j.amc.2012.11.046Search in Google Scholar

[22] W. X. Ma and B. Fuchssteiner. Int. J. Non-Linear Mech. 31, 329 (1996).10.1016/0020-7462(95)00064-XSearch in Google Scholar

[23] W. X. Ma, H. Wu, and J. He. Phys. Lett. A 364, 29 (2007).10.1016/j.physleta.2006.11.048Search in Google Scholar

[24] W. X. Ma. Sci. Chin. Math. 55, 1796 (2012).Search in Google Scholar

[25] W. X. Ma. J. Phys. Conf. Series 411, 012021 (2013).10.1088/1742-6596/411/1/012021Search in Google Scholar

[26] N. Sirendaoreji. Phys. Lett. A 356, 124 (2006).10.1016/j.physleta.2006.03.034Search in Google Scholar

[27] N. Sirendaoreji. Phys. Lett. A 363, 440 (2007).10.1016/j.physleta.2006.11.049Search in Google Scholar

[28] S. Abdoulkary, A. Mohamadou, T. Beda, B. Gambo, S. Y. Doka, Alim, and A. Mahamoudou. Appl. Math. Comput. 233, 109 (2014).10.1016/j.amc.2014.01.109Search in Google Scholar

[29] F. Calogero and A. Degasperis. Nuovo Ciment B 32, 201 (1976).10.1007/BF02727634Search in Google Scholar

[30] A. Bekir, F. Tascan, and O. Unsal. J. Assoc. Arab U Basic Appl. Sci. 17, 1 (2015).10.1016/j.jaubas.2013.12.002Search in Google Scholar

[31] K. Khan and M. A. Akbar. Shams Eng. J. 5, 247 (2014).10.1016/j.asej.2013.07.007Search in Google Scholar

[32] R. Abazari. Appl. Math. Sci. 5, 2943 (2011).Search in Google Scholar

[33] M. Tajiri. J. Phys. Soc. Jpn. 9, 2277 (1983).10.1143/JPSJ.52.2277Search in Google Scholar

[34] H. Zhao. Chaos Soliton. Fract. 36, 359 (2008).10.1016/j.chaos.2006.06.060Search in Google Scholar

[35] A. Bekir. Commun. Nonlinear. Sci. Numer. Simulat. 14, 1069 (2009).10.1016/j.cnsns.2008.05.007Search in Google Scholar

[36] A. Jabbari, H. Kheiri, and A. Bekir. Comput. Math. Appl. 62, 2177 (2011).Search in Google Scholar

[37] A. Bekir, M. Kaplan, and O. Guner. AIP Conf. Proc. 1611, 30 (2014).Search in Google Scholar

[38] A. Bekir and E. Yusufoglu. Comput. Math. Appl. 54, 1147 (2007).10.1016/j.camwa.2006.12.080Search in Google Scholar

[39] K. R. Raslan. Int. J. Nonlinear Sci. 6, 246 (2008).Search in Google Scholar

Received: 2015-3-11
Accepted: 2015-8-28
Published Online: 2015-9-24
Published in Print: 2015-10-1

©2015 by De Gruyter

Downloaded on 14.5.2024 from https://www.degruyter.com/document/doi/10.1515/zna-2015-0122/html
Scroll to top button