Orthogonal polynomials for the de Rham complex on the disk and cylinder

1. Problématique et Contexte

L'article s'attaque à la construction de bases de polynômes orthogonaux adaptés à la structure du complexe de de Rham sur un disque (2D) et un cylindre (3D), en respectant les conditions aux limites et la symétrie de rotation.

Le problème central réside dans la difficulté de discrétiser efficacement les espaces de Sobolev vectoriels ($\mathring{H}^1$, $\mathring{H}(\text{curl})$, $\mathring{H}(\text{div})$) nécessaires pour la résolution numérique d'équations aux dérivées partielles (EDP) comme l'équation de Maxwell ou l'équation de Stokes, dans le cadre du Calcul Extérieur des Éléments Finis (FEEC). Les bases classiques (comme les polynômes de Zernike scalaires) ne capturent pas naturellement la structure vectorielle, les conditions aux limites (normale ou tangentielle au bord), ni les relations de récurrence simples entre les opérateurs différentiels (gradient, rotationnel, divergence).

L'objectif est de créer des bases qui :

1. Décomposent le complexe de de Rham en sous-complexes exacts.
2. Respectent la symétrie de rotation (bases adaptées à la symétrie).
3. Offrent des relations de récurrence simples entre les gradients, rotationnels et divergences des fonctions de base.

2. Méthodologie

L'auteur développe une approche hiérarchique combinant la théorie des représentations, les polynômes orthogonaux multivariés et l'algèbre linéaire matricielle.

A. Bases Adaptées à la Symétrie (Symmetry-Adapted Bases)

Le travail repose sur la définition de fonctions invariantes et équivariantes sous l'action du groupe de rotation $SO(2)$.

• Fonctions scalaires : Une fonction $p$ est adaptée à la symétrie de mode $m$ si une rotation d'angle $\varphi$ multiplie la fonction par $e^{im\varphi}$.
• Fonctions vectorielles et matricielles : L'équivariance impose que le vecteur ou la matrice soit également transformé par la matrice de rotation $\rho(\varphi)$.
• Orthogonalité : Grâce au lemme de Schur (version simplifiée), les fonctions de modes différents sont automatiquement orthogonales par rapport à des poids invariants ou équivariants.

B. Construction des Polynômes

La construction procède en plusieurs étapes :

1. Polynômes homogènes : Construction de bases de polynômes homogènes scalaires, vectoriels et matriciels adaptés à la symétrie, utilisant des harmoniques complexes $y_m$ et une matrice équivariante $\Sigma(x,y)$.
2. Polynômes de Zernike Généralisés : Extension aux polynômes orthogonaux sur le disque avec un poids ultrasphérique $(1-r^2)^\lambda$.
3. Réduction à une dimension (Module $\mathcal{M}$) : Pour gérer les poids matriciels complexes, l'auteur relie les polynômes vectoriels sur le disque à une classe spéciale de polynômes vectoriels univariés dans un module $\mathcal{M}$, orthogonaux par rapport à des poids diagonaux.
4. Factorisation LU et Cholesky : Utilisation de factorisations matricielles (inspirées de Gautschi) pour construire explicitement des polynômes orthogonaux pour des poids modifiés, permettant d'obtenir des expressions fermées en termes de polynômes de Jacobi.

C. Le Poids $N(x,y)$ et la Condition Normale

Un élément clé est l'introduction d'une matrice de poids $N(x,y)$ qui force les vecteurs à être normaux au bord du disque ($r=1$).

• $N(x,y) = \begin{pmatrix} 1-y^2 & xy \\ xy & 1-x^2 \end{pmatrix}$.
• Les polynômes vectoriels orthogonaux par rapport à $N$ (notés $n^\nu_{mj}$) forment une base naturelle pour l'espace des vecteurs normaux au bord.

3. Contributions Clés

