Analytical solutions for the linear and nonlinear shallow-water wave equations are developed for evolution and runup of tsunamis &ndash / long waves&ndash / over one- and two-dimensional bathymetries. In one-dimensional case, the nonlinear equations are solved for a plane beach using the hodograph transformation with eigenfunction expansion or integral transform methods under different initial conditions, i.e., earthquake-generated waves, wind set-down relaxation, and landslide-generated waves. In two-dimensional case, the linear shallow-water wave equation is solved for a flat ocean bottom for initial waves having finite-crest length. Analytical verification of source focusing is presented. The role of focusing in unexpectedly high tsunami runup observations for the 17 July 1998 Papua New Guinea and 17 July 2006 Java Island, Indonesia tsunamis are investigated. Analytical models developed here can serve as benchmark solutions for numerical studies.

Reproducing moisture retention behavior of the upper and lower vadose zone in shallow water table settings provides unique challenges for integrated (combined surface and groundwater) hydrological models. Field studies indicate that moisture retention in shallow water table settings is highly variably affected by antecedent state and air entrapment. The theory and vertical behavior of a recently developed integrated surface and groundwater model (IHM) is examined through comparisons to collected field data in West-Central Florida. The objectives of this study were to (1) Identify important considerations and behavior of the vadose zone for reproducing runoff, ET and recharge in shallow water table settings; (2) Develop a conceptual model that describes vertical soil moisture behavior while allowing for field scale variability; (3) Test the model against observations of the vertical processes; (4) Investigate the sensitivity of model parameters on model vs. observed vertical behavior, and (5) offer recommendations for improvements and parameterization for regional model application. Rigorous testing was made to better understand the robustness and/or limitations of the methodology of the IHM for upper and lower vadose zone. The results are also generally applicable and useful to the upper zone and lower zone conceptualization and parameterization of stand alone HSPF and perhaps other surface water models. Simulation results indicate IHM is capable of providing reasonable predictions of infiltration, depth to water table response, ET distributions from the upper soil, lower soil and water table, and recharge while incorporating field scale variability of soil and land cover properties.

The two-dimensional depth-averaged shallow water equations have attracted considerable attentions as a practical way to solve flows with free surface. Compared to three-dimensional Navier-Stokes equations, the shallow water equations give essentially the same results at much lower cost. Solving the shallow water equations by the Godunov-type finite volume method is a newly emerging area. The Godunov-type finite volume method is good at capturing the discontinuous fronts in numerical solutions. This makes the method suitable for solving the system of shallow water equations. In this dissertation, both the shallow water equations and the Godunov-type finite volume method are described in detail. A new surface flow routing method is proposed in the dissertation. The method does not limit the shallow water equations to open channels but extends the shallow water equations to the whole domain. Results show that the new routing method is a promising method for prediction of watershed runoff. The method is also applied to turbulence modeling of free surface flow. The κ - ε turbulence model is incorporated into the system of shallow water equations. The outcomes prove that the turbulence modeling is necessary for calculation of free surface flow. At last part of the dissertation, a total load sediment transport model is described and the model is tested against 1D and 2D laboratory experiments. In summary, the proposed numerical method shows good potential in solving free surface flow problems. And future development will be focusing on river meandering simulation, non-equilibrium sediment transport and surface flow - subsurface flow interaction.

