High-Order Schemes for Hyperbolic Conservation Laws Using Young Measures

Résumé Technique : Schémas d'ordre élevé pour les lois de conservation hyperboliques utilisant les mesures de Young

Énoncé du problème
L'article traite de la résolution numérique de lois de conservation hyperboliques aléatoires, qui apparaissent dans des applications d'ingénierie, atmosphériques et géophysiques. Ces systèmes sont caractérisés par le développement de motifs d'ondes complexes, incluant des chocs et des ondes de raréfaction, même à partir de données initiales lisses. Les méthodes numériques existantes font face à des défis importants :

• Les méthodes de Monte Carlo sont robustes mais coûteuses en termes de calcul.
• La collocation stochastique permet de réutiliser les solveurs déterministes mais souffre de limitations d'approximation spectrale.
• Les méthodes intrusives (par exemple, la méthode de Galerkin stochastique avec chaos polynomial généralisé) entraînent souvent des systèmes couplés qui perdent l'hyperbolicité et nécessitent une réimplémentation substantielle des solveurs.
• Les méthodes spectrales sont sujettes au phénomène de Gibbs à proximité des discontinuités, nécessitant des techniques supplémentaires de filtrage ou de viscosité.

De plus, les méthodes numériques standards peinent souvent face à des problèmes où les solutions sont véritablement à valeurs de mesures (non atomiques) ou impliquent des flux discontinus, ce qui conduit fréquemment à l'effondrement des distributions vers des états atomiques.

Méthodologie
Les auteurs étendent un cadre numérique de premier ordre précédemment proposé, basé sur les mesures de Young et la programmation linéaire, à des ordres de précision plus élevés. La méthodologie centrale comprend :

1. Formulation par mesure : Au lieu de chercher des solutions faibles classiques, l'approche traite la solution comme une mesure de Young paramétrée $\nu_{t,x,\xi}$. Les équations directrices sont transformées en une équation aux dérivées partielles linéaire concernant la mesure de Young.
2. Clôture par programmation linéaire : Pour fermer le système, la méthode minimise une fonctionnelle d'entropie convexe sous des contraintes de moments. Cela transforme l'étape de reconstruction en un problème de programmation linéaire (LP) de dimension finie. La mesure de Young discrète est approximée par des coefficients $\mu_{i,j,\ell}$ qui minimisent l'entropie tout en satisfaisant la contrainte du premier moment (le vecteur de solution discrète).
3. Extensions d'ordre élevé :
   • Schéma du second ordre : Obtenu via une reconstruction linéaire par morceaux des valeurs ponctuelles unilatérales des inconnues. Pour prévenir les oscillations parasites, un limiteur minmod généralisé est appliqué aux variables caractéristiques locales via une décomposition caractéristique locale (LCD). Le flux spatial est calculé à l'aide d'un flux de Lax–Friedrichs local (LLF).
   • Schéma du cinquième ordre : Développé dans le cadre des méthodes de différences finies A-WENO (Alternative Weighted Essentially Non-Oscillatory). Cela étend le schéma de volumes finis (FV) de bas ordre à un schéma de différences finies (FD) d'ordre élevé. Le flux numérique inclut des termes de correction d'ordre supérieur calculés efficacement à partir des flux FV plutôt que des valeurs ponctuelles. Comme le schéma du second ordre, il utilise la LCD et l'interpolation WENO-Z sur les variables caractéristiques pour maintenir des propriétés non oscillatoires.
4. Intégration temporelle : Les systèmes de systèmes d'équations différentielles ordinaires (ODE) semi-discrets résultants sont résolus à l'aide d'une méthode Runge-Kutta de troisième ordre à stabilité forte préservée (SSP RK3) à trois étapes.

Contributions clés

• Généralisation d'ordre élevé : L'article étend avec succès le schéma de mesure de Young de premier ordre existant aux ordres de précision deux et cinq.
• Formulation de flux efficace : Les auteurs proposent des formulations de flux LLF computationnellement efficaces dans le contexte d'ordre élevé. Spécifiquement, ils démontrent que le calcul direct des flux à partir des moments reconstruits (plutôt que de résoudre le problème LP deux fois pour les états d'interface) réduit considérablement le temps CPU sans sacrifier la précision.
• Gestion des solutions non atomiques : En contrôlant le paramètre $\lambda_F$ dans les contraintes de la programmation linéaire (spécifiquement en fixant $\lambda_F < 1$), le schéma peut imposer des mesures de Young discrètes non atomiques, permettant de capturer avec précision des solutions qui ne peuvent être représentées par des états atomiques.
• Polyvalence : Le cadre est montré capable de gérer à la fois des systèmes paramétrés aléatoires et des systèmes déterministes présentant des caractéristiques difficiles, telles que des flux discontinus et des systèmes faiblement hyperboliques (dynamique des gaz sans pression).

Résultats numériques
Les schémas proposés ont été validés par une série d'expériences numériques sur des problèmes tant aléatoires que déterministes :

• Lois hyperboliques aléatoires : Des tests sur l'équation de Burgers aléatoire et les équations d'Euler isentropiques ont démontré que les schémas de second et cinquième ordre améliorent significativement la résolution des ondes de choc par rapport au schéma de premier ordre. L'analyse d'erreur ($L^1, L^2, L^\infty$) a confirmé les taux de convergence attendus (environ 2ème et 5ème ordre) et a montré que les schémas d'ordre élevé atteignent une précision et une efficacité supérieures pour un coût de calcul donné.
• Flux discontinu : Des expériences avec des lois de conservation scalaires présentant des flux discontinus ont illustré comment le choix de la fonction d'entropie dans l'objectif de la programmation linéaire agit comme un critère de sélection, influençant la solution reconstruite et le support de la mesure de Young.
• Dynamique des gaz sans pression : Les schémas ont capturé avec succès les chocs de type delta et les états de vide dans le système de dynamique des gaz sans pression, qui est faiblement hyperbolique ; les schémas d'ordre élevé produisant des pics plus nets que leur homologue de premier ordre, même sur des mailles plus fines.
• Support non atomique : Dans l'exemple de l'équation de Burgers avec un support non atomique, les schémas ont réussi à capturer des mesures de Young distribuées lorsque $\lambda_F < 1$, évitant l'effondrement vers des solutions atomiques typique des méthodes standards.

Signification et affirmations
L'article affirme que les extensions d'ordre élevé proposées améliorent considérablement la résolution des discontinuités et des structures d'ondes complexes tout en préservant les propriétés de conservation de la structure de l'approche originale par mesure de Young. Les résultats suggèrent que ces schémas offrent une alternative robuste et efficace aux méthodes existantes pour les lois de conservation hyperboliques aléatoires, particulièrement dans les scénarios impliquant des chocs, des raréfactions et des structures de solutions non atomiques. Les auteurs notent que la méthode s'étend naturellement aux systèmes d'équations et élimine le besoin de filtrage spectral ou de viscosité artificielle souvent requis par d'autres méthodes d'ordre élevé.

Le travail conclut en identifiant des directions futures, spécifiquement l'extension de ces approches aux problèmes hyperboliques aléatoires multidimensionnels et l'amélioration de l'efficacité computationnelle par des stratégies de raffinement adaptatif.