Integrability Tests and Some New Soliton Solutions of an Extended Potential Boiti-Leon-Manna-Pempinelli Equation

Miao Li, Wei Tan^*, Houping Dai

College of Mathematics and Statistics, Jishou University, Jishou, China.

DOI: 10.4236/jamp.2022.1010194                

Abstract

This paper is devoted to the study of a (2 + 1)-dimensional extended Potential Boiti-Leon-Manna-Pempinelli equation. Firstly, By means of the standard Weiss Tabor Carnevale approach and Kruskal’s simplification, we prove the painlevé non integrability of the equation. Secondly, A new breather solution and lump type solution are obtained based on the parameter limit method and Hirota’s bilinear method. Besides, some interaction behavior between lump type solution and N-soliton solutions (N is any positive integer) are studied. We construct the existence theorem of the interaction solution and give the process of calculation and proof. We also give a concrete example to illustrate the effectiveness of the theorem, and some spatial structure figures are displayed to reflect the evolutionary behavior of the interaction solutions with the change of soliton number N and time t.

Keywords

BLMP Equation, Lump Type Solution, Interaction Behavior, Parameter Limit Method, Hirota’s Bilinear Method

1. Introduction

Many nonlinear phenomena in nature can be represented by nonlinear partial differential equations (NLPDE). Soliton theory, as one of the important contents in the field of NLPDE, has received extensive attention in recent years. Therefore, it is of great significance to find soliton solutions of nonlinear equations and analyze their structural properties. With the development of soliton theory, many methods of solving soliton are proposed, such as the Darboux transformation method [1] [2], Hirota bilinear method [3] [4], polynomial-expansion method [5] [6], Lie group method [7], wronskian technique [8], etc. A large of exact solutions of nonlinear equations, including soliton solutions, periodic solutions, rational solutions and interaction solutions and so on, have been studied by these methods. In recent years, the research on lump solutions and interaction behavior between solitons, such as rogue wave, the fusion and degeneration of the lump solution, especially the interaction between lump solutions and solitons [9] [10] [11], has become a hot topic because of its application value in nonlinear science.

In this paper, we consider the following (2 + 1)-dimensional extended Potential Boiti-Leon-Manna-Pempinelli (ePBLMP) equation [12],

$u_{yt}+c_1{u}_{xxxy}-3c_1\left(u_y{u}_{xx}+u_x{u}_{xy}\right)-c_2{u}_{xy}-c_3{u}_{yy}=\mathrm{0,}$ (1)

where u is a multivariate function about $x\mathrm{,}y\mathrm{,}t$, and $c_1\mathrm{,}{c}_2\mathrm{,}{c}_3$ are arbitrary real parameters. In fact, when $c_1=1$, $c_2=c_3=0$, Equation (1) will be simplified as the following standard form of the pBLMP equation.

$u_{yt}+u_{xxxy}-3u_y{u}_{xx}-3u_x{u}_{xy}=0.$ (2)

The (2 + 1)-dimensional BLMP equation [13] was firstly proposed by Boiti et al. and based on the relationship of weak Lax pairs. Equation (2) is described the interaction between the long wave propagating along the x-axis and the Riemann wave propagating along the y-axis. Some exact solutions of the BLMP equation were obtained in works [14] [15] [16] [17]. Recently, some research results have been made on eBLMP, such as Ren studying the D’Alembert wave and soliton molecules of the eBLMP equation [12]. Paliathanasis obtained some new periodic solutions of the eBLMP equation [18]. In addition, there are many methods to study similar equations. For example, Luo obtained the parametric form of soliton solutions according to the solution of the constructed Riemann-Hilbert problem [19], and Zeng obtained a class of soliton solutions of the (2 + 1)-dimensional time fractional equation by using the new mapping method [20].

In this paper, we first give detailed proof of the integrability of Equation (1). Moreover, as far as we know, the breather solutions, lump solutions and some mixed solutions of Equation (1) have not been obtained, and we will get these new results. More importantly, we find that most of the papers only study the interaction solutions between the first-order lump solutions and N-solitons, where the number of solitons N is usually less than 3. We obtain the existence theorem of the interaction solutions between the higher-order lump solutions of Equation (1) and different types of N-soliton solutions. This paper will study these new nonlinear phenomena and dynamic behaviors.

