Exact Solutions for (2 + 1)-Dimensional KdV-Calogero-Bogoyavlenkskii-Schiff Equation via Symbolic Computation

Yan Li^1,2, Temuer Chaolu^2*

¹College of Sciences, Inner Mongolia University of Technology, Hohhot, China.
²College of Art and Sciences, Shanghai Maritime University, Shanghai, China.

DOI: 10.4236/jamp.2020.82015                

Abstract

This paper constructs exact solutions for the (2 + 1)-dimensional KdV-Calogero-Bogoyavlenkskii-Schiff equation with the help of symbolic computation. By means of the truncated Painlev expansion, the (2 + 1)-dimensional KdV-Calogero-Bogoyavlenkskii-Schiff equation can be written as a trilinear equation, through the trilinear-linear equation, we can obtain the explicit representation of exact solutions for the (2 + 1)-dimensional KdV-Calogero-Bogoyavlenkskii-Schiff equation. We have depicted the profiles of the exact solutions by presenting their three-dimensional plots and the corresponding density plots.

Keywords

(2 + 1)-Dimensional KdV-Calogero-Bogoyavlenkskii-Schiff Equation, Trilinear Equation, Exact Solutions

1. Introduction

Over the past few decades, the research of the nonlinear evolution equations (NLEEs) is flourishing because of the rich findings of these equations. NLEEs display many interesting nonlinear dynamic behaviors, such as plasma physics, optical systems, ocean, superfluids, hydrodynamics and other nonlinear fields [1] [2] [3] [4]. The solutions of NLEEs can provide much physical information and more insight into the physical aspects and then lead to further applications, so seeking exact solutions is an important problem in nonlinear science. In fact, there are a variety of powerful methods have been used to solve the NLEEs. For instance, the Darboux transformation [5], the Lie group method [6], the Hirota bilinear method [7], the inverse scattering method [8] and other methods. Recent study for constructing lump solution by taking the function f in its bilinear form as a positive quadratic function based on Hirota bilinear method was obtained in [9]. Then, by using this method, more and more the lump solutions and the interaction solutions have been perfectly obtained in NLEEs equations [10] - [26].

In this paper, we focus on the following (2 + 1)-dimensional KdV-Calogero-Bogoyavlenkskii-Schiff (KdV-CBS) equation [27] [28] [29].

$4u_t-\alpha \left(4u{u}_y+2u_x{\partial }_x^{-1}{u}_y+2u_{xxy}\right)-\beta \left(6u{u}_x+u_{xxx}\right)=0,$ (1)

which with two arbitrary constant $\alpha$ and $\beta$, Equation (1) can be reduced to some other equations with physical meanings. If we setting $\alpha \ne 0,\beta =0$, Equation (1) is reduced to Calogero-Bogoyavlenskii-Schiff equation [30].

$4u_t-\alpha \left(4u{u}_y+2u_x{\partial }_x^{-1}{u}_y+u_{xxy}\right)=0.$ (2)

If setting $\left\{\alpha =0,\beta \ne 0\right\}$, Equation (1) is famous KdV equation

$4u_t-\beta \left(u_{xxx}+6u{u}_x\right)=0.$ (3)

Let $v_x=u_y$, Equation (1) is transformed into the following system

$u_t-\alpha \left(4u{u}_y+2u_{x}v+u_{xxy}\right)-\beta \left(6u{u}_x+u_{xxx}\right)=0,$ (4a)

$v_x-u_y=0.$ (4b)

By applying Painlevé truncated extension, some types of explicit solutions (1) were studied in [31]. A direct bilinear Bäcklund transformation on the basis of quadrilinear representation was constructed in [32] [33]. Ref. [34] obtained soliton solutions, quasiperiodic wave solutions of the Equation (1) based on the Riemann-Bäcklunk method. Interaction solutions of Equation (1) were obtained under consistent Riccati expansion in [29].

The outline of this paper is as follows. In Section 2, we convert the original KdV-CBS equation to a trilinear equation by using the truncated Painlevé expansion. In Section 3, we derived exact solutions for the KdV-CBS equation by taking function f as a positive quadratic function; another kind of exact solutions were got by taking function f as a positive quadratic and an exponential function; by taking function f as two exponential terms and the quadratic function, we obtained three kinds of exact solutions respectively via different parameters. A short conclusion is included in the last section.

2. Trilinear Equation

