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.

Cette étude s’inscrit dans le cadre de la sûreté de l’installation ITER. En effet, des études d’accidentologie concernant cette installation ont permis de mettre à jour un risque d’explosion de mélanges à base d’hydrogène et de poussières. L’objectif de la thèse est d’acquérir les données fondamentales caractérisant l’explosion de ces mélanges diphasiques permettant ainsi d’évaluer les pressions générées par une éventuelle explosion des poussières qu’elle soit concomitante on non à celle de l’hydrogène. Pour se faire, des expériences, en bombe sphérique, ont été réalisées concernant des mélanges gazeux hydrogène - oxygène - azote. Les expériences ont été accomplies pour des températures de 303 et 343 K et des pressions de 50 et 100 kPa pour différentes concentrations en hydrogène et différents rapports N2/O2 dans le mélange. Les paramètres de caractérisation de l’explosion de ces mélanges ont été déterminés tels que la pression maximale de combustion (PMAX), l’indice de déflagration (KG ou KST), le temps de combustion (tC), la vitesse fondamentale de flamme (SL°) et la longueur de Markstein (LB). Une modélisation cinétique de la vitesse de flamme utilisant le code COSILAB a été réalisée, permettant de déduire la vitesse fondamentale de flamme sur la base de trois mécanismes cinétiques détaillés de la littérature. Celle-ci a permis l’évaluation de l’énergie d’activation globale sur la base du modèle cinétique présentant le meilleur accord avec l’expérience. De plus des calculs thermodynamiques à l’équilibre ont été réalisés afin de comparer les pressions maximales de combustion aux valeurs théoriques. Pour les mélanges diphasiques, un nouveau système d’introduction des poussières a été mis en place et des expériences caractérisant les paramètres d’explosion de ces mélanges ont été réalisées dans la bombe sphérique. Celles-ci ont permis de mettre en exergue le fait qu’une explosion de poussières, sous certaines conditions peut-être concomitante à une explosion d’hydrogène. / The context of the study is the safety of the ITER installation. Indeed, studies have shown that it exists a risk for two-phase mixtures of hydrogen and dust can explode and create a safety risk for the ITER installation. This aims to obtain the fundamental data which characterize the explosion of these mixtures and to evaluate the pressure loads they can generate. To do so, experiments in spherical bomb have been carried out for hydrogen - oxygen - nitrogen mixtures at two initial temperatures (303 and 343 K) and pressures (50 and 100 kPa) for different hydrogen concentrations and different N2/O2 ratios. Explosion parameters like maximum combustion pressures (PMAX), deflagration indexes (KG or KST), combustion times (tC), fundamental flame speeds (SL°) and Markstein lengths have been determined. A kinetic modelling of the flame speed, using the COSILAB software was performed based on three detailed kinetic models available in the literature and allowed the calculation of the global activation energy on the basis of the kinetic model which showed the best agreement with the experimental data. Moreover equilibrium calculations were achieved to compare PMAX to the theoretical values. For two-phase mixtures, a new introduction device was tested and set up and experiments characterizing the explosions parameters of the two-phase mixtures have been performed in the spherical bomb. They were able to stress out the fact that, under some circumstances, dust explosion can be concomitant to a hydrogen explosion.

