Simulation d’un écoulement de jet de rive par une méthode VOF.

Numerical simulations of turbulent bore run-up and run-down.

• Auteur / Autrice : Sylvain Mauriet
• Direction : Mathieu Mory, Stéphane Abadie
• Type : Thèse de doctorat
• Discipline(s) : Génie civil
• Date : Soutenance en 2009
• Etablissement(s) : Pau

Mots clés

Mots clés contrôlés

Littoraux -- Modifications

Modèles mathématiques

Turbulence

Ressaut hydraulique

Résumé

On présente les résultats de simulations RANS de la propagation d'un mascaret (obtenu par un ''lâcher de barrage'') sur une plage en pente et le run-up et le run-down ainsi générés. Les résultats numériques sont comparés aux résultats expérimentaux de Yeh et al. (1989). Les simulations ont été réalisées avec le code Navier-Stokes diphasique AQUILON. Deux méthodes de suivi d’interface VOF (VOF TVD ET VOF PLIC) sont implémentées. La viscosité turbulente est calculée par un modèle V2-F (Durbin, 1991). Une modélisation VOF-PLIC & V2-F est appliquée pour reproduire les caractéristiques macroscopiques du lâcher de barrage, qui comme on pouvait s’y attendre dépendent peu de la turbulence. Nous étudions aussi l���impact des conditions initiales sur k et  sur l’établissement de l’écoulement turbulent. Après ces validations vis-à-vis de la turbulence, des simulations du cas décrit par Yeh et al. (1989) sont menées pour optimiser le choix des paramètres de calcul. La théorie de Whitham (1958), prédit un effondrement du mascaret au niveau de la ligne de rivage au repos. La théorie de Shen and Meyer (1963) est toujours à l'heure actuelle le modèle de référence. Les résultats expérimentaux de Yeh et al. (1989) montrent clairement un phénomène différent. L’utilisation conjointe de la technique VOF-TVD et du modèle de turbulence V2-F semble apporter les meilleurs résultats par rapport aux expériences de Yeh et al. (1989). Une étude de la transition mascaret/lame de swash est proposée. Nos résultats montrent que la théorie de Whitham décrit de façon assez précise le mécanisme de d'effondrement du mascaret. Les résultats de nos simulations sont utilisés pour décrire la transition entre l'effondrement du mascaret et l'écoulement du run-up. L'analyse des processus de frottement dans le jet de rive met en évidence une forte dissymétrie entre le run-up et le run-down avec cisaillement plus faible lors du run-down.

The results of RANS numerical simulations of a bore propagating on a sloping beach and resulting in a run-up/run-down phenomenon are presented. Numerical results are compared to the results of a laboratory experiment by Yeh et al. (1989). Numerical simulations were carried out using the AQUILON model. Two different methods of interface tracking are implemented in the model to compute the interface between fluids (VOF PLIC and VOF TVD). The eddy viscosity is computed using the V2F model (Durbin, 1991). We study the bore generated from a dam-break condition. Shortly after its generation the bore propagates over a sloping beach. The geometry considered for numerical simulations is the same as for the experiments by Yeh et al. (1989). The theory of Whitham (1958), predicts the collapse of the bore at the shoreline. The theory of Shen and Meyer (1963) is up to date the reference model of the run-up flow. The experimental results of Yeh et al. (1989) clearly show a different phenomenon. It is discussed in the present paper. We examined in more detail the properties which may distinguish the bore collapse phenomenon from the subsequent run-up down. The results of our simulation show that the theory of Whitham does describe with sufficient detail the mechanism of bore collapse, due to the inherent shortcomings of shallow water equations. The results of our simulation are used to describe the transition between bore collapse and run-up flow. We analyse the evolution in time of the shear stress applied on the bottom along the swash water layer. We observe that the bottom shear stress is much lower during run-down as compared to run-up. The flow is highly asymmetrical between run-up and run-down.