Research Article | Open Access

Jagdev Singh, Devendra Kumar, Adem Kılıçman, "Numerical Solutions of Nonlinear Fractional Partial Differential Equations Arising in Spatial Diffusion of Biological Populations", *Abstract and Applied Analysis*, vol. 2014, Article ID 535793, 12 pages, 2014. https://doi.org/10.1155/2014/535793

# Numerical Solutions of Nonlinear Fractional Partial Differential Equations Arising in Spatial Diffusion of Biological Populations

**Academic Editor:**Hassan Eltayeb

#### Abstract

The main aim of this work is to present a user friendly numerical algorithm based on homotopy perturbation Sumudu transform method for nonlinear fractional partial differential arising in spatial diffusion of biological populations in animals. The movements are made generally either by mature animals driven out by invaders or by young animals just reaching maturity moving out of their parental territory to establish breeding territory of their own. The homotopy perturbation Sumudu transform method is a combined form of the Sumudu transform method and homotopy perturbation method. The obtained results are compared with Sumudu decomposition method. The numerical solutions obtained by the proposed method indicate that the approach is easy to implement and accurate. These results reveal that the proposed method is computationally very attractive.

#### 1. Introduction

Nonlinear phenomena have important effects on applied mathematics, biology, and physics and are related to engineering; many such physical phenomena are modeled in terms of nonlinear partial differential equations. For example, the degenerate parabolic nonlinear partial differential equations arising in the spatial diffusion of biological populations are given as with initial condition , where denotes the population density and represents the population supply due to birth and death. It is worth mentioning that (1) with , where , , arises in the study of flow through porous media [1, 2]. Moreover, the problem (1) leads to Malthusian law [3] for and Verhulst law [3] for , where are positive constants. We consider the model as an example in the population of animals. The movements are made generally either by mature animals driven out by invaders or by young animals just reaching maturity moving out of their parental territory to establish breeding territory of their own. In both cases, it is much more plausible to suppose that they will be directed towards nearby vacant territory. Therefore, movement will take place almost exclusively “down” the population density gradient and will be much more rapid at high population densities than at low ones. In an attempt to model this situation, they considered a walk through a rectangular grid, in which at each step an animal may either stay at its present location or may move in the direction of lowest population density. We consider the more general form of time fractional biological equation by taking as with initial condition , where denotes the populations density and represents the population supply due to births and deaths. The most important advantage of using fractional differential equations in mathematical modelling is their nonlocal property. It is well known that the integer order differential operator is a local operator, but the fractional order differential operator is nonlocal. This means that the next state of a system depends not only upon its current state but also upon all of its historical states. This is more realistic and it is one reason why fractional calculus has become more and more popular [4–10]. There are several methods to solve the fractional biological equations, such as Adomian’s decomposition method (ADM) [11] and homotopy analysis method (HAM) [12]. The optimal homotopy analysis method (OHAM) has been applied for solving nonlinear problems [13–15]. The homotopy perturbation method (HPM) has first proposed by He [16–18] for solving various linear and nonlinear initial and boundary value problems. The HPM was also investigated by many researchers to handle nonlinear equations arising in science and engineering [19–22]. The Laplace transform and Sumudu transform are totally incapable of handling nonlinear equations because of the difficulties that are caused by the nonlinear terms. Various ways have been proposed recently to deal with such nonlinearities such as the Laplace decomposition algorithm [23, 24], the homotopy perturbation transform method (HPTM) [25], Sumudu decomposition method (SDM) [26], and the homotopy perturbation Sumudu transform method (HPSTM) [27] to produce highly effective techniques for solving many nonlinear problems.

In this paper, we implement the homotopy perturbation Sumudu transform method (HPSTM) to find the analytical and numerical solutions of the fractional biological population models. The HPSTM is an elegant combination of the Sumudu transform method, the HPM and Adomian’s polynomials. Moreover, we apply the Sumudu decomposition method (SDM) to the same problems to compare the results. The SDM is a combined form of Sumudu transform method and Adomian’s decomposition method (ADM) which was first introduced by Adomian [28, 29]. Adomian’s polynomials are studied by Jafari et al. [30, 31] for handling nonlinear equations and showed that the Adomian’s polynomials and He’s polynomials are same. The numerical results are depicted through the graphical illustrations. The proposed techniques provide the solution in a rapid convergent series which may lead to the solution in a closed form. The advantage of these methods is their capability of combining two powerful methods for obtaining exact and approximate solutions for nonlinear equations.

