# Traveling Wave Solutions of the Nonlinear(3+1)-Dimensional Kadomtsev-Petviashvili Equation Using the Two Variables (G′/G,1/G)-Expansion Method

Authors: E. M. E. Zayed; S. A. Hoda Ibrahim; M. A. M. Abdelaziz
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

The two variables(G′/G,1/G)-expansion method is proposed in this paper to construct new exact traveling wave solutions with parameters of the nonlinear (3+1)-dimensional Kadomtsev-Petviashvili equation. This method can be considered as an extension of the basic (G′/G)-expansion method obtained recently by Wang et al. When the parameters are replaced by special values, the well-known solitary wave solutions and the trigonometric periodic solutions of this equation were rediscovered from the traveling waves.


---

## Body

## 1. Introduction

In the recent years, investigations of exact solutions to nonlinear PDEs play an important role in the study of nonlinear physical phenomena. Many powerful methods have been presented, such as the inverse scattering transform method [1], the Hirota method [2], the truncated Painlevé expansion method [3–6], the Backlund transform method [7, 8], the exp-function method [9–14], the tanh function method [15–18], the Jacobi elliptic function expansion method [19–21], the original (G′/G)-expansion method [22–33], the two variables (G′/G,1/G)-expansion method [34, 35], and the first integral method [36]. The key idea of the original (G′/G)-expansion method is that the exact solutions of nonlinear PDEs can be expressed by a polynomial in one variable (G′/G) in which G=G(ξ) satisfies the second ordinary differential equation G′′(ξ)+λG′(ξ)+μG(ξ)=0, where λ and μ are constants. In this paper, we will use the two variables (G′/G,1/G)-expansion method, which can be considered as an extension of the original (G′/G)-expansion method. The key idea of the two variables (G′/G,1/G)-expansion method is that the exact traveling wave solutions of nonlinear PDEs can be expressed by a polynomial in the two variables (G′/G) and (1/G), in which G=G(ξ) satisfies a second order linear ODE, namely, G′′(ξ)+λG(ξ)=μ, where λ and μ are constants. The degree of this polynomial can be determined by considering the homogeneous balance between the highest order derivatives and nonlinear terms in the given nonlinear PDEs, while the coefficients of this polynomial can be obtained by solving a set of algebraic equations resulted from the process of using the method. According to Aslan [29], the two variables (G′/G,1/G)-expansion method becomes the basic (G′/G)-expansion method if μ=0 in (2.1) and bi=0 in (2.12). Recently, Li et al. [34] have applied the two variables (G′/G,1/G)-expansion method and determined the exact solutions of the Zakharov equations, while Zayed and abdelaziz [35] have applied this method to determine the exact solutions of the nonlinear KdV-mKdV equation.The objective of this paper is to apply the two variables(G′/G,1/G)-expansion method to find the exact traveling wave solutions of the following nonlinear (3+1)-dimensional Kadomtsev-Petviashvili equation:
(1.1)uxt+6(ux)2+6uuxx-uxxxx-uyy-uzz=0.
This equation describes the dynamics of solitons and nonlinear wave in plasma and superfluids. Recently, Zayed [24] has found the exact solutions of (1.1) using the original (G′/G)-expansion method, while Aslan [14] has discussed (1.1) using the exp-function method. Comparison between our results and that obtained in [14, 24] will be discussed later. The rest of this paper is organized as follows. In Section 2, the description of the two variables (G′/G,1/G)-expansion method is given. In Section 3, we apply this method to (1.1). In Section 4, conclusions are obtained.


## 2. Description of the Two Variables(G′/G,1/G)-Expansion Method

