A Numerical Study on Stratified Shear Layers With Relevance to Oil Boom Failure

Abstract

Interface dynamics of twophase flow, with relevance for leakage of oil retained by mechanical oil barriers, is studied by means of a twodimensional (2D) latticeBoltzmann method (LBM) combined with a phasefield model for interface capturing. A multirelaxationtime (MRT) model of the collision process is used to obtain a numerically stable model at high Reynolds number flow. In the phasefield model, the interface is given a finite but small thickness, where the fluid properties vary continuously across a thin interface layer. Surface tension is modeled as a volume force in the transition layer. The numerical model is implemented for simulations with the graphic processing unit (GPU) of a desktop personal computer. Verification tests of the model are presented. The model is then applied to simulate gravity currents (GCs) obtained from a lockexchange configuration, using fluid parameters relevant for those of oil and water. Interface instability phenomena are observed, and obtained numerical results are in good agreement with theory. This work demonstrates that the numerical model presented can be used as a numerical tool for studies of stratified shear flows with relevance to oilboom failure.