Conservation laws and exact solutions of the (3+1)$(3+1)$-dimensional Jimbo–Miwa equation

*Correspondence: j_manafianheris@tabrizu.ac.ir 1Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Tabriz, Tabriz, Iran 2Natural Sciences Faculty, Lankaran State University, 50, H. Aslanov str., Lankaran, Azerbaijan Full list of author information is available at the end of the article Abstract In this article, by using the Herman–Pole technique the conservation laws of the (3 + 1)– Jimbo–Miwa equation are obtained, and then by using the Lie symmetry analysis all of the geometric vector fields of this equation are given. Also, the non-classical symmetries of the Jimbo–Miwa equation have been determined by applying nonclassical schemes. Eventually, the ansatz solutions of the Jimbo–Miwa equations utilizing the tanh technique have been offered.


Introduction
Differential equations play a significant and key role in all sciences and disciplines; by using the analysis of these differential equations, the physical behaviors and manner of interaction and communication with the surrounding world can be discovered. So far, researchers have made great efforts to solve differential equations, including linear and nonlinear equations, and so on. One of the methods of achieving the solution is the conjoint analysis of the differential equations [14]. A study on Lie groups in the late nineteenth century was conducted by Sophus Lie [21], the Norwegian mathematician, inspired by the methodology of variate Galois, the French mathematician, an algebraist who had studies and articles on group theory, and used it for the analysis of polynomial equations. By using symmetries and applying them to the group's work field, Lie was very interested in simplifying and eliminating the ambiguities of partial differential equations (PDEs) and was able to make a great revolution in science. Thus, Lie's group analysis method is considered as one of the systematic methods employed to achieve the nonlinear differential equations' solutions and plays an important role in this regard [26]. This method not only takes a big step toward obtaining differential equations by providing the appropriate tools, but also has approved its applicability by linking concepts such as conservation law and Lie's symmetries in physics, mechanics, and other sciences [24,25].
The conservation law in physics is the constancy of a physical quantity associated with a specific system during the evolution of it. Conservation laws of a differential equation have a very effective role in analyzing the properties of differential equations, more specifically finding exact solutions and expanding numerical methods to obtain more accurate answers and also finding nonlocal dependent devices [15,16]. One of the most popular methods for obtaining conservation laws is Noether's theorem method that expresses the relation between symmetries of a differential equation and its conservation laws by using Euler-Lagrange operator [13].
Although it is a very important result, there are limitations to the existence/nonexistence of any desired variations. Therefore, many efforts have been recently made to resolve this problem, and numerous other methods have been proposed to calculate the conservation law. Regardless of finding a variational problem, Ibragimov [15,16] presented a method which is a generalization of the more fundamental theorem and obtained the conservation law formally based on the symmetries of the given system.
The conservation laws of differential equations can be obtained from different methods, including the classic method of Lie groups using Lie symmetries, in which one or more conservative laws can be constructed using asymmetry obtained from a system of differential equations. To achieve this, different techniques such as the direct technique, the Noether technique, the Boyer technique, the Ibragimov technique, the variational symmetries techniques, and the Herman-Pole technique can be used. The Hermann-Pole method provided in 2010 can largely remove the weaknesses of the direct method and calculate the density and flow of the conservation law in more systematic relations [13].
The JM equation was first introduced by Jimbo and Miwa [18] and then examined by a lot of scholars in various fields for obtaining its solutions [24,25], integrability properties [11], symmetries [27], and so on. Scientists used the Jimbo-Miwa [18] equation to explain certain interesting (3 + 1)-dimensional waves in physics, and it is the second equation in the well-known Painlev'e hierarchy of integrable systems. Most recently, many researchers have investigated this equation. In the [28] Wazwaz obtained multiple-soliton solutions of the JM equation and its extended version.
Tang in [28] received its Pfaffian solution and extended Pfaffian solutions with applying the Hirota bilinear frame. Tang and Liang in [27] gained two kinds of variable separation solutions and plentiful nonlinear coherent structures with the aid of a multi-linear variable separation approach. By using the generalized Riccati equation mapping method, researchers presented rational solutions of (3 + 1)-JM equation [20]. In this paper, we utilize the conservation law and the symmetry method to get the solution of the (3 + 1)dimensional partial differential equation in the following form: which is a so-called Jimbo-Miwa equation in its potential frame. Bibi and co-workers [8] used the generalized Riccati equation mapping method for the space-time conformable Caudrey-Dodd-Gibbon equation by implementing the conformable derivative. Atangana and Araz [5] studied the spread of COVID-19 cases in Turkey and South Africa with an exhaustive statistical analysis. Alkahtani et al. [3] utilized the different operators of fractional differentiation with power law, exponential decay law, and Mittag-Leffler law to the Klein-Gordon equation with mass parameter and also investigated stability and the convergence of the used numerical scheme. Authors of [6] used the Crank-Nicholson scheme and investigated the stability and the convergence to the space fractional variable-order Schrödinger equation. Halder and co-authors [12] worked on Lie symmetry analysis and similarity solutions for the Jimbo-Miwa equation and generalizations. Also, Chauhan and co-authors [10] used the Lie group theoretic method to the similarity reduction and solitary wave solutions of (2 + 1)-dimensional Date-Jimbo-Kashiwara-Miwa equation. Khalique and Moleleki [19] studied a generalized first extended (3 + 1)-dimensional Jimbo-Miwa equation by using the (G /G)-expansion method. In valuable works of scholars the conservation laws for a generalized Ito-type coupled KdV system [23], a new fourth-order integrable nonlinear equation [7], and a (3 + 1)-dimensional B-type Kadomtsev-Petviashvili equation [1], and fractional order coupled KdV system [17] were constructed by increasing the order of partial differential equations. In what follows, we show the Lie symmetry algebra for a JM equation by utilizing the Lie group analysis and find non-classical symmetries of the equation.

