Uniform convergence and order reduction of the fractional implicit Euler method to solve singularly perturbed 2D reaction-diffusion problems

作者:

Highlights:

摘要

In this paper we analyze the uniform convergence of a numerical method designed to approximate efficiently the solution of 2D parabolic singularly perturbed problems of reaction diffusion type. The method combines a modified fractional implicit Euler method to discretize in time, and the classical central finite difference scheme, on a special nonuniform mesh, to discretize in space. The resulting fully discrete scheme is uniformly convergent with respect to the diffusion parameter. The analysis of the convergence is made by using a two step technique, which discretizes first in time and later on in space. We show the order reduction phenomenon associated to the fractional implicit Euler method, which typically appears if the boundary conditions are time dependent and a natural evaluation of them is done. An appropriate choice for the boundary conditions is proposed and analyzed in detail, proving that the order reduction can be removed. Some numerical tests show the practical effects of our method; as well, we compare it with the classical choice for the boundary data in terms of the uniform consistency and the order of uniform convergence of the numerical scheme.

论文关键词:2D parabolic singularly perturbed problems,Fractional implicit Euler method,Nonuniform special meshes,Uniform convergence,Order reduction

论文评审过程:Received 7 August 2015, Accepted 12 April 2016, Available online 31 May 2016, Version of Record 31 May 2016.

论文官网地址:https://doi.org/10.1016/j.amc.2016.04.017