Primitive variable regularization to derive novel Hyperbolic Shallow Water Moment Equations

Cet article propose de nouvelles équations des moments pour les eaux peu profondes en régularisant les variables primitives plutôt que les variables convectives, garantissant ainsi l'hyperbolicité, la conservation de l'équation de la quantité de mouvement et une plus grande précision par rapport aux modèles existants.

L'essentiel

Imaginez que vous essayez de prédire comment une vague va se déplacer dans une rivière ou lors d'une inondation. Pour les scientifiques, c'est un peu comme essayer de comprendre le mouvement de millions de gouttes d'eau individuellement. C'est trop compliqué et prend trop de temps pour faire des prédictions rapides.

Pour simplifier, on utilise souvent une "moyenne". C'est comme dire : "Toute l'eau dans cette rivière va à 5 km/h". C'est rapide à calculer, mais ce n'est pas très précis, car en réalité, l'eau près du fond va plus lentement que celle en surface (à cause du frottement).

Le problème : La recette qui tourne mal

Les chercheurs ont essayé de faire mieux en ajoutant des "détails" à cette moyenne. Ils ont créé des modèles appelés équations des moments. Au lieu de juste une vitesse moyenne, ils ajoutent des coefficients (comme des ingrédients supplémentaires dans une recette) pour décrire comment la vitesse change du fond au sommet.

Le problème, c'est que ces modèles existants avaient deux gros défauts :

1. Ils devenaient fous (instables) : Parfois, les mathématiques derrière ces modèles s'effondraient, donnant des résultats impossibles (comme une vitesse imaginaire). C'est comme si votre GPS vous disait de conduire à travers un immeuble.
2. Ils ne savaient pas s'arrêter : Quand la rivière est calme (un état stable), ces modèles ne pouvaient pas décrire mathématiquement cet état de repos de manière simple.

La solution : Changer de point de vue

Dans cet article, l'auteur, Julian Koellermeier, propose une astuce géniale.

Imaginez que vous essayez de réparer une voiture qui fait des bruits bizarres.

• Les anciennes méthodes regardaient le moteur sous le capot (les variables "convectives"). Elles ont essayé de simplifier les pièces pour que ça marche, mais en simplifiant trop, ils ont perdu la précision ou la stabilité.
• La nouvelle méthode de l'auteur consiste à regarder la voiture depuis l'extérieur, en regardant la route et le conducteur (les variables "primitives").

En changeant de point de vue (en passant des variables "convectives" aux variables "primitives"), l'auteur a pu appliquer une régularisation (une sorte de "stabilisateur" mathématique) qui fonctionne beaucoup mieux.

Les deux nouveaux modèles magiques

L'auteur a créé deux nouveaux modèles, mais le vrai champion est le dernier, le PMHSWME.

1. Le modèle PHSWME : C'est une première version qui stabilise le système. C'est comme mettre un stabilisateur sur une caméra : l'image ne tremble plus. Mais, pour le faire, il a dû modifier légèrement la façon dont on calcule la vitesse moyenne, ce qui a rendu le résultat un peu moins fidèle à la réalité physique.
2. Le modèle PMHSWME (Le gagnant) : C'est la version ultime. L'auteur a eu une idée brillante : "Gardez la recette de la vitesse moyenne exactement telle qu'elle est, ne touchez à rien, et ne régularisez (ne stabilisez) que les détails complexes."

C'est comme si vous aviez une recette de gâteau parfaite pour la base (la vitesse moyenne), mais que la crème fouettée (les détails complexes) était instable. Au lieu de changer la base, vous stabilisez juste la crème. Résultat : vous avez un gâteau stable et délicieux.

Pourquoi c'est important ?

Grâce à cette nouvelle approche, les chercheurs peuvent maintenant :

• Être sûrs que le modèle ne va pas planter (il est "globalement hyperbolique", ce qui signifie mathématiquement stable).
• Calculer des états de calme (quand la rivière ne bouge plus) sans se casser la tête.
• Garder une grande précision, car ils n'ont pas sacrifié la partie la plus importante de l'équation (la conservation de la quantité de mouvement).

