Shape stability of a gas bubble in a soft solid

Predicting the onset of non-spherical oscillations of bubbles in soft matter is a fundamental cavitation problem with implications to sonoprocessing, polymeric materials synthesis, and biomedical ultrasound applications. The shape stability of a bubble in a Kelvin-Voigt viscoelastic medium with nonlinear elasticity, the simplest constitutive model for soft solids, is analytically investigated and compared to experiments. Using perturbation methods, we develop a model reducing the equations of motion to two sets of evolution equations: a Rayleigh-Plesset-type equation for the mean (volume-equivalent) bubble radius and an equation for the non-spherical mode amplitudes. Parametric instability is predicted by examining the natural frequency and the Mathieu equation for the non-spherical modes, which are obtained from our model. Our theoretical results show good agreement with published experiments of the shape oscillations of a bubble in a gelatin gel. We further examine the impact of viscoelasticity on the time evolution of non-spherical mode amplitudes. In particular, we find that viscosity increases the damping rate, thus suppressing the shape instability, while shear modulus increases the natural frequency, which changes the unstable mode. We also explain the contributions of rotational and irrotational fields to the viscoelastic stresses in the surroundings and at the bubble surface, as these contributions affect the damping rate and the unstable mode. Our analysis on the role of viscoelasticity is potentially useful to measure viscoelastic properties of soft materials by experimentally observing the shape oscillations of a bubble.