An efficient technique to solve coupled–time fractional Boussinesq–Burger equation using fractional decomposition method

For this work, a novel numerical approach is proposed to obtain solution for the class of coupled time-fractional Boussinesq–Burger equations which is a nonlinear system. This system under consideration is endowed with Caputo time-fractional derivative. By means of the natural decomposition approach, approximate solutions of the proposed nonlinear fractional system are obtained where the exact solutions are presented in the classical case of fractional order at γ = 1 . Some numerical examples are given to support the theoretical framework and to point out the role and the effectiveness of the intended method. Our results clearly show the approximate analytical solutions eventually will converge quickly to the already known exact solutions. AMS Classification: 35A22, 35C05, 35C10, 35R11, 44A30.


Introduction
It was only the starts when Leibniz sent an amazing notice to l'Hospital in 1695 which led to the theory of fractional derivatives and integrals of arbitrary order such as one-half. Since then, scholars started developing theories of fractional derivative as theoretical field of pure Mathematics for three centuries. For the past three decades, researchers started bringing the attention about non-integer fractional integrals order along with its noninteger derivatives order. Moreover, they clarify the fact that they are very useful and more adequate than classical integer order in description the properties of real materials such as polymers, modelling mechanical and electrical properties of real materials, description of archaeological properties of rocks, theory of fractals. 1 Surveys of the history of fractional derivatives can be found in Miller and Ross, 2 Podlubny, 3 Al-Smadi et al., 4 and Eid et al. 5 We can summarize all the above describing the relation between these two concepts in one sentence: fractional derivative is the generalized form of classical integer derivative. In addition, the subject of fractional derivatives can be widely applied in many real life applications, such as; engineering, mathematical biology, quantum physics, fluid mechanics fields. [6][7][8][9][10][11] Due to the fast development of software programs such as Mat Lab, Mathematica and Maple, many new powerful analytical techniques have been proposed to find new and approximate solutions for fractional linear and nonlinear differential equations such as; the sub-equation method, 12 Exponential function method, 13 first integral method, 14 the expansion method, 15 fractional reproducing kernel method, [16][17][18] fractional Adomian decomposition method, 19 fractional homotopy perturbation, 20 fractional homotopy analysis, 21 fractional residual power series, [22][23][24][25] fractional Laplace decomposition, 26 fractional differential transform method [27][28][29][30] and other advanced numerical methods. [31][32][33][34] For the current work, an efficient method will be explored which we choose to call fractional decomposition method (FDM). It is a combination of the Adomian decomposition method and the natural transform method (NTM). [35][36][37] In fact, finding closed-form solutions for fractional partial differential equations is not an easy task for researchers due to the complexities of involving the fractional order. The goal of our study is to look for exact solution to coupled time-fractional fractional Boussinesq-Burger equation which is a nonlinear system using the FDM. Consider the non-linear coupled time-fractional Boussinesq-Burger equation of the form: along with the I.C's In the equation above, 0\g ł 1, x 2 ½a, b, t ø 0. Also, both functions g 1 (x), g 2 (x) are analytic, and w(x, t), v(x, t) are real-valued functions to be determined. Here D g t is for the Caputo time-fractional derivative of order g. We shall proceed in our paper as follow: Section 2 is intended for background of fractional calculus. In section 3 we introduce the NTM. Section 4 is devoted for the methodology of the FDM. In section 5, we give some examples of a fractional nonlinear system of differential equations. Section 6 is for discussion and conclusion of this paper.

Background of fractional calculus
In recent years, scientists showed some interest in nonlocal field theories and their interest really became more consistence. This development became more clear due to the expectation and the needs to use these theories so that they can treat problems in a more elegant and effective fashion. A particular group of non-local field theories plays an outstanding role and may be described with operators of fractional nature and is specified within the framework of fractional calculus. [4][5][6][7][8][9] Definition 2.1.If m 2 N, x.0 with g 2 C m À1 and m À 1\g ł m. Then, the fractional Caputo derivative of g is defined as: D g g(x) = J mÀg D m g(x) = 1 G(mÀg) Ð x 0 (x À y) mÀgÀ1 g (m) (y)dy: Note that, since the I.C and B.C will be included in the formulation of our applications, then the Caputo fractional derivative will be used through out the paper.

The natural transform method
Here, we refer the reader to Belgacem and Silambarasan 38 to view some of the background and history of the (NTM). Definition 3.1.Suppose H t ð Þ is the Heaviside function. Consider a real-valued function, where the natural transform N + is well-defined on the half plane s.bu for some b.0. Let g t ð Þ be continuous on R. Let K, b.0, then we define Definition 3.2.The natural transform of the function g(t) for t 2 R is given by 38 : where N g(t) ½ is the natural transform of g(t). Note that equation (3) can be written in the form 38 N g(t) ½ = R À (s, u) + R + (s, u): Suppose that g(t)H(t) is defined on R + , and s, u, t 2 R + , where g(t) 2 B. Then, we define the Ntransform as Consequently, if u = 1 we obtain the Laplace transform and if s = 1, we get the Sumudu transform. Also, one can find most of the properties of the N-transforms in Rawashdeh 30 and Belgacem and Silambarasan. 38 For

