Lie Symmetry Reductions and Exact Solutions of a Multidimensional Double Dispersion Equation

Jiali Yu¹, Fuzhi Li², Lianbing She³

¹School of Mathematics and Information Science, Guangzhou University, Guangzhou, China.
²School of Mathematics and Statistics, Southwest University, Chongqing, China.
³Department of Mathematics, Southwest University, Chongqing, China.

DOI: 10.4236/am.2017.85056                

Abstract

In this paper, based on classical Lie group method, we study a multidimensional double dispersion equation, and get its infinitesimal generator, symmetry group and similarity reductions. In particular, similarity solutions and travelling wave solutions of the multidimensional double dispersion equation are derived from the reduction equations.

Keywords

Lie Group, Multidimensional Double Dispersion Equation, Similarity Solutions, Travelling Wave Solutions

1. Introduction

The double dispersion Equation (1) was introduced as a mathematical model of nonlinear dispersive waves in various contexts (see [1] [2] [3] [4] )

$u_{tt}-u_{xx}+a{u}_{xxxx}-b{u}_{xxtt}=f{\left(u\right)}_{xx},$ (1)

where $u=u\left(x,t\right)$ is a real-valued function, $a\mathrm{,}b$ are positive real constants with $a\ge b$ . It also presents the plots of the instability/stability regions of travelling waves for various values of p, where $f\left(u\right)={\left|u\right|}^{p-1}u,p>1$ . The present paper provides an overview of results obtained in [5] and [6] for travelling wave solutions of the Equation (1).

Considering the possibility of energy exchange through lateral surfaces of the waveguide in the physical study of nonlinear wave propagation in waveguide, the longitudinal displacement $u\left(x\mathrm{,}t\right)$ of the rod satisfies the following double dispersion equation (DDE) (see [7] [8] [9] )

$u_{tt}-u_{xx}=\frac{1}{4}{\left(6u^2+a{u}_{tt}-b{u}_{xx}\right)}_{xx},$ (2)

and the general cubic DDE (CDDE)

$u_{tt}-u_{xx}=\frac{1}{4}{\left(c{u}^3+6u^2+a{u}_{tt}-b{u}_{xx}+d{u}_t\right)}_{xx},$ (3)

where $a\mathrm{,}b$ and $c$ are positive constants. The Equations (2) and (3) were studied in some literatures, the travelling wave solutions, depending upon the phase variable $z=x\pm ct$ were studied by Samsonov in [10] [11] .

In [12] [13] , Chen and Wang studied the initial-boundary value problem and the Cauchy problem of the following generalized double dispersion equation which includes above Equation (3) as special cases

$u_{tt}-u_{xx}-a{u}_{xxtt}+b{u}_{xxxx}-d{u}_{xxt}=f{\left(u\right)}_{xx},$ (4)

where $a>0,b>0$ and $d$ are constants.

Recently in [14] , the authors considered the Cauchy problem of the multidimensional nonlinear evolution equation

$u_{tt}-\Delta u-\Delta u_{tt}+{\Delta }^{2}u-k\Delta u_t=\Delta f\left(u\right),x\in R^n,t>0,$ (5)

$u\left(x,0\right)=u_0\left(x\right),u_t\left(x,0\right)=u_1\left(x\right),x\in R^n,$ (6)

where k is an arbitrary real constant. The authors gave the existence of local solution and the existence of global solution. And in [15] , Zhijian Yang et al. studied the existence of global attractor for the generalized double dispersion equation arising in elastic waveguide model. And there have been lots of research studies on the well-posedness, blowup, asymptotic behavior and other properties of solutions for both the IBVP and the IVP of the equation of type (6) (see [16] - [24] ) and references therein). While for the investigation on the global attractor to Equation (6), one can see [25] [26] [27] and references therein.

Symmetry reductions have several important applications in the context of differential equations. Since solutions of partial differential equations asymptotically tend to solutions of equations obtained by symmetry reduction, some of these special solutions will illustrate important physical phenomena (see [28] [29] [30] [31] [32] ). Solitary wave solutions and similarity solutions are usually applied to describe physical phenomena and to check on the accuracy and reliability of numerical algorithm, so deriving travelling wave solutions and similarity solutions has a great significance (see [33] [34] [35] ).

In [36] , the authors applied the method of Lie and the nonclassical method of Bluman and Cole to undertake the following equation

$u_{xxxx}+f{\left(u\right)}_{xx}=u_{tt},$ (7)

In [37] , the authors applied the Lie-group formalism to deduce symmetries of a generalized double dispersion equation

$u_{tt}-u_{xx}+u_{xxxx}-u_{xxtt}-{\left(u^n\right)}_{xx}=0,$ (8)

and they obtained exact solutions which can be expressed by various single and combined nondegenerative Jacobi elliptic function solutions.

In [38] , the authors consider the generalized double dispersion equation

$u_{tt}-u_{xx}-a{u}_{xxtt}+b{u}_{xxxx}-d{u}_{xxt}-f{\left(u\right)}_{xx}=0,$ (9)

