New Exact Solutions for Benjamin-Bona-Mahony-Burgers Equation

Chunhuan Xiang¹, Honglei Wang^2*

¹School of Public Health and Management, Chongqing Medical University, Chongqing, China.
²College of Medical Informatics, Chongqing Medical University, Chongqing, China.

DOI: 10.4236/ojapps.2020.108038                

Abstract

Obtaining the new solutions for the nonlinear evolution equation is a hot topic. Benjamin-Bona-Mahony-Burgers equation is this kind of equation, the solutions are very interest. Several new exact solutions for the nonlinear equation are obtained by using truncated expansion method in this paper. The numerical simulations with different parameters for the new exact solutions of Benjamin-Bona-Mahony-Burgers equation are given.

Keywords

Truncated Expansion Method, Benjamin-Bona-Mahony-Burgers Equation, Evolution, Nonlinear Equation

1. Introduction

Nonlinear partial differential equations are very important in the natural world; the solutions are the important way for us to know the nature. However, the different initial condition means different solutions exit. People had found out many powerful methods to obtain the solutions for nonlinear partial differential equations. Such as the useful methods-Soliton-Like Solutions, complex travelling wave, truncated expansion method, hyperbolic tangent method [1] [2] [3] [4] [5], high-order multi-symplectic schemes, simple fast method, Invariant-conserving finite difference algorithms, stable bound states [6] - [11], which help us deeply study the relation of the nature.

The Nonlinear Benjamin-Bona-Mahony-Burgers equation was proposed by Peregrine [12] and Benjamin [13], long waves on the surface of water in a channel with small-amplitude can be described by the nonlinear equation. The form of the Benjamin-Bona-Mahony-Burgers equation [14] is read as

$u_t+u_x+u{u}_x-m{u}_{xx}-u_{xxt}=0$ (1)

where the subscripts denote the partial derivatives of position x and time t. $u\left(x,t\right)$ is a real-valued function and $u_{xx}$ is considered as dissipative term. The following structure of this work is organized as follows: Section 2 is a brief introduction to the truncated expansion method and its Properties. In Section 3, applying the truncated expansion method, some new exact wave solutions for Nonlinear Benjamin-Bona-Mahony-Burgers equation are given The conclusion is summarized in the final.

2. The Truncated Expansion Method and Its Properties

The nonlinear partial differential equation with independent variables position x and time t is generally in the following form

$Q\left(u,u_x,u_t,u_{xt},u_{xx},u_{tt},\cdots \right)=0$ (2)

The above equation is a function about $u\left(x,t\right)$ variable with position x and time t, the subscripts denote the partial derivatives with x and t, respectively. The wave variable $\xi =ax+bt$ is applied to the Equation (2), which is changed into the following ordinary differential equation

$Q\left(u,u_{\xi },u_{\xi \xi },u_{\xi \xi \xi },\cdots \right)=0$ (3)

where $u_{\xi },u_{\xi \xi },\cdots$ denotes the derivative with respect to the same variable $\xi$.

Generally, the function $u\left(x,t\right)$ in terms of truncated expansion method [15] [16] can be expressed as follow

$u\left(ϵ\right)={\displaystyle {\sum }_{n=0}^N{A}_n{\phi }^n\left(ϵ,t\right)}$ (4)

where A_n is constant parameter; N is determined by balancing the linear term of the highest order derivative with nonlinear term of Equation (3).

The function $\phi \left(ϵ,t\right)$ has the nature

${\phi }_{ϵ}=\lambda +{\phi }^2$ (5)

where the $\lambda$ is a constant. The second and third derivative of $\phi$ are read as the follow

${\phi }_{ϵϵ}=2\lambda \phi +2{\phi }^3$ (6)

${\phi }_{ϵϵϵ}=2{\lambda }^2+8\lambda {\phi }^2+6{\phi }^4$ (7)

The function $\phi \left(ϵ,t\right)$ with the form of Equation (5) has the follow form with different $\lambda$

$\phi =-{\left(-\lambda \right)}^{\frac{1}{2}}\tanh \left[{\left(-\lambda \right)}^{\frac{1}{2}}ϵ\right],\lambda <0$ (8)

$\phi =-{\left(-\lambda \right)}^{\frac{1}{2}}\coth \left[{\left(-\lambda \right)}^{\frac{1}{2}}ϵ\right],\lambda <0$ (9)

$\phi =-\frac{1}{ϵ},\lambda =0$ (10)

$\phi =-{\left(\lambda \right)}^{\frac{1}{2}}\tan \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right],\lambda >0$ (11)

$\phi =-{\left(\lambda \right)}^{\frac{1}{2}}\cot \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right],\lambda >0$ (12)

3. The Truncated Expansion Method for Benjamin-Bona-Mahony-Burgers Equation

Consider the Nonlinear Benjamin-Bona-Mahony-Burgers equation in the form of Equation (1) with $ϵ=ax+bt$, we obtain

$b{u}_{ϵ}+a{u}_{ϵ}+au{u}_{ϵ}-m{a}^2{u}_{ϵϵ}-a^{2}b{u}_{ϵϵϵ}=0$ (13)

