Le développement et la certification de systèmes propulsifs aéronautiques nécessitent une phase d'essais sur banc moteur. Ces essais permettent entre autres de caractériser les limites d'allumage des foyers de combustion de turbomachines, mais sont extrêmement coûteux et générateurs de délais pour l'industriel. Afin de limiter leur recours, il est nécessaire de développer des méthodes permettant de prévoir de la façon la plus fiable possible les performances d'allumage d'un prototype de chambre de combustion, et ce aussi bien pour des conditions de décollage au niveau de la mer que pour le cas critique de la haute altitude.L'objectif de cette thèse est de contribuer au développement et à la validation d’outils numériques pour la prévision de l'allumage des foyers de combustion à partir de données expérimentales obtenues sur le banc MERCATO de l’ONERA. Ces travaux ont été conduits en étroite coopération avec TURBOMECA. Des améliorations ont été apportées à un modèle permettant de simuler l'allumage d'un brouillard de carburant suite à un dépôt d'énergie par bougie à arc électrique. Ce modèle a été couplé à un code de calcul multiphysique (code CEDRE) afin d'établir des cartographies d'allumage à partir d'un champ aérodiphasique moyenné, mais également de simuler la phase de propagation de la flamme à l’ensemble du foyer par approche RANS pseudo-stationnaire. Afin de valider ces deux approches, une caractérisation expérimentale d'une maquette de chambre mono-secteur a été réalisée sur le banc d'essai MERCATO. Ces essais ont permis de constituer une banque de données très fournie pour différents cas tests. La comparaison de la cartographie d’allumage expérimentale à celle fournie par la simulation numérique donne des résultats très satisfaisants et encourageants en vue d’une application à un foyer réel de turbomachine. En parallèle, des résultats très prometteurs ont été obtenus sur une nouvelle maquette de chambre de combustion trisecteur,dérivée d'un foyer industriel, permettant de comparer ses limites d’allumage à celles obtenues sur la chambre réelle. Cette configuration fera l’objet par la suite d’une étude détaillée de la propagation de la flamme intersecteurs, impossible à réaliser sur la géométrie d'un foyer annulaire complet. / Design of aircraft engines requires tests on engine benches. These tests allow characterizing combustor ignition limits, but are extremely expensive and time consuming. In order to limit their number, it is necessary to develop alternative methods enabling to predict the ignition performances of a combustor prototype, for both ground conditions and high altitude conditions, the latter being particularly critical.The purpose of this thesis is to contribute to the development and validation of numerical tools enabling to predict ignition performances of combustor. Validation will be possible using an experimental data base obtained on the MERCATO test bench, at ONERA. This work results from a close cooperation with TURBOMECA. A numerical model, previously developed to predict the ignition of fuel spray following a spark discharge, has been improved. This model has been combined to a multiphysics CFD code (CEDRE) in orderto build ignition map from a mean two-phase flow field, and also to compute the flame propagation stage using RANS approach. In order to validate both methods, experimental characterization of a one-sectorcombustor has been performed on the MERCATO test rig. A complete data base for validation of CFD code isnow available. Comparison of experimental and numerical ignition mapping showed good agreements for the different tests cases, and seems encouraging for an application on an industrial combustor. In the same time,promising results have been obtained with a new experimental three-sector combustor. Its geometry is basedon an industrial combustor, allowing a comparison of ignition performances between the simplified and the industrial combustors

