Fusion of bubbles is a common and significant phenomenon in the nature. In this work, the potential flow theory and the boundary integral method are employed to simulate the fusion of underwater explosion bubbles. Based on the numerical and experimental results in the literature, a three-dimensional model of bubble fusion is established. A good agreement is shown between the numerical results and experimental data. Influences
of characteristic parameters such as distance and depth are specifically investigated by using the developed three-dimensional program, resulting in favorable curves and conclusions. This work provides references for the relevant research on fusion mechanics and dynamic characteristics of bubbles.