De nos jours, la combustion représente encore un 90% de la production totale d'énergie au monde. La plupart des brûleurs de type industriel utilisent comme carburant des hydrocarbures en forme liquide. Cependant, un grand nombre d'études ont été dédiés aux flammes gazeuses et l'impact du spray liquide est encore loin d'être totalement compris. Le but de cet étude est l'amélioration de la modélisation des deux phénomènes principaux qui ont lieu entre l'atomisation du spray et la combustion, i.e. la dispersion des gouttes par la turbulence gazeuse et le procès d'évaporation dans le contexte de la Simulation Aux Grandes Echelles (SGE) des configurations complexes. Premièrement, l'approche Euler-Euler mésoscopique (Février et al. (2005)), basée sur une moyenne d'ensemble conditionnée et implémentée dans AVBP est amélioré. Le modèle de fermeture (Simonin et al. (2001); Kaufmann (2004)) pour les moments de deuxième ordre qui apparait dans les équations de transport résolues échoue quand appliqué à des configurations cisaillées (Riber (2007)). Plusieurs modèles proposés récemment par Masi (2010) et qui ont été valides a priori dans une configuration de nappe chargée de particules sont validés a posteriori dans la même configuration. Un analyse quantitative sur plusieurs cas avec diffèrent nombres de Stokes, nombres de Reynolds de la phase gazeuse et résolutions du maillage ont permit de retenir un modèle non-linéaire nommé 2EASM3, qui utilise le tenseur de déformations de la phase dispersée comme échelle de temps caractéristique. La deuxième partie a pour but l'amélioration du modèle d'évaporation implémenté dans AVBP. Ce modèle suppose une conduction infinie dans la phase liquide et symétrie sphérique dans la phase gazeuse ainsi que des lois simplifiées pour les propriétés thermodynamiques et de transport. Un nouveau modèle prenant en compte la dépendance de la viscosité du mélange gazeux avec la composition locale, et des nombres de Prandtl et Schmidt fixés par les valeurs à l'équilibre obtenus par moyen d'une simulation prenant en compte des lois complèxes pour les propriétés thermodynamiques et de transport est proposé. Cette nouvelle méthode produit des résultats en bon accord avec les mesures expérimentales pour l'évaporation d'une goutte isolé en une atmosphère d'azote au calme sans pourtant augmenter le cout du calcul. Finalement, l'impacte des nouveaux modèles est analysé dans une SGE de la configuration semi-industrielle MERCATO (García-Rosa (2008)). Bien que les données expérimentales ne soient pas suffisantes pour confirmer les résultats, les distributions de gouttes et de carburant gazeux sont significativement affectés par les modèles, ce qui pourrait avoir un impact directe sur le procès d'allumage. / Combustion still represents about 90% of the energy production in the world. Most industrial burners are fuelled with liquid hydrocarbons. However, most studies have been dedicated to gaseous ßames and the impact of liquid spray is still misunderstood. The purpose of this study is to improve the modelisation of two main phenomena occurring between atomization and combustion, i.e. the droplet dispersion in the turbulent gaseous flow and the evaporation process, in the context of Large Eddy Simulation (LES) of complex configurations. First, the mesoscopic Euler-Euler approach (Février et al. (2005)) based on a conditioned ensemble averaging and implemented in AVBP is improved. The closure model (Simonin et al. (2001), Kaufmann (2004)) for the second-order moments appearing in the transport equations solved fails in mean-sheared configurations (Riber (2007)). Several new models proposed by Masi (2010) and a priori tested in a particle-laden slab are tested a posteriori in the same configuration. A quantitative analysis based on several calculations varying the Stokes number, the gaseous Reynolds number and the grid resolution allows to retain a non-linear model using the particle rate-of-strain tensor as timescale and called 2EASM3. The second part consists in improving the evaporation model implemented in AVBP which assumes infinite conduction in the liquid and spherical symmetry in the gas phase along with simplified thermodynamics and transport properties calculation. A new model is proposed, where the dependence of gaseous mixture viscosity on local composition is accounted for, and the Prandtl and Schmidt numbers are fixed by a reference equilibrium calculation using complex thermodynamics and transport properties. This method shows good agreement with experimental measurements in the configuration of an isolated droplet evaporating in quiescent N2 without further increasing the computational cost. Finally, the impact of the new models is analysed in the LES of the MERCATO semi-industrial configuration (García-Rosa (2008)). Although the experimental data are not sufficient to confirm the results, both the droplet distribution and the fuel mass fraction are significantly affected, which would eventually affect the ignition process.

