[Paper Review] A comprehensive bifurcation method to analyze the super-harmonic and ultra-harmonic behavior of the acoustically excited bubble oscillator

AJ Sojahrood, Dana Wegierak+3 more|arXiv (Cornell University)|Sep 24, 2018

Ultrasound and Cavitation Phenomena39 references4 citations

TL;DR

This paper introduces a novel bifurcation analysis method that identifies super-harmonic (SuH) and ultra-harmonic (UH) oscillations in acoustically driven bubbles by tracking oscillation maxima, revealing hidden resonances missed by conventional methods. The approach enables precise mapping of SuH and UH behavior, particularly in lipid-shelled microbubbles, and provides a framework for optimizing ultrasound-based medical and sonochemical applications.

ABSTRACT

Acoustically excited bubbles are involved in a wide range of phenomena and applications ranging from oceanography to sonoluminescence; they have applications in chemistry, medical imaging, and therapeutic ultrasound. The complexity of bubble dynamics and the limited understanding of their behavior restricts the exploration of their full potential. The bubble oscillator is a highly nonlinear system, which makes it difficult to generate a comprehensive understanding of its oscillatory behavior. One method used to investigate such complex dynamical systems is the bifurcation analysis. Numerous investigations have employed the method of bifurcation diagrams to study the effect of different control parameters on the bubble behavior. These studies, however, focused mainly on investigating the subharmonic (SH) and chaotic oscillations of the bubbles. Super-harmonic (SuH) and ultra-harmonic (UH) bubble oscillations remain under-investigated. One reason is that the conventional method used for generating bifurcation diagrams cannot reliably identify features that are responsible for the identification of SuH and UH oscillations. Additionally, the conventional method cannot distinguish between the UHs and SHs. We introduce a simple procedure for the generation of bifurcation diagrams to address this shortcoming. This method selects the maxima of the bubble oscillatory response and plots them alongside the traditional bifurcation points for the corresponding control parameter. Through applying this method, the oscillatory behavior of the bubble oscillator is analyzed, and stable SuH and UH bubble oscillations are investigated. Based on this new analysis, the conditions for the generation and amplification of UH and SuH regimes are discussed.

Motivation & Objective

• To address the lack of comprehensive analysis methods for super-harmonic (SuH) and ultra-harmonic (UH) oscillations in acoustically excited bubble oscillators.
• To overcome the limitations of conventional bifurcation analysis, which samples data once per acoustic period and fails to reveal sub-harmonic and ultra-harmonic dynamics.
• To develop a method that extracts and plots the maxima of bubble oscillations to expose hidden SuH and UH resonances.
• To provide a framework for identifying parameter regimes where SuH and UH oscillations dominate, especially in nonlinear systems like lipid-shelled microbubbles.
• To enable optimization of ultrasound applications such as contrast-enhanced imaging, drug delivery, and histotripsy by revealing previously undetected nonlinear behaviors.

Proposed method

• The method modifies conventional bifurcation diagrams by plotting the local maxima of radial bubble oscillations (R/R₀) instead of sampling at fixed intervals per acoustic cycle.
• It uses the 4th-order Runge-Kutta method to numerically solve the Keller-Miksis equation for bubble dynamics under acoustic excitation.
• The backscattered pressure (PSc) is calculated from radial acceleration and velocity to analyze frequency content via spectral analysis.
• Phase portraits and time-domain R(t) curves are used alongside bifurcation diagrams to validate oscillation regimes and identify periodicity.
• The approach is applied alongside traditional bifurcation analysis to compare and contrast conventional sampling with maximum-tracking methods.
• Control parameters such as initial radius (R₀), driving frequency (f), and pressure amplitude (PA) are systematically varied to map nonlinear dynamics.

Experimental results

Research questions

• RQ1 Why do conventional bifurcation methods fail to reveal super-harmonic and ultra-harmonic oscillations in acoustically driven bubbles?
• RQ2 What conditions enable the emergence of 2nd-order super-harmonic (SuH) and ultra-harmonic (UH) resonances in bubble oscillators?
• RQ3 How do the maxima of radial oscillations correlate with the presence of specific harmonic components in the backscattered pressure spectrum?
• RQ4 Can the proposed method distinguish between period-1 (P1) and period-2 (P2) oscillation regimes when conventional methods obscure underlying sub-harmonic behavior?
• RQ5 How do initial conditions and nonlinear shell dynamics (e.g., lipid shells) influence the generation of SuH and UH components?

Key findings

• The maxima-based bifurcation method successfully reveals 2nd-order super-harmonic (SuH) resonance at f = 2.6 MHz and PA = 275 kPa, where the signal exhibits two distinct maxima per cycle.
• At f = 1.2 MHz and PA = 145 kPa, the method identifies a P2 oscillation with four maxima per cycle, indicating a 2nd-order SuH resonance with strong 5/2 and 7/2 ultra-harmonic components.
• The 2nd-order SuH component was found to be the strongest signal in the backscattered pressure spectrum, surpassing both 1/2-order sub-harmonics and other UH components.
• Conventional bifurcation diagrams misrepresent the dynamics by showing P1 or P2 regimes when the maxima method reveals multiple internal maxima, indicating complex SuH and UH behavior.
• Ultra-harmonic (UH) components such as 5/2 and 7/2 were stronger than 1/2-order sub-harmonics, indicating that UH resonance can dominate under specific parameter settings.
• The method enables identification of parameter ranges where UH and SuH oscillations coexist and dominate, particularly in lipid-shelled microbubbles, which exhibit enhanced nonlinear responses at low pressures.