Cette thèse porte sur l’étude des écoulements de Taylor (ou poche/bouchon) gazliquide en microcanal. Ces écoulements où les effets de tension de surface sont prépondérants ont été étudiées expérimentalement et numériquement pour des géométries rectangulaires avec divers rapports d’aspects. Une première partie expérimentale a consisté à caractériser la formation de bulles (taille, fréquence) en fonction des conditions opératoires, des propriétés des fluides (notamment à travers le nombre capillaire) et du mode de mise en contact des fluides. La dynamique de l’écoulement établi a par la suite été étudiée à l’aide du code JADIM. La simulation de ces écoulements dominés par la tension de surface a nécessité de lever les limitations liées à la prise en compte de la force capillaire. En effet des courants parasites numériques sont créés à proximité de l’interface lors de la simulation d’écoulements capillaires. Une méthode Level Set a été implémentée et comparée à la méthode Volume of Fluid d’origine en termes de courants parasites. Des simulations numériques 3D ont permis l’étude des effets du nombre capillaire et de la géométrie sur la dynamique des bulles de Taylor (vitesse, pression et formes de bulles). Les effets inertiels souvent négligés ont été considérés et leur influence, notamment sur les sauts de pression à l’interface, a été mise en évidence. Le mélange dans le bouchon liquide a également été étudié. / This thesis focuses on the hydrodynamics of gas-liquid Taylor flow (or slug flow) in microchannels. These flows, which are generally dominated by surface tension forces, have been investigated in rectangular channels of various cross-sectional aspect ratios by means of both experimental visualizations and numerical simulations. The first experimental part aims at characterizing the bubble generation process (bubble length and frequency of break-up) depending on the operating conditions, the fluid properties, as well as the junction where both fluids merge. Numerical simulations of fully developed Taylor flow have been carried out with the JADIM code. The computation of such surface tension dominated flows requires an accurate calculation of the surface tension force. Some limitations of the Volume of Fluid method have been highlighted and a Level Set method has been developed in order to improve the calculation of capillary effects. Both methods have been compared in detail in terms of spurious currents. 3D numerical simulations have been performed and the influence of the capillary number, as well as the effects of geometry have been highlighted. Inertial effects have been taken into account and their influence on the pressure drop has been shown to be non-negligible. Mixing in the liquid slug has also been studied.

Le modèle multi-fluide permet de décrire par une approche Eulérienne les sprays polydispersés et apparaît donc comme une méthode indiquée pour les applications de combustion diphasique. Sa pertinence pour la simulation à l’échelle d’applications industrielles est évaluée dans ce travail, par sa mise en oeuvre dans des configurations bi-dimensionnelle et tri-dimensionnelle plus représentatives de ce type de simulations. Cette évaluation couple une étude de faisabilité en terme de coût de calcul avec une analyse de la précision obtenue, par des comparaisons avec les résultats de méthodes de références pour la description des sprays. Afin de définir une telle référence, une hiérarchisation des modèles de spray est proposée dans ce travail, soulignant les niveaux de modélisation associée aux diverses méthodes. Une première configuration d’écoulements tourbillonnaires est utilisée pour caractériser la méthode multi-fluide. L’étude de la structure mathématique du système de lois de conservation permet d’analyser la formation de singularités et de fournir les outils permettant d’évaluer leur impact sur la modélisation. Cette étude permet également de dériver un schéma numérique robuste et efficace pour des configurations bi- et tri-dimensionnelle. La description des dynamiques de gouttes conditionnées par la taille est évaluée dans ces configurations tourbillonnaires au moyen de comparaisons quantitatives, sur des champs instantanés, où le multi-fluide est confronté à une méthode Lagrangienne, ainsi qu’à des résultats expérimentaux. Afin d’évaluer le comportement de la méthode multi-fluide dans des configurations plus représentatives des problématiques industrielles, le solveur MUSES3D est développé, permettant, entre autres, une évaluation fine des méthodes de résolution des sprays. Une implémentation originale de la méthode multi-fluide, conciliant généricité et efficacité pour le calcul parallèle, est réalisée. Le couplage de ce solveur avec le code ASPHODELE, développé au CORIA, permet d’effectuer une évaluation opérationnelle des approches Euler/Lagrange et Euler/Euler pour la description des écoulements diphasiques à inclusions dispersées. Finalement, le comportement de la méthode multi-fluide dans des jets bi-dimensionnels et dans une turbulence homogène isotrope tri-dimensionnelle permet de montrer sa précision pour la description de la dynamique de sprays évaporant dans des configurations plus complexes. La résolution de la polydispersion du spray permet de décrire précisément la fraction massique de combustible en phase vapeur, un élément clé pour les applications de combustion. De plus, l’efficacité du calcul parallèle par décomposition de domaine avec la méthode multi-fluide permet d’envisager son utilisation à l’échelle d’applications industrielles. / The multi-fluid model, providing a Eulerian description of polydisperse sprays, appears as an interesting method for two-phase combustion applications. Its relevance as a numerical tool for industrial device simulations is evaluated in this work. This evaluation assesses the feasibility of multi-fluid simulations in terms of computational cost and analyzes their precision through comparisons with reference methods for spray resolution. In order to define such a reference, the link between the available methods for spray resolution is provided, highlighting their corresponding level of modeling. A first framework of 2-D vortical flows is used to assess the mathematical structure of the multi-fluid model governing system of equations. The link between the mathematical peculiarities and the physical modeling is provided, and a robust numerical scheme efficient for 2-D/3-D configurations is designed. This framework is also used to evaluate the multi-fluid description of evaporating spray sizeconditioned dynamics through quantitative, time-resolved, comparisons with a Lagrangian reference and with experimental data. In order to assess the multi-fluid efficiency in configurations more representative of industrial devices, a numerical solver is designed, providing a framework devoted to spray method evaluation. An original implementation of the multifluid method, combining genericity and efficiency in a parallel framework, is achieved. The coupling with a Eulerian/Lagrangian solver for dispersed two-phase flows, developed at CORIA, is conducted. It allows a precise evaluation of Euler/Lagrange versus Euler/Euler approaches, in terms of precision and computational cost. Finally, the behavior of the multi-fluid model is assessed in 2D-jets and 3-D Homogeneous Isotropic Turbulence. It illustrates the ability of the method to capture evaporating spray dynamics in more complex configurations. The method is shown to describe accurately the fuel vapor mass fraction, a key issue for combustion applications. Furthermore, the method is shown to be efficient in domain decomposition parallel computing framework, a key issue for simulations at the scale of industrial devices.

