Modal analysis of a domain decomposition method for Maxwell's equations in a waveguide
For computational electromagnetics researchers, this work provides theoretical and numerical insights into the scalability of domain decomposition methods for Maxwell's equations, though it is incremental in extending existing analysis to more complex settings.
The paper analyzes the weak scalability of one-level Schwarz domain decomposition methods for Maxwell's equations in waveguides, providing a theoretical framework that extends to complex geometries and transmission conditions. Numerical experiments confirm that the limiting spectrum predicts practical behavior, and the method achieves robustness with respect to wave number under specific parameters.
Time-harmonic wave propagation problems, especially those governed by Maxwell's equations, pose significant computational challenges due to the non-self-adjoint nature of the operators and the large, non-Hermitian linear systems resulting from discretization. Domain decomposition methods, particularly one-level Schwarz methods, offer a promising framework to tackle these challenges, with recent advancements showing the potential for weak scalability under certain conditions. In this paper, we analyze the weak scalability of one-level Schwarz methods for Maxwell's equations in strip-wise domain decompositions, focusing on waveguides with general cross sections and different types of transmission conditions such as impedance or perfectly matched layers (PMLs). By combining techniques from the limiting spectrum analysis of Toeplitz matrices and the modal decomposition of Maxwell's solutions, we provide a novel theoretical framework that extends previous work to more complex geometries and transmission conditions. Numerical experiments confirm that the limiting spectrum effectively predicts practical behavior even with a modest number of subdomains. Furthermore, we demonstrate that the one-level Schwarz method can achieve robustness with respect to the wave number under specific domain decomposition parameters, offering new insights into its applicability for large-scale electromagnetic wave problems.