The physical processes in lakes remain only partially understood despite successful data collection from a variety of sources spanning several decades. Although numerical models are already frequently employed to simulate the physics of lakes, especially in the context of water quality management, improved methods are necessary to better capture the wide array of dynamically important physical processes, spanning length scales from ~ 10 km (basin-scale oscillations) - 1 m (short internal waves). In this thesis, high-order numerical methods are explored for specialized model equations of lakes, so that their use can be taken into consideration in the next generation of more sophisticated models that will better capture important small scale features than their present day counterparts. The full three-dimensional incompressible density-stratified Navier-Stokes equations remain too computationally expensive to be solved for situations that involve both complicated geometries and require resolution of features at length-scales spanning four orders of magnitude. The main source of computational expense lay with the requirement of having to solve a three-dimensional Poisson equation for pressure at every time-step. Simplified model equations are thus the only way that numerical lake modelling can be carried out at present time, and progress can be made by seeking intelligent parameterizations as a means of capturing more physics within the framework of such simplified equation sets. In this thesis, we employ the long-accepted practice of sub-dividing the lake into vertical layers of different constant densities as an approximation to continuous vertical stratification. We build on this approach by including weakly non-hydrostatic dispersive correction terms in the model equations in order to parameterize the effects of small vertical accelerations that are often disregarded by operational models. Favouring the inclusion of weakly non-hydrostatic effects over the more popular hydrostatic approximation allows these models to capture the emergence of small-scale internal wave phenomena, such as internal solitary waves and undular bores, that are missed by purely hydrostatic models. The Fourier and Chebyshev pseudospectral methods are employed for these weakly non-hydrostatic layered models in simple idealized lake geometries, e.g., doubly periodic domains, periodic channels, and annular domains, for a set of test problems relevant to lake dynamics since they offer excellent resolution characteristics at minimal memory costs. This feature makes them an excellent benchmark to compare other methods against. The Discontinuous Galerkin Finite Element Method (DG-FEM) is then explored as a mid- to high-order method that can be used in arbitrary lake geometries. The DG-FEM can be interpreted as a domain-decomposition extension of a polynomial pseudospectral method and shares many of the same attractive features, such as fast convergence rates and the ability to resolve small-scale features with a relatively low number of grid points when compared to a low-order method. The DG-FEM is further complemented by certain desirable attributes it shares with the finite volume method, such as the freedom to specify upwind-biased numerical flux functions for advection-dominated flows, the flexibility to deal with complicated geometries, and the notion that each element (or cell) can be regarded as a control volume for conserved fluid quantities. Practical implementation details of the numerical methods used in this thesis are discussed, and the various modelling and methodology choices that have been made in the course of this work are justified as the difficulties that these choices address are revealed to the reader. Theoretical calculations are intermittently carried out throughout the thesis to help improve intuition in situations where numerical methods alone fall short of giving complete explanations of the physical processes under consideration. The utility of the DG-FEM method beyond purely hyperbolic systems is also a recurring theme in this thesis. The DG-FEM method is applied to dispersive shallow water type systems as well as incompressible flow situations. Furthermore, it is employed for eigenvalue problems where orthogonal bases must be constructed from the eigenspaces of elliptic operators. The technique is applied to the problem calculating the free modes of oscillation in rotating basins with irregular geometries where the corresponding linear operator is not self-adjoint.

Esta dissertação tem por objetivo a implementação de um código para simular problemas hidrodinâmicos, bem como a possibilidade de controlar as elevações de onda resultantes numa determinada região por meio de uma vazão ótima controlada dentro do sistema estudado. O algoritmo implementado é baseado nas equações de águas rasas, as quais são aplicáveis em situações onde a altura d’água é de ordem muito menor do que as dimensões do sistema, que é discretizado espacial e temporalmente pelo Método dos Elementos Finitos e pelo método CBS (Characteristic Based-Split), respectivamente. O método de controle consiste na busca de uma curva de vazão de controle ótima que minimize a função objetivo, a qual compara os valores de altura d’água que se deseja encontrar em uma região especificada com os calculados pela simulação numérica. Para isso, utiliza-se um algoritmo evolutivo SCE-UA (Shuffled Complex Evolution – University of Arizona), que busca otimizar parâmetros de geração das curvas de vazão de controle, podendo estas serem modeladas por NURBS (Non- Uniform Rational B-Splines), que são capazes de encontrar a solução ótima, ou modeladas com curvas de forma triangular (linear) ou parabólica (quadrática) que apresentam uma solução aproximada de fácil implementação. Por fim, várias aplicações são realizadas, tanto para a simples simulação, quanto para o controle de problemas hidrodinâmicos, a fim de validar os algoritmos desenvolvidos e os resultados obtidos mostraram que os objetivos foram alcançados, encontrando uma forma eficiente de se fazer o controle de enchentes. / Implementation of a computational code for the numerical simulation of hydrodynamic problems as well as the ability to control the resulting wave elevations in a specific area, using an optimal flow controlled within the studied system are the aims of this work. The implemented algorithm is based on the shallow waters equations, which are applicable in situations where the water height is much smaller than the system dimensions, and are spatial and temporally discretized by the Finite Element Method and the CBS method (Caractheristic Based-Split), respectively. The control method consists in finding an optimal control flow curve that minimizes the objective function, which compares the objective value of water elevations in a specified region with those calculated by numerical simulation. An evolutionary algorithm called SCE-UA (Shuffled Complex Evolution - University of Arizona), which looks for optimize parameters of control flow curves generation, is used. These curves may be modeled by NURBS (Non-Uniform Rational B-Splines) which are able to find the optimal solution, or by curves of triangular (linear) or parabolic quadratic forms, which are an approximate solution easy to implement. Finally, several applications are performed for both simulation and control of hydrodynamic problems in order to validate the developed algorithms, and the results showed that the aims of this work were reached, finding an efficient way to control floods.

