Spelling suggestions: "subject:"écoulement diphasique,"" "subject:"ecoulement diphasique,""
71 |
Étude et modélisation de l'orientation de fibres dans des thermoplastiques renforcésMegally, Alexandra 13 July 2005 (has links) (PDF)
Cette thèse porte sur l'étude et la modélisation de l'orientation de fibres de verre dans des thermoplastiques renforcés. Dans une optique de maîtrise et de compréhension du comportement des composites renforcés de fibres longues, nous proposons dans ce travail d'étudier les phénomènes d'orientation et la rhéologie de ces suspensions. Une approche expérimentale a permis de caractériser l'orientation et la structure des fibres dans des pièces injectées. L'analyse de la distribution de longueur des fibres dans une pièce a mis en évidence que la casse des fibres se faisait préférentiellement durant la phase de la plastification de la matière. Les mesures d'orientation, réalisées avec des techniques de mesure 2D et 3D montrent que la concentration en fibres est un paramètre déterminant sur la formation et les caractéristiques de la structure coeur-peau. Une approche numérique de modélisation de l'orientation des fibres dans des pièces injectées, basée sur le calcul direct du mouvement d'orientation d'une population de fibres en écoulement, a permis d'accéder à la physique intrinsèque de la suspension. Ce modèle repose sur une formulation éléments finis multi- domaine développée dans le cadre du logiciel Rem3D. La principale originalité du modèle est qu'il n'est pas nécessaire d'exprimer de manière explicite l'ensemble des forces et interactions hydrodynamiques régissant le système. Lorsqu'on applique un cisaillement simple à un volume élémentaire représentatif, le calcul de simulation directe nous permet d'étudier les phénomènes d'interaction entre particules en fonction de la concentration et du rapport de forme des fibres. Ainsi, nous proposons une méthode d'identification du coefficient d'interaction Ci , paramètre déterminant du modèle statistique de Folgar & Tucker. Les valeurs de Ci sont fonction du rapport de forme et de la concentration des particules. On s'intéresse d'autre part à la validité des approximations de fermeture. On montre que l'approximation hybride est satisfaisante dans le cas d'un écoulement de cisaillement simple. Enfin, ce modèle est appliqué à l'étude rhéologique d'une suspension de fibres.
|
72 |
Analyse et approximation numérique de systèmes hyperboliques de lois de conservation avec termes sources. Application aux équations d'Euler et à un modèle simplifié d'écoulements diphasiques.Gosse, Laurent 19 September 1997 (has links) (PDF)
Cette thèse est consacrée à l'étude des systèmes non linéaires hyperboliques de lois de conservation avec des termes sources autorisées à devenir raides. Les applications viennent pour la plupart de problèmes physiques, par exemple l'equation de Burgers ou de Buckley-Leverett et le système d'Euler muni de lois de pression pour un gaz parfait ou une version simplifiée de mélange à deux phases. La première partie traite des aspects théoriques de ce sujet. Nous nous concentrons sur le cas unidimensionnel et considérons l'équation scalaire et d'un cas particulier d'un système d'écoulement diphasique réactif. La question principale à gérer lorsqu'on établit la convergence forte d'une suite d'approximations au moyen de la théorie de compacité par compensation dans de tels cas est de trouver des régions invariantes dans l'espace des phases. Pour le problème scalaire, il suffit que le terme source conserve une borne uniforme sur la norme L ∞ de la séquence. Mais c'est moins évident pour un système car les zones invariantes sont plutôt compliquées et peuvent être perturbées en raison de la forme de certains termes de forçage. La deuxième partie s'intéresse aux critères de stabilité pour les schémas numériques implicites impliquant des quadratures en temps sophistiquées. Nous présentons d'abord l'analyse d'une semi-discrétisation (méthode des lignes) afin d'en tirer des conditions suffisantes pour le système différentiel d'être bien posé. Cela amène des restrictions assez naturelles pour assurer la stabilité d'une méthode à une seule étape implicite: le pas de temps doit être petit près des points de compression et de répulsion du terme source. Une analyse de variation totale est ensuite réalisée pour dériver les régimes du second ordre sans oscillations autour de chocs pour le plus grand pas de temps possible. Des expériences numériques sont présentés sur des problèmes scalaires non linéaires. La troisième partie est une suite directe afin de corriger certaines des difficultés rencontrées auparavant. L'objectif principal est de construire une nouvelle discrétisation qui peut gérer un terme source arbitrairement raide sans aucune influence (sauf la condition habituelle CFL) sur le pas de temps. Suivant les idées de J. Greenberg et A.-Y. LeRoux, nous calculons d'abord un schéma de Godunov scalaire afin d'obtenir un résultat de convergence vers la solution entropique à l'aide des estimations BV. Comme la caractéristique principale est de gérer les sources au moyen des produits non-conservatifs, nous décidons d'utiliser le formalisme proposé par G. Dalmaso, PG LeFloch et F. Murat pour étendre ces idées à des systèmes non-homogènes. Suivant I. Toumi, nous introduisons des matrices de type Roe le long des chemins régularisants afin d'obtenir des approximations qui montrent des propriétés intéressantes. En fait, c'est ce que l'on appelle les régimes bien équilibrés, à savoir les régimes qui préservent les intégrales premières théoriques du mouvement. Des tests numériques sont basés sur le système d'Euler avec termes sources géométriques. Le quatrième (et dernière) partie concerne la mise en place de cette formulation non-conservative dans un schéma de type flux-splitting. L'objectif est de construire une méthode très robuste, bien équilibrée et facile à mettre en oeuvre. Les relations de saut généralisées provenant des sources sont détaillées et le système est testé dans un très large éventail de problèmes: les flux de tuyères, diphasique avec chimie et amortissement, les systèmes de relaxation ... Un chapitre est également consacré à des cas de résonance: en suivant Majda, nous étudions le problème stationnaire pour valider nos résultats. Ce travail conclut avec un écoulement diphasique en deux dimensions qui évolue de manière assez compliquée, les relations non conservatives étant traitées par un processus itératif convergent.
|
73 |
MODELISATION DE LA DISPERSION ET L'EVAPORATION DE SPRAYS DANS DES CHAMBRES DE COMBUSTION AERONAUTIQUESSierra Sanchez, Patricia 23 January 2012 (has links) (PDF)
De nos jours, la combustion représente encore un 90% de la production totale d'énergie au monde. La plupart des brûleurs de type industriel utilisent comme carburant des hydrocarbures en forme liquide. Cependant, un grand nombre d'études ont été dédiés aux flammes gazeuses et l'impact du spray liquide est encore loin d'être totalement compris. Le but de cette étude est l'amélioration de la modélisation des deux phénomènes principaux qui ont lieu entre l'atomisation du spray et la combustion, i.e. la dispersion des gouttes par la turbulence gazeuse et le procès d'évaporation dans le contexte de la Simulation Aux Grandes Echelles (SGE) des configurations complexes. Premièrement, l'approche Euler-Euler mésoscopique (Février et al. (2005)), basée sur une moyenne d'ensemble conditionnée et implémentée dans AVBP est améliorée. Le modèle de fermeture (Simonin et al. (2001); Kaufmann (2004)) pour les moments de deuxième ordre qui apparait dans les équations de transport résolues échoue quand appliqué à des configurations cisaillées (Riber (2007)). Plusieurs modèles proposés récemment par Masi (2010) et qui ont été valides a priori dans une configuration de nappe chargée de particules sont validés a posteriori dans la même configuration. Une analyse quantitative sur plusieurs cas avec diffèrent nombres de Stokes, nombres de Reynolds de la phase gazeuse et résolutions du maillage ont permit de retenir un modèle non-linéaire nommé 2EASM3, qui utilise le tenseur de déformations de la phase dispersée comme échelle de temps caractéristique. La deuxième partie a pour but l'amélioration du modèle d'évaporation implémenté dans AVBP. Ce modèle suppose une conduction infinie dans la phase liquide et symétrie sphérique dans la phase gazeuse ainsi que des lois simplifiées pour les propriétés thermodynamiques et de transport. Un nouveau modèle prenant en compte la dépendance de la viscosité du mélange gazeux avec la composition locale, et des nombres de Prandtl et Schmidt fixés par les valeurs à l'équilibre obtenus par moyen d'une simulation prenant en compte des lois complexes pour les propriétés thermodynamiques et de transport est proposé. Cette nouvelle méthode produit des résultats en bon accord avec les mesures expérimentales pour l'évaporation d'une goutte isolée en une atmosphère d'azote au calme sans pourtant augmenter le cout du calcul. Finalement, l'impact des nouveaux modèles est analysé dans une SGE de la configuration semi-industrielle MERCATO (García-Rosa (2008)). Bien que les données expérimentales ne soient pas suffisantes pour confirmer les résultats, les distributions de gouttes et de carburant gazeux sont significativement affectés par les modèles, ce qui pourrait avoir un impact direct sur le procès d'allumage.
|
74 |
LES combined with statistical models of spray formation closely to air-blast atomizer / Modélisation d'atomisation air-assistée au voisinage de l'injecteur : LES couplée avec les approches stochastiquesDeng, Tian 09 November 2011 (has links)
Cette thèse présente une extension de l'approche stochastique de l'atomisation primaire de type air assisté près d'un injecteur. Cette approche avait déjà été introduite dans les publications de Gorokhovski et al. Dans le cadre de la simulation des grandes échelles, la zone d'atomisation primaire est simulée comme un corps immergé avec une structure stochastique. Ce dernier est défini par la simulation stochastique de la position et de la courbure de l'interface entre le liquide et le gaz. La simulation de la position de l'interface est basée sur l'hypothèse de symétrie d'échelle pour la fragmentation. La normale extérieure à l'interface est modélisée en supposant une relaxation statistique vers l'isotropie. Les statistiques de la force du corps immergé servent de conditions aux limites pour le champ de vitesse issu de la LES ainsi que pour la production des gouttes de l'atomisation primaire. Celles-ci sont ensuite transportées par une approche lagrangienne. Les collisions entre les gouttes dans la zone d'atomisation primaire sont prises en compte par analogie avec l'approche standard de la théorie cinétique des gaz. Une fermeture est proposée pour la température statistique des gouttelettes. Cette approche est validée par des comparaisons avec les mesures expérimentales de la thèse de Hong. Les résultats numériques pour la vitesse et de la taille des gouttes dans le spray à différentes distances du centre du jet et de l'orifice de la buse sont relativement proches des résultats expérimentaux. Différentes conditions d'entrée pour la vitesse sont testées et comparées aux résultats expérimentaux. Par ailleurs, le rôle spécifique de la zone de recirculation devant le dard liquide est soulignée par le battement du dard liquide et la production de gouttelettes. / This thesis introduced an extension to stochastic approach for simulation of air-blast atomization closely to injector. This approach was previously proposed in publications of Gorokhovski with his PHD students. Our extension of this approach is as follows. In the framework of LES approach, the contribution of primary atomization zone is simulated as an immersed solid body with stochastic structure. The last one is defined by stochastic simulation of position-and-curvature of interface between the liquid and the gas. As it was done previously in this approach, the simulation of the interface position was based on statistical universalities of fragmentation under scaling symmetry. Additionally to this, we simulate the outwards normal to the interface, assuming its stochastic relaxation to isotropy along with propagation of spray in the down-stream direction. In this approach, the statistics of immersed body force plays role of boundary condition for LES velocity field, as well as for production of primary blobs, which are then tracked in the Lagrangian way. In this thesis, the inter-particle collisions in the primary atomisation zone are accounted also by analogy with standard kinetic approach for the ideal gas. The closure is proposed for the statistical temperature of droplets. The approach was assessed by comparison with measurements of Hong in his PHD. The results of computation showed that predicted statistics of the velocity and of the size in the spray at different distances from the center plane, at different distances from the nozzle orifice, at different inlet conditions (different gas velocity at constant gas-to-liquid momentum ratio, different gas-to-liquid momentum ratio) are relatively close to measurements. Besides, the specific role of recirculation zone in front of the liquid core was emphasized in the flapping of the liquid core and in the droplets production.
|
75 |
Experimental study of third (HFC) and fourth generation (HFO) refrigerants during flow boiling in singularities / Etude expérimentale de la troisième (HFC) et quatrième génération (HFO) frigorigènes lors de l'écoulement bouillant dans les singularitésPadilla Gomez, Miguel David 18 November 2011 (has links)
La réduction de charge de fluides frigorigènes dans les systèmes de production de froid est un enjeu important s'inscrivant dans les politiques environnementales sur la contribution de ces fluides à l'effet de serre. Une des voies menant à la réduction de charge est l'augmentation de la compacité des échangeurs de chaleur, conduisant inévitablement à la conception d'évaporateurs plus complexes. Néanmoins, si de nombreuses études ont été publiées sur l'analyse thermo-hydraulique de l'ébullition convective de frigorigènes dans les tubes horizontaux, très peu ont été menées sur des écoulements diphasiques dans des géométries non conventionnelles. Cette thèse vise à étudier les caractéristiques des écoulements diphasiques avec changement de phase des frigorigènes de troisième (HFC) et de quatrième génération tels que le HFO-1234yf dans des singularités (coudes de retour, contractions). Ainsi, un banc d'essais expérimental a été spécifiquement conçu et construit pour mener des expériences sur l'ebullition convective. Ce banc d'essais a permis de visualiser les régimes d'écoulement et de mesurer les chutes de pression dans les singularités. Des visualisations de régimes d'écoulement ont été réalisées pour les frigorigènes HFO-1234yf, R-134a et R-410A dans des tubes droits et dans des singularités. Les résultats ont été confrontés avec diverses méthodes de prédiction de la litérature. Les perturbations sur l'écoulement causées par ces singularités et leurs effets sur le comportement hydrodynamique des fluides frigorigènes ont également été etudiés. Une serie d'éxperience a été menée pour déterminer le longueur de perturbation en amont et en aval des singularités. Finalement, les valeurs expérimentales des chutes de pression ont été confrontées à différentes corrélations de la littérature. Des nouvelles méthodes de prédiction de la chute de pression pour des écoulements diphasiques dans des singularités sont présentées. / The refrigerant charge reduction in HVAC\R systems is an important issue because it falls within environmental policies regarding refrigerants contributions to the greenhouse effect. A way to move toward charge reduction is to increase the compactness of heat exchangers, which means more complex designs of the evaporators. Nevertheless, while a large amount of studies have has been published on the thermal and hydraulic analysis of flow boiling of refrigerants in horizontal tubes, very little attention has been given to flow boiling in geometries different from straight tubes. This PhD thesis aims at studying the flow boiling characteristics of third generation (HFC) and fourth generation fluids such as HFO-1234yf in geometries which modify the fluid dynamics and two-phase flow with respect to horizontal straight tubes. To achieve this goal, an experimental test facility was specifically designed and built to conduct refrigerant evaporation experiments. This test facility allowed to perform flow regimes visualizations and pressure drop measurements in singularities (such as sudden contractions and return bends). First, two-phase flow regimes visualizations have been carried out using HFO-1234yf, R-134a and R-410A either in straight tubes or in singularities. A qualitative analysis of the flow behavior and also several comparisons to flow pattern prediction methods from the literature were conducted. The second objective of this work was to characterize the flow disturbances caused by singularities such as sudden contractions and return bends, and to study their effects on the hydrodynamic performance (e.g. pressure drop) of refrigerants.Lastly, pressure drop databases for third and fourth generation refrigerants are presented.
|
76 |
Contribution à l'étude des propriétés thermiques et hydrodynamiques d'un écoulement d'hélium normal (5HeI) diphasique en circulation naturelle pour le refroidissement des aimants supraconducteurs / Contribution to the study of thermal and hydrodynamical properties of HeI two phase natural circulation flow for cooling superconducting magnetsBenkheira, Lahcène 29 June 2007 (has links)
La méthode de refroidissement basée sur le principe thermosiphon présente un grand intérêt en raison de sa simplicité, de sa nature passive et de son coût faible. Elle est adoptée pour le refroidissement à 4,5 K de l’aimant supraconducteur du détecteur de particules CMS auprès du LHC en construction au CERN à Genève. Le travail présenté dans cette thèse étudie expérimentalement les propriétés thermiques et hydrodynamiques d’un écoulement d’He I diphasique en circulation naturelle. Le dispositif expérimental utilisé consiste en une boucle thermosiphon monobranche composée principalement d’un séparateur de phases, d’un tube descendant et d’une section d’essai. Les expériences ont été réalisées en faisant varier plusieurs paramètres tels que le diamètre des sections d’essai (10 mm ou 14 mm) et le flux de chaleur allant jusqu’à l’apparition de la crise d’ébullition. Ces expériences ont permis de déterminer les lois d’évolution des différentes grandeurs caractérisant l’écoulement (le débit massique de circulation, le débit massique vapeur, le titre massique, le coefficient de friction et le coefficient d’échange thermique) en fonction de la densité du flux de chaleur appliquée. Au regard des résultats obtenus, nous discutons la validité des différents modèles classiques existants dans la littérature. Nous montrons que le modèle homogène est le modèle le mieux adapté pour prédire les propriétés hydrodynamiques de ce type d’écoulement dans la gamme de titre massique 0?x?30%. De plus, nous proposons deux modèles pour la prédiction du coefficient de transfert de chaleur diphasique et la densité de flux de chaleur critique. Le premier considère que les effets de la convection forcée et de l’ébullition nucléée agissent simultanément et contribuent au transfert de chaleur. Le deuxième corrèle la densité de flux de chaleur critique mesurée en fonction du rapport altitude sur diamètre / The method of cooling based on the thermosiphon principle is of great interest because of its simplicity, its passivity and its low cost. It is adopted to cool down to 4,5 K the superconducting magnet of the CMS particles detector of the Large Hadron Collider (LHC) experiment under construction at CERN, Geneva. This work studies heat and mass transfer characteristics of two phase He I in a natural circulation loop. The experimental set-up consists of a thermosiphon single branch loop mainly composed of a phase separator, a downward tube, and a test section. The experiments were conducted with varying several parameters such as the diameter of the test section (10 mm or 14 mm) and the applied heat flux up to the appearance of the boiling crisis. These experiments have permitted to determine the laws of evolution of the various parameters characterizing the flow (circulation mass flow rate, vapour mass flow rate, vapour quality, friction coefficient, two phase heat transfer coefficient and the critical heat flux) as a function of the applied heat flux. On the base of the obtained results, we discuss the validity of the various existing models in the literature. We show that the homogeneous model is the best model to predict the hydrodynamical properties of this type of flow in the vapour quality range 0?x?30%. Moreover, we propose two models for the prediction of the two phase heat transfer coefficient and the density of the critical heat flux. The first one considers that the effects of the forced convection and nucleate boiling act simultaneously and contribute to heat transfer. The second one correlates the measured critical heat flux density with the ratio altitude to diameter
|
77 |
Modélisation numérique d’écoulement diphasique compressible et transport réactif en milieux poreux - Applications à l'étude de stockage de CO2 et de réservoir de gaz naturel. / Numerical simulation of compressible two-phase flow and reactive transport in porous media - Applications to the study of CO2 storage and natural gas reservoir.Sin, Irina 08 December 2015 (has links)
Les activités humaines dans la subsurface se développent rapidement (stockage de déchets,nouvelles techniques minières, stockage à haute fréquence de l’énergie), alors que dans le même temps les attentes du public et des autorités s’intensifient. L’évaluation de chaque étape de ces opérations souterraines repose sur des études détaillées de la sûreté et des impacts environnementaux.Elles reposent sur des simulateurs élaborés et sur de la modélisation multiphysique. Avec leur approche orientée processus, les simulations en transport réactifs proposent une méthode efficace pour comprendre et prévoir le comportement de ces systèmes complexes, à différentes échelles de temps et d’espace.Le but de ce travail est d’intégrer la résolution de l’écoulement diphasique compressible dans le cadre de codes de transport réactifs à l’aide d’une méthode de séparation d’opérateurs. Un module multiphasique a été créé dans le code de transport réactif HYTEC. Une nouvelle approche a ensuite été développée pour coupler écoulement multicomposant multiphasique compressible, description de propriétés thermo-dynamiques complexes pour les fluides, avec des codes de transport réactif. La méthode a été intégrée dans HYTEC. Des cas de validation sont proposés, puis des exemples d’application pour la simulation du stockage souterrain de CO2 et des impuretés associées. / Human activity in the subsurface has rapidly been expanding and diversifying (waste disposal, new mining technologies, high-frequency storage of energy), while the public and regulatory expectations keep growing. The assessment of each step of underground operations requires careful safety and environmental impact evaluations. They rely on elaborate simulators and multiphysics modeling. With its process-based approach, reactive transport simulation provides an effective way to understand and predict the behavior of such complex systems at different time and spatial scale.This work aims at incorporating a compressible multiphase flow into conventional reactive transport framework by an operator splitting approach. A multiphase flow module is developed in the HYTEC reactive transport software. A new approach is then developed to fully couple multiphase multicomponent compressible flow, the complex thermodynamic description of the fluid properties, with existing reactive transport codes. The method is implemented in HYTEC. Some validation is provided, before application to the simulation of underground storage of CO2 and associated impurities.
|
78 |
Interaction between Flow Induced Pulsations and Multiphase Flows in Gas Liquid Systems / Oscillations auto-induites en systèmes de conduites en présence d’écoulement diphasique air/eauSanna, Francesco 15 April 2016 (has links)
Le couplage entre un écoulement instable et des résonances acoustiques dans des systèmes de conduites peut conduire à des phénomènes d’oscillations auto-induites. Ce type de phénomènes trouve principalement place dans des conduites latérales fermées, par exemple dans des systèmes de transport ou de compression de gaz. L’objectif de ce travail est d’étudier les oscillations auto-induites dans le cas où le fluide transporté ne se limite pas à un gaz, mais est un mélange de gaz et de liquide. Les pulsations sont mesurées dans des conduites latérales fermées, pour deux types de configurations (en tandem et en croix), avec écoulement d’un mélange variable d’air et d’eau. La position de l’injection d’eau est variable afin d’obtenir plusieurs régimes d’écoulement diphasique. Les résultats indiquent que la présence d’eau a un effet important sur les niveaux de pulsations dans les conduites. Cet effet a pu être attribué à deux mécanismes dus à la présence d’eau : les instabilités de couches de mélange sont modifiées et l’amortissement des ondes acoustiques est amplifié.Le deuxième mécanisme a été quantifié à l’aide de mesures sur un montage expérimental dédié conçu pour avoir un écoulement stratifié. On a observé que, dans tous les cas, la présence d’eau augmente l’amortissement. Cette augmentation a pu être attribuée à la réduction de la section effective de la conduite (due au remplissage partiel par l’eau) et à l’augmentation de la friction turbulente à l’interface entre les phases liquide et gazeuse. / Coupling between flow instabilities and acoustic resonances in ducts with closed side branches leads to Flow Induced Pulsations (FIPs). This is a typical phenomenon in engineering applications (gas transport systems, compressor installations, and chemical plants). The objective of this work is to extend the knowledge about FIPs when the transported medium is not uniquely gas but a combination of gas and (a small quantity of) liquid. For two configurations of double side branches (in tandem and in quasi-cross), the amplitude of pressure pulsations in the side branches was measured for different liquid injection rates. This was repeated with the liquid injection point located at different places to allow different flow regimes at the pipe connections. The results show a strong effect of the water content on the pulsations. On basis of these results and additional measurements, the following hypotheses for the effect of liquid were made: (1) interaction of the liquid with the flow instability and (2) increase of the acoustical damping in the ducts in presence of liquid.The effect of liquid on damping was measured with a dedicated test setup designed to have a stratified flow. It was found that the liquid always increases the acoustical damping, mainly due to the reduction of the effective cross section by the liquid, and because of the increased turbulent friction at the interface between gas and liquid.
|
79 |
Mise à l’échelle d’un écoulement diphasique avec gravité dans un milieu géologique hétérogène : application au cas de la séquestration du CO₂ / Upscaling of a two-phase flow model including gravity effect in geological heterogeneous media : application to CO₂ sequestrationNgo, Tri Dat 26 January 2016 (has links)
Ce travail de thèse porte sur la modélisation mathématique et la simulation numérique de la migration par gravité et capillarité du CO₂ supercritique injecté dans un site de séquestration géologique hétérogène. Les simulations sont réalisées à l'aide du code DuMux. Particulièrement, on s'intéresse à la mise à l'échelle, de l'échelle de la cellule à l'échelle du réservoir, d'un modèle d'écoulement diphasique CO₂ -saumure, au sein d'un milieu stratifié périodique constitué d'un réseau de barrières peu perméables horizontales, continues ou discontinues. La mise à l'échelle est effectuée par la méthode asymptotique à double échelle. Dans un premier temps, on considère le cas d'une colonne verticale parfaitement stratifiée. Un modèle homogénéisé est développé puis validé par simulation numérique pour différentes valeurs du nombre capillaire et du flux incident de CO₂ . La méthode d'homogénéisation est appliquée au cas d'un écoulement dans un milieu bidimensionnel constitué de strates discontinues. Par l'effet de gravité, le CO₂ s'accumule sous les strates peu perméables, ce qui conduit à un problème mathématique local non standard. Cette stratification est modélisée à l'aide de l'approche des courants de gravité. L'approche est étendue au cas des strates semi-perméables et en prenant en compte la capillarité. Le modèle mis à l'échelle est comparé à des simulations numériques effectuées pour différents types de strates, avec ou sans pression capillaire, et sa limite de validité est discutée pour chacun de ces cas. La dernière partie de la thèse est dédiée à l'étude des performances du code DuMux pour simuler par calcul parallèle l'injection et la migration de CO₂ dans des milieux hétérogènes tridimensionnels (milieu périodique stratifié, milieu fluviatile et milieu réservoir SPE10). / This work deals with the mathematical modeling and the numerical simulation of the migration under gravity and capillarity effects of the supercritical CO₂ injected into a geological heterogeneous sequestration site. The simulations are performed with the code DuMux. Particularly, we consider the upscaling, from the cell scale to the reservoir scale, of a two-phase (CO₂ -brine) flow model within a periodic stratified medium made up of horizontal low permeability barriers, continuous or discontinuous. The upscaling is done by the two-scale asymptotic method. First, we consider perfectly layered media. An homogenized model is developed and validated by numerical simulation for different values of capillary number and the incident flux of CO₂ . The homogenization method is then applied to the case of a two-dimensional medium made up of discontinuous layers. Due to the gravity effect, the CO₂ accumulates under the low permeability layers, which leads to a non-standard local mathematical problem. This stratification is modeled using the gravity current approach. This approach is then extended to the case of semi-permeable stratas taking into account the capillarity. The upscaled model is compared with numerical simulations for different types of layers, with or without capillary pressure, and its limit of validity is discussed in each of these cases. The final part of this thesis is devoted to the study of the parallel computing performances of the code DuMux to simulate the injection and migration of CO₂ in three-dimensional heterogeneous media (layered periodic media, fluvial media and reservoir model SPE 10).
|
80 |
Magnetic flux distorsion in two-phase liquid metal flow / Étude d’une méthode de détection de gaz dans du sodium liquide par méthode électromagnétique pour les Réacteurs à Neutrons Rapides au Sodium (RNR Na)Kumar, Mithlesh 23 March 2016 (has links)
Cette thèse se situe dans le cadre du projet ASTRID du CEA. La surveillance de la présence de gaz dans un RNR est importante pour son fonctionnement en toute sécurité. Cette thèse porte sur la détection et la caractérisation de présence de gaz par un débitmètre à distorsion de flux (DDF). Du point de vue technologique, l’objectif est d’évaluer la faisabilité de l’utilisation d’un DDF, et s’il est possible de mesurer le débit et le taux de vide simultanément avec un DDF. Du point de vue de la physique, le DDF mesure la perturbation du flux magnétique due aux courants de Foucault induits par l’hydrodynamique d’un écoulement de liquide conducteur. En présence de vide dans le conducteur, une nouvelle source de perturbation apparaît du fait de la modification, par le taux de vide, de la distribution du champ magnétique. En effet, la présence de vide agit sur les distributions locales de courants électriques dues au couplage des phénomènes d’induction et de forces de Lorentz. Notre objectif est de comprendre la nature de ces couplages et de proposer une méthode qui permette de caractériser la présence de vide à partir du signal de sortie du DDF, de mesurer le taux de vide et d’étudier la sensibilité de cette mesure aux variations des paramètres de l’écoulement et du champ électromagnétique (vitesse et pulsation, en particulier). Dans ce travail, des expériences spécifiques ont été développées pour étudier les effets de la vitesse, du taux de vide et de la pulsation du flux magnétique sur la réponse d’un DDF. Ces expériences modélisent un écoulement diphasique idéal (écoulement piston) consistant en une barre mobile d’aluminium contenant des trous et des cannelures pour simuler le taux de vide. Dans ces expériences la vitesse, le taux de vide et la pulsation sont parfaitement contrôlé et leur domaine de variations sont les suivants : 0<U<1 m/s pour la vitesse, 0< <6.9% pour le taux de vide et 1500<<12000 rad/s pour la pulsation (dans cette gamme de pulsations, on notera que la profondeur de pénétration du champ électromagnétique est de l’ordre de, mais plus petit que, le rayon du barreau d’aluminium). Une méthode d’ellipse-fitting appliqué au signal de sortie du DDF a été utilisée pour caractériser l’effet du taux de vide. Les résultats montrent que le DDF est sensible pour des fractions volumique de vide entre 0.3 % and 6.9% . En outre, la réponse aux variations de taux de vide est insensible à la vitesse du barreau. Cette technique est peu adaptée aux mesures des faibles taux de vide (< 0.1%). Une deuxième approche a été développée, sur la base de l’analyse du signal démodulé du DDF. Cette analyse s’appuie sur un modèle théorique du flux magnétique obtenu par un développement au 1er ordre par rapport à U et a. Ce modèle permet d’interpréter les résultats expérimentaux en termes de contributions de U et a au flux magnétique. Malgré le couplage fort entre l’induction Faraday et les forces de Lorentz, les résultats montrent que les contributions de U et alpha peuvent être séparées correctement pour des petites valeurs du nombre de Reynolds magnétique (Rem < 0.12), et de faible taux de vide ( < 6.9%). Un résultat très important est que la contribution de a sur le flux magnétique est insensible aux variations de vitesse dans cette gamme de Rem . De plus, l’effet de la géométrie du vide a été étudié pour deux configurations : cannelures et trous. Il a été observé que la géométrie du vide n’a pas d’effet sur les variations du flux magnétique avec a. Cette seconde approche est plus sensible aux variations de fraction volumique du vide que la méthode de l’ellipse-fitting. Enfin, des expériences préliminaires avec un liquide métallique (galinstan) contenant des perles de verre, ont été faites et ont montré une bonne sensibilité du signal du DDF avec la vitesse et le taux de vide. En conclusion, ce travail a montré qu’il est possible de mesurer simultanément un taux de vide et un débit pour la gamme de variations des paramètres étudiés. / A Generation IV Sodium cooled Fast Reactor (SFR) is being researched and developed at CEA, Cadarache France under the project named ASTRID. Monitoring gas presence in SFR is important with respect to its safe operation. In accordance with the principles of diversity, techniques based on different measurement principles have been proposed. This thesis concerns the detection and characterization of void using magnetic flux perturbation principle. An Eddy Current Flow Meter (ECFM) device is used for this purpose. From the technological point of view, the objective is to evaluate the feasibility of ECFM as a flow and/or void monitoring/characterizing device; and to determine which parameters are of interest and what are the precision of these measurements; and whether it is possible to measure the flow rate and void fraction simultaneously with the same ECFM device. From the physics point of view, the ECFM system involves the magnetic flux perturbation due to voids in the presence of Faraday induction and Lorentz force effects. Therefore ECFM integrated signal contains informations about the void, Faraday induction and Lorentz force effects based perturbation in magnetic flux and their couplings. Our objective is to understand the nature and extent of these couplings. Specific experiments have been developed to study the effects of flow velocity, void fraction and magnetic flux pulsations on the response of an ECFM. It consists in modeling the two-phase flow by a moving aluminium rod (plug flow) with holes and grooves to simulate voids. Flow velocity range of variation was 0<U<1 m/s, void fraction 0<<6.9% and pulsation 1500<<12000 rad/s (for this range of pulsations the electromagnetic skin depth is of order but smaller than the aluminium rod radius). An ellipse fitting method was proposed to analyze the output signal of the secondary coils. The results show that ECFM is sensitive to void fractions between 0.3 % and 6.9%. Furthermore, the response to void fraction is insensitive to the mean velocity of the twophase medium. A second approach based on demodulation analyses of the secondary coils output signal has been developed. We have proposed a theoretical model based on a first order expansion of magnetic flux in U and a. With this model it was possible to interpret the experimental results in terms of contributions of U and a. Despite the strong coupling between Faraday induction and Lorentz force effects, the results show that the contributions of U and a. can be well separated at low magnetic Reynolds number (Rem < 0.12) and low a values ( < 6.9%). A very important result is that the contribution of _ on magnetic flux is insensitive to variations of velocity in this range of Rem. Moreover, different geometries of void have been studied: grooves and holes. It was observed that the geometry of void do not change the variation of magnetic flux with a. This second approach revealed to be more sensitive to void fraction variations than ellipse fitting method. Finally, preliminary experiments with liquid metal galinstan with glass beads were done, which showed sensitivity of ECFM signal with velocity and void. In conclusion, this work has shown that ECFM can measure simultaneously void and velocity in the range of parameters studied, in particular 0.06%< < 6.9%.
|
Page generated in 0.0896 seconds