Made available in DSpace on 2014-08-13T14:50:37Z (GMT). No. of bitstreams: 0 Previous issue date: 2013-11-01Bitstream added on 2014-08-13T18:01:37Z : No. of bitstreams: 1 000741253_20150101.pdf: 141727 bytes, checksum: bf2e02cb75454b0d8c9536e890144bde (MD5) Bitstreams deleted on 2015-01-05T11:00:41Z: 000741253_20150101.pdf,Bitstream added on 2015-01-05T11:01:44Z : No. of bitstreams: 1 000741253.pdf: 2376805 bytes, checksum: d8aaa4bb45713ee7c22ca08294c74e07 (MD5) / O leito fluidizado é uma importante operação unitária de grande interesse na indústria química, geração de energia, extração e refino do petróleo, entre outros. As muitas aplicações do leito fluidizado na indústria e a demanda constante de melhorias na eficiência do processo de fluidização, tem produzido um aumento no número de trabalhos envolvendo simulações numéricas do escoamento gás-sólido em leito fluidizado através de CFD. Na simulação numérica, ainda não existe consenso na seleção de um modelo de arrasto e na definição de um valor adequado para o coeficiente de restituição. Nesta tese investiga-se a influência de diferentes modelos de arrasto e do coeficiente de restituição; sobre a fluidodinâmica do escoamento gás-sólido em leito fluidizado, através de simulações numéricas. Para isto é usado um modelo baseado na abordagem Euleriana-Euleriana considerando um leito fluidizado bidimensional, o qual será resolvido usando o codigo livre MFIX. Para avaliar os modelos de arrasto, os resultados numéricos são validados por meio da determinação dos dados experimentais de Taghipour et al. (2005), o qual apresenta parâmetros de projeto de leito fluidizado, tais como a queda de pressão, expansão do leito, velocidade mínima de fluidização e distribuição das fases. Estudou-se também o efeito do coeficiente de restituição sobre a fluidodinâmica do escoamento, através do desenvolvimento de uma correlação em função da fração de vazio. Os resultados demonstram que o modelo de arrasto têm uma influência muito grande na determinação das variáveis de projeto de leito fluidizado. No estudo, um grupo de modelos representam adequadamente a fluidodinâmica, sendo que o modelo de arrasto de Gidaspow (1994) foi o que mais se aproximou aos resultados experimentais. Foram desenvolvidas simulações ... ( Resumo completo, clicar acesso eletrônico abaixo) / Fluidized bed is an important unit operation of great interest in petroleum, chemical, metallurgical and energy industries. The large number of applications of fluidized bed in the industry and constant demand for improvements in the efficiency of the fluidization process, has produced an increase in the number of studies involving numerical simulations of gas-solid flow in fluidized bed through CFD. In numerical simulation, there is not yet an agreement on the selection of a drag model and in defining a suitable value for the coefficient of restitution. The purpose of this work was to study the influence of different drag models and restitution coefficient on the hydrodynamics of gas-solid flow in fluidized-bed, through numerical simulations. It is used a model based on the Eulerian-Eulerian approach considering a two-dimensional fluidized bed, which will be resolved using the open source MFIX code. To evaluate drag models, numerical results will be validated through the determination of the experimental data of Taghipour et al. (2005), which presents fluidized bed design parameters, such as fluidization curve, bed expansion, minimum fluidization velocity and void fraction profiles. It was also studied the effect of restitution coefficient on the hydrodynamics, through the development of a correlation in function of void fraction. The results demonstrated that drag model have an great influence in fluidized bed design parameters determination. In this study, a group of models were able to represent adequately the hydrodynamics, however Gidaspow (1994) drag model was the most accurate to represent the experimental data. Simulations were carried out considering restitution coefficient as a function of void fraction and using Gidaspow (1994) drag model. The simulation results obtained was able to improve the prediction of void fraction profile

