**American Journal of Numerical Analysis**

## Numerical Simulation of Time-fractional Fourth Order Differential Equations via Homotopy Analysis Fractional Sumudu Transform Method

**Rishi Kumar Pandey**^{1}, **Hradyesh Kumar Mishra**^{1,}

^{1}Department of Mathematics, Jaypee University of Engineering and Technology, Guna, INDIA

### Abstract

The work provides an incipient analytical technique called the Homotopy Analysis Fractional Sumudu Transform Method (HAFSTM) for solving time-fractional fourth order differential equations with variable coefficients. The HAFSTM is the cumulation of the homotopy analysis method (HAM) and sumudu transform method (STM). The numerical simulation of the proposed method has the sundry applications, it can solve linear and nonlinear boundary value quandaries without utilizing Adomian polynomial, and He’s polynomial, which can be considered a clear advantage of this incipient algorithm. The solutions obtained by proposing technique are very lucid and less computationally implementable.

**Keywords:** Homotopy Analysis Method, Fractional Sumudu Transform Method, Fractional Partial Differential equation, variable coefficients, boundary value problem

**Copyright**© 2015 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Rishi Kumar Pandey, Hradyesh Kumar Mishra. Numerical Simulation of Time-fractional Fourth Order Differential Equations via Homotopy Analysis Fractional Sumudu Transform Method.
*American Journal of Numerical Analysis*. Vol. 3, No. 3, 2015, pp 52-64. http://pubs.sciepub.com/ajna/3/3/1

- Pandey, Rishi Kumar, and Hradyesh Kumar Mishra. "Numerical Simulation of Time-fractional Fourth Order Differential Equations via Homotopy Analysis Fractional Sumudu Transform Method."
*American Journal of Numerical Analysis*3.3 (2015): 52-64.

- Pandey, R. K. , & Mishra, H. K. (2015). Numerical Simulation of Time-fractional Fourth Order Differential Equations via Homotopy Analysis Fractional Sumudu Transform Method.
*American Journal of Numerical Analysis*,*3*(3), 52-64.

- Pandey, Rishi Kumar, and Hradyesh Kumar Mishra. "Numerical Simulation of Time-fractional Fourth Order Differential Equations via Homotopy Analysis Fractional Sumudu Transform Method."
*American Journal of Numerical Analysis*3, no. 3 (2015): 52-64.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### At a glance: Figures

### 1. Introduction

Sumudu integral transform was introduced by Watugala ^{[1]} to facilitate the process of solving differential and integral equations in the time domain, and for the utilization in sundry applications of system engineering and applied physics. Albeit the mathematical properties of the Sumudu transform have been explored in some details [2-9]^{[2]}, to the best of our cognizance, no systematic derivation of the Sumudu transform is available in the open literature.

The fractional calculus deals with arbitrary orders of derivatives and integrals of applied mathematics. In the last decade, the fractional calculus has found applications in numerous ostensibly diverse fields of science and engineering. Fractional differential equations are increasingly used to model quandaries in fluid mechanics, acoustics, biology, electromagnetism, diffusion, signal processing, and many other physical processes.

In sundry fields of science and engineering, it is consequential to obtain exact or numerical simulation of the nonlinear partial differential equations. In science and engineering, to find the exact and numerical solution of nonlinear equations is still challenging, therefore it’s required lucid methods for finding the exact and approximate solutions. Sundry potent mathematical methods such as Adomian decomposition method (ADM) ^{[10, 11, 12, 13, 14]}, homotopy perturbation method (HPM) [15-20]^{[15]}, homotopy analysis method (HAM) ^{[21, 22, 23, 24, 25]}, variational iteration method (VIM) [26-32]^{[26]}, Laplace decomposition method (LDM) ^{[33, 34, 35]}, homotopy perturbation transform method (HPTM) ^{[36]}, homotopy perturbation sumudu transform method (HPSTM) ^{[37]} and homotopy analysis transform method (HATM) ^{[38, 39, 40]} have been proposed to obtain exact and approximate analytical solutions of nonlinear equations.

Inspired by all these perpetual research, we propose HAFSTM for the solution of fourth order differential equations with variable coefficient.

### 2. Sumudu Transform

