An Eulerian scheme for direct numerical simulation of multibubble dynamicsin an acoustic field

An Eulerian numerical scheme for the direct simulation of bubble dynamics i n an acoustic field is proposed. The compressible Navier-Stokes equation wi th a surface tension term is used as a governing equation. The convection t erms in the equation are solved using the hybrid interpolation-extrapolatio n scheme which stably provides a nonsmoothed solution of density interface between bubbles and a surrounding material. The acoustic terms in the equat ion are solved by the generalized Crank-Nicholson method which is an implic it method and has a temporally second-order accuracy under the low-Courant- number condition and is applicable to high-Courant-number computation. Usin g this scheme, some typical results of single- and multibubble phenomena ar e given.