Conselho Nacional de Desenvolvimento Científico e Tecnológico / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Fundação de Amparo a Pesquisa do Estado de Minas Gerais / Petrobrás / Em estudos de escoamentos do tipo colunas de bolhas o uso de ferramentas computacionais desempenham um importante papel por permitir que, através de experimentos virtuais, a comunidade científica agregue mais conhecimento sobre as forças que atuam nesse tipo de escoamento. Considerando a confiabilidade dos resultados, torna-se possível o aperfeiçoamento de modelos constitutivos de forças fluidodinâmicas como de arrasto, sustentação e massa virtual. O código AMR3D trabalha com malhas de refinamento localizado bloco-estruturadas que se adaptam dinamicamente às características de um dado escoamento. Utiliza-se um método de projeção para o acoplamento pressão-velocidade e o método do Volume de Fluido (VOF) para a representação das fases e interfaces do escoamento bifásico. As equações do balanço da quantidade de movimento linear são discretizadas no tempo por um esquema semi-implícito de segunda ordem (IMEX). A presente tese contribuiu com o código com a implementação de um novo critério de refinamento que captura as estruturas turbilhonares do escoamento e com um algoritmo para identificação de bolhas para que seja possível a realização de operações sobre uma mesma bolha em meio à uma população. O código é verificado e validado com problemas da literatura. Resultados de simulações de bolhas ascendentes isoladas em meio quiescente são apresentados e comparados com o diagrama de Grace indicando boa concordância. Análises de forças de resistência fluidodinâmicas atuantes nas bolhas foram apresentadas e coeficientes de arrasto, sustentação e massa virtual calculados. Através de comparações com correlações da literatura, foi possível observar concordância dos resultados. Simulações da interação entre duas bolhas de volumes conhecidos são realizadas e comparações são realizadas com resultados da literatura. Em um dos casos, o algorítimo de identificação é verificado ao rastrear consistentemente os volumes e centroides das bolhas. Simulações com população de bolhas foram apresentadas e, através do algoritmo de identificação, estima-se a força de resistência fluidodinâmica total que age em uma bolha em meio à população. Conclui-se que, com base nos resultados, que o código AMR3D, com as implementações realizadas neste trabalho, possui potencial para a investigação de escoamentos borbulhantes. / In bubble dynamics studies computational tools plays an important role as it allows, through virtual experimentations, the increase of knowledge about forces that act in this kind of flow. Considering results reliability, it is common to observe constitutive force models improvement, as in drag, lift and virtual mass correlations. The AMR3D code, developed in the Fluid Mechanics Laboratory - MFlab, uses locally block-structured meshes that dynamically adapts itself to characteristics of a given flow. A projection method is used for the pressure-velocity coupling, and the Volume of Fluid - VOF method for the phases representation, including interface dynamics. The equations of linear momentum balance are discretized with a second order semi-implicit scheme (IMEX). The present thesis contributes to the code with the implementation of a new refinement criteria that captures the main structures of the flow, and an identification algorithm that allows users to perform operations on a single bubble between a swarm of bubbles. The present code is verified and validated against well known literature problems. Results of isolated ascending bubbles in a quiescent media are presented and compared with Grace’s diagram, showing good agreement. Fluid-dynamic forces are computed and analyzed, as well as fluid-dynamics forces coefficients for drag, lift and virtual mass, showing good agreement with literature. Simulations of a pair of bubbles interacting were presented and comparisons were realized against results of literature, showing good proximity. In one of the cases, the bubble identification algorithm was tested as it calculates bubbles known volumes and centroid positions consistently. Bubble swarm simulations were presented, and the identification algorithm were also applied to track a bubble as it ascend though the domain to expose the fluid-dynamics forces acting on it. With the presented results, one can conclude that the AMR3D code has a potential for investigative approaches in bubbly flows. / Tese (Doutorado)