where $a>0,b>0$ and d are constants. They study the functional forms $f\left(u\right)$ for which Equation (9) with $a\mathrm{,}b\ne 0$ admits classical symmetries.

In this paper, we consider the following multidimensional double dispersion equation

$u_{tt}-\Delta u-\Delta u_{tt}+{\Delta }^{2}u+k\Delta u_t=\Delta f\left(u\right),x\in R^n,t>0,$ (10)

where $n=3,f\left(u\right)={\left|u\right|}^p,p>1$ or $f\left(u\right)=u^{2k},k=1,2,\cdot \cdot \cdot ,$

In this paper, the symmetry group of the n-dimensional double dispersion Equation (10) is obtained by using the classical method in Section 2. In Section 3, we discuss the Lie symmetry group of Equation (10). Finally, we obtain similarity solutions or travelling wave solutions of Equation (10) by using similarity variables to obtain reduction equations, and solving the reduction equations in Section 4.

2. Lie Symmetry Analysis of the Double Dispersion Equation

In this section, we perform Lie symmetry analysis for Equation (10), and obtain its infinitesimal generator.

Theorem 1. [28] Assume

$\underset{\_}{V}={\displaystyle {\sum }_{i=1}^p}{\xi }^i\left(x,u\right)\frac{\partial }{\partial x^i}+{\displaystyle {\sum }_{\alpha =1}^p}{\varphi }_{\alpha }\left(x,u\right)\frac{\partial }{\partial x^{\alpha }},$

be a vector field on $X\times U$ , where $X=\left(x^1,x^2,\cdots ,x^p\right)$ , $U=\left(u^1,u^2,\cdots ,u^q\right)$ . Then its n-st prolongation is defined a vector on $X\times U$

$p{r}^{\left(n\right)}\underset{\_}{V}=\underset{\_}{V}+{\displaystyle {\sum }_{\alpha =1}^p}{\displaystyle {\sum }_{J=1}}{\varphi }_{\alpha }^J\left(x,u^{\left(n\right)}\right)\frac{\partial }{\partial u_J^{\alpha }},$

where, by definition

${\varphi }_{\alpha }^J\left(x,u^{\left(n\right)}\right)=D_J{Q}_{\alpha }+{\displaystyle {\sum }_{i=1}^p}{\xi }^i{u}_{Ji}^{\alpha },J=\left(j_1,\cdots ,j_l\right),0\le k\le n,1\le j_k\le p,$

where, $Q_{\alpha }={\varphi }_{\alpha }\left(x,u\right)-{\displaystyle {\sum }_{i=1}^p}{\xi }^i\left(x,u\right)\frac{\partial }{\partial x^i{u}_i^{\alpha }}$ , $\alpha =1,\cdots ,q$ , $Q\left(x,u\right)=\left(Q_1,\cdots ,Q_q\right)$

is referred to as the characteristic of the vector field $\underset{\_}{V}$ .

Here are four independent variables $x\mathrm{,}y\mathrm{,}z$ being spatial coordinates and t the time. According to the method of determining the infinitesimal generator of nonlinear partial differential equation [28] , we take the infinitesimal generator of equation as follows:

$\begin{array}{c}\underset{\_}{V}=\xi \left(x,y,z,t,u\right)\frac{\partial }{\partial x}+\eta \left(x,y,z,t,u\right)\frac{\partial }{\partial y}+\zeta \left(x,y,z,t,u\right)\frac{\partial }{\partial z}\\ +\tau \left(x,y,z,t,u\right)\frac{\partial }{\partial t}+\varphi \left(x,y,z,t,u\right)\frac{\partial }{\partial u},\end{array}$ (11)

be a vector field on $X\times U$ . Where $\xi \mathrm{,}\eta \mathrm{,}\zeta \mathrm{,}\tau \mathrm{,}\varphi$ are coefficient functions of the infinitesimal generator to be determined. We wish to determine all possible coe- fficient functions $\xi \mathrm{,}\eta \mathrm{,}\zeta \mathrm{,}\tau$ and $\varphi$ so that the corresponding one-parameter group $\exp \left(\epsilon \underset{\_}{V}\right)$ is a symmetry group of the double dispersion equation. Applying the forth prolongation of $\underset{\_}{V}$ to Equation (10), we find the invariance condition ${p{r}^{\left(4\right)}\underset{\_}{V}\left(\Delta \right)|}_{\Delta }=0$ , where $\Delta$ is $u_{tt}-\Delta u-\Delta u_{tt}+{\Delta }^{2}u+k\Delta u_t-\Delta \left(u^p\right)$ and with help of Maple software, we find the following system of symmetry equations

$\begin{array}{l}{\varphi }^{tt}-\left({\varphi }^{xx}+{\varphi }^{yy}+{\varphi }^{zz}\right)-\left({\varphi }^{xxtt}+{\varphi }^{yytt}+{\varphi }^{zztt}\right)\\ +\left({\varphi }^{xxxx}+{\varphi }^{yyyy}+{\varphi }^{zzzz}\right)+k\left({\varphi }^{xxt}+{\varphi }^{yyt}+{\varphi }^{zzt}\right)\\ =p\left(p-1\right)\left(p-2\right)u^{p-3}\left(u_x^2+u_y^2+u_z^2\right)\varphi \\ +p\left(p-1\right)u^{p-2}\left(2\nabla u{\varphi }^x+\Delta u\varphi \right)+p{u}^{p-1}\left({\varphi }^{xx}+{\varphi }^{yy}+{\varphi }^{zz}\right),\end{array}$ (12)

