Étude numérique des interactions multi-échelles écoulement-sédiment-structure par une approche multiphasique

Numerical study of multi-scale flow-sediment-structure interactions using a multiphase approach.

• Auteur / Autrice : Tim Nagel
• Direction : Achim Wirth, Julien Chauchat
• Type : Thèse de doctorat
• Discipline(s) : Mécanique des fluides Energétique, Procédés
• Date : Soutenance le 17/07/2018
• Etablissement(s) : Université Grenoble Alpes (ComUE)
• Ecole(s) doctorale(s) : École doctorale Ingénierie - matériaux mécanique énergétique environnement procédés production (Grenoble ; 2008-....)
• Partenaire(s) de recherche : Laboratoire : Laboratoire des écoulements géophysiques et industriels (Grenoble)
• Jury : Président / Présidente : Philippe Gondret
• Jury : Examinateurs / Examinatrices : Cyrille Bonamy, Tian Jian Hsu, Xiaofeng Liu
• Jury : Rapporteurs / Rapporteuses : Sylvain Guillou, Ming Li

Mots clés

Mots clés contrôlés

Simulation par ordinateur

Sédiments (géologie) -- Transport

Mots clés libres

Transport de sédiment

Turbulence

Écoulement diphasique

Simulation numérique

Sediment transport

Turbulence

Numerical simulation

Two-Phase flow

Résumé

Le travail réalisé dans cette thèse a consisté en le développement et l'utilisation des modèles numériques pour étudier les interactions multi-échelles entre une éolienne offshore et la dynamique locale océanique et sédimentaire. Dans une première partie, les interactions entre le système couplé océan-sédiment et le sillage atmosphérique généré par une turbine éolienne offshore sont étudiées à l'aide d'un modèle numérique 2D développé au cours de la thèse et écrit en fortran. Ce modèle résout les équations de Barré-De-Saint-Venant pour l'océan et l'équation d'Exner pour le sédiment. Dans un seconde partie, le phénomène d'affouillement 3D autour d'un cylindre vertical est étudié à l'aide d'un modèle diphasique eulérien-eulérien, sedFoam, implémenté dans la boîte à outils numériques OpenFOAM. L'approche diphasique permet de tenir compte des processus de petite échelle en s'affranchissant des hypothèses classiquement faites pour la modélisation du transport sédimentaire, notamment la corrélation locale entre le flux de sédiments et la contrainte de cisaillement fluide sur le fond.Concernant l'impact du sillage atmosphérique généré par une turbine, nous avons montré que celui-ci peut générer des allées tourbillonnaires dans l'océan. La dynamique turbulente océanique est alors contrôlée par le paramètre de sillage S=Cd D/H, où D est le diamètre du sillage au point d'impact sur la surface de l'océan, Cd est le coefficient de la loi de friction quadratique entre l'océan et le fond et H la profondeur de l'océan. Une paramétrisation des flux turbulents basée sur S est proposée pour modéliser la dynamique océanique dans des modèles à plus grande échelle de type RANS (Reynolds Averaged Navier-Stokes). Les résultats montrent que la dynamique océanique a une rétro-action sur la puissance du vent disponible. Les résultats montrent également que la dynamique sédimentaire instantanée est couplée à la dynamique océanique. Cependant, les variations de l'élévation du fond marin sont faibles (mm/mois) et l'impact morphodynamique du sillage est négligeable.Concernant la simulation diphasique de l'affouillement, après une validation du modèle sur des configurations 1D et 2D, des simulations tridimensionnelles autour d'une pile cylindrique sont présentées. Dans un premier temps, une configuration sans sédiments est réalisée afin de valider la capacité du modèle de turbulence URANS (Unsteady Reynolds Averaged Navier-Stokes) développé dans ce travail de thèse à reproduire les structures tourbillonnaires responsables de l'affouillement comme le tourbillon en fer à cheval et le lâché tourbillonnaire à l'aval du cylindre. Ensuite, les premières simulations diphasiques 3D de l'affouillement autour du cylindre ont été réalisées en régime de transport de type lit-mobile. Ces simulations constitue un véritable challenge en terme calcul numérique à haute performance. La comparaison favorable des résultats de simulations avec les résultats expérimentaux de la littérature apporte la preuve de concept que l'approche diphasique est pertinente pour étudier des configurations d'écoulements complexes instationnaire et tridimensionnelle. Les résultats de simulation sont ensuite analysés pour étudier la relation entre le flux local de transport de sédiments, la valeur de la contrainte fluide sur le fond et la pente locale du lit sédimentaire. La déviation par rapport aux résultats obtenus en écoulement uniforme permet d'identifier les mécanismes prépondérant de transport associées au tourbillon en fer à cheval, à la pente de fond et aux tourbillons lâchés dans le sillage du cylindre. Les résultats obtenus montrent une sensibilité à la résolution numérique en particulier à l'aval du cylindre illustrant le besoin de réaliser des simulations des grandes échelles turbulentes diphasiques.

