Um dos atributos fundamentais associados aos escoamentos multifásicos é a existência de estruturas características segundo as quais as diferentes fases do líquido escoam. O surgimento de uma dessas estruturas, conhecidas como configurações ou regimes de escoamento, é determinado pelas vazões e propriedades físicas dos componentes, além de parâmetros geométricos como diâmetro e inclinação do conduto. O desenvolvimento de metodologias de caracterização de regimes, bem como a caracterização e o diagnóstico da transição destes regimes de escoamento são de fundamental importância. Este trabalho utiliza a análise tempo-frequência da transformada de Gabor para caracterizar os regimes de escoamento horizontais gás-líquido intermitentes. Mais especificamente, o principal objetivo está em investigar a existência de sub-regimes dentro do regime intermitente, para tanto recorremos à covariância tempo-frequência da transformada de Gabor, que é capaz de detectar transições através da não-estacionaridade associada com as correspondentes transições. Testes experimentais foram conduzidos no circuito TALC em CEA-Grenoble e uma extensiva base de dados foi obtida, cobrindo diversos tipos de escoamento intermitente. Uma sonda de condutividade elétrica, consistindo de dois anéis de eletrodos montados junto à tubulação, produziu sinais dos quais a covariância tempo-frequência foi calculada através da correspondente transformada de Gabor. / One of the main features associated to multiphase flows is the existence of characteristic dynamic structures according to which the different phases of a mixture of immiscible fluids can flow. The manifestation of one of these structures, known a flow pattern or regime, is determined by the flow rates as well as by physical and geometrical properties of the fluids and piping. The development of flow pattern characterization and diagnostic methods, and the associated transitions in between, is of crucial importance for an efficient engineering of such phenomena. Time-frequency analysis based on the Gabor transform is used in this work to characterize horizontal air-water intermittent flow regimes. More specifically, our main objective is to reveal the existence of sub-regimes inside the intermittent regimes region with the help of the corresponding time-frequency covariance based on the Gabor transform, which is capable of detecting transitions by assessing the unstationarity associated with the corresponding transitions. Experimental tests were conducted at the TALC facility at CEA-Grenoble and an extensive database was obtained, covering several types of intermittent flow. A conductivity probe, consisting in two ring electrodes flush mounted to the pipe, delivered signals from which the time-frequency covariance were calculated from the corresponding Gabor transform.

Simulações de escoamentos multifásicos são de grande interesse em aplicações práticas na indústria, em particular na indústria petrolífera, entre outras. Vários processos dependem do entendimento físico de escoamentos envolvendo iteração com partículas, sedimentação e separação de fluidos. Dos muitos métodos existentes para a simulação dos processos acima descritos, há um crescente interesse no aumento de precisão, o que levou ao desenvolvimento de estratégias que utilizam esquemas de elementos finitos discretizados em malhas dinâmicas e adaptativas, usando uma formulação ALE (do inglês, Arbitrary Lagrangian-Eulerian), juntamente com uma representação geométrica da interface. Neste sentido, este trabalho tem o objetivo de estudar e implementar estratégias robustas de controle e adaptação de malhas, em situações onde a malha dinâmica é sujeita a grandes deformações. Uma biblioteca de algoritmos e rotinas foi então desenvolvida para este fim, implementando técnicas de controle e otimização da qualidade dos elementos da malha, técnicas de adaptação da interface entre fluidos com esquemas de conservação de massa, técnicas de mudanças topológicas e preservação de propriedades materiais, além de uma comunicação facilitada destas rotinas com códigos de simulação numérica de escoamentos multifásicos existentes / Multiphase flow simulations are of great interest in practical applications, particularly in the oil industry. Several processes depend on understanding physical aspects of flows with particle interaction, sedimentation and fluid separation. Among the several existing methods to simulate the processes described above, theres a growing interest in achieving higher precision, which led to the development of strategies that use finite element discretization in adaptive, dynamic meshes, using the ALE formatulation along with a geometrical representation of the interface. In this context, this thesis aims to study and implement robust strategies for mesh adaptation, for cases where the dynamic mesh is subject to large deformations. A library of routines and algorithms was developed, implementing mesh elements control and quality optimization techniques, fluid interface adaptation techniques with a mass conservation scheme, topological modifications and material properties preservation techniques, and also a decoupled, simplified communication between these routines with existing multiphase flow numerical simulation code