Cette thèse est motivée par le besoin de compléter les connaissances actuelles des phénomènes ayant lieu lors d’un renvoi d’eau dans un lit de débris radioactifs, opération appelée « renoyage » et qui intervient dans une séquence d’accident grave où un cœur de réacteur nucléaire est dégradé suite à une perte prolongée de refroidissement primaire. Notre étude, de nature expérimentale, vise à mieux caractériser la crise d’ébullition en convection forcée, dans un milieu poreux chauffant. Le cœur du dispositif expérimental est un milieu poreux modèle quasibidimensionnel, composé de 276 cylindres disposés entre deux plaques de céramique distantes de 3 mm, dont l’une, transparente, permet de visualiser les écoulements. Les cylindres, de 2 mm de diamètre, sont des sondes thermo-résistives qui ont une double fonction : elles sont utilisées comme éléments chauffants et comme capteurs de température. Une boucle fluide permet de contrôler le débit d’injection de liquide dans la section test, la température d’injection ainsi que la pression. La section test est placée verticalement, le liquide est injecté par le bas à une température proche de la saturation. Dans une première série d’expériences, la puissance thermique dissipée globalement par un ensemble de cylindres chauffants est augmentée de façon progressive jusqu’à atteindre l’assèchement d’une zone du milieu poreux. Les résultats montrent deux types de phénoménologies dans le déclenchement de la crise d’ébullition. Pour des débits d’injection faibles (densités de flux massique de l’ordre de 4 kg.m^-2.s^-1 maximum), l’atteinte de la puissance d’assèchement se traduit par un lent recul du front diphasique jusqu’à sa stabilisation en haut de la zone chauffée ; en aval de la zone chauffée, l’écoulement est monophasique vapeur. Pour des débits d’injection plus élevés, la crise d’ébullition apparaît autour d’un des éléments chauffants, conduisant à une ébullition en film localisée, tandis qu’un écoulement diphasique liquide-vapeur continue de parcourir l’aval de la section test. Les visualisations de ces expériences permettent d’identifier qualitativement la structure des écoulements. D’autres expériences consistent à mesurer le flux critique local autour d’un cylindre choisi, pour différentes configurations d’écoulements. Le débit d’injection est fixé. Une puissance de chauffe est imposée à une ligne horizontale de cylindres en amont du cylindre choisi. Les résultats montrent que le flux critique sur ce cylindre diminue en fonction de la puissance délivrée à la ligne chauffée. La distance du cylindre étudié à la ligne chauffée semble avoir peu d’influence sur le flux critique. Des visualisations expérimentales sont utilisées pour caractériser l’écoulement diphasique en aval de la ligne chauffée, dans le but de mettre en relation le flux critique local avec des paramètres hydrodynamiques (saturations, vitesses des phases). Les images obtenues sont difficiles à exploiter. Afin de calibrer les paramètres des algorithmes de traitement d’images, nous avons reproduit une cellule d’essai de géométrie identique à l’originale, mais où l’on injecte du gaz par une ligne de cylindres en amont de la section test dans une configuration d’écoulement diphasique isotherme. Dans ce dispositif, le débit d’injection de gaz est contrôlé et mesuré. Les visualisations obtenues servent alors de références auxquelles sont comparées les visualisations d’ébullition convective. / This work is motivated by the need to better understand the phenomena occurring while some water is injected into a heated porous debris bed. This reflooding operation is a part of the planned mitigation procedure used during a Loss Of Coolant Accident (LOCA) that may occur into a nuclear power plant and results into a severe core damage. Our experimental study aims to characterize the boiling crisis that can happen in a boiling flow taking place within a heatgenerating model porous medium. The test section is a two-dimensional model porous medium, composed of an array of 276 cylinders placed between two ceramic plates spaced from one another by 3 mm, one of which is transparent and allows visualizations of the flow. The 2 mm diameter cylinders are Pt100 resistance temperature detectors that perform a dual function: they act as heating elements (heated by Joule effect) and are also used as temperature probes. A fluid loop allows controlling the liquid injection flow rate, its inlet temperature as well as its pressure. The test section is held vertically, the liquid injected from bottom at a temperature close to the saturation temperature. In a first series of experiments, the thermal power applied to a bundle of heating cylinders is progressively increased until a dry zone is detected in the porous medium. Two kinds of phenomenology are observed during these “dryout experiments”. First, at low liquid injection rate (4 kg.m^-2.s^-1 maximum mass flux), reaching the dryout power results into a liquid front receding down to the upper limit of the heated zone, while downstream the heated zone, the porous medium is vapour-saturated. Second, at higher flow rate, the boiling crisis happens at the surface of a single heating element, resulting in a local film boiling, whereas a two-phase flow still go through the whole test section. High-speed visualizations allow characterizing the flow regimes. Other experiments focus on determining the local critical heat flux on a given cylinder, for different upstream flow configurations. The inlet liquid flow rate is fixed. A thermal power is uniformly applied to a line of heating cylinders, upstream the cylinder under study. Results show that the local critical heat flux decreases as the power applied to the heated line increases. The distance from the cylinder under study to the heated line seems not to have a significant effect on the critical heat flux. Visualizations are used to characterize the two-phase flow upstream the heated line, aiming at expressing the critical heat flux as a function of the hydrodynamic parameters (saturations, phase velocities). The image analysis is particularly challenging. In order to calibrate the image processing parameters, we use a second model porous medium with the same geometry as the heat generating one, but where an isothermal two-phase flow is obtained by injecting gas into the liquid flow rather than generated by boiling. The gas injection flow rate is controlled and measured. Isothermal two-phase flow visualizations provide a reference case and are compared to flow boiling visualizations.

