Le travail présenté concerne l'étude des instabilités thermo-acoustiques apparaissant dans une chambre de combustion aéronautique. Le montage expérimental permet de faire varier continument les fréquences de résonances et de passer pour une même condition d'écoulement, d'un régime de combustion stable à un régime instable. La caractérisation complète d'un cas instable comprend une mesure des conditions acoustiques aux limites du banc, une analyse du comportement des phases liquide et évaporée, de celui du dégagement de chaleur instationnaire et une mesure de la fonction de transfert de la flamme. Ces travaux ont mis en évidence trois phénomènes jouant sur l’injection cyclique de carburant liquide. Les temps caractéristiques des différents phénomènes intervenant dans le couplage thermo-acoustique et une distribution spatiale de l'indice de Rayleigh sont présentés. / The purpose of this experimental study was to further our understanding of the fuel spray behavior during combustion instability phenomena in combustion chambers. An aeronautical injection system with dual kerosene lines was mounted on the LOTAR setup, which was equipped with an adjustable exhaust length. Stability maps were generated by varying the global equivalent ratio and the fuel split parameter, for two Inner Exhaust Lengths (IEL). A non-unique multiphase flow condition was found to produce stable and unstable combustion for different IELs. Each configuration was fully characterized. Acoustic boundary conditions were measured using the 2-microphone technique. Different optical techniques were used toanalyze the unsteady behavior of the liquid phase, fuel vapor, and heat release. Moreover, two techniques were exploited to study the Flame Transfer Function using velocity measurement supstream and downstream of the injection device. Altogether, these results highlighted three atomization phenomena occurring during the cycles of thermo-acoustic instabilities. The phase-averaged analysis applied on the different measurements permitted to determine thetime scales associated with each process appearing in the thermo acoustic coupling. This cyclic injection of liquid fuel into the chamber was followed by a vapor phase increase corresponding to a wave of equivalent ratio. The delay between the two phenomena was of10°. In addition, OH* emissions showed a cyclic behavior following these waves. The delay between the wave of equivalent ratio and the unsteady heat release was approximately of 25°.Finally, spatial distribution of the Rayleigh index revealed that the inner recirculation zone contributed to sustain the combustion instability.

