Citation: Huaji Cheng, Yanxia Hu. Exact solutions of the generalized (2+1)-dimensional BKP equation by the G'/G-expansion method and the first integral method[J]. AIMS Mathematics, 2017, 2(3): 562-579. doi: 10.3934/Math.2017.2.562
[1] | Weiping Gao, Yanxia Hu . The exact traveling wave solutions of a class of generalized Black-Scholes equation. AIMS Mathematics, 2017, 2(3): 385-399. doi: 10.3934/Math.2017.3.385 |
[2] | Mustafa Inc, Mamun Miah, Akher Chowdhury, Shahadat Ali, Hadi Rezazadeh, Mehmet Ali Akinlar, Yu-Ming Chu . New exact solutions for the Kaup-Kupershmidt equation. AIMS Mathematics, 2020, 5(6): 6726-6738. doi: 10.3934/math.2020432 |
[3] | Abdulah A. Alghamdi . Analytical discovery of dark soliton lattices in (2+1)-dimensional generalized fractional Kundu-Mukherjee-Naskar equation. AIMS Mathematics, 2024, 9(8): 23100-23127. doi: 10.3934/math.20241123 |
[4] | Yunmei Zhao, Yinghui He, Huizhang Yang . The two variable (φ/φ, 1/φ)-expansion method for solving the time-fractional partial differential equations. AIMS Mathematics, 2020, 5(5): 4121-4135. doi: 10.3934/math.2020264 |
[5] | Amjad Hussain, Muhammad Khubaib Zia, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar, Ilyas Khan . Lie analysis, conserved vectors, nonlinear self-adjoint classification and exact solutions of generalized $ \left(N+1\right) $-dimensional nonlinear Boussinesq equation. AIMS Mathematics, 2022, 7(7): 13139-13168. doi: 10.3934/math.2022725 |
[6] | Abeer S. Khalifa, Hamdy M. Ahmed, Niveen M. Badra, Jalil Manafian, Khaled H. Mahmoud, Kottakkaran Sooppy Nisar, Wafaa B. Rabie . Derivation of some solitary wave solutions for the (3+1)- dimensional pKP-BKP equation via the IME tanh function method. AIMS Mathematics, 2024, 9(10): 27704-27720. doi: 10.3934/math.20241345 |
[7] | M. TarikulIslam, M. AliAkbar, M. Abul Kalam Azad . Traveling wave solutions in closed form for some nonlinear fractional evolution equations related to conformable fractional derivative. AIMS Mathematics, 2018, 3(4): 625-646. doi: 10.3934/Math.2018.4.625 |
[8] | Ghazala Akram, Saima Arshed, Maasoomah Sadaf, Hajra Mariyam, Muhammad Nauman Aslam, Riaz Ahmad, Ilyas Khan, Jawaher Alzahrani . Abundant solitary wave solutions of Gardner's equation using three effective integration techniques. AIMS Mathematics, 2023, 8(4): 8171-8184. doi: 10.3934/math.2023413 |
[9] | M. Hafiz Uddin, M. Ali Akbar, Md. Ashrafuzzaman Khan, Md. Abdul Haque . New exact solitary wave solutions to the space-time fractional differential equations with conformable derivative. AIMS Mathematics, 2019, 4(2): 199-214. doi: 10.3934/math.2019.2.199 |
[10] | M. Ali Akbar, Norhashidah Hj. Mohd. Ali, M. Tarikul Islam . Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics. AIMS Mathematics, 2019, 4(3): 397-411. doi: 10.3934/math.2019.3.397 |
During the past decades, the exact solutions of nonlinear partial differential equations have been investigated by many authors. Meanwhile, many powerful methods have been proposed by them, such as Backlund transformation method [1], multiple exp-function method [2], homogeneous balance principle [3], tanh-sech method [4], G′/G-expansion method [5,6,7], the first integral method [8,9] and so on.
The G′/G-expansion method was first presented by Wang [5] which can be used to deal with all types of nonlinear evolution equations. The first integral method was first proposed by Feng [8] for obtaining the exact solutions of Burgers-KdV equation which is based on the ring theory of commutative algebra. The basic idea of the first integral method is to construct a first integral with polynomial coefficients of an explicit form to an equivalent autonomous planer system by using the division theorem. Both the G′/G-expansion method and the first integral method are powerful methods for computing the exact solutions of nonlinear partial differential equations. They are direct, elementary and effective algebraic methods.
In this paper, we consider the following generalized (2+1)-dimensional BKP equation [10]
{(wn)t+(wm)xxx+(wm)yyy+α(uw)x+β(vw)y=0,uy=wx,vx=wy, | (1.1) |
where α,β are arbitrary constants and α+β≠0,m,n are integers and m,n≥2. In [10], authors studied traveling wave solutions in the parameter space of this system by bifurcation theory of dynamical systems and they obtained some exact explicit parametric representations of periodic cusp wave solutions, solitary wave solutions and compacton solutions. In this paper, we continue to consider the problem of solving system (1.1) by using the G′/G-expansion method and the first integral method and we obtain the rational function solutions, periodic function solutions and the hyperbolic function solutions of (1.1) under some parametric conditions and the values of m, n in several cases.
Specially, when m=1,n=1,α=β=6, (1.1) becomes
{wt+wxxx+wyyy+6(uw)x+6(vw)y=0,uy=wx,vx=wy. |
It is the famous (2+1)-dimensional BKP equation which was introduced by Date et al. [11] and describes the processes of interaction of exponentially localized structures. It is one of a hierarchy of integrable systems emerging from a bilinear identity related to a Clifford algebra which is generated by two neutral fermion fields [12]. This equation has been studied by using many methods, such as the sine-cosine method [13], the G′/G-expansion method [6], the improved G′/G-expansion method [14] and so on.
The aim of this paper is to extract the exact solutions of the generalized (2+1)-dimensional BKP equation by using the G′/G-expansion method and the first integral method. The paper is arranged as follows: In section 2, we apply the G′/G-expansion method to this equation. In section 3, we apply the first integral method to solve this equation. In section 4, we give the conclusion of the paper.
We suppose the wave transformations
w(x,y,t)=w(ξ),u(x,y,t)=u(ξ),v(x,y,t)=v(ξ), ξ=k1x+l1y+λ1t | (2.1) |
where k1,l1,λ1 are constants. By using the wave transformations (2.1), (1.1) can be converted into ODEs
{λ1(wn)′+(k31+l31)(wm)′′′+αk1(u′w+uw′)+βl1(v′w+vw′)=0,l1u′=k1w′,k1v′=l1w′, | (2.2) |
where "′" is the derivative with respect to ξ. Integrating the second and third equation of (2.2) and neglecting integral constants, we obtain
{l1u=k1w,k1v=l1w. |
Substituting the above equations into the first equation of (2.2) and integrating it, then it becomes
λ1wn+(k31+l31)(wm)′′+(αk21l1+βl21k1)w2=g, | (2.3) |
where g is an integral constant. We assume that (2.3) has the following formal solutions [7,15]:
w(ξ)=D(G′G)N, D≠0, | (2.4) |
where D is a constant to be determined later. N is determined by balancing the linear term of the highest order derivatives with the highest order nonlinear term in (2.3) and G satisfies a second order constant coefficient ODE which is
G″(ξ)+λG′(ξ)+μG(ξ)=0, | (2.5) |
where λ, μ are constants and will be determined later. Next, we will obtain the exact solutions of (1.1) by considering the values of m and n in several cases.
Balancing (wm)′′ with wn of (2.3), we have mN+2=nN, i.e., N=2/(n−m). Thus, we assume
w(ξ)=D1(G′G)2n−m, D1≠0 | (2.6) |
where D1 is a constant to be determined later. Then, we have
wn=Dn1(G′G)2nn−m, w2=D21(G′G)4n−m,(wm)″=2mn−mDm1[(2mn−m+1)(G′G)2mn−m+2+(4mn−m+1)λ(G′G)2mn−m+1+2mn−m(2μ+λ2)(G′G)2mn−m+(4mn−m−1)λμ(G′G)2mn−m−1+(2mn−m−1)μ2(G′G)2mn−m−2]. |
Substituting the above formulas into (2.3) and collecting all terms with the same order of G′/G together, we can convert the left-hand side of (2.3) into a polynomial in G′/G. Then, setting each coefficient of each polynomial to zero, we can derive a set of algebraic equation for λ,μ and D1:
(G′G)2mn−m+2 coeff:
(k31+l31)(2mn−m+1)2mn−mDm1+λ1Dn1=0, | (2.7) |
(G′G)2mn−m+1 coeff:
(k31+l31)(4mn−m+1)2mn−mλDm1=0. | (2.8) |
Here, we need to consider the value of 4/(n−m) in the following cases:
Case 1. 4n−m=2mn−m−1
(G′G)2mn−m coeff:
(k31+l31)(2mn−m)2(2μ+λ2)Dm1=0, | (2.9) |
(G′G)2mn−m−1 coeff:
(k31+l31)(4mn−m−1)2mn−mλμDm1+(αk21l1+βl21k1)D21=0, | (2.10) |
(G′G)2mn−m−2 coeff:
(k31+l31)(2mn−m−1)2mn−mμ2Dm1=0. | (2.11) |
Solving the set of (2.7)-(2.11), we obtain
λ=μ=0, g=0, αk21l1+βl21k1=0, D1=(−(k31+l31)(2mn−m+1)2mn−mλ1)1/(n−m). | (2.12) |
Case 2. 4n−m=2mn−m−2
(G′G)2mn−m coeff:
(k31+l31)(2mn−m)2(2μ+λ2)Dm1=0, | (2.13) |
(G′G)2mn−m−1 coeff:
(k31+l31)(4mn−m−1)2mn−mλμDm1=0, | (2.14) |
(G′G)2mn−m−2 coeff:
(k31+l31)(2mn−m−1)2mn−mμ2Dm1+(αk21l1+βl21k1)D21=0. | (2.15) |
Solving the set of (2.7)-(2.8) and (2.13)-(2.15), we get the same results as those of Case 1.
Case 3. 4n−m≠2mn−m−1 } and 4n−m≠2mn−m−2
(G′G)2mn−m coeff:
(k31+l31)(2mn−m)2(2μ+λ2)Dm1=0, | (2.16) |
(G′G)2mn−m−1 coeff:
(k31+l31)(4mn−m−1)2mn−mλμDm1=0, | (2.17) |
(G′G)2mn−m−2 coeff:
(k31+l31)(2mn−m−1)2mn−mμ2Dm1=0, | (2.18) |
(G′G)4n−m coeff:
(αk21l1+βl21k1)D21=0. | (2.19) |
Solving the set of (2.7)-(2.8) and (2.16)-(2.19), we obtain the same results as those of former cases. Substituting (2.12) into (2.5) and (2.6), then, we can get the rational function solutions
w(x,y,t)=(−(k31+l31)(2mn−m+1)2mn−mλ1)1n−m(C1C1(k1x+l1y+λ1t)+C2)2n−m,v(x,y,t)=l1k1(−(k31+l31)(2mn−m+1)2mn−mλ1)1n−m(C1C1(k1x+l1y+λ1t)+C2)2n−m,u(x,y,t)=k1l1(−(k31+l31)(2mn−m+1)2mn−mλ1)1n−m(C1C1(k1x+l1y+λ1t)+C2)2n−m, |
where C1,C2 are arbitrary constants and αk31+βl31=0.
(2.3) becomes
λ1wn+(k31+l31)(w2)′′+(αk21l1+βl21k1)w2=g. | (2.20) |
Balancing (w2)′′ with wn, we have N=2/(n−2). Thus, (2.20) has the following formal solutions
w(ξ)=D2(G′G)2n−2, D2≠0, | (2.21) |
where D2 is a constant to be determined later and G satisfies (2.5). Similarly, we can get a set of algebraic equations:
(G′G)4n−2+2 coeff:
(k31+l31)(4n−2+1)4n−2D22+λ1Dn2=0, | (2.22) |
(G′G)4n−2+1 coeff:
(k31+l31)(8n−2+1)4n−2λD22=0, | (2.23) |
(G′G)4n−2 coeff:
(k31+l31)(4n−2)2(2μ+λ2)D22+(αk21l1+βl21k1)D22=0, | (2.24) |
Ⅰ. The case g=0
(G′G)4n−2−1 coeff:
(k31+l31)(8n−2−1)4n−2λμD22=0, | (2.25) |
(G′G)4n−2−2 coeff:
(k31+l31)(4n−2−1)4n−2μ2D22=0. | (2.26) |
Solving that set of (2.22)-(2.26), we obtain
λ=μ=0, αk21l1+βl21k1=0, D2=(−(k31+l31)(4n−2+1)4n−2λ1)1/(n−2). | (2.27) |
Specially, when 4n−2−1=0, i.e., n=6, we obtain
λ=0, μ=−αk21l1+βl21k12(k31+l31), D2=(−2(k31+l31)λ1)1/4. | (2.28) |
Substituting (2.27) into (2.5) and (2.21), then, we can get the rational function solutions
w(x,y,t)=(−(k31+l31)(4n−2+1)4n−2λ1)1n−2(C3C3(k1x+l1y+λ1t)+C4)2n−2,v(x,y,t)=l1k1(−(k31+l31)(4n−2+1)4n−2λ1)1n−2(C3C3(k1x+l1y+λ1t)+C4)2n−2,u(x,y,t)=k1l1(−(k31+l31)(4n−2+1)4n−2λ1)1n−2(C3C3(k1x+l1y+λ1t)+C4)2n−2, |
where C3,C4 are arbitrary constants and αk31+βl31=0. Substituting (2.28) into (2.5) and (2.21), then, we have
G″+(−αk21l1+βl21k12(k31+l31))G=0. |
Case 1. αk21l1+βl21k12(k31+l31)>0
We obtain the hyperbolic function solutions
w(x,y,t)=(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{C5sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C5cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)}1/2,v(x,y,t)=l1k1(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{C5sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C5cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)}1/2,u(x,y,t)=k1l1(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{C5sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C5cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C6sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)}1/2, |
where C5,C6 are arbitrary constants.
Case 2. αk21l1+βl21k12(k31+l31)<0
We obtain the hyperbolic function solutions
w(x,y,t)=(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{−C7sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C7cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)}1/2,v(x,y,t)=l1k1(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{−C7sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C7cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)}1/2,u(x,y,t)=k1l1(−(αk21l1+βl21k1)22λ1(k31+l31))1/4{−C7sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C7cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C8sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)}1/2, |
where C7,C8 are arbitrary constants.
Case 3. αk21l1+βl21k12(k31+l31)=0
We obtain the rational function solutions
w(x,y,t)=(−2(k31+l31)λ1)1/4{C9C9(k1x+l1y+λ1t)+C10}1/2,v(x,y,t)=l1k1(−2(k31+l31)λ1)1/4{C9C9(k1x+l1y+λ1t)+C10}1/2,u(x,y,t)=k1l1(−2(k31+l31)λ1)1/4{C9C9(k1x+l1y+λ1t)+C10}1/2, |
where C9,C10 are arbitrary constants.
Ⅱ. The case g≠0
When 4n−2−2=0, i.e, n=4.
(G′G) coeff:
6(k31+l31)λμD22=0, | (2.29) |
(G′G)0 coeff:
2(k31+l31)μ2D22=g. | (2.30) |
Solving the set of (2.22)-(2.24) and (2.29)-(2.30), we obtain
λ=0, μ=−αk21l1+βl21k12(k31+l31), D2=(−6(k31+l31)λ1)1/2, αk21l1+βl21k1≠0, −3(αk21l1+βl21k1)2λ1=g. | (2.31) |
Similarly, we can obtain the hyperbolic function solutions and trigonometric function solutions
Case 1. αk21l1+βl21k12(k31+l31)>0
w(x,y,t)=(−3(αk21l1+βl21k1)λ1)1/2{C11sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C11cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)},v(x,y,t)=l1k1(−3(αk21l1+βl21k1)λ1)1/2{C11sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C11cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)},u(x,y,t)=k1l1(−3(αk21l1+βl21k1)λ1)1/2{C11sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)C11cosh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)+C12sinh(αk21l1+βl21k12(k31+l31))1/2(k1x+l1y+λ1t)}, | (2.32) |
where C11,C12 are arbitrary constants and λ1(k31+l31)<0.
Case 2. αk21l1+βl21k12(k31+l31)<0
w(x,y,t)=(−3(αk21l1+βl21k1)λ1)1/2{−C13sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C13cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)},v(x,y,t)=l1k1(−3(αk21l1+βl21k1)λ1)1/2{−C13sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C13cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)},u(x,y,t)=k1l1(−3(αk21l1+βl21k1)λ1)1/2{−C13sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)C13cos(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)+C14sin(−(αk21l1+βl21k1)2(k31+l31))1/2(k1x+l1y+λ1t)}, | (2.33) |
where C13,C14 are arbitrary constants and λ1(k31+l31)<0.
(2.3) becomes
(λ1+αk21l1+βl21k1)w2+(k31+l31)(wm)′′=g. | (2.34) |
Balancing (wm)′′ with w2, we have N=2/(2−m). Thus, (2.34) has the following formal solution
w(ξ)=D3(G′G)22−m, D3≠0, | (2.35) |
where D3 is a constant to be determined later and G satisfies (2.5). Similarly, we can get a set of algebraic equations:
(G′G)2m2−m+2 coeff:
(k31+l31)(2m2−m+1)(2m2−m)Dm3+(λ1+αk21l1+βl21k1)D23=0, |
(G′G)2m2−m+1 coeff:
(k31+l31)(4m2−m+1)2m2−mλDm3=0, |
(G′G)2m2−m coeff:
(k31+l31)(2m2−m)2(2μ+λ2)Dm3=0, |
(G′G)2m2−m−1 coeff:
(k31+l31)(4m2−m−1)2m2−mλμDm3=0, |
(G′G)2m2−m−2 coeff:
(k31+l31)(2m2−m−1)2m2−mμ2Dm3=0. |
Solving the above algebraic equations, we obtain
λ=μ=0, g=0, D3=(−(k31+l31)(2m2−m+1)2m2−mλ1+αk21l1+βl21k1)1/(2−m). | (2.36) |
Substituting (2.36) into (2.5) and (2.35), then, when m≠n=2, we have the rational function solutions
w(x,y,t)=(−(k31+l31)(2m2−m+1)2m2−mλ1+αk21l1+βl21k1)12−m(C15C15(k1x+l1y+λ1t)+C16)22−m,v(x,y,t)=l1k1(−(k31+l31)(2m2−m+1)2m2−mλ1+αk21l1+βl21k1)12−m(C15C15(k1x+l1y+λ1t)+C16)22−m,u(x,y,t)=k1l1(−(k31+l31)(2m2−m+1)2m2−mλ1+αk21l1+βl21k1)12−m(C15C15(k1x+l1y+λ1t)+C16)22−m, |
where C15,C16 are arbitrary constants and λ1+αk21l1+βl21k1≠0.
Now, (2.3) can be converted into a second order ODE
(λ1+αk21l1+βl21k1)w2+(k31+l31)(w2)′′=g. | (2.37) |
Obviously, the characteristic equation of (2.37) is r2+(λ1+αk21l1+βl21k1k31+l31)=0, r is the characteristic value.
Case 1. λ1+αk21l1+βl21k1k31+l31<0
We can obtain the exact solution
w(x,y,t)=(C17e(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+C18e−(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+gλ1+αk21l1+βl21k1)1/2,v(x,y,t)=l1k1(C17e(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+C18e−(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+gλ1+αk21l1+βl21k1)1/2,u(x,y,t)=k1l1(C17e(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+C18e−(−λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)+gλ1+αk21l1+βl21k1)1/2, | (2.38) |
where C17,C18 are arbitrary constants.
Case 2. λ1+αk21l1+βl21k1k31+l31>0
We can obtain the periodic function solutions
w(x,y,t)=(C19cos[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+C20sin[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+gλ1+αk21l1+βl21k1)1/2,v(x,y,t)=l1k1(C19cos[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+C20sin[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+gλ1+αk21l1+βl21k1)1/2,u(x,y,t)=k1l1(C19cos[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+C20sin[(λ1+αk21l1+βl21k1k31+l31)1/2(k1x+l1y+λ1t)]+gλ1+αk21l1+βl21k1)1/2, | (2.39) |
where C19,C20 are arbitrary constants.
Case 3. λ1+αk21l1+βl21k1k31+l31=0
We can obtain the rational function solutions
w(x,y,t)=(g2(k31+l31)(k1x+l1y+λ1t)2+C21(k1x+l1y+λ1t)+C22)1/2,v(x,y,t)=l1k1(g2(k31+l31)(k1x+l1y+λ1t)2+C21(k1x+l1y+λ1t)+C22)1/2,u(x,y,t)=k1l1(g2(k31+l31)(k1x+l1y+λ1t)2+C21(k1x+l1y+λ1t)+C22)1/2, | (2.40) |
where C21,C22 are arbitrary constants.
For simplicity, we let g=0 and propose a transformation w=φ2n−m. Then, (2.3) is converted to
λ1φ4+(k31+l31)(2mn−m−1)2mn−m(φ′)2+(k31+l31)2mn−mφφ″+(αk21l1+βl21k1)φ2−2m−4n−m=0. | (3.1) |
Let x=φ,y=dφdξ, thus (3.1) is equivalent to the two dimensional autonomous system
{x′=y,y′=−(λ1x4+(αk21l1+βl21k1)x2−2m−4n−m+(k31+l31)(2mn−m−1)2mn−my2(k31+l31)2mn−mx). | (3.2) |
Making the transformation dη=dξ(k31+l31)2mn−mx, then, (3.2) becomes
{dxdη=(k31+l31)2mn−mxy,dydη=−(λ1x4+(αk21l1+βl21k1)x2−2m−4n−m+(k31+l31)(2mn−m−1)2mn−my2). | (3.3) |
Then, we will apply the Division Theorem to seek the first integral of system (3.3). Suppose that x=x(η),y=y(η) are the nontrivial solutions to (3.3), and p(x,y)=∑Mi=0ai(x)yi is an irreducible polynomial in C[x,y], where ai(x)(i=0,1...,M) are polynomials of x and ai(x)≠0. Let p(x(η),y(η))=0 be the first integral to system (3.3). dpdη is a polynomial in x, y and dpdη|(3.3)=0. According to the Division Theorem, there exists a polynomial g(x)+h(x)y in C[x,y], such that
dpdη|(3.3)=(∂p∂xdxdη+∂p∂ydydη)|(3.3)=M∑i=0[a′i(x)yi⋅(k31+l31)2mn−mxy]−2∑i=0[iai(x)yi−1(λ1x4+(αk21l1+βl21k1)x2−2m−4n−m+(k31+l31)(2mn−m−1)2mn−my2)]=[g(x)+h(x)y][M∑i=0ai(x)yi]. | (3.4) |
Here, let M=1, thus, p(x,y)=a0(x)+a1(x)y. By comparing with the coefficients of yi of both sides of (3.4), we have
(k31+l31)2mn−mxa′1(x)=h(x)a1(x)+(k31+l31)(2mn−m−1)2mn−ma1(x), | (3.5) |
(k31+l31)2mn−mxa′0(x)=g(x)a1(x)+h(x)a0(x), | (3.6) |
g(x)a0(x)=−(λ1x4+(αk21l1+βl21k1)x2−2m−4n−m)a1(x). | (3.7) |
Since ai(x)(i=0,1) are polynomials, then from (3.5), we deduce that h(x)=−(k31+l31)(2mn−m−1)2mn−m and a1(x) is a constant. For simplicity, take a1(x)=1. Balancing the degrees of g(x) and a0(x), we conclude that deg(g(x))=deg(a0(x)). Then, we derive deg(g(x))=deg(a0(x))=j,(j∈Z+,j≥2).
When 2−2m−4n−m=4 (n=2) and deg(g(x))=deg(a0(x))=2, we suppose that
g(x)=A0+A1x+A2x2,a0(x)=B0+B1x+B2x2, (A2≠0,B2≠0), | (3.8) |
where Ai,Bi,(i=0,1,2) are all constants to be determined. Substituting (3.8) into (3.6), we obtain
g(x)=(k31+l31)2m2−m[(2m2−m−1)B0+2m2−mB1x+(2m2−m+1)B2x2]. |
Substituting a0(x),a1(x) and g(x) into (3.7), and setting all the coefficients of powers x to be zero, we can get a system of nonlinear algebraic equations. After solving it, we can get the following solutions
B0=B1=0, B2=±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2. | (3.9) |
Using the conditions (3.9) in p(x,y)=a0(x)+a1(x)y=0, we obtain
y±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2x2=0. | (3.10) |
Combining (3.3) with (3.10), we find
dxdη=±(k31+l31)2mn−m(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2x3. |
Thus, (3.10) can be reduced to
dφdξ=±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2φ2. |
Then, we have
φ(ξ)=[±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2ξ+C23]−1. |
Thus, we can have the rational function solutions
w(x,y,t)=[±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2(k1x+l1y+λ1t)+C23]2/(m−2),v(x,y,t)=l1k1[±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2(k1x+l1y+λ1t)+C23]2/(m−2),u(x,y,t)=k1l1[±(−λ1+αk21l1+βl21k1(k31+l31)(2m2−m+1)2m2−m)1/2(k1x+l1y+λ1t)+C23]2/(m−2), |
where C23 is an arbitrary constant and k31+l31≠0.
Remark 1: When deg(g(x))=deg(a0(x))=2 and 2−2m−4n−m=i,(i∈Z,i<4), there is no solution for them by using the method as that of 2−2m−4n−m=4.
Remark 2: When deg(g(x))=deg(a0(x))=j,(j∈Z,j>2), there is no exact solution of (1.1) by using the method as that of deg(g(x))=deg(a0(x))=2.
Similarly, we propose a transformation denoted by w=ϕ22−m. Then, (2.3) can be converted to
λ1ϕ2+(αk21l1+βl21k1)ϕ4+(k31+l31)(2m2−m−1)2m2−m(ϕ′)2+(k31+l31)2m2−mϕϕ″−gϕ2−2m2−m=0. | (3.11) |
Let x=ϕ,y=dϕdξ, thus (3.11) is equivalent to the two dimensional autonomous system
{x′=y,y′=−(λ1x2+(αk21l1+βl21k1)x4−gx2−2m2−m+(k31+l31)(2m2−m−1)2m2−my2(k31+l31)2m2−mx). | (3.12) |
Making the transformation dη=dξ(k31+l31)2mn−mx, then, (3.12) becomes
{dxdη=(k31+l31)2m2−mxy,dydη=−(λ1x2+(αk21l1+βl21k1)x4−gx2−2m2−m+(k31+l31)(2m2−m−1)2m2−my2). | (3.13) |
Similarly, let M=1, we have
(k31+l31)2m2−mxa′1(x)=h(x)a1(x)+(k31+l31)(2m2−m−1)2m2−ma1(x), | (3.14) |
(k31+l31)2m2−mxa′0(x)=g(x)a1(x)+h(x)a0(x), | (3.15) |
g(x)a0(x)=−(λ1x2+(αk21l1+βl21k1)x4−gx2−2m2−m)a1(x). | (3.16) |
According to m≠n, we have h(x)=−(k31+l31)(2m2−m−1)2m2−m, a1(x)=1 and deg(g(x))=deg(a0(x))=j,(j∈Z+,j≥2). Considering all cases, only when deg(g(x))=deg(a0(x))=3, i.e., 2−2m2−m=6 (m=n=4), there exist solutions of (1). We suppose that
g(x)=a0+a1x+a2x2+a3x3,a0(x)=b0+b1x+b2x2+b3x3, (a3≠0,b3≠0), | (3.17) |
where ai,bi,(i=0,1,2,3) are all constants to be determined. Substituting (3.17) into (3.15), we obtain
g(x)=4(k31+l31)(5b0+4b1x+3b2x2+2b3x3). |
Substituting a0(x),a1(x) and g(x) into (3.16), and setting all the coefficients of powers x to be zero, we have
b0=b2=0, 16(k31+l31)b21=−λ1, 8(k31+l31)b23=g, 24(k31+l31)b1b3=−(αk21l1+βl21k1). |
Solving it, we find
b0=b2=0, b1=±√−λ116(k31+l31), b3=±√g8(k31+l31), −9λ1g=2(αk21l1+βl21k1)2. | (3.18) |
Using the conditions (3.18) in p(x,y)=a0(x)+a1(x)y=0, we obtain
y=±√−λ116(k31+l31)x±√g8(k31+l31)x3. | (3.19) |
Then, (3.19) can be reduced to
dϕdξ=±√−λ116(k31+l31)ξ±√g8(k31+l31)ξ3. | (3.20) |
Solving (3.20), we obtain
ϕ(ξ)=±(±√−2gλ1+C24e±12√−λ1k31+l31ξ)−1/2. |
Thus, we can have the exact solution
w(x,y,t)=±(±√−2gλ1+C24e±12√−λ1k31+l31(k1x+l1y+λ1t))1/2,v(x,y,t)=±l1k1(±√−2gλ1+C24e±12√−λ1k31+l31(k1x+l1y+λ1t))1/2,u(x,y,t)=±k1l1(±√−2gλ1+C24e±12√−λ1k31+l31(k1x+l1y+λ1t))1/2, |
where C24 is an arbitrary constant and αk31+βl31≠0, λ1(k31+l31)<0.
This paper considered the generalized (2+1)-dimensional BKP equation, by the aid of the G′/G-expansion method and the first integral method. Rational function solutions, periodic function solutions and hyperbolic function solutions are obtained under some parametric conditions and the values of m and n in several cases. In [10], authors gave some exact solutions of system (1.1) under some parametric conditions by using the bifurcation theory of dynamical systems. Here, we make a simple comparison:
1. When m=2,n=3,g=0, in [10], authors gave the exact solution (3.20) in P2443 under the parametric conditions α+β<0,c<0 and in this paper, we get w(x,y,t)=−20(k31+l31)λ1(C1C1(k1x+l1y+λ1t)+C2)2 under the parametric condition αk31+βl31=0.
2. When m=2,n=2(k+1),(k∈Z+),g=0, in [10], authors gave the solitary wave solutions (3.9) in P2441 under the parametric conditions α+β<0,c<0 and in this paper, we get w(x,y,t)=(√−2(k+2)(k31+l31)λ1k2C1C1(k1x+l1y+λ1t)+C2)1/k under the parametric condition αk31+βl31=0.
3. When m=3,n=4,g=0, in [10], authors gave the compacton solution (3.23) in P2443 under the parametric conditions α+β<0,c<0 and in this paper, we get w(x,y,t)=−42(k31+l31)λ1(C1C1(k1x+l1y+λ1t)+C2)2 under the parametric condition αk31+βl31=0.
4. When m=3,n=5,g=0, in [10], authors gave the exact solution (3.27) in P2443 under the parametric conditions c<0 and in this paper, we get w(x,y,t)=√−12(k31+l31)λ1C1C1(k1x+l1y+λ1t)+C2 under the parametric condition αk31+βl31=0.
5. When m=4,n=6,g=0, in [10], authors gave the periodic cusp wave solutions (3.6) in P2441 under the parametric conditions α+β<0,c<0 and in this paper, we get w(x,y,t)=√−20(k31+l31)λ1C1C1(k1x+l1y+λ1t)+C2 under the parametric condition αk31+βl31=0.
6. When m=4,n=2k+1,(k∈Z+),g=0, in [10], authors gave the exact solutions (3.17) and (3.18) in P2442 under the parametric conditions α+β<0,c<0 and in this paper, we get w(x,y,t)=(−8(2k+5)(k31+l31)λ1(2k−3)2)12k−3(C1C1(k1x+l1y+λ1t)+C2)22k−3 under the parametric condition αk31+βl31=0.
7. When m=2,n=4,g≠0, in [10], authors gave the exact solutions (3.30) in P2443 under the parametric conditions g<(α+β)24c,c>0,g>0,α+β>0, (3.33) in P2444 under the parametric conditions g<(α+β)24c,c<0,g>0,α+β>0 or g>(α+β)24c,c<0,g>0,α+β<0 and (3.41), (3.43) in P2445 under the parametric conditions g>(α+β)24c,c<0,g<0,α+β<0 and in this paper, we get (2.32) under the parametric conditions λ1(k31+l31)<0, αk21l1+βl21k12(k31+l31)>0 and (2.33) under the parametric conditions λ1(k31+l31)<0, αk21l1+βl21k12(k31+l31)<0.
8. When m=2,n=2,g≠0, in [10], authors gave the exact solutions (3.36) and (3.38) in P2444 under the parametric conditions α+β−c>0,g>0 and in this paper, we get the exact solutions (2.38) under the parametric conditions λ1+αk21l1+βl21k1k31+l31<0, (2.39) under the parametric conditions λ1+αk21l1+βl21k1k31+l31>0 and (2.40) under the parametric conditions λ1+αk21l1+βl21k1k31+l31=0.
In addition, when let m,n be other values, we have got other exact solutions of (1.1) under some parametric conditions that haven't been given in [10]. Certainly, system (1.1) should be studied further, which will be left to a further discussion.
All authors declare no conflicts of interest in this paper.
[1] | M. R. Miura, Backlund Transformation, New York: Springer-Verlag, Berlin, 1978. |
[2] | W. Ma, T. Huang, and Y. Zhang, A multiple exp-function method for nonlinear differential equations and its application, Phys. Scr., 82 (2010), 065003. |
[3] | W. Ma and J. H. Lee, A transfortiom rational function method and exact solutions to (3+1)-dimensional Jimbo-Miwa equation, Choas Solitons Fractals, 42 (2009), 1356-1363. |
[4] | A. M. Wazwaz, Two reliable methods for solving variants of the KdV equation with compact and noncompact structures, Choas Solitons Fractals, 28 (2006), 454-462. |
[5] | M. L. Wang, X. Z. Li, and J. L. Zhang, The G'/G-expansion method and traveling wave solutions of nonlinear evolutions in mathematical physics, Phys. Lett., 372 (2008), 417-423. |
[6] | E. M. E. Zayed and Khaled A. Gepreel, Some applications of the (G'/G)-expansion method to nonlinear partial differential equations, Appl. Math. Comput., 212 (2009), 1-13. |
[7] | H. Q. Zhang, New application of the G0'/G-expansion method, Commun Nolinear Sci. Numer. Simul., 14 (2009), 3220-3225. |
[8] | Z. S. Feng, The first integral method to study the Burgers-Korteweg-de Vries equation, J. Phys. A., 35 (2002), 343-349. |
[9] | N. Taghizadeh and M. Mirzazadeh, Exact solutions of some nonlinear evolution equations via the first integral method, Ain Shams Engineering Journal, 4 (2013), 493-499. |
[10] | Y. G. Xie, B. W. Zhou, and S. Q. Tang, Bifurcations of traveling wave soluions for the generalized (2+1)-dimensional Boussinesq-Kadomtesv-Petviashvili equation, Appl. Math. Comput., 217 (2010), 2433-2447. |
[11] | E. Date, M. Jimbo, M. Kashiwara, and T. Miwa, Transformation groups for soliton equations. Ⅳ. A new hierarchy of soliton equations of KP type, Physica D., 4 (1982), 343-365. |
[12] | H. C. Ma, Y. Wang, and Z. Y. Qin, New exact complex traveling wave solutions for (2+1)-dimensional BKP equation, Appl. Math. Comput., 208 (2009), 564-568. |
[13] | F. Tascan and A. Bekir, Analytic solutions of the (2+1)-dimensional nonlinear evolution equations using the sine-cosine method, Appl. Math. Comput., 215 (2009), 3134-3139. |
[14] | H. Q. Zhang, A note on exact complex travelling wave solutions for (2+1)-dimensional B-type Kadomtsev-Petviashvili equation, Appl. Math. Comput., 216 (2010), 2771-2777. |
[15] | M. Mirzazadeh and M. Eslami, Topological solitons of resonant nonlinear Schrödinger's equation with dual-power law nonlinearity by G'/G-expansion technique, Optik, 125 (2014), 5480-5489. |
1. | Gui Mu, Yan Zhu, Tingfu Feng, New Localized Structure for (2+1) Dimensional Boussinesq-Kadomtsev-Petviashvili Equation, 2022, 10, 2227-7390, 2634, 10.3390/math10152634 | |
2. | Sanjaya K. Mohanty, Sachin Kumar, Apul N. Dev, Manoj Kr. Deka, Dmitry V. Churikov, Oleg V. Kravchenko, An efficient technique of G′G–expansion method for modified KdV and Burgers equations with variable coefficients, 2022, 37, 22113797, 105504, 10.1016/j.rinp.2022.105504 |