🔬 physics

A Hyperbolic Moment Based Shallow Water Model for Coupled Bedload Suspended Load Morphodynamics with Variable Density

Cet article présente le modèle HSWEMED, une extension hyperbolique des équations de Saint-Venant par moments qui intègre la densité variable, l'érosion, le dépôt et le transport de sédiments en suspension pour simuler avec précision la morphodynamique couplée.

Auteurs originaux : Afroja Parvin, Giovanni Samaey, Julian Koellermeier

Publié 2026-04-21

📖 4 min de lecture☕ Lecture pause café

CC BY 4.0

Auteurs originaux : Afroja Parvin, Giovanni Samaey, Julian Koellermeier

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez une rivière qui coule. Ce n'est pas juste de l'eau qui glisse sur un lit de pierre. C'est un ballet complexe où l'eau transporte des cailloux (le lit) et de la boue fine (la suspension). Parfois, l'eau arrache des cailloux au fond (érosion), parfois elle les dépose (dépôt). Le lit de la rivière change de forme en temps réel, comme un paysage vivant qui respire.

Jusqu'à présent, les mathématiciens et ingénieurs utilisaient des modèles simplifiés pour prédire ces changements. C'était un peu comme essayer de décrire le vent dans une tempête en disant simplement : « Il souffle fort ». C'est utile, mais ça ne vous dit pas si le vent tourne en spirale ou s'il y a des tourbillons dangereux près du sol.

Voici ce que cette nouvelle recherche propose, expliqué simplement :

1. Le problème : La vision "plate" des anciens modèles

Les modèles classiques (appelés Shallow Water) traitent l'eau comme une seule couche plate. Ils calculent la vitesse moyenne de l'eau.

L'analogie : Imaginez que vous regardez un tapis roulant. Vous voyez la vitesse moyenne du tapis. Mais si vous voulez savoir si un objet posé dessus va glisser, vous devez savoir si le tapis bouge plus vite en haut ou en bas.
Le souci : Près du fond de la rivière, l'eau frotte contre les roches et ralentit. Les vieux modèles ne voient pas ce ralentissement. Ils pensent que l'eau va vite partout. Résultat : ils prédisent souvent qu'il y a trop d'érosion (trop de roches arrachées) parce qu'ils surestiment la force de l'eau au contact du sol.

2. La solution : Le modèle "Moment" (HSWEMED)

Les auteurs ont créé un nouveau modèle, le HSWEMED. Au lieu de voir l'eau comme une couche plate, ils la découpent en plusieurs couches imaginaires pour voir comment la vitesse change de haut en bas.

L'analogie : C'est comme passer d'une photo en noir et blanc à une photo en 3D avec des couleurs. Au lieu de dire "l'eau va à 2 m/s", le modèle dit : "En surface, elle va à 3 m/s, au milieu à 2 m/s, et tout contre le fond, elle ne va qu'à 0,5 m/s".
Pourquoi c'est génial ? Parce que c'est cette vitesse tout contre le fond qui décide si un caillou va bouger ou non. En ayant cette vision précise, le modèle calcule beaucoup mieux quand la rivière va creuser son lit ou le combler.

3. La touche de magie : La densité variable et la boue

Ce modèle fait aussi quelque chose de très intelligent : il se rend compte que l'eau n'est pas toujours pure. Quand il y a beaucoup de boue en suspension, l'eau devient plus lourde (plus dense).

L'analogie : C'est comme comparer un verre d'eau claire à un verre de chocolat chaud épais. Le chocolat est plus lourd et se comporte différemment.
Le nouveau modèle intègre cette "lourdeur" de la boue dans ses calculs. Il sait que si l'eau est chargée de sable, elle pousse différemment sur le fond, ce qui change la façon dont les vagues se propagent.

4. La stabilité mathématique : Le "Frein de sécurité"

Un défi majeur en mathématiques est que ces modèles complexes peuvent parfois devenir "fous" et donner des résultats impossibles (comme des vitesses imaginaires !).

L'analogie : C'est comme conduire une voiture de course très puissante. Si vous appuyez trop fort sur l'accélérateur sans freins, vous perdez le contrôle.
Les auteurs ont ajouté un "frein mathématique" (une régularisation hyperbolique) qui garantit que le modèle reste stable et réaliste, même dans les situations les plus extrêmes, comme une crue soudaine (rupture de barrage).

5. Les résultats : Plus proche de la réalité

