Cette thèse est consacrée au calcul et à l'étude du champ acoustique rayonné par un jet fortement supersonique, sur-détendu et chauffé. L'approche numérique utilisée consiste à résoudre simultanément l'écoulement turbulent et le champ proche acoustique par simulation des grandes échelles.Afin de résoudre les équations de Navier-Stokes instationnaires et compressibles, des schémas numériques de dérivation spatiale et d'intégration temporelle peu dissipatifs et peu dispersifs sont utilisés. Pour le système de coordonnées cylindriques, une procédure de déraffinement progressive au voisinage de l'axe du jet est proposée pour augmenter le pas de temps de l'intégration temporelle. Une procédure de capture de choc est finalement développée pour supprimer les oscillations de Gibbs aux voisinages des chocs. Cette procédure est explicite, adaptative et conservative.Un jet rond, sur-détendu et présentant un nombre de Mach d'éjection de 3.30 ainsi qu'une pression et une température de réservoir de 28.6 × 105Pa et 1144K est résolu par simulation des grandes échelles. Ce jet est initialement laminaire et possède un nombre de Reynolds de 0.94 × 105. Les champs aérodynamique et acoustique obtenus sont comparés aux données disponibles dans la littérature. Afin de pouvoir identifier les sources acoustiques,le champ turbulent est caractérisé le long de la couche de mélange, notamment en utilisant une analyse spectrale et une décomposition selon les modes azimutaux. Les champs acoustiques proche et lointain sont enfin étudiés et différentes composantes de bruit sont identifiées à partir de comparaisons avec des mesures et des modèles de la littérature, ainsi qu'en établissant des liens avec les propriétés de l'écoulement turbulent. L'influence des effets de la propagation non linéaire des ondes acoustiques est finalement discutée.

Dans le but d'améliorer le rendement des propulseurs aérospatiaux, on s'intéresse à l'utilisation de ladétonation dans le cycle moteur. Cette thèse porte sur le développement et l'utilisation d'un codepour la compréhension du fonctionnement d'un moteur à détonation continue (CDWE). Le 1erchapitre place le cadre de l'étude, et positionne le CDWE par rapport à différents concepts demoteurs à détonation. Un état des lieux des simulations numériques concernant le fonctionnementd'un CDWE est établi afin de justifier l'approche numérique à utiliser. Cette approche numérique estdétaillée dans le 2e chapitre. Les équations d'Euler, les modèles thermochimiques, ainsi que lesschémas cinétiques utilisés dans cette étude y sont présentés. Le 3e chapitre décrit les méthodesnumériques implémentées dans le code. Le schéma WENO d'ordre 5 est utilisé pour l'évaluation desflux numériques. L'avancement temporel est assuré par le schéma semi-implicite d'ordre 2 ASIRK2Cou explicite d'ordre 3 RK3. Le 4e chapitre est consacré à la technique de raffinement adaptatif demaillage (AMR) et à la bibliothèque choisie. Le code est testé dans le 5e chapitre sur différents cas etappliqué à la simulation d'une onde de détonation afin de préparer les simulations présentées dans ledernier chapitre. Le 6e chapitre présente les résultats des simulations d'un CDWE. La structure 2Dd'une onde de détonation continue est présentée et comparée avec la structure 3D. L'influence durayon de courbure du canal et l'effet d'une injection par une fente sur la structure de l'écoulementsont étudiés.

