Simulation of three-dimensional bubbles using desingularized boundary integral method

Abstract

A very simple model based on the three-dimensional desingularized boundary integral method is applied to study the evolution of bubble(s) with or without the presence of solid structures. The choice of the desingularization parameters, which is crucial to the success of the method, is studied in the context of bubble dynamics. With the proper choice of parameters, the new model is far more efficient than previous models with virtually the same level of accuracy being achieved. This is largely attributed to the simplicity of the desingularization method. Furthermore, the new model offers a simple and attractive way for mesh refinement. Although it has limitations in the sense that, with this model the time stepping tends to slow down as two surfaces approach each other, this can be easily rectified by switching over to a direct method so that the two surfaces can be drawn closer as required in the context of jet impact. After this the new model can be reinstated to treat the complicated doubly connected geometry involving toroidal bubbles that would otherwise be very difficult to deal with.