A comparison of theoretical investigation and numerical simulation of a single bubble oscillation is carried out in this paper.The theoretical research is based on solving Rayleigh-Plesset (R-P) equation and Keller-Miksis (K-M) equation using Runge-Kutta method.The numerical method focuses on the two discrete methods in the volume of fluid (VOF) method,geometric reconstruction(GR) and modified high resolution interface capturing (MHRIC).The results show that the interface captured by MHRIC in the collapse stage is more stable,and the evolution of bubble radius agrees better with the theoretical solution.The R-P equation and K-M equation,which omit the effect of the energy equation,have limitations when the bubble collapses under ultra-high pressure.