L'étude effectuée au cours de cette thèse s'inscrit dans le contexte du contrôle des écoulements en turbomachines. Suite à l'intérêt avéré de l'aspiration de la couche limite en compresseurs axiaux pour en améliorer les performances, elle vise à évaluer si cette méthodologie de contrôle est pertinente en géométrie radiale. Elle s'appuie pour cela sur un compresseur centrifuge transsonique Turbomeca. Le fonctionnement du compresseur est d'abord analysé sur la base de simulations numériques stationnaires et instationnaires chorochroniques, réalisées grâce au code de calcul elsA développé par l'Onera et en utilisant le modèle de turbulence k-l de Smith. La validité des résultats numériques est vérifiée par comparaison avec les données expérimentales disponibles : performances mono-dimensionnelles et signaux de pressions instationnaires. Ces données sont issues à la fois de la caractérisation expérimentale du compresseur effectuée par Turbomeca et de campagnes de mesures menées a l'ISAE, ou le compresseur est monté sur un banc de recherche dédié à l'étude des interactions entre le rouet et le diffuseur radial. L'analyse détaillée des résultats numériques révèle le développement d'un décollement de coin dans le diffuseur, entre la paroi moyeu et la face en dépression des aubages, lorsque le point de fonctionnement se rapproche du pompage. L'étude topologique du spectre du frottement pariétal permet d'en préciser le lieu et la structure, et de localiser les points singuliers ainsi que la ligne de séparation principale. En stationnaire, la croissance du décollement aboutit au décrochage du diffuseur, et compromet les performances de l'étage de compression et finalement la stabilité numérique. L'exploitation des résultats instationnaires permet de préciser le comportement temporel de l'écoulement : l'étendue du décollement fluctue sous l'effet du défilement des ondes de pression générées par l'interaction rouet-diffuseur, mais les trajectoires des particules fluides impliquées dans le décollement sont en accord avec les lignes de courant du champ de l'écoulement moyenné temporellement. Le décollement instationnaire est fixe, et sa topologie correspond à celle prévue par les simulations stationnaires. En particulier, le lieu du décollement est prévu de façon similaire par les modèles stationnaires et instationnaires. Une stratégie d'aspiration est ensuite développée sur la base de ces observations : une fente de prélèvement est positionnée au voisinage du col de séparation principal, dont le lieu correspond à celui du maximum du gradient de pression adverse stationnaire. Cette stratégie est implémentée au sein des modèles numériques stationnaires et instationnaires. En stationnaire, un prélèvement de 1 % de la valeur totale du débit traversant le rouet permet un contrôle total du décollement dans le diffuseur, et conduit à une augmentation significative de la plage de stabilité numérique. La réduction du débit de prélèvement à 0,3 % du débit total permet un contrôle partiel du décollement, et conduit également à une augmentation de la plage de stabilité numérique. En instationnaire, le décollement de coin initial est contrôlé. Mais les résultats mettent en évidence le rôle majeur joue par le défilement des ondes de pression le long des aubages du diffuseur. Celles-ci se renforcent au passage du col du diffuseur, et engendrent l'existence d'un important gradient de pression adverse instantané. Ce maximum du gradient de pression instantané conduit au développement d'un nouveau décollement de coin, en aval de la fente de prélèvement, et la plage de stabilité des calculs instationnaires n'est pas augmentée. Ces résultats mettent en évidence la possibilité d'agir sur les décollements se produisant dans les diffuseurs radiaux à l'aide de la technique d'aspiration de la couche limite. [...]

Ce travail a pour but de comparer les approches de simulation numérique des grandes échelles explicite (LES) et implicite (ILES) pour un jet turbulent non-réactif ou réactif d’hydrogène à grande vitesse dans un co-courant d’air, typique d’un super-statoréacteur. La résolution des calculs va de 32 × 32 × 128 à 256 × 256 × 1024, à l’aide d’un schéma WENO d’ordre 5. Les LES explicites emploient les modèles sous-maille de Smagorinsky et de Fonction de Structure Sélective, associés au transport moléculaire. Les LES implicites sont réalisées avec et sans modèle de transport moléculaire, en résolvant les équations de Navier- Stokes ou d’Euler. Dans le cas non-réactif, le modèle de Smagorinsky est trop dissipatif. Le modèle de Fonction de Structure Sélective améliore les résultats, sans faire mieux que l’approche ILES quelle que soit la résolution. Dans le cas réactif, une coupure physique visqueuse est indispensable pour fixer une épaisseur à la flamme, et assurer la convergence en maillage de l’approche ILES. On montre aussi que les résultats LES/ILES sont moins sensibles aux conditions d’injection que ceux de l’approche RANS. Le premier chapitre est une introduction générale au contexte de l’étude. Au second chapitre, on rappelle les équations générales pour un écoulement réactif et on détaille les modèles thermodynamique et de transport retenus. Au troisième chapitre, les équations de la LES et les modèles sous-maille sont présentés. On examine également quelques propriétés du schéma numérique. Le chapitre 4 est consacré à la méthode numérique et au code de calcul. Enfin, on présente les cas-tests et on discute les résultats au chapitre 5. / This work is intended to compare Large Eddy Simulation and Implicit Large Eddy Simulation (LES and ILES) for a turbulent, non-reacting or reacting high speed H2 jet in co-flowing air, typical of scramjet engines. Numerical simulations are performed at resolutions ranging from 32 × 32 × 128 to 256 × 256 × 1024, using a 5th order WENO scheme. Physical LES are carried out with the Smagorinsky and the Selective Structure Function models associated to molecular diffusion. Implicit LES are performed with and without molecular diffusion, by solving either the Navier-Stokes or the Euler equations. In the nonreacting case, the Smagorinsky model is too dissipative. The Selective Structure Function leads to better results, but does not show any superiority compared to ILES, whatever the grid resolution. In the reacting case, a molecular viscous cut-off in the simulation is mandatory to set a physical width for the reaction zone in the ILES approach, hence to achieve grid-convergence. It is also found that LES/ILES are less sensitive to the inlet conditions than the RANS approach. The first chapter is an introduction to the context of this study. In the second chapter, the governing equations for multispecies reacting flows are presented, with emphasis on the thermodynamic and transport models. In the third chapter, physical LES equations and explicit sub-grid modeling strategies are detailed. Some properties of the numerical scheme are also investigated. In chapter four, the numerical scheme and some aspects of the solver are explained. Finally, non-reacting and reacting numerical experiments are presented and the results are discussed.

