Accurate models to reliably predict sediment and chemical transport in watershed water systems enhance the ability of environmental scientists, engineers and decision makers to analyze the impact of contamination problems and to evaluate the efficacy of alternative remediation techniques and management strategies prior to incurring expense in the field. This dissertation presents the conceptual and mathematical development of a general numerical model simulating (1) sediment and reactive chemical transport in river/stream networks of watershed systems; (2) sediment and reactive chemical transport in overland shallow water of watershed systems; and (3) reactive chemical transport in three-dimensional subsurface systems. Through the decomposition of the system of species transport equations via Gauss-Jordan column reduction of the reaction network, fast reactions and slow reactions are decoupled, which enables robust numerical integrations. Species reactive transport equations are transformed into two sets: nonlinear algebraic equations representing equilibrium reactions and transport equations of kinetic-variables in terms of kinetically controlled reaction rates. As a result, the model uses kinetic-variables instead of biogeochemical species as primary dependent variables, which reduces the number of transport equations and simplifies reaction terms in these equations. For each time step, we first solve the advective-dispersive transport of kinetic-variables. We then solve the reactive chemical system node by node to yield concentrations of all species. In order to obtain accurate, efficient and robust computations, five numerical options are provided to solve the advective-dispersive transport equations; and three coupling strategies are given to deal with the reactive chemistry. Verification examples are compared with analytical solutions to demonstrate the numerical accuracy of the code and to emphasize the need of implementing various numerical options and coupling strategies to deal with different types of problems for different application circumstances. Validation examples are presented to evaluate the ability of the model to replicate behavior observed in real systems. Hypothetical examples with complex reaction networks are employed to demonstrate the design capability of the model to handle field-scale problems involving both kinetic and equilibrium reactions. The deficiency of current practices in the water quality modeling is discussed and potential improvements over current practices using this model are addressed.

Nous analysons dans cette thèse différents aspects associés à la modélisation des écoulements à surface libre en eaux peu profondes (Shallow Water). Nous étudions tout d’abord le système d’équations de Saint-Venant à deux dimensions et leur résolution par la méthode numérique des volumes finis, en portant une attention particulière sur les aspects hyperboliques et conservatifs. Ces schémas permettent de traiter les équilibres stationnaires, les interfaces sec/mouillé et aussi de modéliser des écoulements subcritique, transcritique et supercritique. Nous présentons ensuite la théorie de la méthode d’assimilation variationnelle de données adaptée à ce type d’écoulement. Son application au travers des études de sensibilité est longuement discutée dans le cadre de l'hydraulique à surface libre. Après cette partie à caractère théorique, la partie tests commence par une qualification de l’ensemble des méthodes numériques qui sont implémentées dans le code DassFlow, développé à l’Université de Toulouse, principalement à l’IMT mais aussi à l’IMFT. Ce code résout les équations Shallow Water par une méthode de volumes finis et est validé par comparaison avec les solutions analytiques pour des cas tests classiques. Ces mêmes résultats sont comparés avec un autre code d’hydraulique à surface libre aux éléments finis en deux dimensions, Telemac 2D. Une particularité notable du code DassFlow est de permettre l’assimilation variationnelle de données grâce au code adjoint permettant le calcul du gradient de la fonction coût. Ce code adjoint a été obtenu en utilisant l'outil de différentiation automatique Tapenade (Inria). Nous testons ensuite sur un cas réel, hydrauliquement complexe, différentes qualités de Modèles Numériques de Terrain (MNT) et de bathymétrie du lit d’une rivière. Ces informations proviennent soit d’une base de données classique type IGN, soit d’informations LIDAR à très haute résolution. La comparaison des influences respectives de la bathymétrie, du maillage et du type de code utilisé, sur la dynamique d’inondation est menée très finement. Enfin nous réalisons des études cartographiques de sensibilité aux paramètres du modèle sur DassFlow. Ces cartes montrent l’influence respective des différents paramètres ou de la localisation des points de mesure virtuels. Cette localisation optimale de ces points est nécessaire pour une future assimilation de données efficiente. / We analyze in this thesis various aspects associated with the modeling of free surface flows in shallow water approximation. We first study the system of Saint-Venant equations in two dimensions and its resolution with the numerical finite volumes method, focusing in particular on aspects hyperbolic and conservative. These schemes can process stationary equilibria, wetdry interfaces and model subcritical, transcritical and supercritical flows. After, we present the variational data assimilation method theory fitted to this kind of flow. Its application through sensitivity studies is fully discussed in the context of free surface water. After this theoretical part, we test the qualification of numerical methods implemented in the code Dassflow, developed at the University of Toulouse, mainly at l'IMT, but also at IMFT. This code solves the Shallow Water equations by finite volume method and is validated by comparison with analytical solutions for standard test cases. These results are compared with another hydraulic free surface flow code using finite elements in two dimensions: Telemac2D. A significant feature of the Dassflow code is to allow variational data assimilation using the adjoint method for calculating the cost function gradient. The adjoint code was obtained using the automatic differentiation tool Tapenade (INRIA). Then, the test is carried on a real hydraulically complex case using different qualities of Digital Elevation Models (DEM) and bathymetry of the river bed. This information are provided by either a conventional database types IGN or a very high resolution LIDAR information. The comparison of the respective influences of bathymetry, mesh size, kind of code used on the dynamics of flooding is very finely explored. Finally we perform sensitivity mapping studies on parameters of the Dassflow model. These maps show the respective influence of different parameters and of the location of virtual measurement points. This optimal location of these points is necessary for an efficient data assimilation in the future.