Em virtude da sua facilidade de operação e elevada eficiência, sedimentadores são bastante utilizados pelas indústrias na separação de componentes com diferentes densidades. Porém o desempenho do tanque sedimentador depende do campo de fluxo dentro do equipamento que, por sua vez, é influenciado pelas características dos sólidos em suspensão, da geometria e dimensões do tanque separador. As características hidrodinâmicas em um sedimentador circular vertical foram investigadas no presente trabalho, tanto por abordagens experimentais como numéricas. As experiências foram realizadas por velocimetria de imagem de partículas em um protótipo de sedimentador, construído em acrílico cristal transparente. As simulações foram realizadas, utilizando o pacote ANSYS CFX na versão 16.0, considerando escoamento 3D, turbulento, isotérmico e em regime permanente. O modelo em fluidodinâmica computacional foi criado em forma progressiva. Inicialmente, foram testados modelos para o sistema monofásico (água) e verificados através da comparação aos dados experimentais. Seis modelos de turbulência do tipo RANS foram analisados: epsilon, Renormalization Group (RNG) k-epsilon, k-omega, Shear Stress Transport (SST) e os de tensão de Reynolds (BSL RSM e SSG RSM). Para o estudo bifásico, empregou-se o modelo multifásico Euler-Euler, juntamente com o modelo de partícula e o de Gidaspow. Nesta situação, diversas configurações de condições de contorno e modificações no projeto do tanque original foram analisadas por meio de campo vetorial, linhas de fluxo, planos de concentração volumétrica e de energia cinética turbulenta. Os resultados mostraram que o modelo de turbulência BSL RSM foi o mais adequado para descrever o escoamento no sedimentador contínuo, sendo capaz de prever as recirculações e as variações de velocidade ao longo do equipamento. As zonas de recirculações (posição e intensidade) mostraram-se influenciadas pelas condições de contorno e forma geométrica do tanque de separação. Foi verificado que o aumento do diâmetro do tanque sedimentador (em 40,0%) provocou um favorecimento na separação água e sólido, elevando a eficiência em 17,3%, enquanto o aumento da altura (em 40,0%) do tanque reduziu a eficiência de separação em 27,5% em comparação ao projeto original. / Due to their operation facility and high efficiency, settlers are widely used by industries for separating components of different densities. The performance of sedimentation, however, depends on the flow field is inside the equipment, which in turn is influenced by the characteristics of the solids in suspension and by the geometry and size of the separation tank. The hydrodynamic characteristics of a vertical circular settler were investigated in the present work by means of both experimental and numeric approaches. The experimental analyses were carried out by particle image velocimetry in a settler prototype built in transparent crystal acrylic. The simulations were carried out using the software ANSYS CFX version 16.0, considering a 3-D, turbulent, isothermal and steady-state flow. The computational fluid dynamics model was used in the progressive approach. The models were initially tested for the one-phase process (water) and checked against the experimental data. Six RANS turbulence models were analyzed: k-epsilon, Renormalization Group (RNG) k- epsilon, k-omega, Shear Stress Transport (SST) and Reynolds Stress models (BSL RSM and SSG RSM). For the two-phase flow, the Euler-Euler multiphase model was used together with the particle model and the Gidaspow model. In this context, several configurations of boundary conditions and modifications in the original tank design were analyzed using vector field, flow lines, volumetric concentration and turbulent kinetic energy. The results showed that the BSL RSM turbulence model was more appropriate to explain the flow in the continuous settler, being able to predict the recirculations and velocity variations along the equipment. The recirculation zones (position and magnitude) were influenced by the boundary conditions and by the geometric shape of the separation tank. It was found that an increase in the diameter of the settler tank (in 40.0%) favored the water and solid separation, increasing its efficiency by 17.3%, while an increase in the height of the tank (in 40.0%) reduced the efficiency of separation by 27.5% when compared to the original design.

