2 research outputs found
An implicit Keller Box numerical scheme for the solution of fractional subdiffusion equations
In this work, we present a new implicit numerical scheme for fractional subdiffusion equations. In this approach we use the Keller Box method [1] to spatially discretise the fractional subdiffusion equation and we use a modified L1 scheme (ML1), similar to the L1 scheme originally developed by Oldham and Spanier [2], to approximate the fractional derivative. The stability of the proposed method was investigated by using Von-Neumann stability analysis. We have proved the method is unconditionally stable when and , and demonstrated the method is also stable numerically in the case and . The accuracy and convergence of the scheme was also investigated and found to be of order in time and in space. To confirm the accuracy and stability of the proposed method we provide three examples with one including a linear reaction term