Positional stability and radial dynamics of sonoluminescent bubbles under bi-harmonic driving: Effect of the high-frequency component and its relative phase

Numerical investigation of the translational motion of bubbles: The comparison of capabilities of the time-resolved and the time-averaged methods

In the present study, the accuracies of two different numerical approaches used to model the translational motion of acoustic cavitational bubble in a standing acoustic field are compared. The less accurate but less computational demanding approach is to decouple the equation of translational motion from the radial oscillation, and solve it by calculating the time-averaged forces exerted on the bubble for one acoustic cycle. The second approach is to solve the coupled ordinary differential equations directly, which provides more accurate results with higher computational effort. The investigations are carried out in the parameter space of the driving frequency, pressure amplitude and equilibrium radius. Results showed that both models are capable to reveal stable equilibrium positions; however, in the case of oscillatory solutions, the difference in terms of translational frequency may be more than three fold, and the amplitude of translational motion obtained by the time-averaged method is roughly 1.5 times higher compared to the time-resolved solution at particular sets of parameters. This observation implies that where the transient behaviour is important, the time-resolved approach is the proper choice for reliable results.

GPU accelerated numerical investigation of the spherical stability of an acoustic cavitation bubble excited by dual-frequency

The spherical stability of an acoustic cavitation bubble under dual-frequency excitation is investigated numerically. The radial dynamics is described by the Keller-Miksis equation, which is a second-order ordinary differential equation. The surface dynamics is modelled by a set of linear ordinary differential equation according to Hao and Prosperetti (1999), which takes into account the effect of vorticity by boundary layer approximation. Due to the large amount of investigated parameter combinations, the numerical computations were carried out on graphics processing units. The results showed that for bubble size between R_E=2μm and 4μm, the combination of a low and a high frequency, and the combination of two close but not equal frequencies are important to prevent the bubble losing its shape stability, while reaching the chemical threshold (R_max/R_E=3) (Kalmár et al., 2020). The phase shift between harmonic components of dual-frequency excitation has no effect on the shape stability.

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.

Relationship between the radial dynamics and the chemical production of a harmonically driven spherical bubble

The sonochemical activity and the radial dynamics of a harmonically excited spherical bubble are investigated numerically. A detailed model is employed capable to calculate the chemical production inside the bubble placed in water that is saturated with oxygen. Parameter studies are performed with the control parameters of the pressure amplitude, the forcing frequency and the bubble size. Three different definitions of collapse strengths (extracted from the radius vs.time curves) are examined and compared with the chemical output of various species. A mathematical formula is established to estimate the chemical output as a function of the collapse strength; thus, the chemical activity can be predicted without taking into account the chemical kinetics into the bubble model. The calculations are carried out by an in-house code exploiting the high processing power of professional graphics cards (GPUs). The results shown that chemical activity can be approximated qualitatively from the values of relative expansion. This could be helpful in order to optimise chemical output of sonochemical reactors either from measurement data or simulations as well.

Cross-frequency couplings in non-sinusoidal dynamics of interacting oscillators: Acoustic estimation of the radial position and spatial stability of nonlinear oscillating bubbles

In this work, the analysis of cross-frequency couplings (CFC) is introduced in the context of nonlinear acoustics related to the dynamics of bubble(s)-resonator systems. The results obtained from experiments specifically designed to untangle the causal connection between the CFC patterns observed at the signal level and the underlying physical processes, are discussed. It was found that "causal" amplitude-to-amplitude (AAC) and amplitude-to-phase (APC) couplings emerge in the system dynamics as a consequence of the bubble(s)-resonator mechanistic interaction in the oscillatory steady-state. In these CFC patterns, the amplitude of the fundamental frequency component (f₀) effectively modulates the amplitude and relative phase of the harmonic components (Nf₀). Moreover, these AAC and APC couplings give rise to "epiphenomenal" phase-to-amplitude (PAC) and phase-to-phase (PPC) couplings, in which the link between modulating and modulated parameters represents a correlation rather than a causal connection. It is shown that these CFC patterns can be exploited to determine the presence, spatial stability and radial position of nonlinear oscillating bubble(s) trapped within the acoustic chamber. Potential applications of the proposed techniques are also discussed. Substantial evidence is presented showing that CFC patterns emerging from quasi-periodic non-sinusoidal waveforms are informative on the interaction between underlying oscillators.