Um modelo matemático multifásico para a solidificação equiaxial de ligas binárias capaz de prever o efeito de taxa de resfriamento e da concentração de soluto no tamanho médio final de grãos foi proposto no presente trabalho. O modelo matemático foi desenvolvido através do conceito de envelope envolvendo os grãos, utilizando as equações de conservação de energia, massa e espécies químicas. O modelo de nucleação utilizado possibilita que novos núcleos possam surgir durante todo o período de resfriamento. As equações diferenciais foram obtidas através de uma média volumétrica das equações de conservação em um volume elementar representativo contendo três \"pseudofases\": sólido, líquido interdendrítico e líquido extradendrítico. O efeito de algumas variáveis de processamento sobre o tamanho médio final de grão foi analisado. Os resultados do modelo proposto foram comparados com resultados de alguns modelos disponíveis na literatura. / A multiphase mathematical model for the equiaxed solidification of binary alloys was proposed in the present work to predict the effects of the cooling rate and the average solute concentration on the final average grain size. The mathematical model was based on the concept of the grain envelope and on the conservation of energy, mass and chemical species. A nucleation model was adopted to consider the nucleation of new grains during the whole solidification time. The differential equations were derived from the volume average of conservation equations within a representative elementary volume that consisted of three pseudophases: solid, interdendritic liquid, and extradendritic liquid. The effect of some processing variables on the final average grain size was studied. The results from the proposed model were compared with those available in the literature from other models.

Este trabalho consiste em introduzir uma nova abordagem de representação de superfície no ambiente de simulação Freeflow2D. Consiste em usar Partição da Unidade Implícita para estimar da superfície a geometria, normais e curvatura. Procurando se valer das vantagens de métodos do tipo meshless (sem malha) conservando no entanto a malha Lagrangiana, no interesse de manter o fácil acesso de vizinhança, inserção e eliminação de pontos / The objective of this work is to introduce a new approache of surface representation within the Freeflow system. It consists of using implicit functions by means of Partition of Unit Implicit to estimate surface geometry, normals and curvature. Aiming at the advantages of meshless methods of surface representation whilst keeping the Lagrangian mesh in order to preserve ease of access of geometric vicinity, particle insertion and removal

Escoamentos multifásicos estão presentes em diversas aplicações industriais, principalmente na indústria do petróleo. Um dos casos de aplicação, objeto de estudo desse trabalho, é a determinação da produção de poços de petróleo através de válvulas choke. É apresentada uma revisão dos efeitos físicos e do equacionamento adotado pelos principais modelos multifásicos existentes para tais válvulas. Um estudo de sensibilidade de tais efeitos físicos é realizado, analisando as possíveis alternativas para a diferença de velocidade entre as fases, o mecanismo adotado para a troca de calor entre as fases, a influência da energia cinética a montante da válvula, a influência da área efetiva da garganta e a metodologia de cálculo da recuperação de pressão a jusante da válvula. Conclui-se que dos diversos fatores que influenciam no cálculo da vazão mássica e da condição crítica, a área efetiva da garganta é um parâmetro importante e que é necessário uma maior investigação de como determinar tal parâmetro. / Multiphase flow are present in many industrial applications, mainly at the petroleum industry. One of these application cases, aim of this work study, is to determine a petroleum well production by the choke valves. It is presented a revision of the physical effects and for the adopted equation by the main existing multiphase models for such valves. A sensibility study of such physical effects is performed, analyzing the possible alternatives for the phases velocities difference, the adopted mechanism for the heat transfer between the phases, the upstream kinetic energy influence, the throat effective area influencie and the calculation methodology adopted for the valve downstream pressure recover. It is concluded that from many factors that influence in the mass flow and critical condition calculation, the throat effective area is a important parameter and it is necessary a deeper investigation in how to determine such parameter.

