Title The impact of smooth W-grids in the numerical solution of singular perturbation two-point boundary value problems
Authors Gustaf Söderlind, Arjun Singh Yadaw
Alternative Location http://dx.doi.org/10.1016/j..., Restricted Access
Publication Applied Mathematics and Computation
Year 2012
Volume 218
Issue 10
Pages 6045 - 6055
Document type Article
Status Published
Quality controlled Yes
Language eng
Publisher Elsevier Science Inc
Abstract English This paper develops a semi-analytic technique for generating smooth nonuniform grids for the numerical solution of singularly perturbed two-point boundary value problems. It is based on the usual idea of mapping a uniform grid to the desired nonuniform grid. We introduce the W-grid, which depends on the perturbation parameter epsilon << 1. For problems on 0, 1 with a boundary layer at one end point, the local mesh width h(i) = x(i+1) - x(i), with 0 = x(0) < x(1) < ... < x(N) = 1, is condensed at either 0 or 1. Two simple 2nd order finite element and finite difference methods are combined with the new mesh, and computational experiments demonstrate the advantages of the smooth W-grid compared to the well-known piecewise uniform Shishkin mesh. For small epsilon, neither the finite difference method nor the finite element method produces satisfactory results on the Shishkin mesh. By contrast, accuracy is vastly improved on the W-grid, which typically produces the nominal 2nd order behavior in L(2), for large as well as small values of N, and over a wide range of values of epsilon. We conclude that the smoothness of the mesh is of crucial importance to accuracy, efficiency and robustness. Published by Elsevier Inc.
Keywords Singular perturbation, Boundary value problems, Finite difference, method, Galerkin method, Adaptive grid, W-grid, Grid density, Shishkin, mesh,
ISBN/ISSN/Other ISSN: 0096-3003

Questions: webmaster
Last update: 2013-04-11

Centre for Mathematical Sciences, Box 118, SE-22100, Lund. Telefon: +46 46-222 00 00 (vx)