Scholarly article on topic 'Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form'

Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form Academic research paper on "Mathematics"

CC BY
0
0
Share paper
Academic journal
Mathematical Problems in Engineering
OECD Field of science
Keywords
{""}

Academic research paper on topic "Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form"

Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2013, Article ID 468206, 8 pages http://dx.doi.org/10.1155/2013/468206

Research Article

Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form

Reza Abazari1 and Adem Kilicman2

1 Young Researchers and Elite Club, Ardabil Branch, Islamic Azad University, P.O. Box 5616954184, Ardabil, Iran

2 Department of Mathematics, Institute of Mathematical Research, University Putra Malaysia (UPM), 43400 Serdang, Malaysia

Correspondence should be addressed to Reza Abazari; abazari-r@uma.ac.ir Received 17 July 2013; Revised 5 September 2013; Accepted 6 September 2013 Academic Editor: Guo-Cheng Wu

Copyright © 2013 R. Abazari and A. Kilicman. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents the general case study of previous works on generalized Boussinesq equations, (Abazari, 2011) and (Kilicman and Abazari, 2012), that focuses on the application of (G'/G)-expansion method with the aid of Maple to construct more general exact solutions for the coupled Boussinesq equations. In this work, the mentioned method is applied to construct more general exact solutions of Boussinesq equation and improved Boussinesq equation, which the French scientist Joseph Valentin Boussinesq (1842-1929) described in the 1870s model equations for the propagation of long waves on the surface of water with small amplitude. Our work is motivated by the fact that the (G'/G)-expansion method provides not only more general forms of solutions but also periodic, solitary waves and rational solutions. The method appears to be easier and faster by means of a symbolic computation.

1. Introduction

In the recent five decades, a new direction related to the investigation of nonlinear evolution equations (NLEEs) and processes has been actively developing in various areas of sciences. Nonlinear evolution equations have been the important subject of study in various branches of mathematical-physical sciences such as physics, fluid mechanics, and chemistry. The analytical solutions of NLEEs are of fundamental importance, since many of mathematical-physical models are described by NLEEs. Among the possible solutions to NLEEs, certain special form solutions may depend only on a single combination of variables such as solitons. In mathematics and physics, a soliton is a self reinforcing solitary wave, a wave packet or pulse, that maintains its shape while it travels at constant speed. Solitons are caused by a cancelation of nonlinear and dispersive effects in the medium. The term "dispersive effects" refers to a property of certain systems where the speed of the waves varies according to frequency. Solitons arise as the solutions of a widespread class of weakly nonlinear dispersive partial differential equations describing physical systems. The soliton phenomenon was first described by John Scott Russell (1808-1882) who observed a solitarywave in the

Union Canal in Scotland. He reproduced the phenomenon in a wave tank and named it the "wave of translation" (also known as solitary wave or soliton) [1]. The soliton solutions are typically obtained by means of the inverse scattering transform [2] and be in dept their stability to the integrability of the field equations.

In fluid mechanics, the Boussinesq approximation for water waves is an approximation valid for weakly nonlinear and fairly long waves. The approximation is named after Joseph Valentin Boussinesq (1842-1929), who first derived them in response to the observation by John Scott Russell of the wave of translation [3, 4]. According to the 1872 paper of Boussinesq, for water waves on an incompressible fluid and irrotational flow in the (x, z) plane, the boundary conditions at the free surface elevation z = q(x, t) are

v-1 - w = 0, ot ox

+ 1 (v2 + w2) + an = 0, dt 2 V ' y '

where v is the horizontal flow velocity component, v = dcp/dx, w is the vertical flow velocity component, w = df/dz,

and g is the acceleration by gravity. Now, the Boussinesq approximation for the velocity potential <p, as given previously, is applied in these boundary conditions. Further, in the resulting equations, only the linear and quadratic terms with respect to ^ and vb are retained (with vb = the

horizontal velocity at the bed z = -fo). The cubic and higher order terms are assumed to be negligible. Then, the following partial differential equations are obtained:

^ 9 ui \ i vb

^ + y ^ + ^ = V d3yb

This set of equations has been derived for a flat horizontal bed; that is, the mean depth h is a constant independent of position x. When the right-hand sides of the previous equations are set to zero, they reduce to the shallow water equations. Under some additional approximations, but at the same order of accuracy, (2) can be reduced to a single partial differential equation for the free surface elevation ^(x, i):

, , 92

= 0' (3)

