We propose a Trefftz discontinuous Galerkin (TDG) method for the approximation of plane wave scattering by periodic diffraction gratings, modelled by the two-dimensional Helmholtz equation. The periodic obstacle may include penetrable and impenetrable regions. The TDG method requires the approximation of the Dirichlet-to-Neumann (DtN) operator on the periodic cell faces, and relies on plane wave discrete spaces. For polygonal meshes, all linear-system entries can be computed analytically. Using a Rellich identity, we prove a new explicit stability estimate for the Helmholtz solution, which is robust in the small material jump limit.
Comment: 30 pages, 16 figures