En résumé

C'est comme si on avait un véhicule de course (le modèle mathématique) qui avait tendance à dérailler sur les virages. Les anciens mécaniciens ont essayé de retirer des pièces pour le rendre plus stable, mais il était moins rapide.

L'auteur de cet article a dit : "Non, on ne retire rien ! On change juste la façon dont on regarde la route et on ajoute un système de suspension intelligent."

Le résultat ? Un véhicule qui va vite, qui ne dérape pas, et qui arrive exactement là où on veut, même dans les situations les plus extrêmes. C'est une avancée majeure pour mieux prévoir les inondations et gérer les cours d'eau.

Résumé technique

1. Problématique

Les écoulements à surface libre (tsunamis, inondations, estuaires) sont souvent modélisés par les équations de Saint-Venant (SWE), qui supposent un profil de vitesse constant dans la verticale. Cependant, cette hypothèse est insuffisante pour capturer des effets physiques complexes comme le frottement au fond, la végétation ou le transport de sédiments.

Pour pallier cela, les Équations des Moments des Eaux Peu Profondes (SWME) ont été développées. Elles utilisent un développement de la vitesse verticale sur une base de polynômes de Legendre, générant des équations d'évolution pour les coefficients de ce développement (les moments). Bien que plus précises, les SWME originales souffrent de deux défauts majeurs :

1. Manque d'hyperbolicité globale : Pour des ordres de troncature $N \ge 2$, le système perd son caractère hyperbolique (les valeurs propres deviennent complexes) en dehors de l'équilibre, entraînant des instabilités numériques.
2. Absence d'états stationnaires analytiques : Contrairement aux SWE classiques, il est impossible de dériver des états stationnaires (solutions de type saut hydraulique) à partir des conditions de Rankine-Hugoniot pour les SWME régularisées existantes, ce qui empêche le développement de schémas numériques bien-équilibrés.

Des régularisations antérieures (HSWME, SWLME) ont tenté de résoudre l'un ou l'autre problème, mais aucune n'a réussi à satisfaire simultanément l'hyperbolicité globale, la conservation de l'équation de quantité de mouvement, la présence d'états stationnaires analytiques et la précision physique.

2. Méthodologie

L'auteur propose une nouvelle approche de régularisation basée sur un changement de variables fondamental. Au lieu de régulariser le système dans les variables conservatives (ou convectives) $U_c = (h, hu_m, h\alpha_1, \dots)$, comme cela a été fait précédemment, la méthode proposée opère dans les variables primitives $U_p = (h, u_m, \alpha_1, \dots)$.

La méthodologie se déroule en plusieurs étapes clés :