[pt] A produção de óleo diminui e a produção de água aumenta com o passar do tempo na vida de um reservatório. A mistura de óleo e água é geralmente produzida na forma de uma emulsão. A formação da emulsão começa no escoamento bifásico no interior do reservatório e sua estrutura muda na medida em que os líquidos escoam através de tubulações, bombas e válvulas até as instalações de superfície. Durante a produção, as gotas maiores da fase dispersa quebram-se mudando a distribuição do tamanho das gotas. É importante conhecer a distribuição de tamanho de gota da fase dispersa da emulsão a fim de projetar as unidades de separação e prever as quedas de pressão ao longo do escoamento. O objetivo do presente trabalho é o estudo do processo de quebra de gotas em capilares retos e válvulas agulha a fim de prever o tamanho das gotas resultantes em função das condições de escoamento. O principal desafio é entender como os diferentes parâmetros operacionais de escoamento afetam o processo de quebra. Duas bancadas experimentais foram utilizadas para o estudo do processo. Na primeira, foi realizada uma análise paramétrica de formação de emulsões em um escoamento laminar através de capilares retos. Os experimentos foram realizados utilizando duas seringas conectadas através de um capilar. A emulsão óleo-emágua foi forçada a uma passagem de ida e vinda através do capilar. O diâmetro médio de gota e a superfície específica da fase dispersa foram obtidos em função da vazão de injeção, taxa de cisalhamento, tempo de residência e trabalho dissipado na parede do capilar. Como esperado, o diâmetro médio da fase dispersa diminui com o aumento do tempo de cisalhamento atingindo um valor assintótico e possui uma grande dependência com a taxa de cisalhamento na parede do capilar. Na segunda bancada, foi realizada uma análise paramétrica de um escoamento turbulento de emulsão óleo-em-água através de uma válvula agulha. Os experimentos foram realizados utilizando uma bomba helicoidal para controlar a vazão através da válvula. O diâmetro médio de gota e a superfície específica da fase dispersa a montante e a jusante da válvula foram obtidos em função da queda de pressão, vazão e taxa de dissipação de energia na válvula. O diâmetro médio da fase dispersa diminui e a superfície específica aumenta com o aumento da queda de pressão na válvula agulha até atingir um valor assintótico. / [en] Oil production decreases and water production increases as time goes by in the life of a hydrocarbon reservoir. The mixture of oil and water is usually produced as an emulsion. Emulsion formation starts in the two-phase flow inside the reservoir. The emulsion structure changes as it flows through pipes, pumps and valves up to the surface facilities. During all stages, large drops of the dispersed phase break up leading to smaller drops. It is important to know the droplet size distribution of the dispersed phase in order to design separation units and predict the pressure drop along the flow. The aim of the this work is to study the droplets break-up process that takes place in capillaries and in a needle valve in order to make predictions of the size of the resulting droplets that emerge from this process. The main challenge is to understand how the different operating flow parameters affect the break up process. In order to achieve this goal, two laboratory scale experimental set-ups have been used. In the first experiment, we conducted a parametric analysis of oil-water emulsion formation in laminar flow through straight capillaries. The experiments were carried out using two syringe pumps connected by a double-hubbed capillary pipe. The oil-water emulsion is forced back and forth through the pipe. The mean diameter and the specific surface area of the dispersed phase were obtained as a function of flow rate, shear rate, residence time and rate of energy dissipation at the capillary wall. As expected, keeping all other variables fixed, the dispersed phase mean diameter decreases with the shearing time, reaching an asymptotic value, which was a strong function of the shear rate at the capillary wall. Secondly, we conducted a parametric analysis of turbulent oil-in-water emulsion flow through a needle valve. The experiments were carried out using a helicoidal pump to control the flow rate through the needle valve. The mean diameter and the specific surface area of the dispersed phase upstream and downstream of the valve were obtained as a function of the pressure drop in the valve, flow rate, and rate energy dissipation of the flow. The dispersed phase mean diameter falls and the specific surface area rises with the pressure drop in the valve until reaching an asymptotic value.

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Química, Florianópolis, 2013. / Made available in DSpace on 2014-08-06T17:09:28Z (GMT). No. of bitstreams: 1 323813.pdf: 923726 bytes, checksum: a8015d8b04c2edbedf60057ec1e4ebba (MD5) Previous issue date: 2013 / Esta dissertação de mestrado originou-se da necessidade de avaliação de estratégias numéricas que permitam uma adequada predição in silico dos fenômenos de transferência de calor e massa em sistemas multicomponentes e multifásicos. Em particular, o foco do trabalho consistiu na implementação de modelos robustos em ferramenta de CFD. Para tanto, as equações de Maxwell-Stefan em conjunto com a lei de Fick generalizada foram codificadas na linguagem de programação C e o código gerado foi dinamicamente acoplado ao código comercial ANSYS® CFD (FLUENT®), versão 14.0. Determinou-se, então, a composição e a temperatura de equilíbrio em um sistema vapor-líquido, sendo cada fase composta por uma mistura de quatro hidrocarbonetos puros (metano, n-pentano, n-hexano e n-octano). Levou-se em consideração todas as características inerentes a misturas multicomponentes (a exemplo da correção para altas taxas e da contribuição de todos os gradientes de composição para a taxa de transferência da espécie i), o que introduz uma complexidade considerável ao sistema quando comparada à solução de escoamentos envolvendo misturas binárias. Os resultados obtidos nos estudos com a ferramenta de CFD foram verificados com predições obtidas do simulador de processos comercial em estado estacionário denominado PRO/II®, versão 8.2. Um estudo preliminar com a mistura binária água/ar foi conduzido com o objetivo de validar os resultados obtidos com o código CFD por meio de comparação direta com dados fornecidos pela carta psicrométrica padrão. <br> / Abstract : This dissertation was originated from the need of the evaluation of numerical strategies that may allow an adequate prediction in silico of heat and mass transfer phenomena in multicomponent and multiphase systems. In particular, the aim of the work has been to implement robust models in a CFD tool. Therefore, the Maxwell-Stefan's equations in conjunction with the generalized Fick's law have been codified in the C programming language and the code generated has been dynamically coupled to the commercial code ANSYS® CFD (FLUENT®), version 14.0. It has been determined, then, the equilibrium composition and temperature in a vapor-liquid system, in which each phase has been composed by a mixture of four pure hydrocarbons (methane, n-pentane, n-hexane and octane). All the inherent characteristics of multicomponent mixtures (e.g., the high flux correction and the contribution of all composition gradients to the transport rate of species i) have been taken into account, which introduces considerable complexity to the system when compared to the solution of flow involving binary mixtures. The results obtained in the studies conducted with the CFD code have been verified with predictions obtained from the steady state process simulator named PRO/II®, version 8.2. A preliminary study with the binary mixture water/air has been carried out with the aim to validate the results obtained with the CFD code by means of direct comparison with data provided by standard psychrometric chart.