Before we describe the main steps of this method, we need the following remarks (see [34, 35]).Remark 2.1.
If we consider the second order linear ODE(2.1)G′′(ξ)+λG(ξ)=μ,
and set ϕ=G′/G and ψ=1/G, then we get
(2.2)ϕ′=-ϕ2+μψ-λ,ψ′=-ϕψ.Remark 2.2.
Ifλ<0, then the general solutions of  (2.1) is
(2.3)G(ξ)=A1sinh(ξ-λ)+A2cosh(ξ-λ)+μλ,
where A1 and A2 are arbitrary constants. Consequently, we have
(2.4)ψ2=-λλ2σ+μ2(ϕ2-2μψ+λ),
where σ=A12-A22.Remark 2.3.
Ifλ>0, then the general solutions of (2.1) is
(2.5)G(ξ)=A1sin(ξλ)+A2cos(ξλ)+μλ,
and hence
(2.6)ψ2=-λλ2σ-μ2(ϕ2-2μψ+λ).
where σ=A12+A22.Remark 2.4.
Ifλ=0, then the general solutions of (2.1) is
(2.7)G(ξ)=μ2ξ2+A1ξ+A2,
and hence
(2.8)ψ2=1A12-2μA2(ϕ2-2μψ),Suppose we have the following NLPDEs in the form:(2.9)F(u,ut,ux,uxx,utt,…)=0,
where F is a polynomial in u and its partial derivatives. In the following, we give the main steps of the two variables (G′/G,1/G)-expansion method [34, 35].Step 1.
The traveling wave variable(2.10)u(x,t)=u(ξ),ξ=x-Vt
reduces (2.9) to an ODE in the form
(2.11)P(u,u′,u′′,…)=0,
where V is a constant and P is a polynomial in u and its total derivatives, while {}′=d/dξ.Step 2.
Suppose that the solutions of (2.11) can be expressed by a polynomial in the two variables ϕ and ψ as follows:
(2.12)u(ξ)=∑i=0i=Naiϕi+∑i=1i=Nbiϕi-1ψ,
where ai(i=0,1,…,N) and bi(i=1,…,N) are constants to be determined later.Step 3.
Determine the positive integerN in (2.12) by using the homogeneous balance between the highest order derivatives and the nonlinear terms in (2.11).Step 4.
Substitute (2.12) into (2.11) along with (2.2) and (2.4), the left-hand side of (2.11) can be converted into a polynomial in ϕ and ψ, in which the degree of ψ is not longer than 1. Equating each coefficients of this polynomial to zero yields a system of algebraic equations which can be solved by using the Maple or Mathematica to get the values of ai,bi,V,μ,A1,A2, and λ where λ<0. Thus, we get the exact solutions in terms of the hyperbolic functions.Step 5.
Similar to Step4, substitute (2.12) into (2.11) along with (2.2) and (2.6) for λ>0 (or (2.2) and (2.8) for λ=0), we obtain the exact solutions of (2.11) expressed by trigonometric functions (or by rational functions), respectively.


## 3. An Application