Conservation laws for Jimbo-Miwa equation
There are various methods to obtain conservation laws. Herein, we employ the multiplier method. Taking the equation system via n-independent variables x = (x 1 , . . . , x n ) and m-dependent variables. The conservation law for a system of PDEs (2) is a divergence expression that is defined as follows: It is valid for all of the solutions of (3) that here ℵ i [u] = ℵ i (x, u, ∂u, . . . , ∂ k u) is named the fluxes of the conservation law and the order of conservation law is considered. The maximum degree of derivation in the flux expression.
for arbitrary functions u, then a set of multipliers coefficient function The Euler operator for u J is given by exists for ℵ i (x, u, ∂u, . . . , ∂ s-1 u), i = 1, . . . , n, functions.
by arbitrary degree results in the conservation law for system (2) if and only if the following rela-tionship for every arbitrary function u(x) is maintained: The equations obtained from (2) are a set of linear characteristic equations that can be solved by the set of all coefficients σ [u] N σ =1 . Now, we can obtain the functional coefficients of the conservation law for the JM equation.
We can obtain the characteristic equation (3) for the JM equation from the following relation: where u(x, y, z, t) is an arbitrary function. By solving the characteristic equation obtained from (5) are arbitrary coefficients. Then local multipliers are determined in the following frame:

Each of the functional coefficients
The t, y, and z-components can also be defined similar to (6). In (6), I x u j f is given by where Consider I y u j f , I z u j f , and I t u j f same as I x u j f . Now, we determine conserved quantities ψ, φ, , and which conclude of multiplier ε = a (t)x -3a(t)u x . So we consider By using (6) and (7), we have I t u j f = 3a(uu xyu x u y ) + a t u y x.
So, the conservation law of the JM equation for the case ε = a (t)x -3a(t)u x is D x + D y + D z + D t = 0. Now, by using a similar method for other cases, we can find all the conservation laws. Therefore, we obtain in the case ξ = u x : and in the case ξ = b(y, z): and finally for ξ = c(z, t):

