Study on the Nonlinear Dynamics of the (3+1)-Dimensional Jimbo-Miwa Equation in Plasma Physics

Abstract

The Jimbo-Miwa equation (JME) that describes certain interesting (3+1)-dimensional waves in plasma physics is studied in this work. The Hirota bilinear equation is developed via the Cole-Hopf transform. Then, the symbolic computation, together with the ansatz function schemes, are utilized to seek exact solutions. Some new solutions, such as the multi-wave complexiton solution (MWCS), multi-wave solution (MWS) and periodic lump solution (PLS), are successfully constructed. Additionally, different types of travelling wave solutions (TWS), including the dark, bright-dark and singular periodic wave solutions, are disclosed by employing the sub-equation method. Finally, the physical characteristics and interaction behaviors of the extracted solutions are depicted graphically by assigning appropriate parameters. The obtained outcomes in this paper are more general and newer. Additionally, they reveal that the used methods are concise, direct, and can be employed to study other partial differential equations (PDEs) in physics.

1. Introduction

Complex phenomena in engineering and physics can usually be reduced to PDEs [1,2,3,4,5,6]. The study on the properties of these equations such as the explicit analytical solutions, especially the soliton solutions, is of great significance since they can help us to better understand complex phenomena and their inner nature. Up to now, a series of different effective methods have been developed to construct the exact solutions of PDEs such as the Hirota bilinear method [7,8,9,10], Wang’s Bäcklund transformation-based method [11,12], trial equation method [13,14], Sardar subequation method [15,16,17], exp-function method [18,19], Riccati equation mapping method [20] and so on [21,22,23,24,25,26,27,28]. In this work, we aim to examine the (3+1)-dimensional JME given by [29]:

$${\mathrm{\Pi }}_{xxxy}+3{\mathrm{\Pi }}_x{\mathrm{\Pi }}_{xy}+3{\mathrm{\Pi }}_y{\mathrm{\Pi }}_{xx}+2{\mathrm{\Pi }}_{yt}-3{\mathrm{\Pi }}_{xz}=0,$$

(1)

Equation (1) is derived from the second equation in the well-known KP hierarchy of integrable systems and used widely to describe some interesting (3+1)-dimensional waves in plasma and optics. Up to now, some important research achievements have been developed to deal with Equation (1). In [29], the Kudryashov method is used with the symbolic computation and different solutions are obtained. In [30], four kinds of different wave forms are constructed via the Hirota bilinear method. In [31], the authors employ the direct algebraic method to handle Equation (1) and some different wave forms are constructed. In [32], several closed-form solutions are developed by using the singular manifold method. In [33], the Riccati equation mapping method is adopted. The exp-function method is utilized in [34] and some generalized solutions with parameters are constructed. In [35], the authors carry out the linear superposition principle to seek for multi-resonant solutions of Equation (1). In [36], the authors make use of the generalized Bernoulli equation method to inquire into Equation (1). In this study, we will present th results of a detailed investigation of Equation (1). The rest of the content of this work is given as follows. In Section 2, the Cole-Hopf transform is adopted to establish the Hirota bilinear form, and symbolic computation, combined with the ansatz function schemes, is utilized to search for the MWCS, MWS and PLS. In Section 3, the sub-equation method is used to seek for the TWSs. In Section 4, the physical characteristics and interaction behaviors are presented. Finally, we reach a conclusion in Section 5.

2. The Hirota Bilinear Equation and the Exact Solutions

To obtain the Hirota bilinear form of Equation (1), we adopt the Cole-Hopf transform as:

$$\mathrm{\Pi }=2\ln {\left(\mathrm{\Xi }\right)}_x,$$

(2)

Taking it into Equation (1), we can obtain the bilinear form as:

$$\left(D_x^3{D}_y+2D_y{D}_t-3D_x{D}_z\right)\mathrm{\Xi }\cdot \mathrm{\Xi }=0.$$

(3)

Here, the definition of the operators $D_x^m{D}_{\tau }^n$ is [37,38]:

