Nous nous intéressons à la modélisation et à la simulation numérique d'écoulements gravitaires transitoires à surface libre, pour des fluides visqueux et incompressibles. La loi de comportement est de type viscoplastique avec transition fluide/solide. Plus précisément, nous considérons la loi rhéologique de Drucker-Prager. Nous nous intéressons tout d'abord le cas unidimensionnel d'un écoulement longitudinal cisaillé. Nous étudions un modèle simplifié, avec terme source empirique, pour lequel nous concevons une méthode numérique pour le suivi de la position de l'interface entre la couche solide et la couche fluide. Nous présentons des résultats numériques, avec divers termes sources, et nous comparons ces résultats, lorsque la viscosité est petite, à la solution analytique non visqueuse. Dans le cas visqueux, nous étudions les phases de démarrage et d'arrêt de l'écoulement. Dans un second temps, nous étudions le cas bidimensionnel d'écoulement de Drucker-Prager avec surface libre. La loi de comportement du fluide est traitée par régularisation, et nous utilisons la méthode ALE pour traiter le mouvement du domaine. Nous présentons des résultats numériques pour l'étude de la mise en mouvement d'un talus / This thesis deals with the modeling and numerical simulation of transient free-surface gravity flows, for viscous and incompressible fluids. The constitutive law is viscoplastic, with fluid/solid transition. More precisely, we consider the Drucker-Prager rheological law. We first study the case of a one-dimensional shear flow. We investigate a simplified model, with an empirical source term, for which we develop a numerical method to compute the position of the solid/fluid interface. We present numerical results for various source terms, and compare, in the case of small viscosity, our results to the inviscid analytical solution. In the viscous case, we study the case of a two-dimensional Drucker--Prager flow with free surface. The constitutive law of the fluid is regularized, and the ALE method is used to treat the displacement of the domain. Numerical results are presented for the setting in motion of an embankment

L'océan et l'atmosphère sont des milieux stratifiés, supports d'ondes internes de gravité. Ces ondes peuvent se propager sur de très grandes distances, transportant de l'énergie loin de leurs sources et jouant ainsi un rôle clé dans le bilan énergétique global des écoulements géophysiques. Notre approche expérimentale modélise de façon idéalisée les topographies océaniques et la stratification du milieu afin de comprendre les mécanismes fondamentaux mis en jeu en présence d'ondes internes. Cette démarche a été couplée au développement de techniques d'analyse des données telles que la démodulation complexe. Différents types d'ondes internes ont été utilisées. Pour cela, un point clé a été la caractérisation et le contrôle d'une nouvelle source d'ondes. Notre étude combinant expériences, simulations numériques et prédictions analytiques permet la génération contrôlée de structures spatiales variées (faisceaux localisés, modes verticaux occupant toute la hauteur du uide, ondes planes). L'étude du problème classique de la réflexion d'une onde plane sur un plan incliné a permis de vérifier l'absence d'onde retour pour des pentes proches de la criticalité, tout en mettant en évidence que ce cas d'école est encore un sujet ouvert. Nous avons également étudié la conversion d'un mode vertical modélisant la marée interne lors de l'interaction avec une topographie de grande échelle. Cette conversion est régie par des mécanismes linéaires et non-linéaires. La présence de pentes sous-critiques inuence la conversion du mode incident en de plus petites échelles spatiales, tandis que la hauteur de la topographie joue un rôle notable dans l'amplitude des modes transmis ainsi que l'intensité de processus non-linéaires générant des ondes harmoniques et un courant moyen. Nous avons réalisé la première mise en évidence expérimentale de la génération locale de trains d'ondes non-linéaires (solitons) lors de la réflexion d'un faisceau d'ondes internes intense et localisé au niveau d'une pycnocline. Ce phénomène dépend de la stratification considérée, qui caractérise la nature de la pycnocline. Enfin, une étude de la dynamique couplée entre un bateau évoluant à force constante et le fluide stratifié environnant a été menée. Le régime particulier des dépend d'un unique nombre sans dimension, le nombre de Froude associé à l'onde la plus rapide de la stratification sur laquelle le bateau évolue. Nous avons généralisé le cas classique à deux couches à tout type de stratification constituée d'une couche homogène en surface.

Ce travail numérique porte sur le comportement de la lutocline (interface entre l'eau et la vase fluide) en écoulement cisaillé et vise à une meilleure compréhension des mécanismes de remise en suspension de sédiments cohésifs. La crème de vase, ou vase partiellement solidifiée, est modélisée par un fluide homogène équivalent miscible dans l'eau, de rhéologie newtonienne ou viscoplastique. Une étude de l'hydrodynamique de cet écoulement stratifié en densité ainsi qu'en viscosité est ensuite proposée. Considérant une crème de vase initialement non-turbulente, l'étude se focalise sur le développement des instabilités au niveau de la lutocline et de la transition vers une couche de mélange turbulente. La particularité de cet écoulement réside dans la forte viscosité de vase et son seuil de mise en mouvement lorsqu'elle présente un caractère viscoplastique. Une étude de stabilité linéaire permet d'évaluer l'influence des différents paramètres de l'écoulement, notamment les stratifications en densité et en viscosité. La stratification en viscosité augmente sensiblement le taux de croissance de l'instabilité pour des nombres de Reynolds intermédiaires. L'évolution non-linéaire de l'écoulement est ensuite étudiée en utilisant des simulations numériques directes, la stratification en viscosité entrainant un épaississement de la couche de mélange finale. Enfin, des simulations numériques directes basées sur un modèle de fluide de Bingham régularisé permettent d'étudier l'influence de la contrainte seuil sur le développement de l'instabilité. / This numerical study focuses on the behavior of the lutocline in a shear flow and aims to better understand the mechanism of resuspension of cohesive sediment. Mud flow, or mud partially consolidated, is modeled by an equivalent homogenous fluid miscible in water, with newtonian or viscoplastic rheology. A study of the hydrodynamics of this shear flow, stratified both in density and viscosity, is presented. Considering an initially laminar mud flow, the focus of the study is on the development of instabilities on the lutocline and the transition to a turbulent mixing layer. The specificity of this flow lies on the large viscosity of the mud and its threshold to be put in motion, when it presents a viscoplastic feature. A linear stability study assesses the influence of the various parameters of the flow, especially of density and viscosity stratification. The viscosity stratification slightly increases the growth rate of the instability for intermediate Reynolds numbers. Then, the non linear evolution of the flow is studied by using direct numerical simulations, viscous stratification leading to a thicker mixing layer. At last, direct numerical simulations based on a Bingham regularized model, permits to study the influence of the critical strain on the development of the instability.