Lie Symmetry Analysis, Explicit Solutions and Conservation Laws of a Spatially Two-Dimensional Burgers–Huxley Equation

In this paper, we investigate a spatially two-dimensional Burgers–Huxley equation that depicts the interaction between convection effects, diffusion transport, reaction gadget, nerve proliferation in neurophysics, as well as motion in liquid crystals. We have used the Lie symmetry method to study the vector fields, optimal systems of first order, symmetry reductions, and exact solutions. Furthermore, using the power series method, a set of series solutions are obtained. Finally, conservation laws are derived using optimal systems.


introduction
Mathematical modeling of dynamical systems results in nonlinear partial differential equations (PDEs). In reality, most complicated phenomena-namely diffusion, reaction, conservation, and many more-can be illustrated by means of partial differential equations. Due to their quintessence, PDEs are studied profusely in science and engineering. Various peculiar methods are designed for obtaining their exact and approximate solutions, which, in turn, help us in quantitative and qualitative analysis of these PDEs. The interested reader can see some of these methods in [1,2].
Lie symmetry analysis is a powerful and influential tool for mathematically analyzing partial differential equations. It can be used in securing analytic solutions or in switching PDEs into solvable ordinary differential equations (ODEs). Diverse symmetry vectors are also discovered for the considered system in some cases, however, sometimes there emerges a chance of linear combination of these vectors. To avoid this, an optimal system is constructed. Each member of this system is used in lessening independent variables of the system until analytic solutions are obtained or PDEs are switched to solvable ODEs [3][4][5][6][7]. The system also analyzes solutions of PDEs of different kinds as well as opens many fields [8]. 2

of 13
The generalized Burgers-Huxley equation, which has many utilizations in the fields of biology, metallurgy, chemistry, mathematics, and engineering is of the following type, This is a non-linear equation that has secured much importance due to its appearance in many physical phenomena and its scientific utilization. The parameters a, b ≥ 0 are real constants, n is a positive integer and k ∈ [0, 1]. When a = 0 and n = 1, Equation (1) reduces to the Huxley equation and with b = 0, n = 1, it becomes the Burgers equation. Some exact numerical and traveling wave solutions to (1) were reported in [9][10][11]. However, the spline collocation method for the Burgers-Huxley equation was discussed in a book by Schiesser [12]. In addition, many other analytical and numerical methods for generalized Burgers-Huxley equations have been developed in the past, see for example [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29].
In this paper, we will analyze the Burgers-Huxley equation in two spatial dimensions, which is of the following form, Equation (2) couples both the assets of the Burgers equation (this is one of the basic models in fluid mechanics and is used to catch some of the properties of turbulent flow in a channel, which occurs due to the interaction of the reverse outcome of convection and diffusion and also describes the format of shock waves, traffic flow, and acoustic transmission [30]) and the Huxley equation (which is used for nerve proliferation in neurophysics and wall proliferation in liquid crystals [31]). So, we will designate this equation as simply a two dimensional Burgers-Huxley equation. The combined Burgers-Huxley equation shows a prototypical imitation that specifies the interaction between the reaction gadget, diffusion transport and convection effects, nerve proliferation, and motion in liquid crystals [32].
In Section 2, the vector fields and optimal systems for (2) are obtained by using the Lie symmetry method. In Section 3, we computed the similarity diminution for one and two dimensional subalgebras and hence obtained the group invariant solutions for (2). Employing the power series method, certain power series solutions are achieved in Section 4. Finally, in Section 5, conservation laws are derived using optimal systems.

Lie Symmetry Analysis
In this section, we will study the Lie symmetries and optimal systems of the Burgers-Huxley equation. Consider the one-parameter Lie group of transformation: where is the group parameter. The infinitesimal operator associated with the above transformations is: The coefficient functions ξ(x, y, t; u), η(x, y, t; u), ζ(x, y, t; u), and ϕ(x, y, t; u) are to be determined and the vector field X satisfies the Lie symmetry condition, where The second prolongation of the infinitesimal generator is given by with where the the operator D i is defined as: Coupling Equations (5) and (6), we can obtain the compatible condition for Equation (5). Substituting extended transformations into the obtained compatible conditions and making the coefficients of several monomials in partial derivatives and numerous powers of u equal, we get the following over determining system of PDEs: By solving the over determining system of PDEs (8), we obtain the coefficient functions ζ, ξ, and η as: where c 1 , c 2 , and c 3 are random constants. The Lie algebra of infinitesimal symmetry of Equation (2) with k = 0 is given by,

