L'objectif de ce travail est le développement d'un modèle et d'une méthode numérique adaptés à la simulation du ruissellement d'eau de pluie sur des surfaces agricoles. Pour cela, nous utilisons un système d'équations aux dérivées partielles pour les eaux peu profondes : le système de Saint Venant. La pluie et l'infiltration y sont pris en compte par l'ajout de termes source et le couplage avec un modèle d'infiltration. Le système étant hyperbolique et conservatif nous choisissons d'utiliser un schéma aux volumes finis avec reconstruction hydrostatique. Avec cette reconstruction nous obtenons un schéma permettant de traiter les équilibres stationnaires et les interfaces sec/mouillé dues aux événements pluvieux. Tout d'abord, nous effectuons une comparaison systématique de différents flux numériques, de différentes reconstructions d'ordre deux (MUSCL et ENO modifiées) et de différentes méthodes de traitement des frottements sur des solutions analytiques unidimensionnelles. Ces comparaisons nous permettent de choisir une méthode numérique adaptée à la simulation du ruissellement qui en outre capture les phénomènes de type roll-waves. Ensuite nous nous intéressons au traitement de la pluie et de l'infiltration à l'aide du modèle de Green-Ampt. Cette méthode généralisée en dimension deux est validée sur des résultats expérimentaux obtenus à l'INRA d'Orléans et à l'IRD. La méthode numérique choisie et validée a été implémentée dans FullSWOF_2D un logiciel libre écrit en C++ pour la simulation de ruissellement d'eau de pluie sur des surfaces agricoles.

Le système de Saint Venant est répandu pour modéliser des fluides dont la hauteur est inférieure au domaine d'écoulement. Son écriture nécessite des hypothèses sur le profil de vitesse pour connaître le flux de la quantité de mouvement ainsi que le cisaillement sur le fond. Dans cette thèse, nous nous sommes intéressés à un couplage entre un fluide parfait et une couche visqueuse dans l'esprit des couches limites interactives (IBL) introduites en aéronautique. Cette interaction nous permet de proposer un terme de friction en adéquation avec les attentes physiques au regard de la position du maximum local. Une part importante de cette thèse est donc consacrée à la compréhension de la couche visqueuse dans laquelle la recherche du profil de vitesse est cantonnée. Cette étude se décompose en l'écriture des équations de Prandtl puis en l'établissement de l'équation de von Kármán. Cette dernière met en jeu les quantités nécessaires à la définition du flux recherché et est donc un élément clé de la fermeture du système. Des résultats numériques viennent illustrer le modèle obtenu par le couplage entre le fluide parfait et la couche visqueuse. Le dernier chapitre expose deux formulations alternatives obtenues d'un point de vue d'un écoulement d'un fluide parfait dont les conditions sur les bords du domaine sont modifiées, soit par une condition de transpiration définie sur le fond, soit par une modification du domaine enlien avec une topographie apparente. / Shallow Water system is widely used for flows when the depth is smaller than the longitudinal scale. The establishment needs some hypothesis on the velocity profile in order to describe the moment flux and the shear stress on ground. In this thesis, we present a two layer decomposition of the fluid between an ideal fluid and a viscous layer in the spirit of the Interactive Boundary Layer (IBL) introduced in aeronautics. This interaction leads to obtain in our equations a friction term which fits with the physical expectations for the local maximum. So a major part of this work is interested in the comprehension of the viscous layer where the velocity profile is confined. The study is based on the writing of Prandtl equations then the establishment of the von Kármán equation. The last one contains the necessary quantities for a definition of the researched flux. Also this equation is essential for a closure of the system. Some numerical results illustrate the proposed model with the association of ideal fluid ans viscous layer. A last chapter presents two alternatives formulations of the model based on an ideal fluid with modified boundary conditions. The first one keeps the same domain but has a transpiration boundary.