2. Painlevé Test

In this section, we give out the standard test with Kruskal’s simplification, which is constituted by three steps, the leading order analysis, the resonance determination and the resonance conditions’ verification, respectively. According to the standard WTC approach, we investigated the general solution of Equation (1) of the following form:

$\{\begin{array}{l}u\left(x,y,t\right)={\phi }^{-\alpha }{\displaystyle \underset{j=0}{\overset{N}{\sum }}}{u}_j\left(x,y,t\right){\phi }^j,\\ \phi \left(x,y,t\right)=x+h\left(y,t\right).\end{array}$ (3)

Step 1: Leading item analysis.

We assume the first item of (3) is

$u\left(x\mathrm{,}y\mathrm{,}t\right)={\phi }^{-\alpha }{u}_0\mathrm{.}$ (4)

Substituting (4) into (1), comparing the lowest derivative term and the nonlinear term, we get

$\alpha =1,\mathrm{<mtext>\; \hspace{0.17em}\; </mtext>\; <mtext>\; \hspace{0.17em}\; </mtext>\; <mtext>\; \hspace{0.17em}\; </mtext>\; <mtext>\; \hspace{0.17em}\; </mtext>}{u}_0\left(x,y,t\right)=-2{\phi }_x.$ (5)

Step 2: Determine resonance point and recurrence relation.

Substituting $\alpha =1$ into (3), we obtain the expansion of u, that is

$u\left(x,y,t\right)={\displaystyle \underset{j=0}{\overset{N}{\sum }}}{u}_j{\phi }^{j-1}.$ (6)

Inserting (6) and its derivatives into (1). Comparing the power terms of $\phi$ and letting the coefficient be zero, we can get the recurrence relation about $u_j$,

$\begin{array}{l}{u}_{j-3,t}\left(j-4\right){\phi }_y+u_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y+u_{j-3}\left(j-4\right){\phi }_{yt}\\ -c_2{u}_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y-c_3\left(u_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y+u_{j-3}\left(j-4\right){\phi }_{yy}\right)\\ +c_1{u}_j\left(j-1\right)\left(j-2\right)\left(j-3\right)\left(j-4\right){\phi }_y\\ -6c_1{\displaystyle \underset{i=0}{\overset{j}{\sum }}}{u}_i\left(i-1\right)u_{j-i}\left(j-i-1\right)\left(j-i-2\right){\phi }_y=0.\end{array}$ (7)

Finishing items that contain $u_j$,

$\begin{array}{l}{c}_1{u}_j{\phi }_y\left(j-1\right)\left(j+1\right)\left(j-4\right)\left(j-6\right)\\ =-u_{j-3,t}\left(j-4\right){\phi }_y-u_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y-u_{j-3}\left(j-4\right){\phi }_{yt}\\ +c_2{u}_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y+c_3\left(u_{j-2}\left(j-3\right)\left(j-4\right){\phi }_y+u_{j-3}\left(j-4\right){\phi }_{yy}\right)\\ +6c_1{\displaystyle \underset{i=1}{\overset{j-1}{\sum }}}{u}_i\left(i-1\right)u_{j-i}\left(j-i-1\right)\left(j-i-2\right){\phi }_y.\end{array}$ (8)

It can be found from (8) that the resonance point occurs when

$j=-1,1,4,6.$ (9)

Step 3: Verify compatibility conditions.

In (7), when $j=0$,

$24c_1{u}_0+12c_1{u}_0^2=0.$ (10)

Solve (10) and get $u_0=-2$.

When $j=1$, we have

$-6c_1\left(u_0{u}_1\cdot 0+u_1{u}_0\cdot 0\right)=0,$ (11)

We get that $u_1$ is an arbitrary function, so $j=1$ is the resonance point.

When $j=2$, we get

$u_2=\frac{1-c_2-c_3}{6c_1}.$ (12)

When $j=3$, we get

$u_3=\frac{c_1{\phi }_{yy}-{\phi }_{yt}}{48c_1{\phi }_y}.$ (13)

When $j=4$, we have

$-6c_1\left(-6u_0{u}_4+6u_4{u}_0\right)=0.$ (14)

We can get that $u_4$ is an arbitrary function, so $j=4$ is the resonance point.

Substituting $j=6$ into the left side of (7), we get

$\begin{array}{l}2u_{\mathrm{3,}t}{\phi }_y+6u_4{\phi }_y+2u_3{\phi }_{yt}-6c_2{u}_4{\phi }_y-6c_3\left(u_4{\phi }_y+2u_3{\phi }_{yy}\right)\\ +\text{120}{c}_1{u}_6{\phi }_y-6c_1\left(40u_6{\phi }_y+6u_2{u}_4{\phi }_y+4u_3^2{\phi }_y+10u_0{u}_6\right)\mathrm{.}\end{array}$ (15)

The items containing $u_3$, $u_4$, $u_6$ in (16) are combined respectively, and the sum of these items is zero. But the left side of (15) is $2u_{\mathrm{3,}t}{\phi }_y$, which is reduced to $\frac{-{\phi }_{ytt}}{24c_1}$. Therefore, when ${\phi }_{ytt}=0$, the (10) is integrable.

In summary, $j=-1,1,4$ satisfy the resonance condition, but $j=6$ does not satisfy the resonance condition. Although it does not pass the Painlevé test, it has many undiscovered mathematical properties. For example, we can obtain the lump solutions and N-solitons by introducing bilinear method.

3. Degradation of Breather Solution and Emergence of Lump Solution

In this section, we study the degeneration of the double breathing solutions and the emergence of lump solutions of the ePBLMP equation. According to Ref. [12], we can assume that

$u\left(x,y,t\right)=-2{\left(\ln f\left(x,y,t\right)\right)}_x,$ (16)

where $f=f\left(x,y,t\right)$ represents the real function to be selected. Then substituting Equation (16) into Equation (1), we can obtain the following bilinear form of Equation (1),

$\left(D_y{D}_t+c_1{D}_x^3{D}_y-c_2{D}_x{D}_y-c_3{D}_y^2\right)f\cdot f=0,$ (17)

where D denotes the bilinear operators. In Equation (17), through the extended homoclinic test method [21], we choose a new test function, which is first applied to Equation (1).

$f\left(x,y,t\right)=m_1\cosh \left({\eta }_1\right)+m_2\cos \left({\eta }_2\right)+m_3\cosh \left({\eta }_3\right),$ (18)

where ${\eta }_i=k_i\left(a_{i}x+b_{i}y+d_{i}t+r_i\right)$ and $k_i,a_i,b_i,d_i,r_i\left(i=1,2,3\right)$ are arbitrary parameters. By substituting Equation (18) into Equation (17), the following parameters are obtained:

$k_1=-\frac{a_3{k}_3}{a_1},a_2=\frac{a_3{k}_3^2{b}_3\left(m_1^2-m_3^2\right)}{k_2^2{b}_2{m}_2^2},b_1=-\frac{a_1{b}_3}{a_3},$

$\begin{array}{l}{d}_1=a_1\left(-a_3^2{c}_1{k}_3^2+c_2-\frac{b_3{c}_3}{a_3}+3\alpha \right),\\ d_2=\frac{a_3{b}_3{k}_3^2\left(m_1^2-m_3^2\right)}{b_2{k}_2^2{m}_2^2}\left(-3a_3^3{k}_3^2{c}_1-c_2+\alpha \right)+b_2{c}_3,\\ d_3=a_3(-a_3^2{k}_3^2{c}_1+c_2+3\alpha )+b_3{c}_3,\end{array}$ (19)

where $\alpha =\frac{a_3^3{b}_3^2{k}_3^4{c}_1{\left(m_1^2-m_3^2\right)}^2}{b_2^2{k}_2^2{m}_2^4}$. Plugging Equation (19) into Equation (16), we can obtain the double breathing solution of Equation (1), (See Figures 1(a)-(c)).

$u=-\frac{2(m_1{k}_1{a}_1\sinh \left({\eta }_1\right)-m_2{k}_2{a}_2\sin \left({\eta }_2\right)+m_3{k}_3{a}_3\sinh \left({\eta }_3\right)}{m_1\cosh \left({\eta }_1\right)+m_2\cos \left({\eta }_2\right)+m_3\cosh \left({\eta }_3\right)}$ (20)

where

$\{\begin{array}{l}{\eta }_1=-\frac{a_3{k}_3}{a_1}\left(a_{1}x-\frac{a_1{b}_3}{a_3}y+a_1\left(-a_3^2{c}_1{k}_3^2+c_2-\frac{b_3{c}_3}{a_3}+3\alpha \right)t+r_1\right),\\ {\eta }_2=k_2\left(\frac{a_3{k}_3^2{b}_3\left(m_1^2-m_3^2\right)}{k_2^2{b}_2{m}_2^2}x+b_{2}y+\frac{a_3{b}_3{k}_3^2\left(m_1^2-m_3^2\right)}{b_2{k}_2^2{m}_2^2}\left(-3a_3^3{k}_3^2{c}_1-c_2+\alpha \right)t+b_2{c}_{3}t+r_2\right),\\ {\eta }_3=k_3\left(a_{3}x+b_{3}y+a_3\left(-a_3^2{k}_3^2{c}_1+c_2+3\alpha \right)t+b_3{c}_{3}t+r_3\right).\end{array}$ (21)

When $m_i\left(i=1,2,3\right)$ satisfies $\underset{k_i\to 0}{\lim }\left(m_1+m_2+m_3\right)=0$, $k_2=p{k}_3$, $m_1=\beta \cosh \left(\rho k_3\right)$, $m_3=\gamma \cosh \left(\sigma k_3\right)$, $m_2=-\left(\beta +\gamma \right)\cos \left(n{k}_3\right)$, letting $k_3\to 0$ in Equation (20). We will get the lump type solution of Equation (1), (See Figure 1(d)).

$u=-\frac{a_3\left(\gamma {\theta }_3-\beta {\theta }_1\right)-\left(\beta +\gamma \right)p^2\delta {\theta }_2}{\beta {\theta }_1^2+\left(\beta +\gamma \right)p^2{\theta }_2^2+\gamma {\theta }_3^2+\beta {\rho }^2+\left(\beta +\gamma \right)n^2+\gamma {\sigma }^2}$ (22)

where $\delta =\underset{k_3\to 0}{\lim }{a}_2=\frac{a_3{k}_3^2{b}_3\left(\beta -\gamma \right)}{p^2{b}_2\left(\beta +\gamma \right)}$,

$\{\begin{array}{l}{\theta }_1=-a_{3}x+b_{3}y+a_1\left(c_2-\frac{b_3{c}_3}{a_3}\right)t-\frac{a_3}{a_1}{r}_1,\\ {\theta }_2=\delta x+b_{2}y+\delta c_{2}t+b_2{c}_{3}t+r_2,\\ {\theta }_3=a_{3}x+b_{3}y+a_3{c}_{2}t+b_3{c}_{3}t+r_3.\end{array}$ (23)

In fact, the breather solution (20) represented is double breather solitons which take on the characteristics of a periodic wave with variable ${\eta }_2$, and shows the characteristics of double solitary waves with variable ${\eta }_1$ and ${\eta }_3$, the breather solution with this form is also called double breather solution by some scholars [22] [23]. In Figure 1, we give the evolution diagram of spatial structure for the double breather solutions with the change of time t. However, the period $\frac{2\pi }{p{k}_3}$ of $\cos \left({\eta }_2\right)$ tends to infinity as $k_3\to 0$, the breather and period of solution (20) is degenerate as a lump type solution (22) when the parametric $k_3$ tends to 0.

(a) (b) (c) (d)

Figure 1. The evolution behavior of the breather solution Equation (20) with time: (a) $t=-0.1$ ; (b) $t=0$ ; (c) $t=0.1$, and where the parameters are set as $a_1=-\frac{9}{10}$, $a_3=2$, $b_2=-\frac{1}{2}$, $b_3=3$, $m_1=m_2=\frac{3}{2}$, $m_3=\frac{1}{2}$, $k_2=k_3=\frac{1}{4}$, $c_1=c_3=1$, $c_2=3$. (d) Spatial structures of lump type solution Equation (22): $a_1=a_3=b_3=n=p=2$, $m_2=-\frac{3}{2}$, $m_3=\frac{1}{2}$, $k_2=k_3=\frac{1}{4}$, $m_1=b_2=c_1=c_3=\beta =\rho =\gamma =\sigma =1$, $c_2=3$, $r_1=r_2=r_3=0$, $t=0$.

4. Evolutionary Behavior of Interaction

Now, we study the interaction between lump solutions and N-solitons (where N is any positive integer). Firstly, we give a theorem and prove the calculation process. Here, we choose a new test function:

$\{\begin{array}{l}h\left(x,y,t\right)=a_0+{\displaystyle \underset{i=1}{\overset{M}{\sum }}}{\left(a_{i}x+b_{i}y+d_{i}t+m_i\right)}^2,\\ g\left(x,y,t\right)={\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j{\text{e}}^{\left(p_{j}x+q_{j}y+n_{j}t+s_j\right)},\\ f\left(x,y,t\right)=h\left(x,y,t\right)+g\left(x,y,t\right),\end{array}$ (24)

where $a_0,a_i,b_i,d_i,m_i\left(i=1,2,\cdots ,M\right),{\delta }_j,p_j,n_j,s_j\left(j=1,2,\cdots ,N\right)$ are some free real parameters. Obviously, $h\left(x,y,t\right)$ is a polynomial function, and $g\left(x,y,t\right)$ is an exponential function. They combined to form $f\left(x,y,t\right)$. By substituting (24) into (16), we will get a new interaction solution of ePBLMP.

Theorem 1. Let $f\left(x,y,t\right)$ satisfies ${\displaystyle {\sum }_{i=1}^N}{a}_i{b}_i=0$, $d_i=c_2{a}_i+c_3{b}_i$, ${\displaystyle {\sum }_{i=1}^N}{b}_i{d}_i-c_3{\displaystyle {\sum }_{i=1}^N}{b}_i^2=0$, $q_j=0$, $n_j=c_2{a}_j+c_1{p}_j^3$, then $f\left(x,y,t\right)$ is the solution of Equation (1).

Proof. By direct calculation

$\begin{array}{l}P\left(D_x,D_y{D}_t\right)f\cdot f\\ =\left(D_t+c_1{D}_x^3-c_2{D}_x-c_3{D}_y\right)D_y\left(h\cdot h+h\cdot g+g\cdot h+g\cdot g\right).\end{array}$ (25)

We first consider

$\begin{array}{l}\left(D_y{D}_t+c_1{D}_x^3{D}_y-c_2{D}_x{D}_y-c_3{D}_y^2\right)h\cdot h\\ =2h_{yt}h-h_y{h}_t+2c_1\left(h_{yxxx}h-3h_{yxx}{h}_x+3h_{yx}{h}_{xx}-h_y{h}_{xxx}\right)\\ -2c_2\left(h_{yx}h-h_y{h}_x\right)-2c_3\left(h_{yy}h-h_y{h}_y\right)\\ =6c_1{h}_{yx}{h}_{xx}+2\left(h_{yt}-c_2{h}_{yx}-c_3{h}_{yy}\right)h+2\left(-h_t+c_2{h}_x+c_3{h}_y\right)h_y.\end{array}$ (26)

From $h_{yx}=0$, $h_{yt}-c_2{h}_{yx}-c_3{h}_{yy}=0$, $-h_t+c_2{h}_x+c_3{h}_y=0$, and relational expression

${\displaystyle {\sum }_{i=1}^N}{a}_i{b}_i=0,{\displaystyle {\sum }_{i=1}^N}{b}_i{d}_i-c_3{\displaystyle {\sum }_{i=1}^N}{b}_i^2=0,d_i=c_2{a}_i+c_3{b}_i.$ (27)

Obviously $g_y=0$, so $P\left(D_x,D_y{D}_t\right)\left(g\cdot g\right)=0$. So that

$\begin{array}{l}2\left(D_t+c_1{D}_x^3-c_2{D}_x-c_3{D}_y\right)D_y\left(h\cdot g\right)\\ =2\left(D_t+c_1{D}_x^3-c_2{D}_x-c_3{D}_y\right)D_y\left(h_{y}g-h{g}_y\right)\\ =2\left(h_{yt}g-h_y{g}_t\right)+2c_1\left(h_{yxx}g-3h_{yxx}{g}_x+3h_{yx}{g}_{xx}-h_y{g}_{xxx}\right)\\ -2c_2\left(h_{yx}g-h_y{g}_x\right)-2c_3\left(h_{yy}g-h_y{g}_y\right)\\ =6c_1{h}_{yx}{g}_{xx}+2\left(h_{yt}-c_2{h}_{yx}-c_3{h}_{yy}\right)g-2\left(g_t+c_1{g}_{xxx}-c_2{g}_x\right)h_y.\end{array}$ (28)

From $g_t+c_1{g}_{xxx}-c_2{g}_x=0$, we get ${\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j\left(c_2{p}_j-c_1{p}_j^3\right){\text{e}}^{l_j}+c_1{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j{p}_j^3{\text{e}}^{l_j}-c_2{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j{p}_j{\text{e}}^{l_j}=0$, then

$n_j=\left(c_2{p}_j-c_1{p}_j^3\right)\mathrm{,}$ (29)

where $l_j=p_{j}x+\left(c_2{p}_j-c_1{p}_j^3\right)t+s_j$. Substituting (27), (29) into (25) completes the proof of the theorem.

Corollary 1. Take

$\begin{array}{c}f\left(x,y,t\right)=a_0+{\displaystyle \underset{i=1}{\overset{M}{\sum }}}{\left(a_{i}x+b_{i}y+d_{i}t+m_i\right)}^2\\ +{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j\cosh \left(p_{j}x+\left(c_2{p}_j-c_1{p}_j^3\right)t+s_j\right)\\ +{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\beta }_l\cos \left(p_{l}x+\left(c_2{p}_l-c_1{p}_l^3\right)t+s_l\right),\end{array}$ (30)

when (30) satisfies ${\displaystyle {\sum }_{i=1}^N}{a}_i{b}_i=0$, $d_i=c_2{a}_i+c_3{b}_i$, ${\displaystyle {\sum }_{i=1}^N}{b}_i{d}_i-c_3{\displaystyle {\sum }_{i=1}^N}{b}_i^2=0$, then $h\left(x\mathrm{,}y\mathrm{,}t\right)$ is the solution of (17), and $f\left(x\mathrm{,}y\mathrm{,}t\right)$ is also the solution of (17).

Secondly, we will apply the results of Theorem 1 to the (1). When $M=2,N=0$ in (24), two sets of relations can be obtained by direct calculation.

$a_0=a_0,a_1=a_1,a_2=a_2,b_1=0,b_2=0,d_1=d_1,d_2=d_2,m_1=m_1,m_2=m_2,$ (31)

and

$\begin{array}{l}{a}_0=a_0,a_1=a_1,a_2=-\frac{a_1{b}_1}{b_2},b_1=b_1,b_2=b_2,d_1=a_1{c}_2+b_1{c}_3,\\ d_2=-\frac{a_1{b}_1{c}_2}{b_2}+b_2{c}_3,m_1=m_1,m_2=m_2.\end{array}$ (32)

According to (16), (24) and (32), we obtain the interaction solution between lump type solution and N-soliton solution (See Figure 2 and Figure 3).

$u=-\frac{-4a_1^2\left(1+\frac{b_1^2}{b_2^2}\right)\left(x+c_{2}t\right)+{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\delta }_j{p}_j{\text{e}}^{\left(p_{j}x+q_{j}y+n_{j}t+s_j\right)}}{a_0+\left(1+\frac{b_1^2}{b_2^2}\right)(a_1^2{\left(x+c_{2}t\right)}^2+{\left(x+c_{3}y\right)}^2+g}.$ (33)

Figure 2 shows the evolution of the spatial structure of the interaction solution with the increase in the number of solitons. Here, we select parameters in Figure 2, (a): $a_0=10$, $a_1=-2$, $b_1=2$, $b_2=\frac{1}{2}$, $c_1=1$, $c_2=-2$, $c_3=4$, $m_1=m_2=0$. (b), (c), (d): $a_0=4$, $a_1=-2$, $b_1=2$, $b_2=\frac{1}{2}$, $c_1=1$, $c_2=-2$, $c_3=4$, $p_1=\frac{1}{4}$,

(a) (b) (c) (d)

Figure 2. Spatial evolution behavior of interaction solution (33) with soliton number: (a) $N=0$ ; (b) $N=1$ ; (c) $N=2$ ; (d) $N=3$.

(a) (b) (c) (d)

Figure 3. When $N=2$, the evolution behavior of the interaction solution (24) with the change of time t. (a) t = −20 (b) t = −10 (c) t = 10 (d) t = 20.

$p_2=\frac{1}{2}$, $p_3=-\frac{1}{5}$, ${\delta }_1=1$, $m_1=m_2=s_1=s_2=s_3=t=0$. Besides, we also give the evolution behavior of the interaction solution with time when $N=2$, from Figure 3, we find that with the development of time, the lump type solution is swallowed up by 2-solitons.

5. Conclusion

In summary, according to the bilinear form of ePBLMP (1), we constructed different test functions and obtained the double breather solution of the ePBLMP equation. The lump type solution is obtained from the breather solution by using the parameter limit method. Meanwhile, we also studied the interaction between lump type solutions and N-solitons. The spatial evolution behavior is demonstrated by images. More importantly, we give the theoretical calculation process of the interaction solution between lump type solution and N-solitons, and the correctness of the theory is verified by a concrete example. We believe that there are many nonlinear partial differential systems, although not integrable, also have similar dynamic behavior, such as high-dimensional Sine-Gordon equation, Ablowitz-Kaup-Newell-Segur systems, etc. we will focus on this in the future. At the same time, we also hope that the dynamic behavior obtained in this paper can provide new perspectives and useful information for the field of mathematical physics.

Acknowledgements

This work was partially supported by the Natural Science Foundation of Hunan Province (No. 2021JJ40434), the Scientific Research Project of the Hunan Education Department (No. 21B0510), and Jishou University Natural Science Foundation (No. Jdy21016).

Conflicts of Interest

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

References

[1]	Ma, W.M. (1997) Darboux Transformations for a Lax Integrable System in 2n Dimensions. Letters in Mathematical Physics, 39, 33-49. https://doi.org/10.1007/s11005-997-3049-3
[2]	Lu, Q., Ilhan, O.A., Manafian, J., et al. (2021) Multiple Rogue Wave Solutions for a Variable-Coefficient Kadomtsev-Petviashvili Equation. International Journal of Computer Mathematics, 98, 1457-1473. https://doi.org/10.1080/00207160.2020.1822996
[3]	Hirota, R. and Satsuma, J. (1981) Soliton Solutions of a Coupled Korteweg-de Vries Equation. Physics Letters A, 85, 407-408. https://doi.org/10.1016/0375-9601(81)90423-0
[4]	Hirota, R. (1980) Direct Methods in Soliton Theory. Physics Letters A, 379, 1975-1978.
[5]	Ma, W.X. (2016) Lump-Type Solutions to the (3+1)-Dimensional Jimbo-Miwa Equation. International Journal of Nonlinear Sciences and Numerical Simulation, 17, 355-359. https://doi.org/10.1515/ijnsns-2015-0050
[6]	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
[7]	Jimbo, M. and Miwa, T. (1983) Solitons and Infinite Dimensional Lie Algebras. Publications of the Research Institute for Mathematical Sciences, 19, 943-1001. https://doi.org/10.2977/prims/1195182017
[8]	Wang, D.S., Li, X.G., Chan, C.K. and Zhou, J. (2016) Double Wronskian Solution and Soliton Properties of the Nonisospectral BKP Equation. Communications in Theoretical Physics, 65, 259-265. https://doi.org/10.1088/0253-6102/65/3/259
[9]	Ma, W.X. and Zhou, Y. (2016) Lump Solutions to Nonlinear Partial Differential Equations via Hirota Bilinear Forms. Journal of Differential Equations, 264, 2633-2659.
[10]	Tan, W. (2019) Evolution of Breathers and Interaction between High-Order Lump Solutions and N-Solitons (N→∞) for Breaking Soliton System. Physics Letters A, 383, Article ID: 125907. https://doi.org/10.1016/j.physleta.2019.125907
[11]	Tan, W. and Li, M. (2021) Breather Degeneration and Lump Superposition for the (3+1)-Dimensional Nonlinear Evolution Equation. Modern Physics Letters B, 35, Article ID: 2150250. https://doi.org/10.1142/S021798492150250X
[12]	Ren, B. (2021) Dynamics of a D’Alembert Wave and a Soliton Molecule for an Extended BLMP Equation. Communications in Theoretical Physics, 73, Article ID: 035003. https://doi.org/10.1088/1572-9494/abda17
[13]	Gilson, C., Nimmo, J. and Willox, R. (1993) A (2+1)-Dimensional Generalization of the AKNS Shallow Water Wave Equation. Physics Letters A, 180, 337-345. https://doi.org/10.1016/0375-9601(93)91187-A
[14]	Darvishi, M.T., Najafi, M., Kavitha, L. and Venkatesh, M. (2012) Stair and Step Soliton Solutions of the Integrable (2+1) and (3+1)-Dimensional Boiti-Leon-Manna-Pempinelli Equations. Communications in Theoretical Physics, 58, 785-794. https://doi.org/10.1088/0253-6102/58/6/01
[15]	Shen, Y., Tian, B., Zhang, C.R., Tian, H.Y. and Liu, S.H. (2021) Breather-Wave, Periodic-Wave and Traveling-Wave Solutions for a (2+1)-Dimensional Extended Boiti-Leon-Manna-Pempinelli Equation for an Incompressible Fluid. Modern Physics Letters B, 35, Article ID: 2150261. https://doi.org/10.1142/S0217984921502614
[16]	Song, L.L., Pu, Z.L. and Dai, Z.D. (2017) Spatio-Temporal Deformation of Kink-Breather to the (2+1)-Dimensional Potential Boiti-Leon-Manna-Pempinelli Equation. Communications in Theoretical Physics, 67, 493. https://doi.org/10.1088/0253-6102/67/5/493
[17]	Seadawy, A.R., Ali, A. and Helal, M.A. (2021) Analytical Wave Solutions of the (2+1)-Dimensional Boiti-Leon-Pempinelli and Boiti-Leon-Manna-Pempinelli Equations by Mathematical Methods. Mathematical Methods in the Applied Sciences, 44, 14292-14315. https://doi.org/10.1002/mma.7697
[18]	Paliathanasis, A. (2022) Lie Symmetry Analysis for a (2+1) Extended Boiti-Leon-Manna-Pempinelli Equation. Quaestiones Mathematicae, 1-8. https://doi.org/10.2989/16073606.2022.2035844
[19]	Luo, Y. (2019) New Soliton Solutions to the Initial Value Problem for the Two-Component Short Pulse Equation. Journal of Applied Mathematics and Physics, 7, 13-22. https://doi.org/10.4236/jamp.2019.71002
[20]	Zeng, Z., Liu, X., Zhu, Y., et al. (2022) New Exact Traveling Wave Solutions of (2+1)-Dimensional Time-Fractional Zoomeron Equation. Journal of Applied Mathematics and Physics, 10, 333-346. https://doi.org/10.4236/jamp.2022.102026
[21]	Tan, W., Dai, Z.D., Xie, J.L. and Qiu, D.Q. (2018) Parameter Limit Method and Its Application in the (4+1)-Dimensional Fokas Equation. Computers and Mathematics with Applications, 75, 4214-4220. https://doi.org/10.1016/j.camwa.2018.03.023
[22]	Dai, Z.D., Wang, C.J. and Liu, J. (2014) Inclined Periodic Homoclinic Breather and Rogue Waves for the (1+1)-Dimensional Boussinesq Equation. Pramana—Journal of Physics, 83, 473-480. https://doi.org/10.1007/s12043-014-0811-9
[23]	Tan, W. (2021) Some New Dynamical Behaviour of Double Breathers and Lump-N-Solitons for the Ito Equation. International Journal of Computer Mathematics, 98, 961-974. https://doi.org/10.1080/00207160.2020.1792454