Welcome to visit Communications in Theoretical Physics,
Mathematical Physics

Consistent Burgers equation expansion method and its applications to high- dimensional Burgers-type equations

  • Gangwei Wang , 1 ,
  • Li Li 2 ,
  • A H Kara 3
Expand
  • 1School of Mathematics and Statistics, Hebei University of Economics and Business, Shijiazhuang 050061, China
  • 2Library, Hebei University of Economics and Business, Shijiazhuang 050061, China
  • 3School of Mathematics, University of the Witwatersrand, Private Bag 3, Wits 2050, Johannesburg, South Africa

Received date: 2022-02-20

  Revised date: 2022-04-13

  Accepted date: 2022-06-08

  Online published: 2022-08-01

Copyright

© 2022 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

In this paper, a novel method, named the consistent Burgers equation expansion (CBEE) method, is proposed to solve nonlinear evolution equations (NLEEs) by the celebrated Burgers equation. NLEEs are said to be CBEE solvable if they are satisfied by the CBEE method. In order to verify the effectiveness of the CBEE method, we take (2+1)-dimensional Burgers equation as an example. From the (1+1)-dimensional Burgers equation, many new explicit solutions of the (2+1)-dimensional Burgers equation are derived. The obtained results illustrate that this method can be effectively extended to other NLEEs.

Cite this article

Gangwei Wang , Li Li , A H Kara . Consistent Burgers equation expansion method and its applications to high- dimensional Burgers-type equations[J]. Communications in Theoretical Physics, 2022 , 74(8) : 085004 . DOI: 10.1088/1572-9494/ac76b6

1. Introduction

Nonlinear evolution equations (NLEEs) play an important role in many fields [19], they come up in various science fields, such as nonlinear optics, fluid mechanics, plasma physics, condensed matter physics, biophysics, etc. It is well known that NLEEs is an important part of applied mathematics and mathematical physics. In recent years, a great many of authors have made many contributions in NLEEs from different aspects, which has promoted their development. Especially in recent decades, with the rapid development of mathematical physics and computer science, the subject has developed more rapidly and achieved more results.
Due to the complexity of NLEEs, in general, a large number of existing important NLEEs, are not easy to find analytical solutions. Even if one can find some exact solutions, it requires some skill. Besides, new solutions with physical significance need to be further constructed and discovered. Through the continuous efforts of mathematicians and physicists, a large number of effective methods have been established and developed, for example, the symmetry method [516], inverse scattering transformation [1], ${\rm{B}}\ddot{a}$cklund transformation [17, 18], Darboux transformations [19], the Hirotas bilinear direct method [20], auxiliary equation expansion method [21], generalized multi-symplectic method [2225] homogeneous balance method [26], F-expansion method [27], homotopy perturbation method [28, 29], CKs direct symmetry reduction method [30], consistent Riccati expansion method [31], consistent KdV expansion method [32] and so on.
Motivated by these papers, we propose a new approach, named as consistent Burgers equation expansion (CBEE) method to solve the NLEEs. On the basis of the (CBEE) method, NLEEs can be said to be CBEE solvable if it satisfies the CBEE.
As mentioned earlier, the exact solution of NLEE plays a very important role in explaining complex nonlinear phenomena. In addition to some methods mentioned in the previous references and the references cited therein, it is necessary to find new methods to solve NLEEs through appropriate skills or transformations. In the present paper, we found such a transformation, which directly converts complex research objects into simple classical research objects. For the high-dimensional Burgers equation, we can directly use this method to get their new solutions. In this way, the calculation is greatly simplified. Because compared with high-dimensional and high-order NLEEs, low-dimensional and low-order ones are easier to study. Therefore, not only the calculation is simplified, but also new solutions are obtained. In a word, the advantage of this method is that as long as we study the (1+1)-dimensional Burgers equation, we can get new solutions of high-dimensional burgers and other NLEEs.
The present paper is divided into the following parts. In section 2, the basic idea and steps of CBEE and some conceptions are given. In section 3, the application of this method to the (2+1)-dimensional Burgers equation is described systematically. Symmetries and conservation laws are presented in section 4. A brief summary and discussion are presented in the last section.