Sous la direction de SNECMA, un nouveau moteur cryotechnique pour Ariane 5 est en phase de développement. Ce moteur, VINCI, fonctionne à l'oxygène et à l'hydrogène liquides et devra être rallumable en cours de vol. Le transitoire de remplissage d'une cavité intermédiaire par laquelle transite l'oxygène liquide avant d'arriver dans la chambre de combustion s'avère être une étape critique qu'il faut s'efforcer de maîtriser. Cette cavité, appelée dôme LOX, est directement reliée à la chambre de combustion par les injecteurs à oxygène. Des outils numériques sont actuellement en cours de conception : ils permettront à terme de simuler le remplissage de cette cavité dans l'espace. Afin de valider ces outils numériques, un programme expérimental a vu le jour, impliquant SNECMA, le CNES (Centre National d'Etudes Spatiales) et le LEGI (Laboratoire des Ecoulements Géophysiques et Industriels). Il s'agit de mener à bien des expériences « simples » et reproductibles afin de disposer d'une base de données expérimentales qui servira de cas tests pour les simulations. Un banc d'essais expérimental a été progressivement mis en place et instrumenté au LEGI. Le fluide utilisé en substitution de l'oxygène liquide est de l'eau. On a choisi de respecter une similitude du nombre de Weber entre le cas « réel » et les conditions expérimentales. Ce dernier équivaut au rapport des forces d'inertie sur les forces de tension de surface. Deux campagnes d'essais ont été réalisées, dans lesquelles on s'est centré sur l'étude des aspects hydrodynamique du transitoire de remplissage : variation du débit total et des pressions au cours d'un essai, évaluation du taux de vide dans la cavité, visualisation de l'écoulement dans la cavité et en sortie des injecteurs. L'instrumentation à disposition est constituée d'un débitmètre à effet Coriolis, de capteurs de pression, d'une sonde optique, de caméras rapides et d'un laser pour l'imagerie. La première campagne d'essais a visé le remplissage de la cavité en eau seule. Le paramètre de contrôle principal était la pression génératrice de l'écoulement liquide. Dans la deuxième campagne d'essais on injecte simultanément dans la cavité un écoulement d'eau et un écoulement d'air. Cela se rapproche plus des conditions réelles du transitoire de remplissage, au cours duquel la cavité est balayée par un écoulement d'hélium. Le paramètre de contrôle supplémentaire est le débit de gaz initial. Ces campagnes ont également souligné l'importance du profil d'ouverture de la vanne de l'écoulement liquide sur le transitoire de remplissage de la cavité. Ces campagnes d'essais constituent une première étape dans la compréhension du transitoire de remplissage du dôme LOX. Elles ont permis de visualiser la forme de l'écoulement dans la cavité et en sortie des injecteurs et d'identifier certains phénomènes intervenant dans le remplissage de la cavité. En particulier, nous avons mis en évidence l'existence d'un délai de mise en place de l'écoulement par les injecteurs, qui peut être responsable d'un pic de pression dans la cavité au cours du transitoire. L'influence de la fraction gazeuse sur l'écoulement dans les injecteurs a été soulignée mais reste à quantifier de façon précise. La prochaine étape de l'étude concerne les aspects énergétiques du transitoire de remplissage, notamment les transferts thermiques ayant lieu, dans la réalité, entre l'hélium, l'oxygène et les parois du dôme LOX. / Under the supervision of SNECMA, a new cryotechnic engine is being developed for Ariane 5. This engine, named VINCI, uses liquid oxygen and liquid hydrogen as propellant. It must be re-ignitable in flight. The filling transitory phase of an intermediate tank where the liquid oxygen passes through before entering the combustion chamber, has proved to be a very important stage that must be handled. This tank, called LOX dome, is directly linked to the combustion chamber through the oxygen injectors. Numerical tools are currently under development. They will allow to simulate the filling of this tank in the space. In order to validate those numerical tools, an experimental program has been launched. It involves SNECMA, the CNES (Centre National d'Etudes Spatiales : National Centre for Spatial Studies) and the LEGI (Laboratoire des Ecoulements Géophysiques et Industriels : Laboratory of Geophysical and Industrial Flows). Simple and repeatable experiments must be run. They will allow to gather experimental data that will further be used as test cases for the simulations. A test bench has been brought into service step by step at the LEGI, as well as scientific instruments. Water is used in place of liquid oxygen. A similarity of flows based on the Weber number has been chosen between the real case and the experiment. The Weber number measures the relative importance of the fluid inertia compared to its surface tension. Two experimental campaigns have been realized, that have focused on the dynamic aspects of the filling transitory phase : variations of the total flow and of the pressures measured during an experiment, evaluation of the void fraction in the tank, flow visualization in the tank and at the outlet of the injectors. The scientific instrumentation used is made of a Coriolis flow-meter, pressure probes, an optical probe, and high speed cameras with a laser for the flow visualization. The first experimental campaign has studied the tank filling with water only. The main control parameter is the reference pressure of the liquid flow. In the second campaign, both liquid and air flows are simultaneously injected in the tank. It aims at reproducing the real conditions of the filling transitory phase, where helium is injected in the tank with the liquid oxygen. The additional control parameter is the initial gas flow. Those campaigns have shown as well the importance of the valve opening that controls the liquid flow. Those campaigns are a first step in the understanding of the filling transitory phase of the LOX dome. They have permitted to visualize the flow in the tank and at the outlet of the injectors and to point out some important phenomena occurring during the tank filling. In particular, they have highlighted the existence of a delay before the flow can develop through the injectors. This delay can be responsible for a pressure peak in the tank during the transitory phase. The influence of the gas fraction on the flow through the injectors has been underlined as well but still must be accurately quantified. The next step of the study concerns the energetics of the filling transitory phase, especially the thermal transfers that occur between the helium, the oxygen and the walls of the tank.