In early 90’s, Watugala ^{[1]} introduced an incipient integral transform, designated the sumudu transform and applied it to the solution of mundane differential equation in control engineering quandaries. The sumudu transform is defined over the set of functions

by the following formula

Weerakoon ^{[2, 3, 4]} was established the properties of mentioned transform. Further fundamental properties of this transform were also established by Asiru ^{[4]}. In the same way, this transform was applied to the unidimensional neutron convey equation in ^{[5]} by Kadem. In fact, it was shown that there is a vigorous relationship between Sumudu and other integral transforms refer to Kılıҫman et al. ^{[6]}. In particular the relation between Sumudu transform and Laplace transforms was proved in Kılıҫman and Gadain ^{[7]}. Further, in Eltayeb et al. ^{[9]}, the Sumudu transform was elongated to the distributions and some of their properties were furthermore studied in Kılıc¸man and Eltayeb ^{[9]}. Recently, this transform is applied to solve the system of differential equations ^{[41, 42, 43]}.

### 3. Basic Definition of Fractional Calculus

**Definition 3.1 **A real function is said to be in the space if there exists a real number such that where and it is said to be in the space iff .

**Definition 3.2 **The Riemann Liouville Fractional integral operator of order of a function and is defined as ^{[44, 45]}

and

For the Riemann – Liouville fractional integral, we have

**D****efinition 3****.****3 **The fractional derivative of in the Caputo sense is defined as ^{[46]}

where

For the Riemann – Liouville fractional derivative, we obtain the following relation.

where

### 4. Solution by Homotopy Analysis Fractional Sumudu Transform Method

To illustrate the rudimental conception of the HAFSTM for the fractional partial differential equation, we consider the following fractional partial differential equation as

(4.1) |

where is the linear operation in is the general nonlinear operation in and is a continuous function.

All initial and boundary conditions are ignored for minimalism, which can be treated in a homogeneous way now applying the Sumudu transform first on both sides of the equation (4.1), we get

(4.2) |

Using the differentiation property of the Sumudu transform

(4.3) |

we define nonlinear operator as

(4.4) |

where be an embedding parameter and is a real function of and we construct a homotopy as follow:

(4.5) |

where is a nonzero auxiliary parameter and An auxiliary function is an initial guess of and is an unknown function. It is important that one has great freedom to choose auxiliary parameter in HAFSTM. Obviously, when and it holds

(4.6) |

Thus, as increases from 0 to 1,the solution varies from initial guess to the solution Now, expanding on Taylor’s series with respect to we get

(4.7) |

where

(4.8) |

The convergence of the series solution (4.7) is controlled by If the auxiliary linear operator, the initial guess, the auxiliary parameter and the auxiliary function are properly chosen, the series (4.7) converges at Hence we obtain

(4.9) |

which must be one of the solutions of original nonlinear equations. The above expression provides us with a relationship between the initial guess and the exact solution by means of the terms which are still to be determined.

Define the vectors

(4.10) |

Differentiating the zero order deformation Eq. (4.5) times with respect to embedding parameter and then setting and finally dividing them by we obtain the order deformation equation as follows:

(4.11) |

Operating the inverse Sumudu transform of both sides, we get

(4.12) |

where

(4.13) |

and

In this way, it is easy to obtain for at order, we have

(4.14) |

where we obtain an accurate approximation of the original equation (4.1).

### 5. Numerical Examples

In this section, we apply the HAFSTM developed in Section 4 to solve one and two dimensional initial boundary value quandaries with variable coefficients. The methods may additionally be applicable for higher dimensional spaces. Numerical results reveal that the HAFSTM is facile to implement and reduce the computational work to a tangible level while still maintaining a higher caliber of precision. All the results for the following three applications are calculated by utilizing the symbolic calculus software MATHEMATICA 8.0.

**Example 5****.1****: **Consider the following problem of one-dimensional time-fractional fourth-order PDE ^{[47]}

(5.1) |

subject to the initial and boundary conditions

(5.2) |

where the function is defined as

(5.3) |

We start with initial condition

(5.4) |

Operating the Sumudu transform of both sides in Eq. (5.1) and after using the differentiation property of Sumudu transform for fractional derivative, we get