La réduction de charge de fluides frigorigènes dans les systèmes de production de froid est un enjeu important s'inscrivant dans les politiques environnementales sur la contribution de ces fluides à l'effet de serre. Une des voies menant à la réduction de charge est l'augmentation de la compacité des échangeurs de chaleur, conduisant inévitablement à la conception d'évaporateurs plus complexes. Néanmoins, si de nombreuses études ont été publiées sur l'analyse thermo-hydraulique de l'ébullition convective de frigorigènes dans les tubes horizontaux, très peu ont été menées sur des écoulements diphasiques dans des géométries non conventionnelles. Cette thèse vise à étudier les caractéristiques des écoulements diphasiques avec changement de phase des frigorigènes de troisième (HFC) et de quatrième génération tels que le HFO-1234yf dans des singularités (coudes de retour, contractions). Ainsi, un banc d'essais expérimental a été spécifiquement conçu et construit pour mener des expériences sur l'ebullition convective. Ce banc d'essais a permis de visualiser les régimes d'écoulement et de mesurer les chutes de pression dans les singularités. Des visualisations de régimes d'écoulement ont été réalisées pour les frigorigènes HFO-1234yf, R-134a et R-410A dans des tubes droits et dans des singularités. Les résultats ont été confrontés avec diverses méthodes de prédiction de la litérature. Les perturbations sur l'écoulement causées par ces singularités et leurs effets sur le comportement hydrodynamique des fluides frigorigènes ont également été etudiés. Une serie d'éxperience a été menée pour déterminer le longueur de perturbation en amont et en aval des singularités. Finalement, les valeurs expérimentales des chutes de pression ont été confrontées à différentes corrélations de la littérature. Des nouvelles méthodes de prédiction de la chute de pression pour des écoulements diphasiques dans des singularités sont présentées. / The refrigerant charge reduction in HVAC\R systems is an important issue because it falls within environmental policies regarding refrigerants contributions to the greenhouse effect. A way to move toward charge reduction is to increase the compactness of heat exchangers, which means more complex designs of the evaporators. Nevertheless, while a large amount of studies have has been published on the thermal and hydraulic analysis of flow boiling of refrigerants in horizontal tubes, very little attention has been given to flow boiling in geometries different from straight tubes. This PhD thesis aims at studying the flow boiling characteristics of third generation (HFC) and fourth generation fluids such as HFO-1234yf in geometries which modify the fluid dynamics and two-phase flow with respect to horizontal straight tubes. To achieve this goal, an experimental test facility was specifically designed and built to conduct refrigerant evaporation experiments. This test facility allowed to perform flow regimes visualizations and pressure drop measurements in singularities (such as sudden contractions and return bends). First, two-phase flow regimes visualizations have been carried out using HFO-1234yf, R-134a and R-410A either in straight tubes or in singularities. A qualitative analysis of the flow behavior and also several comparisons to flow pattern prediction methods from the literature were conducted. The second objective of this work was to characterize the flow disturbances caused by singularities such as sudden contractions and return bends, and to study their effects on the hydrodynamic performance (e.g. pressure drop) of refrigerants.Lastly, pressure drop databases for third and fourth generation refrigerants are presented.

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids.<p>We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a subtriangulation of the triangular Pk elements where we apply the distribution used for a P1 element.<p>This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes.<p><p>The first part is subdivided in four chapters:<p>The aim of the second chapter is to present the multidimensional upwind residual distributive schemes and to explain what was the status of their development at the beginning of this work.<p>The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes.<p>The fourth chapter is composed of two parts: we start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation.<p>This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion.<p>Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems.<p>In the second part, we extend the schemes of the first part to systems of equations as follows:<p>The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries.<p>Then, we look at the performance of 2nd and 3rd order schemes on viscous flow.<p>Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Depuis quelques années, la microplasturgie est un secteur en plein développement. Cependant, le comportement rhéologique des matériaux polymères dans des géométries très minces (dimension inférieure à 100 µm) n’est pas bien caractérisé. Peu de travaux ont été entrepris à ce sujet, en particulier en ce qui concerne les écoulements de Poiseuille qui sont pourtant les plus représentatifs des conditions de mise en œuvre usuelles. Ainsi, ce travail a pour but la mise au point d’une méthode expérimentale permettant d’obtenir des données pertinentes afin de caractériser de façon approfondie le comportement des matériaux polymères en écoulement de Poiseuille dans des géométries micrométriques. Afin de décrire au mieux la physique de tels écoulements, nous avons également cherché à les simuler numériquement, soit en utilisant des lois de comportement classiques, soit à l’aide de modèles se référant à la dynamique moléculaire. Pour atteindre ces objectifs, une filière à fente plate instrumentée avec des capteurs de pression et température, d’entrefer variant entre 50 et 200µm, a été conçue afin d’effectuer des mesures à l’aide d’un rhéomètre capillaire. Ce dispositif a été validé en confrontant les mesures à celles obtenues par d’autres méthodes (rhéométrie capillaire en filière classique et rhéométrie dynamique). Le glissement à la paroi a également été étudié, selon la méthode de Mooney. La simulation numérique de l’écoulement a d’abord été réalisée à l’aide de POLYFLOW®. L’effet de la pression sur la géométrie ainsi que sur le matériau polymère a été étudié. Puis, l’écoulement a également été simulé sous MATLAB® en utilisant des lois constitutives de type moléculaire basées sur le modèle du tube de Doï-Edwards ainsi que sur le concept de « Molecular Stress Function » introduit par Wagner pour rendre compte des effets d’orientation des molécules (variation du diamètre du tube) dans le champ de contraintes. L’écart constaté entre ces calculs et les résultats expérimentaux est expliqué et discuté à la lumière des simulations sous POLYFLOW®. Il modifie les perspectives d’étude de ce type d’écoulements. / The rheological behavior of polymer melts in microscale geometries is not really understood yet. In such processes which involve gaps thinner than 100µm (e.g. micro-injection molding), the material behaves differently compared to macroscopic flows. Besides, most polymer processing techniques involve pressure flows and only very few studies can be found about pressure flows in such thin geometries. The aim of this study was, first, to develop an experimental method which can provide relevant data about the rheological behavior of polymer melts in pressure flow taking place in microscale geometries. In order to get better descriptions of the physics involved in such flows, numerical simulation with commercial and home-made softwares was also implemented, especially with molecular dynamics constitutive models. Thus, a modular rheometrical slit die equipped with pressure and temperature transducers was designed to be adapted to a capillary rheometer, with different gap dimensions available, between 50µm and 200µm. The device was assessed by comparing to usual rheological ones, and wall slip was investigated according to Mooney’s method. Then, simulation of the flow was performed with POLYFLOW®. The pressure effect on the geometry and on the polymer material was investigated. Besides, simulation was also conducted with MATLAB® by implementing the Doi-Edwards’ tube model (reptation theory) and the Molecular Stress Function concept of Wagner to take into account the enhanced orientation of the molecules due to the very close vicinity of the die walls. Experimental results were compared to calculations, and the discussion of the discrepancies was supported by POLYFLOW® simulations. The conclusions somewhat modify the prospects for future studies of such flows.

