and -Spectral Element Methods for Two-dimensional Elliptic Boundary Layer Problems
En palabras de los autores
We propose and -spectral element methods for elliptic boundary layer problems on two dimensional rectangular domains. To resolve boundary layers, we use -version with a fixed mesh and an -version with a boundary layer mesh consisting of thin needle-like elements near the boundary layer and coarse elements away from the layer. Stability estimates are derived using non-conforming spectral element functions. The numerical scheme for both the methods is based on minimising a least-squares functional in appropriate Sobolev norms. We construct robust preconditioners to manage the condition number of the normal equations derived from the least-squares formulation. The method is able to approximate boundary layers at a rate for the -version and at the rate , uniformly in , for the -version, where is the boundary layer parameter, is a constant and denotes the degree of the approximating polynomial. Numerical results are provided for model elliptic boundary layer problems for a range of boundary layer parameters. Simulation results demonstrate the efficiency and robustness of the method in capturing boundary layers on rectangular domains.
Apareció: viernes, 25 de septiembre. arXiv. Preprint, todavía sin revisión por pares.
Comentario de los autores: 34 pages, 24 figures (including graphs)