1. Transformation : Le système SWME est transformé des variables convectives vers les variables primitives en utilisant une matrice de Jacobienne explicite.
2. Régularisation Primitive : Une régularisation hyperbolique est appliquée directement sur la matrice du système dans les variables primitives. Cela consiste à évaluer la matrice autour de profils de vitesse linéaires (en posant $\alpha_i = 0$ pour $i \ge 2$ dans les termes de couplage non linéaires des équations de moments supérieurs), tout en préservant la structure des équations de masse et de quantité de mouvement.
3. Retour aux variables convectives : La matrice régularisée dans les variables primitives est transformée de nouveau en variables convectives via une transformation de similarité. Cela permet de prouver que les valeurs propres (et donc l'hyperbolicité) sont préservées.
4. Développement de deux nouveaux modèles :
   • PHSWME : Régularisation primitive complète (modifie l'équation de quantité de mouvement).
   • PMHSWME (Primitive Moment Hyperbolic SWME) : Une variante hybride qui applique la régularisation primitive uniquement aux équations des moments d'ordre supérieur ($i \ge 2$), tout en conservant exactement l'équation de quantité de mouvement originale des SWME.

3. Contributions Clés

L'article introduit et analyse théoriquement deux nouveaux modèles hiérarchiques :

• PHSWME (Primitive Hyperbolic Shallow Water Moment Equations) :

  • Démonstration analytique de l'hyperbolicité globale pour tout ordre $N$. Les valeurs propres sont réelles et distinctes.
  • Dérivation d'états stationnaires analytiques fermés, généralisant ceux des SWE classiques.
  • Cependant, ce modèle simplifie excessivement la physique en modifiant l'équation de quantité de mouvement, ce qui peut réduire la précision.

• PMHSWME (Primitive Moment Hyperbolic Shallow Water Moment Equations) :

  • C'est la contribution principale. Ce modèle combine la régularisation primitive pour les moments d'ordre supérieur avec la préservation stricte de l'équation de quantité de mouvement.
  • Il conserve l'hyperbolicité globale et admet des états stationnaires analytiques (dépendant de tous les moments, contrairement au PHSWME).
  • Il préserve la non-linéarité physique des équations de transport des moments, assurant une meilleure fidélité physique.

Le papier fournit des preuves mathématiques rigoureuses (polynômes caractéristiques, valeurs propres liées aux dérivées de polynômes de Legendre) pour l'hyperbolicité et les états stationnaires de ces nouveaux systèmes.

4. Résultats Numériques

Un test de rupture de barrage (dam-break) a été utilisé pour comparer les nouveaux modèles (PHSWME, PMHSWME) avec les modèles existants (SWME, HSWME, SWLME, MHSWME).

• Précision : Le modèle PMHSWME démontre la meilleure précision globale pour toutes les variables (hauteur d'eau $h$, vitesse moyenne $u_m$, et coefficients de moments $\alpha_i$). Il s'écarte très peu des solutions de référence (SWME non régularisé ou solveur de référence 2D).
• Comparaison avec d'autres modèles :
  • Le SWLME présente les erreurs les plus élevées, confirmant que sa régularisation trop forte (linéarisation des moments) dégrade la physique.
  • Le HSWME et le MHSWME (régularisation dans les variables convectives) montrent des erreurs significatives sur la hauteur d'eau et la vitesse, soulignant l'importance de la régularisation dans les variables primitives.
  • Le PHSWME est précis pour les moments mais moins précis pour la hauteur et la vitesse moyenne, ce qui confirme que modifier l'équation de quantité de mouvement est préjudiciable.
• Conclusion sur la régularisation : Le résultat clé est que la régularisation doit être effectuée dans les variables primitives et que l'équation de quantité de mouvement doit rester inchangée pour maintenir une haute précision tout en garantissant la stabilité mathématique.

5. Signification et Impact

Ce travail résout un problème ouvert majeur dans la modélisation des écoulements à surface libre de haute fidélité : la conciliation entre la stabilité mathématique (hyperbolicité), la capacité à construire des schémas numériques robustes (états stationnaires analytiques) et la précision physique.

• Avancée théorique : La démonstration que la régularisation dans les variables primitives est supérieure à celle dans les variables convectives pour ce type de système ouvre une nouvelle voie pour la dérivation de modèles hyperboliques dans d'autres domaines (dynamique des gaz raréfiés, etc.).
• Application pratique : Le modèle PMHSWME proposé est le premier à offrir un compromis optimal. Il permet des simulations stables et précises de phénomènes complexes (inondations rapides, tsunamis) tout en étant compatible avec des méthodes numériques avancées (schémas bien-équilibrés).
• Perspectives : L'article ouvre la voie à l'extension de ces modèles en dimensions supérieures (2D/3D), à l'inclusion de termes de frottement réalistes et à l'analyse des propriétés d'entropie.

En résumé, l'article propose une méthodologie de régularisation novatrice qui permet de dériver des modèles d'ondes de gravité à surface libre à la fois mathématiquement robustes et physiquement précis, comblant ainsi le fossé entre les modèles simplifiés (SWE) et les modèles complets coûteux (Navier-Stokes).