Le calcul d'écoulements dans les générateurs de vapeur des réacteurs à eau pressurisée est un problème complexe, faisant intervenir différents régimes d'écoulement et plusieurs échelles de temps et d'espace. Un scénario accidentel peut être caractérisé par des variations très rapides pour un nombre de Mach de l'ordre de l'unité. A l'inverse en régime nominal l'écoulement peut être stationnaire, à bas nombre de Mach. De plus quelque soit le régime considéré, la complexité de la géométrie d'un générateur de vapeur conduit à modéliser le faisceau de tubes par un milieu poreux, d'où le problème de couplage à l'interface avec le milieu fluide.Un schéma de correction de pression tout-Mach en volumes finis colocalisés a été introduit pour les équations d'Euler et de Navier-Stokes. L'existence d'une solution discrète, la consistance du schéma au sens de Lax et la positivité de l'énergie interne ont été démontrées. Le schéma a été ensuite étendu aux modèles diphasiques homogènes du code GENEPI développé au CEA. Enfin un algorithme Multigrille-AMR a été adaptée pour permettre de mettre en oeuvre notre schéma sur des maillages adaptatifs.Concernant la seconde problématique, une extension de la loi de Beavers-Joseph a été proposée pour le régime convectif. En introduisant un saut d'énergie cinétique à l'interface, on retrouve une loi de type Beavers-Joseph mais avec un coefficient de glissement non-linéaire, qui dépend de la vitesse fluide à l'interface et de la vitesse Darcy. La validité de cette nouvelle condition d'interface a été évaluée en réalisant des calculs de simulation numérique directe à différents nombres de Reynolds. / The numerical simulation of steam generators of pressurized water reactors is a complex problem, involving different flow regimes and a wide range of length and time scales. An accidental scenario may be associated with very fast variations of the flow with an important Mach number. In contrast in the nominal regime the flow may be stationary, at low Mach number. Moreover whatever the regime under consideration, the array of U-tubes is modelled by a porous medium in order to avoid taking into account the complex geometry of the steam generator, which entails the issue of the coupling conditions at the interface with the free-fluid.We propose a new pressure-correction scheme for cell-centered finite volumes for solving the compressible Navier-Stokes and Euler equations at all Mach number. The existence of a discrete solution, the consistency of the scheme in the Lax sense and the positivity of the internal energy were proved. Then the scheme was extended to the homogeneous two-phase flow models of the GENEPI code developed at CEA. Lastly a multigrid-AMR algorithm was adapted for using our pressure-correction scheme on adaptive grids.Regarding the second issue addressed in this work, an extension to the Beavers-Joseph law was proposed for the convective regime. By introducing a jump in the kinetic energy at the interface, we recover an interface condition close to the Beavers-Joseph law but with a non-linear slip coefficient, which depends on the free-fluid velocity at the interface and on the Darcy velocity. The validity of this new transmission condition was assessed with direct numerical simulations at different Reynolds numbers.

