Simulação numérica de dutos enterrados, submetidos à perda de apoio e elevação localizada.

Os recentes desastres ocorridos no país, como o rompimento da adutora em Campo Grande e os desastres relacionados às enchentes urbanas, mostram a necessidade de desenvolvimento de pesquisas científicas que auxiliem na compreensão e no dimensionamento das estruturas projetadas para atender a demanda da população. Os métodos analíticos e experimentais mais utilizados possuem algumas limitações de ordem teórica ou prática. Por outro lado, os métodos numéricos, capazes de simular etapas construtivas e envolver materiais com diferentes modelos constitutivos numa mesma análise, buscam atender às necessidades práticas dos projetos de geotecnia e, ao mesmo tempo, complementam os modelos analíticos e experimentais. Nesse trabalho foram realizadas comparações entre resultados obtidos em ensaios experimentais e resultados extraídos do modelo computacional, buscando aumentar a compreensão sobre a interação solo-estrutura em relação à distribuição de tensões mobilizadas e aos deslocamentos e deformações provocados. A simulação numérica foi feita com a utilização do PLAXIS/3D, software de análise geotécnica baseado no método dos elementos finitos. Os ensaios foram confeccionados na Escola de Engenharia de São Carlos/USP por Costa (2005) e envolveram dutos enterrados submetidos à perda de apoio ou elevação localizada. O estudo experimental foi realizado através de modelos físicos compostos por um maciço de areia pura, contendo um tubo repousando sobre um alçapão no centro do vão. Os modelos físicos foram equipados com instrumental capaz de medir as deflexões e as deformações específicas ao longo do duto, além das tensões totais no maciço de solo circundante e na base do equipamento.

Análise numérica de dutos enterrados sujeitos a elevação e recalque localizado

Problems associated to longitudinal interactions in buried pipelines are characterized as three-dimensional and can lead to different soil-pipe issues. Despite the progress achieved in research on buried pipelines, little attention has been given to the three-dimensional nature of the problem throughout the last decades. Most of researches simplify the problem by considering it in plane strain condition. This dissertation aims to present a study on the behavior of buried pipelines under local settlement or elevation, using three-dimensional simulations. Finite element code Plaxis 3D was used for the simulations. Particular aspects of the numerical modeling were evaluated and parametric analyzes were performed, was investigated the effects of soil arching in three-dimensional form. The main variables investigated were as follows: relative density, displacement of the elevation or settlement zone, elevated zone size, height of soil cover and pipe diameter/thickness ratio. The simulations were performed in two stages. The first stage was involved the validation of the numerical analysis using the physical models put forward by Costa (2005). In the second stage, numerical analyzes of a full-scale pipeline subjected to a localized elevation were performed. The obtained results allowed a detailed evaluation of the redistribution of stresses in the soil mass and the deflections along the pipe. It was observed the reduction of stresses in the soil mass and pipe deflections when the height of soil cover was decreased on regions of the pipe subjected to elevation. It was also shown for the analyzed situation that longitudinal thrusts were higher than vi circumferential trusts and exceeded the allowable stresses and deflections. Furthermore, the benefits of minimizing stress with technical as the false trench, compressible cradle and a combination of both applied to the simulated pipeline were verified

Análise de critérios de dimensionamento de carga em dutos de concreto instalados em vala na situação de recobrimento mínimo.

Perante diversas situações da engenharia são utilizadas formulações empíricas de dimensionamento baseadas em dados de campo e experiência profissional que definem muito o caráter subjetivo da metodologia padrão de projeto. O presente trabalho de pesquisa aborda os diversos métodos de obtenção dos esforços gerados em dutos enterrados submetidos a cargas dinâmicas e estáticas e sua posterior reavaliação através de modelagem numérica com o programa Plaxis 3D. Os métodos analíticos não convencionais foram comparados com o método padrão de cálculo sendo que o mesmo demonstrou ter uma boa precisão mesmo sem considerar outros fatores importantes como a parcela de resistência devida à coesão do solo e sua deformabilidade. A modelagem numérica demonstrou o conservadorismo do método de Marston e o subdmensionamento do espraiamento em prisma devido aos efeitos locais ocasionados pela adoção do recobrimento mínimo e sobrecarga dinâmica elevada. Também se observou, através da modelagem 3D, que a utilização dos dois métodos clássicos favorecem a obtenção de resultados dentro da razoabilidade.Verificou-se também, como resultado desta pesquisa, que a proposta de um método clássico modificado permite uma melhor aproximação da carga que atinge o duto.

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.

