Academic literature on the topic 'Écoulements surface libre'

Une simulation numérique par des éléments finis des écoulements transitoires à surface libre dans les canaux prismatiques est présentée. Dans cette étude, l'écoulement est supposé unidirectionnel dans un canal de faible pente. Le modèle mathématique est constitué d'un système de deux équations aux dérivées partielles de type hyperbolique résolu numériquement par la méthode des éléments finis. Pour définir les fonctions d'interpolation dans la forme intégrale des résidus pondérés, la méthode de Galerkin a été utilisée. Dans les applications, différentes sections prismatiques sont examinées. Les régimes transitoires étudiés sont dus à des manœuvres de vanne placée en aval du canal, l'extrémité amont étant connectée à un réservoir de niveau constant. Dans ces conditions, le régime transitoire correspond à une évolution de l'écoulement d'un régime permanent initial vers un régime permanent final. Ces deux régimes sont supposés uniformes à débit constant défini par la formule de Manning. Les résultats obtenus concernent l'évolution des paramètres hydrauliques en différentes sections du canal, suite à la manœuvre en aval. Deux cas de manœuvres sont considérés ; le cas d'une ouverture et le cas d'une fermeture. L'étude a permis d'analyser la propagation des ondes de surface et la réflexion de ces ondes sur les deux extrémités du canal. En particulier, les résultats numériques montrent que lorsque la largeur du lit du canal est très petite (cas de la section triangulaire), les fluctuations des profondeurs sont rapidement amorties.

La méthode SPH (Smoothed Particles Hydrodynamics) est une méthode purement Lagrangienne qui a été utilisée depuis 1977 dans la simulation des problèmes astrophysiques et les chocs à grandes vitesses. Dans le présent travail, une adaptation de cette méthode aux écoulements à surfaces libres est présentée. Comme première étape, le problème de bris de barrage en mono-dimensionnel est simulé. Cette étape nous a servi à évaluer et se familiariser avec la méthode SPH. Les résultats obtenus ont été validés et comparés avec la solution analytique ainsi qu'avec les résultats déjà publiés dans la littérature. La deuxième partie de ce mémoire a été consacrée aux cas bidimensionnels. Le problème de bris de barrage dans diverses configurations a été simulé. Les résultats obtenus dépendent directement de chaque cas traité. Pour le cas d'un barrage standard ainsi que le barrage circulaire, les résultats sont très encourageants. En ce qui concerne le cas d'un barrage avec étranglement, le résultat est oscillant. Ces résultats ont été obtenus avec l'approche standard de la méthode SPH. De plus, une nouvelle forumlation basée sur une combinaison EF-SPH a été introduite et testée. Les résultats obtenus sont meilleurs que celle donnés par la formulation standard. Cette amélioration n'est palpable que dans les premiers instants de la simulation. Durant les derniers instants des simulations, un manque de particules à l'emplacement du barrage, cause une concordance moins bonne des courbes trouvées avec celles fournies par la méthode des volumes finis. Le succès de la méthode SPH dans sa version standard ainsi que dans sa nouvelle version combinée EF-SPH, dépendent de plusieurs facteurs. Des études supplémentaires sont par conséquent suggérées pour standardiser les paramètres qui interviennent dans ces formulations.

La méthode des éléments finis (M. E. F. ) a été appliquée avec succès pour représenter l'hydrodynamique à surface libre, mais avec une modélisation de la turbulence très simplifiée. Ce travail a pour but de simuler, à l'aide de la M. E. F. , des écoulements à surface libre avec des représentations de la turbulence plus complexes. Il commence par une présentation générale des modèles de turbulence. Nous avons choisi d'utiliser plus particulièrement deux types de modèle : - le modèle k-e, - des modèles aux longueurs de mélange développés à partir de la simplification du modèle précédent. Une simplification supplémentaire est effectuée en bidimensionnalisant les modèles. Un travail de ce type a déjà été effectué par RODI et al pour le modèle k-e. Nous l'avons complété pour les modèles aux longueurs de mélange. Dans le cas où les courbures de l'écoulement horizontal sont importantes, nous proposons des variantes des deux types de modèle précédents et par ailleurs, une modification des frottements de fond intervenant dans l'équation de l'écoulement (SAINT-VENANT). Cette dernière est notamment importante pour les études sédimentologiques. La méthode de résolution des équations du modèle k-e utilisée dans cette thèse est la méthode des éléments finis ; il s'agit alors de résoudre un système non linéaire complexe ; dans ce travail, nous avons fait un choix d'interpolation particulier et testé plusieurs méthodes de résolution. Le couplage avec le modèle hydrodynamique est effectué par un algorithme explicite. Sont présentés plusieurs critères pour contrôler la convergence du modèle k-e. Nous terminons par une validation des modèles de turbulence utilisés sur de nombreux exemples. Sur ces derniers, nous avons vérifié les méthodes employées pour résoudre les équations du modèle k-e.