In dimensionless quantities, by using the water depth h and gravitational acceleration g for nondimensionalization, (3) leads to the following, after normalization:

a2 i 2 9y\

a? b+ä?. ' = °

where f = 3(^/h), r = ^3(g/h)i, and ^ = V3(x/h). In the recent years, (4) rewrites as follows [5]:

"rt -M:o: = 0,

^2 /XX

where |q| = 1 is a real parameter. Setting q = -1 gives the good Boussinesq equation (GB) or well-posed Boussinesq equation, while by setting q = 1, we get the bad Boussinesq equation (BB) or ill-posed classical Boussinesq equation. Following Bogolubsky's modification [6] in (5) when the term qw^, is replaced with qwrt it gives the so-called improved Boussinesq equation (IBq):

"ft- uvv

Similarily, using an analogous characterization used for Boussinesq equation (5), the IBq equation for q = -1 will give the good or well-posed (GIBq), while for q = 1 the bad or ill-posed (BIBq) equation. The IBq equation appears in studying the transverse motion and nonlinearity in acoustic waves on elastic rods with circular cross-section. In particular, the BIBq is used to discuss the wave propagation at right angles to the magnetic field and also to approach the bad BS equation (see Makhankov [7]) or to study ion-sound(s) waves (see Bogolubsky [6]).

There are some review articles and some collected works that have been focused to study the classical Boussinesq equation from various points of view. The initial boundary value and the Cauchy problem of (5) have been described in [8-11]. Yajima [12] has studied the nonlinear evolution of a linearly stable solution, while the exponentially decaying solution of the spherical Boussinesq equation was obtained by Nakamura [13]. The global existence of the strong solution and the small amplitude solution for the Cauchy problem of the multidimensional equation (5) is proved in [14]. A general approach to construct exact solution to (5) is given by Clarkson [9], and Hirota [10] has deduced conservation laws and has examined N-soliton interaction. Bona and Sachs, in [8], have discussed that the special solitary-wave solutions for (5), when nonlinear term is u , are nonlinearly stable for a range of their wave speeds.