Le transport sédimentaire contrôle l'évolution morphologique des rivières, l'érosion du littoral et l'équilibre des écosystèmes. Il constitue également un facteur de risque pour les populations et les infrastructures. Le sheet-flow, ou charriage intense, est un régime de transport sédimentaire qui s'établi lors de crues dans les fleuves et les rivières ou lors du déferlement des vagues littorales sur les plages sableuses. Le fort taux de transport associé à ce régime le rend très morphogène et une bonne compréhension des processus physiques impliqués est fondamentale pour prédire la morphodynamique. Cependant, les interactions granulaires et les fluctuations turbulentes, qui sont les principaux mécanismes à l'œuvre dans ce phénomène, constituent des verrous scientifiques pour la modélisation du régime de sheet-flow. Cette déficience s'explique essentiellement par le manque de données expérimentales haute résolution. Partant de ce constat, l'objectif de la thèse est de proposer un modèle diphasique et un dispositif expérimental haute résolution permettant de mieux caractériser les mécanismes impliqués. indent Dans un premier temps, le modèle diphasique est présenté et les résultats obtenus sont confrontés aux données de la littérature. L'analyse des résultats montre que la rhéologie des écoulements granulaires denses ($mu(I)/phi(I)$) et l'approche de longueur de mélange utilisées sont des fermetures appropriées pour reproduire les principales caractéristiques du régime de sheet-flow pour une large gamme d'écoulements et de propriétés sédimentaires. La deuxième partie de la thèse est consacrée à la mise en place d'un dispositif expérimental capable de fournir des mesures instantanées de vitesse et de concentration en régime de sheet-flow uniforme. Dans la troisième partie les grandeurs moyennes sont analysées pour décrire la structure verticale de l'écoulement. Les résultats obtenus montrent qu'une formulation en longueur de mélange et un profil de Rouse permettent de décrire la contrainte turbulente et le profil de concentration dans la suspension à condition de fortement modifier le paramètre de von Karman ($kappa approx 0.2$) et le nombre de Schmidt ($sigma_s=0.44$). La rhéologie frictionnelle ($mu(I)/phi(I)$) et la théorie cinétique des écoulements granulaires prédisent qualitativement le comportement observé, mais échoue à reproduire quantitativement les mesures. Le lien étroit existant entre les structures cohérentes turbulentes et la dynamique du lit sédimentaire illustre l'importance des fluctuations et de l'intermittence de l'écoulement. Ce couplage pourrait expliquer l'écart observé entre le comportement prédit par les modèles de contraintes inter-granulaires et les mesures expérimentales. Finalement, la comparaison des analyses statistiques en régime de sheet-flow et en écoulement sur fond fixe rugueux permet de montrer que l'énergie cinétique turbulente est peu affectée par la présence des sédiments mais que le niveau de corrélation entre fluctuations horizontales et verticales est sensiblement diminué, impliquant une diminution de la longueur de mélange et de la viscosité turbulente. Une augmentation significative de la rugosité équivalente induite par le lit mobile est aussi observée. / Sediment transport controls river morphological evolution, coastal erosion and ecosystem equilibrium. It represents a risk factor for populations and infrastructures. The sheet-flow, or intense bed-load, is a regime of sediment transport occurring during river floods or in the coastal wave breaking region above sandy beaches. The large amount of sediment transported in this regime is the main source for morphological evolution in our natural systemswater bodies. A good understanding of the underlying physical processes is a pre-requisite for accurate morphodynamic predictions. However, particle-particle interactions and turbulent flow interactions, which are the main driving mechanisms in this problem, constitute the scientific bottlenecks for sheet-flow modelling. This deficiency is mainly caused by the lack of high resolution experimental data. Based on this observation, the objective of the present thesis is to propose a novel two-phase model and to generate a new set of high resolution experiment data to improve process based sheet-flow modelling. indent First, the two-phase flow model is presented and the obtained results are compared with data from the literature. The result analysis has shown that the dense granular flow rheology ($mu(I)/phi(I)$) combined with a turbulent mixing length concept predicts the main sheet flow characteristics over a wide range of flow and sediment properties. Secondly, the experimental set up providing high-rate measurements of velocity and concentration under a uniform sheet-flow regime is presented. Third, the measured mean flow quantities are analysed to describe the vertical structure of the flow. The obtained results show that a mixing length formulation and a Rouse profile allow to describe the turbulent stress and the concentration profiles in the turbulent suspension layer, provided that the von Karman parameter and the Schmidt number are modified ($kappa approx 0.2$ and $sigma_s=0.44$). The frictional rheology ($mu(I)/phi(I)$) and the kinetic theory of granular flows predict qualitatively the observed behaviour but fail to reproduce measurements quantitatively. The observed link between the turbulent coherent structures and the bed dynamic illustrates the importance of flow fluctuations and intermittency. This coupling could be responsible for the discrepancy found between the predictions from the intergranular stresses models and the measurements. Finally, the comparison between the statistical analysis performed for a sheet-flow regime and for a clear water flow over a rough fixed bed demonstrates that the turbulent kinetic energy is weakly affected by the presence of sediments whereas the turbulent correlation level between horizontal and vertical fluctuations is significantly reduced, leading to a decrease of both the mixing length and the turbulent eddy viscosity. An increase of the equivalent roughness height induced by the moving bed is also observed.

