Interaction of two three-dimensional explosion bubbles

Abstract

The interaction of two underwater explosion bubbles was mathematically analyzed in this paper. Based on the assumption of potential flow, high-order curved elements were used to discretize the boundary integral equation and solve it. Assuming that gas inside the bubble follows the isentropic rule, the Euler-Lagrange method was used to trace the evolution of the bubble, and when calculating the singular integral, the singularity of the double-layer singular integral was eliminated by reconstructing a principal-value integral of double-layer potential so that a more precise result could be obtained. Elastic mesh technique (EMT) was also used when tracing the evolution of the bubble interface, and numerical smoothing wasn’t needed. A comparison of calculations using this three-dimensional model with results of the Reyleigh-Plesset bubble model shows that the three-dimensional model and calculation method in this paper is practical. This three-dimensional model was applied to simulate the interaction of two bubbles under the action of gravity, and the dynamic characteristics of two bubbles near the surface was also analyzed. Bubbles influenced by surface effects and gravity present severe non-linearity. This paper provides a reference for research into the dynamics of multi-bubbles.

Keywords

two bubbles / boundary integral / EMT / three-dimensional model

References

[1] Rayleigh J. W. On the pressure developed in a liquid during the collapse of a spherical cavity[J]. Philos Mag, 1917, 34: 94-98

[2] GIBSON D C. Cavitation adjacent to plane boundaries[A]. Proc 3rd Conf Hydraulic Fluid Mech[C]. Syndey, 1968.

[3] Kling C. L., Hammitt F. G. A photographic study of spark-induced cavitation bubble collapse[J]. Trans ASME, 1972, 94: 825-833

[4] Zhang Y. L., Yeo K. S., B. Khoo C., Chong W. K. Three-dimensional computation of bubbles near a free surface[J]. Comput Phys, 1998, 146: 105-123

[5] Klaseboer E. K., Hung C., Wang C., Wang C. W., Khoo B. C., Boyce P., Debono S., Charlier H. Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient/rigid structure[J]. J Fluid Mech, 2005, 537: 387-413

[6] Wang C., Khoo B. C. An indirect boundary element method for three-dimensional explosion bubbles[J]. Journal of Comput. Phys, 2004, 194: 451-480

[7] Lu C. 3-D Numerical simulation of underwater explosion bubble[J]. Chinese Journal of Aeronautics, 1996, 9(1): 38-42

[8] Lu C. 3-D Numerical simulation of underwater explosion bubble[J]. Chinese Journal of Aeronautics, 1996, 17(1): 92-95

[9] Xiao Z., Reginald B. H. An improved model for bubble formation using the boundary-integral method[J]. Chemical Engineering Science, 2005, 60(1): 179-186

[10] Best J. P., Kucera A. A numerical investigation of non-spherical rebounding bubbles[J]. J Fluid Mech, 1992, 245: 137-154

[11] Wilkerson S. A. A boundary integral approach to three dimensional underwater explosion bubble dynamics[D]. 1990, Baltimore: The Johns Hopkins University

[12] Wang C., Khoo B. C., Yeo K. S. Elastic mesh technique for 3D BIM simulation with an application to underwater explosion bubbles [J]. Computers and Fluids, 2003, 32(9): 1195-1212

[13] Wang Q. X., Yeo K. S., Khoo B. C., Lam K. Y. Strong interaction between a buoyancy bubble and a free surface[J]. Theor Comput Fluid Dyna, 1996, 8: 73-88

[14] Wang Q. X., Yeo K. S., Khoo B. C., Lam K. Y. Nonlinear interaction between gas bubble and free surface[J]. Comput Fluids, 1996, 25(7): 607-628

[15] Rungsiyaphornrat S., Klaseboer E., Khoo B. C., Yeo K. S. The merging of two gaseous bubbles with an application to underwater explosions[J]. Computers & Fluids, 2003, 32: 1049-1074

[16] Cole R. H. Underwater explosion[M]. 1948, Princeton: Princeton University Press

[17] Brebbia C. A. The boundary element method for engineers [M]. 1978, New York: Halstead Press

[18] Yao X. L., Zhang A. Man.Simulation of the motion of three-dimensional underwater explosion bubble using simple Green function method[J]. Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(6): 749-759