which must be satisfied whenever u satisfy Equation (10). Here ${\varphi }^{tt}\mathrm{,}{\varphi }^{xx}$ , etc. are

the coefficients of the second order derivatives $\frac{\partial }{\partial u_{tt}}\mathrm{,}\frac{\partial }{\partial u_{xx}}$ , etc. appearing in $p{r}^{\left(4\right)}\underset{\_}{V}$ .

According to Th. 1, ${\varphi }_{\alpha }^J\left(x\mathrm{,}{u}^{\left(n\right)}\right)=D_J\left({\varphi }_{\alpha }-\Sigma {\xi }^i{u}_i^{\alpha }\right)+\Sigma {\xi }^i{u}_{J\mathrm{,}i}^{\alpha }$ , we have

$\begin{array}{c}{\varphi }^{xx}=D_x^2\varphi -u_x{D}_x^2\xi -u_y{D}_x^2\eta -u_z{D}_x^2\varsigma -u_t{D}_x^2\tau \\ -2u_{xx}{D}_x\xi -2u_{xy}{D}_x\eta -2u_{xz}{D}_x\varsigma -2u_{xt}{D}_x\tau .\end{array}$

Similarly, we can get ${\varphi }^{tt}\mathrm{,}{\varphi }^{xxt}\mathrm{,}{\varphi }^{xxxx}$ , etc. we find the determining equations for the symmetry group of Equation (10) to be the following:

$\{\begin{array}{l}{\xi }_u={\eta }_u={\zeta }_u=0,{\xi }_t={\eta }_t={\zeta }_t=0,\hfill \\ {\varphi }_{uu}=0,\hfill \\ {\tau }_x={\tau }_y={\tau }_z={\tau }_u=0,\hfill \\ \left({\xi }_x+{\eta }_y+{\zeta }_z\right)-{\tau }_t=0,\hfill \\ 2\nabla {\varphi }_u-3\left({\xi }_{xx}+{\eta }_{yy}+{\zeta }_{zz}\right)=0,\hfill \\ 2\nabla {\varphi }_u-\left({\xi }_{xx}+{\eta }_{yy}+{\zeta }_{zz}\right)=0,\hfill \\ \Delta {\varphi }_u-2\left({\xi }_x+{\eta }_y+{\zeta }_z\right)=0,\hfill \\ 2{\varphi }_{tu}-k{\tau }_t-{\tau }_{tt}=0,\hfill \\ \nabla {\varphi }_{tu}+k\left({\xi }_{xx}+{\eta }_{yy}+{\zeta }_{zz}\right)-2k\nabla {\varphi }_u=0,\hfill \\ k\Delta {\varphi }_u-2\Delta {\varphi }_t+2{\varphi }_{tu}-{\tau }_{tt}=0,\hfill \\ u^{p-3}\left[u{\varphi }_u+\left(p-2\right)\varphi +2{\tau }_{t}u\right]=0,\hfill \\ -6\Delta {\varphi }_u+{\varphi }_{ttu}-k{\varphi }_{tu}+p\left(p-1\right)u^{p-2}\eta +2p{u}^{p-1}{\tau }_t=0,\hfill \\ {\Delta }^2\varphi -p{u}^{p-1}\Delta \varphi -\Delta \varphi -\Delta {\varphi }_{tt}+{\varphi }_{tt}+k\Delta {\varphi }_t=0,\hfill \\ 2p\left(p-1\right)u^{p-2}\nabla \varphi -4\nabla \Delta {\varphi }_u+2p{u}^{p-1}\nabla {\varphi }_u+2\nabla {\varphi }_u+2\nabla {\varphi }_{ttu}-2k\nabla {\varphi }_{tu}\hfill \\ -p{u}^{p-1}\left({\xi }_{xx}+{\eta }_{yy}+{\zeta }_{zz}\right)+\left({\xi }_{xxxx}+{\eta }_{yyyy}+{\zeta }_{zzzz}\right)-\left({\xi }_{xx}+{\eta }_{yy}+{\zeta }_{zz}\right)=0,\hfill \end{array}$

Since we have now satisfied all the determining equations, we conclude that general infinitesimal symmetry of Equation (10) has coefficient functions of the following form:

$\xi =c_1,$

$\eta =c_2,$

$\zeta =c_3,$

$\tau =c_4,$

where $c_1\mathrm{,}\cdots \mathrm{,}{c}_4$ are arbitrary constants. Thus the Lie-algebra of infinitesimal of the double dispersion equation is spanned by four vector fields:

$V_1=\frac{\partial }{\partial x},$

$V_2=\frac{\partial }{\partial y},$

$V_3=\frac{\partial }{\partial z},$

$V_4=\frac{\partial }{\partial t},$

so we have

$\underset{\_}{V}=c_1{V}_1+c_2{V}_2+c_3{V}_3+c_4{V}_4.$

3. Symmetry Groups of the Double Dispersion Equation

In this section, in order to get some exact solutions from a known solution of Equation (10), we should find the one-parameter symmetry groups $g_i\mathrm{<mo>\; :\; </mo>}\left(x\mathrm{,}y\mathrm{,}z\mathrm{,}y\mathrm{,}u\right)\to \left(\bar{x}\mathrm{,}\bar{y}\mathrm{,}\bar{z}\mathrm{,}\bar{t}\mathrm{,}\bar{u}\right)$ of corresponding infinitesimal generators. To get the Lie symmetry groups, we should solve the following initial problems of ordinary differential equations:

$\{\begin{array}{l}\frac{\text{d}\left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)}{\text{d}\epsilon }=\left(\xi ,\eta ,\zeta ,\tau ,\varphi \right)\\ {\left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)|}_{\epsilon =0}=\left(x,y,z,t,u\right),\end{array}$ (13)