Classical symmetries of Jimbo-Miwa equation
The symmetry group of equation (1) is made of the vector field of the following form: The fourth order prolongation X is determined by the following vector field: where ϕ ι = D ι Q + ξ u xι + ηu tι , ϕ ιj = D ı (D j Q) + ξ u xıj + ηu tıj , Q = ϕξ u xηu t is the characteristic of X given by (10) and D i indicates total derivative. Thus, equation (1) admits a Lie point symmetry X if X (4) [u x 3 y + 3u x u xy + 3u y u x 2 + 2u yt -3u xz ] = 0, where X (4) denotes the fourth prolongation of X. We can find the determining equations for the symmetry group of the general JM equation as follows: By solving this system of PDEs, we find that ξ 1 = c 1 x + 3tF 1 v + F 2 t + F 3 z, ξ 2 = -4F 1 z + (3c 3 t -2c 1 )y, ξ 3 = 3c 3 z + c 4 , ξ 4 = 3c 1 t + c 2 , η 1 = 3ytF 1 z -yF 3 z + 4xF 2 t -12xF 1 z -2c 2 + F 4 (z, t), Table 1 The commutator table of g where c i , i = 1, . . . , 4, are arbitrary constants and F i (t), i = 1, . . . , 4, are arbitrary smooth functions. Solving the above determining equations, we reach the following Lie point symmetry generators: where f = f (z), g = g(t), h = h(z), and d = d(z, t) are arbitrary smooth functions. Having functional coefficients, these vector fields produce a Lie pseudo-group g = L(G). This Lie pseudo-algebra g has a 4-sub-algebra L 4 af(1) × af(1) made by X 1 , . . . , X 4 and an infinite dimensional ideal L ∞ generated by X f 5 , X g 6 , X h 7 , X d 8 . Therefore L L 4 L ∞ . The commutation relations satisfied by generators above are shown Table 1.
Consider Ad(exp(s.X).Y ) = ∞ k=0 (-s) k (k!) -1 ad k X Y , where s is the group parameter. An adjoint action is considered for Lie algebra g, so we consider the following theorem.
Theorem 3.1 The optimal system of 1-subalgebras for JM is Proof Let L 4 be the symmetry algebra of Eq. (1) by applying adjoint and X = a 1 X 1 + · · · + a 4 X 4 , that is, it is a nonzero vector field of L 4 . We can simplify as many of the coefficients a 1 , . . . , a 4 as possible through proper adjoint applications on X. So we have the following: • If we assume that a 2 = 0, by scaling X 2 and then assuming a 4 = 0 and scaling X 1 , X is reduced to X 3 , X 2 + X 3 , and X 2 -X 3 . In this case if a 4 = 0 and by scaling X 3 , X is reduced to X 4 , X 2 + X 4 and X 2 -X 4 . • If we assume that a 2 = 0 and a 4 = 0, by scaling X 4 and X 3 , X is reduced to X 4 , X 1 + X 4 , and X 1 -X 4 . If we assume that a 2 = a 4 = 0, X is reduced to a 1 X 1 + a 3 X 3 .
Similarly, we can classify 2-sub-algebra and 3-sub-algebra. For convenience, we can use the normal sub-algebra to make the sub-algebra.
we obtain the following reduced equation: where R := r∂ r + ∂ q . The polynomial solution of this equation (11) is f (w) = c 1 w 2 + c 2 w + c 3 . By a similar argument, for X 1 + X 2 + X 3 , we can obtain the reduced form of the equation as follows: Besides, we could use 3-sub-algebras for reducing the equation. We use 1-sub-algebra to obtain 3-sub-algebras because 3-sub-algebras are vertical to 1-sub-algebras. For instance, with the selection of several algebras of 3-sub-algebras, we reduce the equation. By choosing X 4 , X 2 , X 1 , the equation is reduced to where k = yx 2 , and also selecting X 1 , X 3 , X 4 leads to f = ϕ(t/x 3 ) and the equation becomes zero. By selecting X 1 , X 2 , X 3 , the equation is reduced to the following: 8f 3 k 3 (kf + 9) = 6k 4k 2 f + kf -25k -1 f + 48k 2 ff 2 + 3(5kf -20k -3)f .

Non-classical symmetries of JM
Non-classic symmetry is another way to determine some other solutions for a PDE and ODE system. Here, we use the conventional method to obtain the non-classical symmetry of the JM equation according to the compatibility of the evolutionary equations [4,9].
In a non-classical way, first, we add the conditions of the invariance surface to the equation and then apply the classical symmetry method. So, we express X (4) 1 ≡ 0 mod 1 = 0, 2 = 0, where X is defined in (9) and 1 and 2 are given as 1 := u x 3 y + 3u x u xy + 3u y u x 2 + 2u yt -3u xz and 2 := ηu t + ξ u x + βu y + ζ u zϕ. By using non-classical methods on the JM, we consider the following equations: Therefore we have these solutions: where c 1 and c 2 are free amounts and f = f (t), g = g(z), h = h(z), and d = d(z, t) are arbitrary smooth functions.

Some ansatz exact solutions of JM
Herein, we perform one of the most important ansatz methods (the tanh-function method [22]) to gain exact traveling wave solutions of this nonlinear system of PDEs. For this, a new variable τ = tanh(c 1 x + c 2 y + c 3 z + c 4 t + c 5 ), where c i are arbitrary constants, is considered. By placing this expression in equation (1), we get Then, using the ansatz u = A 0 + A 1 τ + A 2 τ 2 + A 3 τ 3 , where A i are arbitrary constants, we obtain the exact solution by using required simplifications and linear algebra: In Figs. 1-4, the solutions have been plotted for some constant coefficients. By choosing suitable amounts, the dynamical structures of hyperbolic wave solutions are presented in Figs. 1 and 2 including three-dimensional, density, and y-curves plot in Fig. 1 and tcurves plot in Fig. 2. Moreover, by choosing suitable amounts, the dynamical structures  Fig. 3 and t-curves plot in Fig. 4. Similarly, ansatz τ = tan(c 1 x + c 2 y + c 3 z + c 4 t + c 5 ), where c i are arbitrary constants, leads to u = -2c 1 tan c 1 x + c 2 y + c 3 z + c 1 (4c 2 1 c 2 + 3c 3 ) 2c 2 t + c 5 + c 0 .

Conclusions
In this article, we obtained the conservation laws for JM equation that give unequivocal expressions of conserved quantities. Then, we attempted to show the Lie symmetry algebra for JM equation by utilizing the Lie group analysis. Further, we found non-classical symmetries of the equation. Finally, we obtained the exact solitary wave solution of JM equation by utilizing the tanh-function technique. The graphes were plotted containing 3D plot, density plot, and 2D plot. The results are beneficial to the study of the mathematics physics, fluid dynamics, and applied mechanics. All calculations in this paper have been made quickly with the aid of the Maple.