Numerical prediction of interfacial instabilities: Sharp interface method (SIM)

We introduce a sharp interface method (SIM) for the direct numerical simulation of unstable fluid-fluid interfaces. The method is based on the level set approach and the structured adaptive mesh refinement technology, endowed with a corridor of irregular, cut-cell grids that resolve the interfacial region to third-order spatial accuracy. Key in that regard are avoidance of numerical mixing, and a least-squares interpolation method that is supported by irregular datasets distinctly on each side of the interface. Results on test problems show our method to be free of the spurious current problem of the continuous surface force method and to converge, on grid refinement, at near-theoretical rates. Simulations of unstable Rayleigh-Taylor and viscous Kelvin-Helmholtz flows are found to converge at near-theoretical rates to the exact results over a wide range of conditions. Further, we show predictions of neutral-stability maps of the viscous Kelvin-Helmholtz flows (Yih instability), as well as self-selection of the most unstable wave-number in multimode simulations of Rayleigh-Taylor instability. All these results were obtained with a simple seeding of random infinitesimal disturbances of interface-shape, as opposed to seeding by a complete eigenmode. For other than elementary flows the latter would normally not be available, and extremely difficult to obtain if at all. Sample comparisons with our code adapted to mimic typical diffuse interface treatments were not satisfactory for shear-dominated flows. On the other hand the sharp dynamics of our method would appear to be compatible and possibly advantageous to any interfacial flow algorithm in which the interface is represented as a discrete Heaviside function.