Rational Solutions and Their Interaction Solutions of the (3 + 1)-Dimensional Jimbo-Miwa Equation

In this paper, we gave a form of rational solution and their interaction solution to a nonlinear evolution equation. The rational solution contained lump solution, general lump solution, high-order lump solution, lump-type solution, etc. Their interaction solution contained the classical interaction solution, such as the lump-kink solution and the lump-soliton solution. As the example, by using the generalized bilinear method and symbolic computation Maple, we obtained abundant high-order lump-type solutions and their interaction solutions between lumps and other function solutions under certain constraints of the (3 + 1)-dimensional Jimbo-Miwa equation. Via three-dimensional plots, contour plots and density plots with the help of Maple, the physical characteristics and structures of these waves are described very well. These solutions have greatly enriched the exact solutions of the (3 + 1)-dimensional Jimbo-Miwa equation on the existing literature.


Introduction
Nonlinear phenomena have a lot of significant applications in different sides of physics with natural and engineering fields. Basically, all the fundamental equations of physics are nonlinear and, generally, such types of nonlinear evolution equations (NLEEs) are often very tough to solve clearly. The exact solutions of NLEEs play a crucial role in the study of nonlinear physical or natural phenomena. In the recent decade, several direct methods for finding the exact solutions to NLEEs have been proposed [1][2][3][4][5][6][7][8][9]. Thousands of examples have shown that these methods are powerful for obtaining exact solutions of NLEEs, such as soliton [10][11][12][13][14], rogue wave [15,16], breather solution [17], periodic wave solution [18][19][20][21], and optical solution [22,23].
The lump solution has attracted a great deal of attention since lump solutions were firstly discovered [24]. The research to lump solution has not been well developed, because it is very complex to solve the lump solution of NLEEs. Recently, based on the Hirota bilinear method, Ma and Zhou introduced a new way to get the lump solution of NLEEs by using symbolic computation and gave a theoretical testimony [25,26]. By using this method, researchers successfully obtained the lump solutions and interaction solutions of NLEEs . In the present paper, we will propose the form of rational solution and their interaction solution to NLEE. The rational solution contains lump solution, general lump solution, high-order lump solutions, lump-type solution, etc. Their interaction solution contains the classical interaction solution, such as the lump-kink solution and the lump-soliton solution.
The rest of the paper is organized as follows. In Section 2, we will give the form of rational solution and their interaction solution to NLEE. In Section 3, by using the generalized bilinear method and symbolic computation Maple, we will obtain the high-order lump-type solutions of the (3 + 1)dimensional Jimbo-Miwa equation. In Section 4, by using the symbolic computation Maple, we will get abundant interaction solutions between the high-order lump-type solution and other function solutions. Via three-dimensional plots, contour plots, and density plots with the help of Maple, the physical characteristics and structures of these waves are described very well. In Section 5, a few of the conclusions and outlook will be given.

Rational Solution and Their Interaction Solution
Consider a Kth order NLEE ðK ≥ 2Þ F x, u, ∂u, ∂ 2 u, ⋯, where x = ðx 1 , x 2 ,⋯,x n Þ are n independent variables and x i ði ≠ 1Þ contain time variable t. u is the dependent variable.

Rational Solution.
In order to get the rational solution of NLEE (1), we take its main steps as follows.
Step 3. By substituting (4) and (5) into Equation (3), collecting all terms with the same order of x i together, the left-hand side of Equation (3) is converted into another polynomial in x i . Equating each coefficient of this different power terms to zero yields a set of nonlinear algebraic equations for a 0 , a ij . With the aid of Maple (or Mathematica), we solve the above nonlinear algebraic equations.

General Interaction Solution.
In order to obtain the general interaction solution, we take its main steps as follows: Step 1. By using transformation (2), Equation (1) is transformed into bilinear form (3).

Bilinear Form. Under the Cole-Hopf transformation,
Equation (8) becomes the generalized Hirota bilinear equation: where p, being an arbitrarily natural number, is often a prime number. D is a generalized bilinear differential operator as follows [3]: where m, n ≥ 0, α s p = ð−1Þ r p ðsÞ , if s ≡ r p ðsÞ mod p. When taking p = 2, we obtain the Hirota bilinear equation: When taking p = 3, we can obtain the generalized bilinear Jimbo-Miwa equation: By using transformation (9), generalized bilinear Jimbo-Miwa equation (13) is transformed into the following form: where u y = v x . Transformation (9) is also a characteristic one in establishing Bell polynomial theories of soliton equations [59], and an accurate relation is Hence, if f solves generalized bilinear Jimbo-Miwa equation (13), Jimbo-Miwa equation (14) will be solved.

