Numerical simulations of stable cavitation bubble generation and primary Bjerknes forces in a three-dimensional nonlinear phased array focused ultrasound field

We present a model developed for studying the generation of stable cavitation bubbles and their motion in a three-dimensional volume of liquid with axial symmetry under the effect of finite-amplitude phased array focused ultrasound. The density of bubbles per unit volume is determined by a nonlinear law which is a threshold-dependent function of the negative acoustic pressure reached in the liquid, in which nuclei are initially distributed. The nonlinear mutual interaction of ultrasound and bubble oscillations is modeled by a nonlinear coupled differential system formed by the wave and a Rayleigh-Plesset equations, for which both the pressure and the bubble oscillation variables are unknown. The system, which accounts for nonlinearity, dispersion, and attenuation due to the bubbles, is solved by numerical approximations. The nonlinear acoustic pressure field is then used to evaluate the primary Bjerknes force field and to predict the subsequent motion of bubbles. In order to illustrate the procedure, a medium-high and a low ultrasonic frequency configurations are assumed. Simulation results show where bubbles are generated, the nonlinear effects they have on ultrasound, and where they are relocated. Despite many physical restrictions and thanks to its particularities (two nonlinear coupled fields, bubble generation, bubble motion), the numerical model used in this work gives results that show qualitative coherence with data observed experimentally in the framework of stable cavitation and suggest their usefulness in some application contexts.