FAST AND PARALLEL RUNGE-KUTTA APPROXIMATION OF FRACTIONAL EVOLUTION EQUATIONS

被引:10
作者
Fischer, Marina [1 ]
机构
[1] Heinrich Heine Univ, Math Inst, D-40225 Dusseldorf, Germany
关键词
convolution quadrature; inverse Laplace transform; Runge-Kutta methods; sub-diffusion equation; parallelizable algorithm; COMPACT FINITE-DIFFERENCE; CONVOLUTION QUADRATURE; OBLIVIOUS CONVOLUTION; LAPLACE TRANSFORM; SCHEME;
D O I
10.1137/18M1175616
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We consider a linear inhomogeneous fractional evolution equation which is obtained from a Cauchy problem by replacing its first-order time derivative with Caputo's fractional derivative. The operator in the fractional evolution equation is assumed to be sectorial. By using the inverse Laplace transform a solution to the fractional evolution equation is obtained which can be written as a convolution. Based on L-stable Runge-Kutta methods a convolution quadrature is derived which allows a stable approximation of the solution. Here, the convolution quadrature weights are represented as contour integrals. On discretizing these integrals, we are able to give an algorithm which computes the solution after N time steps with step size h up to an arbitrary accuracy epsilon. For this purpose the algorithm only requires O(N) Runge-Kutta steps for a large number of scalar linear inhomogeneous ordinary differential equations and the solutions of O(log(N) log(1/epsilon)) linear systems which can be done in parallel. In numerical examples we illustrate the algorithm's performance.
引用
收藏
页码:A927 / A947
页数:21
相关论文
共 39 条
  • [1] Abramowitz M., 1972, HDB MATH FUNCTIONS, DOI DOI 10.1115/1.3625776
  • [2] Arendt W., 2001, VECTOR VALUED LAPLAC, DOI [10.1007/978-3-0348-5075-9, DOI 10.1007/978-3-0348-5075-9]
  • [3] Ashyralyev A., 1994, WELL POSEDNESS PARAB, DOI [10.1007/978-3-0348-8518-8, DOI 10.1007/978-3-0348-8518-8]
  • [4] A KERNEL COMPRESSION SCHEME FOR FRACTIONAL DIFFERENTIAL EQUATIONS
    Baffet, Daniel
    Hesthaven, Jan S.
    [J]. SIAM JOURNAL ON NUMERICAL ANALYSIS, 2017, 55 (02) : 496 - 520
  • [5] Bakaev NY, 2003, MATH COMPUT, V72, P1597, DOI 10.1090/S0025-5718-02-01488-6
  • [6] Runge-Kutta convolution quadrature for operators arising in wave propagation
    Banjai, Lehel
    Lubich, Christian
    Melenk, Jens Markus
    [J]. NUMERISCHE MATHEMATIK, 2011, 119 (01) : 1 - 20
  • [7] BAZHLEKOVA E. G., 2001, THESIS, DOI [10.6100/IR549476, DOI 10.6100/IR549476]
  • [8] On fractional resolvent operator functions
    Chen, Chuang
    Li, Miao
    [J]. SEMIGROUP FORUM, 2010, 80 (01) : 121 - 142
  • [9] Convolution quadrature time discretization of fractional diffusion-wave equations
    Cuesta, E
    Lubich, C
    Palencia, C
    [J]. MATHEMATICS OF COMPUTATION, 2006, 75 (254) : 673 - 696
  • [10] D'ERRICO J., 2012, FRESNELS FRESNELC