Linearly implicit predictor–corrector methods for space-fractional reaction–diffusion equations with non-smooth initial data
Autor: | S.S. Alzahrani, T. A. Biala, A. Q. M. Khaliq, Khaled M. Furati |
---|---|
Rok vydání: | 2018 |
Předmět: |
Predictor–corrector method
Discretization 010103 numerical & computational mathematics Function (mathematics) Derivative 01 natural sciences 010101 applied mathematics Computational Mathematics Matrix (mathematics) Computational Theory and Mathematics Modeling and Simulation Reaction–diffusion system Padé approximant Applied mathematics Matrix exponential 0101 mathematics Mathematics |
Zdroj: | Computers & Mathematics with Applications. 75:2629-2657 |
ISSN: | 0898-1221 |
DOI: | 10.1016/j.camwa.2017.12.033 |
Popis: | In this paper, linearly implicit predictor–corrector methods are proposed for solving space-fractional reaction–diffusion equations with non-smooth initial data. The methods are based on Matrix Transfer Technique for spatial discretization and are shown to be unconditionally stable. It is observed that the linearly implicit predictor–corrector method derived by using (1,1)-Pade approximation to matrix exponential function incurs oscillatory behavior for some time steps. These oscillations are due to high frequency components present in the solution and are diminished as the order of the space-fractional derivative decreases (slow diffusion). We present a priori reliability constraint to avoid these unwanted oscillations and generalize the constraints for all ( m , m ) -Pade approximants, m ∈ Z + , to the matrix exponential functions. These time stepping constraints are seen to be dependent on the order of the space-fractional derivative. The linearly implicit predictor–corrector method based on the (0,2)-Pade approximations to the matrix exponential function is shown to be oscillation-free for any time step. Error estimates are obtained for the methods and are theoretically shown to be second-order convergent. Computational complexity of the algorithms is discussed for solving multidimensional space-fractional reaction–diffusion systems. Several numerical experiments are performed to support our theoretical observations and to show the effectiveness, reliability, and efficiency of the methods. |
Databáze: | OpenAIRE |
Externí odkaz: |