$$D_x^m{D}_t^{n}f\cdot g={\left(\frac{𝜕}{𝜕x}-\frac{𝜕}{𝜕x^{\prime }}\right)}^m{\left(\frac{𝜕}{𝜕t}-\frac{𝜕}{𝜕t^{\prime }}\right)}^{n}f\left(x,t\right)g\left(x^{\prime },t^{\prime }\right){|}_{x=x^{\prime },t=t^{\prime }}.$$

(4)

Additionally, there are

$$D_x\left(f\cdot g\right)=f_{x}g-f{g}_x,$$

$$D_x^2\left(f\cdot g\right)=f_{xx}g-2f_x{g}_x+f{g}_{xx},$$

$$D_x^2\left(f\cdot f\right)=2\left(f_{xx}f-f_x^2\right),$$

$$D_t{D}_x\left(f\cdot g\right)=f_{tx}g-f_t{g}_x-f_x{g}_t+g_{tx}f.$$

2.1. The MWCS

In order to find the MWCS, it is assumed that the solution of Equation (3) is:

$$\mathrm{\Xi }=u_1{e}^p+u_2{e}^{-p}+u_3\sin \left(q\right)+u_4\sin \mathrm{h}\left(\rho \right),$$

(5)

with

$$\{\begin{array}{l}p=x+k_{1}y+k_{2}z+k_{3}t\\ q=x+k_{4}y+k_{5}z+k_{6}t\\ \rho =x+k_{7}y+k_{8}z+k_{9}t\end{array},$$

where $u_i$($i=1,2,3,4.$) and $k_i$($i=1,2,3,4,5,6,7,8,9.$) are constants that can be determined later. Substituting Equation (5) into Equation (3) and setting the coefficients of different terms to zero, an algebraic equation system is attained. Solving it, we derive:

Case 1:

$$\begin{gathered} k_1=\frac{3k_2}{2\left(k_9+2\right)}, k_2=k_2, k_3=k_9, k_4=-\frac{3k_2}{2\left(k_9+2\right)}, k_5=\frac{k_2-k_2{k}_9}{2+k_9}, k_6=k_9+1, k_7=\frac{3k_2}{2\left(2+k_9\right)}, \\ {k}_8=k_2, k_9=k_9, u_1=u_1, u_2=u_2, u_3=u_3, u_4=u_4. \end{gathered}$$

The MWCS is obtained as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[\begin{array}{l}{u}_1{e}^{x+\frac{3k_2}{2\left(k_9+2\right)}y+k_{2}z+k_{9}t}-u_2{e}^{-\left(x+\frac{3k_2}{2\left(k_9+2\right)}y+k_{2}z+k_{9}t\right)}+u_3\cos \left(x-\frac{3k_2}{2\left(k_9+2\right)}y+\frac{k_2-k_2{k}_9}{2+k_9}z+\left(k_9+1\right)t\right)\\ +u_4\cos \mathrm{h}\left(x+\frac{3k_2}{2\left(2+k_9\right)}y+k_{2}z+k_{9}t\right)\end{array}\right]}{\begin{array}{l}{u}_1{e}^{x+\frac{3k_2}{2\left(k_9+2\right)}y+k_{2}z+k_{9}t}+u_2{e}^{-\left(x+\frac{3k_2}{2\left(k_9+2\right)}y+k_{2}z+k_{9}t\right)}+u_3\sin \left(x-\frac{3k_2}{2\left(k_9+2\right)}y+\frac{k_2-k_2{k}_9}{2+k_9}z+\left(k_9+1\right)t\right)\\ +u_4\sin \mathrm{h}\left(x+\frac{3k_2}{2\left(2+k_9\right)}y+k_{2}z+k_{9}t\right)\end{array}}.$$

(6)

For the special case $u_1=-u_2=2u_4$, Equation (6) becomes:

Case 2:

$$\begin{gathered} k_1=k_1, k_2=k_2, k_3=-2+\frac{3k_2}{2k_1}, k_4=-k_1, k_5=2k_1-k_2, k_6=-1+\frac{3k_2}{2k_1}, k_7=k_1, k_8=k_2, \\ {k}_9=-2+\frac{3k_2}{2k_1}, u_1=u_1, u_2=u_2, u_3=u_3, u_4=u_4. \end{gathered}$$