De nombreuses hypothèses sont classiquement utilisées pour décrire le comportement du fluide dans un contact lubrifié : film continu, viscosité constante dans l’épaisseur, film mince, fluide newtonien… Or, certaines s’avèrent erronées dès lors que l’on s’intéresse aux contacts Elasto- HydroDynamiques fortement glissants ou à la répartition du lubrifiant en sortie de contact. Une approche numérique originale, basée sur un retour aux équations de la mécanique des fluides générale et prenant en compte le couplage fluide/solide et les effets thermiques sont proposés ici dans le but d’apporter des éléments physiques supplémentaires aux modélisations usuelles. Dans un premier temps, l’influence des effets thermiques sur l’évolution du frottement dans les contacts Thermo-EHD est mise en évidence. La présence d’un minimum de frottement pour le cas du glissement pur est expliquée par l’analyse des transferts thermiques entre le lubrifiant et les solides. L’origine des modifications locales d’épaisseur de film observées et l’existence même d’une épaisseur de film lubrifiant pour les cas de vitesse d’entraînement nulle sont alors reliées à la présence d’un fort gradient de viscosité dans l’épaisseur de film. Une comparaison qualitative avec des éléments expérimentaux de la littérature est réalisée, validant les tendances obtenues. Dans un second temps, l’écoulement à surface libre du lubrifiant en périphérie du contact est étudié expérimentalement puis numériquement par une méthode à interface diffuse. Le rôle des effets capillaires est analysé et les résultats numériques confrontés à des résultats issus de la littérature. Un bon accord est obtenu tant qualitativement que quantitativement. Validé par l’étude numérique diphasique (air/lubrifiant) réalisée, un modèle analytique simplifié est alors développé, prédisant une loi de répartition du lubrifiant en sortie de contact. La zone de sortie des contacts EHD est ensuite traitée par un modèle de cavitation vaporeuse et la prise en compte nécessaire de l’air environnant est discutée. Enfin, une première modélisation tridimensionnelle de l’écoulement à surface libre du lubrifiant autour d’un contact ponctuel est réalisée mettant en avant l’influence des effets capillaires et la faisabilité d’une telle approche. / Classically, many assumptions are used to model the fluid behaviour in a lubricated contact : continuous film, constant viscosity across the film thickness, film thickness is very thin compared to other contact dimensions, Newtonian lubricant... However, some of them are not well-founded for the study of Elasto-HydroDynamic contacts with high sliding or to estimate the liquid distribution at the exit of the contact. An original numerical approach, based on the general fluid mechanics equations and taking into account the fluid/solid coupling and thermal effects, is developed here in order to give more physical insights to the usual modelling. First of all, the thermal effects are shown on the friction coefficient evolution for Thermo- EHD contacts. A minimum value is found concerning the friction value for the pure sliding case. It is explained by analyzing the heat transfer between the solids and the lubricant. The origin of the resulting local modifications of the film thickness and the existence of a film thickness for zero entrainment velocity cases are related to the presence of a high viscosity gradient through the film. A qualitative comparison is performed with experimental data from literature, validating the results. Second, the free surface flow of the lubricant around the contact is experimentally and numerically studied with a diffuse interface method. The capillary effects on the air/lubricant meniscus position are analyzed and quantitatively compared with experimental data from literature. Good agreements are found. An analytical approach is then developed, based on the numerical study of the two-phase flow. An analytical law predicting the liquid distribution is obtained. The exit area of EHD contacts is then investigated with a vaporous cavitation model highlighting the necessity of taking into account the effects of surrounding air and surface wettabillity. Finally, a first approach of the tri-dimensional two-phase flow is performed, showing the capillary effects on the interface location.