On s’intéresse dans ce travail à l’aspect tridimensionnel (3D) de la propagation sonore en milieux océaniques petits fonds dans le cadre des problèmes inverses. Les problèmes inverses en acoustique sous-marine se basent sur la modélisation bidimensionnelle (2D) de la propagation, ignorant ainsi les effets de réfraction horizontale, qualifiés d’effets 3D. Toutefois, la propagation acoustique en environnements petits fonds, tels le plateau continental, peut être affectée par des effets 3D, leur prise en compte nécessitant l’utilisation de modèles pleinement 3D. Une inversion basée sur un modèle 3D devient inabordable pour plus de deux paramètres à la fois en raison de temps CPU particulièrement élevés. L’objectif de ce travail est d’examiner l’importance des effets 3D sur la performance et la fiabilité des procédures d’inversion habituellement utilisées dans les problèmes de l’acoustique sous-marine. Pour cela, on se place dans un guide d’onde océanique à géométrie inclinée. Des expérimentations à échelle réduite sont menées afin d’identifier et d’interpréter les effets 3D. Une procédure d’inversion par champs d’onde adaptés, formulée dans un cadre Bayesien et basée sur la recherche exhaustive dans l’espace des paramètres, est élaborée. L’inversion s’effectue en comparant des données basse fréquence du champ acoustique, recueillies le long d’antennes linéaires verticales ou horizontales, aux répliques générées par des modèles d’équation parabolique 2D et 3D. Les paramètres importants sont identifiés au moyen d’une étude de sensibilité de la fonction de coût. Dans une étape préliminaire, la performance de l’inversion est étudiée, sur données synthétiques bruitées, dans un cas simple permettant l’utilisation de modèles 2D. Une stratégied’inversion en sous-espaces résultant en une réduction importante des temps CPU pour l’inversion, est examinée. L’inversion est ensuite abordée en présence d’un fond incliné. La possibilité et les limites d’une inversion basée sur un modèle 2D sont explorées. Cette approche, mise en œuvre sur données synthétiques, met en évidence la pertinence de l’utilisation de modèles 2D en champ relativement proche. Une inversion basée sur un modèle 3D n’étant alors nécessaire que pour la pente, des temps CPU raisonnables sont ainsi réalisés. En revanche, en champ lointain, nous sommes confrontés à un risque potentiel d’estimation erronée et le recours à une modélisation 3D devient nécessaire. / This work deals with the three-dimensional (3D) aspect of sound propagation in shallow-water oceanic environments with respect to inverse problems. Inverse problems in underwater acoustics are based on twodimensional (2D) modeling of sound propagation, hence ignoring the effects of horizontal refraction, referred to as 3D propagation effects. However, the acoustic propagation in shallow-water environments, like the continental shelf, may be affected by 3D effects requiring 3D modeling to be accounted for. An inversion based on a 3D model for more than two parameters at a time becomes prohibitive due to dramatically increased CPU times. The aim of this work is to investigate the importance of the 3D effects with respect to the performance and reliability of the inversion procedures typically applied in problems of underwater acoustics. To this aim, we focus on a wedge-shaped oceanic wave guide. Laboratory scale experiments of long-range acoustic propagation are performed to identify and interpret the 3D effects due to a sloping bottom, as predicted by numerical simulations. A matched-field inversion procedure implemented within a Bayesian framework and based on the exhaustive search over the parameter space is elaborated. The inversion is performed by comparing low frequency acoustic field data, collected along vertical or horizontal line arrays, to replica generated from 2D and 3D parabolic equation codes. The recoverable parameters are identified by means of a sensitivity study of the cost function. In a preliminary step, the inversion performance is investigated on noisy synthetic data in a simple waveguide where 2D codes apply. A ubspace inversion strategy providing significant reduction in CPU times is examined. The inversion in the presence of a sloping bottom is then considered. The feasibility and the limits of an inversion matching replica from a 2D code are explored. This approach, applied on synthetic data, highlights the relevance of using 2D codes at relatively short ranges. An inversion based on a 3D code is thus only needed for the slope, and reasonable CPU times are achieved. On the other hand, important mismatch might occur at farther ranges and 3D modeling is required.

