Modelagem dos processos envolvidos nos sistemas de secagem e armazenamento de grãos

A produção de soja é uma das principais atividades econômicas na Região Noroeste do Estado do Rio Grande do Sul. As perdas de produto em condições de comercialização ocasionadas nas atividades de secagem e armazenamento são significativas, justificando a pesquisa e aprimoramento destes processos. Nesta tese foram pesquisados dois problemas: 1. Modelamento matemático dos processos de secagem, utilizando parâmetros conhecidos de soja e 2. Modelamento matemático do problema de aeração para o cálculo da distribuição da pressão e da velocidade do ar na massa de grãos em unidades de armazenamento de soja. No problema de secagem foi desenvolvido um sistema composto de quatro equações diferenciais parciais hiperbólicas acopladas não-lineares, que descreve o comportamento da temperatura e do teor de umidade do ar e dos grãos em função do tempo. Para resolver o sistema foram utilizados os métodos das diferenças finitas (p. ex., métodos de MacCormack e Crank- Nicolson.) e o método dos volumes finitos. A análise dos resultados permitiu recomendar o método mais adequado para cada tipo do problema. Para determinação da intensidade do fluxo de massa e de calor foram utilizados os dados experimentais de camada fina obtidos da literatura e complementados com dados experimentais desta tese. Foi desenvolvido um equipamento para obtenção das curvas de secagem de grãos em secador de leito fixo, a fim de identificar o modelo para secagem em camada espessa. A comparação entre os resultados experimentais e das simulações numéricas mostrou que o modelo descreve razoavelmente a dinâmica de secagem No problema de aeração foi desenvolvido um modelo matemático que descreve o escoamento do ar em sistemas de armazenamento de grãos, baseado em relações experimentais entre velocidade e gradiente de pressão. Para resolver o problema de aeração foi utilizado o método dos elementos finitos e desenvolvido um programa computacional. Um teste realizado com o programa mostrou que os resultados da solução numérica convergem para uma solução analítica conhecida. As simulações realizadas mostraram que o programa computacional pode ser usado como instrumento auxiliar para o projeto de silos, possibilitando o cálculo e a visualização gráfica da distribuição das pressões e das linhas de corrente em diferentes seções do armazém.

Simulação numérica da combustão turbulenta de gás natural em câmara cilíndrica

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

Análise numérica de uma torre de resfriamento de grande porte

Este trabalho tem como objetivo principal determinar a troca de calor em uma torre de resfriamento de grande porte, do tipo seca. A força motriz deste sistema é o empuxo resultante da diferença de temperatura entre o ar interior e exterior da torre. A partir de uma certa velocidade o vento externo também desempenha um papel importante. A convecção natural ocorre em regime turbulento, devido às grandes dimensões. A modelagem numérica é feita pelo Método dos Volumes Finitos com malhas criadas sobre Elementos Finitos, através dos softwares comerciais Fluent e CFX. Os resultados são obtidos para diferentes potências dissipadas no trocador de calor. A partir dos resultados são realizados estudos sobre a viabilidade de mudanças na geometria da torre, a fim de se obter melhorias no escoamento, bem como uma troca de calor mais eficiente. Faz-se ainda uma regressão a fim de obter correlações entre a potência dissipada no trocador de calor e a vazão mássica de ar necessária para a troca. Estas correlaçõoes são de utilidade, para aplicações posteriores, na implantação de um simulador para o processo de troca de calor nesta torre. O trocador de calor também é enfocado a fim de se determinar o coeficiente global de troca térmica, a fim de disponibilizá-lo para aplicação em outros trabalhos

Simulação numérica da transferência simultânea de energia e umidade através do solo em um sistema trocador-armazenador de calor

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Transferência de calor e controle de temperatura em tubos capilares utilizados em eletroforese capilar