Made available in DSpace on 2014-06-11T19:23:39Z (GMT). No. of bitstreams: 0 Previous issue date: 2006-05-26Bitstream added on 2014-06-13T18:51:01Z : No. of bitstreams: 1 dias_jp_me_ilha.pdf: 1702007 bytes, checksum: 8fc663dc91ce05f27bfdce5cb7c35609 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Este trabalho apresenta um estudo numérico do escoamento bifásico com formação de espuma de uma mistura composta por óleo lubrificante sintético Freol a10 com pequenas quantidades de refrigerante R134a, usualmente encontrada em compressores de sistemas de refrigeração e ar condicionado. A geometria considerada é um tubo longo horizontal de seção circular constante. Em razão da queda de pressão inicialmente causada pelas forças de atrito viscoso, a solubilidade do refrigerante no óleo diminui e o refrigerante evapora da mistura líquida formando bolhas de gás que escoam com a fase líquida. A formação de bolhas de refrigerante pode ser tão intensa a ponto de o escoamento assumir a forma de espuma quando a fração de vazio atinge valores superiores à aproximadamente 70%. Um ponto importante com relação ao escoamento de espuma é que o seu comportamento é muito diferente dos escoamentos bifásicos convencionais, apresentando quedas de pressão bem superiores. Os resultados gerados para as distribuições de pressão e temperatura ao longo do escoamento são confrontados com dados experimentais para a validação do modelo numérico. Será verificado que o modelo representa satisfatoriamente os resultados experimentais para várias condições do escoamento. / This work presents a numerical study of the two-phase flow with foam formation of a mixture composed by synthetic lubricant oil Freol a10 with small amounts of refrigerant R134a, usually encountered in compressors of refrigeration and air conditioning systems. A straight horizontal tube with circular cross section was chosen to accomplish the simulation. Due to the pressure drop initially caused by the viscous friction forces, the refrigerant solubility in the oil reduces and the refrigerant evaporates from the liquid mixture (outgassing) forming gas bubbles that flow with the liquid phase. In this type of flow, the bubble formation can be so large that foam is formed as the void fraction reaches values above 70%. An important point with reference to foam is that the foam flow behavior is much different from the conventional two-phase flows, presenting larger pressure drops. Results for pressure and temperature distributions along the flow are compared with experimental data in order to validate the numerical model. It will be seen that the model represents satisfactorily the experimental results for many flow conditions.