De manière à respecter les normes européennes en matière d’émission de gaz polluants, les constructeurs automobiles pratiquent le downsizing. Cette pratique consiste en une réduction de la cylindrée des moteurs tout en maintenant un bon niveau de performance. Il s’en suit des puissances spécifiques moteur importantes, l’objectif cible étant de l’ordre de 100 kW/l. Pour de telles puissances, les températures atteintes au niveau des pistons sont élevées, ce qui pose alors le problème de leur tenue thermomécanique. Le refroidissement du piston devient donc un acteur important de la fiabilité du moteur. Le procédé le plus répandu actuellement est le refroidissement par jet d’huile. Le piston est alors refroidi par un écoulement turbulent diphasique incompressible (Air/Huile) dans un environnement mobile. Il s’agit de phénomènes physiques complexes qui sont pour l’instant mal connus. L’optimisation du refroidissement nécessitant une bonne compréhension des phénomènes physiques concernés, nous nous proposons dans ce mémoire de le modéliser au moyen du modèle 1-Fluide diphasique couplé à l’équation de l’énergie. / The present work aims at studying the cooling system used in cars engine to ensure piston thermo mechanical resistance by numerical simulation. Most of actual engines use an oil jet cooling system coupled with 'cocktail shaking' to extract heat from piston. This cooling method brings into play a two-phase incompressible turbulent flow in a mobile environment, due to motion of pistons in the cylinder. The need today for more effective cooling of pistons involves an accurate understanding of the physical mechanisms which are concerned. Modeling could be a good way to achieve it. The idea is to support the engine design process to account for advanced technologies to improve turbine or engine performances, less fuel burn and green house gases. In the present work, a numerical model dedicated to the simulation at small scale of oil/air two-phase flows and related heat transfers is proposed to characterize the cooling of engine elements under fragmented jet impact.

Cette travail de thèse est lié au besoin de modélisation des écoulements diphasiques en générateurs de vapeur (entrée liquide et sortie vapeur). La démarche proposée consiste à faire le choix d'une modélisation hybride de l'écoulement, en scindant la phase gaz en deux champs, modélisés de manières différentes. Ainsi, les petites bulles sphériques sont modélisées avec une approche dispersée classique avec le modèle eulérien à deux fluides, et les bulles déformées sont simulées à l'aide d'une méthode de localisation d'interface.Le travail effectué porte sur la mise en place, la vérification et la validation du modèle dédié aux larges bulles déformées, ainsi que le couplage entre les deux approches pour le gaz gaz, permettant des premiers calculs de démonstration utilisant l'approche hybride complète. / The present study is related to the need of modeling the two-phase flows occuring in a steam generator (liquid at inlet and vapour at outlet). The choice is made to investigate a hybrid modeling of the flow, considering the gas phase as two separated fields, each one being modeled with different closure laws. In so doing, the small and spherical bubbles are modeled through a dispersed approach within the two-fluid model, and the distorted bubbles are simulated with an interface locating method.The main outcome is about the implementation, the verification and the validation of the model dedicated to the large and distorted bubbles, as well as the coupling of the two approaches for the gas, allowing the presentation of demonstration calculations using the so-called hybrid approach.