(5.5) |

or

(5.6) |

The nonlinear operator is

(5.7) |

and thus

(5.8) |

The order deformation equation is given by

Applying the inverse Sumudu transform, we have

(5.9) |

On solving above equation from we get

(5.10) |

(5.11) |

(5.12) |

and so on, in the same manner the rest of the components of the series can be obtained.

Finally, the solution of (5.1) is given as

(5.13) |

According to Liao ^{[24]}, the accuracy and convergence of the HAM series solution depends on the careful selection of the auxiliary parameter Here we choose then

(5.14) |

For the above expressions are exactly the same as those given by the Decomposition method by N. A. Khan ^{[47]}.

Nevertheless, most of the results given by the ADM, LDM, HPM and HPTM converge to the corresponding numerical solutions in a rather small region. Except, different from these methods, the HAFSTM provides us through a simple way to alter and manage the convergence region of solution series by choosing a proper value for the auxiliary parameter . So, the valid region for where the series converge is the horizontal segment of each curve. When we choose then clearly, we can conclude that the obtained solution converges to the exact solution obtained by Wazwaz ^{[45, 46]}.

**Example 5.2:**** **Consider the following problems of one-dimensional time-fractional fourth-order PDE ^{[47]}

(5.15) |

subject to the initial and boundary conditions

(5.16) |

where the function

We start with initial condition

Operating the Sumudu transform on both sides in (5.15) and after using the differentiation property of Sumudu transform for fractional derivative, we get

(5.17) |

Applying the initial and boundary conditions of (5.15) from (5.16) to (5.17), we obtain

The nonlinear operator is

(5.18) |

and thus

(5.19) |

The - order deformation equation

Applying the inverse Sumudu transform,

(5.20) |

On solving above equation from we get

(5.21) |

(5.22) |

(5.23) |

and so on, in the same manner the rest of the components of the series can be obtained.

Finally, the solution of Eq. (5.15) is given as

(5.24) |

According to Liao ^{[24]}, the accuracy and convergence of the HAM series solution depends on the careful selection of the auxiliary parameter here, we choose then

(5.25) |

For the above expressions are exactly the same as those given by the Decomposition method by N. A. Khan ^{[47]}.

When we choose then clearly, we can conclude that the obtained solution converge to the exact solution obtained by Wazwaz ^{[48, 49]}, and Biazar and Ghavini ^{[50]}.

**Example 5.3: **Consider the following problem of two-dimensional time-fractional fourth-order PDE ^{[47]}

(5.26) |

subject to the initial and boundary conditions:

(5.27) |

we start with initial condition

Operating the Sumudu transform on both sides in (5.26) and after using the differentiation property of Sumudu transform for fractional derivative, we get

Applying the initial and boundary conditions of equation (5.26) from equation (5.27), we obtain

(5.28) |

The nonlinear operator is

(5.29) |

and thus

(5.30) |

The - order deformation equation

Applying the inverse Sumudu transform,

(5.31) |

On solving above equation from we get

(5.32) |

(5.33) |

(5.34) |

and so on, in the same manner the rest of the components of the series can be obtained.

Finally, the solution of Eq. (5.26) is given as

(5.35) |

The accuracy and convergence of the HAM series solution depend on the careful selection of the auxiliary parameter here, we choose them

(5.36) |

For the above expressions are exactly the same as those given by the Decomposition method by N. A. Khan ^{[47]}.

### 6. Numerical Results and Discussion

**Example 1**

**Figure**

**1**

**.**Plot of w.r. to ћ at and

**Figure**

**2**

**.**Plot of w.r.t and at

**Figure**

**3**

**.**Plot of w.r.t and at

**Figure**

**4**

**.**Plot of w.r.t and at

**Figure**

**5**

**.**

**Plot of Exact Solution of w.r.t and**

**Example 2 **

**Figure**

**6**

**.**Plot of w.r. to ћ at and

**Figure**

**7**

**.**Plot of w.r.t and at

**Figure**

**8**

**.**Plot of w.r.t and at

**Figure**

**9**

**.**Plot of w.r.t and at

**Figure**

**10**

**.**Plot of Exact Solution of w.r.t and

**Example 3 **

**Figure**

