A Model for the Dynamics of Stable Gas Bubbles in Viscoelastic Fluids Written in 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 and validated by comparison with theoretical 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 prediction and control of stable cavitation phenomena in complex media.