#### 2. Basic Definitions of Fractional Calculus

In this section, we mention the following basic definitions of fractional calculus and Sumudu transform.

*Definition 1. *The Riemann-Liouville fractional integral operator of order , of a function , is defined as [4]
For the Riemann-Liouville fractional integral we have

*Definition 2. *The fractional derivative of in the Caputo sense is defined as [5]
for .

*Definition 3. *The Sumudu transform [32] is defined over the set of functions:
by the following formula:
Some of the important properties of the Sumudu transform were established by Asiru [33]. Further, fundamental properties of this transform were established by Belgacem et al. [34] and Belgacem and Karaballi [35]. In fact it was shown that there is strong relationship between Sumudu and other integral transform; see Kılıçman et al. [36]. In particular the relationship between Sumudu transform and Laplace transforms was presented in Kılıcman and Gadain [37]. The Sumudu transform has scale and unit preserving properties, so it can be used to solve problems without resorting to a new frequency domain.

*Definition 4. *The Sumudu transform of the Caputo fractional derivative is defined as follows [38]:

#### 3. Homotopy Perturbation Sumudu Transform Method (HPSTM)

To illustrate the basic idea of this method, we consider a general fractional nonlinear nonhomogeneous partial differential equation with the initial condition of the form: where is the Caputo fractional derivative of the function , is the linear differential operator, represents the general nonlinear differential operator, and , is the source term.

Applying the Sumudu transform (denoted in this paper by ) on both sides of (9), we get Using the differentiation property of the Sumudu transform, we have Operating with the Sumudu inverse on both sides of (12) gives where represents the term arising from the source term and the prescribed initial conditions. Now we apply the HPM and the nonlinear term can be decomposed as for some Adomian’s polynomials that are given by Substituting (14) and (15) in (13), we get which is the coupling of the Sumudu transform method, the HPM, and Adomian’s polynomials. Comparing the coefficients of like powers of , the following approximations are obtained: Proceeding in the same manner, the rest of the components can be completely found and the series solution is thus entirely determined. Finally, we approximate the analytical solution by truncated series as

#### 4. Sumudu Decomposition Method (SDM)

We consider a general fractional nonlinear nonhomogeneous partial differential equation with the initial condition of the form: where is the Caputo fractional derivative of the function is the linear differential operator, represents the general nonlinear differential operator, and is the source term.

Taking the Sumudu transform on both sides of (20), we get Using the differentiation property of the Sumudu transform and the above initial conditions, we have Now, applying the inverse Sumudu transform on both sides of (23), we get where represents the term arising from the source term and the prescribed initial conditions.

The second step in SDM is that we represent solution as an infinite series given below: and the nonlinear term can be decomposed as where are Adomian polynomials of and it can be calculated by formula given below: Using (25) and (26) in (24), we get On comparing both sides of (28), we get In general the recursive relation is given by Now first of all applying the Sumudu transform of the right hand side of (30) then applying the inverse Sumudu transform, we get the values of , respectively.

#### 5. Numerical Examples

In this section, we demonstrate the performance and efficiency of HPSTM by giving three examples of time-fractional biological models and compare the results with the results obtained with help of SDM.

*Example 1. *We consider the following time-fractional biological population model:
subject to the initial condition
Operating with the Sumudu transform on both sides of (31), we have
The inverse of the Sumudu transform gives that
According to the HPM, we construct the following homotopy:
where are Adomian’s polynomials that represent the nonlinear terms. The first few components of Adomian’s polynomials are given by
Comparing the coefficient of like power of , we have
and so on.

Making use of (37), the series solution is given as
Setting in (38), we reproduce the solution of the problem as follows:
which converge very rapidly to the exact solution .

According to SDM and to (25) and (34), the iteration formula for (31) is given by
where are Adomian’s polynomials that represent the nonlinear terms. The first few components of Adomian’s polynomials are given by
The recursive relation is given below:
The other components of the solution can be easily found by using above recursive relation:
and so on.

Therefore, the decomposition series solutions is given by
which is the same solution as obtained by using HPSTM and setting *α* = 1, it converge to the exact solution .