Dois testes de imageamento com GPR em problemas de controle ambiental em regiões tropicais: migração de dunas e localização de dutos de óleo enterrados

Because the penetration depth of Ground Penetrating Radar (GPR) signals is very limited in high conductive soils, the usefullness of this method in tropical regions is not yet completly known. The main objective of this researh is to test the usefullness of the method in Brazil. Two typical problems where GPR has been used in Europe and North American were choosed for this test: the first one is to characterize the internal structures of a sand body and the second problem is the localization of old buried pipes lines. The first test was done near the city of São Bento do Norte, in the northern coast of Rio Grande do Norte state, NE Brazil. In this region, there is a sand dune that is migrating very fast in the direction of adjacent settling areas. To characterize the internal structure of the dune and its relationship to the prevailing wind direction, as a preliminary step to understand the dune migration, GPR profiles using the 400 MHz frequency were performed in E-W, N-S, NE-SW, and SE-NW directions over the sand dune intersecting at the top of the dune. The practical resolution of the GPR data is around 30 cm; this was sufficient to distinguish individual foresets inside the dune. After applying the elevation correction to the data, we identified that dips of bedding structures are smallest for the N-S profile, which is perpendicular to the dominant wind direction, largest for the E-W profile, and intermediate for the SW-NE and SE-NW profiles. Foresets in the E-W profile dip with angles varying from 2 to 6 degrees. In the E-W profile, the water table and a horizontal truncation interface separating two generations of dunes were identified, as well as an abrupt directional change in the foreset patterns associated to a lateral contact between two dune generations, the older one extending to the west. The used high frequency of 400 Mhz does not allow a penetration deep enough to map completely these internal contacts. The second test was done near Estreito, a small town near Carnaúbais city, also in Rio Grande do Norte state. In this locality, there are several old pipe lines buried in area covered by plantations where digging should be minimized. Several GPR profiles using the 400 and 200 MHz frequency were performed trying to intercept perpendicularly the possible pipe lines. Because of the high conductivity of the soil, the raw original data can hardly be use to identify the pipe lines. However, after an adequate processing over the 200 MHz profiles, six pipe lines were identified. As a global result of the tests, GPR can be very usefull if the conductivity of the ground is low or, in the case of medium conductivities of the soils, if adequate processing is performed

Avaliação por pig de perfilagem de danos superfciais nos materiais das paredes de dutos de petróleo

The modern technology of materials and structural integrity of pipelines requests the use of inspection tools named inspection pigs to detect, localize and measure the length, width and depth dimensions of the thickness losses of walls of buried and underwater pipelines in service. These tools run them internally, performing and recording measurements, with performance that varies according to the pig s technology. It has been developed recently an instrumented pig technology, called feller pig. This work aims to indicate factors that influence the feller pig technology performance in the detection and in the accuracy of measurement of the length, width and depth dimensions of the thickness losses on the internal surface of an oil pipeline wall under normal conditions of oil pipe inspection with pig. In this work, is made a collection of factors and an analyses of the technology based on the available literature, as well as an experiment to observe the technology and the factors operating. In the experiment, a feeler pig is used in a pipeline built in carbon steel and in operation that flows petroleum, in witch are observed areas with internal thickness losses occurred naturally. Some of these areas and their dimensions taken by automated ultra-sound scanner are compared with the ones indicated by the feller pig. Based on the data collection, on the analysis and on the experiment, the influence of factors object of this research is discussed. It is concluded that, among these, there are factors related to pipe fabrication tolerances, to wear of pig components, to internal adhesive wear of pipeline, to other pipeline damages and to technology characteristics. Finally, actions are suggested to know better, improve and define the applicability of this technology

