On dérive un nouveau modèle d’écoulements cisaillés et turbulents d’eau peu profonde. Les écarts de la vitesse horizontale par rapport à sa valeur moyenne sont pris en compte par une nouvelle variable appelée enstrophie, liée à la vorticité et à l’énergie turbulente. Le modèle comporte trois équations qui sont les bilans de masse, de quantité de mouvement et d’énergie. Le modèle est hyperbolique et peut être écrit sous forme conservative. L’énergie turbulente, dont l’intensité peut être importante, est produite par les ondes de choc qui apparaissent naturellement dans le modèle. Les écoulements rapidement variés étudiés sont caractérisés par l’existence d’une structure turbulente appelée rouleau dans laquelle la dissipation d’énergie turbulente joue un rôle majeur. Cette dissipation, qui détermine notamment le profil de profondeur, est modélisée par l’introduction d’un terme nouveau dans le bilan d’énergie. Le modèle comporte deux paramètres. L’un gouverne la dissipation de l’énergie turbulente du rouleau. L’autre paramètre, l’enstrophie de paroi, liée au cisaillement sur le fond, peut être considéré comme constant dans la partie rapidement variée d’un écoulement, sur laquelle il exerce une influence assez faible. Ce modèle a été appliqué avec succès aux vagues des trains de rouleaux et au ressaut hydraulique classique. Le profil de la surface libre est en très bon accord avec les résultats expérimentaux. L’étude numérique en régime non stationnaire permet notamment de prédire le régime oscillatoire du ressaut hydraulique. La fréquence d’oscillations correspondante est en accord satisfaisant avec les mesures expérimentales de la littérature. / We derive a new model of turbulent shear shallow water flows. The deviation of the horizontal velocity from its average value is taken into account by a new variable called enstrophy, which is related to the vorticity and to the turbulent energy. The model consists of three equations which are the balances of mass, momentum and energy. The model is hyperbolic and can be written in conservative form. The turbulent energy, which can be of high intensity, is produced in shock waves which appear naturally in the model. The rapidly varied flows we studied are characterized by the presence of a turbulent structure called roller in which the turbulent energy dissipation plays a major part. This dissipation, which determines, in particular, the depth profile, is modelled by the introduction of a new term in the energy balance equation. The model contains two parameters. The first one governs the dissipation of the turbulent energy of the roller. The second one, the wall enstrophy, related to the shearing at the bottom, can be considered as constant in the rapidly varied part of the flow on which it does not exert an important influence. This model was successfully applied to roll waves and to the classical hydraulic jump. The free surface profile was found in very good agreement with the experimental results. The numerical study in the non-stationary case can notably predict the oscillations of the hydraulic jump. The corresponding oscillation frequency is in good agreement with the experimental measures found in the literature.

Cette thèse porte sur la résolution numérique du modèle multi-dimensionnel d'écoulement cisaillé en eau peu profonde. Dans le cas d'un mouvement unidimensionnel, ces équations coïncident avec les équations de la dynamique de gaz pour un choix particulier de l'équation d'état. Dans le cas multi-dimensionnel, le système est complètement différent du modèle de la dynamique de gaz. Il s'agit d'un système EDP hyperbolique 2D non-conservatif qui rappelle un modèle de turbulence barotrope. Le modèle comporte trois types d'ondes correspondant à la propagation des ondes de surface, des ondes de cisaillement et à celle de la discontinuité de contact. Nous présentons dans le cas 2D un schéma numérique basé sur une nouvelle approche de ``splitting" pour les systèmes d'équations non-conservatives. Chaque sous-système ne contient qu'une seule famille d'ondes: ondes de surface ou ondes de cisaillement, et discontinuité de contact. La précision d'une telle approche est testée sur des solutions exactes 2D décrivant l'écoulement lorsque la vitesse est linéaire par rapport aux variables spatiales, ainsi que sur des solutions décrivant des trains de rouleaux 1D. Finalement, nous modélisons un ressaut hydraulique circulaire formé dans un écoulement convergent radial d'eau. Les résultats numériques obtenus sont clairement similaires à ceux obtenus expérimentalement: oscillations du ressaut et son rotation avec formation du point singulier. L'ensemble des validations proposées dans ce manuscrit démontre les aptitudes du modèle et de la méthode numérique pour la résolution des problèmes complexes d'écoulements cisaillés en eau peu profonde multidimensionnels. / This thesis is devoted to the numerical modelling of multi-dimensional shear shallow water flows. In 1D case, the corresponding equations coincide with the equations describing non--isentropic gas flows with a special equation of state. However, in the multi-D case, the system differs significantly from the gas dynamics model. This is a 2D hyperbolic non-conservative system of equations which is reminiscent of a generic Reynolds averaged model of barotropic turbulent flows. The model has three families of characteristics corresponding to the propagation of surface waves, shear waves and average flow (contact characteristics). First, we show the ability of the one-dimensional conservative shear shallow water model to predict the formation of roll-waves from unstable initial data. The stability of roll waves is also studied.Second, we present in 2D case a new numerical scheme based on a splitting approach for non-conservative systems of equations. Each split subsystem contains only one family of waves (either surface or shear waves) and contact characteristics. The accuracy of such an approach is tested on exact 2D solutions describing the flow where the velocity is linear with respect to the space variables, and on the solutions describing 1D roll waves. Finally, we model a circular hydraulic jump formed in a convergent radial flow of water. Obtained numerical results are qualitatively similar to those observed experimentally: oscillation of the hydraulic jump and its rotation with formation of a singular point. These validations demonstrate the capability of the model and numerical method to solve challenging multi--dimensional problems of shear shallow water flows.