ABSTRACT

Inside each fluid, ρ and µ are constants. Equations (95) and (96) were solved by using a finite difference method on a fixed two-or three-dimensional grid. The spatial terms were discretized by second-order finite differences on a staggered Eulerian grid. The discretization of time was achieved by an explicit Euler method or a second order

Adams-Bashforth method. The boundary conditions used in their study were either periodic or full slip in the horizontal directions and rigid, stress-free on the top and bottom.