L'objectif de cette thèse est l'étude du problème de Cauchy pour les solutions faibles de trois problèmes (systèmes paraboliques dégénérés et fortement couplés) modélisant des écoulements diphasiques et compressibles en milieu poreux. La motivation de ce travail est un "benchmark" du GNR MoMaS pour l'étude de l'impact de l'écoulement du gaz d\^{u} à la corrosion des matériaux ferreux dans un site de stockage de déchets radioactifs. Cette thèse est divisée en trois chapitres indépendants. Premièrement, on s'intéresse à l'analyse mathématique d'un problème modélisant l'écoulement de deux phases immiscibles et en considérant qu'une phase est compressible et l'autre est incompressible (eau/gaz). Deuxièmement, on traite le cas général du déplacement de deux fluides compressibles et immiscibles dans un milieu poreux. Enfin, le dernier chapitre est consacré à la construction et à la convergence de la méthode des volumes finis pour le système eau-gaz sous l'hypothèse que la densité du gaz est une fonction de la pression globale.

En aérodynamique automobile, les écoulements sont caractérisés par des structures fortement décollées, tridimensionnelles et instationnaires qui sont à l'origine de la trainée du véhicule et de la génération de bruit. Une configuration modèle générique de cette aérodynamique a été proposée. On étudie de manière expérimentale l'interaction du sillage tridimensionnel d'un disque de diamètre D placé à une distance H d'une paroi. Sur la base de cet écoulement, l'objectif est de développer une analyse instationnaire des structures de l'écoulement et de leur couplage avec la pression fluctuante. On a mis en oeuvre pour cela des moyens d'essais modernes tels que la mesure synchronisée pression/vitesse ou la mesure HS-PIV (High Speed - PIV) et développé des moyens d'analyse (analyse corrélatoire, outils lagrangien FTLE, moyenne conditionnelle et estimation stochastique). Pour H/D=0.75, le battement du sillage associé au lâcher tourbillonnaire est mis en évidence. Les moyennes de phase, déterminées à l'aide d'une analyse POD du signal de pression, associent le déplacement vertical du sillage au lâcher de structures cohérentes dans le plan de symétrie. L'utilisation de l'estimation stochastique de type spatio-temporelle, dont la supériorité par rapport à l'approche spatiale a été démontrée, a permis de montrer que les structures cohérentes sont reliées de manière linéaire aux fluctuations de pression. Pour H/D=0.3, l'interaction du sillage du disque avec la paroi a considérablement changé. Dans la zone proche du point de décollement, il a été montré à l'aide de l'outil lagrangien FTLE que la fréquence dominante de 30Hz sur les DSP de pression est associée aux fluctuations du point de décollement.