On the other hand, recently, the (G'/G)-expansion method, firstly introduced by Wang et al. [15], has become widely used to search for various exact solutions of NLEEs [15-19]. The value of the (G'/G)-expansion method is that one treats nonlinear problems by essentially linear methods. The method is based on the explicit linearization of NLEEs for traveling waves with a certain substitution which leads to a second-order differential equation with constant coefficients. Moreover, it transforms a nonlinear equation to a simple algebraic computation. Although many efforts have been devoted to find various methods to solve (integrable or nonintegrable) NLEEs, there is no unified method. The main merits of the (G'/G)-expansion method over the other methods are that it gives more general solutions with some free parameters which, by suitable choice of the parameters, turn out to be some known solutions gained by the existing methods.

Our first interest in the present work is in implementing the (G'/G)-expansion method to show its power in handling nonlinear partial differential equations (PDEs), so that one can apply it to other models of various types of nonlinearity. The next interest is in the determination of exact travelling wave solutions for generalized equations (5) and (6).

2. Description of the (G'/G)-Expansion Method

The objectiveofthissection is to outlinethe use expansion method for solving certain nonlinear PDEs. Suppose that we have a nonlinear PDE for u(x, i), in the form

') = °,

where P is a polynomial in its arguments, which includes nonlinear terms and the highest order derivatives. The transformation w(x, i) = £ = fcx + «i, reduces (7) to the ordinary differential equation (ODE)

P (u, H/', wü', fc V, fcwü", «V, ''') = 0, (8)

where U = and prime denotes derivative with respect to We assume that the solution of (8) can be expressed by a polynomial in (G'/G) as follows:

m ( r'

ordinary differential equation (LODE):

+ «0> «m = 0

+ pG(Ç) = 0,

where a0 and a;, for i = 1,2,..., m, are constants to be determined later and G(£) satisfies a second order linear _I

where A and ^ are arbitrary constants. Using the general solutions of (10), we have

G' (fl)

^Я2 -4p ( Ci sinh ((VA2 -4F/2)fl) + C2 cosh ((^A2 - 4^/2) fl) 2 yCj cosh ((^A2 -4^/2)fl) +C2 sinh ((^A2 -4^/2)fl)

V^-A2 ^ -Cx sin ((^4^-а2/2)а) + С2 cos ((^ - A2/2) fl) Cx cos ((^4^-A2/2)fl) + C2 sin ((^4^-A2/2)fl)

A2 - 4^ > О,

- -, A2 - 4^ < О,

and it follows from (9) and (10), that

\i+i / ^ \ i / ^ \ i-i G^ G \ ( G

и" = £i«i ((1 + 1)(-) + (21 + 1) a( —

I(A2 + 2p)

+ (21 - 1) A^

+ («-1)p2( f

and so on, here the prime denotes the derivative with respective to To determine m explicitly, we take the following four steps.

Step1. Determine the integer m by substituting (9) alongwith (10) into (8) and balance the highest order nonlinear term(s) and the highest order partial derivative.

Step 2. Substitute (9) to give the value of m determined in Step 1 along with (10) into (8) and collect all terms with the same order of (G'/G) together, the left-hand side of (8) is converted into a polynomial in (G'/G). Then, set each coefficient of this polynomial to zero to derive a set of algebraic equations for fc, to, a0, and a;, for i = 1,2,..., rn.

Step 3. Solve the system of algebraic equations obtained in Step 2, for fc, to, a0, and , for i = 1,2,..., rn, by the use of Maple.

Step 4. Use the results obtained in the above steps to derive a series of fundamental solutions м(£) of (8) depending on (G'/G); since the solutions of (10) have been well known for us, then we can obtain exact solutions of (7).

3. Application

In this section, we will demonstrate the (G'/G)-expansion method on three of the well-known Boussinesq type equations (5) and (6).

3.1. Boussinesq Equation. To look for the traveling wave solution of Boussinesq equation (5), we use the gauge transformation:

и (я, i) = U (fl) '

where £ = fcx + wt, and fc and to are constants. We substitute (13) into (5) to obtain the nonlinear ordinary differential equation

(w2 -fc2)u" -fc2(±u2 + )" = О.

According to Step 1, we get m + 4 = 2m + 2; hence, m = 2. We then suppose that (14) has the following formal solutions:

U = a-,

+ a0, a2 = 0,

where a2, ai, and a0 are constants which are unknowns to be determined later. Substituting (15) along with (10) into (14) and collecting all terms with the same order of (G'/G) together, the left-hand sides of (14) are converted into a polynomial in (G'/G). Setting each coefficient of each

polynomial to zero, we derive a set of algebraic equations for k, w, A, m, and a2 as follows:

: (-^A (A2 + 8^) a - 2^2a2 (8^ + 7A2)) ^k4

+ (-2«2^2 - Aa (1 + 2a0) ^ - 2^2 (1 + 2a0) a2) k2 + 2w2a2^2 + = 0,

— ) : ((-22^A2 - 16^2 - A4) a

-30^Aa2 (4^ + A2)) ^k4 + (-6a2A^ + (-2^ - A2 - 2a0A2 - 4a0^ - 12a2^2) a -6^Aa2 (1 + 2a0) ) k2 + w2 (2^ + A2) a

+ 6w2a2A^ = 0,

:(-15A(4^ + A2) a

-8a2 (29^A2 + 2A4 + 17^2)) ijk4 + ((-4A2 - 8^) a2 - 3A (2a0 + 1 + 10^a2) a -4a2 (3a2^2 + 2^ + A2 + 4a0^ + 2a0A2)) k2 + 3w2a1A + 4w2a2 (2^ + A2) = 0,

: ((-40p - 50A2) a - 10Aa2 (13A2 + 44^)) ijk4

+ (-10a2A + (-36^a2 - 18a2A2 - 2 - 4a0) a

-2a2A (14^a2 + 5 + 10a0) ) k2 + 2w2a + 10w2a2A = 0,

: (-60aA-30a2 (11A2 + 8^))ijk4

+ ( - 6«2 - 42 a2aA

-2 a2 (16^a2 + 8a2A2 + 3 + 6a0)) k2

+ 6w a2 = 0,

— ) : (-24a - 336a2A) ^k4

+ (-36a^A - 24a2a) k2 = 0,

— ) : -120^k4a2 - 10fc2«2 = 0.

20 10 0

-10 -20

Figure 1: Hyperbolic function solution (20a) of the Boussinesq equation (5) for q = -1, k = -1, W = 1/2, A = 4, m = 1, and = 0.

Solving the obtained algebraic equations by the use of Maple, we get the following results:

.4.. , „r-412 ,.,2 , ,2

+ gk4A2 - w2 + k kr

a = -12k2^A, a2 = -12k2^

and k, w, A, and m are arbitrary constants. Therefore, substitute the previous case in (15), we get

U= -12k2<î(— ) -12kV(——-

8^kV + gk4A2 - w2 + k2

Substituting the general solutions (11) into (18), we obtain three types of traveling wave solutions of Boussinesq equation (5) in the view of the positive, negative, or zero of A2 - 4^.

When D = A2 - 4^ > 0, we obtain hyperbolic function solution UH of Boussinesq equation (5) (see Figure 1) as follows:

-3fc2q (( (Vd [C sinh ((1/2) Vd.Ç)

+ C2 cosh ((1/2) VdÇ)]) x (c2 sinh ((1/2) Vd^)

+ C cosh ((1/2) VdÇ))-1) - A

■6k2qA(((Vd [Cj sinh ((1/2) Vd^)

+ C2 cosh ((1/2) Vd)]) x (C2 sinh ((1/2) Vd^)

+ Cj cosh ((1/2) Vd^))-1) - A 8qk4p + qk4A2 - c2 + k2

where £ = kx + o>i, and C2 are arbitrary constants. It is easy to see that the hyperbolic solution (19) can be rewritten at Cl > C2 as follows:

(x,f) = -3k2qdtanh2 (1vd.Ç +

-2qk4D - c2 + k2 k2 '

while at C < C2, one can obtain

mh (x,f) = -3k2qDcoth2 +

-2qk4D - c2 + k2 k2 '

where £ = kx + cot, = tanh 1(C1/C2), and k, c, A, and ^ are arbitrary constants. Now, when d = A2 - 4^ < 0, the trigonometric function solutions UT of Boussinesq equation (5) will be

Ut (?)

= -3k2q(((V-d [-C1 sin((1/2) V-dÇ)

+ C2 cos ((1/2)v-d^ )]) x (C2 sin ((1/2) >/=dÇ)

+ C1 cos

-6k2qA(((V-d [-C1 sin ((1/2) V-dÇ)

+C2 cos ((1/2) v/-D£)]) x (C2 sin ((1/2) >/=dÇ)

+C1 cos ((1/2) v-d^))-1) - A)

8qk> + qk4A2 - c2 + k к2

20 0 -20 -40

Figure 2: Trigonometric function solution (22a) of the Boussinesq equation (5) for q = -1, к = -1, c = 1/20, A = 4, ^ = -1, and

Пт =°.

where £ = fcx + «i, and C1, C2 are arbitrary constants (see Figure 2). Similarly, the trigonometric solutions (21) can be rewritten at C > C and C < C2, respectively, as follows:

(x, i) = 3k2qdtan2 ( 1 +

-2qk4D - c2 + k2 k2"

(x,f) = 3k2qdcot2 (^V-d +

-2qk4D - c2 + k2 k2"

where £ = fcx + cot, = tan 1(C1/C2), and fc, w, A, and ^ are arbitrary constants. Finally, when A2 - 4^ = 0, then the rational function solutions to (5) will be

"rat (*,f) = -"

12k2qC2

+ C-1,

(C1 + C2 (kx + ci))2 k2 where C1, C2, k, and cc are arbitrary constants.

3.2. Improved Boussinesq Equation. Similar to the previous section, to obtain the traveling wave solution of improved Boussinesq equation (6) we substitute the gauge transformation (13) into (6) to obtain nonlinear ordinary differential equation

(c2 - k2) U" -k2(2u2 +qc2U") =0.

According to Step 1, we get m + 4 = 2m + 2; hence, m = 2. Then, similar to the previous section, we suppose that (24) has the same formal solutions (15). Substituting (15) along with (10) into (24) and collecting all terms with the same order of (G'/G) together, the left-hand sides of (24) are converted into a polynomial in (G'/G). Setting each coefficient of each

polynomial to zero, we derive a set of algebraic equations for fc, to, A, m, and a2 as follows:

( ) : ((-^to2A (A2 + 8^) aj - 2^2to2a2 (8^ + 7A2)) ^

-2«2- A^ (1 + 2a0) a1 - 2^2a2 (1 + 2a0)) fc

+ 2to2 a2^2 + to2a1A^ = О,

: ((-to2 (22A2^ + 16^2 + A4)a1

^ ) : (-1О«2 - 120^to2a2) fc2 = О,

-30^to2a2A (A2 + 4^)) ^ - 6a12A^ + (-2^ - A2 - 12a2- 4a0^ - 2a0A2) a1 -6A^a2 (1 + 2a0) ) fc2 + to2 (2^ + A2) a1 + 6to2a2A^ = О,

: ((-15to2A(A2 + 4^) a1

-8to2a2 (29A2^ + 2A4 + 17^2)) q + (-4A2 - 8^) a]2 - 3A (1 + 2a0 + 1Оа2^) a1 -4a2 (3a2+ 2^ + A2 + 4a0^ + 2a0A2)) fc2 + 3to2a1A + 4to2a2 (2^ + A2) = О,

: ((-10to2 (4^ + 5A2)a1

-10to2a2A (13A2 + 44^)) q - 1О«2A + (-4a0 - 2 - 36a2^ - 18a2A2) a1 -2a2A (14a2^ + 1Оа0 + 5) ) fc2 + 10to2a2A + 2to2a1 = О,

: ((-60to2a1A - 30to2a2 (11A2 + 8^)) q

- 6«2 - 42a2 a1A

-2a2 (16a2^ + 8a2A2 + 3 + 60,3) ) fc2 + 6to2 a2 = О,

: ((-24to2a1 - 336to2a2A) q - 36a^A - 24a2a1) fc = О,

-10 -3

Figure 3: Hyperbolic function solution (20a) of the improved Boussinesq equation (6) for q = -1, fc = -1, to = 1/2, A = 4, ^ = 1, and -цж = О.

and solving by use of Maple, we get the following results:

«П = -

8^fc2to2^ + ijfcVA2 + fc2 - to2 fc2 ,

a1 = -12to2^A, a2 = -12to2^

and fc, to, A, and ^ are arbitrary constants. Therefore, the solution (15) leads to

U = - 12to2^( ^ ) - 12to2^A ( ^

8^fc2to2^ + gfcVA2 + fc2 - to2 kr

Now, for d = A2-4^ > О and d = A2-4^ < О, the hyperbolic function solution and trigonometric function solution of improved Boussinesq equation (6) are obtained as follows (see Figures 3 and 4), respectively:

-3to2^(((VD [c1 sinh ((1/2) VDfl)

+ C2 cosh ((1/2) VDfl)]) x(c2 sinh ((1/2) VDfl)

+ C1 cosh ((1/2) VDfl))-1 ) - A

-0.9 -0.95 -1 -1.05

Figure 4: Trigonometric function solution (22a) of the improved Boussinesq equation (6) for q = -1, k = -1, c = 1/20, A = 4, p = -1, and чт = 0.

■6c2qA(((Vd [C1 sinh ((1/2) VdÇ)

+ C2 cosh ((1/2) VdÇ)]) x (c2 sinh ((1/2) Vd^)

+ C1 cosh ((1/2)Vd^))-1)-A

8qk2«2^ + qk2c2A2 + k2 - c2 kr

= -3c2q (( ( V-D [-C1 sin ((1/2) V-D)

+ C2 cos ((1/2) V-D)]) x (C2 sin ((1/2) V-DÇ)

+ C1 cos ((1/2) V-DÇ))-1) - A)2 - 6c2qA (( ( V-D [-C1 sin ((1/2) V-D£) +C2 cos ((1/2) V-D)]) x (C2 sin ((1/2) V-DÇ)

+C1 cos((1/2) V-DÇ))-1 )-A) 8qk2«2^ + qk2c2A2 + k2 - c2

where £ = kx+w>i, and C, C2 are arbitrary constants. It is easy to see that the hyperbolic solution (28) and trigonometric

solution (29) can be rewritten at C > C2 as follows:

.(x,i) = -3c2qDtanh2 (^vD; + 4 -2qk2c2D - c2 + k2

(x, i) = 3c2qdtan2 ( 1V-d; + чт )

-2qk2c2D - c2 + k2 k2"

while at C2 < C2, one can obtain

(x,i) = -3c2qdcoth2 (^Vd^ +

-2qk2c2D - c2 + k2 k2 '

(x, i ) = 3c2qdcot2 ( 1V-d + Чт

-2qk 2c2D - c2 + k2 k2

where £ = fcx + «i, = tanh-1(C1/C2) and Ъ, « A, and ^ are arbitrary constants.

Finally, when A2 -4^ = 0, then the rational function solutions of improved Boussinesq equation (6) will be

"rat (X, 0 = -

12c2 qC2

(C1 + C2 (kx + c))2 k2

where C:, C2, k, and w are arbitrary constants.

4. Conclusions

This study shows that the (G'/G)-expansion method is quite efficient and practically well suited for use in finding exact solutions for the Boussinesq equation and improved Boussinesq equations. The reliability of the method and the reduction in the size of computational domain give this method a wider applicability. Though the obtained solutions represent only a small part of the large variety of possible solutions for the equations considered, they might serve as seeding solutions for a class of localized structures existing in the physical systems. Furthermore, our solutions are in more general forms, and many known solutions to these equations are only special cases of them. With the aid of Maple, we have assured the correctness of the obtained solutions by putting them back into the original equation.

Conflict of Interests

This paper does not have any conflict of interest with the authors research topics.

Acknowledgment

The first author would like to thank the Young Researchers

and Elite Club, Islamic Azad University, Ardabil Branch for

its financial support.

References

[1] J. S. Russell, "Report on waves," in Proceedings of the 14th Meeting of the British Association for the Advancement of Science, 1844.

[2] M. J. Ablowtiz and J. F. Ladik, "On the solution of a class of nonliear partial difference equations," Studies in Applied Mathematics, vol. 57, pp. 1-12,1977.

[3] J. Boussinesq, "Thorie de l'intumescence liquide, applele onde solitaire ou de translation, se propageant dans un canal rectangulaire," Comptes Rendus de l'Académie des Sciences, vol. 72, pp. 755-759, 1871.

[4] J. Boussinesq, "Thorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond," Journal de Mathématiques Pures et Appliquées, vol. 17, pp. 55-108,1872.

[5] O. V. Kaptsov, "Construction of exact solutions of the Boussinesq equation," Journal of Applied Mechanics and Technical Physics, vol. 39, no. 3, pp. 389-392,1998.

[6] I. L. Bogolubsky, "Some examples of inelastic soliton interaction," Computer Physics Communications, vol. 13, no. 3, pp. 149155,1977.

[7] V. G. Makhankov, "Dynamics of classical solitons (in non-integrable systems)," Physics Reports, vol. 35, no. 1, pp. 1-128, 1978.

[8] J. L. Bona and R. L. Sachs, "Global existence of smooth solutions and stability of solitary waves for a generalized Boussinesq equation," Communications in Mathematical Physics, vol. 118, no. 1, pp. 15-29, 1988.

[9] P. A. Clarkson, "New exact solutions of the Boussinesq equation," European Journal of Applied Mathematics, vol. 1, no. 3, pp. 279-300, 1990.

[10] R. Hirota, "Solutions of the classical boussinesq equation and the spherical boussinesq equation: the Wronskian technique," Journal of the Physical Society of Japan, vol. 55, no. 7, pp. 21372150, 1986.

[11] Y. Liu, "Instability and blow-up of solutions to a generalized Boussinesq equation," SIAM Journal on Mathematical Analysis, vol. 26, no. 6, pp. 1527-1546, 1995.

[12] N. Yajima, "On a growing mode of the Boussinesq equation," Progress of Theoretical Physics, vol. 69, no. 2, pp. 678-680,1983.

[13] A. Nakamura, "Exact solitary wave solutions of the spherical boussinesq equation," Journal of the Physical Society of Japan, vol. 54, no. 11, pp. 4111-4114, 1985.

[14] S. B. Wang and G. W. Chen, "Small amplitude solutions of the generalized IMBq equation," Journal of Mathematical Analysis and Applications, vol. 274, no. 2, pp. 846-866, 2002.

[15] M. Wang, X. Li, and J. Zhang, "The (frac(G'/G))-expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics," Physics Letters A, vol. 372, no. 4, pp. 417-423, 2008.

[16] R. Abazari, "The (G' /G) -expansion method for the coupled Boussinesq equation," Procedia Engineering,, vol. 10, pp. 28452850, 2011.

[17] E. M. E. Zayed and K. A. Gepreel, "The (G'/G) method for finding traveling wave solutions of nonlinear partial differential equations in mathematical physics," Journal of Mathematical Physics, vol. 50, no. 1, Article ID 013502, 12 pages, 2009.

[18] Y. He, S. Li, and Y. Long, "Exact solutions of the Kudryashov-Sinelshchikov equation using the multiple (G'/G)-expansion method," Mathematical Problems in Engineering, vol. 2013, Article ID 708049, 7 pages, 2013.

[19] A. Kilicman and R. Abazari, "Travelling wave solutions of the Schrodinger-Boussinesq system," Abstract and Applied Analysis, vol. 2012, Article ID 198398,11 pages, 2012.

Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.