Schemes Convergent ε-Uniformly for Parabolic Singularly Perturbed Problems with a Degenerating Convective Term and a Discontinuous Source

Abstract

We consider the numerical approximation of a 1D singularly perturbed convection-diffusion problem with a multiply degenerating convective term, for which the order of degeneracy is 2p + 1, p is an integer with p ≥ 1, and such that the convective flux is directed into the domain. The solution exhibits an interior layer at the degeneration point if the source term is also a discontinuous function at this point. We give appropriate bounds for the derivatives of the exact solution of the continuous problem, showing its asymptotic behavior with respect to the perturbation parameter ε, which is the diffusion coefficient. We construct a monotone finite difference scheme combining the implicit Euler method, on a uniform mesh, to discretize in time, and the upwind finite difference scheme, constructed on a piecewise uniform Shishkin mesh condensing in a neighborhood of the interior layer region, to discretize in space. We prove that the method is convergent uniformly with respect to the parameter ε, i.e., ε-uniformly convergent, having first order convergence in time and almost first order in space. Some numerical results corroborating the theoretical results are showed.