Three-dimensional adaptive, Cartesian grid method for multiphase flow computations

Multiphase flows are marked by the presence of interfacial dynamics, steep jump in fluid-properties across the interface and moving boundaries between different phases and materials. Presence of moving boundaries and nonlinear coupling of phases across interface makes three-dimensional numerical simulations a difficult and expensive task. The presented work uses a continuous interface combining marker points on a stationary Eulerian grid. For an efficient computation an adaptive Cartesian grid is used to resolve local length scales of the flow. The interface location and flow solution are used as the adaptation criteria. The multiphase-interface is tracked using unstructured surface grid representation. A volume preserving interface reorganization procedure has been developed to maintain the desired interface surface grid resolution. Also, a momentum-source computation technique improving the performance of immersed boundary treatment in the governing field equation has been investigated.