By using N to balance the highest order derivative term and the nonlinear term of Equation (13), that means $u{u}_{ϵ}$ and $u_{ϵϵϵ}$ are the nonlinear term and the highest order derivative, $\text{2}N+\text{1}=N+\text{3}$, we have N = 2. Then, Equation (4) reduced as

$u\left(ϵ\right)=A_0+A_1\phi +A_2{\phi }^2$ (14)

Substituting Equation (14) into the equation Equation (13) and collecting all terms with the same power term together and equating each coefficient of $\phi$ to zero, a set of simultaneous algebraic equations are yielded as follows:

$a\lambda A_1+a{A}_0\lambda A_1+\lambda A_{1}b-2a^2{\lambda }^2{A}_{1}b-2m{a}^2{A}_2{\lambda }^2=0$ (15)

$a\lambda A_1^2+2a\lambda A_2+2a\lambda A_0{A}_2+2\lambda A_{2}b-16b{a}^2{A}_2{\lambda }^2-2m{a}^2{A}_1\lambda =0$ (16)

$a{A}_1+a{A}_0{A}_1+3a\lambda A_1{A}_2+A_{1}b-8b{a}^2{A}_1\lambda -8m{a}^2{A}_2\lambda =0$ (17)

$a{A}_1^2+2a{A}_2+2a{A}_0{A}_2+2\lambda A_2^{2}a+2b{A}_2-2m{a}^2{A}_1-40b{a}^2{A}_2\lambda =0$ (18)

$3a{A}_1{A}_2-6a^2{A}_{1}b-6m{a}^2{A}_2=0$ (19)

$2A_2^2-24ab{A}_2=0$ (20)

Solving the above algebraic equations, we get the results:

$\begin{array}{l}m=\frac{\sqrt{24b^2+24ab-100\lambda a^2{b}^2}}{a};A_0=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab};\\ A_1=\frac{12am}{5};A_2=12ab\end{array}$ (21)

By using Equations (8)-(10) and (14), we obtained the solution of Equation (1) as the follow:

$\begin{array}{l}u{\left(x,t\right)}_1=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab}-\frac{12am}{5}{\left(-\lambda \right)}^{\frac{1}{2}}\tanh \left[{\left(-\lambda \right)}^{\frac{1}{2}}\left(ax+bt\right)\right]\\ +12ab{\left({\left(-\lambda \right)}^{\frac{1}{2}}\tanh \left[{\left(-\lambda \right)}^{\frac{1}{2}}\left(ax+bt\right)\right]\right)}^2,\lambda <0\end{array}$ (22)

$\begin{array}{l}u{\left(x,t\right)}_2=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab}-\frac{12am}{5}{\left(-\lambda \right)}^{\frac{1}{2}}\coth \left[{\left(-\lambda \right)}^{\frac{1}{2}}\left(ax+bt\right)\right]\\ +12ab{\left({\left(-\lambda \right)}^{\frac{1}{2}}\coth \left[{\left(-\lambda \right)}^{\frac{1}{2}}\left(ax+bt\right)\right]\right)}^2,\lambda <0\end{array}$ (23)

$u{\left(x,t\right)}_3=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab}-\frac{12am}{5\left(ax+bt\right)}+12ab{\left(\frac{-1}{ax+bt}\right)}^2,\lambda =0$ (24)

$\begin{array}{l}u{\left(x,t\right)}_4=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab}+\frac{12am}{5}{\left(\lambda \right)}^{\frac{1}{2}}\tan \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right]\\ +12ab{\left({\left(\lambda \right)}^{\frac{1}{2}}\tan \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right]\right)}^2,\lambda >0\end{array}$ (25)

$\begin{array}{l}u{\left(x,t\right)}_5=\frac{200\lambda a^2{b}^2-a^2{m}^2-ab-b^2}{25ab}-\frac{12am}{5}{\left(\lambda \right)}^{\frac{1}{2}}\cot \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right]\\ +12ab{\left(-{\left(\lambda \right)}^{\frac{1}{2}}\cot \left[{\left(\lambda \right)}^{\frac{1}{2}}ϵ\right]\right)}^2,\lambda >0\end{array}$ (26)

The rest section is the numerical simulation for the above sol are determined randomly in the corresponding interval. The contour plot of numerical simulation for Equation (22) with $a=\text{1}$ ; $b=0.5$ ; $\lambda =-2$ ; $x\in \left[-5,5\right]$ and $t\in \left[-5,5\right]$ is shown in Figure 1.

The contour plot of numerical simulation for Equation (23) with $a=\text{1}$ ; $b=0.5$ ; $\lambda =-1$ ; $x\in \left[3,4\right]$ and $t\in \left[3,4\right]$ is shown in Figure 2.

The contour plot of numerical simulation for Equation (24) with $a=\text{1}$ ; $b=0.5$ ; $x\in \left[1,2\right]$ and $t\in \left[1,2\right]$ is shown in Figure 3.

The contour plot of numerical simulation for Equation (25) with $a=\text{2}$ ; $b=6$ ; $\lambda =0.05$ ; $x\in \left[-1,1\right]$ and $t\in \left[-1,1\right]$ is shown in Figure 4.

