Polynomial quasi-Trefftz DG for PDEs with smooth coefficients: elliptic problems
Fuente:
arXiv
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2024
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866910045394960384 |
|---|---|
| author | Imbert-Gérard, Lise-Marie Moiola, Andrea Perinati, Chiara Stocker, Paul |
| author_facet | Imbert-Gérard, Lise-Marie Moiola, Andrea Perinati, Chiara Stocker, Paul |
| contents | Trefftz schemes are high-order Galerkin methods whose discrete spaces are made of elementwise exact solutions of the underlying PDE. Trefftz basis functions can be easily computed for many PDEs that are linear, homogeneous, and have piecewise-constant coefficients. However, if the equation has variable coefficients, exact solutions are generally unavailable. Quasi-Trefftz methods overcome this limitation relying on elementwise "approximate solutions" of the PDE, in the sense of Taylor polynomials.
We define polynomial quasi-Trefftz spaces for general linear PDEs with smooth coefficients and source term, describe their approximation properties and, under a non-degeneracy condition, provide a simple algorithm to compute a basis. We then focus on a quasi-Trefftz DG method for variable-coefficient elliptic diffusion-advection-reaction problems, showing stability and high-order convergence of the scheme. The main advantage over standard DG schemes is the higher accuracy for comparable numbers of degrees of freedom. For non-homogeneous problems with piecewise-smooth source term we propose to construct a local quasi-Trefftz particular solution and then solve for the difference. Numerical experiments in 2 and 3 space dimensions show the excellent properties of the method both in diffusion-dominated and advection-dominated problems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2408_00392 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Polynomial quasi-Trefftz DG for PDEs with smooth coefficients: elliptic problems Imbert-Gérard, Lise-Marie Moiola, Andrea Perinati, Chiara Stocker, Paul Numerical Analysis 65N15, 65N30, 35J25, 41A10, 41A25 Trefftz schemes are high-order Galerkin methods whose discrete spaces are made of elementwise exact solutions of the underlying PDE. Trefftz basis functions can be easily computed for many PDEs that are linear, homogeneous, and have piecewise-constant coefficients. However, if the equation has variable coefficients, exact solutions are generally unavailable. Quasi-Trefftz methods overcome this limitation relying on elementwise "approximate solutions" of the PDE, in the sense of Taylor polynomials. We define polynomial quasi-Trefftz spaces for general linear PDEs with smooth coefficients and source term, describe their approximation properties and, under a non-degeneracy condition, provide a simple algorithm to compute a basis. We then focus on a quasi-Trefftz DG method for variable-coefficient elliptic diffusion-advection-reaction problems, showing stability and high-order convergence of the scheme. The main advantage over standard DG schemes is the higher accuracy for comparable numbers of degrees of freedom. For non-homogeneous problems with piecewise-smooth source term we propose to construct a local quasi-Trefftz particular solution and then solve for the difference. Numerical experiments in 2 and 3 space dimensions show the excellent properties of the method both in diffusion-dominated and advection-dominated problems. |
| title | Polynomial quasi-Trefftz DG for PDEs with smooth coefficients: elliptic problems |
| topic | Numerical Analysis 65N15, 65N30, 35J25, 41A10, 41A25 |
| url | https://arxiv.org/abs/2408.00392 |