Calibration of Numerical Simulation Methods for Underwater Explosion with Centrifugal Tests

The centrifugal underwater explosion tests and corresponding numerical simulations were carried out to study the laws of shock wave and bubble pulsation. A semiempirical method to determine JWL state equation parameters was given. The influence on numerical results caused by factors such as state equation of water, boundary condition, and mesh size was analyzed by comparing with the centrifugal underwater explosion test results. The results show that the similarity criterion is also suitable in numerical simulation; the shock wave peak pressure calculated by polynomial state equation is smaller than that of shock state equation. However, the maximum bubble radius and the pulsation period calculated by the two equations are almost the same. The maximum bubble radius is mainly affected by the boundary simulating the test model box, and the pulsation period is mainly affected by the artificial cutoff boundary. With the increase of standoff distance of measuring point, the mesh size required for the numerical calculation decreases; the size of the two-dimensional model is recommended to take 1/30 approximate to 1/10 explosion radius. In three-dimensional models, when mesh size is 2 times larger than explosion radius, the bubble motion change in the second pulsation period is not obvious. When mesh size is smaller than 1 time explosive radius, the full period of bubble pulsation can be well simulated, but calculation errors increase slowly and computation time greatly increases, so the three-dimensional mesh size is suggested to take the charge radius. Shock wave peak pressure is basically unaffected by gravity. As the gravity increases, the bubble maximum radius and the first pulsation period both decrease.