The present research is an experimental and numerical study of aeroacoustic phenomena occurring in large solid rocket motors (SRM) as the Ariane 5 boosters. The emphasis is given to aeroacoustic instabilities that may lead to pressure and thrust oscillations which reduce the rocket motor performance and could damage the payload. The study is carried out within the framework of a CNES (Centre National d'Etudes Spatiales) research program. Large SRM are composed of a submerged nozzle and segmented propellant grains separated by inhibitors. During propellant combustion, a cavity appears around the nozzle. Vortical flow structures may be formed from the inhibitor (Obstacle Vortex Shedding OVS) or from natural instability of the radial flow resulting from the propellant combustion (Surface Vortex Shedding SVS). Such hydrodynamic manifestations drive pressure oscillations in the confined flow established in the motor. When the vortex shedding frequency synchronizes acoustic modes of the motor chamber, resonance may occur and sound pressure can be amplified by vortex nozzle interaction. Original analytical models, in particular based on vortex sound theory, point out the parameters controlling the flow-acoustic coupling and the effect of the nozzle design on sound production. They allow the appropriate definition of experimental tests. The experiments are conducted on axisymmetric cold flow models respecting the Mach number similarity with the Ariane 5 SRM. The test section includes only one inhibitor and a submerged nozzle. The flow is either created by an axial air injection at the forward end or by a radial injection uniformly distributed along chamber porous walls. The internal Mach number can be varied continuously by means of a movable needle placed in the nozzle throat. Acoustic pressure measurements are taken by means of PCB piezoelectric transducers. A particle image velocimetry technique (PIV) is used to analyse the effect of the acoustic resonance on the mean flow field and vortex properties. An active control loop is exploited to obtain resonant and non resonant conditions for the same operating point. Finally, numerical simulations are performed using a time dependent Navier Stokes solver. The analysis of the unsteady simulations provides pressure spectra, sequence of vorticity fields and average flow field. Comparison to experimental data is conducted. The OVS and SVS instabilities are identified. The inhibitor parameters, the chamber Mach number and length, and the nozzle geometry are varied to analyse their effect on the flow acoustic coupling. The conclusions state that flow acoustic coupling is mainly observed for nozzles including cavity. The nozzle geometry has an effect on the pressure oscillations through a coupling between the acoustic fluctuations induced by the cavity volume and the vortices travelling in front of the cavity entrance. When resonance occurs, the sound pressure level increases linearly with the chamber Mach number, the frequency and the cavity volume. In absence of cavity, the pressure fluctuations are damped.

L’objectif est de construire un modèle d’écoulement qui se prête bien à des problèmes de contrôle, en associant un faible nombre de degrés de liberté à la possibilité de décrire la dynamique d’un écoulement relativement complexe. Dans ce travail nous considérons un écoulement bidimensionnel laminaire autour d’un cylindre carré. Des actionneurs placés sur le cylindre permettent un contrôle actif par sou?age et aspiration. Ce contrôle peut être dé?ni par rétroaction, exploitant des mesures de la vitesse dans le sillage du cylindre. Nous construisons un modèle d’ordre réduit (ROM) des équations de Navier-Stokes incompressibles, basé sur la technique de décomposition orthogonale aux valeurs propres (POD). Une façon classique de construire un tel modèle est de réaliser une projection Galerkin des équations sur le sous-espace réduit obtenu par POD. Un tel modèle peut cependant être peu précis, voire instable. Une technique de calibration est alors mise en place pour assurer la bonne représentativité dynamique du modèle. Nous dé?nissons ensuite une stratégie pour mettre à jour le modèle au cours d’un processus d’optimisation. La méthode est en?n appliquée pour réduire la di?érence entre l’écoulement contrôlé et la solution stationnaire instable à Re = 150. / The aim is to build a ?ow model adapted for control applications combining a low number of degrees of freedom with the possibility of describing relatively complex ?ows. In this work a two-dimensional laminar ?ow past a square cylinder is considered. Actuators placed on the cylinder enable active control by blowing and suction. Proportional feedback control can then be applied using velocity measurements taken in the cylinder wake. The proper orthogonal decom- position (POD) approach is used to build a low order model of the incompressible Navier-Stokes equations. A classical way of obtaining a Reduced-Order Model (ROM) is to perform a Galerkin projection of the equations onto the subspace spanned by the POD modes. Such a model can however be inaccurate, even unstable. A calibration technique is therefore applied, leading to a model that is accurate and robust to variations of the control parameters. A strategy is then de?ned to update the model within an optimisation loop. The method is tested at Re = 150 for reducing the di?erence between the actuated ?ow ?eld and the steady unstable solution.

Ce travail est consacré à l'étude expérimentale des écoulements diphasiques et triphasiques d'inclusions (bulles, gouttes, sphères) en milieux tant newtoniens que non newtoniens à l'échelle microscopique et mésoscopique, en utilisant la visualisation par une caméra rapide, la vélocimétrie par images des particules (PIV) ainsi que la micro-vélocimétrie par images des particules. Des bulles et des gouttes ont été étudiées expérimentalement depuis leur formation, en passant par leur déformation jusqu'à leur coalescence. La formation de bulles dans des micro-mélangeurs a été étudiée et caractérisée par l'obtention de champs de vitesses. Différents paramètres, tels que le cisaillement, la géométrie de la zone de formation, les débits ou encore les propriétés physiques ont été testés afin de développer des lois d'échelles. La traversée d'une interface liquide-liquide par une inclusion a été abordée par des expériences originales, permettant de décrire la dynamique du phénomène, de définir des nombres adimensionnels et de mettre en évidence des instabilités interfaciales. L'effet Weissenberg a également été étudié aux différentes échelles afin de comprendre les phénomènes conduisant à son amplification. Enfin, en milieu viscoélastique et rhéofluidifiant, nous avons caractérisé l'écoulement autour d'une inclusion isolée solide par l'obtention de champs de vitesses. Ces résultats ont permis de confirmer l'origine viscoélastique du sillage négatif et de prédire ses caractéristiques / The present work was devoted to the experimental study of the multiphase flow around inclusions in both Newtonian and non-Newtonian media at respectively microscopic and mesoscopic scales, by means of the Particle Image Velocimetry (both PIV and µ-PIV) and fast camera visualization. Bubbles and drops were experimentally studied starting from their formation, the rising and up to their coalescence and fragmentation. Bubble formation in micro-mixers was also investigated and characterized by measuring the liquid velocity fields. Different parameters, such as the shear rate, the geometry of the micro-mixer, the flow rates or the physical properties were tested to develop correlations of power-law kind. The deformation of a liquid-liquid interface due to the passage of an inclusion was investigated by original experiments in order to observe and describe the dynamics of the phenomenon, to define the reliable dimensionless numbers and also to highlight several interfacial instabilities. The Weissenberg effect was also studied at different scales in various non-Newtonian fluids in comparison with Newtonian fluids to understand the amplification phenomena under the combining effects of instability and normal forces. Finally, the flow around an isolated solid inclusion was characterized by performing the measurements of velocity fields in viscoelastic and shear-thinning fluids. These results were used to confirm the viscoelastic origin of the negative wake behind the solid sphere and to model its characteristics

Cette thèse porte sur le développement d’un modèle numérique de l’affouillement causée par des obstacles montes sur le lit, combinant les processus hydrodynamiques et morphologiques. Le modèle numérique est basé sur le solveur de champ d’écoulement polyphasique de l’outil CFD open-source OpenFOAMR qui est distribue par OpenCFD Ltd. Le module hydrodynamique du modèle résout les équations de Navier-Stokes avec moyennes de Reynolds (RANS) et les modèles des turbulences k-ε ou k-ω. Il existe deux interfaces dans le domaine de simulation: la surface libre entre l’eau et l’air, qui est suivi par la méthode de Volume de Fluide (VOF); et l’interface entre l’eau et le lit du sédiment, qui est représentée par un maillage de surface finie déformable construit à partir de la limite en bas du maillage de volume fini. En outre, un module morphologique qui a été développé dans le cadre du projet se compose de trois composantes: un modèle de transport de sédiments comprenant la charge suspendue et le charriage; l’équation d’Exner pour mesurer la déformation du lit; et un mécanisme de glissement du sable pour limiter la pente du lit à être plus petite que l’angle de repos du sédiment. Le changement morphologique est incorporé dans le modèle hydrodynamique par la déformation du maillage. Des conditions limites spéciales et des corrections nécessaires pour le calcul en parallèle sont également ajoutées au modèle. Chaque partie du modèle est validée séparément avec les tests préliminaires correspondants, y compris les fonctions de paroi rugueuse, les performances de la méthode VOF, le modèle de transport de charge suspendu et le mécanisme de glissement de sable. Le modèle numérique est ensuite appliqué pour étudier un affouillent bidimensionnelle cause par un jet immerge provenant d’une ouverture sous écluse. Comparaison des résultats de la simulation avec des données expérimentales prouve la capacité du modèle. Et les limites du modèle sont également discutées. Enfin, le modèle est appliqué à l’étude du champ d’écoulement tridimensionnel et de la formation d’affouillement autour d’un obstacle dans l’écoulement. Tout d’abord, la déformation du lit n’est pas activée. Le tourbillon en fer à cheval devant un obstacle et le champ d’écoulement turbulent autour d’un cylindre sur un lit lisse ou rugueux sont simulés. Deux types de simulation pour le module hydrodynamique sont effectués: une simulation qui utilise une surface fixe et rigide pour représenter l’interface air-eau, et une simulation incluant à la fois les domaines de l’eau et de l’air avec la surface libre suivie par la méthode VOF. Les influences de la surface libre sur le champ d’écoulement sont identifiées et discutées. La comparaison avec les données expérimentales confirme l’importance de la déformation de la surface libre sur le champ d’écoulement. Ensuite, le lit est autorisé à se déformer et le développement temporel de l’affouillement tridimensionnelle autour d’un cylindre sur le lit est simule. Le développement temporel d’affouillement et la profondeur maximale du trou calcule devant et derrière le cylindre conviennent assez bien avec les mesures expérimentales. Les influences de l’affouillement sur le champ d’écoulement sont aussi étudiées et la performance du modèle numérique développé est discutée. / This thesis describes the development of a numerical model for local scour caused by bed-mounted obstacles, combining the hydrodynamic and morphological processes. The basis of the numerical model is the multiphase flow field solver in the open-source CFD toolbox OpenFOAMR which is released by OpenCFD Ltd. The hydrodynamic module of the model solves the Reynolds Averaged Navier-Stokes (RANS) equations with either a k-ε or a k-ω model. There are two interfaces in the simulation domain: the free surface between water and air, which is tracked using the Volume of Fluid (VOF) method, and the interface between the water and the sediment, which is represented by a finite area mesh constructed from the bottom boundary of the finite volume mesh. A morphological module which has been developed as part of the project consists of three components: a sediment transport model which includes suspended load and bed load transport; the Exner equation to compute the bed deformation, and a sand-sliding mechanism to restrict the bed slope angle to be smaller than the angle of repose. The morphological changes are incorporated into the hydrodynamic field through deformation of the computational mesh. Additional boundary conditions and parallel computing corrections are also added into the model. Each individual part of the model has been validated separately with corresponding preliminary test cases including the rough wall functions, the performance of the VOF method, the suspended load transport model and the sand-sliding mechanism. The numerical model is then applied to study two-dimensional scour caused by a submerged jet issuing from an opening under sluice gate. Comparison of the simulation results with the experimental measurements proves the ability of the model for conducting two-dimensional simulations and the limitations of the model are also discussed. Finally, the model is applied to study the three-dimensional flow field and scour formation around an obstacle in flow. Initially, the bed deformation is not activated in the model. The horseshoe vortex formed in front of an obstacle in water and the turbulent flow field around a cylinder on smooth and rough beds are simulated. Two types of simulations for the hydrodynamic module are used: a rigid lid simulation with a slip boundary condition to represent the air-water interface, and a free surface simulation including both the water and air domains with the free surface tracked by the VOF method. The influences of the variation of the water depth on the flow field are identified and discussed. Comparison with the experimental data also confirms the importance of the water surface variation on the flow field. Next, the bed is allowed to deform in the model. The temporal development of three-dimensional scour around a cylinder on live-bed in a steady current is simulated. The development of the scour with time and the computed maximum scour depths in front of and behind the cylinder agree quite well with the experimental measurements. The influences of the scour process on the flow field are also studied and the performance of the numerical model is discussed.