High-Order Lump-Type Solutions.
In the section, we will study the high-order lump-type solutions of (3 + 1)-dimensional Jimbo-Miwa equation (8) by constructing positive quadratic function solutions to the corresponding generalized bilinear equation (13).
Step 2. To get the positive quadratic function solution of generalized bilinear equation (13), we take N = 3, n 1 = 2, n 2 = 1, n 3 = 1 in expression (4), where where a ij ði = 1, 2, 3 ; j = 0, 1, 2, 3, 4Þ are arbitrary real constants, and Step 3. By substituting (16) and (17) into Equation (13), 3 Advances in Mathematical Physics collecting all terms with the same order of x, y, z, t together, the left-hand side of Equation (13) is converted into another polynomial in x, y, z, t. Equating each coefficient of this different power terms to zero yields a set of nonlinear algebraic equations for a 0 , a ij .
Solving the algebraic equations by Maple yields the following sets of solutions, Case 1. where other parameters in Cases 1-3 are arbitrary real constants.
where f i ðx, y, z, tÞ are the positive quadratic function solutions to the generalized bilinear Jimbo-Miwa equation (13), and It is also readily observed that at any given time t, the above high-order lump-type solutions u i → 0 if and only if the corresponding sum of squares ξ 4 1 In order to exhibit the dynamical characteristics of these waves, we plot various three-dimensional, contour, and density plots as follows. We choose the following parameters to illustrate the high-order lump-type solution u 2 ðx, y, z, tÞ for Jimbo-Miwa equation (8), The physical properties and structures for the highorder lump-type solution u 2 ðx, y, z, tÞ are described in Figure 1. Figure 1 shows the three-dimensional dynamic graphs A 1 , B 1 , C 1 , corresponding contour maps A 2 , B 2 , C 2 , and density plots A 3 , B 3 , C 3 in the ðx, yÞ plane when t = −6, 0, 6, respectively. The three-dimensional graphs reflect the localized structures, and the density plots show the energy distribution. Remark 6. By substituting a 0 , a ij in Cases 4-7 to (16) and (17) and using bilinear transformation (9), we obtain the new lump solutions for Jimbo-Miwa equation (8). Due to the lack of space, we omit the expressions of lump solutions.

Interaction Solutions between Lump and Soliton Solutions of the (3 + 1)-Dimensional Jimbo-Miwa Equation
In this section, we will study the general interaction solutions between the high-order lump-type solutions and other function solutions of (3 + 1)-dimensional Jimbo-Miwa equation (8).

Advances in Mathematical Physics
where m j , b jk ðj = 1, 2, 3, 4 ; k = 0, 1, 2, 3, 4Þ are arbitrary real constants. In order to obtain the interaction solution between the high-order lump-type solution and the double exponential function, the trigonometric function, and the hyperbolic function of (3 + 1)-dimensional Jimbo-Miwa equation (8), we suppose The interaction solution of generalized bilinear equation (13) is written the following form: Step 3. By substituting (33) into Equation (13), collecting all terms with the same order of x, y, z, t, e η 1 , e −η 2 , tan η 3 , tanh η 4 together, the left-hand side of Equation (13) is converted into another polynomial in x, y, z, t, e η 1 , e −η 2 , tan η 3 , tanh η 4 . Equating each coefficient of this different power terms to zero yields a set of nonlinear algebraic equations for a 0 , a ij , b jk , m j . Solving the algebraic equations by Maple, yields the following sets of solutions.
where other parameters are arbitrary real constants.
Remark 9. In addition to the above result Case 1, we can also get the same solutions as Cases 1-6 in Section 4.3 when m 3 = m 1 , η 3 = η 4 , respectively.
where other parameters are arbitrary real constants.

Advances in Mathematical Physics
Remark 10. In addition to the above results Cases 1-6, we can also get the same solutions as Cases 1-3 in Section 4.3 when η 3 = η 1 , respectively.
Step 4. By substituting the parameters a 0 , a ij , b jk , m j in the Sections 4.1-4.6 into the solution (33) and using transformation (9), we can obtain abundant interaction solutions of Jimbo-Miwa equation (1).
These sets of solutions for the parameters generate 42 classes of combination solutions f i , 1 ≤ i ≤ 42 to the generalized bilinear Jimbo-Miwa equation (13), and then, the resulting combination solutions present 42 classes of interaction solutions u i , 1 ≤ i ≤ 42 to Equation (8) under transformation (9). Therefore, various kinds of interaction solutions could be constructed explicitly this way.
As the example, substituting (38) into (33), we can get f as follows: where By using transformation (9), we get the interaction solution between the high-order lump-type solutions and a pair of line-soliton solution of (3 + 1)-dimensional Jimbo-Miwa equation (8) where f ðx, y, z, tÞ is given in (62). In order to exhibit the dynamical characteristics of these waves, we plot various three-dimensional, contour, and density plots as follows. We choose the following parameters to illustrate interaction solution (64), The physical properties and structures for interaction solution (64) are shown in Figure 2. Figure 2 shows the three-dimensional dynamic graphs A 1 , B 1 , C 1 , corresponding contour maps A 2 , B 2 , C 2 , and density plots A 3 , B 3 , C 3 in the ðx, yÞ-plane when t = −1,0, and 1, respectively. The threedimensional graphs reflect the localized structures, and the density plots show the energy distribution. We can see that the high-order lump-type wave and the exponential function wave react with each other.
When we choose the following parameters and t = −2, we illustrate interaction solution (64) of (3 + 1)-dimensional Jimbo-Miwa equation (8), The physical properties and structures for interaction solution (64) are shown in Figure 3. Figure 3 shows the three-dimensional dynamic graph D 1 , corresponding contour map D 2 , and density plot D 3 in the ðx, yÞ-plane, respectively.

Conclusion
In this paper, we gave the form of rational solution and their interaction solution to NLEE. The rational solution contained lump solution, general lump solution, high-order lump solution, lump-type solution, etc. The general interaction solution contain the classical interaction solution, such as the lump-kink solution and the lump-soliton solution. As the example, by using the generalized bilinear method and symbolic computation Maple, we successfully constructed 14 Advances in Mathematical Physics  the high-order lump-type solutions and their interaction solutions between lumps and other function solutions under certain constraints of (3 + 1)-dimensional Jimbo-Miwa equation. Three-dimensional plots, contour plots, and density plots of these waves are observed in Figures 1-3, respectively. We can find the physical structure and characteristics of the interactions between the high-order lump-type solutions and the exponential function wave.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that they have no conflicts of interest.