Dans cette thèse, nous nous intéressons à la résolution numérique des équations tridimensionnelles de Navier–Stokes à surface libre, écrites sous formalisme ALE, par la méthode des éléments finis. Ces équations permettent la représentation d’écoulements environnementaux. L’objectif initial et principal porte sur une analyse précise des limitations existantes des méthodes numériques et de possibles voies d’amélioration, mathématiquement justifiées, afin de réviser et d’améliorer le cœur numérique du code hydrodynamique Telemac-3D, développé par le Laboratoire National d’Hydraulique et Environnement d’EDF R&D. Nous décrivons ainsi précisément l’algorithme de ce dernier et nous l’évaluons à l’aune des publications les plus récentes. Ainsi, au sein de ce logiciel,les équations de Navier–Stokes à surface libre dans lesquelles la pression a été décomposée en une partie hydrostatique et une partie non-hydrostatique sont résolues. Une des limitations majeures concerne le fait que la vitesse calculée ne satisfait pas l’équation de continuité. De plus, nous avons relevé l’existence d’une restriction sur le pas de temps.Des voies alternatives sont explorées, et un algorithme très différent est analysé et programmé à des fins de comparaison. La stratégie numérique consiste à convecter la surface libre, à mettre à jour le domaine et à résoudre ensuite les équations de Navier–Stokes.Considérant cette approche, nous analysons un schéma éléments finis explicite d’ordre 1 en temps avec un terme de stabilisation symétrique pour l’équation cinématique de surface libre. Nous montrons que ce schéma permet de satisfaire deux propriétés importantes: la conservation globale de la masse d’eau et la divergence nulle faible. De plus, nous démontrons la stabilité conditionnelle de l’algorithme et nous proposons une nouvelle formulation variationnelle associée aux équations de Navier–Stokes permettant d’aboutir à un schéma éléments finis semi-implicite en temps, qui est stable indépendamment de la vitesse du domaine et d’une vitesse de convection à divergence nulle exacte.Finalement, certains modèles simplifiés utilisés dans d’autres logiciels développés par EDF R&D, notamment le modèle asymptotique en fréquence de Berkhoff utilisé au sein du logiciel Artemis [51], sont également étudiés afin d’acquérir une vue d’ensemble du domaine.L’équation de Berkhoff est dérivée asymptotiquement des équations de la houle linéaire.Une étude systématique de la dérivation de différents modèles simplifiés (Berkhoff, Shallowwater linéaire, Helmholtz) à partir des équations de la houle linéaire est effectuée.Le calcul d’une solution approchée des équations de Berkhoff utilisant des résultats de l’optique géométrique complète ce chapitre
In this thesis, we are interested by Finite Elements methods for the three-dimensional free surface Navier–Stokes equations under the ALE formulation. They enable to simulate geophysical flows. The initial and main goal is to analyse the existing limitations of these numerical methods and to provide perspectives of improvement, justified mathematically.This purpose helps us to present a review and improvement way for Telemac-3D, which is a hydrodynamics industrial software developed by the Laboratoire National d’Hydraulique et Environnement of EDF R&D. Therefore, we analyse precisely and we evaluate this algorithm, with respect to the recent scientific publications. This software solves the free surface Navier–Stokes equations with the decomposition of the pressure through a hydrostatic part and a dynamic part. A major limitation is that the velocity field of the fluid is not divergence-free. Furthermore, we highlight a time restriction on the time step. Moreover, alternative approaches are studied and compared. In particular, we focus on a numerical strategy which consists in advecting the free surface, in updating the domain and in solving the Navier–Stokes equations. Based on this strategy, we analyze a first order explicit scheme in time with a Finite Elements stabilization term. The numerical method allows to ensure important properties : the mass conservation of the water quantity and the weak free divergence condition. We demonstrate that this scheme is conditionally stable in time. Besides, we propose a new variational formulation allowing to obtain a semi-implicit scheme in time combined with the Finite Elements method,which is stable independently from the velocity of the mesh and without an exact free divergence velocity.Finally, in order to expand the hydrodynamic knowledges, some simplified models used in other software developed by EDF R&D are studied. In particular, we focus on the mild-slope equation solved in the software Artemis [51]. It is an asymptotic model derived from the linear water wave equation. As a consequence, we study the hypothesis and the validity of the derivation. An approximate analytical solution is additionally derived for this purpose. Moreover, comparisons with other asymptotic models, such as the linear shallow water equation or the Helmholtz equation, are presented

