Dispersion and dissipation error in high-order Runge-Kutta discontinuous Galerkin discretisations of the Maxwell equations

Different time-stepping methods for a nodal high-order discontinuous Galerkin discretisation of the Maxwell equations are discussed. A comparison between the most popular choices of Runge-Kutta (RK) methods is made from the point of view of accuracy and computational work. By choosing the strong-stability-preserving Runge-Kutta (SSP-RK) time-integration method of order consistent with the polynomial order of the spatial discretisation, better accuracy can be attained compared with fixed-order schemes. Moreover, this comes without a significant increase in the computational work. A numerical Fourier analysis is performed for this Runge-Kutta discontinuous Galerkin (RKDG) discretisation to gain insight into the dispersion and dissipation properties of the fully discrete scheme. The analysis is carried out on both the one-dimensional and the two-dimensional fully discrete schemes and, in the latter case, on uniform as well as on non-uniform meshes. It also provides practical information on the convergence of the dissipation and dispersion error up to polynomial order 10 for the one-dimensional fully discrete scheme

Finite element eigenvalue enclosures for the Maxwell operator

We propose employing the extension of the Lehmann-Maehly-Goerisch method developed by Zimmermann and Mertins, as a highly effective tool for the pollution-free finite element computation of the eigenfrequencies of the resonant cavity problem on a bounded region. This method gives complementary bounds for the eigenfrequencies which are adjacent to a given real parameter. We present a concrete numerical scheme which provides certified enclosures in a suitable asymptotic regime. We illustrate the applicability of this scheme by means of some numerical experiments on benchmark data using Lagrange elements and unstructured meshes.Comment: arXiv admin note: substantial text overlap with arXiv:1306.535

A nodal-based finite element approximation of the Maxwell problem suitable for singular solutions

A new mixed finite element approximation of Maxwell’s problem is proposed, its main features being that it is based on a novel augmented formulation of the continuous problem and the introduction of a mesh dependent stabilizing term, which yields a very weak control on the divergence of the unknown. The method is shown to be stable and convergent in the natural H(curl; ) norm for this unknown. In particular, convergence also applies to singular solutions, for which classical nodal based interpolations are known to suffer from spurious convergence upon mesh refinement.Postprint (published version

Nonstandard Finite Element Methods

[no abstract available