This study numerically investigates the underwater explosion bubble dynamics during the process from detonation to bubble jet with a hybrid algorithm based on Runge Kutta Discontinuous Galerkin-Level Set-Direct Ghost Fluid method (RKDG-LS-DGF) and boundary element method (BEM). RKDG-LS-DGF method is used to simulate the physical process from the detonation of a column charge at one end to the formation of a nonspherical initial bubble and the process of bubble jet. And BEM is adopted to simulate bubble pulsing characteristics. In addition, the numerical results are compared with the experimental results to verify the feasibility of the numerical method. It is found that, during the detonation process of a column charge, the detonation product experiences a shape change from an initial ellipsoid into a sphere during expansion. After the detonation, the bubble experiences expansion and contraction and develops a jet. The jet threads through the bubble in the opposite direction to gravity and induces a high-pressure region. Subsequently, the pressure of this region decreases when the bubble reexpands after being penetrated by the jet. The numerical results agree well with the experimental data, which proves that axisymmetric RKDG-LS-DGF method and BEM are successfully combined to simulate the whole process of underwater explosion.