Volume 54, Number 1, January-February 2020
|Page(s)||335 - 358|
|Published online||31 January 2020|
Long-time behavior of numerical solutions to nonlinear fractional ODEs
Department of Mathematics and Center for Nonlinear Studies, Northwest University, 710075 Xi’an, Shaanxi, P.R. China
2 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University, 411105 Xiangtan, Hunan, P.R. China
3 Deptartment of Mathematics, The Chinese University of Hong Kong, Shatin, NT, Hong Kong
* Corresponding author: email@example.com
Accepted: 24 August 2019
In this work, we study the long time behavior, including asymptotic contractivity and dissipativity, of the solutions to several numerical methods for fractional ordinary differential equations (F-ODEs). The existing algebraic contractivity and dissipativity rates of the solutions to the scalar F-ODEs are first improved. In order to study the long time behavior of numerical solutions to fractional backward differential formulas (F-BDFs), two crucial analytical techniques are developed, with the first one for the discrete version of the fractional generalization of the traditional Leibniz rule, and the other for the algebraic decay rate of the solution to a linear Volterra difference equation. By means of these auxiliary tools and some natural conditions, the solutions to F-BDFs are shown to be contractive and dissipative, and also preserve the exact contractivity rate of the continuous solutions. Two typical F-BDFs, based on the Grünwald–Letnikov formula and L1 method respectively, are studied. For high order F-BDFs, including convolution quadrature schemes based on classical second order BDF and product integration schemes based on quadratic interpolation approximation, their numerical contractivity and dissipativity are also developed under some slightly stronger conditions. Numerical experiments are presented to validate the long time qualitative characteristics of the solutions to F-BDFs, revealing very different decay rates of the numerical solutions in terms of the the initial values between F-ODEs and integer ODEs and demonstrating the superiority of the structure-preserving numerical methods.
Mathematics Subject Classification: 34A08 / 34D05 / 65L07
Key words: Fractional ODEs / contractivity / dissipativity / fractional BDFs
© EDP Sciences, SMAI 2020
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.