Stability of fully discrete schemes with interpolation-type fractional formulas for distributed-order subdiffusion equations (Q2407913)

The authors study the stability of the two numerical schemes developed for the distributed order subdiffusion equation in Capthe Caputo derivative with the midpoint quadrature of the weighted integral. This is referred to as CBDF1 method which preserves the discrete minimum-maximum principle and the non-negative preserving property. In the second scheme the Caputo's BDF2 (L1-2) formula is used to approximate the fractional derivative. The resulting method, denoted by CBDF2, is shown to be stable and convergent. In order to establish the numerical schemes, two numerical examples are considered. It is observed that the solution of distributed order subdiffusion admits a weak initial singularity even for smooth initial data. Several lemmas and theorem are stated and proved.