Thus, we can obtain the MWCS as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[\begin{array}{l}{u}_1{e}^{x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t}-u_2{e}^{-\left(x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t\right)}+u_3\cos \left(x-k_{1}y+\left(2k_1-k_2\right)z+\left(-1+\frac{3k_2}{2k_1}\right)t\right)\\ +u_4\cos \mathrm{h}\left(x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t\right)\end{array}\right]}{\begin{array}{l}{u}_1{e}^{x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t}+u_2{e}^{-\left(x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t\right)}+u_3\sin \left(x-k_{1}y+\left(2k_1-k_2\right)z+\left(-1+\frac{3k_2}{2k_1}\right)t\right)\\ +u_4\sin \mathrm{h}\left(x+k_{1}y+k_{2}z+\left(-2+\frac{3k_2}{2k_1}\right)t\right)\end{array}}.$$

(7)

Case 3:

$$\begin{gathered} k_1=k_1, k_2=\frac{2}{3}\left(2k_1+k_1{k}_9\right), k_3=k_9, k_4=-k_1, k_5=-\frac{2}{3}\left(-k_1+k_1{k}_9\right), k_6=-1+k_9, k_7=k_1, \\ {k}_8=\frac{2}{3}\left(2k_1+k_1{k}_9\right), k_9=k_9, u_1=u_1, u_2=u_2, u_3=u_3, u_4=u_4. \end{gathered}$$

Thus, we obtain the MWCS solution as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[\begin{array}{l}{u}_1{e}^{x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{3}t}-u_2{e}^{-\left(x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{3}t\right)}+u_3\cos \left(x+k_{4}y+k_{5}z+\left(k_9-1\right)t\right)\\ +u_4\cos \mathrm{h}\left(x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{9}t\right)\end{array}\right]}{\begin{array}{l}{u}_1{e}^{x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{3}t}+u_2{e}^{-\left(x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{3}t\right)}+u_3\sin \left(x+k_{4}y+k_{5}z+\left(k_9-1\right)t\right)\\ +u_4\sin \mathrm{h}\left(x+k_{1}y+\frac{2}{3}\left(2k_1+k_1{k}_9\right)z+k_{9}t\right)\end{array}}.$$

(8)

Case 4:

$$\begin{array}{c}{k}_1=-\frac{3k_5}{2\left(k_9-1\right)}, k_2=-\frac{2k_5\left(1+k_9\right)}{k_9-1}, k_3=k_9, k_4=\frac{3k_5}{2\left(k_9-1\right)}, k_5=k_5, k_6=k_9-1, \\ k_7=-\frac{3k_5}{2\left(k_9-1\right)}, k_8=-\frac{2k_5\left(1+k_9\right)}{k_9-1}, k_9=k_9, u_1=u_1, u_2=u_2, u_3=u_3, u_4=u_4.\end{array}$$

Accordingly, the MWCS is:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[\begin{array}{l}{u}_1{e}^{x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t}-u_2{e}^{-\left(x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t\right)}+u_3\cos \left(x+\frac{3k_5}{2\left(k_9-1\right)}y+k_{5}z+\left(k_9-1\right)t\right)\\ +u_4\cos \mathrm{h}\left(x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t\right)\end{array}\right]}{\begin{array}{l}{u}_1{e}^{x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t}+u_2{e}^{-\left(x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t\right)}+u_3\sin \left(x+\frac{3k_5}{2\left(k_9-1\right)}y+k_{5}z+\left(k_9-1\right)t\right)\\ +u_4\sin \mathrm{h}\left(x-\frac{3k_5}{2\left(k_9-1\right)}y-\frac{2k_5\left(1+k_9\right)}{k_9-1}z+k_{9}t\right)\end{array}}.$$

(9)

Case 5:

$$\begin{gathered} k_1=-k_4, k_2=k_8, k_3=-\left(2+\frac{3k_8}{2k_4}\right), k_4=k_4, k_5=-\left(2k_4+k_8\right), k_6=-\left(1+\frac{3k_8}{2k_4}\right), k_7=-k_4, \\ {k}_8=k_8, k_9=-\left(2+\frac{3k_8}{2k_4}\right), u_1=u_1, u_2=u_2, u_3=u_3, u_4=u_4. \end{gathered}$$