Ces travaux de thèse s'incrivent dans le cadre du traitement de gaz acides et captage CO2 dans les colonnes à garnissages structurés. Les gaz à traiter réagissent avec un liquide s'écoulant à contre-courant sur des plaques métalliques dont la compléxité géométrique permet d'accroître l'aire d'échange, et donc l'efficacité du procédé. Dans un contexte de modélisation multi-échelles des contacteurs à garnissages structurés, les écoulements gaz-liquide à la plus petite échelle géométrique des plaques de garnissages (de l'ordre de l'épaisseur du film liquide) sont étudiés, pour améliorer la compréhension et la modélisation des écoulements diphasiques et phénomènes de mouillage dans les garnissages. L'objectif final est de développer une méthodologie CFD pour reproduire des écoulements diphasiques 3D sur des géométries complexes telles que les plaques de garnissages. Pour ce faire, il est nécessaire de progresser en méthodes numériques et de proposer des méthodes expérimentales pour observer des écoulements de film liquide sur des géométries complexes. Ces travaux comprennent une partie numérique et une partie expérimentale. Un écoulement sur une plaque de garnissage structuré peut présenter des zones sèches, et donc des lignes de contact (dynamiques), ce qui présente un défi en simulation numérique à cause des différentes échelles de l'écoulement. La méthodologie employée ici en simulation numérique consiste à résoudre l'écoulement jusqu'à une échelle intermédiaire en modélisant les effets des plus petites échelles. Le code de calcul Two-Phase Level-Set a été utilisé et modifié dans ce but. Différentes méthodes level-set ont d'abord été testées de manière à identifier une méthode satisfaisante quant à la réduction des erreurs de conservation de masse, un problème rencontré en level-set. Il est ici montré que certaines combinaisons de schémas de discrétisation spatiale et temporelle permettent de réduire considérablement ces erreurs de conservation de masse. Après avoir réalisé de nombreux tests de validation, une nouvelle méthode numérique est proposée pour simuler les grandes échelles d'écoulements diphasiques 3D avec ligne de contact dynamique en level-set, dans des conditions réalistes. La méthode est ici validée pour des écoulements axisymétriques de gouttes simulés en 3D, en régime visqueux et en régime inertiel, et pour des écoulements de gouttes sur plan incliné. Les résultats sont en très bon accord avec d'autres travaux numériques et expérimentaux. Afin de faciliter l'utilisation de cette méthodologie pour des applications industrielles, un modèle sous-maille similaire a été implémenté dans un code VOF commercial; les résultats sont aussi en très bon accord avec d'autres travaux. En plus de ces développements numériques, une campagne expérimentale est mise en oeuvre pour observer des écoulements de film liquide sur une plaque de garnissage structuré. Les méthodes expérimentales employées sont d'abord testées et validées pour des écoulements de film plat ou ondulé sur plan incliné, et ensuite utilisées pour observer des écoulements de film sur des plaques de garnissage. L'épaisseur de film liquide est mesurée aux creux et aux crêtes des picots des plaques de garnissages, pour différents débits, par imagerie confocale chromatique. Des lois de puissance de l'épaisseur de film en fonction du Reynolds sont proposées; celles-ci sont très différentes suivant la position des relevés de mesure, aux creux ou aux crêtes des picots. La vitesse à l'interface de l'écoulement gaz-liquide est aussi mesurée, par PIV et PTV, en utilisant des particules hydrophobes. Les résultats montrent que le liquide a tendance à dévier du creux des canaux (corrugations), et la norme de la vitesse semblent présenter des extremums correspondant aux creux et crêtes des picots. [...] / The work described in this thesis is motivated by the use of structured packing columns in acid gas treatment and post-combustion CO2 capture. In a counter-current mode, flue gases react with the liquid that flows down over metal sheets, the geometrical complexity of which allows increasing the specific interfacial area, and thereby the overall efficiency of the process. In the context of multiscale modeling of structured-packing contacting devices, the focus in this work is on the gas-liquid flows at the smallest geometrical scale of packing sheets, of the order of the liquid film thickness, aiming to improve understanding and modeling of two-phase flows and wetting phenomena in structured packings. The ultimate objective is to build up a CFD methodology to reproduce 3D two-phase flows over complex surfaces such as structured packing sheets. For this purpose, progress is necessary both in pertinent computational methods and in the adaptation of experimental methods for observing liquid film flows over complex surfaces. This thesis therefore consists of computational and experimental parts. Flows over structured packing sheets may exhibit dry zones, and hence (moving) contact lines, the numerical simulation of which presents a computational challenge due to the disparity in length scales involved. Here, the methodology for large-scale numerical simulations of flows with moving contact lines consists in resolving the flow down to an intermediate scale and modeling effects of smaller ones. The parallelized freeware Two-Phase Level-Set has been extended for this purpose. First though, because some level-set methods have been reproached to yield mass conservation issues, an assessment is made of the mass conservation properties of a range of level-set methods. It is demonstrated that the combined use of some spatial and temporal discretization schemes allows to drastically reduce mass conservation errors in level-set methods. Having thus implemented a level-set method with satisfactory performance at such tests (and others), a novel numerical method is proposed to perform 3D large-scale simulations of flows with moving contact lines in level-set, under realistic conditions. Validation tests of axisymmetric droplet spreading in a viscous, and in an inertial regime, simulated in 3D, and sliding drops are shown to be in excellent agreement with prior experimental and numerical work. The results show that complex contact-line dynamics observed in prior experimental studies on sliding droplets can be simulated using the present large-scale methodology. To facilitate dissemination of this work in industrial applications, a similar subgrid model has been implemented in a commercial volume-of-fluid code; results of validation tests are shown to be in excellent agreement with other work. These computational developments are accompanied by an experimental campaign to observe liquid film flows over structured packing sheets. All experimental methods used herein are tested and validated for flat and wavy films down an inclined plane before being used for observing liquid film flows over packing sheets. The film thickness is measured at local troughs and crests of small-scale corrugations of the structured packing sheet, for different flow rates, by Chromatic Confocal Imaging. Power laws of the Reynolds number for the mean liquid film thickness are suggested, with significant differences for measurements at crests compared to that at troughs. Interface velocity measurements are also performed by PIV and PTV using hydrophobic particles. Results reveal that the liquid tends to deviate from troughs of large-scale corrugations, and seems to exhibit local extrema of the velocity magnitude corresponding to troughs and crests of small-scale corrugations. [...]