Based on the Painlevé analysis proposed in Ref. [4], Equation (4a) and Equation (4b) possess a truncated Painlevé expansion as follows:

$\begin{array}{l}u=u_0+\frac{u_1}{f}+\frac{u_2}{f^2},\\ v=v_0+\frac{v_1}{f}+\frac{v_2}{f^2},\end{array}$ (5)

with $u_0\mathrm{,}{u}_1\mathrm{,}{v}_0\mathrm{,}{v}_1\mathrm{,}{v}_2\mathrm{,}f$ being the function of $x\mathrm{,}y$ and t. We have

$\begin{array}{l}{v}_1=f_{xy},v_2=-2f_x{f}_y,u_2=-2f_x^2,\mathrm{<mi>\; </mi>}{u}_0=\lambda -\frac{f_{xxx}}{2f_x}+\frac{f_{xx}^2}{4f_x^2},\mathrm{<mi>\; </mi>}{u}_1=2f_{xx},\\ v_0=\frac{-3\beta \lambda }{\alpha }+\frac{2f_t}{\alpha f_x}-2\lambda \frac{f_y}{f_x}+\frac{3\beta }{4\alpha }\frac{f_{xx}^2}{f_x^2}+\frac{f_{xx}{f}_{xy}}{f_x^2}-\frac{f_{xxy}}{f_x}-\frac{\beta }{2\alpha }\frac{f_{xxx}}{f_x},\end{array}$ (6)

the Equation (6) is from Ref. [35]. The function f satisfies the following trilinear equation,

$\begin{array}{l}\alpha f_{xxxy}{f}_x^2+4\alpha \lambda f_x^2{f}_{xy}+\beta f_{xxxx}{f}_x^2-\alpha f{f}_x{f}_{xxx}{f}_{xy}-4\beta f_{xxx}{f}_x{f}_{xx}-3\alpha f_x{f}_{xx}{f}_{xy}\\ -4\lambda \alpha f_{xx}{f}_x{f}_y+3\alpha f_{xy}{f}_{xx}^2+3\beta f_{xx}^3-4f_x^2{f}_{xt}+4\alpha 4f_x^2{f}_x{f}_x=0.\end{array}$ (7)

By solving the Equation (5), we can get

$u=\lambda -\frac{1}{2}\frac{f_{xxx}}{f_x}+\frac{f_{xx}^2}{4f_x^2},$ (8a)

$v=\frac{-3\beta \lambda }{\alpha }+\frac{2f_t}{\alpha f_x}-2\lambda \frac{f_y}{f_x}+\frac{3\beta }{4\alpha }\frac{f_{xx}^2}{f_x^2}+\frac{f_{xx}{f}_{xy}}{f_x^2}-\frac{f_{xxy}}{f_x}-\frac{\beta }{2\alpha }\frac{f_{xxx}}{f_x}.$ (8b)

Based on the idea in Ref. [9], the (2 + 1)-dimensional KdV-Calogero-Bogoy-avlenkskii-Schiff equation has a set of rational solutions determined by

$\begin{array}{l}f={\left(a_{1}x+a_{2}y+a_{3}t+a_4\right)}^2+{\left(a_{5}x+a_{6}y+a_{7}t+a_8\right)}^2+a_9\\ +m\exp \left(k_{1}x+k_{2}y+k_{3}t+k_4\right)+n\exp \left(-\left(k_{1}x+k_{2}y+k_{3}t+k_4\right)\right)\mathrm{,}\end{array}$ (9)

where $a_i\mathrm{,}k\mathrm{,}{k}_j\mathrm{,1}\le i\le \mathrm{9,1}\le j\le 4$ are real parameters to be determined. Substituting function (9) into trilinear Equation (7), we give five kinds of exact solutions of Equation (4a) and Equation (4b) via function (9).

3. Lumps and Interaction Solutions

Case I

Substitution function (9) into trilinear Equation (7), a set of constraining equations for the parameters can be obtained as following

$\begin{array}{l}{a}_1=0,a_2=a_2,a_3=a_2\lambda \alpha ,a_4=a_4,\\ a_5=a_5,a_6=-\frac{a_5\beta }{\lambda },a_7=a_7,a_8=a_8,a_9=a_9,\\ m=0,n=0,k_1=k_1,k_2=k_2,k_3=k_3,k_4=k_4,\beta =\beta ,\lambda =\lambda \end{array}$ (10)