Le travail présenté concerne le développement d’un cadre d’applications pour le traitement d’écoulements dispersés, tenant compte de l’effet de la coagulation sur la distribution des tailles des particules. Nous explorons également quelques techniques de validation expérimentale. Les modèles sont valables pour un écoulement incompressible et isotherme, avec des particules qui ont un temps de relaxation faible en comparaison à celui du fluide. Pour la phase dispersée, une méthode eulérienne est utilisée, ce qui permet d’extrapoler la vitesse des particules de celle du fluide. La distribution des tailles est modélisée à l’aide du « Direct Quadrature Method of Moments ». Cette approche permet de résoudre des équations de transport pour les poids et les abscisses d’une approximation de la distribution à l’aide des fonctions Dirac delta. L’effet de la coagulation est pris en compte à l’aide d’un noyau de collisions qui utilise la vitesse instantanée du fluide. Toutes les équations de transport sont résolues à l’aide de la méthode des éléments finis. Pour le fluide, les stabilisations « Streamline Upwind » et « Pressure Stabilized Petrov-Galerkin » sont utilisées ensemble avec une stabilisation grad-div. Afin de limiter le temps de calcul pour une simulation directe, une formulation utilisant un traitement explicite des termes d’advection est proposée. Avec l’apparition de gradients élevés, les équations de transport pour les particules nécessitent une stabilisation supplémentaire. Tout le travail est disponible dans le projet de logiciel libre Coolfluid 3, en utilisant un langage spécifique permettant une implémentation directe pour des modèles en éléments finis. Le code qui en résulte ressemble à la forme variationnelle des équations utilisées. Le programme est générique en termes de dimensions spatiales et de type d’éléments. Une première validation utilise des résultats trouvés dans la littérature comme référence. La précision des méthodes est vérifiée à l’aide des vortices Taylor-Green. Pour l’écoulement et les concentrations des particules, une simulation directe d’un canal turbulent est effectuée. Le noyau de coagulation est vérifié à l’aide de particules de différentes tailles qui tombent à travers un vortex Burgers. Finalement, quelques techniques de validation expérimentale sont utilisées dans une cellule d’essai. La technique « Particle Image Velocimetry » est utilisée pour les vitesses du fluide, tandis que la distribution des tailles est mesurée à l’aide du « Phase Doppler Anemometry » et « Multiple Wavelength Light Extinction ». La technique d’extinction de lumière est capable de produire des distributions des tailles qui peuvent être comparées facilement avec les résultats numériques. / The present work deals with the development of a framework for the modeling of dispersed flows, including the effect of coagulation on the particle size distribution. We also explore some techniques for experimental validation. Models are developed for incompressible, isothermal flow containing particles that have a small relaxation time compared to the fluid time scale. For the dispersed phase, an equilibrium Eulerian approach is used, extrapolating the particle velocity from the fluid velocity. The size distribution is modeled using the Direct Quadrature Method of Moments. In practice, this results in solving transport equations for the weights and abscissa of a Dirac delta approximation of the size distribution. To model the effect of coagulation, a collision kernel that makes use of the resolved instantaneous velocity is developed. All transport equations are solved using the Finite Element Method. For the fluid, the Streamline Upwind and Pressure Stabilized Petrov-Galerkin method are used, with additional grad-div stabilization. To decrease the solution time for DNS, a segregated formulation with an explicit advection term is proposed. The particle transport equations require cross-wind diffusion in addition to the streamline upwind stabilization when large gradients occur. All work is available in the open source Coolfluid 3 framework, using an Embedded Domain Specific Language we developed for the implementation of finite element models. The resulting code closely resembles the variational form of the equations and is generic in terms of element type and the number of spatial dimensions. A first validation uses literature results as reference. Correctness and accuracy of the methods are verified using the Taylor-Green vortex flow. For the fluid and particle concentration, direct numerical simulation of a turbulent channel flow is performed. The particle coagulation kernel is tested using particles of different sizes falling through a Burgers vortex. Finally, some experimental validation techniques are used on a small test chamber. Particle image velocimetry is used for the fluid motion, while the size distributions are measured using Phase Doppler Anemometry and Multiple Wavelength Light Extinction. The light extinction technique was found to produce size distributions that could provide valuable reference data for our particle model.

