A numerical assessment of shock-droplet interaction modeling including cavitation

A computational study of the wave phenomenon and probable cavitation that occurs when a droplet interacts with a shock is presented. A volume-of-fluid method with and without a cavitation model activated is employed. The model with no cavitation is used to clarify the droplet's internal wave field. The model is benchmarked using experimental data of a 22 mm water column impacted by a shock wave at Mach 2.4. Good agreement between the computational results for a two-dimensional (2D) circular droplet and the experimental results is demonstrated. A simplified 2D, square droplet is also considered in order to highlight the effects of interface curvature on the internal waves. Comparison of the results for the square and circular drops shows that the circular geometry diffuses the initial pressure wave while focusing the reflected wave. The three-dimensional (3D) spherical droplet was then analyzed. The internal wave behavior is similar to the 2D circular drop with enhanced diffusion and focusing. Two cavitation models were then used to explore the probability that cavitation occurs for a nominal sized raindrop interacting with shock waves, Mach 1.5-5. Two cavitation models, full Rayleigh-Plesset (FRP) and Schnerr-Sauer both predicted vapor formation at the same location with the FRP always predicting slightly higher amounts of vapor. The current setting for number of nuclei ensures that the simulation is not swamped by vapor formation but is unrealistically small. Therefore, while the work demonstrates that cavitation will occur, the full nature of the cavitation along with the impulsive pressure waves it should create is not currently captured.