where $k_4\ne 0$. Thus, we can obtain the MWCS as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[\begin{array}{l}{u}_1{e}^{x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t}-u_2{e}^{-\left(x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t\right)}+u_3\cos \left(x+k_{4}y-\left(2k_4+k_8\right)z-\left(1+\frac{3k_8}{2k_4}\right)t\right)\\ +u_4\cos \mathrm{h}\left(x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t\right)\end{array}\right]}{\begin{array}{l}{u}_1{e}^{x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t}+u_2{e}^{-\left(x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t\right)}+u_3\sin \left(x+k_{4}y-\left(2k_4+k_8\right)z-\left(1+\frac{3k_8}{2k_4}\right)t\right)\\ +u_4\sin \mathrm{h}\left(x-k_{4}y+k_{8}z-\left(2+\frac{3k_8}{2k_4}\right)t\right)\end{array}}.$$

(10)

For the special case $u_1=-u_2=2u_4$, Equations (6)–(10) become:

$$\mathrm{\Pi }\left(x,y,z,t\right)=2\cot \left(x-\frac{3k_2}{2\left(k_9+2\right)}y+\frac{k_2-k_2{k}_9}{2+k_9}z+\left(k_9+1\right)t\right).$$

$$\mathrm{\Pi }\left(x,y,z,t\right)=2\cot \left(x-k_{1}y+\left(2k_1-k_2\right)z+\left(-1+\frac{3k_2}{2k_1}\right)t\right).$$

$$\mathrm{\Pi }\left(x,y,z,t\right)=2\cot \left(x+k_{4}y+k_{5}z+\left(k_9-1\right)t\right).$$

$$\mathrm{\Pi }\left(x,y,z,t\right)=2\cot \left(x+\frac{3k_5}{2\left(k_9-1\right)}y+k_{5}z+\left(k_9-1\right)t\right).$$

$$\mathrm{\Pi }\left(x,y,z,t\right)=2\cot \left(x+k_{4}y-\left(2k_4+k_8\right)z-\left(1+\frac{3k_8}{2k_4}\right)t\right).$$

2.2. The MWS

Here, we can use the following ansatz function:

$$\mathrm{\Xi }=u_1\cos \left(p\right)+u_2\cos \mathrm{h}\left(q\right)+u_3\cos \mathrm{h}\left(\rho \right),$$

(11)

with

$$\{\begin{array}{l}p=x+k_{1}y+k_{2}z+k_{3}t\\ q=x+k_{4}y+k_{5}z+k_{6}t\\ \rho =x+k_{7}y+k_{8}z+k_{9}t\end{array},$$

where $u_i$($i=1,2,3.$) and $k_i$($i=1,2,3,4,5,6,7,8,9.$) are constants that can be determined later. In the same manner, substituting Equation (11) into Equation (3) and making the corresponding adjustments, we derive:

Case 1:

$$k_1=k_1, k_2=k_2, k_3=2+\frac{3k_2}{2k_1}, k_7=-k_1, k_8=-2k_1-k_2, k_9=1+\frac{3k_2}{2k_1}, u_1=u_1, u_2=0, u_3=u_3.$$

Then, we obtain the MWS as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[-u_1\sin \left(x+k_{1}y+k_{2}z+\left(2+\frac{3k_2}{2k_1}\right)t\right)+u_3\sin \mathrm{h}\left(x-k_{1}y-\left(2k_1+k_2\right)z+\left(1+\frac{3k_2}{2k_1}\right)t\right)\right]}{u_1\cos \left(x+k_{1}y+k_{2}z+\left(2+\frac{3k_2}{2k_1}\right)t\right)+u_3\cos \mathrm{h}\left(x-k_{1}y-\left(2k_1+k_2\right)z+\left(1+\frac{3k_2}{2k_1}\right)t\right)}.$$

(12)

Case 2:

$$k_1=k_1, k_2=-2k_1-k_5, k_3=-1-\frac{3k_5}{2k_1}, k_4=-k_1, k_5=k_5, k_6=-2-\frac{3k_5}{2k_1}, u_1=u_1, u_2=u_2, u_3=0.$$

