Abstract
We investigate the finite element-Galerkin method for singular self-adjoint second-order differential expressions. The weak formulation of the problem involves integration by parts, which allows the use of the usual piecewise linear functions. Our analysis shows that the method produces the solution corresponding to a particular self-adjoint realization of the differential expression. We also propose two algorithms to approximate the solution of any self-adjoint realization. Numerical examples are given to illustrate the analysis as well as the proposed algorithms.
| Original language | English |
|---|---|
| Pages (from-to) | 735-752 |
| Number of pages | 18 |
| Journal | Journal of Computational and Applied Mathematics |
| Volume | 223 |
| Issue number | 2 |
| DOIs | |
| State | Published - 15 Jan 2009 |
Bibliographical note
Funding Information:We would like to thank the referee for many suggestions which led to improvements in the paper. Research of the first two authors has been funded by King Fahd University of Petroleum and Minerals under Project number MS/Singular ODE/274.
Keywords
- Finite element method
- Galerkin method
- Self-adjoint operators
- Singular differential operators
- Weak formulation
ASJC Scopus subject areas
- Computational Mathematics
- Applied Mathematics