Este trabalho tem por objetivo estudar a transferência de calor em tubos capilares cilíndricos utilizados na técnica de separação de moléculas denominada Eletroforese Capilar. Esta técnica é usada, por exemplo, na análise de biomoléculas e no sequenciamento de DNA, onde o controle da temperatura está diretamente ligado ao desempenho destes métodos e à qualidade dos resultados. Para empregar esta técnica, tensões elétricas da ordem de 20 kV são aplicadas entre as extremidades dos tubos capilares, que possuem normalmente 50 cm de comprimento, 350 µm de diâmetro externo e 50 µm de diâmetro interno, preenchidos por uma solução aquosa. Tais tensões geram uma corrente elétrica na solução, provocando aquecimento distribuído por Efeito Joule. Os tubos capilares são construídos em quartzo amorfo e protegidos por uma camada de material polimérico (poliimida). Para implementar o controle da temperatura, os tubos capilares são colocados em contato com um fluido de resfriamento. Num primeiro momento, os estudos são realizados por simulação numérica, empregando o Método dos Volumes Finitos em rotinas escritas em FORTRAN. São simulados casos onde os tubos são recobertos por camadas cilíndricas de materiais com uma condutividade térmica relativamente boa, com o objetivo de aumentar a superfície de troca de calor com o fluido de resfriamento. Como resultado, obtêm-se curvas da temperatura no centro dos tubos capilares em função do coeficiente de transferência de calor por convecção. Um caso de interesse é quando os tubos capilares são posicionados excentricamente ao recobrimento cilíndrico Num segundo momento, é utilizado o software de simulação numérica ANSYS CFX®, onde é simulado o resfriamento dos mesmos tubos capilares expostos a um escoamento transversal de ar a 15°C. Neste caso, também são aplicados os recobrimentos cilíndricos e, além disso, opta-se por simular o resfriamento de um arranjo de vários tubos (sistema multicapilar) dispostos entre placas de vidro, no formato de um sanduíche. Como resultados mais importantes salientam-se: a) o aumento do raio do recobrimento resulta no aumento da transferência de calor, fazendo com que a temperatura no núcleo do capilar fique estacionada em valores baixos que não comprometem as separações/análises; b) chegou-se a um valor de raio crítico da ordem de 10 mm para a condição de operação mais típicas na área da Eletroforese Capilar; c) as montagens com o tubo capilar concêntrico e excêntrico ao recobrimento não apresentam diferenças significativas no perfil de temperatura da solução tampão; e finalmente d) observa-se que o uso de duas placas de material dielétrico com os capilares posicionados em forma de sanduíche entre elas permite uma eficiente dissipação do calor gerado na solução tampão.

Análise experimental e numérica do comportamento térmico de um coletor solar acumulador

Esta tese apresenta um estudo do comportamento térmico de um coletor solar acumulador e desenvolve uma metodologia para medir a sua eficiência diária. O coletor solar acumulador está instalado na face norte do prédio de Energia Solar da UFRGS e possui cerca de 26 m2. É constituído de uma massa espessa de concreto com uma superfície absorvente feita de tijolos, possuindo uma cobertura dupla de vidros colocada de modo a deixar um espaço para a circulação de ar. Os raios solares atravessam a cobertura de vidro e aquecem a massa absorvente de tijolo, a qual aquece o ar que é introduzido no interior da construção por efeito de termossifão. Uma das principais características do coletor solar acumulador consiste no fato de que a resposta do coletor é defasada no tempo. Este fenômeno permite que o coletor entregue calor ao ambiente mesmo após o término da radiação solar. Essa defasagem dos picos de energia térmica ocorre devido ao baixo valor da difusividade térmica do concreto. Este trabalho foi dividido em duas etapas. A primeira etapa consistiu na montagem de um calorímetro para controle e monitoração das variáveis envolvidas. No interior do calorímetro foram instaladas 36 garrafas com água. As temperaturas dos conteúdos das garrafas, do coletor solar e as radiações envolvidas foram monitoradas através de 26 sensores de temperatura de CI, 8 sensores resistivos PT100 e dois sensores de radiação fotovoltaicos. Para obter as medidas dos sensores instalados foi feita a montagem de um sistema de aquisição de dados interfaceado a um microcomputador A segunda etapa consistiu na produção de um programa computacional, escrito em linguagem Fortran 90, para simular o comportamento térmico dos diversos elementos constituintes do coletor, determinar a potência térmica do coletor solar e sua eficiência diária. Para a simulação numérica do coletor solar acumulador, adotou-se um modelo simplificado bidimensional do mesmo. Foi integrada, através do Método dos Volumes Finitos, a equação de difusão de calor transiente em 2 dimensões. Na formulação das equações lineares optou-se pelo emprego das diferenças centrais no espaço e formulação explícita no tempo. Ao todo foram produzidas 4 malhas computacionais, com distintos refinamentos e foi realizado o estudo da estabilidade numérica das diversas malhas. Através da montagem experimental obtiveram-se várias características térmicas do comportamento do sistema, entre as quais, a transmitância da cobertura, curvas de temperatura do ar fornecido ao calorímetro e curva da eficiência diária do coletor solar . Através da simulação numérica foi possível determinar a potência térmica que o coletor entrega para o laboratório, a eficiência do coletor, os campos de temperatura e a vazão mássica nos diversos canais interiores do coletor solar.

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Simulação da transferência de calor em amostras aquecidas por plasma

