Adaptive Artificial Anti-Diffusion Methods for Hyperbolic Systems of Conservation Laws

Résumé technique : Méthodes adaptatives d'anti-diffusion artificielle pour les systèmes hyperboliques de lois de conservation

Énoncé du problème
L'article aborde le défi fondamental dans la résolution numérique des systèmes hyperboliques de lois de conservation (tels que les équations d'Euler de la dynamique des gaz) : concilier la robustesse face aux oscillations non physiques avec une haute résolution des structures physiques. Bien que les méthodes de capture de choc à haute résolution utilisant des limiteurs non linéaires suppriment efficacement les oscillations parasites, elles introduisent souvent une dissipation numérique excessive. Cette dissipation entraîne un lissage excessif des discontinuités de contact et des structures lisses à petite échelle, en particulier sur des maillages grossiers. À l'inverse, les méthodes classiques de viscosité artificielle (AV), bien que robustes, peinent à maintenir une précision d'ordre élevé dans les régions lisses sans dégrader la résolution. Les auteurs notent que même les méthodes avancées de viscosité artificielle adaptative (AAV) n'atteignent pas les niveaux de résolution des schémas modernes de capture de choc.

Méthodologie
Les auteurs proposent une nouvelle classe de méthodes d'Anti-Diffusion Artificielle Adaptative (AAAD) conçues pour améliorer la résolution des schémas stabilisés existants sans compromettre leurs propriétés non oscillatoires. La méthodologie centrale comprend les composants suivants :

1. Correction d'Anti-Diffusion (AD) : Au lieu d'ajouter de la viscosité, la méthode soustrait un terme anti-diffusif du flux numérique. Ce terme est construit pour agir uniquement dans les champs caractéristiques linéairement dégénérés (par exemple, l'onde de contact dans les équations d'Euler, associée à la valeur propre $\lambda = u$). En restreignant la correction AD à ces champs, la méthode affine les ondes de contact tout en évitant l'introduction d'oscillations indésirables dans les champs véritablement non linéaires (chocs et détentes).
2. Sélection adaptative du coefficient : L'amplitude de la correction AD est déterminée de manière adaptative sur la base d'un Indicateur de Lissage (SI). Les auteurs emploient un SI modifié basé sur minmod (de [7]) pour distinguer les régions « rugueuses » (discontinuités) des régions « lisses ».
   • Ondes de contact : Identifiées par une discontinuité de densité ($\rho$) mais une continuité de pression ($p$). Dans ces régions de « contact rugueux », le coefficient AD est fixé à l'ordre $O(\Delta x)$ pour maximiser la résolution.
   • Autres régions rugueuses : Dans les régions comportant des chocs ou d'autres discontinuités où l'AD n'est pas appliquée (ou appliquée minimalement), le coefficient est fixé à $O(\Delta x^2)$.
   • Régions lisses : Pour préserver l'ordre formel de précision du schéma sous-jacent, le coefficient AD est fixé à une valeur très faible, $O(\Delta x^r)$ (où $r$ est l'ordre formel du schéma), assurant que le terme anti-diffusif s'annule asymptotiquement.
3. Cadres d'implémentation : La stratégie AAAD est implémentée dans deux cadres numériques distincts :
   • Ordre deux : Basé sur les schémas aux volumes finis Central-Upwind (CU).
   • Ordre cinq : Basé sur le cadre aux différences finies Alternatif Weighted Essentially Non-Oscillatory (A-WENO).
   • Des extensions 1-D et 2-D sont fournies, l'implémentation 2-D appliquant la correction AD dimension par dimension.

Contributions clés

• Mécanisme novateur : L'introduction d'un terme anti-diffusif adaptatif ciblant spécifiquement les champs linéairement dégénérés, offrant un mécanisme pour affiner les discontinuités de contact sans déstabiliser la solution.
• Intégration avec des schémas d'ordre élevé : L'extension réussie de cette stratégie aux schémas A-WENO d'ordre cinq, démontrant qu'une précision d'ordre élevé peut être maintenue dans les régions lisses tout en améliorant la résolution dans les zones de contact.
• Procédure de réglage robuste : L'article décrit une procédure pratique pour régler la constante adaptative $C$ sur des maillages grossiers, qui est ensuite transférable à des maillages plus fins, garantissant que la méthode reste robuste à travers différentes résolutions de grille.

Résultats
Les méthodes AAAD proposées (désignées AAAD2 et AAAD5) ont été testées sur une série complète de benchmarks 1-D et 2-D pour les équations d'Euler, incluant :

• Tests de précision : Les deux schémas ont atteint leurs ordres formels de précision conçus (2e et 5e) dans les tests d'écoulement lisse.
• Interactions Choc-Entropie et Choc-Densité : Les schémas AAAD ont produit des résultats nettement plus précis que leurs homologues non adaptatifs CU et A-WENO, résolvant les perturbations de densité à haute fréquence avec beaucoup moins de dissipation.
• Discontinuités de contact : Dans les problèmes impliquant des interactions choc-bulle, les problèmes de Lax et les ondes de choc, les méthodes AAAD ont fourni une résolution plus nette des ondes de contact.
• Écoulements complexes 2-D : Dans les problèmes de Riemann 2-D (Configurations 3, 6 et 12), les problèmes d'explosion/implosion et les tests d'instabilité (Kelvin-Helmholtz et Rayleigh-Taylor), les schémas AAAD ont capturé des structures à petite échelle plus riches, des rues de tourbillons et des couches de mélange.
• Comparaison des performances : Notamment, dans plusieurs exemples 2-D (par exemple, Configurations de Riemann 6 et 12, Implosion), le schéma AAAD2 d'ordre deux a surpassé le schéma A-WENO d'ordre cinq en termes de résolution, soulignant l'efficacité de la stratégie anti-diffusive pour réduire la dissipation numérique.

Signification et affirmations
L'article affirme que les méthodes AAAD proposées offrent une solution robuste et efficace au compromis entre stabilité et résolution dans les lois de conservation hyperboliques. En ajoutant de l'anti-diffusion de manière adaptative uniquement là où c'est nécessaire (champs linéairement dégénérés) et uniquement là où la solution est « rugueuse », la méthode améliore considérablement la résolution des ondes de contact et des structures fines sans détruire le caractère non oscillatoire garanti par les limiteurs non linéaires sous-jacents. Les auteurs affirment que ces méthodes sont « substantiellement moins dissipatives » que leurs homologues non adaptatifs, ce qui les rend particulièrement avantageuses pour les problèmes où la résolution précise des discontinuités de contact et des caractéristiques à petite échelle est critique.