**11**

**.**Plot of w.r. to ћ at and

**Figure**

**12**

**.**

**Plot of w.r.t and at**

**Figure**

**13**

**.**Plot of w.r.t and at

**Figure**

**14**

**.**Plot of w.r.t and at

**Figure**

**15**

**.**Plot of Exact Solution of w.r.t and

Figure 1 Shows that the curve between approximate solution and convergence control parameter for the different values of fractional order viz. The convergence lies between the range at and for example 1.

Figure 2 shows that the three dimensional plot between for independent variables and at Similarly Figure 3 and Figure 4 show the corresponding slight changes for different fractional Brownian motions of respectively.

Figure 5 is plotted for the exact solution of which is equal to the Figure 4.

In the subsequent manner the plot of example 2, curve for lies between in Fig. 6. and Figure 7-Figure 9 are shown the plot of approximate solution correspond to two independent variables and verses Figure 10 is plotted for exact solution which shows the same plot as Figure 9.

Figure 11-Figure 15 show the evaluation results of the approximate analytical solution for the Example 3. These also figures show the behavior of the approximate solution obtained by the proposed method for different fractional Brownian motions and the convergence region for convergence control parameter and approximate solutions.

Solution at integral value at in all above mentioned plots is shown the same as obtained by N.A. Khan et. al ^{[44]} , Wazwaz ^{[45, 46]} and Biazar and Gazvini ^{[47]}.

### 7. Conclusion

In science and engineering, the incipient modification of HAFSTM is potent implement to probe the solution of sundry linear and nonlinear quandaries. The main aim of this article is to provide the approximate solution and additionally analytic approximation utilizing the proposed method for fourth order boundary value quandaries. The analytical results have been given in terms of a potency series with facilely computed terms. The method surmounts the arduousness in different methods because it is efficient and lucid. Three examples were investigated to demonstrate the facileness and multifariousness of our incipient approach. The illustrative examples show that the method is facile to utilize and is an efficacious implement to solve fractional partial differential equations numerically.

### Acknowledgements

The authors are grateful to the referees for their valuable suggestions and comments for the improvement of the paper. The First author acknowledges the financial support provided by the JUET, Guna, India as a teaching assistantship.

### References

[1] | G.K. Watugala, Sumudu transform—a new integral transform to solve differential equations and control engineering problems,” Mathematical Engineering in Industry, vol. 6(4) (1998) 319-329. | ||

In article | |||

[2] | S. Weerakoon, Application of Sumudu transform to partial differential equations, International Journal of Mathematical Education in Science and Technology, vol. 25(2) (1994) 277-283. | ||

In article | View Article | ||

[3] | S. Weerakoon, Complex inversion formula for Sumudu transform, International Journal of Mathematical Education in Science and Technology, vol. 29(4) (1998) 618-621. | ||

In article | |||

[4] | M. A. Asiru, Further properties of the Sumudu transform and its applications, International Journal of Mathematical Education in Science and Technology, vol. 33(3) (2002) 441-449. | ||

In article | View Article | ||

[5] | A. Kadem, Solving the one-dimensional neutron transport equation using Chebyshev polynomials and the Sumudu transform, Analele Universitatii din Oradea, vol. 12 (2005)153-171. | ||

In article | |||

[6] | A. Kılıc¸man, H. Eltayeb, and K. A. M. Atan, A note on the comparison between Laplace and Sumudu transforms, Iranian Mathematical Society, vol. 37(1) (2011)131-141. | ||

In article | |||

[7] | A. Kılıc¸man and H. E. Gadain, On the applications of Laplace and Sumudu transforms, Journal of the Franklin Institute, vol. 347(5) (2010) 848-862. | ||

In article | View Article | ||

[8] | H. Eltayeb, A. Kılıc¸man, and B. Fisher, A new integral transform and associated distributions, Integral Transforms and Special Functions, vol. 21(5-6) (2010) 367-379. | ||

In article | View Article | ||

[9] | A. Kılıc¸man and H. Eltayeb, “A note on integral transforms and partial differential equations,” Applied Mathematical Sciences, vol. 4(1-4) (2010) 109-118. | ||

In article | |||