Le nouveau moteur cryogénique de l'étage supérieur du nouveau lanceur Ariane présente la particularité d'être plusieurs fois réallumable, une fois la mise en orbite du lanceur. Le réallumage d'un moteur est particulièrement difficile durant les conditions de vol spatial. Ce moteur est composé d'un dôme LOX alimenté en oxygène liquide (LOX) qui est approvisionné par une vanne à boisseau positionnée en entrée d'une canne d'alimentation. Le mélange liquide / gaz formé dans le dôme LOX est injecté dans la chambre de combustion à travers des injecteurs reliant le dôme à la chambre. En conséquence, la distribution de l'écoulement diphasique en sortie des injecteurs revêt une importance particulière en terme d'allumage, de l'ouverture à la fermeture de cette vanne. La prise en compte de ces conditions de vol est primordiale pour qualifier le moteur. Cependant ces conditions ne peuvent pas être reproduites de façon représentative au cours d’essais au sol. Dans le cadre de ces études, un programme de recherche a été mis en place par le CNES (Centre National d'Études Spatiales) et SAFRAN Snecma pour étudier le remplissage du dôme LOX, via des études expérimentales et numériques. L'objectif est de connaître les conditions aux limites en sortie des injecteurs qui sont déterminantes pour appréhender la phase d'allumage dans la chambre de combustion. Des expériences ont été menées au LEGI (Laboratoire des Écoulements Géophysiques et Industriels) avec des fluides de substitution (de l'eau et de l'air), sans transfert de masse et de chaleur, sur la maquette du dôme d’alimentation d'un moteur de fusée. Les travaux présentés, menés à l'IMFT (Institut de Mécanique des Fluides de Toulouse), tentent de reproduire les expériences réalisées à l'aide de simulations numériques 3D incompressibles diphasiques. La géométrie du domaine de calcul est représentative de la maquette expérimentale, qui est composée d'une canne d'alimentation, d'un dôme, d'un allumeur et d'un grand nombre d'injecteurs. Le but de cette étude est de démontrer la faisabilité d'un calcul 3D instationnaire diphasique du remplissage du dôme oxygène du moteur-fusée avec le code industriel NEPTUNE_CFD, en prenant en compte la géométrie réelle et les phénomènes physiques prépondérants. La comparaison des prédictions avec les résultats expérimentaux est réalisée afin d'évaluer la capacité du code à prédire l'écoulement à l'aide des modèles de fermeture disponibles. Enfin, plusieurs études de sensibilité sur les modèles de fermeture sont menées pour estimer leur influence sur les résultats des simulations. Un travail important a été effectué pour imposer les mêmes conditions d'entrée que dans les expériences. Des études ont également été conduites sur un injecteur isolé. / New generation cryogenic upper-stage rocket engines are planned to be restartable during the orbit mission. The re-ignition of the engine is particularly difficult in space flight conditions. The engine contains a LOX dome fed with liquid oxygen (LOX) supplied by a bushel valve through a pipe. The gas / liquid mixture forming in the dome is injected into the combustion chamber through a number of injectors. Therefore the two-phase flow distribution at injectors outlet carries a real importance in terms of the ignition from the opening to the closing phases of the main valve. These flight conditions are of paramount importance, however, they are truly difficult to reproduce by experimental ground tests. In the framework of these studies, a research program set up by CNES (the French Space Agency) and SAFRAN Snecma, tries to study the LOX dome filling, through experiments and numerical studies. The aim is to identify the phenomena at sake to know the limit conditions at injectors, which will determine the ignition stage in the combustion chamber. Experiments are carried out at LEGI (Geophysical and Industrial Flows Laboratory in Grenoble) with substitution fluids (air and water), without heat and mass transfer on a rocket engine mockup. The work presented here, conducted at IMFT (Fluid Mechanics Institute in Toulouse), intends to reproduce the experimental results using incompressible two-phase flow simulations. The geometry used is representative of the experimental mockup composed of a feeding pipe, a dome, an igniter pipe and injectors. The aim of this study is to demonstrate the feasibility of a 3D unsteady two-phase flow calculation with the industrial code NEPTUNE_CFD, to simulate the LOX dome filling of the rocket engine, by taking into account the real geometry and the preponderant physical phenomena. The comparison of the predictions with the experimental results is carried out in order to estimate the code capability to predict the flow behavior, according to available closure laws. Finally, several sensitivity studies on the closure laws have been conducted to assess their influence on the numerical results. An important work has been carried out to obtain the proper inlet conditions to be imposed in the code in coherence with the experiments. Studies have equally been conducted on an isolated injector.

