An Approximate Bounded Cochain Projection

Résumé Technique : Une Projection de Cochaîne Bornée Approximative

Énoncé du Problème
L'article traite de la construction d'un projecteur, $\Pi^{(k)}$, qui projette un complexe de Hilbert de formes différentielles de dimension infinie, $H\Lambda^{(k)}(\Omega)$, sur un sous-complexe de polynômes par morceaux de dimension finie, $H\Lambda^{(k)}_h(\Omega)$. L'objectif est de satisfaire trois propriétés spécifiques caractéristiques d'un « projecteur de cochaîne borné » :

1. Idempotence : $\Pi^{(k)}\Pi^{(k)} = \Pi^{(k)}$.
2. Bornitude Uniforme : Le projecteur est borné dans la norme de Sobolev ($H\Lambda^{(k)}$) avec une constante indépendante du paramètre de maillage $h$.
3. Propriété de Commutation : Le projecteur commute avec la dérivée extérieure ($d$), de telle sorte que $d\Pi^{(k)} = \Pi^{(k+1)}d$.

Bien que les projections de cochaînes bornées soient bien établies pour les domaines contractiles (où le complexe de de Rham est exact), l'article souligne une difficulté sur les domaines non contractiles. Sur ces domaines, le complexe possède une cohomologie non triviale (formes harmoniques). Les constructions standards échouent souvent à maintenir une bornitude uniforme ou une commutativité exacte lorsque le sous-espace de dimension finie ne capture pas parfaitement la topologie du domaine (c'est-à-dire lorsque l'écart entre les espaces harmoniques de dimension infinie et de dimension finie est non nul).

Méthodologie
Les auteurs proposent une construction basée sur une formulation mixte augmentée et un complexe discret à deux niveaux.

1. Formulation Mixte Augmentée :
   Au lieu de s'appuyer sur des opérateurs de lissage suivis d'une interpolation (ce qui peut manquer de bornitude uniforme dans la norme $H\Lambda$), les auteurs définissent le projecteur $\Pi^{(k)}$ via un problème variationnel. Étant donné $\phi^{(k)} \in H\Lambda^{(k)}(\Omega)$, le projecteur cherche $\lambda^{(i)} \in H\Lambda^{(i)}_h(\Omega)$ et une variable harmonique $h^{(k+1)} \in H^{(k+1)}_h$ satisfaisant un système d'équations symétriques.

   • Ce système introduit un multiplicateur de Lagrange, $h^{(k+1)}$, qui impose la contrainte que la dérivée projetée appartienne à l'espace des formes exactes.
   • Les auteurs prouvent la bijectivité (bien-posée) de ce système en utilisant des conditions de stabilité, incluant la coercivité de la forme bilinéaire sur les formes fermées et les conditions inf-sup pour la dérivée extérieure et l'espace harmonique.

2. Propriétés du Système Augmenté :

   • Idempotence : Prouvée en montrant que si l'entrée appartient déjà à l'espace de dimension finie, la solution du système mixte est l'entrée elle-même.
   • Bornitude : Prouvée en décomposant la solution en composantes fermées et orthogonales, en utilisant l'inégalité de Poincaré pour borner les normes $L^2$ et de graphe.
   • Propriété de Commutation : L'article démontre que le projecteur commute avec $d$ si le domaine est contractile ou si l'espace de dimension finie est « topologiquement conforme » (c'est-à-dire que l'écart entre les espaces harmoniques discrets et continus est nul). Si l'espace n'est pas topologiquement conforme, le multiplicateur harmonique $h^{(k+1)}$ est non nul, représentant l'incohérence topologique, et la propriété de commutation ne tient que jusqu'à l'ordre de grandeur de cette erreur.

3. Projection de Cochaîne Bornée Approximative (Construction à Deux Niveaux) :
   Pour traiter les domaines non contractiles où le sous-espace discret peut ne pas être topologiquement conforme, les auteurs proposent la composition de deux opérateurs :

   • Étape 1 : Un projecteur $\Pi^{(k)}$ projetant l'espace de dimension infinie vers un espace de dimension finie enrichi $\tilde{H}\Lambda^{(k)}_h(\Omega)$. Cet espace enrichi est construit pour être topologiquement conforme (écart = 0).
   • Étape 2 : Un projecteur canonique $\Pi^{(k)}_C$ projetant l'espace enrichi $\tilde{H}\Lambda^{(k)}_h(\Omega)$ vers l'espace de calcul cible $H\Lambda^{(k)}_h(\Omega)$.
     La composition $\Pi^{(k)}_C \Pi^{(k)}$ sert de projecteur de cochaîne bornée approximatif final.

Contributions Clés et Résultats

• Construction d'un Projecteur Borné : L'article fournit une construction rigoureuse d'un projecteur qui est idempotent et uniformément borné dans la norme de Sobolev pour les domaines contractiles et non contractiles.
• Traitement des Domaines Non Contractiles : Les auteurs montrent que sur les domaines non contractiles, la propriété de commutation est satisfaite exactement si l'espace harmonique discret est topologiquement conforme. Si ce n'est pas le cas, l'erreur de commutativité est quantifiée par le multiplicateur harmonique $h^{(k+1)}$, qui mesure la divergence topologique.
• Stabilité Discrète : Contrairement aux constructions précédentes qui reposent sur des régularisations continues, cette approche utilise une formulation mixte qui évite les difficultés analytiques associées aux inverses à droite continus de la dérivée extérieure.
• Commutativité Approximative : L'article démontre qu'en utilisant une approche à deux niveaux (enrichissement de l'espace pour résoudre la topologie avant de projeter vers le maillage final), on peut obtenir un projecteur qui est borné et qui commute avec la dérivée extérieure avec une précision arbitraire.

Signification et Revendications
L'article affirme présenter une construction de projecteur qui remplit les trois propriétés définissant une projection de cochaîne bornée. Plus précisément, il soutient que :

• Sur les domaines contractiles, le projecteur proposé est idempotent, uniformément borné et commute exactement avec la dérivée extérieure.
• Sur les domaines non contractiles, le projecteur reste idempotent et uniformément borné. La propriété de commutation est satisfaite « jusqu'à une précision arbitraire », ce qui signifie que l'erreur peut être rendue arbitrairement petite en affinant la résolution topologique de l'espace intermédiaire (par exemple, via le raffinement $p$ jusqu'à ce que le multiplicateur harmonique disparaisse à la précision machine).

Ce travail s'inscrit dans le cadre de l'Électro-Calcul de Formes Fines (FEEC - Finite Element Exterior Calculus), en affinant les idées existantes concernant les constructions locales et les conditions aux limites. Il se distingue en offrant une construction « entièrement locale » au sens où elle évite les inverses à droite continus, en s'appuyant plutôt sur des problèmes discrets posés sur un raffinement auxiliaire pour gérer les contraintes topologiques. Les auteurs positionnent cela comme une méthode pour assurer la stabilité et la cohérence des méthodes numériques pour les formes différentielles sur des géométries complexes.