Discontinuous Galerkin method for linear wave equations involving derivatives of the Dirac delta distribution
Linear wave equations sourced by a Dirac delta distribution \(\delta(x)\) and its derivative(s) can serve as a model for many different phenomena. We describe a discontinuous Galerkin (DG) method to numerically solve such equations with source terms proportional to \(\partial^n \delta /\partial x^n\...
Gespeichert in:
| Veröffentlicht in: | arXiv.org |
|---|---|
| Hauptverfasser: | , , , |
| Format: | Paper |
| Sprache: | Englisch |
| Veröffentlicht: |
Ithaca
Cornell University Library, arXiv.org
30.06.2023
|
| Schlagworte: | |
| ISSN: | 2331-8422 |
| Online-Zugang: | Volltext |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Zusammenfassung: | Linear wave equations sourced by a Dirac delta distribution \(\delta(x)\) and its derivative(s) can serve as a model for many different phenomena. We describe a discontinuous Galerkin (DG) method to numerically solve such equations with source terms proportional to \(\partial^n \delta /\partial x^n\). Despite the presence of singular source terms, which imply discontinuous or potentially singular solutions, our DG method achieves global spectral accuracy even at the source's location. Our DG method is developed for the wave equation written in fully first-order form. The first-order reduction is carried out using a distributional auxiliary variable that removes some of the source term's singular behavior. While this is helpful numerically, it gives rise to a distributional constraint. We show that a time-independent spurious solution can develop if the initial constraint violation is proportional to \(\delta(x)\). Numerical experiments verify this behavior and our scheme's convergence properties by comparing against exact solutions. |
|---|---|
| Bibliographie: | SourceType-Working Papers-1 ObjectType-Working Paper/Pre-Print-1 content type line 50 |
| ISSN: | 2331-8422 |
| DOI: | 10.48550/arxiv.2211.14390 |