2. Basic definition of CBEE method and CBEE solvability of NLEEs

In the present section, we give some definitions and concepts of the CBBE and the CBEE solvability. Consider the NLEEs as follows
$\begin{eqnarray}\begin{array}{l}{\boldsymbol{Q}}({\boldsymbol{x}},t,{\boldsymbol{u}})=0,\quad {\boldsymbol{x}}=\left\{{x}_{1},{x}_{2},\ldots ,{x}_{p}\right\},\\ \quad {\boldsymbol{Q}}=\left\{{Q}_{1},{Q}_{2},\ldots ,{Q}_{m}\right\}\\ \quad {\boldsymbol{u}}=\left\{{u}_{1},{u}_{2},\ldots ,{u}_{q}\right\},\end{array}\end{eqnarray}$
while Q is polynomial of functions ui and its derivatives term. In order to achieve this goal, we try to expand equation (1) in the following form
$\begin{eqnarray}{u}_{i}=\sum _{m=0}^{{M}_{i}}{A}_{{im}}({\boldsymbol{x}},t){U}^{m}(\xi ,\tau ).\end{eqnarray}$
The positive integers Mi, i = 1, 2, ..., q should be fixed via the leading order analysis from equation (1). In general, the function U(ξ, τ) is a special known function. Here we require that U(ξ, τ) satisfied the following famous Burgers equation [33]
$\begin{eqnarray}{U}_{\tau }+{{UU}}_{\xi }-{U}_{\xi \xi }=0.\end{eqnarray}$
Inserting equation (2) into equation (1) and using equation (3), it should get
$\begin{eqnarray}Q=\sum {g}_{p}\left({A}_{{im}},\xi ,\tau ,\ldots \right)\displaystyle \frac{{\partial }^{k}U}{\partial {\xi }^{m}}.\end{eqnarray}$
Here ${g}_{p}\left({A}_{{im}},\xi ,\tau ,\ldots \right)$ are functions with respect to ξ(t, X), τ(t, X), Aij(t, X) and their different derivatives. To solve these equations, vanish all the coefficients with different powers of $\tfrac{{\partial }^{k}U}{\partial {\xi }^{m}}$, which leads to determining equations
$\begin{eqnarray}{g}_{p}\left({A}_{{im}},\xi ,\tau ,\ldots \right)=0,p=1,2\ldots q.\end{eqnarray}$
In general, as the number of equations is greater than the number of unknown variables, determining equation (5) is overdetermined. Solving them, we should obtain Aij, ξ, τ. Therefore, some new explicit solutions are presented via the Burgers equation. From the above analysis, we give the following statement:

The expansion (2) is a ‘CBEE’ and the NLEEs (1) is ‘CBEE solvable’, if system (5) is consistent.

3. CBEE solvability and novel explicit solutions of the (2+1)-dimensional Burgers equation

Through the analysis in the previous section, in order to prove the effectiveness of our proposed method, we consider the following (2+1)-dimensional Burgers equation [3438]
$\begin{eqnarray}{\left({u}_{t}+{{uu}}_{x}-{u}_{{xx}}\right)}_{x}+{u}_{{yy}}=0.\end{eqnarray}$
This equation describes weakly nonlinear (1+2)-dimensional shocks which appear in dissipative media. With regard to more descriptions of this equation and its applications, see [3438] and references therein.
Based on the steps of CBEE, by the leading order analysis of equation (6), consider the following transformation
$\begin{eqnarray}\begin{array}{l}u(x,y,t)=\alpha (x,y,t)\\ \quad +\beta (x,y,t)U\left(\xi (x,y,t),\tau (x,y,t)\right),\end{array}\end{eqnarray}$
where α(x, y, y), β(x, y, t), ξ(x, y, t), τ(x, y, t) are needed to be fixed later. Here $U\left(\xi ,\tau \right)$ satisfies the Burgers equation (3), that is to say, $U\left(\xi ,\tau \right)$ is the solution of the Burgers equation (3).
Putting equation (7) into (6) and using equation (3), we get some overdetermined partial differential equations (PDEs) with regard to α, β, ξ, τ. Let the coefficients of U and its different derivatives of U be equal to 0, one should obtain
$\begin{eqnarray}\begin{array}{rcl}\alpha & = & \displaystyle \frac{-{a}_{t}}{a}x+\frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t),\beta =a,\\ \tau & = & {\int }^{t}a{\left(\mu \right)}^{2}{\rm{d}}\mu ,\\ \xi & = & {ax}+{dy}+\int \frac{1}{a}\left(\frac{1}{2}(-{a}_{{tt}}a+2{a}_{t}^{2}){y}^{2}\right.\\ & & \left.-{a}^{2}{gy}-{{ha}}^{2}-{d}^{2}-{{ad}}_{t}y\right){\rm{d}}t+{c}_{1},\end{array}\end{eqnarray}$
where a, d, g, h are arbitrary functions of t, c1 is an integral constant.
From the above analysis, in short, the (2+1)-dimensional Burgers equation (6) is CBEE solvable. Therefore, we get the following statement

