Development of a dispersively accurate conservative level set scheme for capturing interface in two-phase flows

• Authors: Tony W. H. Sheu;C. H. Yu;P. H. Chiu
• Venue: Journal of Computational Physics
• Year: 2009

A two-step interface capturing scheme, implemented within the framework of conservative level set method, is developed in this study to simulate the gas/water two-phase fluid flow. In addition to solving the pure advection equation, which is used to advect the level set function for tracking interface, both nonlinear and stabilized features are taken into account for the level set function so that a sharply varying interface can be stably predicted. To preserve the conservative property, the mapping given by @x@?=x@?-u@?t is performed between two coordinates x@? and @x@? for the advection-diffusion equation in a flow field with velocity u. To capture the interface, the flux term capable of compressing the level set contours is also adopted in the construction of linear inviscid Burgers' equation, which is indispensable in the level set method. To resolve the physically sharp interface without incurring contact discontinuity oscillations, a damping term which is nonlinear in terms of the level set function and can render an adequate artificial diffusion to stabilize the contact surface is added into the reinitialization step in the modified level set method. For accurately predicting the level set function, the advection scheme for solving the linear inviscid Burgers' equation in the advection step of the modified level set method is developed to accommodate the true dispersion relation. The solution computed from the resulting two-dimensional dispersion-relation-preserving advection scheme can minimize the phase error. Less artificial damping is needed to damp the oscillations in the vicinity of contact surface and the interface can be less numerically smeared. For the sake of programming simplicity, the incompressible two-phase flow will be discretized in non-staggered grids without incurring checkerboard oscillations by the developed explicit compact scheme for the approximation of pressure gradient terms. For the verification of the proposed two-dimensional dispersion-relation-preserving scheme and the non-staggered incompressible flow solver, three benchmark problems have been chosen in this study. The proposed conservative level set method for capturing the interface in incompressible fluid flows is also verified by solving the dam-break, bubble rising in water, droplet falling in water and Rayleigh-Taylor instability problems. Good agreements with the referenced solutions are demonstrated for all the four investigated problems.