Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method

Résumé technique : Couplage discret FEM-BEM avec la méthode de Schwarz optimisée généralisée

Énoncé du problème
L'article traite de la résolution numérique des problèmes de propagation d'ondes acoustiques modélisés par l'équation de Helmholtz sur des domaines bornés et non bornés. L'objectif est de développer un solveur polyvalent et efficace capable de gérer des matériaux hétérogènes et des géométries complexes, incluant celles présentant des points de croisement (points partagés par plus de deux sous-domaines). Les auteurs se concentrent sur le couplage de la méthode des éléments finis (FEM), adaptée aux régions bornées hétérogènes, avec la méthode des éléments de frontière (BEM), optimale pour les domaines par morceaux homogènes et non bornés. Un défi spécifique abordé est l'éventuelle malposée des formulations classiques FEM-BEM aux « résonances parasites » (nombres d'onde où la formulation variationnelle admet un noyau non trivial malgré le fait que le problème physique soit bien posé).

Méthodologie
Les auteurs proposent une approche de décomposition de domaine (DD) non chevauchante basée sur la méthode de Schwarz optimisée généralisée (GOSM). La méthodologie comprend :

1. Substructuration : Le domaine est partitionné en sous-domaines non chevauchants. Les inconnues du solveur itératif sont associées aux interfaces (le squelette) plutôt qu'au volume.
2. Cadre discret : L'analyse étend la théorie continue présentée dans [20] à un cadre entièrement discret utilisant des approximations de Galerkin (FEM et BEM).
3. Formulation par opérateur : La méthode reformule le problème global en une équation de sous-structure de la forme $(Id + \Pi S)q = \text{rhs}$, où :
   • $S$ est un opérateur de diffusion qui résout des problèmes locaux avec des conditions aux limites d'impédance optimisées.
   • $\Pi$ est un opérateur d'échange non local (remplaçant l'opérateur d'échange standard utilisé dans les méthodes de Schwarz optimisées classiques) conçu pour gérer les points de croisement.
   • $q$ représente les traces d'impédance sur le squelette.
4. Couplages FEM-BEM : L'étude analyse spécifiquement trois couplages FEM-BEM classiques : Costabel, Johnson-Nédélec et Bielak-MacCamy. Elle examine l'« hypothèse de signe » (liée à la dissipation physique) pour les opérateurs de bord $A_\Gamma$ associés à ces couplages.
5. Solveur itératif : Lorsque l'hypothèse de signe est vérifiée (spécifiquement pour le couplage Costabel), l'opérateur $\Pi S$ est une contraction. Cela permet l'utilisation d'un schéma itératif de Richardson avec une convergence géométrique garantie. Pour les couplages où l'hypothèse de signe ne tient pas (Johnson-Nédélec, Bielak-MacCamy), un solveur GMRes est proposé.
6. Opérateurs de transmission : Les expériences numériques utilisent divers opérateurs de transmission, notamment des opérateurs locaux de type Després, des opérateurs hypersinguliers de type Yukawa non locaux, et des opérateurs de complément de Schur dérivés d'une application de Dirichlet-vers-Neumann positive.

Principales contributions

• Bien poséeté discrète : L'article établit que les formulations de sous-structure des couplages FEM-BEM classiques sont bien posées pour des nombres d'onde différents des résonances parasites classiques.
• Relation de dimension du noyau : Une relation explicite est établie entre les dimensions des noyaux de la formulation variationnelle initiale, des problèmes locaux et de la formulation de sous-structure (Théorème 7.3).
• Robustesse aux résonances parasites : Un résultat théorique important est la preuve que la formulation de sous-structure du couplage FEM-BEM de Costabel reste bien posée même aux résonances parasites. Ceci est dû au fait que la dimension du noyau de l'opérateur de sous-structure correspond à celle du problème local, permettant un noyau trivial dans le système de sous-structure sous certaines conditions.
• Théorie de la convergence : Les auteurs prouvent que pour le couplage Costabel (qui satisfait l'hypothèse de signe), l'opérateur de sous-structure prend la forme « identité + contraction ». Cela garantit la convergence géométrique du solveur de Richardson.
• Gestion des points de croisement : La méthode est démontrée comme étant robuste aux points de croisement sans nécessiter d'inconnues supplémentaires par point de croisement ou le couplage de problèmes locaux en ces points, une difficulté courante dans d'autres méthodes de DD.
• Robustesse en $h$ : Des estimations théoriques et des expériences numériques démontrent que la convergence géométrique peut être obtenue de manière robuste vis-à-vis de la taille de maillage ($h$-robustesse), à condition d'utiliser des opérateurs de transmission non locaux (comme les types Yukawa ou complément de Schur).

Résultats
Des expériences numériques ont été menées en 2D pour des problèmes homogènes et hétérogènes, ainsi que pour des problèmes impliquant des obstacles impénétrables et des cavités ouvertes (introduisant des points de croisement).

• Nombre d'itérations : Le nombre d'itérations requis pour la convergence a été analysé par rapport à la taille de maillage ($h$) et au nombre d'onde ($\kappa$).
  • Les configurations utilisant des opérateurs de transmission locaux (Després) ont montré une croissance quasi-quadratique des itérations lors du raffinement du maillage ($O(h^{1.69} \dots O(h^{1.78})$).
  • Les configurations utilisant des opérateurs de transmission non locaux (Yukawa, complément de Schur) ont présenté des nombres d'itérations quasi-constants par rapport au raffinement du maillage, confirmant la robustesse en $h$.
  • Par rapport au nombre d'onde, les configurations non locales ont généralement montré une croissance quasi-constante ou sous-linéaire, tandis que les configurations locales ont montré une croissance sous-linéaire.
• Comparaison des couplages : Bien que les garanties théoriques de convergence géométrique aient été principalement dérivées pour le couplage Costabel, les résultats numériques indiquent que les couplages Johnson-Nédélec et Bielak-MacCamy produisent des nombres d'itérations similaires dans les scénarios testés.
• Robustesse aux points de croisement : Dans le cas de test de la cavité rectangulaire impliquant des points de croisement, la méthode a maintenu sa robustesse en $h$, le nombre d'itérations restant stable lors du raffinement du maillage.
• Sélection de l'opérateur de transmission : La configuration « Y-S » (opérateur Yukawa pour BEM, complément de Schur pour FEM) a généralement produit le plus faible nombre d'itérations. L'étude souligne qu'il est bénéfique d'adapter le type d'opérateur de transmission à la nature du sous-domaine (creux vs dense) pour l'efficacité de l'implémentation.

Signification et affirmations
L'article affirme fournir le premier algorithme itératif $h$-robuste pour le couplage FEM-BEM avec des bornes de convergence garanties. En reformulant les couplages FEM-BEM classiques en une forme de sous-structure, les auteurs démontrent que ces méthodes classiques peuvent être résolues efficacement à l'aide de techniques de décomposition de domaine, même en présence de points de croisement et de résonances parasites. Ce travail comble le fossé entre la théorie de la décomposition de domaine continue et les implémentations discrètes FEM-BEM, offrant une stratégie systématique pour coupler ces méthodes qui soit à la fois théoriquement solide et numériquement efficace pour les applications de calcul haute performance. Les auteurs soulignent que leur approche évite la complexité du couplage des problèmes locaux aux points de croisement, ce qui constitue un avantage distinct par rapport aux autres formulations FEM-BEM multi-domaines.