An analytical treatment to fractional Fornberg–Whitham equation

In this paper, an analytical technique, namely the new iterative method (NIM), is applied to obtain an approximate analytical solution of the fractional Fornberg–Whitham equation. The obtained approximate solutions are compared with the exact or existing numerical results in the literature to verify the applicability, efficiency, and accuracy of the method.


Introduction
In recent years, the fractional calculus used in many phenomena in engineering, physics, biology, fluid mechanics, and other sciences [1][2][3][4][5][6][7][8] can be described very successfully by models using mathematical tools from fractional calculus. Fractional derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes [9]. The fractional derivative has been occurring in many physical and engineering problems such as frequency-dependent damping behavior materials, signal processing and system identification, diffusion and reaction processes, creeping and relaxation for viscoelastic materials.
The new iterative method (NIM), proposed first by Gejji and Jafari [10], has proven useful for solving a variety of nonlinear equations such as algebraic equations, integral equations, ordinary and partial differential equations of integer, and fractional order and systems of equations as well. The NIM is simple to understand and easy to implement using computer packages and yield better results than the existing Adomain decomposition method [11], homotopy perturbation method [12], and variational iteration method [13].
In the present paper, we have to solve the nonlinear time-fractional Fornberg-Whitham equation by the NIM. This equation can be written in operator form as: with the initial condition where u(x, t) is the fluid velocity, a is constant and lies in the interval (0, 1), t is the time and x is the spatial coordinate. Subscripts denote the partial differentiation unless stated otherwise. Fornberg and Whitham obtained a peaked solution of the form u(x, t) = Ae -1/2(|x-4t/3|) where A is an arbitrary constant.

Preliminaries and notations
In this section, we set up notation and review some basic definitions from fractional calculus [14,15].

Definition 2.3
Let f 2 C a and a C -1, then the (leftsided) Riemann-Liouville integral of order l, l [ 0 is given by: Note that

Basic idea of new iterative method (NIM)
To describe the idea of the NIM, consider the following general functional equation [10,[16][17][18][19][20]: where N is a nonlinear operator from a Banach space B ? B and f is a known function. We are looking for a solution u of (3) having the series form: The nonlinear operator N can be decomposed as follows: From Eqs. (4) and (5), Eq. (3) is equivalent to: We define the recurrence relation: Then: If N is a contraction, i.e., then: and the series P i=0 ? u i absolutely and uniformly converges to a solution of (3) [21], which is unique, in view of the Banach fixed point theorem [22]. The k term approximate solution of (3) and (4) is given by P i=0 k-1 u i .

Convergence analysis of the new iterative method (NIM)
Now, we introduce the condition of convergence of the NIM, which is proposed by Daftardar-Gejji and Jafari in (2006) [10], also called (DJM) [23]. From (5), the nonlinear operator N is decomposed as follows: Let G 0 = N(u 0 ) and Then N(u) = P i=0 ? G i . Set: Then: is a solution of the general functional Eq. (3). Also, the recurrence relation (7) becomes Using Taylor series expansion for G i , i = 1, 2, …, n, we have In general: In the following theorem, we state and prove the condition of convergence of the method.
for any n and for some real L [ 0 and u i k k M\ 1 e ; i ¼ 1; 2; . . .; then the series P n=0 ? G n is absolutely convergent, and moreover, G n k k LM n e nÀ1 ðe À 1Þ; n ¼ 1; 2; . . .: Proof In view of (18) h Thus, the series P n=1 ? kG n k is dominated by the convergent series LM(e -1) P n=1 ? (Me) n-1 , where M\ 1 e . Hence, P n=0 ? G n is absolutely convergent, due to the comparison test. For more details, see [23].
Reliable algorithm of new iterative method (NIM) for solving the Linear and Nonlinear fractional partial differential equations After the above presentation of the NIM, we introduce a reliable algorithm for solving nonlinear fractional PDEs using the NIM. Consider the following nonlinear fractional PDE of arbitrary order: with the initial conditions where A is a nonlinear function of u and qu (partial derivatives of u with respect to x and t) and B is the source function.
In view of the integral operators, the initial value problem (22) is equivalent to the following integral equation and where I t n is an integral operator of n fold. We get the solution of (24) by employing the algorithm (7).

Solution of the problem
We first consider the following time-fractional Fornberg-Whitham equation [24,25]: with the initial condition: Then, the exact solution is given by: Note that Eq. (27) is equivalent to the integral equation and NðuÞ ¼ I a t D a t u À D xxt u þ D x uþ Â uD xxx u À uD x u þ 3D x uD xx u, using (7) we get . . .
Therefore, the NIM series solution is: Numerical results and discussion In this section, we calculate numerical results of the displacement u(x, t) for different time-fractional Brownian motions a = 2/3, 3/4, 1 and for various values of t and x. The numerical results for the approximate solution (30) obtained using NIM and the exact solution for various values of t, x, and a are shown by Figs. 1, 2, 3, 4 and those for different values of t and a at x = 1 are depicted in Fig. 5. It is observed from Figs. 1, 2, 3, 4 that u(x, t) increases with the increase in both x and t for a = 2/3, 3/4 and a = 1. Figure 4 clearly show that, when a = 1, the approximate solution (30) obtained by the present method is very near to the exact solution. It is also seen from Fig. 5 that as the value of a increases, the displacement u(x, t) increases but afterward its nature is opposite. Finally, we remark that the approximate solution (30) is in full agreement with the results obtained homotopy perturbation method [24] and homotopy perturbation transform method [25].
In Table 1, we compute the absolute errors for differences between the exact solution (29) and the approximate solution (30) obtained by the NIM at some points.

Conclusion
In this paper, the new iterative method (NIM) has been applied for approximating the solution for the nonlinear fractional Fornberg-Whitham equation. The accuracy of the NIM for solving nonlinear fractional Fornberg-Whitham equation is good compared to the literature; however, it has the advantage of reducing the computations complexity presented in other perturbation techniques. In fact, in NIM, nonlinear problems are solved without using Adomian's polynomials or He's polynomials that appear in the decomposition methods. The numerical results show that the proposed method is reliable and efficient technique in finding approximate solutions for nonlinear differential equations.