If $U(\xi ,\tau )$ is a solution of the Burgers equation (3), then

$\begin{eqnarray}u=\displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}+g(t)y+h(t)+{aU}(\xi ,\tau )\end{eqnarray}$
is a solution of the (2+1)-dimensional Burgers equation (6). Here $\xi ,\tau $ are given by equation (8).

Next, from the existing literature, we give some new solutions of the (2+1)-dimensional Burgers equation (6). In paper [38], the authors give some explicit solutions of the (1+1)-dimensional Burgers equation (3)
$\begin{eqnarray}{U}_{1}=\pm 2,\end{eqnarray}$
$\begin{eqnarray}{U}_{2}=-2\tanh (\xi ),\end{eqnarray}$
$\begin{eqnarray}{U}_{3}=-2\coth (\xi ),\end{eqnarray}$
$\begin{eqnarray}{U}_{4}=2\tan (\xi ),\end{eqnarray}$
$\begin{eqnarray}{U}_{5}=-2\cot (\xi ),\end{eqnarray}$
$\begin{eqnarray}{U}_{6}=\displaystyle \frac{-2}{\xi },\end{eqnarray}$
$\begin{eqnarray}{U}_{7}=\displaystyle \frac{-2}{1\pm \exp (-\tau -\xi )},\end{eqnarray}$
$\begin{eqnarray}{U}_{8}=\displaystyle \frac{2}{1\pm \exp (-\tau +\xi )},\end{eqnarray}$
$\begin{eqnarray}{U}_{9}=\displaystyle \frac{-2\sinh \xi }{\cosh \xi \pm \exp (-\tau )},\end{eqnarray}$
$\begin{eqnarray}{U}_{10}=\displaystyle \frac{-2\cosh \xi }{\sinh \xi \pm \exp (-\tau )},\end{eqnarray}$
$\begin{eqnarray}{U}_{11}=\displaystyle \frac{2\sin \xi }{\cos \xi \pm \exp (-\tau )},\end{eqnarray}$
$\begin{eqnarray}{U}_{12}=\displaystyle \frac{-2\cos \xi }{\sin \xi \pm \exp (-\tau )},\end{eqnarray}$
$\begin{eqnarray}{U}_{13}=\displaystyle \frac{-2}{l+\xi },\end{eqnarray}$
$\begin{eqnarray}{U}_{14}=\displaystyle \frac{\xi }{\tau },\end{eqnarray}$
$\begin{eqnarray}{U}_{15}=\displaystyle \frac{\xi /\tau }{1\pm {\tau }^{\tfrac{1}{2}}\exp (\tfrac{{\xi }^{2}}{4\tau })},\end{eqnarray}$
$\begin{eqnarray}{U}_{16}=\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\tanh (\displaystyle \frac{\xi }{\tau }),\end{eqnarray}$
$\begin{eqnarray}{U}_{17}=\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\coth (\displaystyle \frac{\xi }{\tau }),\end{eqnarray}$
$\begin{eqnarray}{U}_{18}=\displaystyle \frac{\xi }{\tau }+\displaystyle \frac{2}{\tau }\tan (\displaystyle \frac{\xi }{\tau }),\end{eqnarray}$
$\begin{eqnarray}{U}_{19}=\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\cot (\displaystyle \frac{\xi }{\tau }).\end{eqnarray}$
Therefore, based on the theorem (1), we get new explicit solutions of the (2+1)-dimensional Burgers equation as follows:
$\begin{eqnarray}\begin{array}{rcl}{u}_{1} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)\pm 2a,\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{2} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)-2a\tanh (\xi ),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{3} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)-2a\coth (\xi ),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{4} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+2a\tan (\xi ),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{5} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)-2a\cot (\xi ),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{6} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a}{\xi },\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{7} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a}{1\pm \exp (-\tau -\xi )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{8} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{2a}{1\pm \exp (-\tau +\xi )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{9} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a\sinh \xi }{\cosh \xi \pm \exp (-\tau )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{10} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a\cosh \xi }{\sinh \xi \pm \exp (-\tau )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{11} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{2a\sin \xi }{\cos \xi \pm \exp (-\tau )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{12} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a\cos \xi }{\sin \xi \pm \exp (-\tau )},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{13} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{-2a}{l+\xi },\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{14} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{a\xi }{\tau },\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{15} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+\displaystyle \frac{a\xi /\tau }{1\pm {\tau }^{\tfrac{1}{2}}\exp (\tfrac{{\xi }^{2}}{4\tau })},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{16} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+a\left(\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\tanh (\displaystyle \frac{\xi }{\tau })\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{17} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+a\left(\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\coth (\displaystyle \frac{\xi }{\tau })\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{18} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+a\left(\displaystyle \frac{\xi }{\tau }+\displaystyle \frac{2}{\tau }\tan (\displaystyle \frac{\xi }{\tau })\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{rcl}{u}_{19} & = & \displaystyle \frac{-{a}_{t}}{a}x+\displaystyle \frac{{a}_{{tt}}a-2{a}_{t}^{2}}{2{a}^{2}}{y}^{2}\\ & & +g(t)y+h(t)+a\left(\displaystyle \frac{\xi }{\tau }-\displaystyle \frac{2}{\tau }\cot (\displaystyle \frac{\xi }{\tau })\right).\end{array}\end{eqnarray}$

4. Symmetries and conservation laws

A nontrivial conservation law of equation (6) exists if there exists a vector (Tt, Tx, Ty) whose divergence
$\begin{eqnarray}{D}_{t}{T}^{t}+{D}_{x}{T}^{x}+{D}_{y}{T}^{y},\end{eqnarray}$
vanishes on the solutions of the PDE equation (6). We adopt the ‘multiplier approach’ [39] to construct the conserved flows where, by multipliers, we mean the differential functions Q(x, y, t, u, ux, …,uxx, …) such that the Euler–Lagrange operator (variational derivative) on
$\begin{eqnarray}Q({\left({u}_{t}+{{uu}}_{x}-{u}_{{xx}}\right)}_{x}+{u}_{{yy}}),\end{eqnarray}$
vanishes; each, such Q leads to a conserved flow. It turns out that for Q up to second-order in derivatives, we obtain only derivative independent, infinitely many forms for Q, viz.,
$\begin{eqnarray}Q={{xf}}_{1}\left(t\right)y+{{xf}}_{2}\left(t\right)-\displaystyle \frac{1}{6}\,{f}_{1}^{\prime} {y}^{3}-\displaystyle \frac{1}{2}\,{f}_{2}^{\prime} {y}^{2}+{f}_{3}\left(t\right)y+{f}_{4}\left(t\right).\end{eqnarray}$
The conserved flow may then be determined by a homotopy formula [40] or by substituting into the definition (48). We state some special cases.
i

(i) Q = xy + y + t:

$\begin{eqnarray}\begin{array}{rcl}{T}^{t} & = & \displaystyle \frac{1}{2}{u}_{x}{xy}+\displaystyle \frac{1}{2}\,{{yu}}_{x}+\displaystyle \frac{1}{2}\,{u}_{x}t-\displaystyle \frac{1}{2}\,{{yu}}_{},\\ {T}^{x} & = & {u}_{x}{u}_{}{xy}+{u}_{x}{{yu}}_{}+{u}_{x}{u}_{}t-\displaystyle \frac{1}{2}\,{{yu}}_{}^{2}\\ & & -{u}_{{xx}}{xy}-{u}_{{xx}}y-{u}_{{xx}}t+\displaystyle \frac{1}{2}\,{u}_{t}{xy}\\ & & +\displaystyle \frac{1}{2}\,{u}_{t}y+\displaystyle \frac{1}{2}\,{u}_{t}t+{{yu}}_{x}-\displaystyle \frac{1}{2}\,{u}_{},\\ {T}^{y} & = & -[-{u}_{y}{xy}-{u}_{y}y-{u}_{y}t+{u}_{}x+{u}_{}],\end{array}\end{eqnarray}$

ii

(ii) $Q=-\tfrac{1}{2}{y}^{2}+{xy}+{xt}$:

$\begin{eqnarray}\begin{array}{rcl}{T}^{t} & = & -\displaystyle \frac{1}{4}\,{u}_{x}{y}^{2}+\displaystyle \frac{1}{2}\,{u}_{x}{xy}+\displaystyle \frac{1}{2}\,{u}_{x}{xt}-\displaystyle \frac{1}{2}\,{{yu}}_{}-\displaystyle \frac{1}{2}\,{u}_{}t,\\ {T}^{x} & = & -\displaystyle \frac{1}{2}\,{u}_{x}{u}_{}{y}^{2}+{u}_{x}{u}_{}{xy}+{u}_{x}{u}_{}{xt}-\displaystyle \frac{1}{2}\,{{yu}}_{}^{2}\\ & & -\displaystyle \frac{1}{2}\,{u}_{}^{2}t+\displaystyle \frac{1}{2}\,{u}_{{xx}}{y}^{2}-{u}_{{xx}}{xy}-{u}_{{xx}}{xt}\\ & & -\displaystyle \frac{1}{4}\,{u}_{t}{y}^{2}+\displaystyle \frac{1}{2}\,{u}_{t}{xy}+\displaystyle \frac{1}{2}\,{u}_{t}{xt}+{{yu}}_{x}+{u}_{x}t-\displaystyle \frac{1}{2}\,{u}_{}x,\\ {T}^{y} & = & -[\displaystyle \frac{1}{2}\,{u}_{y}{y}^{2}-{u}_{y}{xy}-{u}_{y}{xt}-{{yu}}_{}+{u}_{}x],\end{array}\end{eqnarray}$

iii

(iii) $Q=1-\tfrac{1}{2}{y}^{2}+{xt}+y$:

$\begin{eqnarray}\begin{array}{rcl}{T}^{t} & = & \displaystyle \frac{1}{2}\,{u}_{x}-\displaystyle \frac{1}{4}\,{u}_{x}{y}^{2}+\displaystyle \frac{1}{2}\,{u}_{x}{xt}+\displaystyle \frac{1}{2}\,{{yu}}_{x}-\displaystyle \frac{1}{2}\,{u}_{}t,\\ {T}^{x} & = & -\displaystyle \frac{1}{2}\,{u}_{x}{u}_{}{y}^{2}+{u}_{x}{u}_{}{xt}+{u}_{x}{{yu}}_{}+{u}_{x}{u}_{}-\displaystyle \frac{1}{2}\,{u}_{}^{2}t-{u}_{{xx}}\\ & & +\displaystyle \frac{1}{2}\,{u}_{{xx}}{y}^{2}-{u}_{{\text{}}{xx}}{xt}-{u}_{{xx}}y+\displaystyle \frac{1}{2}\,{u}_{t}\\ & & -\displaystyle \frac{1}{4}\,{u}_{t}{y}^{2}+\displaystyle \frac{1}{2}\,{u}_{t}{xt}+\displaystyle \frac{1}{2}\,{u}_{t}y+{u}_{x}t-\displaystyle \frac{1}{2}\,{u}_{}x,\\ {T}^{y} & = & -[-{u}_{y}+\displaystyle \frac{1}{2}\,{u}_{y}{y}^{2}-{u}_{y}{xt}-{u}_{y}y-{{yu}}_{}+{u}_{}].\end{array}\end{eqnarray}$

On the basis of the group method [5, 7], for a one-parameter group of infinitesimal transformation
$\begin{eqnarray}\begin{array}{rcl}\overline{t} & = & t+\epsilon {\xi }_{t}(x,y,t,u)+O({\epsilon }^{2}),\\ \overline{x} & = & x+\epsilon {\xi }_{x}(x,y,t,u)+O({\epsilon }^{2}),\\ \overline{y} & = & y+\epsilon {\xi }_{y}(x,y,t,u)+O({\epsilon }^{2}),\\ \overline{u} & = & u+\epsilon \eta (x,y,t,u)+O({\epsilon }^{2}),\end{array}\end{eqnarray}$
where ε is a group parameter. The corresponding vector field is given by
$\begin{eqnarray}\begin{array}{rcl}V & = & {\xi }_{t}(x,y,t,u)\displaystyle \frac{\partial }{\partial t}+{\xi }_{x}(x,y,t,u)\displaystyle \frac{\partial }{\partial x}\\ & & +{\xi }_{y}(x,y,t,u)\displaystyle \frac{\partial }{\partial y}+{\eta }_{u}(x,y,t,u)\displaystyle \frac{\partial }{\partial u}.\end{array}\end{eqnarray}$
As this equation is second order, to solve this equation, the second prolongation Pr(2)V is required. Meanwhile, the invariant condition is given by:
$\begin{eqnarray}\begin{array}{l}{\eta }_{u}^{{tx}}+2{u}_{x}{\eta }_{u}^{x}+u{\eta }_{u}^{{xx}}+{u}_{{xx}}{\eta }_{u}\\ \quad -{\eta }_{u}^{{xxx}}+{\eta }_{u}^{{yy}}=0,\end{array}\end{eqnarray}$
where the total derivative operators is showed
$\begin{eqnarray}\begin{array}{rcl}{\eta }^{t} & = & {D}_{t}(\eta )-{u}_{t}{D}_{t}({\xi }_{t})-{u}_{x}{D}_{t}({\xi }_{x})-{u}_{y}{D}_{t}({\xi }_{y}),\\ {\eta }^{x} & = & {D}_{x}(\eta )-{u}_{t}{D}_{x}({\xi }_{t})-{u}_{x}{D}_{x}({\xi }_{x})-{u}_{y}{D}_{x}({\xi }_{y}),\\ {\eta }^{y} & = & {D}_{y}(\eta )-{u}_{t}{D}_{y}({\xi }_{t})-{u}_{y}{D}_{y}({\xi }_{y})-{u}_{y}{D}_{y}({\xi }_{y}),\\ & & \cdots \end{array}\end{eqnarray}$
Dx, Dt, Dy are functions of x, y and t, respectively.
The Lie point symmetry generators, as would be expected, also form an infinite algebra generated by the vector field
$\begin{eqnarray}\begin{array}{rcl}X & = & \left(-\displaystyle \frac{1}{2}{g}_{1}^{\prime} (t)y+{g}_{2}(t)\right){\partial }_{x}+{f}_{1}(t){\partial }_{y}+{\partial }_{t}\\ & & +\left(-\displaystyle \frac{1}{2}{f}_{1}^{\prime\prime} (t)y+{f}_{2}^{\prime} (t)\right){\partial }_{u},\\ Y & = & \left(\displaystyle \frac{2}{3}\,x-\displaystyle \frac{1}{2}{f}_{3}^{\prime} (t)y+{f}_{4}(t)\right){\partial }_{x}\\ & & +(y+{f}_{3}(t)){\partial }_{y}+\left(\displaystyle \frac{4}{3}t\right){\partial }_{t}\\ & & +\left(-\displaystyle \frac{1}{2}{f}_{3}^{\prime\prime} (t)y+{f}_{4}^{\prime} (t)-\displaystyle \frac{2}{3}u\right){\partial }_{u},\end{array}\end{eqnarray}$
from which we may conclude some ‘polynomial’ generators
$\begin{eqnarray}\begin{array}{l}{\partial }_{x},\qquad {\partial }_{y},\qquad {\partial }_{t},\qquad t{\partial }_{x}+{\partial }_{u},\\ y{\partial }_{x}-2t{\partial }_{y},\\ {t}^{2}{\partial }_{x}+2t{\partial }_{u},\qquad {ty}{\partial }_{x}-{t}^{2}{\partial }_{y}+y{\partial }_{u},\\ \quad \displaystyle \frac{1}{2}x{\partial }_{x}+\displaystyle \frac{3}{4}{\partial }_{y}+t{\partial }_{t}-\displaystyle \frac{1}{2}u{\partial }_{u}.\end{array}\end{eqnarray}$
The variety of solutions obtained above may be attributed to the richness of the symmetries and conservation laws obtained here.

5. Conclusions

In the present paper, on the basis of the (1+1)-dimensional Burgers equation, we have proposed a simple, direct, and efficient method, named as CBEE method, for solving NLEEs. The convenience, simplicity and effectiveness of the method is illustrated by solving the (2+1)-dimensional Burgers equation. The main results in the present paper map solutions of the (1+1)-dimensional Burgers equation onto the (2+1)-dimensional Burgers equation. By choosing appropriate parameters, the transformation yields new solutions for the (2+1)-dimensional Burgers equation. Although not all the NLEEs are CBEE solvable, it is may shed further light on the solutions of some high-dimensional NLEEs.
Also, we note that a variety of solutions were obtained above and this can be predicted and attributed to the large number of symmetries and conservation laws that the PDE under investigation generated.
It is should be noted that this method, of course, should also be easily extended to other more (2+1)-dimensional, (3+1)-dimensional, and even more high-dimensional NLEEs. The method, in this paper, might become a useful, promising and powerful technique for solving other NLEEs.

We are very grateful to the editors and reviewers for their comments, which greatly improved the quality of this article. This work is supported by Natural Science Foundation of Hebei Province, China (No. A2018207030), Youth Key Program of Hebei University of Economics and Business (2018QZ07), Key Program of Hebei University of Economics and Business (2020ZD11), Youth Team Support Program of Hebei University of Economics and Business, Study on system dynamics of scientific and technological innovation promoting the expansion and quality of residents’ consumption in Hebei Province(20556201D), Youth Top-notch Talent Support Program of Higher Education of Hebei Province of China (BJ2020011).

1
Ablowitz M J Segur H 1981 Solitons and the Inverse Scattering Transform Philadelphia SIAM

2
Yang Y Y Song J 2021 On the generalized eigenvalue problem of Rossby waves vertical velocity under the condition of zonal mean flow and topography Appl. Math. Lett. 121 107485

DOI

3
Yu D Zhang Z G Dong H H 2021 Bäcklund transformation, infinite number of conservation laws and fission properties of an integro-differential model for ocean internal solitary waves Commun. Theor. Phys. 73 035005

DOI

4
Bi Y Zhang Z Liu Q 2021 Research on nonlinear waves of blood flow in arterial vessels Commun. Nonlinear Sci. Numer. Simul. 102 105918

DOI

5
Bluman G W Kumei S 1989 Symmetries and Differential Equations New York Springer

6
Lou S Y 2018 Alice-Bob systems, $\hat{P}$-$\hat{T}$-$\hat{C}$ symmetry invariant and symmetry breaking soliton solutions J. Math. Phys. 59 083507

DOI

7
Olver P J 1986 Application of Lie Group to Differential Equation New York Springer

8
Bluman G W Cheviakov A Anco S 2010 Applications of Symmetry Methods to Partial Differential Equations New York Springer

9
Tian C 2001 Lie Groups and its Applications to Differential Equations Beijing Science Press (in Chinese)

10
Wang G W Liu X Q Ying Y Y 2013 Lie symmetry analysis to the time fractional generalized fifth-order KdV equation Commun. Nonlinear Sci. Numer. Simul. 18 2321 2326

DOI

11
Wang G W 2021 A novel (3+1)-dimensional sine-Gorden and sinh-Gorden equation: derivation, symmetries and conservation laws Appl. Math. Lett. 113 106768

DOI

12
Wang G W 2020 A (2+1)-dimensional sine-Gordon and sinh-Gordon equations with symmetries and kink wave solutions Nucl. Phys. B 953 114956

DOI

13
Wang G W Kara A H 2019 A (2+1)-dimensional KdV equation and mKdV equation: Symmetries, group invariant solutions and conservation laws Phys. Lett. A 383 728 731

DOI

14
Wang G W Wazwaz A M 2022 On the modified Gardner type equation and its time fractional form Chaos Solitons Fractals 155 111694

DOI

15
Wang G W 2021 A new (3+1)-dimensional Schröodinger equation: Derivation, soliton solutions and conser- vation laws Nonlinear Dyn. 104 1595 1602

DOI

16
Wang G W 2020 Symmetry analysis for a seventh-order generalized KdV equation and its fractional version in fluid mechanics Fractals 28 2050044

DOI

17
Rogers C Schief W K 2002 Bäcklund and Darboux Transformations, Geometry and Morden Applications in Soliton Theory Cambridge Cambridge University Press

18
Hu X B Zhu Z N Bäcklund A 1998 Transformation and nonlinear superposition formula for the Belov-Chaltikian lattice J. Phys. A: Math. Gen. 31 4755 4716

DOI

19
Gu C H 2005 Darboux Transformations in Integrable Systems Theory and their Applications to Geometry Berlin Springer

20
Hirota R 2004 The Direct Method in Soliton Theory Cambridge Cambridge University Press

21
Wang Q Chen Y Zhang H Q 2005 A new Jacobi elliptic function rational expansion method and its application to (1+1)-dimensional dispersive long wave equation Chaos Solitons Fractals 23 477 483

DOI

22
Hu W P 2020 Symmetry breaking of infinite-dimensional dynamic system Appl. Math. Lett. 103 106207

DOI

23
Hu W 2021 Coupling dynamic behaviors of flexible stretching hub-beam system Mech. Syst. Signal Process. 151 107389

DOI

24
Hu W Ye J Deng Z 2020 Internal resonance of a flexible beam in a spatial tethered system J. Sound Vib. 475 115286

DOI

25
Hu W 2022 Symplectic analysis on dynamic behaviors of tethered tug-debris system Acta Astronaut. 192 182 189

DOI

26
Wang M L Zhou Y B Li Z B 1996 Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics Phys. Lett. A 216 67 75

DOI

27
Wang M L Li X 2005 Extended F-expansion method and periodic wave solutions for the generalized Zakharov equations Phys. Lett. A 343 48 54

DOI

28
He J H 2005 Application of homotopy perturbation method to nonlinear wave equations Chaos Solitons Fractals 26 695 700

DOI

29
Liao S J 1998 Homotopy analysis method: A new analytic method for nonlinear problems Appl. Math. Mech. 19 957 962

DOI

30
Lou S Y Ma H C 2005 Non-Lie symmetry groups of (2+1)-dimensional nonlinear systems obtained from a simple direct method J. Phys. A: Math. Gen. 38 L129 L137

DOI

31
Lou S Y 2015 Consistent Riccati expansion for integrable systems Stud. Appl. Math. 134 372 402

DOI

32
Xian D Q Jiang Y Kang X R 2020 Consistent KdV expansion method and its applications to the KP equation J. Comput. Appl. Math. 366 112423

DOI

33
Burgers J M 1948 A mathematical model illustrating the theory of turbulence Adv. Appl. Mech. 1 171 199

34
Webb G M Zank G P 1990 Painlevé analysis of the two-dimensional Burgers equation J. Phys. A: Math. Gen. 23 5465 5477

DOI

35
Webb G M 1990 Painlevé analysis of the three-dimensional Burgers equation Phys. Lett. A 150 14 22

DOI

36
Webb G M Zank G P 1992 On the Painlevé analysis of the two-dimensional Burgers’ equation Nonlinear Anal.-Theor. 19 167 176

DOI

37
Taniuti T 1974 Reductive perturbation method and far fields of wave equations Prog. Theor. Phys. Suppl. 55 1 35

DOI

38
Bartuccelli M Pantano P Brugarino T 1983 Two-dimensional Burgers equation Lett. Al Nuovo Cimento Ser. 37 433 438

DOI

39
Bluman G W Anco S C 2002 Symmetry and Integration Methods for Differential Equations New York Springer

40
Poole D Hereman W 2010 The homotopy operator method for symbolic integration by parts and inversion of divergences with applications Appl. Anal. 89 433 455

DOI

Outlines

/