The processing of materials through plasma has been growing enough in the last times in several technological applications, more specifically in surfaces treatment. That growth is due, mainly, to the great applicability of plasmas as energy source, where it assumes behavior thermal, chemical and/or physical. On the other hand, the multiplicity of simultaneous physical effects (thermal, chemical and physical interactions) present in plasmas increases the complexity for understanding their interaction with solids. In that sense, as an initial step for the development of that subject, the present work treats of the computational simulation of the heating and cooling processes of steel and copper samples immersed in a plasma atmosphere, by considering two experimental geometric configurations: hollow and plane cathode. In order to reach such goal, three computational models were developed in Fortran 90 language: an one-dimensional transient model (1D, t), a two-dimensional transient model (2D, t) and a two-dimensional transient model (2D, t) which take into account the presence of a sample holder in the experimental assembly. The models were developed based on the finite volume method and, for the two-dimensional configurations, the effect of hollow cathode on the sample was considered as a lateral external heat source. The main results obtained with the three computational models, as temperature distribution and thermal gradients in the samples and in the holder, were compared with those developed by the Laboratory of Plasma, LabPlasma/UFRN, and with experiments available in the literature. The behavior showed indicates the validity of the developed codes and illustrate the need of the use of such computational tool in that process type, due to the great easiness of obtaining thermal information of interest

Simulação computacional do escoamento em bombas de cavidades progressivas

The use of Progressing Cavity Pumps (PCPs) in artificial lift applications in low deep wells is becoming more common in the oil industry, mainly, due to its ability to pump heavy oils, produce oil with large concentrations of sand, besides present high efficiency when compared to other artificial lift methods. Although this system has been widely used as an oil lift method, few investigations about its hydrodynamic behavior are presented, either experimental or numeric. Therefore, in order to increase the knowledge about the BCP operational behavior, this work presents a novel computational model for the 3-D transient flow in progressing cavity pumps, which includes the relative motion between rotor and stator, using an element based finite volume method. The model developed is able to accurately predict the volumetric efficiency and viscous looses as well as to provide detailed information of pressure and velocity fields inside the pump. In order to predict PCP performance for low viscosity fluids, advanced turbulence models were used to treat, accurately, the turbulent effects on the flow, which allowed for obtaining results consistent with experimental values encountered in literature. In addition to the 3D computational model, a simplified model was developed, based on mass balance within cavities and on simplification on the momentum equations for fully developed flow along the seal region between cavities. This simplified model, based on previous approaches encountered in literature, has the ability to predict flow rate for a given differential pressure, presenting exactness and low CPU requirements, becoming an engineering tool for quick calculations and providing adequate results, almost real-time time. The results presented in this work consider a rigid stator PCP and the models developed were validated against experimental results from open literature. The results for the 3-D model showed to be sensitive to the mesh size, such that a numerical mesh refinement study is also presented. Regarding to the simplified model, some improvements were introduced in the calculation of the friction factor, allowing the application fo the model for low viscosity fluids, which was unsuccessful in models using similar approaches, presented in previous works

Estudo numérico de escoamento turbulento em padrão anular gás-líquido em dutos verticais

Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Modelagem dinâmica do escoamento de um sistema de elevação por plunger lift

A critical problem in mature gas wells is the liquid loading. As the reservoir pressure decreases, gas superficial velocities decreases and the drag exerted on the liquid phase may become insufficient to bring all the liquid to the surface. Liquid starts to drain downward, flooding the well and increasing the backpressure which decreases the gas superficial velocity and so on. A popular method to remedy this problem is the Plunger Lift. This method consists of dropping the "plunger"to the bottom of the tubing well with the main production valve closed. When the plunger reaches the well bottom the production valve is opened and the plunger carry the liquid to the surface. However, models presented in literature for predicting the behavior in plunger lift are simplistic, in many cases static (not considering the transient effects). Therefore work presents the development and validation of a numerical algorithm to solve one-dimensional compressible in gas wells using the Finite Volume Method and PRIME techniques for treating coupling of pressure and velocity fields. The code will be then used to develop a dynamic model for the plunger lift which includes the transient compressible flow within the well

Simulação numérica do escoamento de jatos impingentes axissimétricos confinados

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Modelagem do escoamento e da transferência de calor ao longo de tubos capilares não adiabáticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Implementação e teste do método da fronteira imersa para a simulação do escoamento em torno de cilindros estacionários e rotativos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)