The numerical results for the time-fractional biological population model (31) obtained by using the HPSTM, SDM, and the exact solution for various values of , and with and are shown by Figures 1(a)–1(d) and those for various values of at , , and are depicted in Figure 2 and those for different values of and at , and are shown in Figure 3. It is observed from Figure 1 that increases with the increase in both and for , and 1 with and . It is seen from Figure 3 that as the value of *α* decreases, the numerical value of increases. It can also be seen from Figure 1 that the solution obtained by using the present method is nearly identical with the exact solution. From Figure 2, we can easily observe that the results obtained with help of HPSTM and SDM are in a full agreement with the exact solution. It is to be noted that only the seventh order term of the HPSTM and SDM was used in evaluating the approximate solutions for Figures 1 and 2. It is evident that the efficiency of the present method can be dramatically enhanced by computing further terms of when the HPSTM and SDM are used.

**(a)**

**(b)**

**(c)**

**(d)**

*Example 2. *Next, we consider the time-fractional biological population model in the following form:
subject to initial condition
By applying the HPSTM, we have
Equating the like powers of , (47) yields
and so on.

Taking use of (48), the series solution is presented as
Putting in (49), we reproduce the solution of the problem as follows:
which converge very rapidly to the exact solution .

According to SDM, we have
The recursive relation is given below:
The other components of the solution can be easily found by using above recursive relation:
and so on.

Thus, the decomposition series solutions is given by
which is the same solution as obtained by employing HPSTM and setting ; it converses to the exact solution .

The numerical results for the time-fractional biological population model (45) obtained with the help of HPSTM, SDM, and the exact solution for various values of , and with are described through Figures 4(a)–4(d) and those for various values of at , and are given in Figure 5 and those for different values of and *α* at and are depicted in Figure 6. It is observed from Figure 4 that increases with the increase in both and for , and 1 with . It is to be observed from Figure 6 that as the value of decrease, the value of increases. It can also be seen from Figure 4 that the solution obtained by using the present method is nearly identical with the exact solution. From Figure 5, we can easily observe that the results obtained with help of HPSTM and SDM are in a full agreement with the exact solution. It is to be noted that only the seventh order term of the HPSTM and SDM was used in evaluating the approximate solutions for Figures 4 and 5.

**(a)**

**(b)**

**(c)**

**(d)**

*Example 3. *Finally, we consider the following time-fractional biological population model:
subject to the initial condition
By applying the HPSTM, we have
Equating the terms with identical powers of , we have
and so on.

Using (58), the series solution is
Taking in (59), we reproduce the solution of the problem as follows:
which converge very rapidly to the exact solution .

According to SDM, we have
The recursive relation is given below:
The other components of the solution can be easily found by using above recursive relation:
and so on.

So, the decomposition series solutions is given by
which is the same solution as obtained by the application of HPSTM and setting *α* = 1; it converges to the exact solution .

The numerical results for the time-fractional biological population model (55) obtained with the help of HPSTM, SDM, and the exact solution for various values of , and with are described through Figures 7(a)–7(d) and those for various values of at , and are depicted in Figure 8 and those for different values of and at and are presented in Figure 9. It is observed from Figure 7 that increases with the increase in and decreases with increase in for , and 1 with . It is to be observed from Figure 9 that as the value of increase, the value of increases, but afterward its nature is opposite. It can also be seen from Figure 7 that the solution obtained by using the present technique is nearly identical with the exact solution. From Figure 8, we can easily observe that the results obtained with help of HPSTM and SDM are in an excellent agreement with the exact solution. It is to be noted that only the seventh order term of the HPSTM and SDM was used in evaluating the approximate solutions for Figures 7 and 8.

**(a)**

**(b)**

**(c)**

**(d)**

#### 6. Conclusions

In this paper, the HPSTM and SDM are successfully employed to solve the time-fractional biological population models. The proposed methods provide the solutions in terms of convergent series with easily computable components in a direct way without using linearization, perturbation, or restrictive assumptions. The numerical results obtained with the proposed techniques are in an excellent agreement with the exact solution. It is worth mentioning that both methods are capable of reducing the volume of the computational work as compared to the classical methods while still maintaining the high accuracy of the numerical result; the size reduction amounts to an improvement of the performance of the approach for solving the time-fractional biological population models. Hence, we conclude that the HPSTM and SDM are very powerful and efficient in finding analytical as well as numerical solutions for wide classes of nonlinear fractional partial differential equations arising in science engineering and finance.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors are extending their heartfelt thanks to the reviewers for their valuable suggestions for the improvement of the paper. The authors also gratefully acknowledge that this research was partially supported by the University Putra Malaysia under the GP-IBT Grant Scheme having Project no. GP-IBT/2013/9420100.

