Simulação numérica direta de alta ordem para escoamentos bifásicos água-ar

Abstract

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.