A time splitting projection scheme for compressible two-phase flows. Application to the interaction of bubbles with ultrasound waves

This paper is focused on the numerical simulation of the interaction of an ultrasound wave with a bubble. Our interest is to develop a fully compressible solver in the two phases and to account for surface tension effects.As the volume oscillation of the bubble occurs in a low Mach number regime, a specific care must be paid to the effectiveness of the numerical method which is chosen to solve the compressible Euler equations. Three different numerical solvers, an explicit HLLC (Harten–Lax–van Leer-Contact) solver [48], a preconditioning explicit HLLC solver [14] and the compressible projection method [21], [53] and [55], are described and assessed with a one dimensional spherical benchmark. From this preliminary test, we can conclude that the compressible projection method outclasses the other two, whether the spatial accuracy or the time step stability are considered.Multidimensional numerical simulations are next performed. As a basic implementation of the surface tension leads to strong spurious currents and numerical instabilities, a specific velocity/pressure time splitting is proposed to overcome this issue. Numerical evidences of the efficiency of this new numerical scheme are provided, since both the accuracy and the stability of the overall algorithm are enhanced if this new time splitting is used. Finally, the numerical simulation of the interaction of a moving and deformable bubble with a plane wave is presented in order to bring out the ability of the new method in a more complex situation.