Hindawi Publishing Corporation Abstract and Applied Analysis Volume 2012, Article ID 198398,11 pages doi:10.1155/2012/198398

Research Article

Travelling Wave Solutions of the Schrodinger-Boussinesq System

Adem Kilicman1 and Reza Abazari2

1 Department of Mathematics and Institute of Mathematical Research, Universiti Putra Malaysia, 43400 UPM, Serdang, Selangor, Malaysia

2 Young Researchers Club, Ardabil Branch, Islamic Azad University, P.O. Box. 5616954184, Ardabil, Iran

Correspondence should be addressed to Adem Kilicman, akilicman@putra.upm.edu.my Received 16 August 2012; Accepted 4 October 2012 Academic Editor: Mohammad Mursaleen

Copyright © 2012 A. Kilicman and R. Abazari. 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.

We establish exact solutions for the Schrodinger-Boussinesq System iut + uxx - auv = 0, vtt - vxx + vxxxx - b(\u\2)xx = 0, where a and b are real constants. The (G'/G)-expansion method is used to construct exact periodic and soliton solutions of this equation. Our work is motivated by the fact that the (G'/G)-expansion method provides not only more general forms of solutions but also periodic and solitary waves. As a result, hyperbolic function solutions and trigonometric function solutions with parameters are obtained. These solutions may be important and of significance for the explanation of some practical physical problems.

1. Introduction

It is well known that the nonlinear Schrodinger (NLS) equation models a wide range of physical phenomena including self-focusing of optical beams in nonlinear media, the modulation of monochromatic waves, propagation of Langmuir waves in plasmas, and so forth. The nonlinear Schrodinger equations play an important role in many areas of applied physics, such as nonrelativistic quantum mechanics, laser beam propagation, Bose-Einstein condensates, and so on (see [1]). Some properties of solutions for the nonlinear Schrodinger equations on R" have been extensively studied in the last two decades (e.g., see [2]).

The Boussinesq-type equations are essentially a class of models appearing in physics and fluid mechanics. The so-called Boussinesq equation was originally derived by Boussinesq [3] to describe two-dimensional irrotational flows of an inviscid liquid in a uniform rectangular channel. It also arises in a large range of physical phenomena including the propagation of ion-sound waves in a plasma and nonlinear lattice waves. The study on the

soliton solutions for various generalizations of the Boussinesq equation has recently attracted considerable attention from many mathematicians and physicists (see [4]). We should remark that it was the first equation proposed in the literature to describe this kind of physical phenomena. This equation was also used by Zakharov [5] as a model of nonlinear string and by Falk et al. [6] in their study of shape-memory alloys.

In the laser and plasma physics, the Schrodinger-Boussinesq system (hereafter referred to as the SB-system) has been raised. Consider the SB-system

where t > 0, x e [0,L], for some L > 0, and a, b are real constants. Here u and v are, respectively, a complex-valued and a real-valued function defined in space-time [0, L] x R. The SB-system is considered as a model of interactions between short and intermediate long waves, which is derived in describing the dynamics of Langmuir soliton formation and interaction in a plasma [7] and diatomic lattice system [8]. The short wave term u(x,t) : [0,L] x R ^ C is described by a Schrodinger type equation with a potential v(x,t) : [0,L] x R ^ R satisfying some sort of Boussinesq equation and representing the intermediate long wave. The SB-system also appears in the study of interaction of solitons in optics. The solitary wave solutions and integrability of nonlinear SB-system has been considered by several authors (see [7, 8]) and the references therein.

In the literature, there is a wide variety of approaches to nonlinear problems for constructing travelling wave solutions, such as the inverse scattering method [9], Bcklund transformation [10], Hirota bilinear method [11], Painleve expansion methods [12], and the Wronskian determinant technique [13].

With the help of the computer software, most of mentioned methods are improved and many other algebraic method, proposed, such as the tanh/coth method [14], the Exp-function method [15], and first integral method [16]. But, most of the methods may sometimes fail or can only lead to a kind of special solution and the solution procedures become very complex as the degree of nonlinearity increases.

