Convergence of parallel overlapping domain decomposition methods with impedance boundary conditions for time-harmonic Maxwell equations in heterogeneous media

Résumé Technique : Convergence des méthodes de décomposition de domaine parallèles avec chevauchement pour les équations de Maxwell stationnaires

Énoncé du problème
Cet article traite de la résolution numérique des équations de Maxwell stationnaires dans des milieux hétérogènes, un problème central pour des applications telles que le radar, les communications sans fil et l'imagerie médicale. À haute fréquence, ces problèmes présentent des défis importants : la nécessité de maillages fins pour résoudre les oscillations conduit à des systèmes linéaires de grande taille, non hermitiens et hautement indéfinis. Bien que les méthodes de décomposition de domaine (DD) à deux niveaux aient établi des théories de convergence (reposant souvent sur des espaces grossiers qui résolvent la propagation des ondes), il existe un manque de théorie de convergence rigoureuse pour les méthodes DD à un seul niveau avec chevauchement pour les équations de Maxwell, particulièrement dans des contextes hétérogènes. Les auteurs notent que les arguments standards basés sur le champ de valeurs utilisés pour les méthodes à deux niveaux échouent pour les méthodes à un seul niveau considérées ici, car le champ de valeurs contient l'origine.

Méthodologie
Les auteurs analysent une méthode de Schwarz parallèle avec chevauchement utilisant des conditions aux limites d'impédance sur les interfaces des sous-domaines. La méthodologie se déroule en deux étapes principales : l'analyse continue et l'analyse discrète.

1. Analyse au niveau continu :

   • Bien posérescence : L'article établit la bien posérescence de l'itération de Schwarz parallèle dans des espaces de fonctions appropriés ($X_{imp}$ et $U$) pour des sous-domaines Lipschitziens généraux. Cela repose sur des résultats d'analyse harmonique concernant la régularité des solutions de Maxwell. Un obstacle technique majeur abordé est la régularité des fonctions de partition de l'unité et les termes d'erreur résultants, ce qui nécessite des hypothèses spécifiques (Hypothèse 2.11) concernant la géométrie des décompositions en bandes.
   • Propagation de l'erreur : L'opérateur de propagation de l'erreur est caractérisé à l'aide de cartes d'impédance à impédance. Ces cartes décrivent comment les données d'impédance (traces tangentielles du rotationnel et du champ lui-même) sur une interface de sous-domaine se propagent vers les interfaces voisines.
   • Décompositions en bandes : Pour les décompositions de domaine en bandes, les auteurs dérivent des estimations de convergence explicites. Ils montrent que la convergence de la $N$-ième puissance de l'opérateur de propagation de l'erreur dépend des normes de cartes d'impédance à impédance spécifiques (notées $\rho$ et $\gamma$), qui représentent les interactions « de gauche à droite » et « de gauche à gauche » (ou « de droite à droite »).

2. Analyse au niveau discret :

   • Discrétisation : La méthode continue est discrétisée à l'aide d'éléments finis de Nédélec. La méthode de Schwarz parallèle discrète est identifiée comme une méthode de Schwarz additive restreinte (RAS-imp) préconditionnant une itération de Richardson.
   • Propagation de l'erreur discrète : Les auteurs formulent un opérateur de propagation de l'erreur discret en termes de cartes d'impédance à impédance discrètes.
   • Héritage de la convergence : Sous l'hypothèse que les cartes d'impédance à impédance discrètes convergent vers leurs homologues continus lorsque la taille de la maille $h \to 0$, l'article prouve que la méthode discrète hérite du comportement de convergence de la méthode continue.

Principales contributions