Dans l’industrie chimique, la majorité des réacteurs est équipée d’un dispositif « évent de sécurité » permettant d’éviter leur éclatement en cas de surpression accidentelle conséquente à un emballement. La nature physique du rejet (gazeux ou diphasique gaz-liquide) influe fortement sur la taille requise de l’évent. L’objectif de notre étude est de connaitre mieux le comportement de l’écoulement (désengagement) en cas d’un emballement de réaction et de prédire la nature du rejet (mono ou diphasique) notamment à l’échelle industrielle.Une étude expérimentale d’une réaction d'estérification fut réalisée dans un réacteur en verre de 0,5 l afin de visualiser la nature de l’écoulement et de mesurer la fraction de vide moyenne (ᾱ) lors de l’emballement de réaction. Une étude paramétrique fut effectuée de façon à identifier, pour le régime d'écoulement et la transition, les paramètres clés : la viscosité et l’agitation. Un calorimètre pseudo-adiabatique (VSP2) fut utilisé pour étudier la thermodynamique de l’emballement d’estérification. Les bilans massique et énergétique ont été utilisés pour calculer la vitesse superficielle de la vapeur (jg,max). En outre, le jg, max fut calculé aussi à partir de différentes corrélations, issues d’études en colonne à bulles. Une comparaison entre les résultats obtenus a montré l’inadéquation de ces corrélations pour un système réactif. Une carte de régime pour un système réactif a été construite pour la première fois (jg,max versus ᾱ, combinée à l'observation du régime d'écoulement dans le réacteur), avec une méthodologie pour la prédiction de la nature du rejet en cas d’emballement. / In chemical industry, most reactors are equipped with an emergency relief vent to prevent bursting in case of accidental overpressure due to a runaway reaction scenario. The physical nature of the vent release (gas phase or gas-liquid) strongly influences the necessary vent size. The objective is to enhance the knowledge on the flow behavior (disengagement) during a runaway reaction and to be able to predict the nature of the vent flow (1 or 2-phase) namely at industrial scale.Experiments of esterification were done in a 0.5 l glass reactor in order to visualize the flow pattern (hydrodynamic) and to measure the average void fraction (ᾱ) during the runaway. A parametric study was carried out to identify the key parameters onto the flow pattern and transition, they were: the viscosity and the stirring. A pseudo-adiabatic calorimeter (VSP2) was used to study the runaway P and T profile and the kinetic of the esterification. The thermodynamic data and the mass & heat balances were used to calculate the superficial vapor velocity (jg,max). Furthermore, this jg,max was also calculated from different correlations, obtained from bubble column research. A comparison showed the incorrectness of using these correlations for a reactive system. A flow pattern map (jg,max versus ᾱ combined with the observation of the flow pattern,) for a reactive system was built up for the first time together with a methodology to predict the vent flow nature.

