Ingeniero Civil Matemático / El objetivo de este trabajo es revisar la historia de un problema formulado hace ya más de 250 años, y que todavía no ha abandonado el terreno de la conjetura. De hecho, las ecuaciones de Euler y de Navier - Stokes en tres dimensiones espaciales constituyen hoy en día un desafío para matemáticos, físicos e ingenieros; aunque mucho se ha descubierto, la naturaleza de las soluciones sigue siendo un gran misterio. Precisamente, se ignora si las soluciones de la ecuación de Euler o de Navier - Stokes tridimensional (en el caso incompresible), partiendo desde condiciones iniciales regulares, mantienen esta propiedad en todo tiempo posterior, o bien desarrollan en tiempo finito una singularidad. La investigación comienza con un repaso de aquellos conceptos esenciales de la mecánica de medios continuos que se consideran indispensables para un futuro estudio de la formación de singularidades en las soluciones de la ecuación de Euler o de Navier - Stokes tridimensional. Dentro de este repaso se otorga particular atención a la ecuación de evolución de la vorticidad, una de las herramientas fundamentales en el tratamiento matemático de los fluidos, ya sean ideales o viscosos. Posteriormente se revisan los resultados clásicos concernientes a la existencia, unicidad y regularidad de soluciones de la ecuación de Euler y Navier - Stokes incompresible (en los casos bidimensional y tridimensional). A partir de estos teoremas surge naturalmente el fenómeno del quiebre, en tiempo finito, de la regularidad de dichas soluciones. Este misterio ha sido parcialmente desvelado por el ya famoso criterio de Beale - Kato - Majda, que establece que si una solución inicialmente suave de la ecuación de Euler o de Navier - Stokes 3D desarrolla una singularidad en el instante $T^{*} > 0$, entonces su campo de vorticidad $\omega(t)$ se acumula tan rápidamente en el tiempo de modo tal que: $$ \lim\limits_{t \nearrow T^{*}} \int\limits_{0}^{t} \Vert \omega(s) \Vert_{L^{\infty}} \,ds = \infty. $$ Luego de elaborar un recuento histórico sobre algunos de los intentos que han sido llevados a cabo con la intención de poner fin a esta polémica (en el caso de la ecuación de Euler 3D), en el capítulo N°6 se describe detalladamente un experimento numérico del año 2014, diseñado por Thomas Hou y Guo Luo con el propósito de hallar potenciales soluciones singulares y axisimétricas de la ecuación de Euler 3D. La principal novedad de este trabajo de memoria está en el estudio del \textit{ansatz} auto - similar propuesto por Hou y Luo para formalizar sus observaciones numéricas: se demuestra analíticamente que dicho \textit{ansatz} no conduce hacia una solución singular de la ecuación de Euler incompresible y tridimensional. El trabajo de memoria concluye con la exposición de algunos resultados que son aplicables únicamente a la ecuación de Navier - Stokes incompresible y tridimensional, tales como la estimación de Caffarelli - Kohn - Nirenberg de la medida de Hausdorff del conjunto de puntos singulares, o bien, diversos teoremas del tipo Liouville en este contexto. / Este trabajo ha sido parcialmente financiado por el proyecto Fondecyt N° 1150066

Publicación a texto completo no autorizada por el autor / Realiza un estudio de dos modelos fisicoquímicos de la dinámica de fluidos contaminantes en medios porosos a escala geológica, específicamente se estudia la interacción del modelo de flujo de aguas subterráneas y el modelo de transporte reactivo, orientados particularmente a los sistemas porosos como son los acuíferos o reservorios geológicos de aguas subterráneas. Complementando el estudio de los fenómenos en medios porosos de la presente tesis se presenta un ejemplo de simulación numérica de la ecuación del transporte de contaminantes en medios porosos. Debido la vital importancia de la simulación computacional se aborda también en el último capítulo una guía básica descriptiva de la simulación computacional de las ecuaciones de flujo y transporte. Básicamente se estudia un ejemplo de simulación computacional de la evolución espacio-temporal de una pluma contaminante reactiva inmersa en un medio poroso constituido por tres tipos de acuíferos. Esta simulación está basada en el programa computacional MODFLOW. Este programa es un modelador de flujo por diferencias finitas desarrollado por el Servicio Geológico de los Estados Unidos, el cual está basado en un código fuente que resuelve mediante interacciones la ecuación de flujo del agua subterránea. En cuanto al modelo de transporte éste se encuentra soportado por el programa MT3D, el cual es un modelo diseñado en computar las ecuaciones del transporte de solutos en 3D para la simulación de la advección, dispersión y reacciones geoquímicas de los constituyentes de solutos contaminantes disueltos en los sistemas de aguas subterráneas. / Tesis