Dans de nombreuses applications aérospatiales, on peut trouver des films liquides cisaillés, c'est-à-dire une fine couche liquide qui ruisselle sur une paroi entrainée par le gaz. Par exemple, une couche de liquide peut se développer sur la voilure des avions, givrer et dégrader les performances. Des vagues peuvent se développer à l'interface liquide-gaz, et l'analyse correcte de ces instabilités devient très importante pour modéliser ce phénomène physique. En effet, la présence d'instabilités modifie les échanges liquide-gaz, notamment les transferts de masse et chaleur. Le but de cette thèse est de développer une technique permettant de coupler la phase gazeuse afin de reproduire les interactions à l'interface. La couche de liquide étant beaucoup plus mince que celle du gaz, une approche intégrale sur l'épaisseur est utilisée pour la modélisation. Enfin, deux cas d'un écoulement diphasique se développant dans une conduite confinée et dans une conduite plus large sont étudiés. les résultats sont ensuite comparés à des autres méthodes de référence, plus coûteuses en temps de calcul. / In many aerospace applications one can find liquid films sheared by a gas flow. In example, these liquid sheets can develop on aircraft wings, freeze and then destroy the aerodynamics performances. Waves can develop at the liquid-gas interface, and the correct analysis of such instabilities becomes very important to model this physical phenomenon. Indeed, instabilities mdify liquid-gas exchanges, such as mass and heat transfers. The aim of the present work consists in developing a technique to couple the liquid phase to the gas phase in order to reproduce the interactions at the interface. Since the liquid layer is much thinner then the gas, anintegral approach is used for modeling. Finally, two cases of a two-phase flow developing in a strictly confined channel and in a large channel are studied. Results are then compared to other reference methods which are more expensive in terms of computational cost.

Esta dissertação tem por objetivo a implementação de um código para simular problemas hidrodinâmicos, bem como a possibilidade de controlar as elevações de onda resultantes numa determinada região por meio de uma vazão ótima controlada dentro do sistema estudado. O algoritmo implementado é baseado nas equações de águas rasas, as quais são aplicáveis em situações onde a altura d’água é de ordem muito menor do que as dimensões do sistema, que é discretizado espacial e temporalmente pelo Método dos Elementos Finitos e pelo método CBS (Characteristic Based-Split), respectivamente. O método de controle consiste na busca de uma curva de vazão de controle ótima que minimize a função objetivo, a qual compara os valores de altura d’água que se deseja encontrar em uma região especificada com os calculados pela simulação numérica. Para isso, utiliza-se um algoritmo evolutivo SCE-UA (Shuffled Complex Evolution – University of Arizona), que busca otimizar parâmetros de geração das curvas de vazão de controle, podendo estas serem modeladas por NURBS (Non- Uniform Rational B-Splines), que são capazes de encontrar a solução ótima, ou modeladas com curvas de forma triangular (linear) ou parabólica (quadrática) que apresentam uma solução aproximada de fácil implementação. Por fim, várias aplicações são realizadas, tanto para a simples simulação, quanto para o controle de problemas hidrodinâmicos, a fim de validar os algoritmos desenvolvidos e os resultados obtidos mostraram que os objetivos foram alcançados, encontrando uma forma eficiente de se fazer o controle de enchentes. / Implementation of a computational code for the numerical simulation of hydrodynamic problems as well as the ability to control the resulting wave elevations in a specific area, using an optimal flow controlled within the studied system are the aims of this work. The implemented algorithm is based on the shallow waters equations, which are applicable in situations where the water height is much smaller than the system dimensions, and are spatial and temporally discretized by the Finite Element Method and the CBS method (Caractheristic Based-Split), respectively. The control method consists in finding an optimal control flow curve that minimizes the objective function, which compares the objective value of water elevations in a specified region with those calculated by numerical simulation. An evolutionary algorithm called SCE-UA (Shuffled Complex Evolution - University of Arizona), which looks for optimize parameters of control flow curves generation, is used. These curves may be modeled by NURBS (Non-Uniform Rational B-Splines) which are able to find the optimal solution, or by curves of triangular (linear) or parabolic quadratic forms, which are an approximate solution easy to implement. Finally, several applications are performed for both simulation and control of hydrodynamic problems in order to validate the developed algorithms, and the results showed that the aims of this work were reached, finding an efficient way to control floods.

