Modeling stable cavitation of coated microbubbles: A framework integrating smoothed dissipative particle dynamics and the Rayleigh-Plesset equation

Coated microbubbles are widely used in medical applications, particularly in enhanced drug and gene delivery. One of the mechanisms underlying these applications involves the shear stress exerted on the cell membrane by acoustic microstreaming generated through cavitation bubbles. In this study, we develop a novel simulation approach that combines the smooth dissipative particle dynamics (SDPD) simulation method with numerical modeling of the Rayleigh-Plesset-like equation in an ad hoc manner to simulate stable cavitation of microbubbles at microsecond and micrometer scales. Specifically, the SDPD method is utilized to model fluid dynamics, while the Rayleigh-Plesset-like equation is employed to describe bubble dynamics. Adopting a 1.5 mu m coated microbubble driven by ultrasound with a frequency of 2 MHz and a pressure of 500 kPa as a representative example, we observe a high-velocity microstreaming pattern emerging around the bubble on a very small scale of a few micrometers after only a few microseconds. These spatiotemporal scales may pose challenges for experimental observation. The formation of this microstreaming arises from the opposing motion of the fluid layer next to the bubble and the fluid layers further away. Furthermore, our simulations reveal high shear stress levels of thousands of Pascals exerted on a wall located a few micrometers from the bubble. This contrasts with the shear stress values of a few Pascals calculated from theoretical models in the literature, which do not incorporate radial streaming into their theories. The implications of our results for bubble cavitation-induced pore formation on the cell membrane are discussed in some details.