La compréhension des écoulements multiphasiques dans les évaporateurs à mini-canaux est primordiale pour la performance des boucles de climatisation dans le secteur automobile notamment. Cette thèse s’est principalement intéressée à l’écoulement d’entrée de tels évaporateurs ainsi qu’à la répartition des phases dans les mini-canaux. Dans un premier temps, l’écoulement adiabatique diphasique en entrée d’évaporateur a été étudié. Un dispositif expérimental transparent, respectant au mieux la géométrie d’entrée de l’évaporateur, a été réalisé afin de reproduire l’écoulement diphasique d’entrée en eau-air mais en respectant les régimes d’écoulement rencontrés avec du R134a. Plusieurs techniques de caractérisation ont été mises en œuvre (visualisation, conductimétrie, tube de Pitot et prises de pression) afin de quantifier les pertes de pression, les épaisseurs de film et les vitesses du gaz dans un régime principalement annulaire. Suivant le même principe, un autre module en acier-inox a été développé pour caractériser l’écoulement directement en entrée d’évaporateur avec du réfrigérant R134a. Dans un second temps, nous avons étendu l’étude au cas d’un évaporateur compact à mini-canaux. Dans deux situations adiabatiques : monophasique (eau) et diphasique (eau-air), les pertes de pression, la répartition des phases le long de l’évaporateur et le régime d’écoulement dans les mini-canaux ont été étudiés sur un échangeur fabriqué en polycarbonate dont la géométrie s’approche au mieux de celle d’un échangeur réel. Les nombreux résultats ainsi obtenus constituent une base de données conséquente utile à la simulation numérique de ce type d'écoulements diphasiques / Understanding of multiphase flows in mini-channel evaporators is essential for the performance of air-conditioning systems, particularly in automotive sector. This thesis is mainly interested in behavior of inlet flow and phase distribution in the mini-channels. Initially, an adiabatic two-phase flow at the evaporator's inlet was studied. A transparent experimental apparatus with the same geometry as an evaporator's inlet has been designed. This test section helped us to reproduce the same flow regimes with air-water as flow regimes encountered with R134a in an evaporator. Several characterization techniques were used (visualization, conductance probes, Pitot tube and pressure taps) to determine pressure losses, liquid film thickness and gas velocity in a predominantly annular flow regime. Following the same principle, another experimental facility in stainless steel was developed to directly characterize the R134a flow at the evaporator's inlet. Finally, we have extended the study to the case of a compact evaporator in two adiabatic situations: single-phase (water) and two-phase (air-water). Pressure losses, phase distribution along the evaporator and flow regime in mini-channels were studied on an evaporator made of transparent materials (polycarbonate) with a close geometry to that of a real evaporator. Numerous results were obtained to provide a consistent database that would be useful for numerical simulation of this type of two-phase flows

Ce travail de thèse présente l'implémentation de la simulation d'écoulements diphasiques dans des conditions de réacteurs nucléaires à caloporteur eau, à l'échelle de bulles individuelles. Pour ce faire, nous étudions plusieurs modèles d'écoulements thermohydrauliques et nous focalisons sur une technique de capture d'interface mince entre phases liquide et vapeur. Nous passons ainsi en revue quelques techniques possibles de maillage adaptatif (AMR) et nous fournissons des outils algorithmiques et informatiques adaptés à l'AMR par patchs dont l'objectif localement la précision dans des régions d'intérêt. Plus précisément, nous introduisons un algorithme de génération de patchs conçu dans l'optique du calcul parallèle équilibré. Cette approche nous permet de capturer finement des changements situés à l'interface, comme nous le montrons pour des cas tests d'advection ainsi que pour des modèles avec couplage hyperbolique-elliptique. Les calculs que nous présentons incluent également la simulation du système de Navier-Stokes incompressible qui modélise la déformation de l'interface entre deux fluides non-miscibles. / This PhD work presents the implementation of the simulation of two-phase flows in conditions of water-cooled nuclear reactors, at the scale of individual bubbles. To achieve that, we study several models for Thermal-Hydraulic flows and we focus on a technique for the capture of the thin interface between liquid and vapour phases. We thus review some possible techniques for Adaptive Mesh Refinement (AMR) and provide algorithmic and computational tools adapted to patch-based AMR, which aim is to locally improve the precision in regions of interest. More precisely, we introduce a patch-covering algorithm designed with balanced parallel computing in mind. This approach lets us finely capture changes located at the interface, as we show for advection test cases as well as for models with hyperbolic-elliptic coupling. The computations we present also include the simulation of the incompressible Navier-Stokes system, which models the shape changes of the interface between two non-miscible fluids.