Figure 1. The numerical simulation for Equation (22) with $x\in \left[-5,5\right]$ and $t\in \left[-5,5\right]$.

Figure 2. The numerical simulation for Equation (23) with $x\in \left[3,4\right]$ and $t\in \left[3,4\right]$.

Figure 3. The numerical simulation for Equation (24) with $x\in \left[1,2\right]$ and $t\in \left[1,2\right]$.

Figure 4. The numerical simulation for Equation (25) with $x\in \left[-1,1\right]$ and $t\in \left[-1,1\right]$.

Figure 5. The numerical simulation for Equation (26) $x\in \left[2,3\right]$ and $t\in \left[2,3\right]$.

The contour plot of numerical simulation for Equation (26) with $a=10$ ; $b=6$ ; $\lambda =0.01$ ; $x\in \left[2,3\right]$ and $t\in \left[2,3\right]$ is shown in Figure 5.

4. Conclusion

The Nonlinear Benjamin-Bona-Mahony-Burgers equation in the given form Equation (1) have been further studied, some new exact solutions for the equation are obtained by means of the truncated expansion method, which are shown in Equations (22)-(26). The numerical simulation results with contour plot are appended with Figures 1-5. The evolutions of the travelling wave are clear shown.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Zhu, Z.N. (1992) Soliton-Like Solutions of Generalized KdV Equation with External Force Term. Acta Physica Sinca, 41, 1561-1566.
[2]	Demiray, H. (2005) A Complex Travelling Wave Solution to the KdV-Burgers Equation. Physics Letters A, 344, 418-422. https://doi.org/10.1016/j.physleta.2004.09.087
[3]	Zhang, J.F. and Chen, F.Y. (2001) Truncated Expansion Method and New Exact Soliton-Like Solution of the General Variable Coefficient KdV Equation. Acta Physica Sinca, 50, 1648-1650.
[4]	Halim, A.A. and Leble, S.B. (2004) Analytical and Numerical Solution of a Coupled KdV-MKdV System. Chaos, Solitons & Fractals, 19, 99-108. https://doi.org/10.1016/S0960-0779(03)00085-7
[5]	Malfliet, M. and Wieers, E. (1996) A Nonlinear Theory of Charged-Particle Stopping in Non-Ideal Plasmas. Journal of Plasma Physics, 56, 441-443. https://doi.org/10.1017/S0022377800019401
[6]	Wang, Y.S. and Wang, B. (2005) High-Order Multi-Symplectic Schemes for the Nonlinear Klein-Gordon Equation. Applied Mathematics and Computation, 166, 608-632. https://doi.org/10.1016/j.amc.2004.07.007
[7]	Gegenhasi, Hu, X.B. and Wang, H.Y. (2007) A (2+1)-Dimensional Sinh-Gordon Equation and Its Pfaffian Generalization. Physics Letters A, 360, 439-447. https://doi.org/10.1016/j.physleta.2006.07.031
[8]	Liu, S.K., Fu, Z.T. and Liu, S.D. (2001) A Simple Fast Method in Finding Particular Solutions of Some Nonlinear PDE. Applied Mathematics and Mechanics, 22, 326-331. https://doi.org/10.1007/BF02437971
[9]	Zhang, J.L. and Wang, Y.M. (2003) Exact Solutions to Two Nonlinear Equations. Acta Physica Sinca, 52, 1574-1578.
[10]	Vu-Quoc, L. and Li, S. (1993) Invariant-Conserving Finite Difference Algorithms for the Nonlinear Klein-Gordon Equation. Computer Methods in Applied Mechanics and Engineering, 107, 341-391. https://doi.org/10.1016/0045-7825(93)90073-7
[11]	Bonanno, C. (2010) Existence and Multiplicity of Stable Bound States for the Nonlinear Klein-Gordon Equation. Nonlinear Analysis: Theory, Methods & Applications, 72, 2031-2046. https://doi.org/10.1016/j.na.2009.10.004
[12]	Peregrine, D.H. (1966) Calculations of the Development of an Undular Bore. Journal of Fluid Mechanics, 25, 321-330. https://doi.org/10.1017/S0022112066001678
[13]	Benjamin, T.B., Bona, J.L. and Mahony, J.J. (1972) Model Equations for Long Waves in Nonlinear Dispersive Systems. Philosophical Transactions of the Royal Society A, 272, 47-78. https://doi.org/10.1098/rsta.1972.0032
[14]	Vaneeva, O., Posta, S. and Sophocleous, C. (2017) Enhanced Group Classification of Benjamin-Bona-Mahony-Burgers Equations. Applied Mathematics Letters, 65, 19-25.
[15]	Yan, Z.Y. and Zhang, H.Q. (1999) Exact Soliton Solutions of the Variable Coefficient KdV-MKdV Equation with Three Arbitrary Functions. Acta Physica Sinica, 48, 1957-1961.
[16]	Chen, F.J. and Zhang, J.F. (2003) Soliton-Like Solution for the (2+1)-Dimensional Variable Coefficient Kadomtsev-Petviashvili Equation. Acta Armamentarii, 24, 389- 391.