78.3NAJun 2
Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshesSergio Gómez, Chiara Perinati, Paul Stocker et al.
We present and analyze an embedded Trefftz discontinuous Galerkin method for reaction-diffusion problems on anisotropic meshes. The method is constructed by imposing a relaxed local Trefftz condition via an embedding into a tensor-product DG space, yielding a reduced global system while preserving the approximation properties of the underlying high-order discretization. We prove stability and quasi-optimality on anisotropic, possibly curved, quadrilateral elements, and derive anisotropic a priori error estimates. Numerical experiments for $h$- and $hp$-refinement, including curved-domain examples, validate the theoretical results.
97.7NAApr 8
A discontinuous Galerkin method for elliptic-hyperbolic equationsChiara Perinati, Lise-Marie Imbert-Gérard, Andrea Moiola et al.
We present and analyze a discontinuous Galerkin method for the numerical solution of a class of second-order linear mixed-type partial differential equations, i.e. equations that change their nature from elliptic to hyperbolic through the computational domain. Well-posedness of the discrete problem is established via coercivity in an energy norm, achieved through the Morawetz multiplier technique. We derive $hp$-a priori error estimates in the energy norm, which we use to prove convergence rates for standard and quasi-Trefftz polynomial spaces. Numerical experiments validate the theoretical results.
35.0NAMar 13
Embedded Trefftz DG method for the Helmholtz equationPaul Stocker, Igor Voulis
We study an embedded Trefftz discontinuous Galerkin method for the Helmholtz equation. The method starts from a polynomial DG space and enforces the Trefftz property through local constraints, avoiding an explicit construction of Trefftz basis functions. For the global coupling we use a simple symmetric interior penalty DG bilinear form. Since the resulting formulation is not coercive, stability is proved by a $T$-coercivity argument combined with a Schatz-type duality technique. This yields wavenumber-explicit stability, quasi-optimality, and convergence estimates in standard DG norms under an explicit mesh resolution condition.