Cette thèse porte sur l’étude des ondes et instabilités dans les écoulements tournants à surface libre. On s’intéresse en particulier à la situation générique d’un récipient cylindrique partiellement rempli d’un fluide incompressible et dont le fond est en rotation rapide. En particulier, la formation de motifs connus sous le nom de polygones tournants ainsi que des phénomènes d’alternance temporelle nommés switching et sloshing est considérée. Dans ce but, des analyses de stabilité globale sont présentées en considérant des champs d’écoulement de base modèles tels que la rotation solide (seau de Newton), la rotation potentielle et le tourbillon de Rankine. Egalement, quelques résultats expérimentaux originaux obtenus lors d’un séjour au DTU (Copenhague) sont présentés. Les analyses de stabilité globales permettent de montrer que diverses familles d’ondes sont présentes dans ces configurations, à savoir des ondes de surface gravitaires et centrifuges, des ondes inertielles, des ondes de Rossby et des ondes de Kelvin-Kirchhoff. Il est constaté que le seau de Newton est toujours stable dans la gamme des paramètres considérés. Au contraire, des instabilités sont obtenues pour des nombres d’onde azimutaux m>=2 pour le tourbillon potentiel. Le diagramme de stabilité obtenu dans l’espace des paramètres (a, F), avec a le rapport d’aspect et F un nombre de Froude, reproduit qualitativement les zones d’observation expérimentales des motifs polygonaux. Les instabilités sont dues à une résonance entre des ondes gravitaires et des ondes centrifuges. Le mécanisme peut être interprété en terme d’onde à énergie négative à l’aide d’un modèle simplifié à deux couches, ou en terme d’over-reflection en utilisant une approche de type WKBJ pour de grandes valeurs de m et dans l’approximation de faibles profondeurs. Par ailleurs, les mécanismes d’un phénomène d’hysteresis observé expérimentalement ainsi que le phénomène de switching sont appréhendés à l’aide d’une étude faiblement non-linéaire effectuée sur le modèle à deux couches. Enfin, le modèle de Rankine, plus représentatif de l’expérience pour de faibles nombres de Froude est étudié et révèle la présence de nouvelles interactions d’ondes conduisants à des instabilités. Ces nouvelles instabilités capturent le phénomène de sloshing et donnent des éléments additionnels pour la comparaison avec les expériences.