where

$\xi =\xi \left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$ ,

$\eta =\eta \left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$ ,

$\zeta =\zeta \left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$ ,

$\tau =\tau \left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$ ,

$\varphi =\varphi \left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$

and $\epsilon$ is a group parameter.

For the infinitesimal generator $\underset{\_}{V}=c_1{V}_1+c_2{V}_2+c_3{V}_3+c_4{V}_4,$ we will take the following different values to obtain the corresponding infinitesimal generators:

Case 1. $C_1=1,C_2=C_3=C_4=0$ , the infinitesimal generator is $V_1=\frac{\partial }{\partial x},$

Case 2. $C_2=1,C_1=C_3=C_4=0$ , the infinitesimal generator is $V_2=\frac{\partial }{\partial y},$

Case 3. $C_3=1,C_1=C_2=C_4=0$ , the infinitesimal generator is $V_3=\frac{\partial }{\partial z},$

Case 4. $C_4=1,C_1=C_2=C_3=0$ , the infinitesimal generator is $V_4=\frac{\partial }{\partial t},$

Case 5. $C_1=C_2=C_3=1,C_4=0$ , the infinitesimal generator is $V_5=V_1+V_2+V_3=\frac{\partial }{\partial x}+\frac{\partial }{\partial y}+\frac{\partial }{\partial z},$

Case 6. $C_1=C_2=C_3=C_4=1$ , the infinitesimal generator is $V_6=V_1+V_2+V_3+V_4=\frac{\partial }{\partial x}+\frac{\partial }{\partial y}+\frac{\partial }{\partial z}+\frac{\partial }{\partial t},$

Case 7. $C_1=C_2=C_3=\lambda ,C_4=\beta$ , the infinitesimal generator is $V_7=V_1+V_2+V_3+V_4=\lambda \left(\frac{\partial }{\partial x}+\frac{\partial }{\partial y}+\frac{\partial }{\partial z}\right)+\beta \frac{\partial }{\partial t},$

The one-parameter groups $G_i$ generated by the $V_i$ . The entries give the transformed point $\exp \left(\epsilon V_i\right)\left(x,y,z,t,u\right)=\left(\bar{x},\bar{y},\bar{z},\bar{t},\bar{u}\right)$ :

$G_1\mathrm{<mo>\; :\; </mo>}\left(x+\epsilon \mathrm{,}y\mathrm{,}z\mathrm{,}t\mathrm{,}u\right)\mathrm{,}$

$G_2\mathrm{<mo>\; :\; </mo>}\left(x\mathrm{,}y+\epsilon \mathrm{,}z\mathrm{,}t\mathrm{,}u\right)\mathrm{,}$

$G_3\mathrm{<mo>\; :\; </mo>}\left(x\mathrm{,}y\mathrm{,}z+\epsilon \mathrm{,}t\mathrm{,}u\right)\mathrm{,}$

$G_4\mathrm{<mo>\; :\; </mo>}\left(x\mathrm{,}y\mathrm{,}z\mathrm{,}t+\epsilon \mathrm{,}u\right)\mathrm{,}$

$G_5\mathrm{<mo>\; :\; </mo>}\left(x+\epsilon \mathrm{,}y+\epsilon \mathrm{,}z+\epsilon \mathrm{,}t\mathrm{,}u\right)\mathrm{,}$

$G_6\mathrm{<mo>\; :\; </mo>}\left(x+\epsilon \mathrm{,}y+\epsilon \mathrm{,}z+\epsilon \mathrm{,}t+\epsilon \mathrm{,}u\right)\mathrm{,}$

$G_7\mathrm{<mo>\; :\; </mo>}\left(x+\lambda \epsilon \mathrm{,}y+\lambda \epsilon \mathrm{,}z+\lambda \epsilon \mathrm{,}t+\beta \epsilon \mathrm{,}u\right)\mathrm{,}$