where

$\lambda \ne 0,a_9>0$ (11)

to guarantee that the corresponding f is positive. By substituting Equation (10) into function (9), then we get the function f reads

$f={\left(a_{2}y+a_2\lambda \alpha t+a_4\right)}^2+{\left(a_{5}x-\frac{a_5\beta }{\lambda }y+a_{7}t+a_8\right)}^2+a_9\mathrm{.}$ (12)

By means of Equation (8a) and Equation (8b), the solutions of (4a) and (4b) can be obtained

$u=\lambda +\frac{a_5^2}{2A_1^2},$ (13a)

$v=-\frac{3\beta \lambda }{\alpha }+\frac{2B_1}{A_1}-\frac{2\lambda C_1}{A_1}-\frac{\beta a_5^2}{\alpha A_1^2},$ (13b)

where

$\begin{array}{l}f=g^2+h^2+a_9,\\ g=a_{2}y+a_2\lambda \alpha t+a_4,\\ h=a_{5}x-\frac{a_5\beta }{\lambda }y+a_{7}t+a_8,\\ A_1=2a_{5}f,B_1=2a_2\alpha \lambda f+2a_{7}g,C_1=2a_{2}f-2\frac{a_5\beta }{\lambda }g.\end{array}$ (14)

The solutions of u via Equation (13a) which can be seen in Figure 1(a); Figure 1(b) is the corresponding density plot. The solutions of v via Equation (13b) can be seen in Figure 1(c); Figure 1(d) is the corresponding density plot.

Case II

$\begin{array}{l}{a}_1=a_1,a_2=-\frac{a_1\beta }{\alpha },a_3=a_3,a_4=a_4,\\ a_5=0,a_6=a_6,a_7=a_6\alpha \lambda ,a_8=a_8,a_9=a_9,\\ m=m,k_1=0,k_2=k_2,k_3=k_2\alpha \lambda ,k_4=k_4,\end{array}$ (15)

which should satisfy the constraint conditions

$\alpha \ne 0,a_9>0,m>0,$ (16)

to ensure that the corresponding f is positive and well defined. By substituting Equation (15) into function (9), then the function f reads

$\begin{array}{l}f={\left(a_{1}x-\frac{a_1\beta }{\alpha }y+a_{3}t+a_4\right)}^2+{\left(a_{6}y+a_6\alpha \lambda t+a_8\right)}^2\\ +m\exp \left(k_{2}y+k_2\alpha \lambda t+k_4\right)+a_9.\end{array}$ (17)

Using Equation (8a) and Equation (8b), the solutions for (4a) and (4b) can be obtained

$u=\lambda +\frac{a_1^2}{4g^2},$ (18a)

$v=-\frac{3\beta \lambda }{\alpha }+\frac{A_2\alpha \lambda B_2}{\alpha a_{1}g}+\frac{3\beta a_1^2}{16\lambda g^2}-\frac{a_1^3\beta }{\alpha g^2},$ (18b)

where

$\begin{array}{l}f=g^2+h^2+m\exp \left(\xi \right)+a_9,\\ g=a_{1}x-\frac{a_1\beta }{\alpha }y+a_{3}t+a_4,\\ h=a_{6}y+a_6\alpha \lambda t+a_8,\\ A_2=2a_{3}g+2a_6\alpha \lambda h+k_{2}m\alpha \lambda \exp \left(\xi \right),\\ B_2=-2\frac{a_1\beta g}{\alpha }+2a_{6}h+k_{2}m\exp \left(\xi \right),\\ \xi =k_{2}y+k_2\alpha \lambda t+k_4,\end{array}$ (19)

Figure 1. Profiles of the solution u via Equation (13a) and profiles of the solution v via Equation (13b) with $a_2=1$, $a_4=1$, $a_5=1$, $a_7=1$, $a_8=0$, $a_9=0.5$, $\beta =1$, $\alpha =1$, $\lambda =2$, $t=0$. (a) and (c) are the corresponding 3-dimensional plots; (b) and (d) are the corresponding density plots.

The solutions of u via Equation (18a) can be seen in Figure 2(a); Figure 2(b) is the corresponding density plot. The solutions of v via Equation (18b) can be seen in Figure 2(c); Figure 2(d) is the corresponding density plot.

Case III