Orientador: José Luiz Gasche / Banca: André Luiz Seixlack / Banca: Adriano da Silva / Resumo: Este trabalho apresenta um estudo numérico do escoamento bifásico com formação de espuma de uma mistura composta por óleo lubrificante sintético Freol a10 com pequenas quantidades de refrigerante R134a, usualmente encontrada em compressores de sistemas de refrigeração e ar condicionado. A geometria considerada é um tubo longo horizontal de seção circular constante. Em razão da queda de pressão inicialmente causada pelas forças de atrito viscoso, a solubilidade do refrigerante no óleo diminui e o refrigerante evapora da mistura líquida formando bolhas de gás que escoam com a fase líquida. A formação de bolhas de refrigerante pode ser tão intensa a ponto de o escoamento assumir a forma de espuma quando a fração de vazio atinge valores superiores à aproximadamente 70%. Um ponto importante com relação ao escoamento de espuma é que o seu comportamento é muito diferente dos escoamentos bifásicos convencionais, apresentando quedas de pressão bem superiores. Os resultados gerados para as distribuições de pressão e temperatura ao longo do escoamento são confrontados com dados experimentais para a validação do modelo numérico. Será verificado que o modelo representa satisfatoriamente os resultados experimentais para várias condições do escoamento. / Abstract: This work presents a numerical study of the two-phase flow with foam formation of a mixture composed by synthetic lubricant oil Freol a10 with small amounts of refrigerant R134a, usually encountered in compressors of refrigeration and air conditioning systems. A straight horizontal tube with circular cross section was chosen to accomplish the simulation. Due to the pressure drop initially caused by the viscous friction forces, the refrigerant solubility in the oil reduces and the refrigerant evaporates from the liquid mixture (outgassing) forming gas bubbles that flow with the liquid phase. In this type of flow, the bubble formation can be so large that foam is formed as the void fraction reaches values above 70%. An important point with reference to foam is that the foam flow behavior is much different from the conventional two-phase flows, presenting larger pressure drops. Results for pressure and temperature distributions along the flow are compared with experimental data in order to validate the numerical model. It will be seen that the model represents satisfactorily the experimental results for many flow conditions. / Mestre