The process industry is concerned with the processing of crude resources into other products. Such crudes consist of multiphase components that introduce major challenges to the operators; hence the need for efficient instrumentations that address such challenges is highly desirable. One major need is an early deposit detection system that detects deposit before it builds-up in a pipeline or equipment to prevent any possible hazard. Another critical requirement is the need to continuously monitor the flow and deduce the flow rate of every individual phase in order to study and analyse the produced product. Hence, in order to ensure safety, increase profits, optimize production and ensure production quality, the multiphase flow must be adequately monitored and controlled. This thesis demonstrated the efficiency of novel ECT algorithms for early deposit detection and multiphase flow measurement in order to measure the flow rate of all separate phases. This thesis focuses on developments in ECT image reconstruction specifically the inverse solutions and is divided into three main studies where they all build up to complete each other. In the first study, ECT is used for the first time with a narrowband pass filter to focus on targeted locations in a pipe where dielectric contaminants are expected to deposit in order to enhance the resolution of the produced images. The experimental results showed that different deposit regimes and accumulated fine deposits could be detected with high resolution. The second study allowed a better understanding of how conductive material could be imaged using a conventional ECT device and how state of the art algorithms such as iterative total variation regularisation method and the level set method could enhance this application. Also, absolute ECT imaging is presented for the first time where the level set algorithm uses only one set of ECT measurement data. This study gives a novel solution for detecting conductive deposits as well as paves the way to use the new level set algorithm for multiphase flow measurement. In the third study, the novel narrowband level set algorithm was modified to image multiphase media in order to correctly determine the number, location and concentration of the present phases. The innovative absolute ECT imaging using level set method is tested with high contrast and low contrast multiphase data, which adds more to the challenge.

Direct numerical simulations, where both the large and small scales in the flow are fully resolved, provide an excellent instrument to validate multiphase flow processes and also further our understanding of it. Three multiphase systems are studied using a finite difference/front-tracking method developed for direct numerical simulations of time-dependent system¬¬s. The purpose of these studies is to demonstrate the benefit in developing accurate sub-grid models that can be coupled with the direct numerical simulations to reduce the computational time. The primary reason to use the models is that the systems under consideration are sufficiently large that resolving the smallest scales is impractical. The processes that are examined are: (1) droplet motion and impact (2) nucleate boiling and (3) convective mass transfer. For droplet impact on solid walls and thin liquid films, the splash characteristics are studied. The collision of a fluid drop with a wall is examined and a multiscale approach is developed to compute the flow in the film between the drop and the wall. By using a semi-analytical model for the flow in the film we capture the evolution of films thinner than the grid spacing reasonably well. In the nucleate boiling simulations, the growth of a single vapor from a nucleation site and its associated dynamics are studied. The challenge here is the accurate representation of the nucleation site and the small-scale motion near the wall. To capture the evaporation of the microlayer left behind as the base of the bubble expands we use a semi-analytical model that is solved concurrently with the rest of the simulations. The heat transfer from the heated wall, the evolution of the bubble size and the departure diameter are evaluated and compared with the existing numerical results. The mass transfer near the interface, without fully resolving the layer by refining the grid is accommodated by using a boundary layer approximation to capture it. The behavior of the concentration profile is taken to be self-similar. A collection of potential profiles is tested and the accuracy of each of these models is compared with the full simulations.