The work undertaken in this PhD thesis was to develop and use numerical models to investigate the multi-scale interactions between an offshore wind turbine and the local ocean and sediment dynamics. First, the interactions between the coupled ocean-sediment system and the atmospheric wake generated by an offshore wind turbine are investigated using an idealized two-dimensional model developed during this Phd thesis and written in fortran. The model integrates the shallow water equations for the ocean together with the Exner equation for the sediment bed. In a second part, the 3D scour phenomenon around a vertical cylinder in a steady current is studied using a two-phase flow eulerian-eulerian solver, sedFoam, written within the framework of the numerical toolbox OpenFOAM. The two-phase flow approach accounts for small-scale processes by avoiding the traditional assumptions made for sediment transport modeling, such as a local corre- lation between the sediment flux and the fluid bed shear stress.Regarding the atmospheric wake generated by a turbine, the results shows that its impact on the ocean’s surface can generate vortices. The resulting turbulent ocean dynamics is controlled by the wake parameter S = CdD/H, where D is the wake diameter at the impact location on the ocean surface, Cd is the quadratic friction coefficient between the ocean and the sediment and H is the oceanic layer depth. A turbulence parameterization based on S is proposed, allowing for upscaling simulations in larger scales Reynolds Averaged Navier-Stokes (RANS) models. It is shown that the ocean dynamics has an effect on the available wind power. The results also show that the instantaneous sediment dynamics is strongly coupled with the ocean one but that the overall seabed elevation variations remain small (a few millimeters/month). The morphodynamic impact of the wake is thus negligible.Concerning the two-phase flow simulation of scour, sedFoam is first validated on 1D and 2D configurations. Then, 3D simulations around a vertical cylindrical pile are presented. At first, a validation of the Unsteady Reynolds Averaged Navier-Stokes (URANS) turbulence model developed in this work is performed on a configuration without sediment. The results show that the vortices structures responsible for scouring, the Horse Shoe Vortex (HSV) and the vortex-shedding in the lee of the cylinder are correctly reproduced. Then, 3D two-phase flow simulations of the scour around a cylindrical pile have been carried out in a live-bed configuration. This work is the first attempt to model 3D scour phenomenon using the two-phase flow approach. Such simulations represent a real challenge in terms of high performance computing. The good agreement between the numerical predictions and the literature experimental results provide the proof of concept that the two-phase flow approach can be used to study complex 3D and unsteady flow configurations. The relationship between the local bed shear stress, the sediment flux and the local sediment bed slope is further investigated. The deviation of the results from a uniform flow configuration is further analyzed to identify the relevant sediment transport mechanisms associated with the HSV, the slope in the scour mark and the vortex-shedding downstream of the cylinder. Finally, the numerical results show a grid sensitivity of the morphological predictions in the lee of the cylinder that are most probably related to small-scale resolved vortical structures. This highlights the need for two-phase flow Large Eddy Simulations on this configuration in the future.