Esta tesis establece un procedimiento de análisis de las tensiones generadas por la dilatación térmica en tres arreglos de tuberías: una tubería simple enterrada, una tubería con un cambio de dirección de 90° y una tubería enterrada con un cambio de dirección de 90°. Para darle solución a cada uno de los casos se analizan los esfuerzos térmicos por dos vías diferentes: la simulación numérica mediante el programa ANSYS y por teorías de resistencia de materiales. El planteamiento de estos dos métodos da la ventaja de corroborar los valores obtenidos verificando que sean similares entre sí y a la vez analizar los comportamientos complejos de ciertos parámetros en la tubería. Para plantear el análisis por resistencia de materiales se explica los siguientes conceptos: dilatación térmica, esfuerzo térmico, longitud virtual entre apoyos, modelo analítico del terreno, flexión en vigas, vigas sobre fundación elástica y teorías de falla. Junto con los conocimientos anteriores es posible, para cada caso, simplificar el modelo teórico y así obtener valores de tensión crítica por medio del procedimiento realizado con los conceptos estudiados. Luego se plantea el estudio por simulación numérica usando modelos lineales para las tuberías y de resortes para el terreno. Se confirma la validez de la simulación como método de análisis al contratar valores con el planteamiento teórico y obtener porcentajes de similitud del grado del 83 – 98 %. Además es posible analizar con mayor profundidad cada caso ya que la simulación da las herramientas suficientes para ver el comportamiento de la dilatación y fuerza a lo largo de la tubería y obtener los esfuerzos máximos (que varían de 35 a 166 MPa según el arreglo de tubería) para realizar un análisis de fallas.

En el presente trabajo, se estudian los efectos de diferentes fenómenos físicos sobre las frecuencias naturales de un agitador 4PBT45 (modificado en estudios preliminares para obtener un bajo consumo de potencia y alto grado de mezcla) utilizando la herramienta de simulación computacional ANSYS Mechanical, a partir de escenarios simplificados en diferentes condiciones de operación. El proyecto comienza con el análisis de los álabes aislados, para luego estudiar el comportamiento del impulsor por separado y verificar qué sucede con los modos de vibración. Se realiza la simplificación del impulsor a un rotor axisimétrico, con la misma masa y momento de inercia polar, para comparar los efectos en contraste con el agitador original, a través de escenarios estáticos y dinámicos en vacío de cada uno de los modelos. Además, se realiza el modelamiento numérico del agitador estático sumergido y se prueban diferentes propiedades de fluidos. Prácticamente, las variables que se monitorean son: cómo gobiernan las frecuencias locales de los álabes sobre el impulsor, en base al factor de participación de masa; cuáles modos de vibración se ven afectados por la rotación en vacío y cuánto varían las frecuencias debido a dicho fenómeno y por qué; cuán factible es reemplazar el impulsor por un disco y en qué casos da buenos resultados; y cuánto repercute el efecto del fluido envolvente sobre las frecuencias naturales en el sistema estático y cómo se comportan para diferentes tipos de fluidos. En general, se determina que las frecuencias naturales varían en diferentes cantidades dependiendo del fenómeno estudiado. Algunas frecuencias tienden a disminuir o a aumentar en diferentes grados conforme se incrementan las revoluciones, mientras que otras se mantienen constantes. También, se determina que casi todas las frecuencias disminuyen entre 0.5% y 20% por el efecto del fluido dependiendo principalmente de la densidad del mismo y del modo de vibrar. Por otro lado, se encuentra una expresión analítica que permite calcular la frecuencia natural bajo el efecto de un fluido B, conociendo la frecuencia dentro de un fluido A, a partir de la relación de sus densidades, bajo ciertas hipótesis y condiciones. Además, se consigue una variación entre 1% y 6% entre las frecuencias obtenidas mediante simulación numérica y las calculadas analíticamente, reduciendo el sistema a un grado de libertad, para los modos más importantes. Por último, se conocen los alcances de la herramienta computacional en este tipo de simulaciones y sus limitaciones para desarrollar escenarios más complejos.