1. Construction de nouvelles bases orthogonales :

   • $w_{mj}$ : Polynômes scalaires pondérés (dérivés des Zernike généralisés) qui s'annulent au bord (base pour $\mathring{H}^1$).
   • $n^\nu_{mj}$ : Polynômes vectoriels orthogonaux par rapport à la matrice $N$, qui sont normaux au bord (base pour $\mathring{H}(\text{curl})$).
   • $z_{mj}$ : Polynômes de Zernike standards (base pour $L^2$).

2. Relations de récurrence simples entre opérateurs :
   L'article établit des relations explicites et simples reliant ces bases via les opérateurs différentiels :

   • Gradient : $\nabla w_{mj} = -(j+1) n^+_{m, j+1}$. Le gradient d'un polynôme scalaire pondéré est un multiple exact d'un polynôme vectoriel normal.
   • Rotationnel (Curl) : $\text{curl } n^-_{mj} = \kappa_{mj} z_{mj}$. Le rotationnel d'une combinaison spécifique des polynômes vectoriels donne un multiple exact des polynômes de Zernike scalaires.
   • Divergence : Des relations analogues sont établies pour les polynômes tangents $t_{mj}$ (obtenus par rotation de $n_{mj}$).

3. Décomposition du Complexe de de Rham :
   Ces bases permettent de décomposer le complexe de de Rham continu en une somme directe de sous-complexes exacts de petite dimension (souvent de longueur 3 ou 4), indexés par les modes de Fourier $m$ et les degrés $j$. Cela découple le problème global en de nombreux petits problèmes indépendants.

4. Extension aux Cylindres (3D) :
   La méthode est généralisée aux cylindres (périodiques et finis) en combinant les bases 2D avec des séries de Fourier (pour les cylindres périodiques) ou des polynômes de Legendre intégrés (pour les cylindres finis), préservant la structure du complexe de de Rham 3D.

4. Résultats Principaux

• Théorème 6 & 7 : Démonstration que le gradient des polynômes $w_{mj}$ et le rotationnel des polynômes $n^\nu_{mj}$ s'expriment linéairement et simplement en fonction des autres bases, avec des constantes de proportionnalité explicites.
• Exactitude : La séquence $0 \to \text{span}\{w\} \xrightarrow{\nabla} \text{span}\{n\} \xrightarrow{\text{curl}} \text{span}\{z\} \to 0$ est exacte pour chaque sous-espace de mode, reproduisant la topologie du disque (nombres de Betti $b_0=1, b_1=0, b_2=0$).
• Efficacité Computationnelle : La structure de blocs tridiagonales des matrices de récurrence (matrices de Jacobi) permet une implémentation efficace pour les solveurs spectraux ou $hp$-FEM. Les matrices de masse peuvent être déduites de ces récurrences.

5. Signification et Perspectives

• Applications Numériques : Ces bases offrent un cadre idéal pour la résolution spectrale ou $hp$-FEM d'équations de Maxwell, de Stokes et de l'élasticité sur des géométries cylindriques. Elles garantissent la stabilité et la conservation des propriétés topologiques (exactitude du complexe) au niveau discret.
• Complexité Optimale : L'utilisation de ces bases permet de construire des solveurs à complexité optimale pour les équations de type Maxwell dans des cylindres périodiques, y compris avec des coefficients variables équivariants.
• Extensions Futures :
  • Annulaires : L'extension aux anneaux (disques avec un trou) est mentionnée comme un défi, nécessitant peut-être des propriétés de bi-orthogonalité.
  • Sphères et Caps : Généralisation aux boules 3D et aux caps sphériques en utilisant des harmoniques sphériques et des représentations matricielles.
  • Complexes autres : Application potentielle au complexe de Koszul et au complexe d'élasticité.

En résumé, cet article fournit une construction mathématique rigoureuse et pratique de bases de polynômes qui « respectent » la géométrie et la topologie du problème, facilitant grandement la simulation numérique de phénomènes physiques régis par le calcul vectoriel sur des domaines à symétrie de rotation.