Simulation numérique d'un essai de surverse

Nous considérons une digue en terre soumise à un écoulement de surverse. La modélisation physique et numérique d’un tel phénomène est en général un compromis entre pertinence physique, précision numérique et temps de calcul. Nous nous sommes attachés à l’analyse de l’écoulement à surface libre caractéristique d’une surverse, éludant momentanément les phénomènes d’érosion et de transport sédimentaire qui seront inclus ultérieurement par couplage. L’hydrodynamique complexe associée à la surverse est traitée grâce à un modèle bi-fluide Eulérien à faible nombre de Mach. Ce modèle bi-fluide, non visqueux et faiblement compressible est physiquement plus pertinent que des modèles de type Saint-Venant généralement utilisés. En effet, on observe dans l’écoulement des jets et des recirculations plus facilement appréhendés par un modèle diphasique. Ce modèle est évidemment moins riche qu’un modèle de Navier-Stokes avec ou sans turbulence, mais sa formulation permet des simulations tridimensionnelles beaucoup plus économes en temps de calcul pour une représentation hydrodynamique réaliste. En effet, la formulation hyperbolique du problème à résoudre autorise l’utilisation d’un solveur volumes finis explicite parallèle rapide. Ce solveur bénéficie d’un outil de raffinement dynamique de maillage par blocs sur critère de production numérique d’entropie, qui permet d’optimiser le nombre d’inconnues, et donc le temps de calcul, tout en augmentant la précision de la simulation dans les zones d’intérêt. Ce modèle est appliqué à l’étude de la surverse d’une digue expérimentale. Le parement aval de la digue est considéré lisse ou bien avec un profil en marches représentatif des observations expérimentales. Grâce au raffinement dynamique de maillage, l’établissement de l’écoulement est simulé rapidement. Puis, lorsque que l’écoulement stationnaire est atteint, les hauteurs et champs de vitesse sont confrontés avec succès à quelques mesures expérimentales. Ces résultats sont également comparés avec des simulations réalisées avec un code commercial utilisant un modèle Navier-Stokes diphasique avec turbulence k-w SST