[pt] Diversos processos práticos apresentam uma seqüência de diferentes configurações de escoamentos bifásicos, com diferentes padrões de escoamento. Dependendo das razões entre as velocidades superficiais do gás e do líquido, diferentes padrões de escoamento bifásico podem ser encontrados. O padrão de escoamento slug pode ser formado a partir do padrão estratificado, devido ao crescimento das instabilidades hidrodinâmicas ou devido às ondulações no terreno onde se instala a tubulação. O presente trabalho consiste na análise do escoamento bifásico no padrão slug, ao longo de tubulações horizontais, mediante a aplicação do modelo de dois fluidos em sua forma transiente e unidimensional. Através de uma análise de estabilidade de Kelvin-Helmholtz para escoamento estratificado, estima-se a possível faixa de operação para a obtenção de um modelo matemático bem posto. Diferentes técnicas de discretização foram implementadas e testadas. Comparações com dados disponíveis na literatura foram realizadas. / [en] Several practical processes show a sequence of different configurations of two-phase flows, with different flows pattern. Depending upon the relation between the superficial velocities of gas and liquid, different two- phase flows patterns can be found. Slug flow can be formed from the stratified flow, due to growth of hydrodynamic instabilities or due to undulations of the surface where the pipes are installed. The present work consists in the analysis of two-phase flow in slug flow through horizontal pipes using the two-fluid model in its transient and one-dimensional form. By means of a Kelvin-Helmholtz stability analysis for stratified flows, the operation range for obtaining of a well-posed mathematical model is estimated. Different techniques of discretization were implemented and tested. Comparisons with result obtained in the literature were done.

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.

