Numerical method for second order singularly perturbed delay differential equations with fractional order in time via fitted computational method

Abstract

The numerical solution of time fractional parabolic differential equations with singular perturbations and delay is the subject of this article. An arbitrarily small perturbation parameter ɛ(0<ɛ<<1) is multiplied the highest order derivative term. The problem’s solution shows an exponential boundary layer on the right side of the spatial domain for small values of ɛ. The solution and its derivatives are discussed together with their properties and bounds. The Crank–Nicolson approach for time discretization and the fitted operator non-standard finite difference method for space discretization are used to solve the time fractional singularly perturbed time delay problem under consideration. The comparison principle and solution bound are used to examine and analyze the scheme’s stability. The scheme’s uniform convergence is examined and demonstrated. The proposed scheme has a uniform convergence order of Mx−1+(Δt)2−γ. Two numerical illustrations for various values of ɛ are taken into consideration to validate the theoretical analysis of the scheme. The constructed scheme provides a better order of convergence and more accurate results than methods found in the literature.

Keywords

• Singular perturbed problem
• Caputo derivative
• Non-standard finite difference scheme
• Time fractional
• Partial differential equation