Dans la production d'électricité, un des leviers centraux pour réduire les détériorations et les pertes causées par l'humidité dans les turbines à vapeur est l'étude des films liquides. Ces films minces, sont créés par la déposition de gouttes et sont fortement cisaillés. Des gouttes peuvent ensuite être arrachées du film. A l'heure actuelle, aucun modèle complet et valide n'existe pour décrire ce phénomène. Un modèle 2D à formulation intégrale associé à des lois de fermetures a été dérivé pour représenter ce film. Comparé aux équations classiques de Saint-Venant, le modèle prend en compte davantage d'effets : le transfert de masse, l'impact des gouttes, le cisaillement à la surface libre, la tension de surface, le gradient de pression et la rotation. Une analyse des propriétés du modèle (hyperbolicité, entropie, conservativité, analyse de stabilité linéaire, invariance par translation et par rotation) est réalisée pour juger de la pertinence du modèle. Un nouveau code 2D est implémenté dans un module de développement libre du code EDF Code Saturne et une méthode de volumes finis pour un maillage non-structure a été développée. La vérification du code est ensuite effectuée avec des solutions analytiques dont un problème de Riemann. Le modèle, qui dégénère en modèle classique de Saint-Venant pour le cas d'un film tombant sur un plan inclinée, est validé par l'expérience de Liu and Gollub, 1994, PoF et comparé à des modèles de références (Ruyer-Quil and Manneville, 2000, EPJ-B et Lavalle, 2014, PhD thesis). Un autre cas d'étude met en scène un film cisaillé en condition basse-pression de turbine à vapeur et, est validé par l'expérience de Hammitt et al., 1981, I. Enfin, le code film est couplé aux données 3D du champ de vapeur autour d'un stator d'une turbine basse-pression du parc EDF, issues de Blondel, 2014, PhD thesis. Cette application industrielle montre la faisabilité d'une simulation d'un film en condition réelle du turbine à vapeur. / In the electricity production, one central key to reduce damages and losses due to wetness in steam turbines is the study of liquid films. These thin films are created by the deposition of droplets and are highly sheared. This film may then be atomized into coarse water. At the moment, no comprehensive and validated model exists to describe this phenomenon. A 2D model based on a integral formulation associated with closure laws is developed to represent this film. Compared to classical Shallow-Water equation, the model takes into account additional effect : mass transfer, droplet impact, shearing at the free surface, surface tension, pressure gradient and the rotation. The model properties (hyperbolicity, entropy, conservativity, linear stability, Galilean invariance and rotational invariance) has been analyzed to judge the pertinence of the model. A new 2D code is implemented in a free module of the code EDF Code Saturne and a finite volume method for unstructured mesh has been developed. The verification of the code is then carried out with analytical solutions including a Riemann problem. The model, which degenerates into classical Shallow-Water equations for the case of a falling liquid film on a inclined plane, is validated by the experiment of Liu and Gollub, 1994, PoF and compared to reference models (Ruyer-Quil and Manneville, 2000, EPJ-B et Lavalle, 2014, PhD thesis). Another study depicts a sheared film under low-pressure steam turbine conditions and is validated by the experiment of Hammitt et al., 1981, FiI. Lastly, the code film is coupled to 3D steam data around a fixed blade of a BP100 turbine, from Blondel, 2014, PhD thesis. This industrial application shows the feasibility of liquid film's simulation in real steam turbine condition.