[pt] Escoamentos intermitentes bifásicos são caracterizados pela sucessão de golfadas de líquido separadas por bolhas de diferentes tamanhos, escoando com frequência variável. A previsão deste tipo de escoamento é muito importante, uma vez que o mesmo pode ser encontrado em diversas aplicações industriais. No presente trabalho, é realizado um estudo do regime intermitente em tubulações verticais, utilizando o Modelo de Dois Fluidos uni-dimensional, juntamente como o método de Volumes Finitos. O Modelo de Dois Fluidos uni-dimensional requer a inclusão de diversas relações de fechamento para as fases líquida e gasosa. Estas relações devem representar realisticamente o fenômeno físico de interesse, e as equações resultantes devem ser bem postas. No presente trabalho, duas relações de fechamento foram consideradas: salto de pressão devido à curvatura da interface e parâmetro de distribuição de fluxo de quantidade de movimento. A presença do salto de pressão não se mostrou significativa, enquanto que o parâmetro relacionado a distribuição da fase líquida ampliou a região em que o sistema de equações é bem posto. Para os casos analisados, o parâmetro de distribuição da fase líquida apresentou influência nas características do escoamento. As principais grandezas do regime estatisticamente permanente, como comprimento e velocidade de translação da golfada e da bolha, foram comparadas com dados experimentais disponíveis na literatura, apresentando boa concordância. / [en] Intermittent two phase flows are characterized by a succession of liquid slugs separated by bubbles of different sizes, with varying frequency. The prediction of this type of flow is very important, since it can be found in several industrial applications. In the present work, studies of intermittent flows in vertical pipes are performed, employing the Two Fluid Model along with the Finite Volume method. The one-dimensional Two Fluid Model requires the inclusion of several closure relations for both gas and liquid phases. These relations should realistically represent the physical phenomenon of interest, and the resulting equations should be well posed. In this work, two closure relations were considered: pressure jump due to the interface curvature and momentum flux distribution parameter. The inclusion of the pressure jump did not present a significant influence, whereas the momentum flux parameter of the liquid phase increased the region in which the system of equations is well posed. For the cases analyzed, the momentum flux parameter of the liquid phase affected the flow characteristics. The main statistical slug quantities, such as length and translational speed of the bubble and slug, were compared with experimental data available in the literature, showing good agreement.

A simulação numérica de escoamento tem se desenvolvido em diferentes aspectos para melhorar sua exatidão e se aproximar cada vez mais da representação de fenômenos reais. O desenvolvimento científico sobre escoamentos turbulentos multifásicos dependem da simulação numérica, onde o custo computacional destas ainda é um fator determinante, principalmente quando se visa resultados mais precisos. Uma das grandes dificuldades é que códigos numéricos de alta ordem, utilizados para simular escoamentos bifásicos, produzem erros de dispersão numérica presentes principalmente na interface, já que nesta posição se apresenta uma descontinuidade nas propriedades físicas dos fluidos envolvidos. O presente trabalho tem o objetivo de propor uma formulação de alta precisão para escoamentos bifásicos. Para isso, realiza-se modificações em um código já existente, denominado Incompact3d que possui a possibilidade de paralelização para até 100:000 núcleos computacionais, utilizando uma decomposição 2D. Entre estas modificações destacam-se a implementação de uma nova formulação do método Level Set, além de pequenas modificações nas equações de Navier-Stokes, como a adição do termo de gravidade, a segregação do termo de difusividade e a adição de um termo referente a força de tensão superficial. Ainda propõem-se uma nova consideração do termo de pressão, que é separado em pressão não hidrostática e pressão hidrostática. Diversos testes de verificação e validação foram realizados, apresentando a capacidade de se realizar simulações numéricas de alta ordem para escoamentos bifásicos. Com esta metodologia, foi identificar características de fenômenos turbulentos bifásicos. / Flow numerical simulation has been developed in di erent ways to improve its results accuracy and have a better approximation to real phenomenon. To scientific development turbulent two-phase flow depends of numerical simulation, where computational cost is an important factor, mainly when accurate results are required. High order two-phase flow codes issue is that they produce numerical dispersion errors that are amplified with the presence of interface, since in its position exist physical proprieties discontinuity. This research main objective is to propose a high precision formulation to two-phase flow. For that, we modified Incompact3d code that has a 2D decomposition making parallelized calculation up to 100:000 computational cores. Some modifications are the implementation of a new Level Set Method formulation; some minor Navier-Stokes equation modifications as the addition of gravity term; the segregation of di usive term; and the addition of surface tension forces. We also proposed a new pressure term consideration where we separate it in hydrostatic and non-hydrostatic pressure. Many verification and validation tests where developed presenting the ability to perform high order numerical simulations for two-phase flows. Also it was possible to identify multiphase turbulent structures in the phenomena simulated.