Nous avons étudié des méthodes de volumes finis pour la simulation numérique d'un flux impliquant deux phases incompressibles ou deux phases générales compressibles en déséquilibre mécanique. Les principales difficultés du régime où il y a une apparition de phase ou une disparition de phase est la singularité de la vitesse. Nous montrons que l'utilisation du l'entropie correction améliorer beaucoup ces problèmes. Enfin, nous simulons certains tests numériques importants pour vérifier les méthodes numériques, telles que la séparation de phase par gravité ou un canal bouillant. / We investigated some finite volume methods for the numerical simulation of a flow involving two incompressible phases or general two compressible phases in mechanical disequilibrium. The main difficulties of the regime where there is either a phase appearance or a phase disappearance is the singularity of the velocity. We show that using the entropy fix will much improve these problems. Finally, we perfom some important numerical tests to verify the numerical methods, such as a phase separation by gravity or a boiling channel.

En géosciences, les applications associées au calage de modèles d'écoulement nécessitent d'appeler plusieurs fois un simulateur au cours d'un processus d'optimisation. Or, une seule simulation peut durer plusieurs heures et l'exécution d'une boucle complète de calage peut s'étendre sur plusieurs jours. Diminuer le temps de calcul global à l'aide des techniques de bases réduites (RB) constitue l’objectif de la thèse.Il s'agit plus précisément dans ce travail d'appliquer ces techniques aux écoulements incompressibles diphasiques eau-huile en milieu poreux. Ce modèle, bien que simplifié par rapport aux modèles utilisés dans l'industrie pétrolière, constitue déjà un défi du point de vue de la pertinence de la méthode RB du fait du couplage entre les différentes équations, de la forte hétérogénéité des données physiques, ainsi que du choix des schémas numériques de référence.Nous présentons d'abord le modèle considéré, le schéma volumes finis (VF) retenu pour l'approximation numérique, ainsi que différentes paramétrisations pertinentes en simulation de réservoir. Ensuite, après un bref rappel de la méthode RB, nous mettons en oeuvre la réduction du problème en pression à un instant donné en suivant deux démarches distinctes. La première consiste à interpréter la discrétisation VF comme une approximation de Ritz-Galerkine, ce qui permet de se ramener au cadre standard de la méthode RB mais n'est possible que sous certaines hypothèses restrictives. La seconde démarche lève ces restrictions en construisant le modèle réduit directement au niveau discret.Enfin, nous testons deux stratégies de réduction pour la collection en temps de pressions paramétrées par les variations de la saturation. La première considère le temps juste comme un paramètre supplémentaire. La seconde tente de mieux capturer la causalité temporelle en introduisant les trajectoires en temps paramétrées. / In geosciences, applications involving model calibration require a simulator to be called several times with an optimization process. However, a single simulation can take several hours and a complete calibration loop can extend over serval days. The objective of this thesis is to reduce the overall simulation time using reduced basis (RB) techniques.More specifically, this work is devoted to applying such techniques to incompressible two-phase water-oil flows in porous media. Despite its relative simplicity in comparison to other models used in the petroleum industry, this model is already a challenge from the standpoint of reduced order modeling. This is due to the coupling between its equations, the highly heterogeneous physical data, as well as the choice of reference numerical schemes.We first present the two-phase flow model, along with the finite volume (FV) scheme used for the discretization and relevant parameterizations in reservoir simulation. Then, after having recalled the RB method, we perform a reduction of the pressure equation at a fixed time step by two different approaches. In the first approach, we interpret the FV discretization as a Ritz-Galerkine approximation, which takes us back to the standard RB framework but which is possible only under severe assumptions. The second approach frees us of these restrictions by building the RB method directly at the discrete level.Finally, we deploy two strategies for reducing the collection in time of pressuresparameterized by the variations of the saturation. The first one simply considers time as an additional parameter. The second one attempts to better capture temporalcausality by introducing parameterized time-trajectories.