Ce travail de thèse, propose un schéma numérique d'ordre élevé, distribuantle résidu (RD) pour l'approximation d'un problème hydro-sédimentairehyperbolique non conservatif, couplant les modèles de Grass et de Saint-Venant. Il fait appel à des méthodes de Runge-Kutta à variation totale diminuanteet de stabilisation (méthode de décentrement amont, dit Upwind),avec ou sans adjonction de limiteurs et présente de bonnes propriétés.L'une des facettes importantes de ce qui a été réalisée, repose sur la conceptionet le développement d'un programme Python 2D-espace, sous la formed'un logiciel faisant appel à un ensemble de modules créés pour l'occasion.Le développement du code de calcul, qui se propose d'approcher la solutiondu problème hydro-sédimentaire, a été e_ectué avec une orientation Objetet pour être e_cace sur calculateur parallèle (utilisant le parallélisme multithreadsOpenMP). L'une des particularités du schéma numérique dans cecadre, est liée à son application à des quadrangles.Un programme 1D-espace, qui se présente également sous forme de logiciel,a aussi été mis en place. Pour des raisons de portabilité et d'e_catité, il aété écrit multilangages (Python-Fortran : via numpy.ctypes pour Python etvia l'interface standard de Fortran pour C). Le schéma RD avec ou sansadjonction de limiteurs de _ux, a été implémenté à la manière d'un schémaprédicteur-correcteur. Des comparaisons avec d'autres schémas ont été e_ectuées a_n de montrer son e_cacité, son ordre de précision élevé a été mis enévidence, et la C-propriété a été testée. Les tests ont révélé que, pour le casd'un transport d'un pro_l sédimentaire parabolique, c'est le limiteur de _uxMUSCL MinMod, qui est le plus adapté parmi ceux testés.Dans le cas scalaire, des tests numériques ont été réalisés a_n de validerle second ordre de précision. / The present work, proposes a high order Residual Distribution (RD) numericalscheme to solve the non conservative hyperbolic problem, coupling Shallow Water and Grass equations. It uses Total Value Diminishing Runge Kutta and stabilisation Upwind methods, with or without limiters. It also has some good properties.A part of the work realised in this thesis, is about the conception and the developpement of a 2D-space Python program, under the form of a software,using a set of moduls created for the occasion. the code developpement, whichis said to approach the _uid-sediment model, coupling Shallow-Water and sedimentequations, has been made with an Object orientation and in orderto be e_cient on parallel architecture (using multithreads OpenMP parallelism). One of the features of the scheme in this case, is due to its application on quadrangles.A 1D-space program, also writen as a software, has been estabished. In order to be portable and e_cient, It has been developped multilinguals (Python- Fortran : by numpy.ctypes for Python and by standart interface FORTRAN for C). The RD scheme with or without Flux Limiters, has been implemented like predictor-corrector one. Comparisons with other schemes results have been realised, in order to show its e_ciency, moreover its high order accuracy has been focus on, and the C-proprerty has been tested. The tests show that MUSCL MinMod _ux limiters, is the most adaptated for a dune test case, between all tested.In the scalar case, numerical tests have been realised, for validating the secondorder of accuracy.

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.

There are two main subjects in this thesis: the first one deals with the numerical simulation of shallow water equations, the other one is the resolution of some problems in image processingThe first part of this dissertation is devoted to the shallow waters. We propose a combined scheme which uses the Marquina's double flux decomposition (extended to the non homogeneous case) when adjacent states are not close and a single decomposition otherwise. This combined scheme satisfies the exact C property. Furthermore, we propose a special treatment of the numerical scheme at dry zones.The second subject is the digital simulation of the Day for Night (or American Night in Europe). The proposed algorithm simulates a night image coming from a day image and considers some aspects of night perception. In order to simulate the loss of visual acuity we introduce a partial differential equation that simulates the spatial summation principle of the photoreceptors in the retina.The gap restoration (inpainting) on surfaces is the object of the third part. For that, we propose some geometrical approaches based on the mean curvature. Then, we also use two interpolation methods: the resolution of the Laplace equation, and an Absolutely Minimizing Lipschitz Extension (AMLE). Finally, we solve the restoration problem of satellite images. The variational problem that we propose manages to do irregular to regular sampling, denoising, deconvolution and zoom at the same time. / Los temas tratados en esta tesis son, por un lado, la simulación numérica de las ecuaciones de aguas someras ("shallow waters") y por otro, la resolución de algunos problemas de procesamiento de imágenes. En la primera parte de la tesis, dedicada a las aguas someras, proponemos un esquema combinado que usa la técnica de doble descomposición de flujos de Marquina (extendida al caso no homogéneo) cuando los dos estados adyacentes no están próximos y una única descomposición en caso contrario. El esquema combinado verifica la propiedad C exacta. Por otro lado, proponemos un tratamiento especial en las zonas secas.El segundo tema tratado es la simulación digital de la Noche Americana ("Day for Night"). El algoritmo propuesto simula una imagen nocturna a partir de una imagen diurna considerando varios aspectos de la percepción visual nocturna. Para simular la pérdida de agudeza visual se propone una ecuación en derivadas parciales que simula el principio de sumación espacial de los fotoreceptores situados en la retina.La restauración de agujeros ("inpainting") en superficies es objeto de la tercera parte. Para ello se proponen varios enfoques geométricos basados en la curvatura media. También se utilizan dos métodos de interpolación: la resolución de la ecuación de Laplace y el método AMLE (Absolutely Minimization Lipschitz Extension).Por último, tratamos la restauración de imágenes satelitales. El método propuesto consigue obtener una colección de muestras regulares a partir de un muestreo irregular, eliminando a la vez el ruido, deconvolucinando la imagen y haciendo un zoom.

L’objectif global de ce travail est d’étudier une modélisation multi échelle et de développer une méthode adaptée pour la simulation numérique du processus d’érosion à l’échelle du bassin versant. Après avoir passé en revue les différents modèles existants, nous dérivons une solution analytique non triviale pour le système couplé modélisant le transport de sédiments par charriage. Ensuite, nous étudions l’hyperbolicité de ce système avec diverses lois de sédimentation proposées dans la littérature. Concernant le schéma numérique, nous présentons le domaine de validité de la méthode de splitting, pour les équations modélisant l’écoulement et celle décrivant l’évolution du fond. Pour la modélisation du transport en suspension à l’échelle de la parcelle, nous présentons un système d’équations couplant les mécanismes d’infiltration, de ruissellement et le transport de plusieurs classes de sédiments. L’implémentation et des tests de validation d’un schéma d’ordre élevé et de volumes finis bien équilibré sont également présentés. Ensuite, nous discutons sur l’application et la calibration du modèle avec des données expérimentales sur dix parcelles au Niger. Dans le but d’aboutir la simulation à l’échelle du bassin versant, nous développons une modélisation multi échelle dans laquelle nous intégrons le taux d’inondation dans les équations d’évolution afin de prendre en compte l’effet à petite échelle de la microtopographie. Au niveau numérique, nous étudions deux schémas bien équilibrés : le schéma de Roe basé sur un chemin conservatif, et le schéma avec reconstruction hydrostatique généralisée. Enfin, nous présentons une première application du modèle avec les données expérimentales du bassin versant de Ganspoel qui nécessite la parallélisation du code. / The overall objective of this thesis is to study a multiscale modelling and to develop a suitable method for the numerical simulation of soil erosion on catchment scale. After reviewing the various existing models, we derive an analytical solution for the non-trivial coupled system modelling the bedload transport. Next, we study the hyperbolicity of the system with different sedimentation laws found in the literature. Relating to the numerical method, we present the validity domain of the time splitting method, consisting in solving separately the Shallow-Water system (modelling the flow routing) during a first time step for a fixed bed and updating afterward the topography on a second step using the Exner equation. On the modelling of transport in suspension at the plot scale, we present a system coupling the mechanisms of infiltration, runoff and transport of several classes of sediment. Numerical implementation and validation tests of a high order wellbalanced finite volume scheme are also presented. Then, we discuss on the model application and calibration using experimental data on ten 1 m2 plots of crusted soil in Niger. In order to achieve the simulation at the catchment scale, we develop a multiscale modelling in which we integrate the inundation ratio in the evolution equations to take into account the small-scale effect of the microtopography. On the numerical method, we study two well-balanced schemes : the first one is the Roe scheme based on a path conservative, and the second one is the scheme using a generalized hydrostatic reconstruction. Finally, we present a first model application with experimental data of the Ganspoel catchment where the parallel computing is also motived.

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.

Ce travail porte sur les processus hydrodynamiques en zone littorale.Deux principaux thèmes sont abordés. Le premier concerne la couchelimite oscillante provoquée par l’interaction entre les vagues et le fondà l’approche des côtes. Le second traite de l’évolution de la circulationet la vorticité induite par la bathymétrie et/ou le forçage des vagues.Un modèle de couche limite turbulente a été élaboré et utilisé pourobserver l’évolution de la couche limite oscillante sous l’effet de vaguesnon-linéaires, en s’appuyant sur une modélisation physique menéedans le canal à houle du LEGI. Les profils expérimentaux de vitesseet positions du fond fixe instantanés permettent de définir l’évolutiondes non-linéarités induites par les vagues au sein de la couche limite.Le modèle numérique couplé à une modélisation du mouvement dulit mobile est capable de reproduire l’évolution de ces non-linéarités,et explique que la diffusion verticale observée expérimentalementest principalement due au mouvement vertical du lit causé induitpar les vagues. Pour l’étude de la circulation et de la vorticité enzone côtière, un modèle numérique 2D moyenné sur la verticale detype Shallow Water est validé avec les données d’une expériencemenée dans le basin à vagues du Laboratoire Hydraulique de France(ARTELIA). La formation de courants sagittaux a été forcée parun front de vagues avec un déficit d’énergie au centre du bassin. Lemodèle numérique est validé par des mesures de surface libre, devitesse, ainsi que de circulation et vorticité. En utilisant ensuitel’équation de vorticité potentielle comme outil de diagnostic, avec unforçage monochromatique on prédit un équilibre entre la générationde vorticité et son advection par l’écoulement moyen. / This work is about the hydrodynamic processes in the nearshorezone. They are of great importance to estimate the overall dynamicsof the coastal zone. This thesis is divided into two main parts; thefirst one investigates the coastal bottom boundary layer induced bythe interaction of the waves and the bottom when approaching thecoast; the second one is about the evolution of the mean circulationand vorticity induced by an inhomogeneity in the bathymetry orthe wave forcing. A turbulent boundary layer numerical model hasbeen developed and used to simulate the evolution of the oscillatingboundary layers under non-linear waves, of a flume experiment at theLaboratoire des Ecoulements Géophysiques et Industriels (LEGI) inGrenoble, France. The experimental instantaneous velocity profilesand still bed positions, allow defining the non-linear velocity distributionsinduced by the waves within the boundary layer. The numericalmodel coupled with a ad-hoc modeling of the mobile bed motionis able to reproduce the vertical distribution of the non-linearities,and also indicates that the vertical diffusion observed experimentallyis mainly caused by the mobile bed motion induced by the passingwaves. A 2D depth-averaged nonlinear shallow water numericalmodel is used to study the circulation and vorticity in the nearshorezone. This model is validated on a mobile bed experiment in thewave basin of the Laboratoire Hydraulique de France (ARTELIA).The formation of rip currents is forced by a damped wave forcing inthe middle of the wave basin. The numerical model is validated withfree surface and velocity measurements, and by the circulation andvorticity. Using the potential vorticity balance as a diagnosis tooland with a monochromatic wave forcing, an equilibrium between thevorticity generation and advection is observed in the nearshore zone. / Este trabajo trata de los procesos hidrodinámicos en la zona litoral,de grande importancia para la dinámica global del flujo costero. Dostemas principales son estudiados. El primero trata de la capa límiteoscilante provocada por la interacción entre el oleaje y el fondo alacercarse a la costa. El segundo tema trata de la evolución de lacirculación y la vorticidad inducida por la batimetría y/o el oleaje.Un modelo de capa límite turbulenta ha sido elaborado y validadopara analizar la evolución de la capa límite oscilante bajo la influenciade oleaje no-lineal, apoyándose en una modelación física, realizada enel canal de olas del LEGI. Los perfiles experimentales instantáneos develocidad y posición del fondo fijo, permiten definir la evolución delas no-linealidades inducidas por las olas dentro de la capa límite. Elmodelo numérico acoplado a una modelación del movimiento del fondomóvil es capaz de reproducir la evolución de estas no-linealidades, yexplica también que la difusión vertical observada experimentalmentees principalmente debida al movimiento vertical del fondo inducidopor el oleaje. El estudio de la circulación y de la vorticidad en zonascosteras se hace mediante un modelo numérico 2D promediado enla vertical de tipo Shallow Water que es validado con los datos deuna experiencia llevada a cabo en la piscina de olas del LaboratoireHydraulique de France (ARTELIA). La formación de corrientes ripse realiza a través de frentes de olas con un déficit de energía en elmedio de la piscina. El modelo numérico es validado con medicionesde superficie libre, de velocidades, y de circulación y vorticidad.Utilizando la ecuación de vortcidad potencial como herramienta dediagnóstico, con un oleaje monocromático se predice un equilibrioentre la generación de vorticidad y su advección por las corrientes.