Cavitation dynamics of nuclei are largely governed by the Rayleigh-Plesset equation. This research focuses on improving the efficiency of numerical solutions to the Rayleigh-Plesset equation using “singularity handling,” which significantly increases the stability of solution methods while retaining solution accuracy. The singularity handling approach is based detection of collapse events that cause solution instabilities, then recovering the solution with a reflected condition that conserves momentum across a collapse event. In this paper, various algorithms are explored and assessed within the context of schemes based on a fixed-time-step-size. The results indicate that solution approaches are achievable that recover unmodified Rayleigh-Plesset solution in the fine-time-step limit, while having the ability to run 980% faster with only a 7% error at much larger time steps. This increase in efficiency and accuracy allows the program to provide useful solutions in the field of fluids engineering, particularly, in the study underwater explosions, optimization methods, cavitation, and   by enabling cavitation simulations that incorporate Rayleigh-Plesset physics.

cavitation dynamics, Rayleigh-Plesset, numerical methods.