Thus, we obtain the MWS as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[-u_1\sin \left(x+k_{1}y-\left(2k_1+k_5\right)z-\left(1+\frac{3k_5}{2k_1}\right)t\right)+u_2\sin \mathrm{h}\left(x-k_{1}y+k_{5}z-\left(2+\frac{3k_5}{2k_1}\right)t\right)\right]}{u_1\cos \left(x+k_{1}y-\left(2k_1+k_5\right)z-\left(1+\frac{3k_5}{2k_1}\right)t\right)+u_2\cos \mathrm{h}\left(x-k_{1}y+k_{5}z-\left(2+\frac{3k_5}{2k_1}\right)t\right)}.$$

(13)

2.3. The PLS

The solution of Equation (3) is assumed as:

$$\mathrm{\Xi }=u_1\sin \left(p\right)+u_2\cos \mathrm{h}\left(q\right)+k_7,$$

(14)

with

$$\{\begin{array}{l}p=x+k_{1}y+k_{2}z+k_{3}t\\ q=x+k_{4}y+k_{5}z+k_{6}t\end{array},$$

where $u_i$($i=1,2.$) and $k_i$($i=1,2,3,4,5,6,7.$) are constants to be determined later. In the same manner, substituting Equation (14) into Equation (3) and making the corresponding adjustments, we derive:

Case 1:

$$\begin{array}{l}{k}_1=k_1, k_2=\frac{2}{3}{k}_1\left(k_6-1\right), k_3=k_6+1, k_4=-k_1, k_5=-\frac{2}{3}{k}_1\left(2+k_6\right), k_6=k_6, k_7=0, u_1=u_1,\\ u_2=u_2.\end{array}$$

The PLS to Equation (1) is:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[u_1\cos \left(x+k_{1}y+\frac{2}{3}{k}_1\left(k_6-1\right)z+\left(k_6+1\right)t\right)+u_2\sin \mathrm{h}\left(x-k_{1}y-\frac{2}{3}{k}_1\left(2+k_6\right)z+k_{6}t\right)\right]}{u_1\sin \left(x+k_{1}y+\frac{2}{3}{k}_1\left(k_6-1\right)z+\left(k_6+1\right)t\right)+u_2\cos \mathrm{h}\left(x-k_{1}y-\frac{2}{3}{k}_1\left(2+k_6\right)z+k_{6}t\right)}.$$

(15)

Case 2:

$$\begin{array}{l}{k}_1=\frac{3k_2}{2\left(k_6-1\right)}, k_2=k_2, k_3=k_6+1, k_4=-\frac{3k_2}{2\left(k_6-1\right)}, k_5=-\frac{2k_2\left(k_6+1\right)}{k_6-1}, k_6=k_6, k_7=0, u_1=u_1, \\ u_2=u_2.\end{array}$$

Thus, we obtain the PLS of Equation (1) as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[u_1\cos \left(x+\frac{3k_2}{2\left(k_6-1\right)}y+k_{2}z+\left(k_6+1\right)t\right)+u_2\sin \mathrm{h}\left(x-\frac{3k_2}{2\left(k_6-1\right)}y-\frac{2k_2\left(k_6+1\right)}{k_6-1}z+k_{6}t\right)\right]}{u_1\sin \left(x+\frac{3k_2}{2\left(k_6-1\right)}y+k_{2}z+\left(k_6+1\right)t\right)+u_2\cos \mathrm{h}\left(x-\frac{3k_2}{2\left(k_6-1\right)}y-\frac{2k_2\left(k_6+1\right)}{k_6-1}z+k_{6}t\right)}.$$

(16)

Case 3:

$$\begin{array}{l}{k}_1=-\frac{3k_5}{2\left(k_3+1\right)}, k_2=\frac{k_5\left(2-k_3\right)}{k_3+1}, k_3=k_3, k_4=\frac{3k_5}{2\left(k_3+1\right)}, k_5=k_5, k_6=k_3-1, k_7=0, u_1=u_1, \\ u_2=u_2\end{array}$$