El desarrollo de una herramienta computacional basada en el método de elementos finitos (FEM) para modelar numéricamente el proceso de transporte de pulpas minerales en tuberías circulares horizontales es el tema principal en este trabajo. Las pulpas minerales son consideradas aquí como un fluido monofásico con propiedades reológicas de fluidos no Newtonianos. La herramienta computacional es implementada siguiendo un paradigma de programación orientada a objetos, usando C++ como lenguaje de programación principal, y una metodología de descomposición de dominio para su paralelización. La referida herramienta es validada realizando simulaciones numéricas de múltiples casos de estudio involucrando flujos de fluidos Newtonianos y no Newtonianos. En esta etapa de validación, errores máximos de 3% en tuberías circulares horizontales transportando pulpas minerales reales son obtenidos. Luego de ser validada, la herramienta computacional es utilizada para estudiar parámetros importantes en el diseño de sistemas de transporte de pulpas minerales. Más específicamente, parámetros relativos a la viscosidad, grado de concentración de partículas, velocidad de la fase continua y caída de presión son particularmente analizados. El principal aporte de este trabajo es el desarrollo de una nueva herramienta computacional basada en FEM para modelar adecuadamente flujos viscosos no Newtonianos, permitiendo resolver problemas prácticos de diseño de tuberías circulares horizontales transportando pulpas minerales reales. Además, comparado con otros enfoques paralelos reportados en la literatura, el utilizado aquí es relativamente fácil de implementar en herramientas de dinámica de fluidos computacional. Así, otra de las contribuciones de este trabajo es el desarrollo de un algoritmo en paralelo adecuado para herramientas basadas en FEM que modelen pulpas minerales. Finalmente, otra de las contribuciones de este trabajo está asociada con el estudio numérico por primera vez de pulpas minerales reales recientemente caracterizadas en la literatura.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo. / The objective of this work is to develop robust and computationally efficient frameworks devoted to the solution of the problem of paraffin deposition, from the point of view of the solid-liquid equilibrium. In this scenario, several thermodynamic models will be evaluated to the description of the behavior of the liquid phase: the Peng-Robinson cubic equation of state and the activity coefficient models of Ideal Solution, Wilson, UNIQUAC and UNIFAC. The solid phase is characterized by the Multisolid model. The formation of the solid phase is initially predicted by a thermodynamic stability test. Subsequently, the nonlinear algebraic system that describes the thermodynamic equilibrium and the material balance equations is solved by three numerical approaches: the multivariable Newton method, the Broyden method and the Newton-Armijo method. Several numerical experiments were conducted in order to evaluate computation times and robustness facing several scenarios of initial estimates of the numerical methods for various models and mixtures. The results indicate the possibility of construction of efficient and robust computational frameworks which can be used coupled with fluid flow simulators, for instance.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards. / Geotechnical Engineering is the area of Civil Engineering that studies the interaction between constructions carried out by man or natural phenomena with geological environment, which most of times is partially saturated soil. In this sense, work developing as stabilization, dam containing, retaining walls, foundations and highways are conditioned to a right prediction of water flow into the soil. However, considering the water flow, the studied region areas are commonly on the order of square kilometers, mathematical models solutions require computational meshes of large proportions, causing serious limitations linked to computational memory requirements and processing time. In order to overcome these limitations, efficient numerical methods must be used in the solution of the considered problem. Hence iterative methods for solving nonlinear and large sparse linear systems must be used in this type of application. In short, this study approached a solution to the Richard partial differential equation by the two dimensions finite volume method, bringing Picard and Newton method with greater efficiency. Linear system resolution iterative techniques based on Krylov space with pre-conditioners matrix were used. Portable Extensible Toolkit for Scientific Computation (PETSc) numerical library was a tool used during the task. The results indicate when a Richards equation is solved considering thr PICARD-KRYLOV method, no matter the soil evaluation model, the best combination for solving linear systems is the stabilized double gradient method and the SOR preconditioning. On the other hand, when the van Genuchten equations are used the gradients methods with the SOR preconditioning must be chosen. Adopting the NEWTON-KRYLOV method, the stabilized double gradient method is more efficient in soling Newton linear system, in relation to the preconditioning it must be giving preference to the Jacob block. Finally, there are strong indications that the PICARDKRYLOV method can be more effective than the NEWTON-KRYLOV one, when used for solving Richards partial differential equation.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo. / The objective of this work is to develop robust and computationally efficient frameworks devoted to the solution of the problem of paraffin deposition, from the point of view of the solid-liquid equilibrium. In this scenario, several thermodynamic models will be evaluated to the description of the behavior of the liquid phase: the Peng-Robinson cubic equation of state and the activity coefficient models of Ideal Solution, Wilson, UNIQUAC and UNIFAC. The solid phase is characterized by the Multisolid model. The formation of the solid phase is initially predicted by a thermodynamic stability test. Subsequently, the nonlinear algebraic system that describes the thermodynamic equilibrium and the material balance equations is solved by three numerical approaches: the multivariable Newton method, the Broyden method and the Newton-Armijo method. Several numerical experiments were conducted in order to evaluate computation times and robustness facing several scenarios of initial estimates of the numerical methods for various models and mixtures. The results indicate the possibility of construction of efficient and robust computational frameworks which can be used coupled with fluid flow simulators, for instance.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards. / Geotechnical Engineering is the area of Civil Engineering that studies the interaction between constructions carried out by man or natural phenomena with geological environment, which most of times is partially saturated soil. In this sense, work developing as stabilization, dam containing, retaining walls, foundations and highways are conditioned to a right prediction of water flow into the soil. However, considering the water flow, the studied region areas are commonly on the order of square kilometers, mathematical models solutions require computational meshes of large proportions, causing serious limitations linked to computational memory requirements and processing time. In order to overcome these limitations, efficient numerical methods must be used in the solution of the considered problem. Hence iterative methods for solving nonlinear and large sparse linear systems must be used in this type of application. In short, this study approached a solution to the Richard partial differential equation by the two dimensions finite volume method, bringing Picard and Newton method with greater efficiency. Linear system resolution iterative techniques based on Krylov space with pre-conditioners matrix were used. Portable Extensible Toolkit for Scientific Computation (PETSc) numerical library was a tool used during the task. The results indicate when a Richards equation is solved considering thr PICARD-KRYLOV method, no matter the soil evaluation model, the best combination for solving linear systems is the stabilized double gradient method and the SOR preconditioning. On the other hand, when the van Genuchten equations are used the gradients methods with the SOR preconditioning must be chosen. Adopting the NEWTON-KRYLOV method, the stabilized double gradient method is more efficient in soling Newton linear system, in relation to the preconditioning it must be giving preference to the Jacob block. Finally, there are strong indications that the PICARDKRYLOV method can be more effective than the NEWTON-KRYLOV one, when used for solving Richards partial differential equation.

