ICONE19-43114 ON A GLOBAL MASS PRESERVING INCOMPRESSIBLE TWO-PHASE FLOW SOLVER

Abstract

An interface capturing scheme, implemented within the framework of the conservative phase-field method, is applied to simulate the incompressible two-phase flows with large density ratio and surface tension. For accurately predicting the solution, the finite difference dispersionrelation-preserving advection schemes are utilized to accommodate the true dispersion relation. The solutions computed from the resulting dispersion-relation-preserving advection schemes can minimize the dispersion error. The surface tension force is calculated by the continuum surface tension force (CSF) formulation of Brackbill et al.. (Brackbill, J. U., et al.., 1992) For the sake of programming simplicity, the proposed incompressible two-phase flow solver will be discretized on the semi-staggered grids without incurring velocity-pressure checkerboard oscillations. To verify the proposed method, four benchmark problems, including Rayleigh-Taylor instability, bubble rising and coaxial/oblique bubble merging problems, are numerically investigated. All predicted results have been shown to compare fairly well with the experimental and/or other numerical results. It also shows that the satisfactions of mass conservations are guaranteed for the proposed two-phase flow solver. Finally, the two-phase flow solver is applied to numerically predict the Kevin-Helmholtz instability in a tilted horizontal channel. The predicted results show the applicability to the stratified two-phase flow analysis in the nuclear engineering field, like the slug flow, droplet impacting or counter-current flow problems