$\begin{array}{l}{a}_1=a_1,a_2=a_2,a_3=a_3,a_4=a_4,a_5=0,\\ a_6=a_6,a_7=a_6\alpha \lambda ,a_8=a_8,a_9=a_9,\alpha =\alpha ,\beta =\frac{a_2\alpha }{a_1},\\ k_1=0,k_2=k_2,k_3=k_2\alpha \lambda ,k_4=k_4,\lambda =\lambda ,m=m,n=n,\end{array}$ (20)

where

$a_1\ne 0,a_9>0,m>0,n>0.$ (21)

to ensure that the corresponding f is positive and well defined. By substituting Equation (20) into function (9),

$\begin{array}{c}f={\left(a_{1}x+a_{2}y+a_{3}t+a_4\right)}^2+{\left(a_{6}x+a_6\alpha \lambda t+a_8\right)}^2+a_9\\ +m\exp \left(k_{2}y+k_2\alpha \lambda t+k_4\right)+n\exp \left(-\left(k_{2}y+k_2\alpha \lambda t+k_4\right)\right)\mathrm{,}\end{array}$ (22)

where $a_i\mathrm{,}{k}_j\mathrm{,1}\le i\le \mathrm{9,1}\le j\le 4$ are all real parameters to be determined. Then we get the solutions of the Equation (4a) and Equation (4b)

Figure 2. Profiles of the solution u via Equation (18a) and the solution v via Equation (18b) with $a_1=0.5$, $a_3=2$, $a_4=1$, $a_5=0$, $a_6=0$, $a_9=0.5$, $a_8=1$, $\beta =0$, $\alpha =1$, $\lambda =1$, $m=1$, $k_1=0$, $k_2=1$, $k_4=1$, $t=1$. (a) and (c) are the corresponding 3-dimensional plots; (b) and (d) are the corresponding density plots.

$u=\lambda +\frac{a_1^2}{4g^2},$ (23a)

$v=-\frac{3a_2\lambda }{a_1}+2\frac{B_3-\lambda \alpha C_3}{\alpha A_3}+\frac{a_1^3{a}_2}{A_3^2},$ (23b)

with

$\begin{array}{l}f=h^2+g^2+a_9+m\exp \left(\xi \right)+n\exp \left(-\xi \right),\\ g=a_{1}x+a_{2}y+a_{3}t+a_4,\\ h=a_{6}x+a_6\alpha \lambda t+a_8,\\ \xi =k_{2}y+k_2\alpha \lambda t+k_4,\\ A_3=2a_{1}g,C_3=2a_{2}g+2a_{6}h+m{k}_2\exp \left(\xi \right)-n{k}_2\exp \left(-\xi \right)\\ B_3=2a_{3}g+2a_6\alpha \lambda h+m{k}_2\alpha \lambda \exp \left(\xi \right)-n{k}_2\alpha \lambda \exp \left(\xi \right)\end{array}$ (24)

The exact solutions for Equation (23a) and Equation (23b) can be seen in Figure 3(a) and Figure 3(c), Figure 3(b) and Figure 3(d) are the corresponding density plots.

Case IV

Figure 3. Profiles of the solution u via Equation (23a) and v via Equation (23b) with $a_1=1$, $a_2=1$, $a_3=2$, $a_4=2$, $a_5=0$, $a_6=1$, $a_8=1$, $a_9=1$, $k_1=0$, $k_2=1$, $k_4=1$, $m=1$, $n=1$, $\alpha =1$, $\lambda =1$, $t=0$. (a) and (c) are the corresponding 3-dimensional plots; (b) and (d) are the corresponding density plots.

$\begin{array}{l}{a}_1=a_1,a_2=-\frac{a_1\beta }{\alpha },a_3=-\frac{a_1{a}_7}{\alpha a_6},a_4=-\frac{a_1{a}_8}{a_6\alpha },a_5=-\frac{\alpha a_6}{\beta },\\ a_6=a_6,a_7=a7,a_8=a_8,a_9=a_9,m=m,n=n,\\ \alpha =\alpha ,\beta =\beta ,k_1=0,k_2=k_2,k_3=k_2\alpha \lambda ,k_4=k_4,\lambda =\lambda \end{array}$ (25)

where

$\alpha \beta a_6\ne 0,a_9>0,m>0,n>0.$ (26)

to guarantee that the corresponding f is positive and well defined. We substitute Equation (25) into Equation (9), hence, reinstall function f as the following formula:

