Optimized Schwarz methods for curl-curl time-harmonic Maxwell's equations

Dolean, Victorita and Gander, Martin J. and Lanteri, Stéphane and Lee, Jin-Fa and Peng, Zhen; Erhel, Jocelyne and Gander, Martin J. and Halpern, Laurence and Pichot, Géraldine and Sassi, Taoufik and Widlund, Olof, eds. (2014) Optimized Schwarz methods for curl-curl time-harmonic Maxwell's equations. In: Domain Decomposition Methods in Science and Engineering XXI. Lecture Notes in Computational Science and Engineering, 98 . Springer-Verlag, FRA, pp. 587-595. ISBN 9783319057897 (https://doi.org/10.1007/978-3-319-05789-7_56)

[thumbnail of Dolean-etal-DD21-2014-Optimized-schwarz-methods-for-curl-curl-time-harmonic-maxwell-equations]
Text. Filename: Dolean_etal_DD21_2014_Optimized_schwarz_methods_for_curl_curl_time_harmonic_maxwell_equations.pdf
Accepted Author Manuscript

Download (480kB)| Preview


Like the Helmholtz equation, the high frequency time-harmonic Maxwell's equations are difficult to solve by classical iterative methods. Domain decomposition methods are currently most promising: following the first provably convergent method in [4], various optimized Schwarz methods were developed over the last decade [1–3, 5, 8, 10, 11, 13, 14, 16]. There are however two basic formulations for Maxwell’s equation: the first order formulation, for which complete optimized results are known [5], and the second order, or curl-curl formulation, with partial optimization results [1, 13, 16]. We show in this paper that the convergence factors and the optimization process for the two formulations are the same. We then show by numerical experiments that the Fourier analysis predicts very well the behavior of the algorithms for a Yee scheme discretization, which corresponds to Nedelec edge elements on a tensor product mesh, in the curl-curl formulation. When using however mixed type Nedelec elements on an irregular tetrahedral mesh, numerical experiments indicate that transverse magnetic (TM) modes are less well resolved for high frequencies than transverse electric (TE) modes, and a heuristic can then be used to compensate for this in the optimization.