The numerical solution of the three-dimensional time-harmonic Maxwell equations using high order methods such as discontinuous Galerkin formulations require efficient solvers. A domain decomposition strategy is introduced for this purpose. This strategy is based on optimized Schwarz methods applied to the first order form of the Maxwell system and leads to the best possible convergence of these algorithms. The principles are explained for a 2D model problem and numerical simulations confirm the predicted theoretical behavior. The efficiency is further demonstrated on more realistic 3D geometries including a bioelectromagnetism application.
- discontinuous Galerkin methods
- domain decomposition methods
- optimized interface conditions