Numerical study of underwater shock wave by a modified method of characteristics

This paper introduces a modified method of characteristics as well as its application to simulation of a 1D spherical underwater explosion. To check the performance of the modified method, corresponding codes for computer calculation are developed to simulate the underwater explosion problem which is a typical isentropic flow problem. In applying the modified method, shock wave is calculated based on the Rankine-Hugoniot conservation relations. Artificial viscosity is not used in the simulation, and thus the corresponding influence of artificial viscosity is not introduced into the simulation. The work is mainly focused on underwater shock wave and secondary shock wave. The results simulated with the modified method are compared with other results from experiment and AUTODYN software, and the comparisons show that the modified method results are coincident with the experimental results in acceptable accuracy. Compared with the AUTODYN results, the modified method results are consistent with the experiment...