The PLS of Equation (1) is obtained as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[u_1\cos \left(x-\frac{3k_5}{2\left(k_3+1\right)}y+\frac{k_5\left(2-k_3\right)}{k_3+1}z+k_{3}t\right)+u_2\sin \mathrm{h}\left(x+\frac{3k_5}{2\left(k_3+1\right)}y+k_{5}z+\left(k_3-1\right)t\right)\right]}{u_1\sin \left(x-\frac{3k_5}{2\left(k_3+1\right)}y+\frac{k_5\left(2-k_3\right)}{k_3+1}z+k_{3}t\right)+u_2\cos \mathrm{h}\left(x+\frac{3k_5}{2\left(k_3+1\right)}y+k_{5}z+\left(k_3-1\right)t\right)}.$$

(17)

Case 4:

$$\begin{array}{l}{k}_1=-k_4, k_2=-\frac{2}{3}{k}_4\left(k_3-2\right), k_3=k_3, k_4=k_4, k_5=\frac{2}{3}{k}_4\left(k_3+1\right), k_6=k_3-1, k_7=0, u_1=u_1, \\ u_2=u_2.\end{array}$$

Thus, the PLS of Equation (1) is attained as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[u_1\cos \left(x-k_{4}y-\frac{2}{3}{k}_4\left(k_3-2\right)z+k_{3}t\right)+u_2\sin \mathrm{h}\left(x+k_{4}y+\frac{2}{3}{k}_4\left(k_3+1\right)z+\left(k_3-1\right)t\right)\right]}{u_1\sin \left(x-k_{4}y-\frac{2}{3}{k}_4\left(k_3-2\right)z+k_{3}t\right)+u_2\cos \mathrm{h}\left(x+k_{4}y+\frac{2}{3}{k}_4\left(k_3+1\right)z+\left(k_3-1\right)t\right)}.$$

(18)

Case 5:

$$k_1=k_1, k_2=k_2, k_3=2+\frac{3k_2}{2k_1}, k_4=-k_1, k_5=-2k_1-k_2, k_6=1+\frac{3k_2}{2k_1}, k_7=0, u_1=u_1, u_2=u_2.$$

We obtain the PLS of Equation (1) as:

$$\mathrm{\Pi }\left(x,y,z,t\right)=\frac{2\left[u_1\cos \left(x+k_{1}y+k_{2}z+\left(2+\frac{3k_2}{2k_1}\right)t\right)+u_2\sin \mathrm{h}\left(x-k_{1}y-\left(2k_1+k_2\right)z+\left(1+\frac{3k_2}{2k_1}\right)t\right)\right]}{u_1\sin \left(x+k_{1}y+k_{2}z+\left(2+\frac{3k_2}{2k_1}\right)t\right)+u_2\cos \mathrm{h}\left(x-k_{1}y-\left(2k_1+k_2\right)z+\left(1+\frac{3k_2}{2k_1}\right)t\right)}.$$

(19)

3. The TWS

This section aims to study the TWS using the sub-equation method [39,40]. For this end, we apply the following variable transformation to Equation (1):

$$\mathrm{\Pi }\left(x,y,z,t\right)=\Im \left(\chi \right), \chi =mx+ny+kz+st,$$

(20)

where m, n, k, and s are non-zero constants. Equation (1) can be converted as:

$$m^{3}n{\Im }^{\left(4\right)}+6m^{2}n{\Im }^{\prime }{\Im }^{″}+\left(2ns-3mk\right){\Im }^{″}=0,$$

(21)

where ${\Im }^{\left(4\right)}=\frac{d^4\Im }{d{\chi }^4}$, ${\Im }^{″}=\frac{d^2\Im }{d{\chi }^2}$, ${\Im }^{\prime }=\frac{d\Im }{d\chi }$. Integrating Equation (21) with respect to $\chi$ once and setting the integral constant to zero, we derive:

$$m^{3}n{\Im }^{‴}+3m^{2}n{\left({\Im }^{\prime }\right)}^2+\left(2ns-3mk\right){\Im }^{\prime }=0.$$

(22)

Based on the sub-equation method, the solution of Equation (22) can be assumed as:

$$\Im \left(\chi \right)={\displaystyle \sum _{i=0}^c{\epsilon }_i{\aleph }^i}\left(\chi \right).$$

(23)

where ${\epsilon }_i\left(i=0,1,2,\dots ,c.\right)$ are constants that can be determined later. Additionally, there is:

$${\aleph }^{\prime }\left(\chi \right)=\sigma +{\aleph }^2\left(\chi \right).$$