• Cadre au niveau des EDP : Développement d'un cadre pour les méthodes de Schwarz parallèles avec chevauchement et conditions de transmission d'impédance pour les équations de Maxwell à coefficients variables. Cela inclut la preuve de la bien posérescence pour des sous-domaines Lipschitziens généraux et la dérivation d'équations d'erreur exprimées via des cartes d'impédance à impédance.
• Estimations de convergence pour les bandes : Dérivation de conditions suffisantes pour la convergence dans les décompositions en bandes. L'analyse réduit le problème de la borne de l'opérateur complet de propagation de l'erreur à la détermination des normes des cartes d'interface ($\rho$ et $\gamma$).
• Théorie de la convergence discrète : Démonstration que la discrétisation de Nédélec (RAS-imp) hérite des propriétés de convergence de la méthode continue, à condition que les cartes d'interface discrètes approximent les continues.
• Validation numérique : Expériences numériques approfondies validant la théorie. Celles-ci incluent :
  • La vérification de la contractivité de puissance de l'opérateur de propagation de l'erreur discret.
  • Le calcul des normes des cartes d'impédance à impédance, confirmant qu'elles sont bornées et cohérentes avec les prédictions théoriques.
  • Des tests de performance sur des décompositions en bandes et en damier (checkerboard) dans des milieux homogènes et hétérogènes, montrant que RAS-imp agit efficacement à la fois comme une itération stationnaire et comme un préconditionneur de Krylov (pour GMRES).

Résultats

• Bornes théoriques : L'article établit que pour les décompositions en bandes, l'opérateur de propagation de l'erreur $T^N$ devient contractif si les normes des cartes d'impédance à impédance sont suffisamment petites. Plus précisément, la norme de la $N$-ième puissance de l'opérateur est bornée par une fonction de $\rho$ et $\gamma$.
• Observations numériques :
  • Dans les décompositions en bandes, l'historique des résidus présente une « décroissance par paliers », correspondant à la contractivité de puissance où la norme de l'erreur chute significativement après environ $N$ itérations (où $N$ est le nombre de sous-domaines).
  • Les normes des cartes d'impédance à impédance ($\rho$ et $\gamma$) ont été calculées numériquement. $\rho$ (représentant la transmission à travers le domaine) s'est avéré être bien en dessous de 1, tandis que $\gamma$ (représentant la réflexion/interaction dans la même direction) est resté proche de 1 mais borné.
  • La méthode a démontré une robustesse vis-à-vis du nombre d'onde $k$ (pour un $N$ fixé) et une scalabilité vis-à-vis du nombre de sous-domaines.
  • Pour les décompositions en damier (qui dépassent le cadre théorique strict de l'analyse en bandes), la méthode a tout de même présenté une contractivité de puissance, bien que le seuil de contraction se soit déplacé vers environ $2\sqrt{N}$ itérations.
  • Dans les milieux hétérogènes, la convergence a été généralement plus lente que dans les cas homogènes, en raison des réflexions internes causées par la variation spatiale de la permittivité, mais la méthode est restée efficace.

Signification et affirmations
L'article affirme fournir la première théorie de convergence pour les méthodes de décomposition de domaine avec chevauchement à un seul niveau pour les équations de Maxwell stationnaires, tant au niveau continu que discret. Les auteurs soulignent que ce travail comble une lacune laissée par les théories existantes à deux niveaux, qui reposent souvent sur des espaces grossiers qui occultent le comportement de propagation crucial pour les méthodes à un seul niveau.

Les auteurs font preuve de modestie quant à l'étendue de leurs résultats par rapport aux travaux analogues sur l'équation de Helmholtz (cités comme [20, 21, 27]). Ils notent explicitement trois limitations :

1. Le résultat de bien posérescence pour l'itération continue nécessite une hypothèse de régularité (Hypothèse 2.11) qui n'est pas encore prouvée pour les sous-domaines Lipschitziens généraux dans le cas de Maxwell, contrairement au cas de Helmholtz.
2. La convergence des cartes d'impédance à impédance discrètes vers les continues est supposée (Hypothèse 6.5) plutôt que prouvée, en raison de la régularité réduite des solutions de Maxwell et de la complexité des arguments de dualité nécessaires.
3. L'article ne fournit pas de bornes de haute fréquence pour les cartes d'impédance à impédance au niveau continu, car la théorie de la mesure de défaut semi-classique requise pour ces bornes n'est pas encore établie pour les équations de Maxwell stationnaires.

Malgré ces limitations, l'article affirme que le cadre développé caractérise avec succès la propagation de l'erreur via les cartes d'interface et que les expériences numériques soutiennent fortement les prédictions théoriques, validant RAS-imp comme un solveur robuste pour ces problèmes électromagnétiques complexes.