Stability and Convergence of an Implicit Numerical Method for the Non-Linear Fractional Reaction-Subdiffusion Process

P. Zhuang,F. Liu,V. Anh,I. Turner
DOI: https://doi.org/10.1093/imamat/hxp015
2009-01-01
IMA Journal of Applied Mathematics
Abstract:In this paper, we consider the following non-linear fractional reaction-subdiffusion process (NFR-SubDP):partial derivative u/partial derivative t = D-0(t)1-gamma [K gamma partial derivative(2)u/partial derivative x(2) + f(u, x, t)] g(u, x, t),where f (u, x, t) is a linear function of u, the function g (u, x, t) satisfies the Lipschitz condition and D-0(t)1-gamma is the Riemann-Liouville time fractional partial derivative of order 1 - gamma. We propose a new computationally efficient numerical technique to simulate the process. Firstly, the NFR-SubDP is decoupled, which is equivalent to solving a non-linear fractional reaction-subdiffusion equation (NFR-SubDE). Secondly, we propose an implicit numerical method to approximate the NFR-SubDE. Thirdly, the stability and convergence of the method are discussed using a new energy method. Finally, some numerical examples are presented to show the application of the present technique. This method and supporting theoretical results can also be applied to fractional integrodifferential equations.
What problem does this paper attempt to address?