A gas bubble in an acoustic field obeys the Keller-Miksis equation, a compressible version of Rayleigh-Plesset: (1 - V/c) R dV/dt + (3/2)(1 - V/(3c)) V^2 = (1 + V/c)(p_L - p_inf)/rho + (R/(rho c)) d(p_L - p_inf)/dt, where V = dR/dt, c = 1481 m/s, p_L = p_g - 4 mu V/R - 2 sigma/R and the acoustic pressure is p_inf = p0 - P_a sin(2 pi f t). The gas follows p_g = (p0 + 2 sigma/R0) [(R0^3 - h^3)/(R^3 - h^3)]^kappa with a van der Waals hard core h = R0/8.5 and kappa = 1.4.

For a weak drive the bubble is a mass-spring oscillator with the Minnaert resonance f0 = (1/(2 pi R0)) sqrt((3 kappa p_g0 - 2 sigma/R0)/rho), about 0.8 MHz for R0 = 4.5 um, far above the 26.5 kHz drive, so the response is quasi-static. At about 1.3 atm the bubble grows to roughly 9 R0 during the rarefaction and then collapses almost in free fall; the wall approaches the speed of sound and the gas is heated adiabatically, T = T0 (Rmax/Rmin)^(3(kappa-1)), which is where the light flash originates in real experiments.

Simplifications: no heat or mass transfer, no gas dissociation, no shape instabilities, no radiated shock; the temperature shown is an idealised upper-bound-style estimate, not a measurement, and the real flash physics is still debated. The simulation runs one bubble in a closed flask and is purely educational.