"Swirling multiphase flows are found in a wide range of industrial processes. Such flows are used for separation of flows containing phases of different densities and for devices such as the spinning tensiometer. These flows are challenging to predict computationally, due to the presence of a phase boundary and the large pressure gradient generated by the swirl. In the present work the applicability of the front tracking method to swirling multi-phase flows is demonstrated by studying the evolution of a bubble in spinning tensiometer. Previous studies show that the evolution of a bubble in the spinning drop tensiometer can be used to measure the interfacial tension and other rheological properties. The front tracking method is applied to the spinning tensiometer problem to study several cases and verify the convergence of the solutions. The results are validated with other computational methods, theoretical models and experimental results. The length scales obtained from the front tracking method are in agreement with the corresponding values from experiments and other computational studies. The shape of the end of the elongated bubble obtained from the simulations is found to be similar to that suggested by a theoretical expression from previous studies. The simulations predict that the relaxation of bubble radius is exponential with time, at a rate that is found to be slightly greater than that predicted by the theoretical model."

La méthode Smoothed Particle Hydrodynamics (SPH) est une méthode lagrangienne, sans maillage développée initialement pour des simulations de phénomènes astrophysiques. Depuis, elle a connu de nombreuses applications, notamment pour la simulation des écoulements des fluides. Contrairement aux méthodes utilisant un maillage, la méthode SPH peut gérer de manière naturelle et sans traitement spécifique les simulations des écoulements à sur- face libre et multiphasiques avec interface subissant de grandes déformations. Dans cette thèse, une modélisation SPH des écoulements des fluides multiphasiques a été réalisée en tenant compte de différentes complexités (écoulements à surface libre et multiphasiques interfacials) et de natures d'écoulement (si- mulation des fluides, des sols et les deux en interactions). Un modèle SPH faiblement compressible (WCSPH) a été proposé pour simuler les écoulements des fluides multiphasiques avec interface comprenant plus de deux phases de fluide. Ce modèle inclut le développement d’une nouvelle formulation de force de tension de surface en utilisant un opérateur SPH consistant de premier ordre. Une modification de condition généralisée aux parois solides a été apportée pour qu’elle soit appliquée sur les écoulements des fluides multiphasiques avec des rapports de densité et de viscosité élevés. Une nouvelle loi de comportement dépendant de la pression nommée RBMC-αμ ( Regularized Bingham Mohr Coulomb où αμ est un paramètre libre) a également été développée. Cette loi peut simuler les fluides (Newtonien, Binghamien), les sols (cohésif, frictionnel) et les deux en interactions. La loi précédente étant sensible à la pression, une extension du terme diffusif δ-SPH a été faite pour le cas des écoulements des fluides multiphasiques afin de réduire les oscillations de pression à haute fréquence qui sont dues à l’utilisation d’une équation d’état. La validation et l’application des modèles développés dans cette thèse sont montrées à travers plusieurs cas tests de difficulté croissante. / Smoothed Particle Hydrodynamics (SPH) is a Lagrangian gridless method developed initially to simulate astrophysical phenomena, and since it has been known for a large number of applications, especially for fluid flow simulations. Contrary to the grid-based method, the SPH method can handle free surface and interfacial fluid flow simulation including large deformations naturally and without the need for any specific treatment. In this thesis a SPH modeling of multiphase fluid flows has been achieved with consideration of different complexities ( free surface and interfacial fluid flows) and natures (simulation of fluids, soil and both in interactions). A consistent weakly compressible SPH model (WCSPH) has been proposed to simulate interfacial multiphase fluid flows with more than two fluid phases. This model includes a new expression of the surface tension force using a first order consistency SPH operator. A modification to the well known generalized wall boundary condition have been brought in order to be applied to multiphase fluid flow with large density and viscosity ratios. A new pressure-based constitutive law named RBMC-αμ (Regularized Bingham Mohr Coulomb with αμ is free parameter) has been developed in this thesis. This model can simulate fluids (Newtonian, Binghamton), soils (cohesive, frictional) and both in interactions. Because the previous model is pressure sensitive, an extension of δ-SPH diffusive term has been proposed for multiphase fluid flows to overcome the hight frequency pressure oscillations due to the determination of pressure from an equation of state. The validation and application of the developed models have been shown in this thesis through several test-cases of increasing difficulty.