Projections, Embeddings and Stability
Fuente:
arXiv
Salvato in:
| Autore principale: | |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2024
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866914629927567360 |
|---|---|
| author | Olsson, Pelle |
| author_facet | Olsson, Pelle |
| contents | In the present work, we demonstrate how the pseudoinverse concept from linear algebra can be used to represent and analyze the boundary conditions of linear systems of partial differential equations. This approach has theoretical and practical implications; the theory applies even if the boundary operator is rank deficient, or near rank deficient. If desired, the pseudoinverse can be implemented directly using standard tools like Matlab. We also introduce a new and simplified version of the semidiscrete approximation of the linear PDE system, which completely avoids taking the time derivative of the boundary data. The stability results are valid for general, nondiagonal summation-by-parts norms. Another key result is the extension of summation-by-parts operators to multi-domains by means of carefully crafted embedding operators. No extra numerical boundary conditions are required at the grid interfaces. The aforementioned pseudoinverse allows for a compact representation of these multi-block operators, which preserves all relevant properties of the single-block operators. The embedding operators can be constructed for multiple space dimensions. Numerical results for the two-dimensional Maxwell's equations are presented, and they show very good agreement with theory. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2401_02197 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Projections, Embeddings and Stability Olsson, Pelle Numerical Analysis 65M06, 65M12 In the present work, we demonstrate how the pseudoinverse concept from linear algebra can be used to represent and analyze the boundary conditions of linear systems of partial differential equations. This approach has theoretical and practical implications; the theory applies even if the boundary operator is rank deficient, or near rank deficient. If desired, the pseudoinverse can be implemented directly using standard tools like Matlab. We also introduce a new and simplified version of the semidiscrete approximation of the linear PDE system, which completely avoids taking the time derivative of the boundary data. The stability results are valid for general, nondiagonal summation-by-parts norms. Another key result is the extension of summation-by-parts operators to multi-domains by means of carefully crafted embedding operators. No extra numerical boundary conditions are required at the grid interfaces. The aforementioned pseudoinverse allows for a compact representation of these multi-block operators, which preserves all relevant properties of the single-block operators. The embedding operators can be constructed for multiple space dimensions. Numerical results for the two-dimensional Maxwell's equations are presented, and they show very good agreement with theory. |
| title | Projections, Embeddings and Stability |
| topic | Numerical Analysis 65M06, 65M12 |
| url | https://arxiv.org/abs/2401.02197 |