Desenvolvimento e implementação de algoritmo para análise de dutos fabricados ou reparados com material compósito carregados com pressão interna utilizando o método dos elementos finitos.

A integridade de dutos ganha importância à medida em que o desenvolvimento da indústria expande a malha para transportar líquidos e gases num contexto de aumento das exigências ambientais. Um importante aliado para manutenção da integridade de dutos são reparos de materiais compósitos. Estes materiais apresentam baixa densidade, capacidade de direcionar resistência de acordo com as direções das solicitações, execução de reparo a frio sem necessidade de interromper produção ou grande maquinário. Este trabalho tem como objetivo desenvolver e implementar um algoritmo de elementos finitos que permita avaliar os esforços e a resistência das paredes de um tubos fabricados ou reparados com laminados de material compósito carregados com pressão interna. Entre as vantagens de desenvolver um programa tem-se: agilidade de avaliação, menor custo com licença, menores exigências computacionais, possibilidade de desenvolver o programa e o melhor entendimento da modelagem dos fenômenos. Utiliza-se como entrada do programa o diâmetro do duto, pressão interna e parâmetros do laminado. A modelagem em elementos finitos é realizada a partir da teoria clássica de laminados. Aplicando o carregamento resultante da pressão interna, determina-se os deslocamentos e são calculadas as tensões e aplicado o critério de falha de Tsai-Hill em cada camada. Estudos experimentais e numéricos encontrados na literatura foram simulados com o programa gerado e os resultados para propriedades do laminado, tensões nos dutos e pressão de ruptura apresentam concordância com os resultados da literatura.O programa ainda tem sua estrutura modificada para encontrar a pressão de falha a partir dos dados do laminado. O programa implementado permite uma avaliação rápida de resistência do reparo e possibilita avaliar rapidamente a resposta a mudanças nos parâmetros de projeto do laminado.

Simulação de escoamentos incompressíveis não newtonianos em dutos com expansão brusca

O objetivo deste trabalho é a simulação numérica de escoamentos incompressíveis bidimensionais em dutos com expansão brusca, considerando o raio de expansão de 3 : 1. As equações governantes implementadas são as de Navier, que junto com relações constitutivas para a tensão visam representar comportamentos não newtonianos. A integração temporal é feita usando o esquema explícito de Runge-Kutta com três estágios e de segunda ordem; as derivadas espaciais são aproximadas pelo método de diferenças finitas centrais. Escoamentos em expansões bruscas para fluidos newtonianos apresentam um número de Reynolds crítico, dependente do raio de expansão, na qual três soluções passam a ser encontradas: uma solução sim étrica instável e duas soluções assimétricas rebatidas estáveis. Aumentando o número de Reynolds, a solução passa a ser tridimensional e dependente do tempo. Dessa forma, o objetivo é encontrar as diferenças que ocorrem no comportamento do fluxo quando o fluido utilizado possui características não newtonianas. As relações constitutivas empregadas pertencem à classe de fluidos newtonianos generalizados: power-law, Bingham e Herschel-Bulkley. Esses modelos prevêem comportamentos pseudoplásticos e dilatantes, plásticos e viscoplásticos, respectivamente. Os resultados numéricos mostram diferenças entre as soluções newtonianas e não newtonianas para Reynolds variando de 30 a 300. Os valores de Reynolds críticos para o modelo power-law não apresentaram grandes diferenças em comparação com os da solução newtoniana. Algumas variações foram percebidas nos perfis de velocidade. Entretanto, os resultados obtidos com os modelos de Bingham e Herschel-Bulkley apresentaram diferenças significativas quando comparados com os newtonianos com o aumento do parâmetro adimensional Bingham; à medida que Bingham é aumentado, o tamanho dos vórtices diminui. Além disso, os perfis de velocidade apresentam diferenças relevantes, uma vez que o fluxo possui regiões onde o fluido se comporta como sólido.

Transferência de calor combinando radiação e convecção no interior de dutos de geradores de vapor fumotubulares