$\begin{array}{l}f={\left(a_{1}x-\frac{a_1\beta }{\alpha }y-\frac{a_1{a}_7}{\alpha a_6}t-\frac{a_1{a}_8}{a_6\alpha }\right)}^2+{\left(a_{5}x+a_{6}y+a_{7}t+a_8\right)}^2+a_9\\ +m\exp \left(k_{2}y+k_2\alpha \lambda t+k_4\right)+n\exp \left(-\left(k_{2}y+k_2\alpha \lambda t+k_4\right)\right)\mathrm{,}\end{array}$ (27)

where $a_i\mathrm{,}{k}_j\mathrm{,1}\le i\le \mathrm{9,1}\le j\le 4$ are all real parameters to be determined. Then, we get

$u=\lambda +\frac{D_4^2}{A_4^2},$ (28a)

$v=-\frac{\beta \lambda }{\alpha }+\frac{2C_4}{\alpha A_4}+\frac{2\lambda B_4}{A_4}+\frac{D_4\left(3\beta D_4+4\alpha E_4\right)}{4\alpha A_4^2},$ (28b)

where

$\begin{array}{l}f=h^2+g^2+\alpha \exp \left(\xi \right)+\beta \exp \left(-\xi \right)+a_9,\\ h=a_{1}x-\frac{a_1\beta }{\alpha }y-\frac{a_1{a}_7}{\alpha a_6}t-\frac{a_1{a}_8}{a_6\alpha },\\ g=a_{5}x+a_{6}y+a_{7}t+a_8,\\ \xi =k_{2}y+k_2\alpha \lambda t+k_4,\end{array}$

$\begin{array}{l}{A}_4=2\left(a_{1}h+a_{5}g\right),\\ B_4=-2\frac{a_1\beta }{\alpha }h+2a_{6}g+k_2\left(m\exp \left(\xi \right)-n\exp \left(-\xi \right)\right),\\ C_4=-2\frac{a_1{a}_7}{\alpha a_6}h+2a_{7}g+k_2\alpha \lambda \left(m\exp \left(\xi \right)-n\exp \left(-\xi \right)\right),\\ D_4=2\left(a_1^2+a_5^2\right),\\ E_4=2\left(a_5{a}_6-\frac{a_1^2\beta }{\alpha }\right).\end{array}$ (29)

The exact solution for Equation (28b) can be seen in Figure 4(a), Figure 4(b) Figure 4(c) the corresponding density plots with the different time.

Figure 4. Profiles of the solution v via Equation (28b) with $a_1=1$, $a_4=0$, $a_5=2$, $a_6=4$, $a_8=0$, $a_9=1$, $\alpha =1$, $\beta =1$, $k_1=1$, $k_2=1.5$, $k_3=2$, $k_4=0$. (a) 3-dimensional plot with the time t = 0; (b) (c) (b) the corresponding density plot with the different time.

Case V.

$\begin{array}{l}{a}_1=0,a_2=a_2,a_3=\alpha \lambda a_2,a_4=a_4,a_5=\frac{a_6{k}_1}{k_2}\\ a_6=a_6,a_7=a_7,a_8=a_8,a_9=a_9,\alpha =-\frac{\beta k_1}{k_2},\beta =\beta ,\\ k_1=k_1,k_2=k_1,k_3=\frac{a_7{k}_2}{a_6},k_4=k_4,m=m,n=n.\end{array}$ (30)

with

$k_2{a}_6\ne 0,a_9>0,m>0,n>0.$ (31)

to guarantee that the corresponding f is positive and well defined. We substitute Equation (30) into Equation (9), then function f reads

$\begin{array}{c}f={\left(a_{2}y+\alpha \lambda a_{2}t+a_4\right)}^2+{\left(\frac{a_6{k}_1}{k_2}x+a_{6}y+a_{7}t+a_8\right)}^2\\ +a_9+m\exp \left(k_{1}x+k_{2}y+\frac{a_7{k}_2}{a_6}t+k_4\right)\\ +n\exp \left(-\left(k_{1}x+k_{2}y+\frac{a_7{k}_2}{a_6}t+k_4\right)\right)\mathrm{,}\end{array}$ (32)

where $a_i\mathrm{,}{k}_j\mathrm{,1}\le i\le \mathrm{9,1}\le j\le 4$ are all real parameters to be determined. Then, the rational solution of system (4a) and (4a) can be got again.

