Research Papers

Numerical Simulation of Noninteger Order System in Subdiffusive, Diffusive, and Superdiffusive Scenarios

[+] Author and Article Information
Kolade M. Owolabi

Institute for Groundwater Studies,
Faculty of Natural and Agricultural Sciences,
University of the Free State,
Bloemfontein 9300, South Africa
e-mails: mkowolax@yahoo.com;

Abdon Atangana

Institute for Groundwater Studies,
Faculty of Natural and Agricultural Sciences,
University of the Free State,
Bloemfontein 9300, South Africa

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received May 3, 2016; final manuscript received October 28, 2016; published online December 5, 2016. Assoc. Editor: Dumitru Baleanu.

J. Comput. Nonlinear Dynam 12(3), 031010 (Dec 05, 2016) (7 pages) Paper No: CND-16-1219; doi: 10.1115/1.4035195 History: Received May 03, 2016; Revised October 28, 2016

In this work, we investigate both the mathematical and numerical studies of the fractional reaction–diffusion system consisting of spatial interactions of three components’ species. Our main result is based on the analysis of the model for linear stability. Mathematical analysis of the main equation shows that the dynamical system is both locally and globally asymptotically stable. We further propose a theorem which guarantees the existence and permanence of the three species. We formulate a viable numerical methods in space and time. By adopting the Fourier spectral approach to discretize in space, the issue of stiffness associated with the fractional-order spatial derivatives in such system is removed. The resulting system of ordinary differential equations (ODEs) is advanced with the exponential time-differencing method of ADAMS-type. The complexity of the dynamics in the system which we discussed theoretically are numerically presented through some numerical simulations in 1D, 2D, and 3D to address the points and queries that may naturally arise.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Garvie, M. , 2007, “ Finite-Difference Schemes for Reaction-Diffusion Equations Modeling Predator-Pray Interactions in MATLAB,” Bull. Math. Biol., 69(3), pp. 931–956. [CrossRef] [PubMed]
Murray, J. D. , 2003, Mathematical Biology II: Spatial Models and Biomedical Applications, Springer-Verlag, Berlin.
Chang, Y. , Feng, W. , Freeze, M. , and Lu, X. , 2015, “ Permanence and Coexistence in a Diffusive Complex Ratio-Dependent Food Chain,” Int. J. Dyn. Control, 3(3), pp. 262–274. [CrossRef]
Holmes, W. R. , 2014, “ An Efficient, Nonlinear Stability Analysis for Detecting Pattern Formation in Reaction Diffusion Systems,” Bull. Math. Biol., 76(1), pp. 157–183. [CrossRef] [PubMed]
Yu, H. , Zhong, S. , Agarwal, R. P. , and Xiong, L. , 2010, “ Species Permanence and Dynamical Behavior Analysis of an Impulsively Controlled Ecological System With Distributed Time Delay,” Comput. Math. Appl., 59(12), pp. 3824–3835. [CrossRef]
Baleanu, D. , Golmankhaneh, A. , Nigmatullin, R. , and Golmankhaneh, A. , 2010, “ Fractional Newtonian Mechanics,” Open Phys., 8, pp. 120–125. [CrossRef]
Uchaikin, V. V. , 2013, Fractional Derivatives for Physicists and Engineers, Springer, Berlin.
Bhrawy, A. H. , Zaky, M. A. , Baleanu, D. , and Abdelkawy, M. A. , 2015, “ A Novel Spectral Approximation for the Two-Dimensional Fractional Sub-Diffusion Problems,” Rom. J. Phys., 60(3–4), pp. 344–359. http://www.nipne.ro/rjp/2015_60_3-4/0344_0359.pdf
Bhrawy, A. H. , Doha, E. H. , Baleanu, D. , and Ezz-eldein, S. S. , 2015, “ A Spectral Tau Algorithm Based on Jacobi Operational Matrix for Numerical Solution of Time Fractional Diffusion-Wave Equations,” J. Comput. Phys., 293, pp. 142–156. [CrossRef]
Bhrawy, A. H. , 2016, “ A Jacobi Spectral Collocation Method for Solving Multi-Dimensional Nonlinear Fractional Sub-Diffusion Equations,” Numer. Algorithms, 73(1), pp. 91–113.
Bhrawy, A. H. , 2016, “ A New Spectral Algorithm for a Time-Space Fractional Partial Differential Equations With Subdiffusion and Superdiffusion,” Proc. Rom. Acad. A, 17(1), pp. 39–47. http://www.acad.ro/sectii2002/proceedings/doc2016-1/06-Bhrawy.pdf
Owolabi, K. M. , and Atangana, A. , 2016, “ Numerical Solution of Fractional-In-Space Nonlinear Schrodinger Equation With the Riesz Fractional Derivative,” Eur. Phys. J. Plus, 131, p. 335. [CrossRef]
Bainov, D. D. , and Simeonov, P. C. , 1993, Impulsive Differential Equations: Asymptotic Properties of the Solutions, World Scientific, Singapore.
Bhrawy, A. H. , 2014, “ An Efficient Jacobi Pseudospectral Approximation for Nonlinear Complex Generalized Zakharov System,” Appl. Math. Comput., 247, pp. 30–46.
Doha, E. H. , Bhrawy, A. H. , Baleanu, D. , and Hafez, R. M. , 2014, “ A New Jacobi Rational-Gauss Collocation Method for Numerical Solution of Generalized Pantograph Equations,” Appl. Numer. Math., 77, pp. 43–54. [CrossRef]
Owolabi, K. M. , and Patidar, K. C. , 2014, “ Higher-Order Time-Stepping Methods for Time-Dependent Reaction-Diffusion Equations Arising in Biology,” Appl. Math. Comput., 240, pp. 30–50.
Boyd, J. P. , 2001, Chebyshev and Fourier Spectral Methods, Dover, Mineola, NY.
Kassam, A. , and Trefethen, L. N. , 2005, “ Fourth-Order Time-Stepping for Stiff PDEs,” SIAM J. Sci. Comput., 26(4), pp. 1214–1233. [CrossRef]
Owolabi, K. M. , and Patidar, K. C. , 2016, “ Numerical Simulations of Multicomponent Ecological Models With Adaptive Methods,” Theor. Biol. Med. Modell., 13, p. [CrossRef]
Owolabi, K. M. , 2017, “ Robust and Adaptive Techniques for Numerical Simulation of Nonlinear Partial Differential Equations of Fractional Order,” Commun. Nonlinear Sci. Numer. Simul., 44, pp. 304–317. [CrossRef]
Cox, S. M. , and Matthews, P. C. , 2002, “ Exponential Time Differencing for Stiff Systems,” J. Comput. Phys., 176(2), pp. 430–455. [CrossRef]
Du, Q. , and Zhu, W. , 2005, “ Analysis and Applications of the Exponential Time Differencing Schemes and Their Contour Integration Modifications,” BIT Numer. Math., 45(2), pp. 307–328. [CrossRef]
Owolabi, K. M. , 2015, “ Robust IMEX Schemes for Solving Two-Dimensional Reaction-Diffusion Models,” Int. J. Nonlinear Sci. Numer. Simul., 16(6), pp. 271–284.
Pindza, E. , and Owolabi, K. M. , 2016, “ Fourier Spectral Method for Higher Order Space Fractional Reaction-Diffusion Equations,” Commun. Nonlinear Sci. Numer. Simul., 40, pp. 112–128. [CrossRef]
Owolabi, K. M. , and Patidar, K. C. , 2015, “ Existence and Permanence in a Diffusive KiSS Model With Robust Numerical Simulations,” Int. J. Differ. Equations, 2015, p. 485860. [CrossRef]


