Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
8902660 | Applied Numerical Mathematics | 2018 | 16 Pages |
Abstract
Fully diagonalized spectral methods using Sobolev orthogonal/biorthogonal basis functions are proposed for solving second order elliptic boundary value problems. We first construct the Fourier-like Sobolev polynomials which are mutually orthogonal (resp. bi-orthogonal) with respect to the bilinear form of the symmetric (resp. unsymmetric) elliptic Neumann boundary value problems. The exact and approximation solutions are then expanded in an infinite and truncated series in the Sobolev orthogonal polynomials, respectively. An identity is also established for the a posterior error estimate with a simple error indicator. Further, the Fourier-like Sobolev orthogonal polynomials and the corresponding Legendre spectral method are proposed in parallel for Dirichlet boundary value problems. Numerical experiments illustrate that our Legendre methods proposed are not only efficient for solving elliptic problems but also equally applicable to indefinite Helmholtz equations and singular perturbation problems.
Related Topics
Physical Sciences and Engineering
Mathematics
Computational Mathematics
Authors
Qing Ai, Hui-yuan Li, Zhong-qing Wang,