Laplace Decomposition Method for Time-Fractional Fornberg-Whitham Type Equations

In this study, the closed form of series solutions of the original and modified nonlinear time-fractional Fornberg-Whitham equations are derived by means of the Laplace decomposition method (LDM). The fractional order derivatives are expressed in the sense of Caputo. For the specific choice of parame-ters, the obtained solutions are compared with the exact solutions to validate the accuracy of this method. Numerical solutions are represented graphically which illustrate the behavior of the solutions. Further, the computations express that the above method is straightforward, and it desires the smaller size of computation.


Introduction
The fractional order nonlinear differential equations (NDE) have been the focus of many studies due to their numerous applications in the areas of mathematical physics and engineering [1]. In recent years, it has been a great interest to find the analytical methods for solving fractional NDE. Exact solutions of most of the fractional NDE cannot be found easily, consequently, analytical and approximate methods have been used. Some of the recent approximate methods for solving NDE include the Variational iteration method (VIM) [2], Differential transformation method (DTM) [3], Homotopy-perturbation method (HPM) [4], Adomian decomposition method (ADM) [5] [6] [7] and Homotopy analysis method (HAM) [8]. tion is one such type of equation which was first introduced for illustrating the qualitative behaviour of wave breaking [9]. Recently, various types of methods are successfully applied to find extract travelling wave solutions and numerical solutions to the nonlinear Fornberg-Whitham equation [10] [11] [12].

Preliminaries
We provided some definitions and basic properties of fractional calculus theory and then we give the basics of Laplace transformation which will be applied later in this work. Definition 2.1 Let a function , 1 f C µ µ ∈ ≥ − and 0 α ≥ , then the Riemann-Liouville fractional integral of order α is defined by Some basic properties of the operator J α are as following [21]: Caputo and Mainardi [22] modified the fractional differentiation operator which overcomes the discrepancy of Riemann-Liouville derivative.
Definition 2.2 The Caputo fractional derivative of order 0 α > is defined by . Also the derivative satisfies the following properties: whenever the improper integral converges. Some properties of the Laplace transform are given below:

The Laplace Decomposition Method
In order to explain the central idea of this method, first we consider a general fractional nonlinear PDE with the initial condition of the following form: Journal of Applied Mathematics and Physics , u x t , R is a linear operator, N is a nonlinear operator and g is a non-homogeneous term.
We take the Laplace transform on both sides of Equation (8) and then using the differentiation property of Laplace transform, we get, Now, taking the inverse Laplace transform to Equation (9) where ( ) , G x t is source term and give initial condition.
The series solution of the Adomian decomposition method is given by, with the nonlinear term ( ) , Nu x t can be decomposed as Adomian polynomials n A can be derived by particular algorithms constructed by Adomian [24] or alternatively by Wazwaz [25]. The formula for n A is given Substituting (11) and (12) into Equation (10), we get, From the equation (14) we get, From (15), the series solution of (8) can be found In generally, the above series is convergence rapidly [   can be calculated using the procedure stated in [24] as follows: and so on. Using the above recursive relationships with Adomian polynomials the decomposition series (16) can be found as follows:

Implementation of LDM
Mathanaranjan, D. Vijayakumar and so on. In the following two subsection, we utilize the above results for two initial value problems associated with the original and modified time-fractional Fornberg-Whitham equations in order to find the series solutions.

Solution of the Original Time-Fractional Fornberg-Whitham Equation
For 1 n = , the Equation (17) To verify the accuracy of the LDM, we compare the numerical values of approximate solutions with the exact solutions in Table 1. These numerical values show that LDM yields a good approximate solution, even for the lowest order approximate solution. The accuracy of solution can be improved by taking higher-order approximate solutions.
In Figure 1, the approximate solution and the exact solution are plotted. its leads to get the behavior of the approximate solutions are the same behavior with the exact solution. Consequently, the approximate series solutions are rabidly convergence as the exact solution. The 3D and 2D plots for different values of α are presented in Figure 2 and Figure 3.

Solution of Modified Time-Fractional Fornberg-Whitham Equation
For 2 n = , the Equation (17)       Similarly as in the previous case, the approximate solutions are the same behavior with the exact solution. This is illustrated by Figure 4. The 3D and 2D plots for some different values of α are presented in Figure 5 and Figure 6.

Conclusion
In this study, the Laplace-Adomian decomposition method has been successfully applied to find a closed form of series solution of the original and modified time-fractional Fornberg-Whitham equations. For some specific choice of parameters, the approximate solutions are compared with the exact solutions. Further, the behaviours of the solutions are presented by 3D and 2D graphs. The graphs showed that the behaviors of the approximate solutions are the same behaviors of the exact solutions. As a result, we realize that the approximate solutions are rapidly convergent series as the exact solutions. Further, the computations show that the described method is easy to apply and it needs the smaller size of computation than the existing methods. Therefore, it is obvious that the LDM is very powerful, efficient and accurate for solving the nonlinear fractional differential equations.