Solitary wave solutions of the fourth order Boussinesq equation through the exp(–Ф(η))-expansion method

Abstract The exp(–Ф(η))-expansion method is an ascending method for obtaining exact and solitary wave solutions for nonlinear evolution equations. In this article, we implement the exp(–Ф(η))-expansion method to build solitary wave solutions to the fourth order Boussinesq equation. The procedure is simple, direct and useful with the help of computer algebra. By using this method, we obtain solitary wave solutions in terms of the hyperbolic functions, the trigonometric functions and elementary functions. The results show that the exp(–Ф(η))-expansion method is straightforward and effective mathematical tool for the treatment of nonlinear evolution equations in mathematical physics and engineering. Mathematics subject classifications 35C07; 35C08; 35P99


Background
The world around us is inherently nonlinear (He 2009) and nonlinear evolution equations (NLEEs) are widely used as models to describe complex physical phenomena in various fields of science and engineering, especially in solid-state physics, plasma physics, fluid mechanics, biology etc. One of the fundamental problems for these models is to obtain their travelling wave solutions as well as solitary wave solutions. In particular, various methods have been utilized to explore different kinds of solutions of physical problems described by nonlinear evolution equations. In the numerical methods, stability and convergence should be considered, so as to avoid divergence or inappropriate results. However, in recent times, a variety of analytical and semi-analytical methods have been developed and use for solving NLEEs, for instance, the inverse scattering transform (Ablowitz and Clarkson 1991), the complex hyperbolic function method (Chow 1995;Zayed et al. 2006), the rank analysis method (Feng 2000), the ansatz method (Hu 2001a, b), the (G′/G)-expansion method (Wang et al. 2008;Bekir 2008;Neyrame et al. 2012;Akbar et al. 2012;Alam and Akbar 2013;Alam et al. 2014), the Exp-functions method (He and Wu 2006), the modified simple equation method (Jawad et al. 2010;, the Jacobi elliptic function method (Chen and Wang 2005;Liu 2005), the Adomian decomposition method (Adomian 1994;Wazwaz 2002), the homogeneous balance method (Wang 1995;Zayed et al. 2004), the F-expansion method (Wang and Zhou 2003;Wang and Li 2005), the Backlund transformation method (Miura 1978), the Darboux transformation method (Matveev and Salle 1991), the homotopy perturbation method (Mohyud-Din 2007; Mohyud-Din and Noor 2009), the generalized Riccati equation method (Yan and Zhang 2001), the tanh-function method (Wazwaz 2005), the Hirota's bilinear method (Hirota 2004), the auxiliary equation method (Sirendaoreji 2007), the exp(-Ф(η))-expansion method  etc.
The objective of this article is to implement the potential exp(-Ф(η))-expansion method to search solitary wave solutions for nonlinear evolution equations via the fourth order Boussinesq equation. In former literature, the solitary wave solutions to the Boussinesq equation have not been studied by this method.
The article is organized as follows: In Methodology Section, we give the description of the exp(-Ф(η))-expansion method. We apply this method to the fourth order Boussinesq equation in the Application Section. 2D and 3D graphs are given in the Graphical Representations of the Solutions Section. Finally, in Conclusion Section, we draw our conclusions.

Methodology
Let us consider the nonlinear evolution equation in the form where u = u(x, t) is an unknown function, F is a polynomial in u(x, t) and its derivatives in which highest order derivatives and nonlinear terms are involved and the subscripts indicate partial derivatives. In order to investigate solitary wave solutions of (1) by using the exp (-Ф(η))-expansion method, we have to perform the following important steps: Step 1. We combine the real variables x and t by a compound variable η where V is the celerity of the traveling wave. By means of traveling wave transformation (2), Eq. (1) switch into an ordinary differential equation (ODE) for u = u(η): where H is a polynomial of u and its derivatives and the superscripts refer to the ordinary derivatives with respect to η.
(3) and then we take into consideration the function exp(-Ф(η)). In consequence of this substitution, we obtain a polynomial in exp(-Ф(η)). We collect all the coefficients of identical power of exp(-Ф(η)) and equalize to zero delivers a system of algebraic equations whichever can be solved to find A N , ⋯⋯,V, λ, μ. The values of A N , ⋯⋯,V, λ, μ along with general solutions of Eq. (5) complete the determination of the solution of Eq. (1).

Application
In this section, we will use the exp(-Ф(η))-expansion method to construct the exact solutions and then the solitary wave solutions to the fourth order Boussinesq equation. Let us consider the equation The above model (11) was introduced by Boussinesq to illustrate the propagation of long waves in shallow water (Lai et al. 2008), where u(x, t) is the elevation of the free surface of the fluid, where the subscripts denoting partial derivatives. The equation also arises in many other physical applications, such as, nonlinear lattice waves, iron sound waves in plasma, and vibrations in a nonlinear string. It was also applied to the study of the percolation of water in porous subsurface strata.
Equation (11) possesses solitary waves, extract from traveling wave solutions and Boussinesq was the first who gave a scientific explanation of their existence. We utilize the traveling wave variable u(η) = u(x, t), η = x − Vt and this operation changes (11) to the following ODE: Integrating Eq. (12) twice with respect to η yields: where C is an integration constant to be determined. Balancing the highest order nonlinear term u 2 and linear term of the highest order u ′′ appearing in (13), yields N = 2. Therefore, the solution of Eq. (13) takes the form where A 0 , A 1 , A 2 are arbitrary constants such that A 2 ≠ 0. We substitute Eq. (14) into Eq. (13) and taking consideration Eq. (5), it generates a polynomial and then setting the coefficients of exp(-Ф(η)) to zero, yields where λ, μ and A 0 are arbitrary constants.
Substituting the values of V, A 0 , A 1 , A 2 into Eq. (14), yields where η ¼ x∓ q t: By using the solutions of Eq. (5) into Eq. (21), we obtain the succeeding traveling wave solutions of the fourth order Boussinesq equation: Type 1: When μ ≠ 0, λ 2 − 4μ > 0, t and E is an arbitrary constant.
Type 5: When μ = 0, λ = 0, and λ 2 − 4μ = 0, t and E is an arbitrary constant. Solitary wave solutions represent an important type of solutions for nonlinear partial differential equations (PDEs) as many nonlinear partial differential equations have been found to have a variety of solitary wave solutions. It is familiar that searching of exact solutions of nonlinear partial differential equations plays a significant role in the study of nonlinear physical phenomena. Exact traveling wave solutions are useful for verifying the accuracy and stability of popular numerical schemes such as the finite difference and finite element methods. The solitary wave solutions obtained in this article are encouraging, applicable, and could be helpful in analyzing long wave propagation on the surface of a fluid layer under the action of gravity, iron sound waves in plasma, and vibrations in a nonlinear string.

Graphical representation of the solutions
Solitary waves can be obtained from each traveling wave solution by setting particular values to its unknown parameters. By adjusting these parameters, one can get an internal localized mode. In this section, we have presented some graphs of solitary waves constructed by taking suitable values of the involved unknown parameters to visualize the underlying mechanism to the original physical phenomena. Using mathematical software Mathematica, two and three-dimensional plots of the obtained solutions have been shown in Figures 1, 2, 3, 4 and 5.
The obtained solutions of the fourth order Boussinesq equation incorporate four types of explicit solutions namely hyperbolic, trigonometric, exponential, and rational function solutions. From these explicit results we observe that solutions u 1 (η) and u 2 (η) are soliton and the rest of the three solutions are cuspon. The above solitary wave solutions might be useful in analyzing the propagation of long waves in shallow water, iron sound waves in plasma, and vibrations in a nonlinear string.

Conclusions
In this article, we have successfully formulated solitary waves solutions from the traveling wave solutions to the fourth order Boussinesq equation through the exp (-Ф(η))-expansion method. The procedure is simple, direct and constructive with the help of a computer algebra system. The method is quite efficient and practically well suited to be used in finding solitary wave solutions of NLEEs and the attained solutions demonstrated the competence of the exp(-Ф(η))-expansion method. We also observed that the method is straightforward and can be applied to many other nonlinear evolution equations.