$\begin{array}{l}u=\lambda -\frac{k_1^3\left(m\exp \left(\xi \right)-n\exp \left(-\xi \right)\right)}{2A_5}+\frac{D_5}{4A_5^2},\\ v=-\frac{3\beta \lambda }{\alpha }+2\frac{B_5}{A_5}+2\lambda \frac{C_5}{A_5},+\frac{3\beta D_5^2}{4\alpha A_5^2}+\frac{D_5{E}_5}{A_5^2}-\frac{F_5}{A_5}-\frac{\beta G_5}{2\alpha A_5},\end{array}$ (33)

and where

$\begin{array}{l}f=g^2+h^2+m\exp \left(\xi \right)+n\exp \left(-\xi \right)+a_9,\\ g=a_{2}y+\alpha \lambda a_{2}t+a_4,\\ h=\frac{a_6{k}_1}{k_2}x+a_{6}y+a_{7}t+a_8,\\ \xi =k_{1}x+k_{2}y+\frac{a_7{k}_2}{a_6}t+k_4,\\ A_5=2\frac{a_6{k}_1}{k_2}+k_{1}h+m{k}_1\exp (\xi )-n{k}_1\exp \left(-\xi \right),\end{array}$

$\begin{array}{l}{B}_5=2\alpha \lambda a_{2}g+2a_{7}h+\frac{a_7{k}_2}{a_6}\left(m\exp \left(\xi \right)+n\exp \left(-\xi \right)\right),\\ C_5=2a_{2}g+2a_{6}h+k_2\left(m\exp \left(\xi \right)+n\exp \left(-\xi \right)\right),\\ D_5=2\frac{a_6^2{k}_1^2}{k_2^2}+k_1^2\left(m\exp \left(\xi \right)+n\exp \left(-\xi \right)\right),\\ E_5=k_1{k}_2\left(m\exp \left(\xi \right)+n\exp \left(-\xi \right)\right)+2a_6\frac{a_6{k}_1}{k_2},\\ F_5=k_1^2{k}_2\left(m\exp \left(\xi \right)+n\exp \left(-\xi \right)\right).\end{array}$ (34)

Figure 5. Profiles of the solution u and v via Equation (33) with $a_2=1$, $a_4=1$, $a_6=1$, $a_7=1$, $a_8=2$, $a_9=1$, $k_1=1$, $k_2=1$, $k_4=1$, $m=1$, $n=1$, $\lambda =1$, $\beta =1$, $t=0$. (a) and (c) are the 3-dimensional plots; (b) and (d) are the corresponding density plots.

The exact solutions for Equation (33) can be seen in Figure 5(a) and Figure 5(c), Figure 5(b) and Figure 5(d) are the corresponding density plots.

4. Conclusion

In this paper, by using the truncated Painlev expansion, the (2 + 1)-dimensional KdV-CBS equation can be changed to a trilinear Equation (7), and by means of symbolic computations, we presented five types of exact solutions ((13a), (13b), (18a), (18b), (23a), (23b), (28a), (28b), (33)) to the (2 + 1)-dimensional KdV-CBS equation. Three-dimensional plots and the corresponding density plots of these exact solutions are given respectively in Figure 1 and Figure 5. This work shows the power of the methods, the exact solutions of some nonlinear equations can be obtained by using this means. It is expected that our results can enrich the exact solutions of the nonlinear equations.

Acknowledgements