Orientador: Antonio Carlos Bannwart / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica, Instituto de Geociências / Made available in DSpace on 2018-08-26T04:52:38Z (GMT). No. of bitstreams: 1 Biazussi_JorgeLuiz_D.pdf: 11756361 bytes, checksum: e163630c268ec8b0c1e2a20bc1bdebc9 (MD5) Previous issue date: 2014 / Resumo: A bomba centrífuga submersa (BCS) é uma bomba de múltiplos estágios que tem recebido muita atenção nos últimos anos, devido à sua importância para a elevação artificial de petróleo em altas vazões. Como uma parte do sistema de Elevação Artificial, a BCS é geralmente instalada no interior do poço, a fim de aumentar a vazão ou até mesmo viabilizar a produção. A presença de gás livre no fluido na entrada da bomba provoca uma diminuição do ganho de pressão fornecido pela BCS e pode conduzir a instabilidades na curva de ganho de pressão versus vazão. O objetivo deste trabalho é desenvolver e avaliar o desempenho de um modelo de deslizamento para representar o ganho de pressão em uma ampla faixa de vazão de líquido de uma BCS. Para este propósito foram realizados testes experimentais em laboratório para três BCS diferentes que operam com misturas ar - água. Especificamente, o ganho de pressão e a potência de eixo foram medidos em diferentes vazões de líquido, de gás, de pressões de entrada e rotação. Os resultados dos testes de água monofásicos foram interpretados por ajuste dos dados com um modelo de correlação genérico para o ganho de pressão e potência, tentando descrever todos os fenômenos físicos envolvidos no escoamento nos canais rotativos. Os resultados dos testes bifásicos foram discutidos em termos dos efeitos da fração de vazio, da pressão de admissão, da velocidade de rotação e da geometria da bomba. Um modelo de deslizamento para escoamento em bolhas foi proposto para representar o ganho de pressão e também foi utilizado para expressar a potência de eixo. Os dois parâmetros envolvidos no modelo, ou seja, C0 e kb?, foram ajustados aos dados e mostrou a capacidade desta abordagem em capturar as principais tendências das curvas experimentais. O parâmetro de distribuição C0 foi, em todos os casos, significativamente menor do que um, confirmando os resultados de outros autores, em escoamento bifásico descendente em tubos. Além disso, este resultado indica que os perfis de velocidade e de distribuição de fase são distorcidos pelo campo centrífugo e de Coriolis. O parâmetro kb? foi significativo apenas para a bomba radial de menor vazão, o que confirma a ideia de que para altas vazões de líquido, as bolhas de gás são dispersadas através do líquido e o deslizamento torna-se insignificante / Abstract: The Electrical Submersible Pump (ESP) is a multiple stage pump which has been receiving a lot of attention in recent years in due to its importance for the Artificial Lift of petroleum at high flow rates. As a part of the AL system, the ESP is often installed inside the well in order to either viabilize the production or increase its flow rate. The presence of free gas within the fluid entering the pump causes a decrease in the head provided by the ESP and may lead to instabilities in the head-capacity curve. The aim of this work is to develop and evaluate the performance of a drift flux model to represent the head in the entire liquid flow rate range of an ESP. For this purpose, experimental performance tests to determine the characteristic curves were performed in laboratory for three different ESPs operating with water and water-air mixtures. Specifically, the head and brakehorsepower were measured at different flow liquid and gas flow rates, inlet pressures and rotation speeds. The results from the single-phase water tests were interpreted by fitting generic correlation models for the head and power to the data, trying to describe all basic phenomena involved in the rotating channel flow. The results from the two-phase runs were discussed in terms of the effects of the mixture composition, intake pressure, rotation speed and pump geometry. A drift flux model for bubbly flow was proposed to represent the head and also used to express the power. The two parameters involved in the model, namely C0 and kb? , were fitted to the data and showed good capability of this approach to capture the main trends of the experimental curves. The distribution parameter C0 was in all cases significantly lower than one, confirming the findings by other authors in two-phase downward pipe flow. Also, this result indicates that the velocity and phase distribution profiles are distorted by both the centripetal and Coriolis fields. The drift parameter kb? was significant for the smallest capacity radial pump only, confirming the idea that at sufficiently high liquid flow rates, the gas bubbles are dispersed through the liquid and drift becomes negligible / Doutorado / Explotação / Doutor em Ciências e Engenharia de Petróleo