An inverse method to fast-track the calculation of phase diagrams for sonoluminescing bubbles

A sound driven air bubble can be transformed into an argon bubble emitting light pulses stably. The very foundation to investigate the sonoluminescing bubble is to accurately determine the ambient radius and gas composition in the interior. The conventional approach is to model the air-to-argon transformation process through a large number of bubble dynamics simulations to obtain the physical parameters of the ultimate argon bubble. In this paper, we propose a highly efficient method to pinpoint this information in a phase diagram. The method is based on the diffusive equilibrium for each species inside the bubble and derives the ambient radius and composition inversely. To calculate the former parameter, the bisection algorithm is employed to consecutively narrow down the searching range until the equilibria is approached. Afterward, several cycles of full dynamics simulations are conducted to refine the composition. The method is validated using published experimental data. The calculated ambient radii deviate from the test results by less than 1 μm, which falls within the margin of measurement error. The advantages of this method over the semi-analytical approach reported by Hilgenfeldt et al. [J. Fluid Mech. 365 (1998)] are also discussed. Our study provides a standard procedure to calculate the ambient radius and composition and is beneficial for the numerical simulation of sonoluminescing bubbles.