#### References

- J. Bear,
*Dynamics of Fluids in Porous Media*, American Elsevier, New York, NY, USA, 1972. - A. Okubo,
*Diffusion and Ecological Problem: Mathematical Models*, vol. 10 of*Biomathematics*, Springer, Berlin, Germany, 1980. View at: MathSciNet - M. E. Gurtin and R. C. MacCamy, “On the diffusion of biological populations,”
*Mathematical Biosciences*, vol. 33, no. 1-2, pp. 35–49, 1977. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - I. Podlubny,
*Fractional Differential Equations*, Academic Press, New York, NY, USA, 1999. View at: MathSciNet - M. Caputo,
*Elasticita e Dissipazione*, Zani-Chelli, Bologna, Italy, 1969. - R. Gorenflo, F. Mainardi, E. Scalas, and M. Roberto, “Fractional calculus and continuous time finance III: the diffusion limit,” in
*Mathematical Finance*, Trends in Mathematics, pp. 171–180, Birkhäuser, Basel, Switzerland, 2001. View at: Publisher Site | Google Scholar - R. Hilfer, Ed.,
*Applications of Fractional Calculus in Physics*, World Scientific Publishing Company, Singapore, 2000. - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, John Wiley & Sons, New York, NY, USA, 1993. View at: MathSciNet - K. B. Oldham and J. Spanier,
*The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order*, Academic Press, New York, NY, USA, 1974. - A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and Applications of Fractional Differential Equations*, Elsevier Science, Amsterdam, The Netherlands, 2006. View at: MathSciNet - A. M. A. El-Sayed, S. Z. Rida, and A. A. M. Arafa, “Exact solutions of fractional-order biological population model,”
*Communications in Theoretical Physics*, vol. 52, no. 6, pp. 992–996, 2009. View at: Publisher Site | Google Scholar | MathSciNet - A. A. M. Arafa, S. Z. Rida, and H. Mohamed, “Homotopy analysis method for solving biological population model,”
*Communications in Theoretical Physics*, vol. 56, no. 5, pp. 797–800, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. Basiri Parsa, M. M. Rashidi, O. Anwar Bég, and S. M. Sadri, “Semi-computational simulation of magneto-hemodynamic flow in a semi-porous channel using optimal homotopy and differential transform methods,”
*Computers in Biology and Medicine*, vol. 43, no. 9, pp. 1142–1153, 2013. View at: Publisher Site | Google Scholar - H. N. Hassan and M. M. Rashidi, “Analytical solution for three-dimensional steady flow of condensation film on inclined rotating disk by optimal homotopy analysis method,”
*Walailak Journal of Science and Technology*, vol. 10, no. 5, pp. 479–498, 2013. View at: Google Scholar - M. M. Rashidi, E. Momoniat, and M. Ferdows, “Lie group solution for free convective flow of a nanofluid past a chemically reacting horizontal plate in a porous media,”
*Mathematical Problems in Engineering*, vol. 2014, Article ID 239082, 21 pages, 2014. View at: Publisher Site | Google Scholar | MathSciNet - J. H. He, “Homotopy perturbation technique,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 178, no. 3-4, pp. 257–262, 1999. View at: Publisher Site | Google Scholar | MathSciNet - J. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,”
*International Journal of Non-Linear Mechanics*, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site | Google Scholar | MathSciNet - J.-H. He, “Some asymptotic methods for strongly nonlinear equations,”
*International Journal of Modern Physics B*, vol. 20, no. 10, pp. 1141–1199, 2006. View at: Publisher Site | Google Scholar | MathSciNet - D. D. Ganji, “The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer,”
*Physics Letters A*, vol. 355, no. 4-5, pp. 337–341, 2006. View at: Publisher Site | Google Scholar | MathSciNet - H. Jafari, M. Zabihi, and M. Saidy, “Application of homotopy-perturbation method for solving gas dynamics equation,”
*Applied Mathematical Sciences*, vol. 2, no. 48, pp. 2393–2396, 2008. View at: Google Scholar | MathSciNet - M. M. Rashidi, D. D. Ganji, and S. Dinarvand, “Explicit analytical solutions of the generalized BURger and BURger-Fisher equations by homotopy perturbation method,”
*Numerical Methods for Partial Differential Equations*, vol. 25, no. 2, pp. 409–417, 2009. View at: Publisher Site | Google Scholar | MathSciNet - M. M. Rashidi and D. D. Ganji, “Homotopy perturbation combined with Padé approximation for solving two dimensional viscous flow in the extrusion process,”
*International Journal of Nonlinear Science*, vol. 7, no. 4, pp. 387–394, 2009. View at: Google Scholar | MathSciNet - S. A. Khuri, “A Laplace decomposition algorithm applied to a class of nonlinear differential equations,”
*Journal of Applied Mathematics*, vol. 1, no. 4, pp. 141–155, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Khan and M. Hussain, “Application of Laplace decomposition method on semi-infinite domain,”
*Numerical Algorithms*, vol. 56, no. 2, pp. 211–218, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using He's polynomials,”
*Computers & Mathematics with Applications*, vol. 61, no. 8, pp. 1963–1967, 2011. View at: Publisher Site | Google Scholar | MathSciNet - D. Kumar, J. Singh, and S. Rathore, “Sumudu decomposition method for nonlinear equations,”
*International Mathematical Forum*, vol. 7, no. 11, pp. 515–521, 2012. View at: Google Scholar | MathSciNet - J. Singh, D. Kumar, and A. Kılıçman, “Homotopy perturbation method for fractional gas dynamics equation using Sumudu transform,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 934060, 8 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - G. Adomian,
*Solving Frontier Problems of Physics: The Decomposition Method*, Kluwer Academic Publishers, Boston, Mass, USA, 1994. View at: Publisher Site | MathSciNet - G. Adomian and S. E. Serrano, “Stochastic contaminant transport equation in porous media,”
*Applied Mathematics Letters*, vol. 11, no. 1, pp. 53–55, 1998. View at: Google Scholar - H. Jafari, S. Ghasempour, and C. M. Khalique, “Comments on “He's homotopy perturbation method for calculating adomian polynomials”,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 14, no. 6, pp. 339–343, 2013. View at: Google Scholar - H. Jafari, S. Ghasempoor, and C. M. Khalique, “A comparison between Adomian's polynomials and He's polynomials for nonlinear functional equations,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 943232, 4 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - G. K. Watugala, “Sumudu transform—a new integral transform to solve differential equations and control engineering problems,”
*International Journal of Mathematical Education in Science and Technology*, vol. 24, no. 1, pp. 35–43, 1993. View at: Publisher Site | Google Scholar | MathSciNet - M. A. Asiru, “Sumudu transform and the solution of integral equations of convolution type,”
*International Journal of Mathematical Education in Science and Technology*, vol. 32, no. 6, pp. 906–910, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. B. M. Belgacem, A. A. Karaballi, and S. L. Kalla, “Analytical investigations of the Sumudu transform and applications to integral production equations,”
*Mathematical Problems in Engineering*, vol. 2003, no. 3, pp. 103–118, 2003. View at: Publisher Site | Google Scholar | MathSciNet - F. B. M. Belgacem and A. Karaballi, “Sumudu transform fundamental properties investigations and applications,”
*Journal of Applied Mathematics and Stochastic Analysis*, vol. 2006, Article ID 91083, 23 pages, 2006. View at: Publisher Site | Google Scholar | MathSciNet - A. Kılıçman, H. Eltayeb, and M. A. K. Ariffin, “A note on the comparison between Laplace and Sumudu transforms,”
*Bulletin of the Iranian Mathematical Society*, vol. 37, no. 1, pp. 131–141, 2011. View at: Google Scholar - A. Kılıcman and H. E. Gadain, “On the applications of Laplace and Sumudu transforms,”
*Journal of the Franklin Institute*, vol. 347, no. 5, pp. 848–862, 2010. View at: Publisher Site | Google Scholar | MathSciNet - V. B. L. Chaurasia and J. Singh, “Application of Sumudu transform in Schrödinger equation occurring in quantum mechanics,”
*Applied Mathematical Sciences*, vol. 4, no. 57–60, pp. 2843–2850, 2010. View at: Google Scholar | MathSciNet

#### Copyright

Copyright © 2014 Jagdev Singh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.