Made available in DSpace on 2014-06-11T19:23:38Z (GMT). No. of bitstreams: 0 Previous issue date: 2007-12-27Bitstream added on 2014-06-13T18:50:59Z : No. of bitstreams: 1 mendes_faa_me_ilha.pdf: 2673818 bytes, checksum: 120558e9283b734ad04dda5b07ec06d0 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / O princípio de funcionamento de um sistema gas-lift é bastante simples, baseando-se na injeção de um gás inerte próximo à extremidade submersa de um tubo utilizado para ascensão do líquido (riser). Esse sistema tem sido utilizado em diversos setores da indústria, dentre os quais se destaca o setor petroquímico. Apesar de ter sido concebido há mais de um século, e a despeito do considerável número de artigos publicados sobre o assunto, vários aspectos relativos às características do escoamento bifásico dentro do riser e ao funcionamento do sistema são, ainda, pouco conhecidos. No presente trabalho, um sistema air-lift de pequeno porte foi inteiramente desenvolvido e testado, permitindo observar o comportamento do sistema quando submetido a variações na vazão de ar, na razão de submersão e na geometria do injetor – contendo três (3F) ou quinze furos (15F). Ensaios de visualização foram conduzidos, possibilitando a caracterização do escoamento bifásico no interior do riser. Medições experimentais foram, também, realizadas para a obtenção da eficiência e das curvas características do sistema. Os resultados quantitativos foram interpretados em associação com imagens estáticas do escoamento, capturadas com uma câmera fotográfica digital, e com imagens dinâmicas, tomadas com o auxílio de uma filmadora de alta velocidade. A partir destas imagens, foi possível, ainda, a estimativa de alguns importantes parâmetros característicos do escoamento pistonado, associados ao deslocamento da chamada bolha de Taylor. Os dados assim obtidos foram comparados com valores oriundos de correlações empíricas propostas por outros autores, apresentando boa concordância. Com relação aos resultados quantitativos, observou-se que, sob determinadas condições de operação, o injetor 3F é capaz de atribuir ao sistema um rendimento significativamente maior... / The operation principle of a gas-lift system is very simple: it is based on the injection of an inert gas close to the immersed edge of a tube, which is used for the rising of the liquid (riser). This system has been used in several industrial areas and more often in the petrochemical industry. Despite the fact that it was conceived over a century ago and that there have been a considerable amount of articles published on the subject, several aspects concerning the two-phase flow characteristics inside the riser as well as its operational system are still not very well known. In the present work, a small size air-lift system was thoroughly developed and tested, which facilitated the visualization of the system behavior when submitted to variations in the air flow, the immersion ratio and the injector geometry – containing three (3H) or fifteen holes (15H). Observation experiments were carried out making it possible to characterize the two-phase flow inside the riser. Experimental measurements were also taken in order to obtain the system characteristic efficiency and slopes. The quantitative results were analyzed together with the statistic images of the flow - taken by a digital photo camera, and with the dynamic images – taken by a high-speed filming camera. Through these images, it was also possible to estimate some important characteristic parameters of the slug flow associated to the displacement of the so-called Taylor bubble. The data obtained this way were compared to the figures originated from experimental correlations proposed by other authors, and they showed a fine agreement. Regarding the quantitative results, it was observed that under certain operation conditions the 3F injector is able to give the system a significantly greater efficiency than when the 15F injector is used

