The Interaction of the Shock Wave with the Bubble and the Effect of Computational Grid Size on the Problem Simulation with A Fully Coupled Pressure-Based Algorithm

When a shock wave propagates through a flow field that has nonlinear thermodynamic properties, different processes occur simultaneously. Wave compression, wave refraction, and vortex generation are examples of these processes that cause the waveform and thermodynamic properties of the fluid to change. The interaction of a shock wave with a cylindrical bubble is an example of a wave-bubble collision problem in which all of the above processes are observed. Due to the high computational cost of density-based algorithms in solving compressible interfacial flow problems such as shock wave interaction with the two-phase flow, using a fully coupled pressure-based algorithm is a good solution that will solve the problem with proper accuracy while reducing computation time. In this paper, using this algorithm, the interaction of the shock wave with the bubble is investigated; while validating the results, the effect of the computational grid size and the method of discretization of the governing equations are determined. It was observed that by increasing the number of computational grids according to the first-order upwind method, the simulation results become more accurate, and the numerical diffusion amount decreases. Also, by changing the discretization method to second-order upwind, the instabilities on the interface of the two phases increase due to spurious fluctuations, and the shape of the interface obtained from the numerical solution moves away from the experimental results.