Transformed Solutions
One can acquire the group transformation initiated by the Lie point symmetry operator X i (i = 1, 2, 3) by solving the following ODEs The one-parameter Lie symmetry groups generated by infinitesimals X 1 , X 2 , and X 3 are given by where 1 , 2 , and 3 are group parameters. Depending on the values of r 1 , r 2 , and r 3 , if f (x, y, t) is any confidential solutions of Equation (2), then the new solutions can be given by

Optimal System of Subalgebras
For the optimal system, we will first construct the tables for commutation relations and adjoint action of the obtained symmetries. For the sake of obtaining the adjoint representation we will use the Lie series in the form: The commutation relations between basis elements satisfies: whereas, the adjoint representation Table 1 is given as: Table 1. Adjoint representation table. Ad Following the method given in [8], consider the vector X with random coefficients a 1 , a 2 , a 3 , such that suppose a 3 = 0 and set up a 3 = 1, so that The scheme involves simplifying the coefficients as much as possible. To abolish the coefficient of X 2 , we will use it to act on X. It is easy to see that the vector form cannot be reduced much more because commutation relations are zero.
Further, suppose a 3 = 0 and establish a 2 = 1 so that, Repeating the same process and normalizing the coefficients we have the following one-dimensional optimal system of subalgebras, where a and b are arbitrary constants.

Symmetry Reductions for Optimal System
Symmetry diminution and explicit solutions have numerous utilizations in the connection of differential equations. Solutions of partial differential equations asymptotically tend to the solutions of lower dimensional equations, which are obtained by symmetry diminution. Especially, exact solutions originating from symmetry technique can be utilized effectively to study properties such as asymptotics etc. In the previous section, we obtained the point symmetries and optimal system for two-dimensional Burgers-Huxley equation, now we will interrogate the symmetry reductions and exact solutions using one-dimensional subalgebras.
(i) For the linear combination aX 1 + X 2 , the invariants can be obtained by solving the characteristic Using, u = Ψ(s, r), we obtained the following PDE with two independent variables, By the application of the similarity transformation method on reduced Equation (11) again, we have where c 1 and c 2 are random constants and the characteristic Equation for (12) is Hence, Ψ can be written as: Substituting Equation (14) into (11), we have the following ODE: (ii) For the linear combination aX 1 + bX 2 + X 3 , the invariants can be obtained by solving the characteristic equation Using, u = Ψ(s, r), we obtained the following PDE with two independent variables: By the application of the similarity transformation method on reduced Equation (17) again, we obtain where c 1 and c 2 are random constants. Choosing c 1 and c 2 equal to 1, the characteristic Equation for (18) is: Therefore, Ψ can be written as: Substituting Equation (20) into (17), we have following ODE: Reductions corresponding to the remaining vectors occuring in the optimal system can be obtained in a similar way, hence, we omit the details here.

Explicit Power Series Solution
After reducing partial differential equations into ordinary differential equations in the previous section, we will now explore the exact analytic solution of the reduced equations using the power series method [8].

Series Solution of Reduced Equation (15)
Consider the power series solution of the form: Incorporating (22) into Equation (15), we get Inspecting the coefficients for n = 0 and n ≥ 1, we have So, we find the solution of Equation (15) in the form of the power series given by Back substituting the values of α and β, we have Thus, in terms of original variables, we obtain the following solution of the two dimensional Burgers-Huxley equation: where b 0 and b 1 are random constants and the remaining constants can be calculated using Equations (24) and (25), respectively. The convergence of Equation (28) is shown in the Appendix A.

Series Solution of Reduced Equation (21)
Now, consider the power series solution of (21) in the form (22). These equations together yield Comparing the coefficients of like powers of α, we obtain Next, we use values of α and β from (20) into (22) to have where b 0 and b 1 are random constants. Hence, we have the following solution in terms of original variables: with b 2 is given in (29) and b n+2 is given in (30).

Conservation Laws
Conservation laws are of fundamental importance in the study of partial differential equations as they provide conserved quantities for all the solutions, can detect integrability and linearization, and can be used to check the precision of numerical solution methods [33]. Conservation laws provide one of the basic rules in formulating models in mathematics and certain times a partial differential equation having a large number of conservation laws depicts a strong implication of its integrability.