Esta dissertação de mestrado considera a transferência de calor combinando convecção e radiação térmica no escoamento de gases participantes em dutos de seção circular. Partindo de uma metodologia geral, o trabalho enfoca principalmente os casos típicos de aplicação em geradores de vapor fumotubulares de pequeno e médio porte, em que gases em alta temperatura escoam através de um tubo mantido em temperatura uniforme. O escoamento é turbulento e o perfil de velocidade é plenamente desenvolvido desde a entrada do duto. A temperatura do gás, contudo, é uniforme na entrada, considerando-se a região de desenvolvimento térmico. Duas misturas de gases são tratadas, ambas constituídas por dióxido de carbono, vapor d’água e nitrogênio, correspondendo a produtos típicos da combustão estequiométrica de óleo combustível e metano. As propriedades físicas dos gases são admitidas uniformes em todo o duto e calculadas na temperatura de mistura média, enquanto que as propriedades radiantes são modeladas pela soma-ponderada-de-gases-cinzas. O campo de temperatura do gás é obtido a partir da solução da equação bidimensional da conservação da energia, sendo os termos advectivos discretizados através do método de volumes de controle com a função de interpolação Flux-Spline; as trocas de energia radiantes são avaliadas por meio do método das zonas, onde cada zona de radiação corresponde a um volume de controle. Em um primeiro passo, a metodologia é verificada pela comparação com resultados apresentados na literatura para a transferência de calor envolvendo apenas convecção e combinando convecção com radiação. Em seguida, discutem-se alguns efeitos da inclusão da radiação térmica, por exemplo, no número de Nusselt convectivo e na temperatura de mistura do gás. Finalmente, são propostas correlações para o número de Nusselt total, que leva em conta tanto a radiação quanto a convecção. Essa etapa exige inicialmente uma análise dos grupos adimensionais que governam o processo radiante para redução do número elevado de parâmetros independentes. As correlações, aplicáveis a situações encontradas em geradores de vapor fumotubulares de pequeno e médio porte, são validadas estatisticamente pela comparação com os resultados obtidos pela solução numérica.

Avaliação de heterociclo do tipo mesoiônico solubilizado em sistema microemulsionado para aplicação em dutos

The efficiency of inhibition to corrosion of steel AISI 1018 of surfactant coconut oil saponified (SCO) and heterocyclic type mesoionics (1,3,4-triazólio-2-tiolato) in systems microemulsionados (SCO-ME and SCO-ME-MI) Of type O/A (rich in water emulsion) region with the work of Winsor IV. The systems microemulsionados (SCO-ME and SCO-ME-MI) were evaluated with a corrosion inhibitor for use in saline 10,000 ppm of chloride enriched with carbon dioxide (CO2). The assessment of corrosion inhibitors were evaluated by the techniques of linear polarization resistance (LPR) and loss of weight (MW) in a cell instrumented given the gravity and electrochemical devices. The systems were shooting speed of less than 60 minutes and efficiency of inhibition [SCO-ME (91.25%) and SCO-ME-MI (98.54%)]

Implementação de um Módulo de Supervisão para um Sistema de Detecção de Vazamentos em Dutos de Petróleo

The transport of fluids through pipes is used in the oil industry, being the pipelines an important link in the logistics flow of fluids. However, the pipelines suffer deterioration in their walls caused by several factors which may cause loss of fluids to the environment, justifying the investment in techniques and methods of leak detection to minimize fluid loss and environmental damage. This work presents the development of a supervisory module in order to inform to the operator the leakage in the pipeline monitored in the shortest time possible, in order that the operator log procedure that entails the end of the leak. This module is a component of a system designed to detect leaks in oil pipelines using sonic technology, wavelets and neural networks. The plant used in the development and testing of the module presented here was the system of tanks of LAMP, and its LAN, as monitoring network. The proposal consists of, basically, two stages. Initially, assess the performance of the communication infrastructure of the supervisory module. Later, simulate leaks so that the DSP sends information to the supervisory performs the calculation of the location of leaks and indicate to which sensor the leak is closer, and using the system of tanks of LAMP, capture the pressure in the pipeline monitored by piezoresistive sensors, this information being processed by the DSP and sent to the supervisory to be presented to the user in real time

Avaliação da corrosão em dutos por técnica gravimétrica e de resistência elétrica

