A level set approach for computing solutions to incompressible two-phase flow is presented. The interface between the two fluids is considered to be sharp and is described as the zero level set of a smooth function. A new treatment of the level set method allows us to efficiently maintain the level set function as the signed distance from the interface. We never have to explicitly reconstruct or find the zero level set. Consequently, we are able to handle arbitrarily complex topologies, large density and viscosity ratios, and surface tension, on relatively coarse grids. We use a second order projection method along with a second order upwinded procedure for advecting the momentum and level set equations. We consider the motion of air bubbles and water drops. We also compute flows with multiple fluids such as air, oil, and water. |