The thesis deals with the numerical solution of partial differential equati- ons describing the flow of the so-called shallow water neglecting the flow in the vertical direction. These equations are of hyperbolical type of the first or- der with a reactive term representing the bottom topology. We discretize the resulting system of equations by the implicit space-time discontinuous Ga- lerkin method (STDGM). In the literature, the explicit techniques are used most of the time. The implicit approach is suitable especially for adaptive methods, because it allows the usage of different meshes for different time niveaus. In the thesis we derive the corresponding method and an adaptive algorithm. Finally, we present usage of the method in several examples. 1

Cette thèse est une partie du projet ONEMA pour le dimensionnement des passes à poissons et pour l’amélioration de la continuité écologique des cours d’eau. Ce travail s’est concentré sur les passes à poissons naturelles qui présentent des avantages de coût et paysager. Il s’agit d’un écoulement à forte pente autour des blocs (macro-rugosités) régulièrement repartis en quinconce avec des grands nombres de Froude. Les conditions hydrodynamiques sont alors très diverses, et peuvent être franchissables par un nombre élargi d’espèces de poisson. Ce mémoire présente les travaux réalisés à l’Institut de Mécanique des Fluides de Toulouse (IMFT). Afin d’étudier l’écoulement dans ces passes, on va mener des expériences sur des canaux réduits ainsi que des simulations numériques à l’aide du modèle Telemac 2D. L’objectif est de mieux connaître la structure de l’écoulement en fonction des conditions hydrauliques et géométriques comme le nombre, la forme et la taille des macro-rugosités. Plus particulièrement, la compréhension de l’interaction de phénomènes physiques généralement étudiés séparément, tels que le passage en régime torrentiel, l’interaction de sillage ou l’écoulement autour de macro-rugosités, a été recherchée. Des relations hauteur-débit ont été établies permettant une aider au dimensionnement des passes naturelles. Elles fournissent des critères de franchissement comme les vitesses maximales, la puissance dissipée ou la hauteur d’eau minimale. Pour atteindre une description plus locale de l’écoulement, des mesures de Vélocimétrie Acoustiques Doppler ont été conduites. Elles ont aussi permis de définir la plage de validité du modèle numérique 2D (Telemac). Ce modèle a alors était utilisé pour extrapoler les critères de franchissement pour des configurations non testées expérimentalement. Finalement, les connaissances sur l’écoulement ont été synthétisées pour définir des préconisations générales de dimensionnement. La précision des relations établies en laboratoire a pu aussi être vérifiée sur des passes réelles. L’hydrodynamique de ces passes est maintenant suffisamment connue pour savoir si un poisson peut remonter le courant et se reposer. Il restera à s’assurer que leur attractivité soit bonne et que des phénomènes liés aux échelles de longueurs de la turbulence ne présentent un obstacle au franchissement. / This thesis is supported by an ONEMA project for the design of fishways and improve ecological continuity of rivers. This work focused on nature-like fish passes that have cost and landscaped appearance advantages. There is a steep flow around the blocks (macro-roughness) regularly distributed in a staggered configuration with large Froude numbers. The hydrodynamic conditions are sufficiently different to be passable by an expanded number of fish species. This thesis presents the work carried out at the Institute of Fluid Mechanics of Toulouse (IMFT). To study the flow in these passes, experiments are conducted on physical models and numerical simulations using the Telemac 2D model. The goal is to better understand the flow structure based on hydraulic and geometric conditions such as the number, shape and size of macro-roughness. In particular, we considered the interaction of physical phenomena usually studied separately such that the passage in supercritical regime, the interaction of wake or flow around macro-roughness. The stage-discharge relationships were established to assist in the design of nature-like passes. They provide criteria useful for passability such as maximum speeds, power dissipation or minimal water height. To reach a local description of the flow measurements with an Acoustic Doppler Velocimeter were conducted. They also helped to define the valid range of the 2D model. This model was then used to extrapolate passability criteria for untested experimental configurations. Finally, the knowledge on flow has been synthesized to define general recommendations sizing. It has also been able to verify on real scale passes, the accuracy of the relations established in the laboratory. The hydrodynamics of these passes is now sufficiently described to know if a fish can swim upstream and rest. It will remain to ensure that their attractiveness is good and that phenomena related to the turbulence length scales do not present an obstacle to the fish passage.

