A numerical algorithm to computationally solve the Hemker problem using Shishkin meshes (Q2114424)

A numerical algorithm is presented to solve a benchmark problem proposed by \textit{P. W. Hemker} [J. Comput. Appl. Math. 76, No. 1--2, 277--285 (1996; Zbl 0870.35020)]. The numerical algorithm is composed of several different Shishkin meshes defined across different coordinate systems aligned to three overlapping subdomains. It is constructed based on parameter explicit pointwise bounds on how the continuous solutions decay away from the circle. Using upwinding in all co-ordinate directions, numerical solutions exhibit no spurious oscillations. Several layer-adapted Shishkin meshes are utilized and these grids are aligned both to the geometry of the domain and to the dominant direction of decay within the boundary/interior layer functions. Numerical experiments indicate that the method is producing accurate approximations over an extensive range of the singular perturbation parameter. The numerical approximations are converging to the continuous solution for each value of the parameter, but not uniformly in the singular perturbation parameter