Investigation on polynomial integrators for time-domain electromagnetics using a high-order discontinuous Galerkin method

Hassan Fahs (1)

1 XLIM-OSA - OSA

In this work, we investigate the application of polynomial integrators in a high-order discontinuous Galerkin method for solving the time-domain Maxwell equations. After the spatial discretization, we obtain a time-continuous system of ordinary differential equations of the form, ∂tY(t)=HY(t), where Y(t) is the electromagnetic field, H is a matrix containing the spatial derivatives, and t is the time variable. The formal solution is written as the exponential evolution operator, exp(tH), acting on a vector representing the initial condition of the electromagnetic field. The polynomial integrators are based on the approximation of exp(tH) by an expansion of the form ∑ _m=0^\infinity gm(t) Pm(H), where gm(t) is a function of time and Pm(H) is a polynomial of order m satisfying a short recursion. We introduce a general family of expansions of exp(tH) based on Faber polynomials. This family of expansions is suitable for non-Hermitian matrices, and consequently the proposed integrators can handle absorbing media and conductive materials. We discuss the efficient implementation of this technique, and based on some test problems, we compare the virtues and shortcomings of the algorithm. We also demonstrate how this scheme provides an efficient alternative to standard explicit integrators.

Mots clés: Maxwell's equations; Discontinuous Galerkin method; Time-stepping schemes; Faber polynomials

Domaines: Physique mathématique [math-ph]; Analyse numérique [math.NA]; Variables complexes [math.CV]; Analyse classique [math.CA]; Electromagnétisme; Modélisation et simulation

Hassan Fahs. Investigation on polynomial integrators for time-domain electromagnetics using a high-order discontinuous Galerkin method. Applied Mathematical Modelling, 2012, 36 (11), pp.Pages 5466--5481. ⟨10.1016/j.apm.2011.12.055⟩. ⟨hal-00722230⟩