Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

Résumé technique : Analyse d'erreur d'un schéma aux éléments finis mixtes préservant la divergence pour les équations de Hall–MHD incompressibles

Énoncé du problème
Le système de magnétohydrodynamique à effet Hall (Hall–MHD) incompressible modélise l'interaction entre le mouvement des fluides et les champs magnétiques dans les fluides électriquement conducteurs, en intégrant l'effet Hall qui devient significatif à des échelles comparables à la profondeur de peau ionique. Cet effet introduit un terme raide et hautement non linéaire ($\eta \text{curl}(\text{curl} B \times B)$) dans l'équation d'induction, conduisant à un comportement dispersif, à des gradients raides et à des oscillations haute fréquence (ondes siffleuses). Ces caractéristiques posent des défis analytiques et computationnels substantiels. Plus précisément, l'existence globale de solutions fortes en trois dimensions reste un problème ouvert, et les méthodes numériques standard éprouvent souvent des difficultés de stabilité et face à l'exigence stricte que le champ magnétique reste solénoïdal ($\text{div } B = 0$). Les violations de cette condition de divergence nulle peuvent induire des monopôles magnétiques non physiques et déclencher des instabilités numériques sévères. De plus, les approches numériques existantes, telles que celles basées sur les méthodes de Variable Auxiliaire Scalaire (SAV) ou les techniques spectrales, dissipent souvent des fonctionnels d'énergie modifiés plutôt que la véritable énergie physique, échouent à imposer la condition de divergence nulle exactement au niveau discret, ou manquent d'une analyse d'erreur rigoureuse et de flexibilité géométrique.

Méthodologie
Pour répondre à ces défis, les auteurs proposent une méthode aux éléments finis mixtes, linéaire, entièrement discrète et préservant la structure, pour un système de Hall–MHD régularisé de Voigt. La régularisation de Voigt remplace la dérivée temporelle dans l'équation de la quantité de mouvement par une version filtrée spatialement $(I - \alpha_1 \Delta)\partial_t u$ et conserve un terme d'inertie électronique $\alpha_2 \partial_t J$ dans la loi d'Ohm généralisée. Cette modification améliore le bien-posé et la stabilité numérique tout en préservant la structure physique essentielle et les lois d'énergie du système original.

Le schéma numérique est construit en utilisant les composants suivants :

1. Formulation mixte : Le problème est reformulé en utilisant le champ électrique $E$ et la densité de courant $J$ comme variables primaires, aboutissant à un système couplé d'équations de quantité de mouvement, d'induction, de loi d'Ohm et de loi d'Ampère.
2. Discrétisation spatiale : Le schéma utilise le Calcul Extérieur aux Éléments Finis (FEEC) pour garantir que la condition de divergence nulle est satisfaite exactement.
   • La vitesse $u$ et la pression $p$ sont approchées à l'aide de la paire d'éléments finis MINI (inf-sup stable d'ordre le plus bas).
   • Le champ magnétique $B$ et la densité de courant $J$ sont approchés respectivement à l'aide des espaces de Raviart–Thomas ($RT_0$) et de Nédélec ($X_0$) d'ordre le plus bas.
   • Ces espaces forment un complexe de de Rham discret, garantissant que $\text{div } B_h = 0$ est préservé exactement au niveau discret.
3. Discrétisation temporelle : Un schéma de pas de temps linéairement implicite et antisymétrique est utilisé. Les termes convectifs et de Lorentz sont traités semi-implicitement en utilisant les valeurs de l'étape de temps précédente pour maintenir la linéarité et la stabilité énergétique inconditionnelle.
4. Analyse d'erreur : Les auteurs établissent des taux de convergence optimaux en décomposant l'erreur en erreurs d'approximation (en utilisant des projections de Stokes et des projections orthogonales) et en erreurs discrètes. Ils utilisent des inégalités de Gaffney discrètes et des estimations d'énergie pour traiter les termes de couplage non linéaires.

Contributions clés

• Schéma préservant la structure : L'article présente un FEM linéaire et entièrement discret qui impose la contrainte solénoïdale sur le champ magnétique exactement, évitant le besoin d'étapes de projection ou de méthodes de pénalisation qui pourraient dégrader la précision ou la stabilité.
• Stabilité énergétique inconditionnelle : Le schéma proposé satisfait une loi de dissipation d'énergie discrète de manière inconditionnelle, reflétant la loi d'énergie physique continue, même en présence du terme de Hall raide.
• Estimations d'erreur rigoureuses : Les auteurs fournissent une analyse d'erreur complète pour le problème régularisé de Voigt, établissant des taux de convergence optimaux de l'ordre $O(h^s + \tau)$ sous des hypothèses de régularité appropriées.
• Extension au système non régularisé : Crucialement, l'analyse s'étend au système de Hall–MHD non régularisé (où $\alpha_1 = \alpha_2 = 0$). Les auteurs démontrent que la convergence peut toujours être prouvée pour le système non régularisé, à condition que la résistivité $\sigma > 0$. Dans le régime non résistif ($\sigma = 0$), la régularisation de Voigt s'avère jouer un rôle crucial dans les estimations d'erreur.
• Implémentation 2,5D et 3D : Le cadre est validé à la fois dans des configurations 2,5D (invariance translationnelle) et pleinement 3D, démontrant l'applicabilité de la méthode à des géométries complexes et à des phénomènes physiques.

Résultats
Les simulations numériques en 2,5D et 3D confirment les résultats théoriques :

• Convergence : Le schéma atteint les taux de convergence spatiaux et temporels prédits pour la vitesse, le champ magnétique et la densité de courant. Bien que les taux de convergence pour la vitesse et la densité de courant soient légèrement réduits dans certains cas (attribués à une régularité de solution limitée due au terme de Hall), le champ magnétique maintient une convergence optimale.
• Préservation de la divergence : L'erreur de divergence maximale ($\max | \text{div } B_h |$) reste au niveau de la tolérance du solveur tout au long des simulations, confirmant la préservation exacte de la contrainte solénoïdale.
• Stabilité énergétique : L'énergie totale du système décroît de manière monotone au cours du temps, conformément à la loi de dissipation d'énergie théorique.
• Fidélité physique : Les simulations de l'écoulement d'Arnold–Beltrami–Childress (ABC), du vortex d'Orszag–Tang et de la reconnexion pilotée par l'effet Hall dans une feuille de Harris capturent avec succès des phénomènes physiques complexes, notamment la reconnexion magnétique, la formation d'îlots magnétiques, de cordes de flux et les motifs de courant quadrupolaires caractéristiques associés à la dynamique Hall–MHD. Le schéma résout ces caractéristiques avec précision même sur des maillages relativement grossiers.

Signification
L'article affirme que ce travail fournit un cadre mathématiquement rigoureux et computationnellement robuste pour la simulation des systèmes Hall–MHD. En combinant la régularisation de Voigt avec une formulation aux éléments finis mixtes préservant la divergence, les auteurs surmontent les défis duaux de la non-linéarité du terme de Hall et de l'exigence stricte de divergence nulle. La méthode offre une approximation physiquement cohérente qui est stable, convergente et capable de gérer des géométries complexes, comblant ainsi une lacune laissée par les méthodes précédentes qui manquaient d'une analyse d'erreur rigoureuse ou d'une préservation exacte de la divergence. Les résultats suggèrent que cette approche est un outil viable pour l'étude des phénomènes de plasma astrophysiques et de laboratoire où les effets Hall sont dominants.