Simulation of bubble dynamics in three-dimensional potential flows on heterogeneous computing systems using the fast multipole and boundary element methods

The bubble dynamics in potential flows of incompressible liquid is studied. The proposed approach is based on the boundary element method for the Laplace equation, which is especially efficient for the 3D bubble dynamics. In order to increase the problem size and to accelerate computations, an efficient numerical algorithm is developed and implemented. Depending on the problem size, for the acceleration of computations we used the direct matrix-vector multiplication on graphics processors (GPU) or the fast multipole method (FMM) implemented on heterogeneous computing systems (multicore CPUs and graphics processors). For the simulation of bubble surfaces, a new method based on the filtration of spherical harmonics is proposed. The proposed approach allows one to directly calculate the 3D dynamics of a single bubble, two and three interacting bubbles as well as a bubble cluster with a high degree of surface discretization. The developed method can be used to solve a wide range of problems related to the potential flow of bubble liquids.