where $G_1\mathrm{,}{G}_2\mathrm{,}{G}_3$ are space translations, $G_4$ is a time translation. $\epsilon$ is an arbitrary constant.

Theorem 2. If $u=f\left(x,y,z,t\right)$ is a known solution of Equation (10), then by using the symmetry groups $G_i\left(i=1,2,3,4\right)$ , so are the functions

$u_1=f\left(x-\epsilon ,y,z,t\right),$

$u_2=f\left(x,y-\epsilon ,z,t\right),$

$u_3=f\left(x,y,z-\epsilon ,t\right),$

$u_4=f\left(x,y,z,t-\epsilon \right),$

$u_5=f\left(x-\epsilon ,y-\epsilon ,z-\epsilon ,t\right),$

$u_6=f\left(x-\epsilon ,y-\epsilon ,z-\epsilon ,t-\epsilon \right),$

$u_7=f\left(x-\lambda \epsilon ,y-\lambda \epsilon ,z-\lambda \epsilon ,t-\lambda \epsilon \right),$

where $\epsilon$ is any real constant.

4. Symmetry Reductions and Exact Solutions of the Double Dispersion Equation

In the previous sections, we obtained the infinitesimal generators $V_i\left(i=1,2,\cdots ,7\right)$ . In this section, we will get similarity variables and its reduction equations, and we obtain similarity solutions and travelling wave solutions of Equation (10) by solving the reduction equations.

Case 1. For the infinitesimal generator $V_1=\frac{\partial }{\partial x}$ , the similarity variables are $r=t,F\left(r\right)=\frac{1}{3}x-tu$ , and the group-invariant solution is $u=\frac{\frac{1}{3}x-F\left(r\right)}{t}$ , subs-

tituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$F_r=0,$ (14)

Obviously, $F=c_1$ is a solution of Equation (14), where $c_1$ is an arbitrary constant. Therefore, Equation (10) has a similarity solution as follows:

$u=\frac{\frac{1}{3}x-c_1}{t}$ (15)

where $p=3$ .

Case 2. For the infinitesimal generator $V_2=\frac{\partial }{\partial y}$ , the similarity variables are $r=t,F\left(r\right)=\frac{1}{3}y-tu$ , and the group-invariant solution is $u=\frac{\frac{1}{3}y-F\left(r\right)}{t}$ , subs-

tituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$F_r=0,$ (16)

Obviously, $F=c_2$ is a solution of Equation (16), where $c_2$ is an arbitrary constant. Therefore, Equation (10) has a similarity solution as follows:

$u=\frac{\frac{1}{3}y-c_2}{t}$ (17)

where $p=3$ .

Case 3. For the infinitesimal generator $V_3=\frac{\partial }{\partial z}$ , the similarity variables are $r=t,F\left(r\right)=\frac{1}{3}z-tu$ , and the group-invariant solution is $u=\frac{\frac{1}{3}z-F\left(r\right)}{t}$ , subs-

tituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$F_r=0,$ (18)

Obviously, $F=c_3$ is a solution of Equation (18), where $c_3$ is an arbitrary constant. Therefore, Equation (10) has a similarity solution as follows:

$u=\frac{\frac{1}{3}z-c_1}{t}$ (19)

where $p=3$ .

Case 4. For the infinitesimal generator $V_4=\frac{\partial }{\partial t}$ , the similarity variables are

$r=x+y+z,F\left(r\right)=u$ , and the group-invariant solution is $u=F\left(r\right)$ , substituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$\Delta F\left(r\right)-{\Delta }^{2}F\left(r\right)+\Delta \left(F{\left(r\right)}^p\right)=0,$ (20)

Case 5. For the infinitesimal generator $V_5=\frac{\partial }{\partial x}+\frac{\partial }{\partial y}+\frac{\partial }{\partial z}$ , the similarity variables are $F\left(r\right)=\frac{1}{3}\left(x+y+z\right)-tu$ , and the group-invariant solution is $u=\frac{\frac{1}{3}\left(x+y+z\right)-F\left(r\right)}{t}$ , substituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$F_r=0,$ (21)

Obviously, $F=c_5$ is a solution of Equation (21), where $c_5$ is an arbitrary constant. Therefore, Equation (10) has a similarity solution as follows:

$u=\frac{\frac{1}{3}\left(x+y+z\right)-c_5}{t}$ (22)

where $p=3$ .

Case 6. For the infinitesimal generator $V_6=\frac{\partial }{\partial x}+\frac{\partial }{\partial y}+\frac{\partial }{\partial z}+\frac{\partial }{\partial t}$ , the similarity

variables are $r=x+y+z-\lambda t,F\left(r\right)=u$ , and the group-invariant solution is $u=F\left(r\right)$ , If we assume k = 0 of Equation (10), substituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$\left(1-{\lambda }^2\right)F_{rrrr}-p\left(p-1\right)F^{p-2}{F}_r^2-p{F}^{p-1}{F}_{rr}+\left({\lambda }^2-1\right)F_{rr}=0,$ (23)

So that Equation (23) is solvable in terms of Jacobi elliptic function, following the method described in [37] .

