Résumé Technique : Advection par Éléments Spectraux Mimétiques

Énoncé du Problème

Le document traite de la discrétisation numérique de l'équation d'advection linéaire pour les formes différentielles sur des domaines bornés. Le problème classique d'advection scalaire, $\partial_t \rho + \nabla \cdot (v\rho) = 0$, est généralisé dans un cadre de géométrie différentielle sous la forme :
$$\frac{\partial \alpha^{(k)}}{\partial t} + \mathcal{L}_v \alpha^{(k)} = 0$$
où $\alpha^{(k)}$ est une $k$-forme différentielle (représentant des scalaires, des vecteurs ou des quantités de dimensions supérieures selon $k$), et $v$ est un champ de vitesse Lipschitz continu, et $\mathcal{L}_v$ est la dérivée de Lie. Le défi consiste à développer une discrétisation qui préserve les structures géométriques sous-jacentes des lois physiques, assure la conservation locale de la masse et atteint une précision spectrale tout en distinguant les quantités évaluées à des instants temporels spécifiques de celles intégrées sur des intervalles de temps.

Méthodologie

La méthode proposée, l'Advection par Éléments Spectraux Mimétiques, étend le cadre mimétique établi dans les travaux antérieurs [4] en incorporant la dérivée de Lie via la formule d'homotopie de Cartan. L'approche repose sur trois piliers fondamentaux :

1. Fondement de Géométrie Différentielle
La méthode utilise le langage des formes différentielles. La dérivée de Lie $\mathcal{L}_v$ est exprimée à l'aide de la formule d'homotopie de Cartan :
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
Un aspect crucial de la formulation est la relation d'adjoint entre le produit intérieur et le produit extérieur, définie via le produit scalaire $L^2$ :
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
Cette dualité permet de mapper les quantités physiques représentées par des produits intérieurs vers leurs 1-formes duales.

2. Discrétisation Spatiale (Éléments Spectraux Mimétiques)
Le domaine spatial est discrétisé à l'aide d'un complexe de cellules composé de points, de segments de ligne, de surfaces et de volumes. L'espace des $k$-formes lisses est approximé par un espace de dimension finie $\Lambda_h^k$ engendré par des formes de base $\epsilon_i^{(k)}$.

• Construction de la Base : Les fonctions de base sont construites à l'aide d'expansions polynomiales par morceaux (produits tensoriels) sur des éléments quadrilatéraux. En 1D, les 0-formes sont représentées par des polynômes de Lagrange évalués aux nœuds de Gauss-Lobatto, tandis que les 1-formes sont représentées par des "polynômes d'arête" dérivés pour satisfaire des propriétés intégrales spécifiques.
• Projection Commutante : Un opérateur de projection $\pi_h$ est défini de telle sorte qu'il commute avec la dérivée extérieure ($\pi_h d = d \pi_h$). Cela garantit que la dérivée extérieure discrète $d_h$ est représentée par des matrices d'incidence ne contenant que des valeurs $\{0, 1, -1\}$, préservant ainsi la structure topologique des opérateurs continus.

3. Discrétisation Temporelle (Intégrateur Mimétique Canonique)
L'intégration temporelle utilise un intégrateur symplectique d'ordre arbitraire dérivé des intégrateurs de collocation de Gauss canoniques [5].

• Nœuds Temporels Décalés : La méthode distingue deux ensembles de nœuds temporels :
  • Nœuds de Gauss-Lobatto ($t_k$) : Où les variables de la solution (par exemple, $\rho$) sont définies.
  • Nœuds de Gauss ($\tilde{t}_q$) : Où les dérivées temporelles (flux) sont évaluées.
• Évolution Discrète : Ce décalage aboutit à un intégrateur discret où le changement de la solution sur un intervalle est égal au flux évalué aux nœuds de Gauss internes. Cette structure reflète les méthodes de type leap-frog et la règle du point milieu implicite, assurant des propriétés symplectiques.

4. Produit Intérieur Discret
Le produit intérieur $\iota_v$ est discrétisé en imposant la dualité de couplage (Éq. 16) dans le cadre discret. Cela conduit à un système où les flux discrets sont calculés en résolvant un système linéaire impliquant les produits scalaires des formes de base et le champ de vitesse.

Contributions Clés

• Cohérence Géométrique : Le schéma incorpore explicitement la dérivée de Lie en utilisant la formule de Cartan, garantissant que la discrétisation respecte la représentation sans métrique des opérateurs différentiels et de leurs adjoints de Hilbert.
• Précision Spectrale : La méthode utilise des fonctions de base polynomiales d'ordre élevé, permettant une convergence spectrale dans l'espace.
• Conservation Locale de la Masse : L'utilisation de matrices d'incidence et la structure spécifique de la dérivée extérieure discrète garantissent la conservation locale de la masse.
• Distinction de l'Échelonnement Temporel : Le cadre sépare rigoureusement les quantités évaluées à des instants temporels de celles intégrées sur des intervalles, utilisant une grille temporelle décalée pour obtenir une intégration symplectique d'ordre élevé.

Résultats Numériques

Les auteurs présentent des expériences numériques sur des domaines 2D avec des conditions aux limites périodiques, testant l'advection d'ondes sinusoïdales et de cloches sinusoïdales dans des champs de vitesse constants et de vortex de Rudman.

• Convergence : La méthode démontre une convergence $h$-algébrique d'ordre $(p+1)$ et une convergence $p$-spectrale, à condition que l'erreur d'intégration temporelle ne domine pas l'erreur spatiale.
• Précision de l'Intégration Temporelle : Il est montré que l'erreur de la solution dépend de l'ordre du schéma d'intégration temporelle ($p_t$). Si le schéma temporel est suffisamment précis, l'erreur de discrétisation initiale est préservée ; sinon, une dispersion artificielle augmente avec le temps.
• Conservation de la Masse : L'erreur de masse totale reste à l'ordre de la précision machine pour les $10^3$ premiers pas de temps et reste inférieure à $10^{-12}$ même après $2 \times 10^4$ étapes, démontifiant d'excellentes propriétés de conservation.
• Réversibilité : Un test impliquant l'advection d'une onde sinusoïdale dans un vortex de Rudman, suivie d'une inversion de la direction du flux, montre que la méthode peut récupérer la solution initiale, prouvant la réversibilité de la méthode d'intégration.
• Dispersion Artificielle : L'article note que la dispersion artificielle est une fonction de l'ordre de l'intégration temporelle ; les schémas temporels d'ordre inférieur introduisent des erreurs de dispersion qui dépendent de la fréquence de l'onde advectée.

Signification et Revendications

L'article affirme que le schéma dérivé combine avec succès la précision spectrale avec la conservation locale de la masse au sein d'un cadre compatible avec la physique (mimétique). En étendant le cadre mimétique pour inclure la dérivée de Lie via la formule d'homotopie de Cartan, les auteurs fournissent une discrétisation qui clarifie la structure géométrique de l'équation d'advection. La méthode est présentée comme une approche robuste pour résoudre des problèmes d'advection où la préservation des propriétés géométriques et topologiques sous-jacentes des lois physiques est essentielle. Les auteurs notent modestement que, bien que la discrétisation spatiale soit hautement précise, la précision globale et les caractéristiques de dispersion dépendent de l'ordre du schéma d'intégration temporelle utilisé.