A Model for the Dynamics of Stable Gas Bubbles in Viscoelastic Fluids Based on Bubble Volume Variation

We present a novel formulation of the Rayleigh–Plesset equation to describe stable gas bubble dynamics in viscoelastic media, using bubble volume variation, rather than radius, as the primary variable of the resulting nonlinear ordinary differential equation. This formulation incorporates the linear Kelvin–Voigt model as the constitutive relation for the surrounding fluid, capturing both viscous and elastic contributions, to track the oscillations of a gas bubble subjected to an ultrasonic field over time. The proposed model is solved numerically, subjected to a convergence analysis, and validated by comparisons with theoretical and experimental results from the literature. We systematically investigate the nonlinear oscillations of a single spherical gas bubble in various viscoelastic environments, each modeled with varying levels of rheological complexity. The influence of medium properties, specifically shear elasticity and viscosity, is examined in detail across both linear and nonlinear regimes. This work improves our understanding of stable cavitation dynamics by emphasizing key differences from Newtonian fluid behavior, resonance frequency, phase shifts, and oscillation damping. Elasticity has a pronounced effect in low-viscosity media, whereas viscosity emerges as the dominant factor modulating the amplitude of oscillations in both the linear and nonlinear regimes. The model equation developed here provides a robust tool for analyzing how viscoelastic properties affect bubble dynamics, contributing to improved the prediction and control of stable cavitation phenomena in complex media.