(24)

Here, $\sigma$ is a constant. Equation (24) has the following different solutions:

$$\aleph \left(\chi \right)=\{\begin{array}{l}-\sqrt{-\sigma }\tan \mathrm{h}\left(\sqrt{-\sigma }\chi \right),\begin{array}{cc}& \sigma <0\end{array}\\ -\sqrt{-\sigma }\cot \mathrm{h}\left(\sqrt{-\sigma }\chi \right),\begin{array}{cc}& \sigma <0\end{array}\\ \sqrt{\sigma }\tan \left(\sqrt{\sigma }\chi \right),\begin{array}{cc}& \sigma >0\end{array}\\ -\sqrt{\sigma }\cot \left(\sqrt{\sigma }\chi \right),\begin{array}{cc}& \sigma >0\end{array}\\ -\frac{1}{\zeta +\mathrm{\Lambda }},\begin{array}{cccc}\mathrm{\Lambda }& \mathrm{is}& \mathrm{a}& \mathrm{constant},\begin{array}{cc}& \sigma =0\end{array}\end{array}\end{array}.$$

(25)

We can determine the value of c in Equation (23) via balancing ${\Im }^{‴}$ and ${\left({\Im }^{\prime }\right)}^2$ in Equation (22) as:

$$c=1.$$

(26)

Then, Equation (23) becomes:

$$\Im \left(\chi \right)={\epsilon }_0+{\epsilon }_1\Im \left(\chi \right).$$

(27)

Substituting Equation (27) with Equation (24) into Equation (22) and setting their coefficients of the different powers of $\Im \left(\chi \right)$ to zero, it yields:

Solving them, we derive:

$$\begin{array}{l}{\Im }^0\left(\chi \right):-3km\sigma {\epsilon }_1+2ns\sigma {\epsilon }_1+2m^{3}n{\sigma }^2{\epsilon }_1+3m^{2}n{\sigma }^2{\epsilon }_1^2=0,\\ {\Im }^2\left(\chi \right):-3km{\epsilon }_1-ns{\epsilon }_1+8m^{3}n\sigma {\epsilon }_1+6m^{2}n\sigma {\epsilon }_1^2=0,\\ {\Im }^4\left(\chi \right): 6m^{3}n{\epsilon }_1+3m^{2}n{\epsilon }_1^2=0.\end{array}$$

Case 1:

$${\epsilon }_0={\epsilon }_0, {\epsilon }_1={\epsilon }_1, m=-\frac{{\epsilon }_1}{2}, n=n, k=k, s=-\frac{{\epsilon }_1\left(3k+n\sigma {\epsilon }_1^2\right)}{4n}, \sigma =\sigma .$$

Thus, the TWS of Equation (1) can be obtained as:

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0-{\epsilon }_1\sqrt{-\sigma }\tan \mathrm{h}\left[\sqrt{-\sigma }\left(-\frac{{\epsilon }_1}{2}x+ny+kz-\frac{{\epsilon }_1\left(3k+n\sigma {\epsilon }_1^2\right)}{4n}t\right)\right],\sigma <0.$$

(28)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0-{\epsilon }_1\sqrt{-\sigma }\cot \mathrm{h}\left[\sqrt{-\sigma }\left(-\frac{{\epsilon }_1}{2}x+ny+kz-\frac{{\epsilon }_1\left(3k+n\sigma {\epsilon }_1^2\right)}{4n}t\right)\right],\sigma <0.$$

(29)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0+{\epsilon }_1\sqrt{\sigma }\tan \left[\sqrt{\sigma }\left(-\frac{{\epsilon }_1}{2}x+ny+kz-\frac{{\epsilon }_1\left(3k+n\sigma {\epsilon }_1^2\right)}{4n}t\right)\right],\sigma >0.$$

(30)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0-{\epsilon }_1\sqrt{\sigma }\cot \left[\sqrt{\sigma }\left(-\frac{{\alpha }_1}{2}x+ny+kz-\frac{{\epsilon }_1\left(3k+n\sigma {\epsilon }_1^2\right)}{4n}t\right)\right], \sigma >0.$$

(31)

Case 2:

$${\epsilon }_0={\epsilon }_0, {\epsilon }_1=-2m, m=m, n=n, k=-\frac{2\left(2m^{3}n\sigma -ns\right)}{3m}, s=s, \sigma =\sigma .$$

Thus, the TWS of Equation (1) can be obtained as:

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0+2m\sqrt{-\sigma }\tan \mathrm{h}\left[\sqrt{-\sigma }\left(mx+ny-\frac{2\left(2m^{3}n\sigma -ns\right)}{3m}z+st\right)\right],\sigma <0.$$

(32)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\alpha }_0+2m\sqrt{-\sigma }\cot \mathrm{h}\left[\sqrt{-\sigma }\left(mx+ny-\frac{2\left(2m^{3}n\sigma -ns\right)}{3m}z+st\right)\right],\sigma <0.$$

