Monolithic Multi-level Overlapping Schwarz Solvers for Fluid Problems

Résumé technique : Solveurs de Schwarz monolithiques à niveaux multiples avec recouvrement pour les problèmes de fluides

Énoncé du problème
L'article traite de la simulation numérique de la dynamique des fluides incompressibles, en particulier la résolution des équations de Stokes discrétisées, qui aboutissent à des systèmes de points de selle à grande échelle. À mesure que la complexité géométrique augmente et que les résolutions de maillage deviennent plus fines pour résoudre les couches limites, les solveurs directs deviennent irréalisables en raison de coûts de calcul prohibitifs et d'une consommation mémoire excessive. Les auteurs se concentrent sur les limitations d'évolutivité des méthodes de décomposition de domaine standard à deux niveaux, en particulier l'approche généralisée Dryja–Smith–Widlund (GDSW). Dans un schéma à deux niveaux, à mesure que le nombre de sous-domaines augmente, le problème grossier ($A_0$) augmente en taille, rendant finalement sa factorisation directe computationnellement intraitable.

Méthodologie
Les auteurs proposent et implémentent un préconditionneur de Schwarz monolithique à niveaux multiples avec recouvrement pour surmonter les limites d'évolutivité de l'approche GDSW à deux niveaux.

• Formulation mathématique : La méthode étend le préconditionneur GDSW à deux niveaux en l'appliquant de manière récursive. Au lieu d'inverser directement le problème grossier $A_0$, les auteurs remplacent $A_0^{-1}$ par une approximation dérivée d'un niveau plus grossier (une hiérarchie à trois niveaux ou à niveaux multiples). Le préconditionneur est défini comme suit :
  $$B_{GDSW, 3-level} = \Phi^{(2)} \left( \Phi^{(3)} A_0^{(3)-1} \Phi^{(3)T} + \sum_{i=1}^{N^{(2)}} R_i^{(2)T} A_i^{(2)-1} R_i^{(2)} \right) \Phi^{(2)T} + \sum_{i=1}^{N^{(1)}} R_i^{(1)T} A_i^{(1)-1} R_i^{(1)}$$
  Cette structure permet de résoudre le problème grossier via un niveau supplémentaire de décomposition de domaine plutôt que par un solveur direct.

• Implémentation : Le solveur est implémenté dans la bibliothèque FROSch (Fast and Robust Overlapping Schwarz), faisant partie du package ShyLU dans l'écosystème Trilinos.

  • Couplage : Une interface évolutive connecte FROSch à la bibliothèque d'éléments finis FEATFLOW. FEATFLOW gère l'assemblage de la matrice du système et des vecteurs du second membre en utilisant des éléments $Q_2$ pour la vitesse et des éléments $P_1$ discontinus pour la pression.
  • Nature algébrique : L'implémentation est purement algébrique, ne nécessitant que des informations sur la distribution parallèle des degrés de liberté (dofs). Elle construit automatiquement les matrices locales de sous-domaines, les opérateurs de restriction/interpolation et le problème grossier.
  • Construction de l'espace grossier : Pour la variable de pression, qui est discontinue à travers les sous-domaines, les auteurs utilisent des fonctions constantes sur tout le sous-domaine comme fonctions de base grossières, évitant ainsi le besoin d'une interface de pression $\Gamma_p$. Pour la vitesse, des extensions minimisant l'énergie des restrictions de l'espace nul sont utilisées.

Contributions clés

1. Extension aux niveaux multiples : Le travail étend avec succès le préconditionneur GDSW monolithique à deux niveaux (introduit précédemment dans [7]) à un cadre à trois niveaux (et potentiellement à niveaux multiples) spécifiquement pour les problèmes de fluides de type point de selle.
2. Couplage évolutif : Le développement d'une interface entre la bibliothèque FEATFLOW et le cadre Trilinos/FROSch, permettant un couplage efficace de l'assemblage géométrique avec des solveurs algébriques de décomposition de domaine.
3. Résultats parallèles à grande échelle : Les auteurs démontrent des résultats d'évolutivité faible sur des architectures de calcul haute performance (HPC), atteignant une évolutivité jusqu'à 32 768 rangs MPI pour un écoulement de Poiseuille sur un cube unité et jusqu'à 4 096 rangs MPI pour une géométrie complexe de filière d'extrusion.

Résultats numériques
Les auteurs ont évalué le solveur en utilisant deux cas de test : un écoulement de Poiseuille sur un cube unité et un écoulement de Poiseuille sur une géométrie complexe de filière d'extrusion.

• Cube unité (décomposition structurée) :

  • Solveur à deux niveaux : À 32 768 sous-domaines, le solveur à deux niveaux a échoué car le problème grossier (407 837 degrés de liberté) ne pouvait pas être factorisé. Le temps de configuration pour la factorisation du problème grossier est devenu le goulot d'étranglement.
  • Solveur à trois niveaux : A résolu avec succès le problème à 32 768 sous-domaines. Bien que le nombre d'itérations de Krylov ait augmenté de 41 à 66 alors que le nombre de sous-domaines augmentait d'un facteur 64, le solveur est resté robuste. L'approche à trois niveaux a considérablement réduit le temps de configuration par rapport à l'approche à deux niveaux pour un grand nombre de sous-domaines, évitant la factorisation irréalisable du problème grossier massif.

• Filière d'extrusion (décomposition non structurée) :

  • Solveur à deux niveaux : À 4 096 sous-domaines, le temps de configuration était dominé par la factorisation du problème grossier (84,3 s).
  • Solveur à trois niveaux : A démontré un gain de performance. Le temps de configuration numérique a diminué d'un facteur d'environ 6 par rapport au solveur à deux niveaux, et le temps de résolution a diminué d'un facteur presque 2. L'approche à trois niveaux a géré avec succès la décomposition non structurée là où l'approche à deux niveaux rencontrait des limites d'évolutivité.

Signification et affirmations
L'article affirme que la signification principale de ce travail est la démonstration réussie que les préconditionneurs de Schwarz monolithiques à niveaux multiples avec recouvrement peuvent surmonter les goulots d'étranglement d'évolutivité inhérents aux méthodes GDSW à deux niveaux pour les problèmes de fluides. Plus précisément :

• L'approche à trois niveaux permet de résoudre des problèmes avec un nombre significativement plus élevé de sous-domaines (jusqu'à 32 768) où le problème grossier à deux niveaux serait autrement irréalisable à factoriser.
• La méthode fournit un gain de performance à la fois dans les temps de configuration et de résolution pour les problèmes à grande échelle, en particulier lorsque l'espace grossier devient grand.
• Le travail représente une étape vers des simulations de dynamique des fluides computationnelle (CFD) efficaces sur les futures architectures exascale avec des systèmes hétérogènes CPU-GPU, dans le cadre du projet allemand "StroemungsRaum".

Les auteurs restent modestes concernant l'augmentation des itérations de Krylov à des nombres de sous-domaines très élevés, notant que le préconditionneur n'est peut-être pas encore dans sa plage asymptotique, mais soulignent que la croissance des itérations est acceptable étant donné l'augmentation massive du parallélisme.