If ${\lambda }^2=\frac{4}{3}\sqrt{\frac{c_6}{5}}+1,F\left(r\right)=\sqrt{5c_6}d{n}^2\left(r,\sqrt{\frac{8}{5}}\right)$ , where $dn\left(r\mathrm{,}\sqrt{\frac{8}{5}}\right)$ is the Jacobi

elliptic of the third kind function. Therefore, Equation (10) has a travelling wave solution as follows:

$u\left(x,y,z,t\right)=\sqrt{5c_6}d{n}^2\left(x+y+z-\sqrt{\frac{4}{3}\sqrt{\frac{c_6}{5}}+1}t,\sqrt{\frac{8}{5}}\right)$ (24)

where $p=2$ , $c_6$ is an arbitrary constant.

If ${\lambda }^2=\frac{4}{3}\sqrt{\frac{c_6}{5}}+1,F\left(r\right)=\sqrt{5c_6}c{n}^2\left(r,\sqrt{\frac{8}{5}}\right)$ , where $cn\left(r\mathrm{,}\sqrt{\frac{8}{5}}\right)$ is the Jacobi

elliptic cosine function. Therefore, Equation (10) has a travelling wave solution as follows:

$u\left(x,y,z,t\right)=\sqrt{5c_6}c{n}^2\left(x+y+z-\sqrt{\frac{4}{3}\sqrt{\frac{c_6}{5}}+1}t,\sqrt{\frac{8}{5}}\right)$ (25)

where $p=2$ , $c_6$ is an arbitrary constant.

Case 7. For the infinitesimal generator $V_7=\left(\lambda \frac{\partial }{\partial x}+\frac{\partial }{\partial x}+\frac{\partial }{\partial x}\right)+\beta \frac{\partial }{\partial t}$ , following

the method described in [38] , the similarity variables are $r=\beta \left(x+y+z\right)-\lambda t,$ $F\left(r\right)=u$ , and the group-invariant solution is $u=F\left(r\right)$ , Substituting the group-invariant solution into Equation (10), we obtain the following reduction equation

$\begin{array}{l}\left({\lambda }^2-{\beta }^2\right){\beta }^2{F}_{rrrr}+k\lambda {\beta }^2{F}_{rrr}-\left({\lambda }^2-{\beta }^2\right){\beta }^2{F}_{rr}\\ -{\beta }^{2}p{F}^{p-1}{F}_{rr}-{\beta }^{2}p\left(p-1\right)F^{p-2}{F}_r^2=0,\end{array}$ (26)

Integrating twice with respect to r, we get

$\left({\lambda }^2-{\beta }^2\right){\beta }^2{F}_{rr}+k\lambda {\beta }^2{F}_r-\left({\lambda }^2-{\beta }^2\right)F+{\beta }^2{F}^p=0,$ (27)

Let us assume that Equation (27) has solution of the form

$F=c_7{H}^{c_8}\left(r\right),$ (28)

where $c_7\mathrm{,}{c}_8$ are arbitrary constants, and $H\left(r\right)$ is a solution of the Jacobi equation

$H_r^2=\mu +r{H}^2+H^4,$ (29)

If $H\left(r\right)=sn\left(r,m\right)$ , where $sn\left(r,m\right)$ is Jacobi elliptic sine function, we obtain $F\left(r\right)=rsn\left(r,m\right)$ is a solution of Equation (27). By substituting $F\left(r\right)$ $=rsn\left(r,m\right)$ into (27), we obtain the equation

$\begin{array}{l}-{\beta }^{2}rd{n}^2\left(r,m\right){\lambda }^2-m{\beta }^{2}rc{n}^2\left(r,m\right)s{n}^2\left(r,m\right){\lambda }^2\\ -rsn\left(r,m\right){\lambda }^2+{\beta }^{2}rc{n}^2\left(r,m\right)d{n}^2\left(r,m\right)s{n}^{-1}\left(r,m\right){\lambda }^2\\ -{\beta }^{2}rc{n}^2\left(r,m\right)d{n}^2\left(r,m\right)s{n}^{-1}\left(r,m\right){\lambda }^2+k{\beta }^{2}rcn\left(r,m\right)dn\left(r,m\right)\\ +{\beta }^{2}rsn\left(r,m\right)+m{\beta }^{4}rc{n}^2\left(r,m\right)n\left(r,m\right)\\ -{\beta }^{4}rc{n}^2\left(r,m\right)d{n}^2\left(r,m\right)s{n}^{-1}\left(r,m\right)\\ +{\beta }^{4}rc{n}^2\left(r,m\right)d{n}^2\left(r,m\right)s{n}^{-1}\left(r,m\right)\\ +{\beta }^{2}rs{n}^2\left(r,m\right)=0,\end{array}$ (30)

If ${\lambda }^2=2{\beta }^2$ , Equation (10) has a travelling wave solution as follows:

$u\left(x\mathrm{,}y\mathrm{,}z\mathrm{,}t\right)=\left[\beta \left(x+y+z\right)-\lambda t\right]sn\left(\beta \left(x+y+z\right)-\lambda t\mathrm{,}-1\right)\mathrm{,}$ (31)

