On the l2 Stability of Crank-Nicolson Difference Scheme for Time Fractional Heat Equations

2018 ◽  
Vol 4 (2) ◽  
pp. 111-122
Author(s):  
Serife R. Bayramoglu Erguner ◽  
Ibrahim Karatay
2019 ◽  
Vol 2019 (1) ◽  
Author(s):  
Xiaozhong Yang ◽  
Xu Dang

Abstract The fractional reaction–diffusion equation has profound physical and engineering background, and its rapid solution research is of important scientific significance and engineering application value. In this paper, we propose a parallel computing method of mixed difference scheme for time fractional reaction–diffusion equation and construct a class of improved alternating segment Crank–Nicolson (IASC–N) difference schemes. The class of parallel difference schemes constructed in this paper, based on the classical Crank–Nicolson (C–N) scheme and classical explicit and implicit schemes, combines with alternating segment techniques. We illustrate the unique existence, unconditional stability, and convergence of the parallel difference scheme solution theoretically. Numerical experiments verify the theoretical analysis, which shows that the IASC–N scheme has second order spatial accuracy and $2-\alpha $ 2 − α order temporal accuracy, and the computational efficiency is greatly improved compared with the implicit scheme and C–N scheme. The IASC–N scheme has ideal computation accuracy and obvious parallel computing properties, showing that the IASC–N parallel difference method is effective for solving time fractional reaction–diffusion equation.


Author(s):  
Ibrahim Karatay ◽  
Nurdane Kale ◽  
Serife Bayramoglu

AbstractIn this paper, we consider the numerical solution of a time-fractional heat equation, which is obtained from the standard diffusion equation by replacing the first-order time derivative with the Caputo derivative of order α, where 0 < α < 1. The main purpose of this work is to extend the idea on the Crank-Nicholson method to the time-fractional heat equations. By the method of the Fourier analysis, we prove that the proposed method is stable and the numerical solution converges to the exact one with the order O(τ 2-α + h 2), conditionally. Numerical experiments are carried out to support the theoretical claims.


Sign in / Sign up

Export Citation Format

Share Document