Orientador: André Luiz Seixlack / Resumo: Neste trabalho apresenta-se um modelo distribuído para análise de condensadores do tipo parede aquecida, hot-wall, usados em refrigeradores domésticos. Nesse modelo considera-se o escoamento do fluido refrigerante no interior do tubo do condensador e a transferência de calor entre a placa de revestimento externo do refrigerador e o ar ambiente. No modelo, inclui-se também a transferência de calor ao longo da fita de alumínio, usada para fixação do tubo na placa de revestimento externo e a transferência de calor através do isolamento térmico para o interior do compartimento refrigerado. O escoamento no interior do tubo é considerado unidimensional e dividido em três regiões: região monofásica de vapor superaquecido, região bifásica líquido-vapor e região monofásica de líquido sub-resfriado. O escoamento bifásico é analisado segundo o modelo homogêneo, ou seja, considerando condições de equilíbrio térmico e hidrodinâmico entre as fases. Considera-se também a queda de pressão do escoamento no interior do tubo. A fita de alumínio e a placa de revestimento externo são consideradas como aletas com extremidades adiabáticas. As equações do escoamento do fluido refrigerante: conservação da massa, quantidade de movimento e conservação de energia, são resolvidas por integração numérica e a equação de conservação de energia ao longo da fita de alumínio é resolvida pelo método de Volumes Finitos. O sistema de equações algébricas, decorrente do processo de discretização das equações difere... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: This work presents a distributed model for the analysis of hot-wall condensers used in domestic refrigerators. In this model, the flow of refrigerant inside the condenser tube and the heat transfer between the outer coating plate of the refrigerator and the ambient air are considered. Also included in this model the heat transfer along the aluminum tape used to hold the tube to the outer coating plate and the heat transfer through the thermal insulation into the refrigerated compartment. The flow inside the tube is taken as one-dimensional and divided into three regions: superheated vapor region, two-phase liquid-vapor region and sub-cooled liquid region. The two-phase flow is analyzed according to the homogeneous model, that is, considering conditions of thermal and hydrodynamic equilibrium between the phases. It is also considered the pressure drop of the flow inside the tube. The aluminum tape and the outer coating plate are considered as fins with adiabatic ends. The equations of the refrigerant flow: mass conservation, momentum and energy conservation, are solved by numerical integration and the energy conservation equation along the aluminum tape is solved by the Finite Volume method. The set of algebraic equations, resulting from the discretization process of the differential equations, is solved iteratively by successive iterations, point-to-point along the condenser tube and the aluminum tape. The results obtained are compared with experimental data available in the ... (Complete abstract click electronic access below) / Mestre