Recently, the (G'/G)-expansion method, firstly introduced by Wang et al. [17], has become widely used to search for various exact solutions of NLEEs [17-19].

The main idea of this method is that the traveling wave solutions of nonlinear equations can be expressed by a polynomial in (G'/G), where G = G(£) satisfies the second order linear ordinary differential equation G"(£) + AG'(£) + nG(Q = 0, where I = kx + wt and k, w are arbitrary constants. The degree of this polynomial can be determined by considering the homogeneous balance between the highest order derivatives and the non-linear terms appearing in the given non-linear equations.

Our first interest in the present work is in implementing the (G'/G)-expansion method to stress its power in handling nonlinear equations, so that one can apply it to models of various types of nonlinearity. The next interest is in the determination of exact traveling wave solutions for the SB-system (1.1).

iut + uxx - auv = 0,

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

The objective of this section is to outline the use of the (G'/G)-expansion method for solving certain nonlinear partial differential equations (PDEs). Suppose that a nonlinear equation, say in two independent variables x and t, is given by

P(U, Ux, Ut, Uxx, Ux,t, utt,...) = 0, (2.1)

where u(x, t) is an unknown function, P is a polynomial in u(x, t) and its various partial derivatives, in which the highest order derivatives and nonlinear terms are involved. The main steps of the (G'/G)-expansion method are the following:

Step 1. Combining the independent variables x and t into one variable £ = kx+wt, we suppose that

u(x, t) = U(£), £ = kx + wt. (2.2)

The travelling wave variable (2.2) permits us to reduce (2.1) to an ODE for u(x,t) = U(£), namely,

p(u,kU',wU',k2U", kwU", w2U",..= 0, (2.3)

where prime denotes derivative with respect to £.

Step 2. We assume that the solution of (2.3) can be expressed by a polynomial in (G'/G) as follows:

m / G'\ i

U(£) = E^yc) + a0/ am /0, (2.4)

where m is called the balance number, a0, and ai, are constants to be determined later, G(£) satisfies a second order linear ordinary differential equation (LODE):

d2G(£) , XdG(£) 0

~d£2 + X~!T + *G(£) =0, (2.5)

where X and * are arbitrary constants. The positive integer m can be determined by considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in ODE (2.3).

Step 3. By substituting (2.4) into (2.3) and using the second order linear ODE (2.5), collecting all terms with the same order of (G'/G) together, the left-hand side of (2.3) is converted into another polynomial in (G'/G). Equating each coefficient of this polynomial to zero yields a set of algebraic equations for k,w,X,*,a0,a\,.. .,am.

Step 4. Assuming that the constants k,w,X,^,a0,a1, ...,am can be obtained by solving the algebraic equations in Step 3, since the general solutions of the second order linear ODE (2.5) is well known for us, then substituting k,w,X, ¡i,a0, ...,am and the general solutions of (2.5) into (2.4) we have more travelling wave solutions of the nonlinear evolution (2.1).

In the subsequent section, we will illustrate the validity and reliability of this method in detail with complex model of Schrodinger-Boussinesq System (1.1).

3. Application

To look for the traveling wave solution of the Schrodinger-Boussinesq System (1.1), we use the gauge transformation:

u(x,t) = U(t,)ein,

v(x,t) = V (I),

where I = kx + wt, n = px + qt, and p, q, w are constants and i = V-T. We substitute (3.1) into

(1.1) to obtain nonlinear ordinary differential equation

k2U" - i(2kp + w)U' - aUV - (p2 + qjU = 0, (3.2)

(w2 - k2)V + k4V(4) - bk2 (u2)" = 0. (3.3)

In order to simplify, integrating (3.3) twice and taking integration constant to zero, the system

(3.2)-(3.3) reduces to the following system:

k2U" - i(2kp + w)U' - aUV - (p2 + q\U = 0,

(w2 - k2j V + k4V" - bk2U2 = 0.

Suppose that the solution of the nonlinear ordinary differential system (3.4) can be expressed by a polynomial in (G'/G) as follows:

m / G'\ i

U(g) = £ ai\~Q/ + a0, am = 0,

'GAi G > + P0,

n /G'\'

V(D = ZP\G) + ß0, ßn / 0,

where m, n are called the balance number, ai, (i = 0,1,...,m) and fy, (j = 0,1,...,n) are constants to be determined later, G(£) satisfies a second order linear ordinary differential equation (2.5). The integers m, n can be determined by considering the homogeneous balance

between the highest order derivatives and nonlinear terms appearing in nonlinear ordinary differential system (3.4) as follow:

m + n = m + 2,

2m = n + 2,

so that m = n = 2. We then suppose that (3.4) has the following formal solutions:

G'\2 ZG

U(£) = -^J + «W )+ «2 ^ 0,

, 2 , (3.7)

V (£) = ß^ G)Z + §) + ß0, ß2 = 0.

Substituting (3.7) along with (2.5) into (3.4) and collecting all the terms with the same power of (G'/G) together, equating each coefficient to zero, yields a set of simultaneous algebraic equations for k, w, X, /, aj, and ßj, (j = 0,1,2), as follows:

2k a2^ + iwa1 n + k a1X^ - aa0fa0 - p a0 - qa0 + 2ikpa1n = 0,

w2po + 2k4p2H2 - k2 (-^Xpik2 + bao2 + fa) = 0, k2a1X2 - p2a1 + (4ika2^ + 2ika1 X)p + (2ia2n + ia1X)w - qa1 - aa0fa + 2k2a1n - aa1fa + 6k2a2X^ = 0, + 6k4p2X^ - k2(-k2X2fai - 2k2^fai + fai + 2baia^ = 0, (4ika2X + 2ika1)p + (2ia2X + ia{)w - p2a2 + ^4X2k2 - q - afa0 + 8^k2^a2 + 3k2a1X - aa1fa1 - aa0fa2 = 0, ^2fa2 + k2(4X2k2 + 8^k2 - l) fa2 - k2(lba2a0 + bai2 - 3Xfaik2) = 0, 4ikpa2 + 2iwa2 + ^10Xk2 - afa^a2 - a1 ^-2k2 + afa2) = 0, 10k4fa2X - 2k2 (-k2fai + ba2ai) = 0, (6k2 - afaija2 = 0, 6k4fa2 - bk2a22 = 0.

Then, explicit and exact wave solutions can be constructed through our ansatz (3.7) via the associated solutions of (2.5).

In the process of constructing exact solutions to the Schrodinger-Boussinesq system (1.1), (2.5) is often viewed as a key auxiliary equation, and the types of its solutions determine

the solutions for the original system (1.1) indirectly. In order to seek more new solutions to (1.1), we here combine the solutions to (2.5) which were listed in [18,19]. Our computation results show that this combination is an efficient way to obtain more diverse families of explicit exact solutions.

Solving (3.8) by use of Maple, we get the following reliable results:

, W3«2 - 3k2 1 k2 - w2 1 w 1 = ±3--' * = 6 P = - 2 -

1 24k4 + 30-4 - 54k2-2 q = 24-k2(k2 - -2)-^ = 0

ai = ±21/-1- V3-2 - 3k2, a2 = 6k2V-T, V ab V ab

a _ „ V3-2 - 3k2 6k2

ßo = 0, ßi = ±2-, ß2 = —

where k and w are free constant parameters and k / ± w. Therefore, substitute the above case in (3.7), we get

u®=6klivXG) ±2iO-b'^^i§

V ({) = * ( §>2

(3.10)

V3-2 - 3k v G

Substituting the general solutions of ordinary differential equation (2.5) into (3.10), we obtain two types of traveling wave solutions of (1.1) in view of the positive and negative of I2 - 4^.

When D = l2-4¡i = (w2-k2)/k4 > 0, using the general solutions of ordinary differential equation (2.5) and relationships (3.10), we obtain hyperbolic function solutions uH(x, t) and vH(x, t) of the Schrodinger-Boussinesq system (1.1) as follows:

uH(x,t) = 2^

f C1 sinh((VD/2^ ^ + C2 cosh

2 I Ci cosh ( (■VD/2) ^ + C2 sinh

1 V3-2 - 3k2

T 6 k2

± V3—2 - 3k2

C1 sinh^Vä/2) t) + Cicosh((VD/2) t)

) t) + C2 si

1 V3—2 - 3k2 T 6 k2

^ (x t) = 6k! f ^ ^^(Q^2) 0 + Ciœsh((vW2) M ^ 1 ^—2 - 3k2 a I 2 I C1 cosh^ (vD/2) t) + C2 sinh ^ (vD/2) t)/ 6

V3—2 - 3k2

^VD^ C1sinh((vWl) t) + Cicosh((vD/l) t

^ 2 ^Qcosh((VD/2) t) + C2sinh((VD/2) t) J

1 V3—2 - 3k2 T 6 k2

(3.11)

where D = (w2 - k2)/k4 > 0, £ = kx + wt, rç = -(1/2)(wx/k) + (1/24)((24k4 + 30w4 -54k2w2)/k2(k2 - w2))t, and k,w,C1, and C2 are arbitrary constants and k = ± w.

It is easy to see that the hyperbolic solutions (3.11) can be rewritten at C2 > C2, as follows:

/ s ^/1 —2 - k2 /„ o/l —2 - k2 . „ . in

"H(x'f) = 2 V ^b"fc!" ( 3tanh2( 2V"F" t + Ph )- 1

VH(x,t) =

1 —2 - k2

2 ak2 V3tanh2( 11/ "t + Ph) - 1

(3.12a)

(3.12b)

while at C2 < C2, one can obtain

. x 1J1 —2 - k2 / ,1(1. —2 - k2 . \ „ \ in UH(x,t) = 5Vk2-- { 3coth2{ ^l—z— t + PH ) - 1 )ein, (3.12c)

VH(x,t) = 1—O--2 {3coth2 {l^/-k-2 t + Ph) - 1), (3.12d)

where £ = kx + wt, rç = -(1/2)(wx/k) + (1/24)((24k4 + 30w4 -54k2w2)/k2(k2 - w2))t, p^ = tanh-1(C1/C2), and k, w are arbitrary constants.

Now, when D = I2 - = (w2 - k2)/k4 < 0, we obtain trigonometric function solutions uT and vT of Schrodinger-Boussinesq system (1.1) as follows:

uT (x,t) = 21

( V-Df~Clsin^V-D/2) + C2cos((V-D/2) ¿)\

I 2 V Cicos((V-D/2) + C2sin( ( V-D/2

i V3w2 - 3k2

± V3w2 - 3k2

-Cisin((V-D/2) + C2cos((VD/2)

2 \ C1 cos ((V-D/2) + C2sin((V-D/2)

i V3w2 - 3k2 T 6 k2

, , 6kYy-D/^ism^v^2) + Qcos^v^W2) M i yw

(x,t)= ---// ,__^-// ,__T" 1^7-T

~D/2) 0

^ 2 \ Cicos((V-D/2) + C2

V3^2 - 3k2f V-Df-^s^V^n) + ^((v^2) ¿) fl I 2 I Ci cos((V=D/2) + C2sin( (V-D/2)

(3.i3)

where D = (w2 - k2)/k4 < 0,1 = kx + wf, rç = -(1/2)(wx/k) + (1/24)((24k4 + 30w4 -54k2w2)/k2(k2 - w2))f, and k, w, C1, and C2 are arbitrary constants and k= ± w. Similarity, the trigonometric solutions (3.13) can be rewritten at C2 > C^, and C2 < C2, as follows:

. x i.fT k2 - w2 / 2/i. Ik2 - w2 > \ „ x

UT(x,t) = 2V Vb,— ^3tan2^ 2V — * + ) + i)en, (3.i4a)

. x i k2 - w2 / 2/i. Ik2 - w2 „ \ \ /„.,,, x

^T(x,t) = 2( n2 ( 2V* + PT) + i ) (3.i4b)

Figure 1: Soliton solution |uH (x, t) | of the Schrodinger-Boussinesq system, (3.12a), for a = 1, p = 16, k = 1, w = -2, and pH = tanh-1(1/4).

Figure 2: Soliton solution |uH (x, t) |2 of the Schrodinger-Boussinesq system, (3.12a), for a = 1, p = 16, k = 1, w = -2, and pH = tanh-1(1/4).

. x 1. IT" k2 - w2f Jl. Ik2 - w2 „ \ „ \ ,„

«t(*,0 = ^-b— ^3cot2^2V^ ^ + PTJ + 1 j^, (3.14c)

, N 1 k2 - w2/ Ik2 - w2 , \ \ .„1ilJN

^T(^f) = 2 -mr ( 3cot2 ( 2 V £ + PTJ + 1J, (3.14d)

where £ = kx + wt, n = -(1/2)(wx/k) + (1/24)((24k4 + 30w4 - 54k2w2)/k2(k2 - w2))t, = tan-1(Ci/C2), and k, w are arbitrary constants.

I? 400

Figure 3: Soliton solution \vH(x,t)\ of the Schrodinger-Boussinesq system, (3.12b), for a = 1, p = 16, k = 1,w = -2, and pH = tanh-1(1/4).

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 Schrodinger-Boussinesq system. With the aid of Maple, we have assured the correctness of the obtained solutions by putting them back into the original equation. To illustrate the obtained solutions, the hyperbolic type of obtained solutions, (3.12a) and (3.12b), are attached as Figures 1, 2, and 3. We hope that they will be useful for further studies in applied sciences.

Acknowledgments

The authors would like to express their sincere thanks and gratitude to the reviewers for their valuable comments and suggestions for the improvement of this paper.

References

[1] C. Sulem and P.-L. Sulem, The Nonlinear Schrodinger Equation. Self-Focusing and Wave Collapse, vol. 139 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1999.

[2] T. Kato, "On nonlinear Schrodinger equations," Annales de l'Institut Henri Poincare. Physique Theorique, vol. 46, no. 1, pp. 113-129,1987.

[3] J. Boussinesq, "Theorie des ondes et des remous qui se propagent le long dun canal rectangulaire horizontal, en communiquant au liquide continu dans 21 ce canal des vitesses sensiblement pareilles de la surface au fond," Journal de Mathematiques Pures et Appliquees, vol. 17, no. 2, pp. 55-108,1872.

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

[5] V. Zakharov, "On stochastization of one-dimensional chains of nonlinear oscillators," Journal of Experimental and Theoretical Physics, vol. 38, pp. 110-108.

[6] F. Falk, E. Laedke, and K. Spatschek, "Stability of solitary-wave pulses in shape-memory alloys," Physical Review B, vol. 36, no. 6, pp. 3031-3041,1978.

[7] V. Makhankov, "On stationary solutions of Schrodinger equation with a self-consistent potential satisfying Boussinesqs equations," Physics Letters A, vol. 50, pp. 42-44,1974.

N. Yajima and J. Satsuma, "Soliton solutions in a diatomic lattice system," Progress of Theoretical Physics, vol. 62, no. 2, pp. 370-378,1979.

M. J. Ablowitz and H. Segur, Solitons and the Inverse Scattering Transform, vol. 4, SIAM, Philadelphia, Pa, USA, 1981.

M. R. Miurs, Backlund Transformation, Springer, Berlin, Germany, 1978.

R. Hirota, "Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons," Physical Review Letters, vol. 27, pp. 1192-1194,1971.

J. Weiss, M. Tabor, and G. Carnevale, "The Painleve property for partial differential equations," Journal of Mathematical Physics, vol. 24, no. 3, pp. 522-526,1983.

N. C. Freeman and J. J. C. Nimmo, "Soliton solutions of the Korteweg-de Vries and the Kadomtsev-Petviashvili equations: the Wronskian technique," Proceedings of the Royal Society of London Series A, vol. 389, no. 1797, pp. 319-329, 1983.

W. Malfliet and W. Hereman, "The tanh method—I. Exact solutions of nonlinear evolution and wave equations," Physica Scripta, vol. 54, no. 6, pp. 563-568,1996.

F. Xu, "Application of Exp-function method to symmetric regularized long wave (SRLW) equation," Physics Letters A, vol. 372, no. 3, pp. 252-257, 2008.

F. Tascan, A. Bekir, and M. Koparan, "Travelling wave solutions of nonlinear evolution equations by using the first integral method," Communications in Nonlinear Science and Numerical Simulation, vol. 14, pp. 1810-1815, 2009.

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

R. Abazari, "The (G'/G)-expansion method for Tzitzeica type nonlinear evolution equations," Mathematical and Computer Modelling, vol. 52, no. 9-10, pp. 1834-1845, 2010.

R. Abazari, "The (G'/G )-expansion method for the coupled Boussinesq equations," Procedia Engineering, vol. 10, pp. 2845-2850, 2011.

Copyright of Abstract & Applied Analysis 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.