Nous présentons la généralisation d'un modèle de dépôt de particules aux écoulements non isothermes. Ce travail comporte deux parties. La première partie présente l'approche de modélisation lagrangienne stochastique adoptée pour simuler le dépôt de particules colloïdales et traite du raccord entre le Modèle de Langevin Généralisé qui calcule le transport des particules dans le cœur de l'écoulement et le Modèle de Proche Paroi destiné au transport des particules au voisinage de la paroi. Le raccord est assuré par l'équilibre du flux de particules traversant l'interface entre les deux modèles. La validation du couplage s'effectue par le passage en particules fluides, il s'agit de vérifier que le raccord est valable dans le cas limite du fluide. En effet, des accumulations non physiques de particules fluides peuvent se produire autour de l'interface. Le travail de couplage a donné lieu à des développements théoriques et numériques sur les statistiques conditionnées d'une particule traversant une interface. Des calculs de dépôt de particules montrent que le modèle reproduit correctement les résultats expérimentaux pour des écoulements isothermes. Cependant, dans un écoulement non isotherme, les particules peuvent faire l'expérience du phénomène de thermophorèse (qui est le mécanisme de déplacement des particules induit par un gradient de température) qui peut considérablement influencer le dépôt de particules. Nous avons donc choisi d'intégrer la modélisation de la thermophorèse dans le Modèle de Proche Paroi afin d'évaluer au mieux le taux de dépôt de particules dans un écoulement non isotherme. Nous commençons par étudier la thermophorèse en milieu gazeux car elle est très bien expliquée par la théorie cinétique des gaz. L'étude de la thermophorèse en milieu liquide est plus ardue car elle ne dispose pas encore de théorie unifié ni de travaux expérimentaux qui étudient ce phénomène. La dernière partie de ce document est donc consacrée à la prise en compte du phénomène de thermophorèse en liquide dans le Modèle de Proche Paroi. Cette étude fournit de nouveaux résultats numériques de dépôt de particules en présence d'un gradient de température dans un milieu liquide.

Cette thèse est consacrée aux écoulements homogènes de bulles, ainsi qu'à leur couplage avec le transport d'un scalaire et la turbulence. Elle s'intéresse plus spécifiquement aux effets de taille finie, des interactions hydrodynamiques et de la microstructure de la suspension qui sont étudiés à l'aide de simulations numériques directes à l'échelle d'une seule bulle. La dynamique d'une suspension laminaire de bulles induite par la seule gravité est d'abord revisitée. L'influence de la fraction volumique sur la vitesse de dérive des bulles est établie analytiquement et numériquement pour une suspension parfaitement ordonnée, puis des ressemblances entre suspensions ordonnées et suspensions désordonnées sont mises en évidence. Ces résultats sont ensuite mis à profit pour la modélisation du transport d'un scalaire passif au sein d'une suspension laminaire, tel que décrit par une diffusivité effective tensorielle, et des différences essentielles entre systèmes ordonnés et systèmes désordonnés concernant le transport de scalaire sont mises en exergue. Enfin, la turbulence est prise en compte dans les simulations et son interaction avec une bulle de taille finie est caractérisée. Il est montré que le comportement dynamique d'une bulle de taille comparable à la microéchelle de Taylor ressemble qualitativement à celui d'une microbulle, avec, notamment, une préférence pour certaines régions caractéristiques de l'écoulement. Une définition de l'écoulement vu par la bulle compatible avec les modèles standards de masse ajoutée et de portance est finalement proposée / This thesis is devoted to the study of homogeneous bubbly flows and their coupling with scalar transport and turbulence. It focuses on the effects of finite size, hydrodynamic interactions, and suspension microstructure, which are investigated using direct numerical simulations at the bubble scale. The dynamics of laminar buoyancy-driven bubbly suspensions is first revisited. More specifically, the effect of volume fraction on the bubble drift velocity is clarified by connecting numerical results to theory for dilute ordered systems, and similarities between perfectly ordered and free disordered suspensions are evidenced. These results are then used for the modeling of passive scalar transport in laminar suspensions as described by an effective diffusivity tensor, and crucial differences between ordered and disordered systems with respect to scalar transport are highlighted. Lastly, turbulence is included in the simulations, and its interaction with a finite-size bubble is characterized. The behavior of a bubble as large as Taylor microscale is shown to share a number of common features with that of a microbubble, most notably, the flow sampled by the bubble is biased. A definition of the liquid flow seen by the bubble, as it enters in usual models for the added mass and the lift forces, is finally proposed