Les cavités latérales sont des zones mortes à surface libre situées sur le côté d’un écoulement fluvial ou côtier. Les vitesses caractéristiques au sein de la cavité étant beaucoup plus faibles que celles de l’écoulement, une couche de mélange se développe à l’interface entre ces deux régions. Cette couche de mélange peut alors transférer de la quantité de mouvement de l’écoulement vers la cavité et ainsi mettre en mouvement la cavité et peut aussi transférer de la masse entre les deux régions, telle une pollution venant de l’écoulement amont. L’étude de cette thèse a alors consisté à étudier les caractéristiques de la couche de mélange, qui est rendue spécifique par le fait qu’elle se développe entre deux coins géométriques formés par l’intersection entre les parois de la cavité et celles de l’écoulement principal. Nous avons alors pu identifier l’origine et l’alternance des mouvements de fluide dans la direction transverse: de la cavité vers l’écoulement et inversement. Concernant la mise en mouvement de la cavité, le choix a été fait de considérer un écoulement principal fixé et de modifier l’extension de la cavité dans la direction perpendiculaire à l’écoulement, passant ainsi d’une cavité rectangulaire alignée avec l’écoulement principal à une cavité allongée dans le sens opposé. La mesure de champ de vitesse par PIV 2D a alors montré une forte évolution de la forme de l’écoulement à mesure que la géométrie de la cavité évolue : un système avec deux cellules alignées dans le sens de l’écoulement à un système à une seule cellule, puis un système à deux cellules et enfin un système complexe 3D ont ainsi été observés pour une cavité de plus en plus allongée. Ensuite, une modification du dispositif expérimental a permis de mesurer de deux façons différentes le transport de scalaire de l’écoulement principal vers la cavité, de comprendre les processus associés à ce transfert et enfin de quantifier cette capacité de transfert pour différents écoulements principaux et différentes géométries de cavités. Nous avons notamment montré que la géométrie de la cavité a peu d’effet alors que le nombre de Reynolds et la profondeur d’eau normalisée ont un effet majeur sur cette capacité de transfert de masse entre les deux régions
Lateral cavities are free-surface dead-zones located on the side of a fluvial or coastal main flow. As the typical velocities are much larger in the main flow than in the cavity, a mixing layer appears at the interface between both regions. This mixing layer is able to transfer between the main flow and the cavity momentum which then sets the fluid in the cavity in motion and also passive scalar, such as a pollution coming from upstream. The objective of this work was then to investigate the characteristics of the mixing layer, which specificity comes from the fact that it is constrained between the upstream and downstream geometrical corners. It was possible to observe the origin and alternation of the transversal fluid motions: from the cavity towards the main flow and conversely. Regarding the motion in the cavity, the choice was made to keep a constant main flow and to measure the 2D horizontal velocity field using PIV as the extension of the cavity increases. The flow pattern then passes from a 2-cell patterns aligned in the direction of the main flow to a single-cell pattern, then a 2-cells patterns aligned along the direction perpendicular to the main flow and finally a complex 3D pattern for the widest cavity. Then a modification of the experimental set-up permitted to investigate the passive scalar exchanges from the main stream towards the cavity. It was possible to understand the processes responsible for such transfer and to quantify the transfer capacity. The analysis dimensional revealed that in the present subcritical, smooth simplified geometry cavity, the three parameters possible responsible for the modification of the transfer capacity are the geometrical aspect ratio of the cavity, the Reynolds number of the main flow and finally the normalized water depth. It was then shown that the impact of the cavity geometry remains negligible but that the Reynolds number and the normalized water do impact this passive scalar transfer capacity

L’objectif de cette thèse est de développer un modèle numérique 3D afin d’étudier le phénomène de laves torrentielles ; écoulements visqueux fortement chargés en matière solide, surgissant en montagne lors d’orages violents. Aujourd’hui, la prévision des zones de vulnérabilité s’appuie sur des outils de calcul 0D, 1D ou 2D. Or ces outils ne peuvent représenter pleinement le comportement à surface libre des écoulements du fait de nombreuses approximations et hypothèses. Ainsi cette thèse met en œuvre un code numérique 3D pour étudier ces écoulements. Ce travail aboutit à la construction d’un modèle 3D à partir de données réelles de terrain. Plusieurs scénarios ont été étudiés et comparés à des résultats issus d’un modèle 2D. Les résultats mettent en évidence les apports non négligeables de la modélisation 3D : zones d’étalement et de dépôt, phénomènes d’encombrement, modélisation fine des écoulements dans les zones chenalisées. Préalablement, le modèle 3D a été validé en comparant les résultats numériques à des données expérimentales issues de la littérature, pour des typologies d’écoulement représentatives de celles observées sur des sites grandeur nature
The objective of this thesis is to develop a 3D numerical model to assess debris flow. These viscous flows, heavily loaded with solid matter, form when heavy rain occurs in mountains. Today, forecasts of potentially impacted areas are based on 0D, 1D and 2D numerical tools. However, these tools cannot fully represent the free surface behaviour of debris flows due to the approximations and assumptions on which they are based. Thus, this work utilises a 3D numerical code to study this phenomenon. A specific model is built with real field data. Several flow scenarios are studied and compared with a 2D numerical model. The results highlight the significant benefits of a 3D approach by providing information on the fine representation of flow dynamics over the catchment area. The model also predicts the impact of debris flow (overflowing on a road bridge) and the zones of deposition and spreading. It highlights possible congestion phenomena and reproduces flows in the channels by fully accounting for parietal friction, capabilities not provided by 2D models. Prior to this application, the 3D model was evaluated with five sets of experimental data to validate its ability to represent viscoplastic flows. Different types of flows are studied and are representative of those observed on real sites when debris flow occur