Les espèces de poissons migrateurs (saumon atlantique, Salmo salar, en particulier) requièrent des conditions thermiques bien spécifiques. Ils sont très sensibles aux températures de l’eau et aux fortes variations estivales. Sur les trente dernières années, l’étude menée sur la Garonne (France) révèle une augmentation des températures estivales associée à un allongement de la durée des périodes chaudes. L’impact de cette modification du régime thermique sur la survie et la reproduction des espèces migratoires est également mis en évidence. Cette étude est menée sur un tronçon de Garonne, situé entre l’amont de Toulouse et l’amont de la retenue deMalause. Ce secteur est fortement touché par cette problématique avec en moyenne 2°C d’écart entre l’amont et l’aval et des températures supérieures à 25°C régulièrement atteintes. Le régime hydrologique de ce tronçon est fortement déficitaire (selon le SDAGE Adour-Garonne), la sensibilité au flux de surface est forte à cause de son lit large et peu profond, les pressions anthropiques sont importantes, ce sont autant de pistes dont l’impact sur le régime thermique est étudié. Une large collection de données hydrologiques et climatiques est exploitée afin de déterminer les processus en jeu dans l’évolution du régime thermique de ce tronçon de fleuve. Des études en tendances et corrélations et des modèles statistiques permettent de mettre en évidence d’une part la relation forte qui existe entre les températures de l’air et les températures de l’eau et d’autre part l’importance des faibles débits durant les périodes estivales. L’estimation des moyennes journalières de température de l’eau à Malause au moyen de modèles statistiques et déterministes donne de bons résultats pour les températures élevées ainsi que pour les franchissements de seuils liés aux conditions de migrations des amphilalins.Enfin un modèle numérique monodimensionnel de résolution de l’équation de transport thermique et des équations de St-Venant est développé. La physique du modèle tant au niveau hydraulique (prise en compte de fortes variabilités de pente, d’ouvrages, etc.) que thermique (apports latéraux, flux de surfaces, flux de conduction avec le lit) permet d’analyser l’évolution des différents flux qui participent au réchauffement du cours d’eau. Une évolution future à l’aide des sorties des modèles de l’IPCC est explorée et des méthodes éventuelles de restauration des conditions de températures favorables pour les espèces piscicoles sont analysées. / Fish species with strong thermal requirements (i.e. Atlantic salmon) are very sensitive to temperature evolution and particularly to large increases. An investigation conducted on the Garonne River (France) during the last three decades revealed global water warming along with an increase of the high temperature period duration. Large impact of this evolution on the survival and breeding of migrating fish species was also reported. Study was thus conducted on a specific reach of the Garonne River located between the immediate upstream of Toulouse and the upstream of the Malause dam. The issue of water temperature warming is particularly relevant on this reach, with an average increase of 2°C between upstream and downstream and temperatures above 25°C frequently reported. Potential causes are numerous: drastic low-flow regime (quoting SDAGE Adour-Garonne), impacts of surface fluxes that are important due to bed shape (wide and shallow), anthropogenic impacts, etc. Large amount of climatic and hydraulic data are used to make a clear determination of the processes involved in the thermal regime evolution of this reach. Trend and correlation analyses and use of statistical models indicate the strong relation between stream temperature and climate. Low flows also seem to be related to water temperatures during summer periods. Statistic and deterministic models give good results in estimating high daily mean water temperatures (RMSE ranging from 0.99°C to 1.22°C) and predicting water temperatures threshold crossings related to the migrating conditions of Atlantic salmon.Finally, a one-dimensional numerical model that solves both shallow water and thermal equations is developed. Both the formulation of the St-Venant equations (high variability in slope, gates …) and the phenomena taken into account in the water temperature model (lateral influx, surface fluxes, bed conduction …) allows studying the evolution of fluxes driving water temperature evolution. Future evolution of the water temperature at the 2050 horizon is also evaluated using IPCC models output and potential solutions to restore favorable stream temperatures conditions for fishes are analyzed.