Characterization of Acoustic Streaming Beyond 100 MHz

High-speed and acceleration micrometric jets induced by GHz streaming: A numerical study with direct numerical simulations

Gigahertz acoustic streaming enables the synthesis of localized microjets reaching speeds of up to meters per second, offering tremendous potential for precision micromanipulation. However, theoretical and numerical investigations of acoustic streaming at these frequencies remain so far relatively scarce due to significant challenges including: (i) the inappropriateness of classical approaches, rooted in asymptotic development, for addressing high-speed streaming with flow velocities comparable to the acoustic velocity; and (ii) the numerical cost of direct numerical simulations generally considered as prohibitive. In this paper, we investigate high-frequency bulk streaming using high-order finite difference direct numerical simulations. First, we demonstrate that high-speed micrometric jets of several meters per second can only be obtained at high frequencies, due to diffraction limits. Second, we establish that the maximum jet streaming speed at a given actuation power scales with the frequency to the power of 3/2 in the low attenuation limit and linearly with the frequency for strongly attenuated waves. Last, our analysis of transient regimes reveals a dramatic reduction in the time required to reach the maximum velocity as the frequency increases (power law in -5/2), leading to characteristic time on the order of μs at gigahertz frequencies, and hence accelerations within the Mega-g range.

GPU accelerated study of a dual-frequency driven single bubble in a 6-dimensional parameter space: The active cavitation threshold

The active cavitation threshold of a dual-frequency driven single spherical gas bubble is studied numerically. This threshold is defined as the minimum intensity required to generate a given relative expansion (R_max-R_E)/R_E, where R_E is the equilibrium size of the bubble and R_max is the maximum bubble radius during its oscillation. The model employed is the Keller-Miksis equation that is a second order ordinary differential equation. The parameter space investigated is composed by the pressure amplitudes, excitation frequencies, phase shift between the two harmonic components and by the equilibrium bubble radius (bubble size). Due to the large 6-dimensional parameter space, the number of the parameter combinations investigated is approximately two billion. Therefore, the high performance of graphics processing units is exploited; our in-house code is written in C++ and CUDA C software environments. The results show that for (R_max-R_E)/R_E=2, the best choice of the frequency pairs depends on the bubble size. For small bubbles, below 3μm, the best option is to use just a single frequency of a low value in the giant response region. For medium sized bubbles, between 3μm and 6μm, the optimal choice is the mixture of low frequency (giant response) and main resonance frequency. For large bubbles, above 6μm, the main resonance dominates the active cavitation threshold. Increasing the prescribed relative expansion value to (R_max-R_E)/R_E=3, the optimal choice is always single frequency driving with the lowest value (20kHz here). Thus, in this case, the giant response always dominates the active cavitation threshold. The phase shift between the harmonic components of the dual-frequency driving (different frequency values) has no effect on the threshold.