[10] | G. Adomian, A review of the decomposition method in applied mathematics, J. Math. Anal. Appl., vol. 135(1988) 501-544. | ||

In article | View Article | ||

[11] | N.T. Shawagfeh, Analytical approximate solutions for nonlinear fractional differential equations, Appl. Math. Comput., vol. 131 (2002)517-529. | ||

In article | View Article | ||

[12] | S.S. Ray, R.K. Bera, An approximate solution of a nonlinear fractional differential equation by Adomian decomposition method, Appl. Math. Comput., vol. 167(2005) 561-571. | ||

In article | View Article | ||

[13] | H. Jafari, V. Daftardar-Gejji, Revised Adomian decomposition method for solving systems of ordinary and fractional differential equations, Appl. Math. Comput., vol. 181(1) (2006)598-608. | ||

In article | View Article | ||

[14] | H. Jafari, V. Daftardar-Gejji, Positive solutions of nonlinear fractional boundary value problems using Adomian decomposition method, Appl. Math. Comput., vol. 180 (2) (2006) 700-706. | ||

In article | View Article | ||

[15] | A. Golbabai, M. Javidi, Application of homotopy perturbation method for solving eighth-order boundary value problems, Appl. Math. Comput., vol. 191(1) (2007) 334-346. | ||

In article | View Article | ||

[16] | J. H. He, A coupling method of a homotopy technique and a perturbation technique for nonlinear problems, Int. J. Non-Linear Mech., vol. 35(1)(2000) 37-43. | ||

In article | View Article | ||

[17] | J. H. He, The homotopy perturbation method for non-linear oscillators with discontinuities, Appl. Math. Comput., vol. 151(1) (2004) 287-292. | ||

In article | View Article | ||

[18] | J. H. He, Application of homotopy perturbation method to nonlinear wave equations, Chaos Solitons Fractals, vol. 26(3) (2005) 695-700. | ||

In article | View Article | ||

[19] | J. H. He, Asymptotology by homotopy perturbation method, Appl. Math. Comput., vol. 156 (3) (2004) 591-596. | ||

In article | View Article | ||

[20] | J. H. He, Homotopy perturbation method for solving boundary problems, Phys. Lett. A, vol. 350(12) (2006) 87-88. | ||

In article | View Article | ||

[21] | Sh. S. Behzadi, Iterative methods for solving nonlinear Fokker-Plank equation, Int. J. Industrial Mathematics, vol. 3 (2011) 143-156. | ||

In article | |||

[22] | Sh. S. Behzadi, Numerical solution of Sawada-Kotera equation by using iterative methods, Int. J. Industrial Mathematics, vol. 4 (2012)269-288. | ||

In article | |||

[23] | M. Ghanbari, Approximate analytical solutions of fuzzy linear Fredholm integral equations by HAM, Int. J. Industrial Mathematics, vol. 4 (2012)53-67. | ||

In article | |||

[24] | S.J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems, Ph.D. Thesis, Shanghai Jiao Tong University, 1992. | ||

In article | |||

[25] | M. Zurigat, S. Momani, A. Alawneh, Analytical approximate solutions of systems of fractional algebraic–differential equations by homotopy analysis method, Comput. Math. Appl., vol. 59 (3) (2010)1227-1235. | ||

In article | View Article | ||

[26] | G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Acad. Publ., Boston, 1994. | ||

In article | |||

[27] | Sh. S. Behzadi, Iterative methods for solving nonlinear Fokker-Plank equation, Int. J. Industrial Mathematics, vol. 3 (2011) 143-156. | ||

In article | |||

[28] | J.H. He, Variational iteration method- a kind of nonlinear analytical technique: some examples, International Journal of Nonlinear Mechanics, vol. 34 (1999) 699-708. | ||

In article | View Article | ||

[29] | J.H. He, X.H. Wu, Variational iteration method: new development and applications, Computers & Mathematics with Applications, vol. 54 (2007) 881-894. | ||

In article | View Article | ||

[30] | J.H. He, G.C.Wu, F. Austin, The variational iteration method which should be followed, Nonlinear Science Letters A, vol. 1 (2009) 1-30. | ||

In article | |||