[pt] A produção de petróleo offshore tem requerido cada vez mais a utilização de equipamentos de bombeio no leito marinho para a produção de poços a longas distâncias das unidades de produção. A grande maioria destes equipamentos tem dificuldade de operar com a presença de gás, ainda mais em condições de golfadas, onde a fração de vazio na sucção pode variar significativamente. Isto limita enormemente a janela operacional dos equipamentos. Uma forma de minimizar os efeitos das golfadas e, assim, ampliar a janela de operação destes equipamentos é a instalação de misturadores de fluxo na sucção das bombas, que atuam de forma a garantir uma maior sobrevida ao equipamento. Portanto, o projeto e a previsão do comportamento deste tipo de equipamento são fundamentais para bom funcionamento dos sistemas de bombeamento submarino. O presente trabalho investiga numericamente, com o modelo de dois fluidos, a eficiência de um tipo de misturador em uniformizar as frações de vazio. Diferentes configurações geométricas são analisadas e a influência das vazões de admissão de ambas as fases é discutida. A metodologia é validada através da comparação da queda de pressão através do misturador com dados experimentais disponíveis, onde se observou boa concordância. O modelo numérico mostra ainda que uma reserva de líquido é mantida no interior do equipamento e, o gás é forçado a se dispersar na corrente de líquido para sair do misturador. O comportamento do misturador frente à presença de golfadas também é investigado. Ao final do trabalho é proposta uma nova geometria, que se mostra mais tolerante a intermitências no escoamento. / [en] The offshore oil production is increasingly using pumping equipment at the seabed for production of wells at long distances from production units. Most of these devices have difficulty to operate in the presence of gas, especially in slug flow conditions, where the pump suction s void fraction can vary significantly. This problem limits the equipment s range of operation. One way to minimize the effects of slugging and thus enlarge the operation range of these devices is to use flow mixers at the pump s suction, which helps to ensure the operation even in the presence of large gas bubbles. Therefore, the design and prediction of the behavior of this type of equipment are essential for proper operation of subsea pumping systems. At the present work, the efficiency of a type of flow mixer to provide uniform void fractions at its outlet is numerically investigated, with the two fluid model. Different geometric configurations are analyzed and the influence of intake flow of both phases is discussed. The methodology is validated through comparison of the pressure drop through the flow mixer with available experimental data, where good agreement was observed. The numerical model also shows that a buffer of liquid is maintained inside the equipment and the gas is forced to disperse in the liquid stream to exit the mixer. The mixer s behavior in the presence of slugs is also investigated. Finally, a new geometry is proposed, which seems to be more efficient to avoid flow s intermittencies.

A simulação numérica de escoamento tem se desenvolvido em diferentes aspectos para melhorar sua exatidão e se aproximar cada vez mais da representação de fenômenos reais. O desenvolvimento científico sobre escoamentos turbulentos multifásicos dependem da simulação numérica, onde o custo computacional destas ainda é um fator determinante, principalmente quando se visa resultados mais precisos. Uma das grandes dificuldades é que códigos numéricos de alta ordem, utilizados para simular escoamentos bifásicos, produzem erros de dispersão numérica presentes principalmente na interface, já que nesta posição se apresenta uma descontinuidade nas propriedades físicas dos fluidos envolvidos. O presente trabalho tem o objetivo de propor uma formulação de alta precisão para escoamentos bifásicos. Para isso, realiza-se modificações em um código já existente, denominado Incompact3d que possui a possibilidade de paralelização para até 100:000 núcleos computacionais, utilizando uma decomposição 2D. Entre estas modificações destacam-se a implementação de uma nova formulação do método Level Set, além de pequenas modificações nas equações de Navier-Stokes, como a adição do termo de gravidade, a segregação do termo de difusividade e a adição de um termo referente a força de tensão superficial. Ainda propõem-se uma nova consideração do termo de pressão, que é separado em pressão não hidrostática e pressão hidrostática. Diversos testes de verificação e validação foram realizados, apresentando a capacidade de se realizar simulações numéricas de alta ordem para escoamentos bifásicos. Com esta metodologia, foi identificar características de fenômenos turbulentos bifásicos. / Flow numerical simulation has been developed in di erent ways to improve its results accuracy and have a better approximation to real phenomenon. To scientific development turbulent two-phase flow depends of numerical simulation, where computational cost is an important factor, mainly when accurate results are required. High order two-phase flow codes issue is that they produce numerical dispersion errors that are amplified with the presence of interface, since in its position exist physical proprieties discontinuity. This research main objective is to propose a high precision formulation to two-phase flow. For that, we modified Incompact3d code that has a 2D decomposition making parallelized calculation up to 100:000 computational cores. Some modifications are the implementation of a new Level Set Method formulation; some minor Navier-Stokes equation modifications as the addition of gravity term; the segregation of di usive term; and the addition of surface tension forces. We also proposed a new pressure term consideration where we separate it in hydrostatic and non-hydrostatic pressure. Many verification and validation tests where developed presenting the ability to perform high order numerical simulations for two-phase flows. Also it was possible to identify multiphase turbulent structures in the phenomena simulated.