The methodology of FDM
In this section, we give the methodology of the FDM which also can be found in Rawashdeh 30 and Rawashdeh and Darweesh. 35 Theorem 4.1.If m 2 Z + , where m À 1 ł g\m and N + g(t) ½ = R + (s, u), then the natural transform of c D g g(t) is given by For the sake of explanation of the method algorithm, let us consider a nonlinear fractional system in the general form: along with the following initial conditions where and v x, t ð Þ respectively, L 1 and L 2 are the linear differential operators, g(x, t) and h(x, t) are the non-homogeneous terms and Now, by applying the N + to equation (5) and theorem (4.1), we have Thus, we apply N À1 to the above equation and we get In equation (8), we have K(x, t) and M(x, t) are counted for both the non-homogeneous and initial conditions. So, consider the series solutions v(x, t) = X ' From the above equation, then equation (8) becomes: Looking at both sides of equation (10), one can obtain If we proceed as before one can obtain this recursive relation Hence, our intended approximate solutions are as follows: v(x, t) = X '

Numerical examples
It has been demonstrated that the FDM deals efficiently with the fractional nonlinear system of differential equations when compared with the other methods that exist in literature. 35 This section provides some applications of nonlinear coupled time-fractional PDEs using the FDM, including coupled Boussinesq-Burger equations which is an application of dynamical system. 39 Example 1. Consider the coupled time-fractional Boussinesq-Burger equation: along with the I.C's According to the FDM algorithm, one can get: So equation (14) becomes Plug equation (13) into equation (14) to conclude Now, taking the N À1 of the above equation, we get Consequently, let us assume the infinite series solutions of the unknown functions w(x, t), v(x, t) and the nonlinear terms ww x , wv can be expanded in the forms From equations (17) and (18), we conclude that Now if we look at both sides of above equation, one can come up with: Thus, we have By continuing in the same way, we conclude: = w 0 (x, t) + w 1 (x, t) + w 2 (x, t) + ::: ð Þ + ::: À Á y 2g 4G 2g + 1 ð Þ + ::: : With the help of Taylor expansion, the intended approximate solutions for g = 1 are given by which coincides with the exact solutions of coupled time-fractional Boussinesq-Burger equations (12). In order for us, to see our proposed method (FDM) is reliable and efficient, the fractional behaviour of the approximate solutions fw n (x, t), v n (x, t)g of Example 5.1 above is discussed when n = 5 by utilizing the 2D and 3D graphs as follows: The 3D plots of w(x, t) and approximate solutions w n (x, t) together with absolute error are presented in Figure 1. While the behaviour of approximate solution w n (x, t) are plotted in Figure 2 for different values of fractional order g such that g 2 f0:6, 0:75, 0:9g. The 2D graphs of fractional levels of approximate solution w n (x, y) are presented in Figure 3. Similar graphical representations of approximate solution v n (x, t) are provided in Figures 4 to 6. From these graphs, it can be seen that the FDM approximations are in closed agreement with each other for various values of g and with exact solution. Further, we can illustrate the validity of the FDM by looking at both, the exact and approximate solutions for different values of fractional order. In Tables 1 and 2, some numerical results of the approximate solutions w(x, t) and v(x, t) are considered at various values of x and t such that x = fÀ5, 0, 5g and t = f0:01, 0:03, 0:05g. Example 2. Consider the fractional system below:  D a x t ð Þ = 2y 2 t ð Þ, 0\a ł 1, along with the I.C's Employ the fractional N-transform algorithm on equation (20), one can conclude: Figure 3. Fractional curves of w 5 (x, t) for t = 1, x 2 ( À 2, 2) (left) and x = 2, t 2 ½0, 1 (right).   So equation (22) becomes Combine equations (23) and (21), to accomplish Now, we take the N À1 of the above equation to obtain Consequently, let us assume the infinite series solutions for x(t), y(t) and z(t) which can be written in the forms and the nonlinear terms y 2 t ð Þ and y t ð Þz t ð Þ as If we combine equations (26) and (27), then equation (25) implies Now looking at equation (28), one can come up with the recursive relation x 0 0 ð Þ = 0, y 0 0 ð Þ = 1, z 0 0 ð Þ = 1: Thus, Subsequently, the 3rd approximate solutions can be obtained as Proceeding in this way, one comes up with these approximate solutions ð Þ + :::: For graphical representation, 2D plots of approximate solutions x(t), y(t) and z(t) are respectively presented in

Conclusion
We successfully employed the The fractional decomposition method and we obtained analytical approximate and exact solutions for two time-fractional order nonlinear systems. We were being able to find exact solutions to nonlinear system of coupled time-fractional Boussinesq-Burger equation. To the best of our knowledge, we are the first to find such exact solutions for   the proposed systems. Since the exact solutions of most FDE's cannot be found easily, then analytical and numerical methods like (FDM) can be used more often. In all cases, the (FDM) provided us with exact solutions in the case when g = 1. The results showed that (FDM) is simple and easy mathematical technique to accomplish exact and numerical solutions of nonlinear time-fractional equations. Finally, one can conclude the FDM can be employ to investigate and study numerous applications of fractional differential equations which usually shows up in many areas of Physics and engineering.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.