where $m=-1,p=2$ . As $sn\left(\beta \left(x+y+z\right)-\lambda t\mathrm{,1}\right)=tanh\left(\beta \left(x+y+z\right)-\lambda t\right)$ , we obtain $F\left(r\right)=\frac{1}{4}tanh\left(r\right)$ is a solution of Equation (27), yb substituting $F\left(r\right)$ $=\frac{1}{4}tanh\left(r\right)$ into Equation (27), we obtain a travelling wave solution as follows:

$u\left(x\mathrm{,}y\mathrm{,}z\mathrm{,}t\right)=\frac{1}{4}tanh\left(x+y+z-\frac{t}{2}\right)\mathrm{,}$ (32)

where $k=8,p=3$

If $F\left(r\right)=rcn\left(r\mathrm{,}m\right)$ is a solution of Equation (27). By substituting $F\left(r\right)$ $=rcn\left(r\mathrm{,}m\right)$ into (27), when ${\lambda }^2=2{\beta }^2$ , Equation (10) has a travelling wave solution as follows:

$u\left(x\mathrm{,}y\mathrm{,}z\mathrm{,}t\right)=\left[\beta \left(x+y+z\right)-\lambda t\right]cn\left(\beta \left(x+y+z\right)-\lambda t\mathrm{,}-1\right)\mathrm{,}$ (33)

where $m=1,p=2$

5. Conclusion

In this paper, we study the symmetry reductions and explicit solutions of a multidimensional double dispersion equation by means of classical Lie group method. First, we get the symmetry groups and the infinitesimal generators of Equation (10). Then, we discuss the Lie symmetry groups of the multidimensional double dispersion equation and obtain the group-invariant solution. Finally, we obtain similarity solutions and travelling wave solutions of Equation (10) using similarity variables.

Acknowledgments