[31] | E. Hesameddini, H. Latifizadeh, An optimal choice of initial solutions in the homotopy perturbation method, International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10 (2009) 1389-1398. | ||

In article | View Article | ||

[32] | E. Hesameddini, H. Latifizadeh, Reconstruction of variational iteration algorithms using the Laplace transform, International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10 (2009) 1377-1382. | ||

In article | View Article | ||

[33] | Yasir Khan, An effective modification of the Laplace decomposition method for nonlinear equations, International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10 (2009) 1373-1376. | ||

In article | View Article | ||

[34] | S.A. Khuri, A Laplace decomposition algorithm applied to a class of nonlinear differential equations, Journal of Applied Mathematics, vol. 1 (2001) 141-155. | ||

In article | View Article | ||

[35] | AR. Vahidi, Gh. A. Cordshooli, On the Laplace transform decomposition algorithm for solving nonlinear differential equations, Int. J. Industrial Mathematics, vol. 3(2011) 17-23. | ||

In article | |||

[36] | Y. Khan, Qingbiao. Wu, Homotopy perturbation transform method for nonlinear equations using He’s polynomials, Computers & Mathematics with Applications, vol.61(8) (2011) 1963-1967. | ||

In article | View Article | ||

[37] | J. Singh, D. Kumar, Sushila, Homotopy perturbation sumudu transform method for nonlinear equations, Adv. Theor. Appl. Mech., vol. 4(2011) 165-175. | ||

In article | |||

[38] | S. Kumar, A new analytical modelling for fractional telegraph equation via Laplace transform, Applied Mathematical Modelling vol. 38(13) (2014) 3154-3163. | ||

In article | View Article | ||

[39] | M. S. Mohamed, F. Al-malki, M. Al-humyani, Homotopy Analysis Transform Method for Time-Space Fractional Gas Dynamics Equation, Gen. Math. Notes, vol. 24(1) (2014) 1-16. | ||

In article | |||

[40] | S. Kumar, M. M. Rashidi, New analytical method for gas dynamics equation arising in shock fronts, Computer Physics Communications, vol. 185 (7) (2014)1947-1954. | ||

In article | View Article | ||

[41] | A. Kılıc¸man, H. Eltayeb, and R. P. Agarwal, “On Sumudu transform and system of differential equations,” Abstract and Applied Analysis, Article ID598702, 11 pages, 2010. | ||

In article | |||

[42] | J. Zhang, “A Sumudu based algorithm for solving differential equations,” Academy of Sciences of Moldova, vol. 15(3) (2007) 303-313. | ||

In article | |||

[43] | V. B. L. Chaurasia and J. Singh, “Application of Sumudu transform in Sch¨odinger equation occurring in quantum mechanics,” Applied Mathematical Sciences, vol. 4(57-60), (2010)2843-2850. | ||

In article | |||

[44] | N. A. Khan, N.U. Khan, M. Ayaz,A. Mahmood, N. Fatima, Numerical study of time- fractional fourth-order differential equations with variable coefficients, Journal of King Saud University – Science, vol. 23(1)(2011) 91-98. | ||

In article | |||

[45] | Y. Luchko, R. Gorenflo, An operational method for solving fractional differential equations with the Caputo derivatives, Acta Math. Vietnamica, 24 (1999) 207-233. | ||

In article | |||

[46] | O. L. Moustafa, On the Cauchy problem for some fractional order partial differential equations, Chaos Solitons& Fractals, 18 (2003) 135–140. | ||

In article | View Article | ||

[47] | I. Podlubny, Fractional Differential Equations, Academic, New York, 1999. | ||

In article | |||

[48] | A.M. Wazwaz, Analytical treatment of variable coefficients fourth-order parabolic partial differential equations. Applied Mathematics Computation, vol. 123(2001) 219–227. | ||

In article | View Article | ||

[49] | A.M. Wazwa, Exact solutions for variable coefficient fourth order parabolic partial differential equations in higher dimensional spaces, Applied Mathematics Computation vol. 130 (2002) 415-424. | ||

In article | View Article | ||

[50] | J. Biazar, H. Ghazvini, He’s variational iteration method for fourth-order parabolic equation, Computers Mathematics with Applications vol.54 (2007) 1047-1054. | ||

In article | View Article | ||