Le contexte général de cette thèse est la simulation numérique de l’injection de carburant dans un moteur à combustion interne, afin d’améliorer son rendement et de limiter la production de polluants. Intrinsèquement, il est possible de simuler l’ensemble de l’écoulement avec les équations classiques de la dynamique des fluides sans avoir recours à des outils de modélisation supplémentaires liés au caractère diphasique. Mais, les tailles des structures générées pendant l’injection (gouttes de diamètre < à 10 μm) conduisent à des temps de calculs prohibitifs pour une application industrielle. C’est pourquoi il est nécessaire d’introduire une modélisation diphasique. C’est dans ce contexte que deux régions sont formellement distinguées: le cœur liquide dense proche de l’injecteur, appelé écoulement à phases séparées, et le spray constitué d’une population de gouttes polydisperse générées après le processus d’atomisation en aval de l’injecteur. Ce travail de thèse étudie les modèles Eulériens pour la description de spray évaporants et polydisperses, en vue d’applications industrielles. Ils représentent une alternative potentielle aux modèles Lagrangiens qui sont majoritairement utilisés en industrie mais présentant des inconvénients majeurs. Ainsi, le modèle multi-fluide est étudié dans un premier temps. Bien que prometteur, deux difficultés sont soulignées: le coût requis pour une description précise de la polydispersion, et son incapacité à décrire les croisements de gouttes (particle trajectory crossing, PTC). La thèse propose des solutions à ces deux limitations. Elles reposent sur des méthodes de moments. Premièrement, le modèle appelé Eulerian Size Multi Size Moment (EMSM) permet de résoudre des sprays évaporants et polydisperses de manière bien plus efficace que le modèle multi-fluide. Des outils mathématiques sont utilisés pour fermer le système d’équations associé au modèle, et combinés à des schémas de types volumes finis appelés schémas cinétiques, afin de préserver la réalisabilité du vecteur de moments, pour le transport et l’évaporation. Une réponse à la seconde limitation est apportée avec le modèle appelé Eulerian Multi Velocity Moment (EMVM) basé sur le transport de moments en vitesse d’ordre élevé. Une distribution bimodale peut être localement reconstruite à partir des moments en utilisant une méthode de quadrature de moments ( QMOM) en une ou plusieurs dimensions d’espace. De la même manière, l’utilisation de schémas cinétiques permet de préserver la réalisabilité du vecteur de moment. De plus, une étude mathématique approfondie de la dynamique du système en une dimension d’espace en révèle toute la complexité et représente une étape indispensable en vue de l’élaboration de schémas de transport d’ordre élevé (supérieur ou égal à 2).Afin de les tester, ces deux modèles ainsi que les outils numériques associés sont implémentés dans MUSES3D, un code académique de simulation numérique directe (DNS) dédié à l’évaluation des modèles de spray. Des résultats de grande qualité démontrent le potentiel des modèles. L’extension du modèle EMSM dans un contexte industriel est ensuite considérée, avec son implémentation dans IFP-C3D, un code résolvant des écoulements réactifs sur des maillages non structurés et mobiles dans un formalisme RANS (Reynolds Averaged Navier Stokes) en présence de sprays. Le formalisme ALE (Arbitrary Lagrangian Eulerian) est utilisé et le modèle EMSM réécrit dans ce formalisme afin de mener des calculs en maillage mobile. De plus, une étude numérique a permis d’étendre les propriétés de précision et de stabilité obtenues en maillage fixe. La robustesse du modèle EMSM est alors démontrée avec succès dans IFP-C3D sur un cas impliquant un mouvement de piston, ainsi que dans le cadre d’une comparaison avec le code MUSES3D. Enfin, des résultats très encourageants prouvent la faisabilité d’un calcul d’injection dans une chambre de combustion d’un spray polydisperse avec le modèle EMSM. / The general context of the PhD is the simulation of fuel injection in an internal combustion engine, in order to improve its thermal and ecological efficiency. This work more generally concerns any industrial device involving a multiphase flow made of liquid fuel injected in a chamber filled with gaz: automotive or aircraft engines, or turbo machines. In and of itself, it is possible to simulate this flow without any modeling. However the small structures created during injection (droplets of diameter until 10 μm or less) lead to a prohibitive computational cost for any industrial application. Therefore modeling is necessary. In this context, two areas are formally distinguished: the dense liquid core close to the injector called separate-phase flow, and the spray made of a polydisperse droplet population (i.e. droplets with different sizes) generated after the atomization processes downstream of the injector. This PhD work investigates Eulerian models for the description of polydisperse evaporating sprays, for industrial computations. They represent a potential alternative to Lagrangian models, widely used at present, yet suffering from major drawbacks. Thus, the Multi-Fluid model is assessed. Although it is very promising, two difficulties are highlighted: its cost for a precise description of polydispersity, and its inability to describe particle trajectory crossing (PTC). Solutions to these two limitations are considered. Both rely on high order moment methods. First, the Eulerian Multi Size Moment (EMSM) proposes a much more efficient resolution of polydisperse evaporating sprays than the Multi-Fluid model does. Mathematical tools are used to close the model and combined with original finite volume kinetic-based schemes in order to preserve the moment-set integrity, for evaporation and advection. An answer to the second limitation is provided with the Eulerian Multi Velocity Moment (EMVM) based on high order velocity moments. A bimodal velocity distribution can be locally reconstructed for the moments using the quadrature method of moments (QMOM), in one or multi-dimensions. Here also, finite volume kinetic-based schemes are studied in order to preserve the moment set integrity. Moreover, a mathematical study of the one-dimensional dynamic system highlights its peculiarity and constitutes a necessary basis for the design of high order numerical schemes. In order to assess them, both the models and their numerical tools are implemented in the MUSES3D code, an academic DNS solver that provides a framework devoted to spray method evaluation. Achievements of the EMSM and the EMVM models are presented. The extension of the EMSM model to an industrial context is then considered, with its implementation in the IFP-C3D code, a 3D unstructured reactive flow solver with spray. In order to perform computations within a moving domain (due to the piston movement) the Arbitrary Lagrangian Eulerian (ALE) formalism is used. A numerical study has been achieved, in order to extent to this formalism the properties of accuracy and stability of the EMSM model, which already induces strong stability condition in an Eulerian approach. The robustness of the EMSM model in the IFP-C3D code has been successfully demonstrated on a case involving a moving piston, and also on a comparison with the MUSES3D code. Moreover, very encouraging results demonstrate the feasibility of the EMSM model for spray injection.

