Numerical solution of Maxwell's equations in axisymmetric domains with the Fourier Singular Complement Method

Patrick Ciarlet (1) , Simon Labrunie (2, 3)

1 POEMS - Propagation des Ondes : Étude Mathématique et Simulation
2 IECN - Institut Élie Cartan de Nancy
3 CALVI - Scientific computation and visualization

We present an efficient method for computing numerically the solution to the time-dependent Maxwell equations in an axisymmetric domain, with arbitrary (not necessarily axisymmetric) data. The method is an extension of those introduced in~\cite{CJK+06} for Poisson's equation, and in~\cite{ACLS03} for Maxwell's equations in the fully axisymmetric setting (i.e., when the data is also axisymmetric). It is based on a Fourier expansion in the azimuthal direction, and on an improved variant of the Singular Complement Method in the meridian section. Thus, it can take into account the lack of regularity of the solution when the domain features non-convex edges or vertices. Moreover, it can handle noisy or approximate data which fail to satisfy the continuity equation, by using either an elliptic correction method or a mixed formulation. We give complete convergence analyses for both mixed and non-mixed formulations. Neither refinements near the reentrant edges or vertices of the domain, nor cut-off functions are required to achieve the desired convergence order in terms of the mesh size, the time step and the number of Fourier modes used.

Mots clés: Maxwell's equations; axisymmetric domains; finite element methods; mixed problems; singularities

Domaines: Equations aux dérivées partielles [math.AP]

Patrick Ciarlet, Simon Labrunie. Numerical solution of Maxwell's equations in axisymmetric domains with the Fourier Singular Complement Method. 2009. ⟨hal-00365391v1⟩