Cette thèse a pour objectif la construction de schémas d’ordre élevé et peu diffusifs pour le transport d’un scalaire dans les écoulements à surface libre, en deux ou trois dimensions. On souhaite en particulier obtenir des schémas robustes, qui gardent au niveau discret les propriétés mathématiques de l’équation de transport avec une faible diffusion numérique, et les utiliser sur des cas industriels. Dans ce travail deux méthodes numériques sont envisagées : une méthode aux volumes finis (VF) et une méthode aux résidus distribués (RD). Dans les deux cas, l’équation de transport est résolue avec une approche découplée, qui est la solution la plus avantageuse en termes de précision et de coûts de calcul. Pour ce qui concerne la méthode aux volumes finis, les équations de Saint-Venant couplées à l’équation du transport sont d’abord résolues avec un schéma dit vertex-centred où le flux numérique est approximé avec un solveur de Riemann appelé Harten-Lax-Van Leer-Contact [135]. A partir de cette approche, une formulation découplée est proposée. Cette dernière permet de résoudre l’équation du transport avec un pas de temps plus grand que celui de la formulation couplée. Cette idée a été d’abord proposée pour d’autres schémas dans [13]. Pour augmenter l’ordre de précision en espace, la technique MUSCL [89] est utilisée en combinaison avec l’approche découplée. Finalement, la problématique des zones sèches est abordée. Dans le cas de la méthode aux résidus distribués, les équations de Saint-Venant sont résolues avec une méthode éléments finis, et la méthode RD est utilisée seulement pour discrétiser l’équation du transport, en focalisant l’attention sur les problèmes non stationnaires. L’équation de continuité du fluide discrétisée est employée pour garantir la conservation de la masse et le principe du maximum. Pour obtenir des schémas d’ordre deux dans les problèmes non stationnaires, un schéma prédicteur-correcteur [112] est utilisé, en l’adaptant au cas de concentration moyennée sur la verticale. Une version d’ordre 1 mais peu diffusive, est aussi présentée dans ce travail. De plus, un schéma localement implicite, complètement nouveau, est aussi formulé pour pouvoir traiter le problème des bancs découvrant. Les deux techniques sont validées d’abord sur des cas simples, pour évaluer l’ordre de précision des schémas et ensuite sur des cas plus complexes pour vérifier aussi les autres propriétés numériques. Les résultats montrent que les nouveaux schémas sont à la fois précis et conservatifs, tout en gardant la monotonie comme le prévoient les démonstrations. Un cas d’application industriel est aussi présenté en conclusion. Le schéma prédicteur-correcteur RD est adapté aussi au cas 3D, sans aucun problème théorique nouveau, par rapport au cas 2D. Les propriétés de base des schémas sont validées sur des cas test préliminaires
The purpose of this thesis is to build higher order and less diffusive schemes for pollutant transport in shallow water flows or 3D free surface flows. We want robust schemes which respect the main mathematical properties of the advection equation with relatively low numerical diffusion and apply them to environmental industrial applications. Two techniques are tested in this work: a classical finite volume method and a residual distribution technique combined with a finite element method. For both methods we propose a decoupled approach since it is the most advantageous in terms of accuracy and CPU time. Concerning the first technique, a vertex-centred finite volume method is used to solve the augmented shallow water system where the numerical flux is computed through an Harten-Lax-Van Leer-Contact Riemannsolver [135]. Starting from this solution, a decoupled approach is formulated and is preferred since it allows to compute with a larger time step the advection of a tracer. This idea was inspired by [13]. The Monotonic Upwind Scheme for Conservation Law [89], combined with the decoupled approach, is then used for the second order extension in space. The wetting and drying problem is also analysed and a possible solution is presented. In the second case, the shallow water system is entirely solved using the finite element technique and the residual distribution method is applied to the solution of the tracer equation, focusing on the case of time-dependent problems. However, for consistency reasons the resolution of the continuity equation must be considered in the numerical discretization of the tracer. In order to get second order schemes for unsteady cases a predictor-corrector scheme [112] is used in this work. A first order but less diffusive version of the predictor-corrector scheme is also introduced. Moreover, we also present a new locally semi-implicit version of the residual distribution method which, in addition to good properties in terms of accuracy and stability, has the advantage to cope with dry zones. The two methods are first validated on academical test cases with analytical solution in order to assess the order of the schemes. Then more complex cases are addressed to test the robustness of the schemes and their performance under different flow conditions. Finally a real test case for which real data are available is carried out. An extension of the predictor-corrector residual distribution schemes to the 3D case is presented as final contribution. Even in this case the RD technique is completely compatible with the finite element framework used for the Navier-Stokes equations, thus its extension to the 3D case does not present any extra theoretical problem. The method is tested on preliminary cases