Made available in DSpace on 2014-06-11T19:23:39Z (GMT). No. of bitstreams: 0 Previous issue date: 2006-11-27Bitstream added on 2014-06-13T18:10:04Z : No. of bitstreams: 1 ferreira_ad_me_ilha.pdf: 1079421 bytes, checksum: dcb82443da323ec46bf3bd207a17001c (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Devido à solubilidade mútua entre o fluido refrigerante e o fluido lubrificante usados em sistemas de refrigeração por compressão de vapor, eles formam uma mistura homogênea que influencia tanto os processos de transferência de calor no evaporador e no condensador, como os processos de lubrificação e de selagem de vazamentos no interior do compressor. O vazamento de refrigerante através da folga radial de compressores rotativos de pistão rolante é de particular importância para o bom desempenho do compressor, uma vez que ele influencia significativamente a eficiência volumétrica do compressor, chegando a somar cerca de 30% das perdas totais de refrigerante. No presente trabalho foram desenvolvidos modelos de escoamento de misturas óleo-refrigerante através desta folga, incluindo a mudança de fase do refrigerante devida à variação da sua solubilidade no lubrificante. A solução da equação da energia constitui uma evolução do processo de modelagem deste escoamento em relação aos modelos até então desenvolvidos. Quatro modelos diferentes foram usados para simular o escoamento: modelo de escoamento bifásico homogêneo isotérmico, modelo de escoamento bifásico isotérmico com formação de espuma, modelo de escoamento bifásico homogêneo não-isotérmico e modelo de escoamento bifásico homogêneo não-isotérmico com termo de força inercial. O estudo foi realizado para três misturas óleo-refrigerante: óleo éster Freol a10 e refrigerante R134a, óleo éster EMKARATE RL10H e refrigerante R134a e óleo mineral SUNISO 1 GS e refrigerante R12. Para todos os modelos e misturas, realizou-se um estudo paramétrico envolvendo as principais variáveis do problema: pressão de entrada, temperatura de entrada, vazão de mistura e valor da folga mínima... / Due to the mutual solubility between the refrigerant and lubricant of refrigeration systems using mechanical compression of vapor, they form a homogeneous mixture which influences the heat transfer processes in the evaporator and condenser as well as the compressor lubrication and refrigerant leakage. The refrigerant leakage through the radial clearance of rolling piston compressors plays an important role to the volumetric efficiency in this type of compressor, in which it represents about 30% of the total refrigerant loss. In the present work several models to predict the lubricant-refrigerant mixture flow through this clearance, including the refrigerant phase change due to the reduction of the refrigerant solubility in the lubricant, are developed. Four different models were developed to simulate the flow: isothermal homogeneous two-phase flow, isothermal two-phase flow with foam formation, non-isothermal homogeneous two-phase flow and non-isothermal homogeneous two-phase flow containing inertial force. The simulations were performed for three mixtures: ester oil Freol a10-refrigerant R134a, ester oil EMKARATE RL10H-refrigerant R134a, and mineral oil SUNISO 1 GS-refrigerant R12. The influences of the inlet pressure, inlet temperature, mixture mass flow rate, and minimal clearance were analyzed for all mixtures. The results showed that it is important to consider the foam formation, the inertial force, and the solution of the energy equation in the flow modeling. Concerning the volumetric efficiency of the compressor the ester oil Freol a10-R134a was the best mixture because it produced the lowest refrigerant leakage.

Submitted by Maike Costa (maiksebas@gmail.com) on 2017-05-24T14:26:07Z No. of bitstreams: 1 arquivototal.pdf: 2237346 bytes, checksum: 088238fab84060042b385a40d43c4c3b (MD5) / Made available in DSpace on 2017-05-24T14:26:07Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 2237346 bytes, checksum: 088238fab84060042b385a40d43c4c3b (MD5) Previous issue date: 2016-02-04 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / The present paper makes an analysis of multiphase flow up (oil-gas) in laminar regime, with hidrodynamically developed profiles and interacting heat. The fluids analyzed are considered immiscible, developing into separated speed and temperature fields. The energy equations are applied to obtain of heat fields in the thermal entry region of the fluids flowing in the same direction. Although the results of parameters of practical interest has been obtained by purely numerical procedures, the system governing equations involving partial differential and have not been resolved successfully analytically yet. Thus, the use of Generalizes Integral Transform Technique, as powerful mathematical tool, is proposed to solve the energy equations associated with the flow and to producing hybrid results, analytical-numerical, more extension e analysis flexibility to a large number of parameter used in oil and gas engineering. Auxiliary problems with the greatest possible number of original problem information are presented in an attempt of improving the numerical convergence. Lastly, a computational routine is developed in FORTRAN to obtain the numerical solution. The numerical results obtained for the fields are presented in tables and graphs for then be compared with the maximum care to the parameters of special cases available in the literature. / O presente trabalho realiza uma análise do escoamento multifásico (óleo-gás) ascendente em regime laminar, com perfis hidrodinamicamente desenvolvidos e interagindo termicamente. Os fluidos em análise são considerados imiscíveis, desenvolvendo-se em campos separados de velocidade e temperatura. A equação da energia é aplicada para a obtenção dos campos térmicos em desenvolvimento na região de entrada térmica dos fluidos escoando na mesma direção. Embora resultados de parâmetros de interesse prático já tenham sido obtidos através de procedimentos puramente numéricos, o sistema de equações que regem esse tipo de escoamento envolve diferenciais parciais e ainda não foram solucionadas analiticamente com sucesso. Dessa forma, o uso da Técnica de Transformada Integral Generalizada, como ferramenta matemática, é proposto para resolver as equações de energia associada ao escoamento produzindo resultados híbridos, analítico-numérico, com maior flexibilidade de análise e extensão para um grande número de parâmetros utilizados na engenharia de petróleo e gás. Problemas auxiliares com o maior número possível de informações do problema original são apresentados na busca da melhoria da convergência numérica. Por fim, uma rotina computacional é elaborada em FORTRAN para a obtenção da solução numérica. Os resultados numéricos obtidos para os campos são apresentados na forma de tabelas e gráficos para em seguida serem comparados com o máximo de esmero aos parâmetros de casos especiais disponíveis na literatura especializada.

Resumo: O princípio de funcionamento de um sistema gas-lift é bastante simples, baseando-se na injeção de um gás inerte próximo à extremidade submersa de um tubo utilizado para ascensão do líquido (riser). Esse sistema tem sido utilizado em diversos setores da indústria, dentre os quais se destaca o setor petroquímico. Apesar de ter sido concebido há mais de um século, e a despeito do considerável número de artigos publicados sobre o assunto, vários aspectos relativos às características do escoamento bifásico dentro do riser e ao funcionamento do sistema são, ainda, pouco conhecidos. No presente trabalho, um sistema air-lift de pequeno porte foi inteiramente desenvolvido e testado, permitindo observar o comportamento do sistema quando submetido a variações na vazão de ar, na razão de submersão e na geometria do injetor - contendo três (3F) ou quinze furos (15F). Ensaios de visualização foram conduzidos, possibilitando a caracterização do escoamento bifásico no interior do riser. Medições experimentais foram, também, realizadas para a obtenção da eficiência e das curvas características do sistema. Os resultados quantitativos foram interpretados em associação com imagens estáticas do escoamento, capturadas com uma câmera fotográfica digital, e com imagens dinâmicas, tomadas com o auxílio de uma filmadora de alta velocidade. A partir destas imagens, foi possível, ainda, a estimativa de alguns importantes parâmetros característicos do escoamento pistonado, associados ao deslocamento da chamada bolha de Taylor. Os dados assim obtidos foram comparados com valores oriundos de correlações empíricas propostas por outros autores, apresentando boa concordância. Com relação aos resultados quantitativos, observou-se que, sob determinadas condições de operação, o injetor 3F é capaz de atribuir ao sistema um rendimento significativamente maior... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The operation principle of a gas-lift system is very simple: it is based on the injection of an inert gas close to the immersed edge of a tube, which is used for the rising of the liquid (riser). This system has been used in several industrial areas and more often in the petrochemical industry. Despite the fact that it was conceived over a century ago and that there have been a considerable amount of articles published on the subject, several aspects concerning the two-phase flow characteristics inside the riser as well as its operational system are still not very well known. In the present work, a small size air-lift system was thoroughly developed and tested, which facilitated the visualization of the system behavior when submitted to variations in the air flow, the immersion ratio and the injector geometry - containing three (3H) or fifteen holes (15H). Observation experiments were carried out making it possible to characterize the two-phase flow inside the riser. Experimental measurements were also taken in order to obtain the system characteristic efficiency and slopes. The quantitative results were analyzed together with the statistic images of the flow - taken by a digital photo camera, and with the dynamic images - taken by a high-speed filming camera. Through these images, it was also possible to estimate some important characteristic parameters of the slug flow associated to the displacement of the so-called Taylor bubble. The data obtained this way were compared to the figures originated from experimental correlations proposed by other authors, and they showed a fine agreement. Regarding the quantitative results, it was observed that under certain operation conditions the 3F injector is able to give the system a significantly greater efficiency than when the 15F injector is used / Orientador: Sérgio Said Mansur / Coorientador: Edson Del Rio Vieira / Banca: André Luiz Seixlack / Banca: Ricardo Augusto Mazza / Mestre