Detailed Two-Dimensional Numerical Analysis of Multiphase Flows

Conselho Nacional de Desenvolvimento Científico e Tecnológico The mathematical modeling of multiphase flows involves the interaction between deformable and moving geometries with the fluid in which they are dispersed (immersed). This kind of interaction is present in many practical applications. A common approach to handle these problems is the so called Front-Tracking/Front-Capturing Hybrid Methods. This methodology consists in separating the problem into two domains: an eulerian domain, which is kept fixed and is used to discretize the fluid equations of both phases, and a lagrangian domain, which is used to solve the equations of motion of the interface. Since there is no geometric dependence between these two domains, the method can easily handle moving and deformable interfaces that are dispersed in the flow. Following this line of research, the present work aims to capture accurately details of such flows by resolving adequately the relevant physical scales in time and in space. This can be achieved by applying locally refined meshes which adapt dynamically to cover special flow regions, e.g. the vicinity of the fluid-fluid interfaces. To obtain the required resolution in time, a semi-implicit second order discretization to solve the Navier-Stokes equations is used. The turbulence modeling is introduced in the present work through Large Eddy Simulation. The eficiency and robustness of the methodology applied are verified via convergence analysis, as well as with simulations of one-phase and two-phase flows for several Reynolds numbers. The results of two-phase flows, with one bubble and with multiple bubbles, are presented. The results obtained for a single bubble case are compared with Clift's shape diagram (Clift et al., 1978). A modelagem matemática de escoamentos multifásicos envolve a interação de geometrias móveis e deformáveis com o meio fluido que as envolve. Este tipo de interação faz parte de uma extensa lista de aplicações. Uma linha proposta para o tratamento num érico deste tipo de problema são os métodos híbridos Front-Tracking/Front-Capturing. Esta abordagem leva à separação do problema em dois domínios distintos (líquido/gás e líquido/líquido), um domínio fixo, euleriano, utilizado para discretizar as equações de ambas as fases, e outro móvel, lagrangiano, usado para as interfaces. Para o presente trabalho, na metologia utilizada, ambos os domínios são geometricamente independentes e não apresentam restrição quanto ao movimento e à deformação da fase dispersa. Seguindo esta linha, no presente trabalho propõe-se capturar detalhes deste de tipo escoamento, resolvendo adequadamente as escalas físicas do tempo e do espaço, utilizando malhas bloco estruturada refinadas localmente, as quais se adaptam dinamicamente para recobrir as regiões de interesse do escoamento (como, por exemplo, ao redor da interface fluido-fluido). Para se obter a resolução necessária no tempo, é usada uma discretização semi-implícita de segunda ordem para solucionar as equações de Navier-Stokes. A modelagem da turbulência é introduzida no presente trabalho via Simulação de Grandes Escalas. A eficiência e a robustez da metodologia implementada são verificadas via análise de convergência do método, bem como a simulação de escoamentos monofásicos e bifásicos para diferentes números Reynolds. São também apresentados resultados para escoamentos bifásicos com uma só bolha assim como para múltiplas bolhas. Os resultados de escoamentos mono-bolhas são comparados com o diagrama de forma de Clift et al. (1978). Doutor em Engenharia Mecânica