This work was supported by the Natural Science Foundation of Education of Guizhou Province: KY[2016]103.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Samsonov, A.M. (2001) Strain Solitons in Solids and How to Construct Them. Chapman and Hall, Boca Raton. https://doi.org/10.1201/9781420026139
[2]	Porubov, A.V. (2003) Amplification of Nonlinear Strain Waves in Solids. World Scientific, Singapore. https://doi.org/10.1142/5238
[3]	Christov, C.I., Maugin, G.A. and Porubov, A.V. (2007) On Boussinesq’s Paradigm in Nonlinear Wave Propagation. Comptes Rendus Mecanique, 335, 521-535. https://doi.org/10.1016/j.crme.2007.08.006
[4]	Engelbrecht, J., Salupere, A. and Tamm, K. (2011) Waves in Microstructured Solids and the Boussinesq Paradigm. Wave Motion, 48, 717-726.
[5]	Erbay, H.A., Erbay, S. and Erkip, A. (2015) Existence and Stability of Traveling Waves for a Class of Nonlocal Nonlinear Equations. Journal of Mathematical Analysis and Applications, 425, 307-336.
[6]	Erbay, H.A., Erbay, S. and Erkip, A. (2016) Instability and Stability Properties of Traveling Waves for the Double Dispersion Equation. Nonlinear Analysis, 133, 1-14.
[7]	Pasternak, N.L. (1954) New Method for Calculation of Foundation on the Elastic Basement. Gosstroiizdat, Moscow. (In Russian)
[8]	Samsonov, A.M. and Sokurinskaya, E.V. (1989) Energy Exchange between Nonlinear Waves in Elastic Waveguides and External Media. In: Engelbrecht, P.J., Ed., Nonlinear Waves in Active Media, Springer, Berlin, 99-104. https://doi.org/10.1007/978-3-642-74789-2_13
[9]	Samsonov, A.M. (1994) Nonlinear Strain Waves in Elastic Waveguide. In: Jeffrey, A. and Engelbrecht, J., Eds., Nonlinear Waves in Solids, CISM Courses and Lectures, Vol. 341, Springer, Wien, 349-382. https://doi.org/10.1007/978-3-7091-2444-4_6
[10]	Samsonov, A.M. (1988) On Existence of Longitudinal Strain Solitons in an Nonlinearly Elastic Rod. Soviet Physics-Doklady, 4, 298-300.
[11]	Samsonov, A.M. (1993) On Some Exact Travelling Wave Solutions for Nonlinear Hyperbolic Equation. In: Fusco, D. and Jeffrey, A., Eds., Nonlinear Waves and Dissipative Effects, Pitmann Research Notes in Mathematics Series, Vol. 227, Longman Scientific Technical, Longman, 123-132.
[12]	Chen, G., Wan, Y. and Wang, S. (2004) Initial Boundary Value Problem of the Generalized Cubic Double Dispersion Equation. Journal of Mathematical Analysis and Applications, 299, 563-577.
[13]	Wang, S. and Chen, G. (2006) Cauchy Problem of the Generalized Double Dispersion Equation. Nonlinear Analysis, 64, 159-173.
[14]	Polat, N. and Ertas, A. (2009) Existence and Blow-Up of Solution of Cauchy Problem for the Generalized Damped Multidimensional Boussinesq Equation. Journal of Mathematical Analysis and Applications, 349, 10-20.
[15]	Yang, Z., Feng, N. and Ma, T.T. (2015) Global Attractor for the Generalized Double Dispersion Equation. Nonlinear Analysis, 115, 103-116.
[16]	Babaoglu, C., Erbay, H.A. and Erkip, A. (2013) Global Existence and Blow-Up of Solutions for a General Class of Doubly Dispersive Nonlocal Nonlinear Wave Equations. Nonlinear Analysis, 77, 82-93.
[17]	Chen, G.W. and Xue, H.X. (2008) Periodic Boundary Value Problem and Cauchy Problem of the Generalized Cubic Double Dispersion Equation. Acta Mathematica Scientia, 28, 573-587.
[18]	Kato, M., Wang, Y.Z. and Kawashima, S. (2013) Asymptotic Behavior of Solutions to the Generalized Cubic Double Dispersion Equation in One Space Dimension. Kinetic and Related Models, 6, 969-987. https://doi.org/10.3934/krm.2013.6.969
[19]	Kawashima, S. and Wang, Y.Z. (2015) Global Existence and Asymptotic Behavior of Solutions to the Generalized Cubic Double Dispersion Equation. Analysis and Applications, 13, 233-254.
[20]	Liu, Y.C. and Xu, R.Z. (2008) Potential Well Method for Cauchy Problem of Generalized Double Dispersion Equations. Journal of Mathematical Analysis and Applications, 338, 1169-1187.
[21]	Liu, Y.C. and Xu, R.Z. (2008) Potential Well Method for Initial Boundary Value Problem of the Generalized Double Dispersion Equations. Communications on Pure and Applied Mathematics, 7, 63-81.
[22]	Wang, S.B. and Da, F. (2013) On the Asymptotic Behaviour of Solution for the Generalized Double Dispersion Equation. Applied Analysis, 92, 1179-1193. https://doi.org/10.1080/00036811.2012.661044
[23]	Xu, R.Z., Liu, R.Z. and Liu, Y.C. (2009) Global Existence and Nonexistence of Solution for Cauchy Problem of Multidimensional Double Dispersion Equations. Journal of Mathematical Analysis and Applications, 359, 739-751.
[24]	Xu, R.Z., Liu, Y.C. and Yu, T. (2009) Global Existence of Solution for Cauchy Problem of Multidimensional Generalized Double Dispersion Equations. Nonlinear Analysis: Theory, Methods & Applications, 71, 4977-4983.
[25]	Dai, Z.D. and Guo, B.L. (1999) Long Time Behavior of Nonlinear Strain Waves in Elastic Waveguides. Journal of Partial Differential Equations, 12, 301-312.
[26]	Dai, Z.D. and Guo, B.L. (2000) Global Attractor of Nonlinear Strain Waves in Elastic Waveguides. Acta Mathematica Scientia, 20, 322-334.
[27]	Yang, Z.J. and Li, K. (2013) Longtime Dynamics for a Elastic Waveguide Model. Discrete and Continuous Dynamical Systems Series S, 797-806.
[28]	Tian, C. (2001) Applications of Lie Groups to Differential Equations. Science Press, Beijing.
[29]	Wang, M.L., Zhou, Y. and Li, Z. (1996) Application of a Homogeneous Balance Method to Exact Solutions of Nonlinear Equations in Mathematical Physics. Physics Letters A, 216, 67-75.
[30]	Fan, E. (2000) Extended Tanh-Function Method and Its Applications to Nonlinear Equations. Physics Letters A, 277, 212-218.
[31]	Bluman, J.C.W. and Cole, J.D. (1974) Similarity Methods for Differential Equations. Springer, Berlin. https://doi.org/10.1007/978-1-4612-6394-4
[32]	Olver, P.J. (2000) Applications of Lie Groups to Differential Equations. Springer Verlag, Berlin.
[33]	Liu, H., Li, J. and Liu, L. (2009) Lie Group Classifications and Exact Solutions for Two Variable-Coefficient Equations. Applied Mathematics and Computation, 215, 2927-2935.
[34]	Stephani, H. and MacCallum, M. (1989) Differential Equations: Their Solution Using Symmetries. Cambridge University Press, Cambridge.
[35]	Cantwell, B. (2002) Introduction to Symmetry Analysis. Cambridge University Press, Cambridge.
[36]	Clarkson, P.A. and Priestley, T.J. (1996) Symmetries of a Generalised Boussinesq Equation. Electronic Edition.
[37]	Bruzón, M.S. and Gandarias, M.L. (2008) Similarity Reductions of a Generalized Double Dispersion Equation. 79th Annual Meeting of the International Association of Applied Mathematics and Mechanics, 8, 10587-10588,. https://doi.org/10.1002/pamm.200810587
[38]	Bruzón, M.S. and Gandarias, M.L. (2009) Travelling Wave Solutions for a Generalized Double Dispersion Equation. Nonlinear Analysis, 71, 2109-2117.