Long-time behavior of numerical solutions to nonlinear fractional ODEs

In this work, we study the long time behaviors, 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 mens 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 Grunwald-Letnikov formula and L1 method respectively, are studied. For high order F-BDFs, including some second order F-BDFs and $3-\alpha$ order method, 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.

[1]  Da Xu,et al.  Decay Properties for the Numerical Solutions of a Partial Differential Equation with Memory , 2014, Journal of Scientific Computing.

[2]  C. Lubich,et al.  On the Stability of Linear Multistep Methods for Volterra Convolution Equations , 1983 .

[3]  G. Dahlquist.  Error analysis for a class of methods for stiff non-linear initial value problems , 1976 .

[4]  Germund Dahlquist,et al.  G-stability is equivalent toA-stability , 1978 .

[5]  Bangti Jin,et al.  Numerical Analysis of Nonlinear Subdiffusion Equations , 2017, SIAM J. Numer. Anal..

[6]  A. T. Hill,et al.  Global Dissipativity for A-Stable Methods , 1997 .

[7]  Wanrong Cao,et al.  Time-Splitting Schemes for Fractional Differential Equations I: Smooth Solutions , 2015, SIAM J. Sci. Comput..

[8]  N. Ford,et al.  A Predictor-Corrector Approach for the Numerical Solution of Fractional Differential Equations , 2013 .

[9]  Xu Da,et al.  Uniform l1 Behavior for Time Discretization of a Volterra Equation with Completely Monotonic Kernel II: Convergence , 2007, SIAM J. Numer. Anal..

[10]  Roberto Garrappa,et al.  On Multistep Methods for Differential Equations of Fractional Order , 2006 .

[11]  C. Lubich,et al.  Fractional linear multistep methods for Abel-Volterra integral equations of the second kind , 1985 .

[12]  Bangti Jin,et al.  Correction of High-Order BDF Convolution Quadrature for Fractional Evolution Equations , 2017, SIAM J. Sci. Comput..

[13]  Eduardo Cuesta,et al.  ASYMPTOTIC BEHAVIOUR OF THE SOLUTIONS OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS AND SOME TIME DISCRETIZATIONS , 2007 .

[14]  Changpin Li,et al.  High-order approximation to Caputo derivatives and Caputo-type advection-diffusion equations (III) , 2016, J. Comput. Appl. Math..

[15]  John C. Butcher,et al.  A stability property of implicit Runge-Kutta methods , 1975 .

[16]  Bangti Jin,et al.  On nonnegativity preservation in finite element methods for subdiffusion equations , 2015, Math. Comput..

[17]  Changpin Li,et al.  High-order approximation to Caputo derivatives and Caputo-type advection-diffusion equations (II) , 2015 .

[18]  Ernst Hairer,et al.  Examples of Stiff Equations , 1996 .

[19]  J. Cermák,et al.  On explicit stability conditions for a linear fractional difference system , 2015 .

[20]  H. Srivastava,et al.  Theory and Applications of Fractional Differential Equations , 2006 .

[21]  Roberto Garrappa,et al.  Trapezoidal methods for fractional differential equations: Theoretical and computational aspects , 2015, Math. Comput. Simul..

[22]  Bangti Jin,et al.  An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data , 2015, 1501.00253.

[23]  I. Petráš.  Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation , 2011 .

[24]  I. Podlubny.  Fractional differential equations , 1998 .

[25]  Andrew M. Stuart,et al.  Runge-Kutta methods for dissipative and gradient dynamical systems , 1994 .

[26]  Chuanju Xu,et al.  Error Analysis of a High Order Method for Time-Fractional Diffusion Equations , 2016, SIAM J. Sci. Comput..

[27]  Changpin Li,et al.  A survey on the stability of fractional differential equations , 2011 .

[28]  Minghua Chen,et al.  Fourth Order Accurate Scheme for the Space Fractional Diffusion Equations , 2014, SIAM J. Numer. Anal..

[29]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[30]  D. Matignon.  Stability results for fractional differential equations with applications to control processing , 1996 .

[31]  J. Hale.  Asymptotic Behavior of Dissipative Systems , 1988 .

[32]  Igor Podlubny,et al.  Mittag-Leffler stability of fractional order nonlinear dynamic systems , 2009, Autom..

[33]  Eduardo Cuesta,et al.  Convolution quadrature time discretization of fractional diffusion-wave equations , 2006, Math. Comput..

[34]  Dongling Wang,et al.  Dissipativity and Contractivity Analysis for Fractional Functional Differential Equations and their Numerical Approximations , 2019, SIAM J. Numer. Anal..

[35]  István Győri,et al.  On exact convergence rates for solutions of linear systems of Volterra difference equations , 2006 .

[36]  N. Ford,et al.  Analysis of Fractional Differential Equations , 2002 .

[37]  R. Temam.  Infinite Dimensional Dynamical Systems in Mechanics and Physics Springer Verlag , 1993 .

[38]  Jiwei Zhang,et al.  Sharp Error Estimate of the Nonuniform L1 Formula for Linear Reaction-Subdiffusion Equations , 2018, SIAM J. Numer. Anal..

[39]  Yubin Yan,et al.  A higher order numerical method for time fractional partial differential equations with nonsmooth data , 2018, J. Comput. Phys..

[40]  Yan Xu,et al.  A Local Discontinuous Galerkin Method for the Camassa-Holm Equation , 2008, SIAM J. Numer. Anal..

[41]  Chuanju Xu,et al.  Finite difference/spectral approximations for the time-fractional diffusion equation , 2007, J. Comput. Phys..

[42]  A. Alikhanov.  A priori estimates for solutions of boundary value problems for fractional-order equations , 2010, 1105.4592.

[43]  Roberto Garrappa,et al.  On linear stability of predictor–corrector algorithms for fractional differential equations , 2010, Int. J. Comput. Math..

[44]  Zhimin Zhang,et al.  Fast evaluation of the Caputo fractional derivative and its applications to fractional diffusion equations , 2015, 1511.03453.

[45]  J. Butcher.  Thirty years of G-stability , 2006 .

[46]  Yubin Yan,et al.  An Analysis of the Modified L1 Scheme for Time-Fractional Partial Differential Equations with Nonsmooth Data , 2018, SIAM J. Numer. Anal..

[47]  Zhi-Zhong Sun,et al.  A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications , 2014, J. Comput. Phys..

[48]  Olavi Nevanlinna,et al.  On the numerical solutions of some Volterra equations on infinite intervals , 1976 .

[49]  A. Xiao,et al.  Dissipativity and contractivity for fractional-order systems , 2015 .

[50]  Wanrong Cao,et al.  Implicit-Explicit Difference Schemes for Nonlinear Fractional Differential Equations with Nonsmooth Solutions , 2016, SIAM J. Sci. Comput..

[51]  P. Eggermont.  Uniform error estimates of Galerkin methods for monotone Abel-Volterra integral equations on the half-line , 1989 .

[52]  Natalia Kopteva,et al.  Error analysis of the L1 method on graded and uniform meshes for a fractional-derivative problem in two and three dimensions , 2017, Math. Comput..

[53]  Zhi‐zhong Sun,et al.  A fully discrete difference scheme for a diffusion-wave system , 2006 .