Les chercheurs ont testé leur modèle avec des expériences de rupture de barrage (l'eau qui dévale une digue soudainement).

Le verdict : Le nouveau modèle (HSWEMED) colle beaucoup mieux aux expériences réelles faites en laboratoire que les anciens modèles.
Pourquoi ? Parce qu'il voit mieux la vitesse au fond de l'eau et comprend mieux comment la boue et l'eau interagissent. Il prédit mieux où l'eau va creuser et où elle va déposer du sable.

En résumé

Cette recherche est comme donner des lunettes à haute définition aux ingénieurs qui étudient les rivières.

Avant : Ils voyaient la rivière de loin, de manière floue, et faisaient des erreurs sur l'érosion.
Maintenant : Avec le modèle HSWEMED, ils voient les détails de la vitesse près du fond et la densité de la boue. Cela permet de mieux prévoir les inondations, l'évolution des lits des rivières et la sécurité des barrages.

C'est un pas de géant vers une compréhension plus précise et plus sûre de la façon dont l'eau sculpte notre planète.

1. Problématique et Contexte

La modélisation du transport sédimentaire en eau peu profonde est cruciale pour la gestion des risques d'inondation, l'ingénierie côtière et fluviale. Les sédiments sont transportés sous deux formes principales :

Charge de fond (Bedload) : particules grossières roulant ou sautant près du lit.
Charge en suspension (Suspended load) : particules fines maintenues en suspension par la turbulence.

Les modèles classiques (Shallow Water Exner) reposent souvent sur des hypothèses simplificatrices limitantes :

Ils utilisent des profils de vitesse verticaux idéalisés (logarithmiques) et des lois de frottement empiriques (Manning, Chézy) calibrées pour des écoulements stationnaires.
Ils négligent la structure verticale réelle de la vitesse, ce qui rend l'estimation des contraintes de cisaillement au fond (cruciales pour l'érosion) imprécise dans des régimes non stationnaires.
Ils traitent souvent la densité du mélange eau-sédiment comme constante, ignorant les effets de flottabilité.
Les modèles à moments (SWM) récents, bien qu'ils reconstruisent le profil de vitesse vertical, perdent l'hyperbolicité (stabilité numérique) pour les ordres supérieurs et négligent souvent les processus d'érosion/dépôt et la charge en suspension.

Objectif : Développer un modèle mathématiquement bien posé, hyperbolique et capable de simuler couplément la charge de fond et la charge en suspension, en tenant compte de la densité variable et des processus d'érosion/dépôt.

2. Méthodologie

Les auteurs proposent le modèle HSWEMED (Hyperbolic Shallow Water Exner Moment model with Erosion and Deposition). La démarche se décompose en plusieurs étapes clés :

A. Dérivation à partir des équations de Navier-Stokes

Le modèle part des équations de Navier-Stokes incompressibles pour un mélange eau-sédiment, où la densité $\rho$ varie selon la concentration volumique de sédiments $c$ .

Approximation hydrostatique : La pression est hydrostatique, mais la densité du mélange est prise en compte dans la pression et les termes de flottabilité.
Approche par moments : La vitesse horizontale $u(z)$ est développée en série de polynômes de Legendre sur la coordonnée verticale normalisée $\zeta \in [0,1]$ . Cela permet de capturer la structure verticale de la vitesse sans discrétiser le domaine en plusieurs couches (comme dans les modèles multi-couches).
Équations couplées : Le système intègre :
1. Les équations de Saint-Venant (masse et quantité de mouvement moyennées).
2. Des équations d'évolution pour les moments de vitesse (coefficients du développement polynomial).
3. Une équation de transport pour la concentration moyenne de sédiments en suspension ( $c_m$ ).
4. L'équation d'Exner pour l'évolution du lit ( $h_b$ ).

B. Gestion de l'érosion et du dépôt

Contrairement aux modèles précédents, HSWEMED inclut explicitement les termes sources d'érosion ( $E$ ) et de dépôt ( $D$ ) dans les équations de masse et de quantité de mouvement.

Reconstruction verticale : Bien que la variable transportée soit la concentration moyenne $c_m$ , le modèle reconstruit un profil vertical quadratique de la concentration pour estimer la concentration près du fond ( $c_b$ ), nécessaire pour calculer les flux d'érosion/dépôt.
Densité variable : La densité du mélange $\rho(c_m)$ couple la charge en suspension à la dynamique hydrodynamique via la pression hydrostatique et les termes de flottabilité.