L’effet thermocapillaire (ou Marangoni) est la résultante mécanique d’un gradient de tensioninterfaciale induit par la présence d’un gradient de température sur une interface fluide. Il semanifeste par (i) la migration d’un objet fini (goutte, bulle) immergé, et (ii) une déflexion del’interface. Sa nature interfaciale le rend particulièrement pertinent à petite échelle, notammenten microfluidique diphasique. Ce travail de thèse montre comment un effet thermocapillaireinduit localement par chauffage laser peut être utilisé pour produire des composants optofluidiquesélémentaires (vanne, aiguillage, échantillonneur), et en présente une étude quantitative.La déstabilisation d’un jet microfluidique forcée par laser, conduisant à sa rupture, est égalementprésentée et caractérisée. Cette « boîte à outils » optique fournit ainsi une approche sans contact,pour produire et manipuler des gouttes en microfluidique digitale sans nécessité d’une microfabricationdédiée. Par ailleurs, afin de caractériser sur des temps longs les gouttes produites,et ainsi considérer des populations statistiquement significatives, un dispositif optoélectroniquesimple pour mesurer les gouttes et leur vitesse en temps réel a également été développé. / The thermocapillary (or Marangoni) effect is the mechanical result of an interfacial tension gradientinduced by a temperature gradient on a fluid interface. This effect manifests itself byinducing (i) the migration of an immersed finite-size object (droplet, bubble), and (ii) a deflexionof the interface. Due to its interfacial nature, the Marangoni effect is particularly relevantat small length scales, especially in the context of two-phase microfluidics. This thesis aims atapplying the thermocapillary effect locally induced by laser heating, in order to create some basicoptofluidic actuators (valve, switch, sampler). A quantitative study of these actuators is presented.The laser-forced destabilization of a co-flowing microfluidic jet, leading to its breakup,is also investigated. This “optical toolbox” represents a non-contacting, and microfabricationfreeapproach for the production and handling of droplets in digital microfluidics. Moreover, tocharacterize these droplet over long times, thus considering statistically significant populations,a simple optoelectronic device has been developed for measuring the size and velocity of thedroplets in real time.