In this section, we apply the method described in Section2, to find the exact traveling wave solutions of the nonlinear (3+1)-dimensional Kadomtsev-Petviashvili equation (1.1). To this end, we see that the traveling wave variables ξ=x+y+z-Vt reduce (1.1) to the following ODE:
(3.1)-(2+V)u′′+6(u′)2+6uu′′-u′′′′=0.
Balancing u′′′′ with uu′′ in (3.1) we get N=2. Consequently, we get
(3.2)u(ξ)=a0+a1ϕ(ξ)+a2ϕ2(ξ)+b1ψ(ξ)+b2ϕ(ξ)ψ(ξ),
where a0,a1,a2,b1, and b2 are constants to be determined later. There are three cases to be discussed as follows.Case 1.
Hyperbolic function solutions(λ<0).
Ifλ<0, substituting (3.2) into (3.1) and using (2.2) and (2.4), the left-hand side of (3.1) becomes a polynomial in ϕ and ψ. Setting the coefficients of this polynomial to zero yields a system of algebraic equations in a0,a1,a2,b1,b2,μ,σ, and λ which can be solved by using the Maple or Mathematica to find the following results:
(3.3)a0=a0,a1=0,a2=1,b1=-μ,b2=±-(λ2σ+μ2)λ,V=6a0-5λ-2.
From (2.3), (3.2), and (3.3), we deduce the traveling wave solution of  (1.1) as follows:
(3.4)u(ξ)=a0-μA1sinh(ξ-λ)+A2cosh(ξ-λ)+(μ/λ)-A1cosh(ξ-λ)+A2sinh(ξ-λ)(A1sinh(ξ-λ)+A2cosh(ξ-λ)+μ/λ)2×[A1λcosh(ξ-λ)+A2λsinh(ξ-λ)∓λ2σ+μ2],
where
(3.5)ξ=x+y+z-(6a0-5λ-2)t.
In particular, by setting A1=0,A2>0 and μ=0 in (3.4), we have the solitary solution
(3.6)u(ξ)=a0-λtanh(ξ-λ)[tanh(ξ-λ)∓isech(ξ-λ)],
where i=-1, while if A2=0,A1>0, and μ=0, then we have the solitary solution
(3.7)u(ξ)=a0-λcoth(ξ-λ)[coth(ξ-λ)∓csch(ξ-λ)].Case 2.
Trigonometric function solutions(λ>0).
Ifλ>0, substituting (3.2) into (3.1) and using (2.2) and (2.6), we get a polynomial in ϕ and ψ. Vanishing each coefficient of this polynomial to get the algebraic equations which can be solved by using the Maple or Mathematica to find the following results:
(3.8)a0=a0,a1=0,a2=1,b1=-μ,b2=±λ2σ-μ2λ,V=6a0-5λ-2.
From (2.5), (3.2), and (3.8), we deduce the traveling wave solution of (1.1) as follows:
(3.9)u(ξ)=a0-μA1sin(ξλ)+A2cos(ξλ)+(μ/λ)+A1cos(ξλ)-A2sin(ξλ)(A1sin(ξλ)+A2cos(ξλ)+μ/λ)2×[A1λcos(ξλ)-A2λsin(ξλ)±λ2σ-μ2],
where ξ has the same form (3.5).
In particular, by settingA1=0,A2>0, and μ=0 in (3.9), we have the periodic solution
(3.10)u(ξ)=a0+λtan(ξλ)[tan(ξλ)∓sec(ξλ)],
while if A2=0,A1>0, and μ=0, then we have the periodic solution
(3.11)u(ξ)=a0+λcot(ξλ)[cot(ξλ)±csc(ξλ)].Case 3.
Rational function solutions(λ=0).
Ifλ=0, substituting (3.2) into (3.1) and using (2.2) and (2.8), we get a polynomial in ϕ and ψ. Setting each coefficients of this polynomial to be zero to get the algebraic equations which can be solved by using the Maple or Mathematica to find the following results:
(3.12)a0=a0,a1=0,a2=1,b1=-μ,b2=±A12-2μA2,V=6a0-2.
From (2.7), (3.2), and (3.12), we deduce the traveling wave solution of (1.1) as follows:
(3.13)u(ξ)=a0-μ(μ/2)ξ2+A1ξ+A2+(μξ+A1)(μξ+A1±A12-2μA2)((μ/2)ξ2+A1ξ+A2)2,
where
(3.14)ξ=x+y+z-(6a0-2)t.Remark 3.1.
All solutions of this paper have been checked with Maple by putting them back into the original equation (1.1).


## 4. Conclusions

The two variables(G′/G,1/G)-expansion method has been used in this paper to discuss (1.1) and obtain the exact traveling wave solutions (3.4), (3.9), and (3.13) of Section 3. As the two parameters A1 and A2 take special values, we obtain the solitary wave solutions (3.6) and (3.7) and the trigonometric periodic solutions (3.10) and (3.11). On comparing these solutions with the result (11) of [14] obtained by Aslan using the exp-function method as well as the results (3.28)–(3.31) of [24] obtained by Zayed using the basic (G′/G)-expansion method, we conclude that all these solutions of (1.1) are different and satisfying that equation. The advantage of the two variables (G′/G,1/G)-expansion method over the basic (G′/G)-expansion method is that the first method is an extension of the second one.



---

*Source: 560531-2012-09-04.xml*