C. Hyperbolicité et Régularisation

Un défi majeur des modèles à moments est la perte d'hyperbolicité (apparition de valeurs propres complexes) pour les ordres supérieurs ( $N \ge 2$ ), entraînant une instabilité numérique.

Régularisation : Les auteurs appliquent une technique de régularisation hyperbolique (inspirée de [36]) pour garantir que la matrice Jacobienne du système possède des valeurs propres réelles.
HSWEMED : Le modèle régularisé résultant est nommé HSWEMED. L'analyse spectrale montre que le système reste hyperbolique pour tout ordre $N$ , assurant la stabilité des schémas numériques.

D. Analyse énergétique

Pour les modèles d'ordre inférieur (N=0 et N=1), les auteurs dérivent des lois de bilan d'énergie mécanique dissipative, prouvant que le système respecte les principes thermodynamiques (dissipation d'énergie par frottement et turbulence).

3. Contributions Clés

Cadre couplé unifié : Première extension d'un modèle à moments hyperbolique pour simuler simultanément la charge de fond et la charge en suspension avec échange de masse et de quantité de mouvement entre le lit et la colonne d'eau.
Densité variable : Intégration rigoureuse de la densité du mélange eau-sédiment dans les équations de quantité de mouvement, permettant de capturer les effets de flottabilité.
Reconstruction du profil de concentration : Utilisation d'un profil quadratique pour relier la concentration moyenne à la concentration près du fond, améliorant la précision des lois de fermeture d'érosion/dépôt sans alourdir le coût computationnel.
Preuve d'hyperbolicité : Démonstration analytique que le système régularisé HSWEMED conserve l'hyperbolicité, garantissant la stabilité numérique.
Amélioration de la vitesse au fond : Le modèle calcule dynamiquement la vitesse au fond ( $u_b$ ) à partir des moments, évitant l'utilisation de profils de vitesse moyens ou logarithmiques empiriques pour le calcul du frottement et du transport solide.

4. Résultats Numériques

Les auteurs ont validé le modèle HSWEMED sur plusieurs cas tests, comparés aux modèles classiques SWEED (Shallow Water Exner avec érosion/dépôt, profil de vitesse uniforme) et HSWEM (modèle à moments sans érosion/dépôt ni suspension).

Test académique de rupture de barrage :
- HSWEMED prédit une érosion intermédiaire entre SWEED (qui surestime l'érosion en utilisant la vitesse moyenne) et HSWEM (qui sous-estime l'érosion car il ignore le dépôt).
- Le modèle capture correctement la forme parabolique du profil de vitesse près du fond.
Expériences de laboratoire (Rupture de barrage sur lit érodible) :
- Trois configurations ont été testées avec deux matériaux : billes de PVC (faible densité) et sable grossier (haute densité).
- PVC : Les particules légères restent plus longtemps en suspension. HSWEMED reproduit mieux l'évolution du lit et la position du saut hydraulique que SWEED et HSWEM, grâce à une meilleure estimation de la vitesse au fond et des processus d'érosion/dépôt.
- Sable : Les différences entre les modèles sont moins marquées car la densité élevée et la vitesse de sédimentation rapide limitent la charge en suspension. Cependant, HSWEMED reste plus précis, surtout lorsque le coefficient de frottement est augmenté.
- Concentration en suspension : HSWEMED prédit des concentrations en suspension plus réalistes que SWEED, qui tend à surestimer l'érosion et donc la concentration.

5. Signification et Conclusion

Le modèle HSWEMED représente une avancée significative dans la modélisation morphodynamique :

Précision physique : Il surpasse les modèles classiques en résolvant la structure verticale de la vitesse, ce qui est critique pour le calcul précis des contraintes de cisaillement au fond et des taux d'érosion.
Robustesse mathématique : La régularisation assure que le modèle reste stable numériquement même pour des ordres de moments élevés, évitant les instabilités des modèles à moments non régularisés.
Efficacité : Il offre un compromis optimal entre la précision des modèles 3D complets (coûteux) et la simplicité des modèles 1D classiques (imprécis), permettant des simulations morphodynamiques réalistes à un coût computationnel modéré.

Ce travail ouvre la voie à des simulations plus fiables des écoulements complexes en présence de sédiments, notamment dans des régimes transitoires et non équilibrés, et pose les bases pour des développements futurs incluant des schémas numériques d'ordre supérieur et des analyses de stabilité.