L’objectif de ce travail est l’interpolation du débit à partir de données issues de vélocimètre à effet Doppler. Pour ce faire une modélisation précise des écoulements turbulents doit être implémentée dans un code conçu pour les conduites d’assainissement. Du fait de l’uniformité de l’écoulement dans le sens du flux, un maillage 2D est utilisé pour effectuer le calcul du champ de vitesse en 3D. Un modèle RSM des écoulements turbulents uniformes à surface libre a donc été développé en version bas-Reynolds afin de pouvoir intégrer les équations en zone proche paroi incluant la sous-couche visqueuse. Une combinaison de (Shir, 1973), (Gibson & Launder, 1978) et (Cokljat, 1993) a été additionnée au terme de corrélation fluctuation de pression/champ moyen de vitesse afin de prendre en compte l’effet d’amortissement de la surface libre et le caractère redistributif d’énergie cinétique turbulente des parois. Ce modèle a été implémenté (sous Matlab) dans un code aux volumes finis, tout comme un modèle k-ε à titre de comparaison. On montre ainsi le rôle prépondérant du caractère anisotrope de la turbulence dans la génération des courants secondaires. Les résultats présentés sont comparés avec ceux de (Kang & Choi, 2005) qui ont développé un modèle RSM utilisant des lois de paroi et avec les résultats expérimentaux de (Nezu & Rodi, 1985). Il y a une bonne corrélation entre les résultats du modèle et les mesures. Les effets sur les tensions de Reynolds de l’annulation de la contrainte tangentielle aux coins sont reproduits. Par ailleurs, la projection du champ moyen de vitesse sur une section transverse met en évidence des courants secondaires, notamment des petites structures au niveau des angles formés par la surface libre avec les parois. Enfin, des cas expérimentaux issus de la littérature sont modélisés afin de comparer les débits mesurés avec ceux calculés par le modèle présenté, là encore il y a une bonne corrélation. Pour finir, un banc d’essai a été conçu pour la mesure de vitesse en canal. Un banc de déplacement automatisé a été assemblé autour d’un canal de 20m de long, il est équipé des systèmes de mesure par ADV et PIV
The aim of this work is the interpolation of a mean flow stream from discreet data supplied by a Doppler velocimeter. This is why a precise flow modeling needs to be implemented in a code designed for sewer net conduit. A 2D mesh is used to calculate 3D velocity flow. This work deals with the description of a Reynolds stress model for the numerical modeling of uniform 3D turbulent open-channel flows. The turbulence model is low-Reynolds one in order to enable the viscosity-affected region to be resolved including the viscous sublayer. Therefore, in order to consider the damping effects at the free surface and the redistribution of turbulent kinetic energy effects at the wall, a combinations of (Shir, 1973), (Gibson & Launder, 1978) and (Cokljat, 1993) models are added to the pressure-strain term. The finite volume method is used for the numerical solution (with Matlab) of the flow equations and transport equations of the Reynolds stress components. K-ε model has been implanted to compare with the developed model; it is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. The results presented here are compared with results predicted by the RSM model of (Kang & Choi, 2005) employing wall functions and with experimental data of (Nezu & Rodi, 1985). There is a good agreement between the measurements and the results predicted. Effects of low corner’s wall shear stress on Reynolds stress are underlined. The computed flow vectors in the plane normal to the main flow direction show secondary currents and moreover a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface. Then experimental cases are modeled, good agreement of flow rate predicted with measurement was found. Finally test bench has been created for velocity channel measurements. An automated portal structure has been assembled all around a 20m long channel to fix onto the ADV and PIV sensors

Depuis une vingtaine d'années se succèdent des inondations dans nos pays industrialisés causant de graves pertes humaines et des dégâts matériels consi¬dérables. Cependant le recours à des outils de simulation des écoulements à l'intérieur des rues, apporterait une meilleure connaissance des processus de submersion et de propagation des ondes de crues. C'est ainsi que des solutions d'aménagement pourront être et faire l'objet de comparaison. Ces variantes permettront de mettre en évidence les contraintes et de proposer des solutions visant à diminuer les risques. Pour cela nous avons conçu un modèle numérique 1D/2D pour les écoulements de crues dans les tissus urbains qui prend en considération les éléments suivants : le choix du schéma d'intégration numérique afin d'assurer la stabilité et la précision du modèle, la capacité de capturer les chocs notamment lors du passage d'un régime torrentiel à un régime fluvial avec apparition de ressauts hydrauliques, la capture du front raide des écoulements sur fond sec en présence ou non des précipitations avec conservation de débit et enfin la capacité de traiter les écoulements dans les réseaux à jonctions tout en assurant une distribution correcte du débit dans les branches associées.