(33)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0-2m\sqrt{\sigma }\tan \left[\sqrt{\sigma }\left(mx+ny-\frac{2\left(2m^{3}n\sigma -ns\right)}{3m}z+st\right)\right],\sigma >0.$$

(34)

$$\mathrm{\Pi }\left(x,y,z,t\right)={\epsilon }_0+2m\sqrt{\sigma }\cot \left[\sqrt{\sigma }\left(mx+ny-\frac{2\left(2m^{3}n\sigma -ns\right)}{3m}z+st\right)\right], \sigma >0.$$

(35)

4. The Physical Interpretations

The obtained solutions will be presented by the 3D plot and 2D contour in this section by taking the reasonable parameters.

By assigning the parameters as $k_2=1$, $k_9=2$, $u_1=1$, $u_2=1$, $u_3=1$, the multi-wave complexiton solution given by Equation (6) for the different time is illustrated in Figure 1 in the form of the 3D plot and 2D contour. Obviously, we can find there is a collision phenomenon between the singular periodic wave and the lump in the outline. As $t$ increases, the waveform propagates in the negative direction of the x axis and y axis.

We illustrate the dynamic behavior of Equation (12) by selecting $k_1=1$, $k_2=1$,$u_1=0.6$, $u_2=0.4$ in Figure 2. From this, collision phenomena between the breather waves and singular periodic waves are revealed. We can observe that the waveform travels along the negative direction of the x axis and positive direction of y axis.

Selecting $k_1=1$, $k_2=2$, $u_1=1$, $u_2=1$, $u_3=1$, $u_4=1$, we present the performance of Equation (15) in Figure 3. Here, it can be found the waveform propagates along the negative direction of the x axis and positive direction of y axis. Additionally, the outline of the wave can be explained as the interaction between lump solution and trigonometric function solution.

By using the parameters as ${\epsilon }_0=1$, ${\epsilon }_1=1$, $n=1$, $k=1$, $\sigma =-1$, the dynamic characteristics of Equation (28) are revealed in Figure 4, where Figure 4a is the 3D plot, Figure 4b is the 2D contour and Figure 4c represents the 2D curve. In our observation, it is a dark wave. With the same parameters, Figure 5 illustrates the behaviors of Equations (3) and (10), which is a bright-dark wave.

The performances of Equations (29) and (30) are presented in Figure 6 and Figure 7, respectively with ${\epsilon }_0=1$, ${\epsilon }_1=1$, $n=1$, $k=1$, $\sigma =1$. We find that the profiles are both singular periodic waves.

5. Conclusions and Future Recommendation

In this article, we obtained multi-wave complexiton solutions, multi-wave solutions and periodic lump solutions of the (3+1)-dimensional Jimbo-Miwa equation with the help of the Hirota bilinear method. Besides, we also construct its diverse travelling wave solutions like the dark, bright-dark and singular periodic wave solutions by applying the sub-equation method. The evolution phenomenon of these different solutions are described graphically. From these descriptions, the physical behavior and the interaction are presented. The obtained results in this work are all new and have not been reported elsewhere. Additionally, they show that the methods adopted are effective and direct, and can moreover be used to study the other PDEs arising in physics.

In recent years, the interest in fractal and fractional calculus [41,42,43,44,45,46,47,48,49] has intensified in different fields due to their strong ability to describe complex phenomena. Applying the fractal and fractional calculus to Equation (1) and obtaining the exact solutions will animate our future research.