This work is supported by the National Natural Science Foundation of China under grant number 11571008.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Gear, J.A. and Grimshaw, R. (1984) Weak and Strong Interactions between Internal Solitary Waves. Studies in Applied Mathematics, 70, 235-258. https://doi.org/10.1002/sapm1984703235
[2]	Ruggieri, M. and Speciale, M.P. (2013) Similarity Reduction and Closed Form Solutions for a Model Derived from Two-Layer Fluids. Advances in Difference Equations, 1, 355-363. https://doi.org/10.1186/1687-1847-2013-355
[3]	Christodoulides, R.G.P. (2010) Steady Gap Solitons in a Coupled Korteweg-de Vries System: A Dynamical Systems Approach. Physica D: Nonlinear Phenomena, 239, 635-639. https://doi.org/10.1016/j.physd.2010.01.016
[4]	Weoss, J., Tabor, M. and Carnevale, G. (1983) Painlevé Property for Partial Differential Equations. Journal of Mathematical Physics, 24, 522-526. https://doi.org/10.1063/1.525721
[5]	Matveev, V.B. and Salle, M.A. (1991) Darboux Transformmations and Solitons. Springer, Berlin. https://doi.org/10.1007/978-3-662-00922-2
[6]	Bluman, G.W., Checiakov, A.F. and Anco, S.C. (2010) Applications of Symmetry Methods to Partial Differential Equations. Springer, New York. https://doi.org/10.1007/978-0-387-68028-6
[7]	Hirota, R. (2004) The Direct Method in Soliton Theory. Spring, Berlin. https://doi.org/10.1017/CBO9780511543043
[8]	Ablowitz, M.J. and Clarkkson, P.A. (1991) Solitons, Nonlinear Evolution Equations and Inverse Scattering. Cambridge University Press, Cambrige. https://doi.org/10.1017/CBO9780511623998
[9]	Ma, W.X. (2015) Lump Solutions to the Kadomtsev-Petviashvili Equation. Physics Letters A, 379, 1975-1978. https://doi.org/10.1016/j.physleta.2015.06.061
[10]	Ma, W.X., Zhou, Y. and Dougherty, R. (2016) Lump-Type Solutions to Nonlinear Differential Equantions Derived from Generalized Bilinear Equations. International Journal of Modern Physics B, 30, Artical ID: 1640018. https://doi.org/10.1142/S021797921640018X
[11]	Zhang, H.Q. and Ma, W.X. (2016) Lump Solution to the (2 + 1)-Dimensional Sawada-Kotera Equation. Nonlinear Dynamics, 87, 2305-2310. https://doi.org/10.1007/s11071-016-3190-6
[12]	Hua, Y.F., Guo, B.L. and Ma, W.X. (2016) Interaction Behavior Associated with a Generalized (2 + 1)-Dimensional Hirota Bilinear Equation for Nonlinear Waves. Applied Mathematical Modelling, 74, 184-198. https://doi.org/10.1016/j.apm.2019.04.044
[13]	Zhang, H.Q. and Ma, W.X. (2017) Mixed Lump-Kink Solutions to the KP Equation. Computers and Mathematics with Applications, 74, 1399-1405. https://doi.org/10.1016/j.camwa.2017.06.034
[14]	Yang, J.Y. and Ma, W.X. (2017) Abundant Lump-Type Solutions of the Jimbo Miwa Equation in (2 + 1)-Dimensions. Computers and Mathematics with Applications, 73, 220-225. https://doi.org/10.1016/j.camwa.2016.11.007
[15]	Wang, Y.H., Wang, H., Zhang, H.S. and ChaoLu, T.M. (2017) Exact Interaction Solutions of an Extended (2 + 1)-Dimensinal Shallow Water Wave Equation. Commnications in Theoretical Physics, 68, 165. https://doi.org/10.1088/0253-6102/68/2/165
[16]	Yang, J.Y., Ma, W.X. and Qin, Z.Y. (2017) Lump and Lump-Soliton Solutions to the (2 + 1) Dimensional Ito Equation. Analysis and Mathematical Physics, 8, 427-436.
[17]	Wang, H., Wang, Y.H. and Ma, W.X. (2018) Lump Solutions of a New Extended (2 + 1)-Dimensional Boussinesq Equation. Modern Physics Letters B, 32, Article ID: 1850376. https://doi.org/10.1142/S0217984918503761
[18]	Ma, W.X. (2018) Abundant Lumps and Their Interaction Solutions of (2 + 1)-Dimensional Linear PDEs. Journal of Geometry and Physics, 133, 10-16. https://doi.org/10.1016/j.geomphys.2018.07.003
[19]	Sumayah, B. and Ma, W.X. (2018) A Study of Lump-Type and Interaction Solutions to a (2 + 1)-Dimensional Jimbo-Miwa-Like Equation. Computers and Mathematics with Applications, 76, 1576-1582. https://doi.org/10.1016/j.camwa.2018.07.008
[20]	Ma, W.X. and Zhou, Y. (2018) Lump Solutions to Nonlinear Partial Differential Equations via Hirota Bilinear Forms. Journal of Differential Equations, 264, 2633-2659. https://doi.org/10.1016/j.jde.2017.10.033
[21]	Ma, W.X., Yong, X. and Zhang, H.Q. (2018) Diversity of Interaction Solutions to the (2 + 1)-Dimensinal Ito Equation. Computers and Mathematics with Applications, 75, 289-295. https://doi.org/10.1016/j.camwa.2017.09.013
[22]	Huang, L.L., Yue, Y.F. and Chen, Y. (2018) Localized Waves and Interaction Solutions to a (2 + 1)-Dimensional Generalized KP Equation. Computers and Mathematics with Applications, 76, 831-844. https://doi.org/10.1016/j.camwa.2018.05.023
[23]	Ma, W.X. (2019) Lump and Interaction Solutions to Linear (2+1)-Dimensional PDEs. Acta Mathematica Scientia, 39, 498-508.
[24]	Ma, W.X. (2019) Lump and Interaction Solutions of Linear PDEs in (2 + 1)-Dimensions. Aast Asian Journal on Applied Mathematics, 9, 185-194. https://doi.org/10.4208/eajam.100218.300318
[25]	Ren, B., Ma, W.X. and Yu, J. (2019) Characteristics and Interactions of Solitary and Lump Waves of a (2 + 1)-Dimensional Coupled Nonlinear Partial Differential Equation. Nonlinear Dynamics, 96, 717-727. https://doi.org/10.1007/s11071-019-04816-x
[26]	Zhou, Y., Manukure, S. and Ma,W.X. (2019) Lump and Lump-Soliton Solutions to the Hirota-Satsuma-Ito Equation. Communications in Nonlinear Science and Numerical Simulation, 68, 56-62. https://doi.org/10.1016/j.cnsns.2018.07.038
[27]	Wang, Y.H. and Chen, Y. (2013) Binary Bell Polynomial Manipulations on the Integrability of a Generalized (2 + 1)-Dimensional Korteweg-de Vries Equation. Journal of Mathematical Analysis and Applications, 400, 624-634. https://doi.org/10.1016/j.jmaa.2012.11.028
[28]	Peng, Y.Z. (2010) A New (2 + 1)-Dimensional KdV Equation and Its Localized Structures. Communications in Theoretical Physics, 54, 863-865. https://doi.org/10.1088/0253-6102/54/5/17
[29]	Chen, J.C. and Ma, Z.Y. (2017) Consistent Riccati Expansion Solvability and Soliton-Cnoidal Wave Interaction Solution of a (2 + 1)-Dimensional Korteweg-de Vries Equation. Applied Mathematics Letters, 64, 87-93. https://doi.org/10.1016/j.aml.2016.08.016
[30]	Toda, K. and Yu, S.J. (2000) The Investigation into the Schwarz-Korteweg-de Vries Equation and the Schwarz Derivative in (2 + 1)-Dimensional. Journal of Mathematical Physics, 41, 4747-4751. https://doi.org/10.1063/1.533374
[31]	Zhang, Y. and Xu, G.Q. (2014) Intergrability and Exact Solutions for a (2 + 1)-Dimensional Variable-Coefficient KdV Equation. Applications of Mathematics, 9, 646-658.
[32]	Lü, X., Lin, F.H. and Qi, F.H. (2015) Analytical Study on a Two Dimensional Korteweg-de Vries Model with Bilinear Representation, Bäcklund Transformaation and Soliton Solutions. Applications of Mathematical Models, 39, 3221-3226. https://doi.org/10.1016/j.apm.2014.10.046
[33]	Lü, X., Ma, W.X. and Khalique, C.M. (2015) A Direct Bilinear Bäcklund Transformaation of a (2 + 1)-Dimensional Korteweg-de Vries Equation. Applied Mathematics Letters, 50, 37-42. https://doi.org/10.1016/j.aml.2015.06.003
[34]	Zhao, Z.L. and Han, B. (2017) The Riemann-Bäcklund Method to a Quasiperiodic Wave Solvable Generalized Variable-Coefficient (2 + 1)-Dimensional KdV Equation. Nonlinear Dynamics, 87, 2661-2676. https://doi.org/10.1007/s11071-016-3219-x
[35]	Huang, L.L. and Chen, Y. (2018) Nonlocal Symmetry and Similarity Reductions for a (2 + 1)-Dimensional Korteweg-de Vries Equation. Nonlinear Dynamics, 92, 221-234. https://doi.org/10.1007/s11071-018-4051-2