Preliminaries
Here, we recall the Ibragimov scheme of constructing conservation laws corresponding to the given symmetries of any system of PDEs, provided that the number of equations in the system is equal to the number of dependent variables [34]. Consider a p-th order system of partial differential equations of p independent variables and q dependent variables, respectively, as x = (x 1 , x 2 , ..., x p ) and u = (u 1 , u 2 , ..., u q ), G Λ (x, u, u (1) , ..., u (p) ) = 0, Λ = 0, 1, ..., q.
The infinitesimal generator for Equation (33) is: The p-th prolongation of the point symmetry operator is: Here m, n = 1, 2, ..., p, Λ = 1, 2, ..., q, and D n is the total derivative operator. A vector T = (T 1 , ..., T n ) is a conserved vector of (33) if is satisfied for all solutions of (33). Equation (34) is called a local conservation law provided that T i are free of integral terms. with, and the formal Lagrangian is defined as In order to find the adjoint equation we have, where δ δu is the variational derivative and (v Λ , Λ = 0, 1, ..., q) is the adjoint variable. (2) is not strictly self-adjoint.

Theorem 2. Equation
Proof. The formal Lagrangian for Equation (2) can be written as, The adjoint Equation of (2) is given by: where δ δu is the Euler-Lagrange operator given by, Using the Euler-Lagrange operator in Equation (2), we have the following adjoint equation, after transforming in the original variable u, ∂u ∂t which is not strictly self-adjoint.

Conservation Laws of a Spatially Two-Dimensional Burgers-Huxley Equation
Now, we will construct conservation laws for each element in the optimal system obtained in Section 2.
(i) For X 1 , we have W = −u t and ξ t = 1. Substituting these values in Equation (35), we find (ii) For X 2 , W = −u x and ξ x = 1. Substituting in Equation (35), we obtain (iii) For the generator X 3 , W = −u y and ξ y = 1. So, in this case, we have the following conserved vector: (iv) For the generator aX 1 + X 2 , W = −au t − u x and ξ t = a, ξ x = 1. Substituting into (35), we have C y = auvu t − au t v y + uvu x − u x v y + avu yt + vu xy .
(v) For the generator aX 1 + bX 2 + X 3 , W = −au t − bu x − u y and ξ t = a, ξ x = b and ξ y = 1. Using these values in Equation (35), we have The conserved vectors comprise random solutions of the adjoint equation, thereby implying the interminable number of conservation laws. Conservation laws play a compelling role in the solution process of an equation or system of equations.

Concluding Remarks
In this work, we have investigated a spatially two-dimensional Burgers-Huxley equation, which combines both the properties of Burgers (convective phenomenon) and Huxley equation (nerve proliferation and motion in liquid crystals). We have computed the Lie point symmetries of (2), and constructed an optimal system of one and two dimensional subalgebras using these symmetries. Symmetry reductions are performed for one-dimensional subalgebras. Furthermore, some power series solutions are calculated using the power series method. Finally, conservation laws are derived using the Ibragimov theory and involvement of random solutions predicted that conservation laws consisted of an interminable number of conserved vectors. Now define a power series γ(α) = ∑ ∞ n=0 a n α n and consider a i = |b i | for i = 0, 1, 2, .... Then we have a n+2 ≤ M n ∑ m=0 a m a n−m + a n + n ∑ m=0 m ∑ j=0 a j a m−j a n−m n ∑ m=0 a m a n−m+1 − a n+1 .

(A3)
It can be easily seen that |b i | ≤ a i and hence the assumed series is majorant series. Now, we will show that the series γ(α) = ∑ ∞ n=0 a n α n is convergent. Through some calculations, we have γ(α) = a 0 + a 1 α + M Now, taking into account the implicit functional system with respect to α, we have We consider the point (0, a 0 ), as γ is analytic in the neighbourhood of the point (0, a 0 ), as γ(0, a 0 ) = 0 and ∂ ∂γ (0, a 0 ) = 0. Hence, by the implicit function theorem, the series γ(α) = ∑ ∞ n=0 a n α n is analytic in the neighbourhood of the point (0, a 0 ) and has a positive radius of convergence. Hence, convergence is proven.