Grahic Jump Location
Fig. 1

Surface plot showing the distribution of species for reaction–diffusion system (1). The first and second rows areobtained at t = 10 and t = 20, respectively. The parameter values and initial data: η=1.5, D1=0.7,D2=0.5,D3=0.2, u0(x)=0.53x+0.47 sin(−1.5πx),v0(x)=1+sin(2πx),w0(x)=ŵ+10−8(x−1200)(x−2800), ŵ=1/30. Other parameter values as in Eq.(16).

Grahic Jump Location
Fig. 2

Different chaotic behaviors of the species arising from the perturbation of the cropping rate αi,i=1,2, at different instants of time with fractional order η=0.5. Other parameters and initial condition are as used in Fig. 1. Take note of variation in the amplitudes.

Grahic Jump Location
Fig. 3

Contour plots of Eq. (1) in 2D indicating fractional power effects at different instants of η=0.5,1.5 which correspond to subdiffusive and superdiffusive scenarios for the upper and lower-rows, respectively. Other parameters are as fixed in Eq. (18).

Grahic Jump Location
Fig. 4

Numerical simulations of system (1) with kinetics (2) in 3D showing the effects of fractional power η at different instances of η=0.5, η=2.0, and η=1.5 which correspond to the subdiffusive, diffusive, and superdiffusive cases in rows 1–3, respectively. Other parameters are given in Eq. (18).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In