Ces travaux portent sur la modélisation et la simulation numérique des effets du déplacement d'un corps indéformable dans un écoulement multiphasique compressible. Ils se placent dans le cas où plusieurs objets sont en mouvement ou dans le cas où un objet est en mouvement dans un milieu aux géométries complexes. L'étude ne peut alors pas être placée dans le référentiel lié à l'objet en mouvement. Le modèle est basé sur une méthode multiphasique à interfaces diffuses où les différentes phases sont en équilibre mécanique. Le système régissant l'écoulement fluide est augmenté d'une équation d'advection. Cette dernière s'applique sur une fonction Level Set dont le niveau zéro permet de localiser le mobile dans l'espace. Des termes de couplage sont ajoutés au membre de droite des équations d'évolution de la quantité de mouvement et de l'énergie totale. Ces termes sont composés d'un facteur du type pénalisation et d'un facteur du type relaxation de vitesses. Cette nouvelle méthode permet de simuler des cas complexes où peuvent interagir des mobiles à hautes vitesses, des ondes de choc et des interfaces liquide/gaz. / This work deals with modelling and the numerical simulation of the effects of a moving rigid body on a multiphase flow. Here more than one object is moving, or an object is moving in a complex geometry domain. So the reference frame linked to the moving body can't be used. The model is build on a multiphase diffuse interface method with mechanical equilibrium. An advection equation is added. It applies on a Level Set function used to track the moving body. Coupling terms are added to the momentum equation and to the total energy equation. These terms are made of a penalization factor and a velocity relaxation factor. This new method allows to simulate complex cases where can interact high velocity objects, shock waves and liquid / gas interfaces.

Parmi les nombreuses problématiques liées à la PEMFC, l'eau est un point crucial car elle est à la fois indispensable et indésirable dans le coeur de pile. En effet, les flux gazeux qui alimentent le système en réactifs sont humidifiés pour que la membrane soit convenablement hydratée afin de faciliter le transport des protons de l'anode vers la cathode. Mais à cet apport initial vient s'ajouter l'eau produite par la réaction électrochimique. Des gouttes d'eau liquide peuvent alors se former dans les pores des couches actives et peuvent grossir jusqu'à boucher les canaux des plaques d'alimentation en gaz. Ce travail de thèse est centré sur l'étude des écoulements gaz-liquide dans ces canaux, dont le diamètre hydraulique est de l'ordre du millimètre. Dans un premier temps, des micro-capteurs ont été mis au point pour détecter la présence d'eau liquide. Après une phase de tests en maquette froide, ces micro-capteurs ont été insérés dans une cellule de 25 cm2 de surface active, identique à celles étudiées au LRGP. Ceci a permis de montrer que dans une pile en fonctionnement, les fluctuations de tension sont liées à la présence d'eau liquide. Par ailleurs, des corrélations existantes ont été utilisées pour prédire les pertes de charge engendrées par la présence de deux phases dans les mini-canaux. Ces calculs ont été comparés aux mesures expérimentales, mettant en évidence l'influence de la couche de diffusion des gaz, située entre les canaux et la membrane, et l'importance des conditions opératoires.Ces méthodes originales viennent compléter les outils de diagnostic existants sans demander d'appareillage complexe et coûteux, ni nécessiter le développement de matériel spécifique, non optimisé pour une utilisation commerciale / Water management is a critical issue in the operation of proton exchange fuel cells (PEMFC). On one hand, the membrane must be sufficiently hydrated to allow the transport of protons from the anode side to the cathode side. But on the other hand, supersaturation of water vapour in the gas phase results in liquid water formation, which can be detrimental to the fuel cell operation. Liquid water can clog up the porous structure of the gas diffusion layer (GDL) and hinder the transport of gases in the catalyst layer. Liquid accumulation can also lead to the formation of water columns inside the gas flow channels, thus preventing the fuel gas from flowing into the reaction area.Micro-sensors have been designed to detect the presence of liquid water. Three of them were inserted inside the gas distribution channels of a 25 cm2 cell. The tests run with this equipped cell highlighted the link between the presence of liquid water and the variations of the cell voltage. Besides, existing correlations were used to calculate the pressure drop between the inlet and the outlet of the cathode channel due to two-phase flow during the cell operation. The results of this calculation were compared with experimental results. This part of the study highlighted the influence of the GDL and of the operating conditions on the flows inside the channels