Some new exact solutions of (3+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(3+1)$\end{document}-dimensional Burgers system via Lie symmetry analysis

In this paper, by using the Lie symmetry analysis, all of the geometric vector fields of the (3+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(3+1)$\end{document}-Burgers system are obtained. We find the 1, 2, and 3-dimensional optimal system of the Burger system and then by applying the 3-dimensional optimal system reduce the order of the system. Also the nonclassical symmetries of the (3+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(3+1)$\end{document}-Burgers system will be found by employing nonclassical methods. Finally, the ansatz solutions of BS equations with the aid of the tanh method has been presented. The achieved solutions are investigated through two- and three-dimensional plots for different values of parameters. The analytical simulations are presented to ensure the efficiency of the considered technique. The behavior of the obtained results for multiple cases of symmetries is captured in the present framework. The outcomes of the present investigation show that the considered scheme is efficient and powerful to solve nonlinear differential equations that arise in the sciences and technology.


Introduction
The Burgers system describes the propagation processes for nonlinear waves in fluid mechanics such as diverse non-equilibrium, nonlinear phenomena in turbulence, and interface dynamics [12]. Also, this system is used in solitary wave theory to expand integrable models with the extending of famous physical equations. The physicists and mathematicians are looking to study the (1 + 1)-dimensional and the (2 + 1)-dimensional integrable models. Because of the real physical space-time being (3 + 1)-dimensional, researchers have been attempting to detect higher-dimensional models in some ways [29,31,32,51]. Researchers have used many effective techniques to discover the different solution of the Burgers system [10,48]. One of the Sophus Lie's significant discoveries in differential equation is to indicate that transforming nonlinear conditions is possible by through infinitesimal invariants which can correspond to the generators of the symmetry group of the system [28]. Having the symmetry group of a system of differential equations has many advantages including the ability to classify the solutions of the differential equations system. This classification is such that we consider both answers in a category that can be con-verted by some transformational equation. Researchers interested in further reading may refer to recently published articles, including [4,16,24,49]. Another use of these groups is that the differential equations can be classified according to the desired parameter or function. Ovsiannikov [43] gave the strategy for halfway invariant solutions. His technique depends on the idea of an equivalence group, which is a Lie symmetry group working in the extended space of free variables, functions and their subordinates, and saving the class of incomplete differential conditions. The investigation of the correct arrangements has a noteworthy influence in the perusing of nonlinear physical frameworks [37]. Probably the main strategies are the inverse scattering method [18], the Hirota bilinear method [19,27], Lie symmetry analysis [11,40], Darboux and Bäcklund transformations [27], the tanh-function method (Duffy and Parkes [15]). In [13], Chao-Qieg and Yue have obtained a new analytical solution of the (3 + 1)-dimensional Burgers system relied on the Riccati equation. And In [12] another mapping equation is utilized to find variable separation solutions of the (3 + 1)-dimensional Burgers system and three class of variable separation solutions are driven. Also in [22,23], Ibragimovs new technique for finding conservation laws and the idea of nonlinear self-adjointness were explained that widely implemented to find conservation laws of equations (for example see [1,5]). Especially In [3], the nonlinear self-adjointness and conservation laws of the (3 + 1)-dimensional Burgers equation has been obtained by the aid of Ibragimovs method, This equation describes the propagation processes for nonlinear waves in fluid mechanics such as diverse non-equilibrium, nonlinear phenomena in turbulence and inter-face dynamics. Lie symmetries of BS will present some additional results which are also obtained from the Lie algebra structure of Lie symmetry group. We, therefore, plan to make an optimal system of 1-subalgebras of the BS which is useful for classifying of group invariant solutions and to construct an optimal system of 2-subalgebras for the BS, also for 3-subalgebra of the BS which are effectively helping to classify the group invariant solutions. After that, the reduced equation for each element of the optimal system is obtained. Then, using the method of Lie symmetry group a solution will be presented for reduced equations. Afterward, the group invariant solutions of BS is achieved. In [30], the Lie point symmetries of (3 + 1)-dimensional Burgers system is obtained and the system is reduced by them. Finally, we found the nonclassical symmetries that were first discussed in Bluman and Cole [8] in their treatment of generalized self-similar solutions of the linear heat equation [42]. This method and its correlation to direct reduction methods of Clarkson [14] and Galaktionov [17] have become the focus of much research and many applications to physically significant partial differential equations. Obviously, other related points, such as partially invariant solutions differentially partially invariant solutions of group foliation, will give rise to effective and delicate methods of treating differential equations. In [33,50], they have been able to find a powerful solution method, By using the tanh method for the computation of traveling wave solutions. First an ansatz of power series in tanh was used to obtain solutions of traveling wave of certain nonlinear evolution equations [34].
Recently, Fan and Hon [16] introduced a generalized tanh method to investigate special types of nonlinear equations.
Hereafter, the symmetry reductions, explicit solutions, convergence analysis and conservation laws to the Chen-Lee-Liu model in nonlinear optics were been studied by the authors in [4]. Inc and co-authors surveyed the time fractional generalized Burgers-Huxley equation with Riemann-Liouville derivative via Lie symmetry analysis and the power series expansion method [24]. Tchier and et al. presented an investigation and analysis for the space-time Carleman equation (STCE) in nonlinear dynamical system [49]. Also, some researchers utilized the Lie symmetry analysis for various nonlinear PD equations. See references [2][3][4][5]. Some authors studied the stochastic influenza model with constant vaccination strategy [6], the stochastic meme epidemic model with investigate effect of threshold number [47], a reliable numerical analysis of a stochastic HIV/AIDS model in a two-sex population considering counseling and anti-retroviral therapy [46]. Some research as regards nonlinear study of equations has been made covering nonlinear vibrations of Euler-Bernoulli beams [44], parametrically excited nonlinear oscillators [45], and nonlinear free vibration analysis of tapered beams [7]. This article is organized as follows. In Sect. 2, the Lie symmetries of BS equations have been obtained. By the aid of obtained symmetries, the invariant groups of BS equations are given in Sect. 3. In Sects. 4, 5, and 6 the optimal systems of order one, two, and three of BS equations are presented, respectively. By considering the new coordinates the BS equations are reduced in Sect. 7. In Sect. 8, the nonclassical symmetries of the BS equation are driven. In Sect. 9, by implementing of tanh-function method ansatz solution of BS equations has been found and plotted. Also, in Sect. 10 the results and discussion of graphs and their behavior are given.

Lie symmetry of BS
The method which is conventionally applied to determine the classical symmetries of a partial differential equation is standard and is explained in [9,40,41] we consider an infinitesimal generator of Lie symmetry from the following form in order to get the Lie symmetry: where ξ t , ξ x , ξ y , ξ z are functions of t, x, y, z, and η u , η v , η w are functions of t, x, y, z, u, v, w. By using the invariant condition, such as applying pr (2) X the second prolongation of X to BS, we have the following system of equations: Therefore, the Lie symmetry group of BS includes a Lie algebra generator in the form of the vector field υ having these functional coefficients: where c i , i = 1, . . . , 5 are arbitrary constants and F i (t), i = 1, 2, 3 are arbitrary smooth functions. In the case the above system is solved, the following theorem is introduced.

Theorem 2.1 The Lie algebra L(G) of projectable Lie symmetries of BS is spanned by
where f (t), g(t), h(t) are arbitrary smooth functions.
Having functional coefficients, these vector fields produce a Lie pseudo-group L(G). This Lie pseudo-algebra L(G) has a 5-subalgebra h generated by v 1 , . . . , v 5 , and an infinite dimensional ideal i generated by A f , B g , C h . Therefore L(G) h × i. To having a reduction in BS, a Lie subgroup of the above pseudo-group is chosen. For executing this chose, simpler forms for every one of the coefficients in the previously mentioned vector fields are chosen.

Group invariant solutions of BS
In order to have the group of transformations which are generated by vector fields v i , i = 1, . . . , 5, we need to at first solve first order system including the first order equation which is in agreement whit each of the same generators. If k (s) is the parametric group represented by v k , k = 1, . . . , 5, then P k = k (s)(t, x, y, z, u, v, w) are, respectively, and if h , g , f are the 1-parameter group generated by A h , B g , C f , respectively, then It should be mentioned that, generally, for each, a parameter of the subgroup of complete symmetry group's system offers a set of invariant solutions [35,40,41].
x, y, z)) be a solution of the BS and i = 1, . . . , 5, there are the functions Furthermore, if h , g , f are the 1-parameter group generated by A f , B g , c h , respectively, and u = u(t, x, y, z), v = v(t, x, y, z) and w = w(t, x, y, z) are a solution of the BS, so are the This theorem has several useful consequences. For example, by using the command PDETools[TWSolutions] of MAPLE 2016© for the system BS, we find a seven parametrized set of tanh-solutions for BS as where c 1 , . . . , c 7 ∈ IR are arbitrary numbers with c 4 = 0. Now, we find a large set of solutions of BS by using Theorem 3.1. For example, by the ψ h of Theorem 3.1 for (1), we find where c 1 , . . . , c 7 ∈ R are arbitrary numbers with c 4 = 0, and h is an arbitrary smooth function of t.

Optimal system of 1-subalgebras
In this part, we take the advantage of symmetry group in order to obtain the OS (optimal system) of 1-subalgebras of BS. Regarding the fact that every linear combination of symmetries is a symmetry, there will be an endless number of 1-subgroups for G. Therefore, determining the subgroups which give different types of solutions is emphasized. So, we need to look symmetry transformations for invariant solutions which are unable to be to convert to each other in the full symmetry groups. This, in turn, leads to the notional of an OS of subalgebra. The problem of classifying this 1-subalgebra is identical to the problem which exists in the classification of the orbits of the adjoint representation [36,40]. The optimal set of subalgebras is achieved by selecting just one representative from any class of equivalent subalgebras. By using a general member in the Lie algebra and simplifying it via diverse adjoint transformations, it is possible to solve the problem of orbits classifications. The following Lie series are considered to include the adjoint representation [38]: where s is the group parameter and ad X i X j = [X i , X j ] is the Lie algebra communicator and i, j = 1, . . . , 9 An adjoint action is considered for Lie algebra L(G) so we consider the following theorem, We can simplify X as follows:

Theorem 4.1 The OS of 1-subalgebras for BS is
• If a 5 = a 3 = a 2 = a 1 = 0, then X is decreased to the case (1).
• If a 5 = a 2 = a 1 = 0 and a 3 = 0, so we can get the coefficient of X 4 vanish by using F 2a 4 /a 3 4 . In this case is reduced to the case (9). • If a 5 = a 3 = a 2 = 0 and a 1 = 0, so we can get the coefficient of X 4 vanish by using F 2a 4 /a 1 5 . In this case is reduced to the case (8). • If a 5 = a 4 = a 3 = 0 and a 2 = 1, then coefficient of X 1 can be vanished or be ±1 by F s 3 3 for s 3 = -ln |1/a 1 |. So X is reduced to the case (2, 6). • If a 5 = a 3 = a 1 = 0, then we can make the coefficient of X 4 vanish or be ±1 by using F s 3 3 for s 3 = (2/3) ln |a 2 /a 4 |. In this case is reduced to the case (3, 7). • If a 1 = a 3 = a 2 = 0 and a 5 = 0, then we can assume that a 5 = 1 so we can get the coefficient of X 4 vanish by using F 2a 4 1 . In this case it is reduced to the case (5). • If a 5 = a 2 = a 1 = 0, then we can suppose that a 1 = a 5 = 1, so we can get the coefficient of X 4 and X 3 vanish by using F 2a 4 1 and F -a 3 2 . In this case is to the case (4). • If a 1 = a 3 = 0 and a 5 = a 2 = 0, then we can make the coefficient of X 4 vanish by using F 2a 4 1 , and X 5 vanish by using F 1/a 2 4 . In this case is decreased to the case (10).

OS of 2-subalgebras
Here, we get OS of 2-subalgebras for BS, choose X 1 or X 2 , as an element of the OS of 2subalgebras in Theorem 3.1 and regard X = c 1 v 1 +· · ·+c 5 v 5 as an optimal vector field where c 5 i s are smooth functions of (t, x, y, z, u, v, w). By such selection, we have [v i , X] = λX 1 + MX.
The following system (2) is a computation of both sides of these equations: The elements of the OS of 2-subalgebras are reached by solving the system of linear equations for each choice of the OS of 1-subalgebras from the Theorem 3.1. When these elements are in the combination form, they could be simplified as a 1-dim case by acting the adjoint matrices to each of them and the following theorem is suggested.

Theorem 5.1 An OS of 2-subalgebras from the BS is
X 3 , X 1 + X 5 , (7) X 3 , X 4 + X 2 , (8) X 1 , X 3 + X 4 , (9) X 2 , X 3 + X 5 , (10) X 3 , X 4 -X 2 , (11) X 1 , X 3 -X 4 , (12) X 3 , X 5 , (13) X 4 , X 2 + X 5 , (14) X 3 , X 1 + X 4 , (15) X 1 , X 4 , Proof Each 2-subalgebra needs two generators. By choosing one of the generators of the OS of 1-subalgebras that have explained in the previous theorem and the second one arbitrary, 2-subalgebras can be classified. let h = Span{X, Y } is a 2-subalgebra where X is a 1subalgebra that is selected from the 1-subalgebra list and Y is an arbitrary vector described by Y = b 1 v 1 + · · · b 5 v 5 . Now, we need to simplify h as much as possible by implementing various adjoint transformations on it and proceed algebraically [38]. Every adjoint transformation is a linear map Fl s i : L(G) → L(G) defined by X − → Ad(exp(sv i )X) for i = 1, . . . , 5, we only illustrate one of the cases in the following: So, we have: for any r, s ∈ RR. Since h, is not closed under the Lie bracket, So there are no two-dimensional subalgebras in this case. • If b 5 = 0 and b 4 = 0 so we can get the coefficient of X 2 vanish by F s 2 2 ; By setting s 3 = (-2/3) ln |b 4 /b 2 |. Then we have h = X 1 , b 3 X 3 + b 4 X 4 and [X 1 , b 3 X 3 + b 4 X 4 ] = -b 3 X 1 . Therefore h is closed under the Lie bracket and we have the case X 1 , So h is closed under the Lie bracket and we have the case X 1 , X 2 + X 3 . • If b 3 = 0 and b 5 = 0 and b 5 = 1 so we can get the coefficient of X 4 vanish by F s 5 5 ; By setting s 5 = -2a 4 also, we have scaling if necessary, we can assume that b 5 = 1. Then we have h = X 1 , b 2 X 2 + X 5 and [X 1 , b 2 X 2 + X 5 ] = b 2 X 3 + b 5 X 4 . So, h is not closed under the Lie bracket and we have no two-dimensional subalgebras in this case. It should be noted that there is no other way to check, and in each case the h is simplified as much as possible. In the same way, two-dimensional algebras can be obtained for other states.

OS of 3-subalgebras
In this section, we get OS of 3-subalgebras for BS. To find 3-subalgebras, we must consider one 3-subalgebra as h = X, Y , Z of symmetry group, so that Z = 5 i=1 c i v i .

Similarity reduction of BS
If expressed in the new coordinates, BS is reduced. The BS is presented in the coordinates of (t, x, y, z, u, v, w) and we have to look for this equation's form in the appropriate coordinates in order to make it reduced. These new coordinates will come in hand by searching independent invariants (p, ϕ, ψ, η) which correspond to the generators of the symmetry group. Hence, If we use the new coordinates and apply chain derivative role into account, we can have the reduced equation. We introduce the present procedure for one of the infinitesimal generators in the OS of Theorem 6.1 and provide the list of results for some other instances. This equation has an independent variable, named p, and three dependent variables labeled as (ϕ, ψ, η). Similarly, it is possible to compute all of the identically reduced equation which correspond to the infinitesimal symmetries that were mentioned in Theorem 6.1. Some of them are listed in Tables 1 and 2. Then we can reduce some equations obtained in the previous section to ODEs, ϕ + p 2 ϕ + 2pϕ (ηψ + 2) + 2ϕ(1ψ) = 0, pϕ + ϕ = 0, ϕ = 0.
Therefore, ϕ = 0, ψ = ψ(z/x), η = η(z/x). We have a non-trivial solution of equation BS: By applying Theorem 3.1 we are able to obtain a new solution of the equation as follows.
are solutions of the equation BS.

Nonclassical symmetries of BS
A kth order system of differential equations is naturally treated as a submanifold E ⊂ J k of the kth order jet space on the space of independent and dependent variables. Consider a kth order system E of differential equations [21]: involving x = (x 1 , . . . , x n ) and u = (u 1 , . . . , u q ) as independent and dependent variables, respectively, and with u (k) denoting the derivatives of the u , s with respect to the x , s up to order k. Suppose that V is a vector field on the space R n × R q of independent and dependent variables: A n-dimensional submanifold of RR n × RR q of the space of independent and dependent variables is defined by a map of the solution u α = f α (x 1 , . . . , x n ) = 0, α = 1, . . . , q, to the system. Then we must add the invariant surface conditions of these equations: where Q = (Q 1 , . . . , Q n ) is known as the characteristic of the vector field (5). The kth prolongation of the invariant surface condition (6) will be denoted by E k Q . If the kth prolongation of the V (k) vector field V is tangent to the intersection E ∩ E k Q then the system (3) and (6) are compatible: The vector field (4) is called a nonclassical infinitesimal symmetry of the system (3) If Eqs. (7) are satisfied [30]. For finding the nonclassical symmetries, according to the system E of BS we do as follows.
Case 3: Again we suppose that the coefficient of ∂ x in (4) is equal to zero and try to obtain the infinitesimal nonclassical symmetries of the form V = ∂ y + ζ ∂ z + φ∂ x + ψ∂ v + θ∂ w . So we obtain the solutions φ = -1/2 and ψ = θ = ζ = 0.

Case 4:
In the end, we assume that the coefficient of ∂ y in (4) is equal to zero and try to obtain the infinitesimal nonclassical symmetries of the form V = ∂ z + φ∂ x + ψ∂ v + θ∂ w . Then we have φ = ψ = 0, and θ = F(t, x, z, w), where F is an arbitrary function.

Conflict of interest
The authors declare that they have no conflict of interest.