Au cours de ce travail, nous nous sommes attaché à la résolution numérique des équations de la Magnétohydrodynamique (MHD) auxquelles s'ajoute une loi hyperbolique de transport des erreurs de divergence.La première étape consista à symétriser le nouveau système de la MHD idéale afin d'en étudier le système propre, ce qui fut l'occasion de rappeler le rôle de l'entropie au niveau de ce calcul comme à celui de l'inégalité de Clausius-Duhem. La suite de cette thèse eut pour objectif la résolution de ces équations idéales à l'aide de schémas distribuant le résidu (notés RD). Les quatre principaux schémas connus furent testés, et nous avons montré entre autres que le schéma N, qui a fait ses preuves sur les équations d'Euler en mécanique des fluides, n'était pas adapté aux équations de la MHD. Les stratégies classiques de limitation et de stabilisation purent être revisitées à ce moment. Les équations étant instationnaires, il fallut intégrer une discrétisation en temps et une distribution spatiale des termes d'évolution (et d'éventuelles sources). Nous avons d'emblée opté pour une approche implicite permettant d'être performant sur les simulations longues des expériences de tokamaks, et de traiter la correction de la divergence d'une manière originale et efficace. Les problèmes de convergence de la méthode de Newton-Raphson n'ayant pas été pleinement résolus, nous nous sommes tournés vers une alternative explicite de type Runge-Kutta. Enfin, nous avons réétabli les principes de la montée en ordre (en théorie, jusqu'à des ordres arbitraires, en prenant en compte le phénomène de Gibbs) à l'aide de tout type d'élément fini (bien construit) 2D ou 3D, sans avoir pu valider tous ces aspects. Nous avons également pris en compte les équations complètes de la MHD réelle classique (i.e. sans effet Hall) à l'aide d'un couplage RD/Galerkin. / During this thesis, we worked on the numerical resolution of the Magnetohydrodynamic (MHD) equations, to which we added a hyperbolic transport equation for the divergence errors of the magnetic field.The first step consisted in symmetrizing the new ideal MHD system in order to study its eigensystem, which was the opportunity to remind the role of the entropy in this calculation as well as in the Clausius-Duhem inequality. Next, we aimed at solving these ideal equations by the mean of Residual Distribution (RD) schemes.The four main schemes were tested, and we showed among other things that the N scheme (although it has been proven very efficient with Euler equations in Fluid Mechanics) could not give satisfying results with the MHD equations. Classical strategies for the limitation and the stabilization were revisited then. Moreover,since we dealt with unsteady equations, we had to formulate atime discretization and a spatial distribution of the unsteady terms (as well as possible sources). We first choosed an implicit approach allowing us to be powerful on the long simulations needed for tokamak experiments, and to treat the divergence cleaning part in an original and efficient way. The convergence problems of our Newton-Raphson algorithm having not been fully resolved, we turned to an explicit alternative (Runge-Kutta type).Finally, we discussed about the principles of higher order schemes (theoretically, up to arbitrary orders, taking into account the Gibbs phenomenon) thanks to any type of 2D or 3D finite element (properly defined), without having been able to to validate all these aspects. We also implemented the dissipative part of the full MHD equations (in the classical sense, i.e. omitting the Hall effect) by the use of a RD/Galerkin coupling.

La thèse est consacrée à la recherche des catastrophes marines (tsunamis, ondes de tempête) dans les Antilles françaises, en utilisant des modèles analytiques et numériques de mécanique des fluides. L'accent est mis sur le développement de modèles de mouvement de glissements de terrain et des vagues causées par ces glissements. Le modèle le plus simple du glissement de terrain «solide block» est utilisé pour évaluer les caractéristiques des flux pyroclastiques du volcan Soufrière-Hills (Montserrat). Un modèle plus complexe de modélisation de glissement de terrain (modèle de Savage-Hutter) a été étudié analytiquement, donnant ainsi une nouvelle famille de solutions exactes décrivant le mouvement de l'écoulement par gravité non linéaire dans une vallée ou dans un canyon sous-marin. Le modèle comprend des ondes simples (Riemann waves),le cas d'un barrage qui cède (dam-break problem), des solutions auto-similaires dans la M-fonne et « chapeaux paraboliques ». Grâce à la théorie linéaire de l'eau peu profonde nous avons étudié le processus de génération de tsunamis par des glissements de terrain de volume variable, se déplaçant à une vitesse variable dans un bassin de profondeur variable. Dans le cas d'un fond marin particulier (cas sans "réflexion"), les phénomènes de résonnance ont été étudiés dans un bassin à profondeur variable. Nous avons utilisé des méthodes numériques pour la résolution non-linéaire des équations des eaux peu profondes afin d'analyser des catastrophes marines réelles (les ondes de tempête causées par le cyclone Lili en 2002, le tsunami volcanique de 2003 à Montserrat) et probables (un tsunami prés des côtes de la Martinique). Des données sur les catastrophes / The dissertation is devoted to research in the field of marine natural hazards (tsunamis, storm surges) in the French West Indies, using analytical and numerical models of fluid mechanics. Emphasis is placed on the development of models of landslide motion and generated tsunami waves. The simple "solid block" model is used to evaluate the characteristics of pyroclastic flow Soufriere Hills volcano (Montserrat). The "fluid model" of a landslide (so called Savage-Hutter model) is studied analytically; within this model a new family of exact solutions that describe the motion of nonlinear gravity flow in a valley or underwater canyon is found: nonlinear Riemann wave, dan break problem, self-similar solutions (M - wave and parabolic cap). In the framework of the linear shallow water theory the process of generation of tsunami waves by landslides of variable volume moving with variable velocity above the basin of variable depth is studied. For the specific bottom profile ("reflectionless" one) the resonant phenomena is investigated in the basin of variable depth. Numerical methods are used to analyze marine hazards: historical (storm surges, caused by Cyclone LILI in 2002; volcanic tsunami 2003 on Montserrat) and possible events (possible tsunami of the coast of Martinique). Various data on marine natural disasters are obtained during field surveys (volcanic tsunami in 2003, stonn surges caused by Hurricane Dean in 2007). Designed catalogs of tsunamis and storm surges are created based on results of numerical modeling and field studies; some statistical analysis is perfomed

Le travail de cette thèse concerne l’étude des écoulements diphasiques dans les mini-canaux d’une pile à combustible. Ces canaux, de taille millimétrique, ont un double rôle d’alimentation de la pile en combustibles gazeux et d’évacuation de l’eau produite. Il a été montré expérimentalement que la configuration d’écoulement a un impact direct sur les performances de la pile. L’objectif est donc de progresser dans la compréhension des mécanismes physiques présents dans ces canaux. Pour cela, l’outil de simulation numérique JADIM est développé pour modéliser la physique de la ligne triple afin de permettre la description des écoulements diphasiques et les transitions entre configurations, où la mouillabilité joue un rôle déterminant. Deux phénomènes physiques sont essentiellement étudiés : la migration dans le canal de gouttelettes formées par l’introduction d’eau dans le canal, et l’évolution et la stabilité de la répartition spatiale des phases lors du remplissage progressif des canaux.