Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem
An elliptic one-dimensional second order boundary value problem involving discontinuous coefficients, with or without transmission conditions, is considered. For the former case by a direct sum spaces method we show that the eigenvalues are real, geometrically simple and the eigenfunctions are orth...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Publishing House of the Romanian Academy
2019-12-01
|
| Series: | Journal of Numerical Analysis and Approximation Theory |
| Subjects: | |
| Online Access: | https://www.ictp.acad.ro/jnaat/journal/article/view/1201 |
| _version_ | 1811224413275684864 |
|---|---|
| author | Calin-Ioan Gheorghiu Bertin Zinsou |
| author_facet | Calin-Ioan Gheorghiu Bertin Zinsou |
| author_sort | Calin-Ioan Gheorghiu |
| collection | DOAJ |
| description |
An elliptic one-dimensional second order boundary value problem involving discontinuous coefficients, with or without transmission conditions, is considered. For the former case by a direct sum spaces method we show that the eigenvalues are real, geometrically simple and the eigenfunctions are orthogonal.
Then the eigenpairs are computed numerically by a local linear finite element method (FEM) and by some global spectral collocation methods.
The spectral collocation is based on Chebyshev polynomials (ChC) for problems on bounded intervals respectively on Fourier system (FsC) for periodic problems.
The numerical stability in computing eigenvalues is investigated by estimating their (relative) drift with respect to the order of approximation. The accuracy in computing the eigenvectors is addressed by estimating their departure from orthogonality as well as by the asymptotic order of convergence. The discontinuity of coefficients in the problems at hand reduces the exponential order of convergence, usual for any well designed spectral algorithm, to an algebraic one.
As expected, the accuracy of ChC outcomes overpasses by far that of FEM outcomes.
|
| first_indexed | 2024-04-12T08:48:58Z |
| format | Article |
| id | doaj.art-f9ecaaacea024f92ae06a9e58fb7c37c |
| institution | Directory Open Access Journal |
| issn | 2457-6794 2501-059X |
| language | English |
| last_indexed | 2024-04-12T08:48:58Z |
| publishDate | 2019-12-01 |
| publisher | Publishing House of the Romanian Academy |
| record_format | Article |
| series | Journal of Numerical Analysis and Approximation Theory |
| spelling | doaj.art-f9ecaaacea024f92ae06a9e58fb7c37c2022-12-22T03:39:38ZengPublishing House of the Romanian AcademyJournal of Numerical Analysis and Approximation Theory2457-67942501-059X2019-12-01482Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblemCalin-Ioan Gheorghiu0Bertin ZinsouTiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy An elliptic one-dimensional second order boundary value problem involving discontinuous coefficients, with or without transmission conditions, is considered. For the former case by a direct sum spaces method we show that the eigenvalues are real, geometrically simple and the eigenfunctions are orthogonal. Then the eigenpairs are computed numerically by a local linear finite element method (FEM) and by some global spectral collocation methods. The spectral collocation is based on Chebyshev polynomials (ChC) for problems on bounded intervals respectively on Fourier system (FsC) for periodic problems. The numerical stability in computing eigenvalues is investigated by estimating their (relative) drift with respect to the order of approximation. The accuracy in computing the eigenvectors is addressed by estimating their departure from orthogonality as well as by the asymptotic order of convergence. The discontinuity of coefficients in the problems at hand reduces the exponential order of convergence, usual for any well designed spectral algorithm, to an algebraic one. As expected, the accuracy of ChC outcomes overpasses by far that of FEM outcomes. https://www.ictp.acad.ro/jnaat/journal/article/view/1201Sturm-Liouville eigenproblemdiscontinuous coefficienttransmission conditionspectral collocationFEMFinite Element Method |
| spellingShingle | Calin-Ioan Gheorghiu Bertin Zinsou Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem Journal of Numerical Analysis and Approximation Theory Sturm-Liouville eigenproblem discontinuous coefficient transmission condition spectral collocation FEM Finite Element Method |
| title | Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem |
| title_full | Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem |
| title_fullStr | Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem |
| title_full_unstemmed | Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem |
| title_short | Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem |
| title_sort | analytic vs numerical solutions to a sturm liouville transmission eigenproblem |
| topic | Sturm-Liouville eigenproblem discontinuous coefficient transmission condition spectral collocation FEM Finite Element Method |
| url | https://www.ictp.acad.ro/jnaat/journal/article/view/1201 |
| work_keys_str_mv | AT calinioangheorghiu analyticvsnumericalsolutionstoasturmliouvilletransmissioneigenproblem AT bertinzinsou analyticvsnumericalsolutionstoasturmliouvilletransmissioneigenproblem |