In this paper, it is assumed that the flow field isirrotational and incompressible at the underwater explosion pulsation phase. Thepotential theory is applied to deduce the governing equation of bubbleboundary motion. The underwater explosion bubble motion is studied byusing the boundary integral method with high order curve triangular-elements tomodel the three-dimensional bubble surface. The computational resultsare comparedwith Rayleigh-Plesset bubble model and experimental data, which showthat the high order surface element can precisely simulate the motion ofunderwater explosion bubble and is superior to linear elements.The gravity, non-gravity and rigid wall have an effect upon the bubblemotion. The paper has simulated these cases, and predicted the dynamicbehavior of underwater explosion bubble in a flow field, such as expansion,collapse, jetting, and also simulated kinetic characteristics ofthree-dimensional bubble around cylinder by establishing a three-dimensionalmodel, which reflects interactions between cylinder, free surface and theunderwater explosion bubble motion.