With the increasing of demand for natural gas and the consequent growth of the pipeline networks, besides the importance of transport and transfer of oil products by pipeline, and when it comes to product quality and integrity of the pipeline there is an important role regarding to the monitoring internal corrosion of the pipe. This study aims to assess corrosion in three pipeline that operate with different products, using gravimetric techniques and electrical resistance. Chemical analysis of residues originated in the pipeline helps to identify the mechanism corrosive process. The internal monitoring of the corrosion in the pipelines was carried out between 2009 and 2010 using coupon weight loss and electrical resistance probe. Physico-chemical techniques of diffraction and fluorescence X-rays were used to characterize the products of corrosion of the pipelines. The corrosion rate by weight loss was analyzed for every pipeline, only those ones that has revealed corrosive attack were analyzed located corrosion rate. The corrosion potential was classified as low to pipeline gas and ranged from low to severe for oil pipelines and the pipeline derivatives. Corrosion products were identified as iron carbonate, iron oxide and iron sulfide

Sistema inteligente para detecção de vazamentos em dutos de petróleo usando transformada Wavelet e redes neurais

This work consists in the use of techniques of signals processing and artificial neural networks to identify leaks in pipes with multiphase flow. In the traditional methods of leak detection exists a great difficulty to mount a profile, that is adjusted to the found in real conditions of the oil transport. These difficult conditions go since the unevenly soil that cause columns or vacuum throughout pipelines until the presence of multiphases like water, gas and oil; plus other components as sand, which use to produce discontinuous flow off and diverse variations. To attenuate these difficulties, the transform wavelet was used to map the signal pressure in different resolution plan allowing the extraction of descriptors that identify leaks patterns and with then to provide training for the neural network to learning of how to classify this pattern and report whenever this characterize leaks. During the tests were used transient and regime signals and pipelines with punctures with size variations from ½' to 1' of diameter to simulate leaks and between Upanema and Estreito B, of the UN-RNCE of the Petrobras, where it was possible to detect leaks. The results show that the proposed descriptors considered, based in statistical methods applied in domain transform, are sufficient to identify leaks patterns and make it possible to train the neural classifier to indicate the occurrence of pipeline leaks

Desenvolvimento e caracterização de nanocompósito de resina epóxi com nanopartículas de sílica para revestimento de dutos para transporte de petróleo

The use of polymer based coatings is a promising approach to reduce the corrosion problem in carbon steel pipes used for the transport of oil and gas in the oil industry. However, conventional polymer coatings offer limited properties, which often cannot meet design requirements for this type of application, particularly in regard to use temperature and wear resistance. Polymer nanocomposites are known to exhibit superior properties and, therefore, offer great potential for this type of application. Nevertheless, the degree of enhancement of a particular property is greatly dependent upon the matrix/nanoparticle material system used, the matrix/nanoparticle interfacial bonding and also the state of dispersion of the nanoparticle in the polymer matrix. The objective of the present research is to develop and characterize polymer based nanocomposites to be used as coatings in metallic pipelines for the transportation of oil and natural gas. Epoxy/SiO2 nanocomposites with nanoparticle contents of 2, 4, and 8 wt % were processed using a high-energy mill. Modifications of the SiO2 nanoparticles‟ surfaces with two different silane agents were carried out and their effect on the material properties were investigated. The state of dispersion of the materials processed was studied using Scanning and Transmission Electron Microscopy (SEM and TEM) micrographs. Thermogravimetric analysis (TG) were also conducted to determine the thermal stability of the nanocomposites. In addition, the processed nanocomposites were characterized by dynamic mechanical analysis (DMA) to investigate the effect of nanoparticles content and silane treatment on the viscoelastic properties and on the glass transition temperature. Finally, wear tests of the pin-on-disc type were carried out to determine the effects of the nanoparticles and the silane treatments studied. According to the results, the addition of SiO2 nanoparticles treated with silane increased the thermal stability, the storage modulus and Tg of the epoxy resin and decreased wear rate. This confirms that the interaction between the nanoparticles and the polymer chains plays a critical role on the properties of the nanocomposites

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