Nonlinear ultrasonic waves in bubbly liquids with nonhomogeneous bubble distribution: Numerical experiments

This paper deals with the nonlinear propagation of ultrasonic waves in mixtures of air bubbles in water, but for which the bubble distribution is nonhomogeneous. The problem is modelled by means of a set of differential equations which describes the coupling of the acoustic field and bubbles vibration, and solved in the time domain via the use and adaptation of the SNOW-BL code. The attenuation and nonlinear effects are assumed to be due to the bubbles exclusively. The nonhomogeneity of the bubble distribution is introduced by the presence of bubble layers (or clouds) which can act as acoustic screens, and alters the behaviour of the ultrasonic waves. The effect of the spatial distribution of bubbles on the nonlinearity of the acoustic field is analyzed. Depending on the bubble density, dimension, shape, and position of the layers, its effects on the acoustic field change. Effects such as shielding and resonance of the bubbly layers are especially studied. The numerical experiments are carried out in two configurations: linear and nonlinear, i.e. for low and high excitation pressure amplitude, respectively, and the features of the phenomenon are compared. The parameters of the medium are chosen such as to reproduce air bubbly water involved in the stable cavitation process.     

Numerical simulation of gas bubble motion in a flow-through resonator

The problem of gas bubble motion in an acoustic resonator with a fluid flow is solved using numerical methods. It is shown that the distribution of the bubble concentration, which is nonuniform over the resonator length, is formed upon homogeneous introduction of bubbles. The problem on the bubble concentration distribution the along the resonator axis (with fluctuations of the bubble introduction period taken into account) is considered, and the fluctuation parameters are determined at which the periodic structure of the concentration distribution is preserved. The distribution of bubbles with different sizes over the resonator length is determined. It is shown that a resonator with a fluid flow accomplishes bubble selection by size (the average bubble concentration in the resonator increases with an increase in bubble size). The field in the resonator was calculated taking into account the effect of bubbles on sound velocity and damping.     

Spatial-temporal dynamics of cavitation bubble clouds in 1.2 MHz focused ultrasound field

Cavitation bubbles have been recognized as being essential to many applications of ultrasound. Temporal evolution and spatial distribution of cavitation bubble clouds induced by a focused ultrasound transducer of 1.2 MHz center frequency are investigated by high-speed photography. It is revealed that at a total acoustic power of 72 W the cavitation bubble cloud first emerges in the focal region where cavitation bubbles are observed to generate, grow, merge and collapse during the initial 600 μs. The bubble cloud then grows upward to the post-focal region, and finally becomes visible in the pre-focal region. The structure of the final bubble cloud is characterized by regional distribution of cavitation bubbles in the ultrasound field. The cavitation bubble cloud structure remains stable when the acoustic power is increased from 25 W to 107 W, but it changes to a more violent form when the acoustic power is further increased to 175 W.     

调频超声脉冲驱动微气泡运动偏移的研究

提出调频超声辐射力技术驱动微泡群,以加强微泡的吸附效率.基于改进的RP方程及粒子轨迹方程研究了微泡群整体的运动位移与调频信号的中心频率、调频范围、信号声压,以及微泡半径分布关系.研究结果表明调频信号在驱动半径具有宽泛分布的气泡群,以及半径分布远离谐振半径的气泡群时,作用效果好于传统正弦波信号.例如中心频率1 MHz、调频范围0.75 MHz的调频脉冲作用高斯分布（平均半径3.5 μm、均方差为1）的微泡群200 μs,可比同等声压的正弦波多约12%的微气泡产生位移30 μm. 关键词： 超声辐射力 调频波 高斯分布     

气泡线性振动时近海面气泡群的声散射

海洋中的不同成因的气泡群是常见的水下声学目标及声呐混响源,因此对水下气泡群进行声学建模意义重大。利用有效媒质理论描述气泡群内部的相速度及声衰减变化,并考虑到海洋中气泡群往往产生于不同界面附近,进一步利用球面波叠加原理描述海面对气泡群散射声波的再辐射,导出了平海面作用下气泡群声散射截面的一般表达式,建立了其声散射模型,研究了单一尺寸及混合尺寸气泡群的声学特性。数值分析表明,气泡群的谐振频率会随其半径或孔隙率增加而降低;由于海面的存在,气泡群声散射截面会随频率进行周期性变化,且随气泡群远离海面,这一变化逐渐加剧。此外,若气泡的黏滞阻尼项在全部阻尼项中占比较高,气泡群声散射强度会在谐振频率附近存在起伏振荡。该模型可为近海面鱼群、气泡羽流及海底泄漏的甲烷气体的声学建模提供一定的理论基础。     

Sound field measurement in a double layer cavitation cluster by rugged miniature needle hydrophones

During multi-bubble cavitation the bubbles tend to organize themselves into clusters and thus the understanding of properties and dynamics of clustering is essential for controlling technical applications of cavitation. Sound field measurements are a potential technique to provide valuable experimental information about the status of cavitation clouds. Using purpose-made, rugged, wide band, and small-sized needle hydrophones, sound field measurements in bubble clusters were performed and time-dependent sound pressure waveforms were acquired and analyzed in the frequency domain up to 20 MHz. The cavitation clusters were synchronously observed by an electron multiplying charge-coupled device (EMCCD) camera and the relation between the sound field measurements and cluster behaviour was investigated. Depending on the driving power, three ranges could be identified and characteristic properties were assigned. At low power settings no transient and no or very low stable cavitation activity can be observed. The medium range is characterized by strong pressure peaks and various bubble cluster forms. At high power a stable double layer was observed which grew with further increasing power and became quite dynamic. The sound field was irregular and the fundamental at driving frequency decreased. Between the bubble clouds completely different sound field properties were found in comparison to those in the cloud where the cavitation activity is high. In between the sound field pressure amplitude was quite small and no collapses were detected.     

鼓泡床中超声驻波的模拟及其对气泡的调制机理

采用计算流体动力学（CFD）的方法，数值生成了鼓泡床中一对声换能器以16kHz高频振动引发的超声场。数值计算是基于包括粘性影响的可压缩流体基本守恒方程，并耦合了水的状态方程。模拟结果表明，在本研究所用的几何布置和换能器与时间相关的速度入口边界条件下，反应器中形成了一个稳定的驻波声场；由于波的非线性以及水的粘性，压力波节点呈现出轻微的时间漂移性。模拟结果与前人的实验结果定性吻合。在模拟的声压分布的基础上，分析了驻波声场调制气泡的机理。如比较熟知，气泡在驻波声场作用下或者向压力波节点运动或者向压力波腹点运动，取决于气泡尺寸与共振尺寸的关系。     

Thresholds for cavitation produced in water by pulsed ultrasound

The threshold for transient cavitation produced in water by pulsed ultrasound was measured as a function of pulse duration and pulse repetition frequency at both 0.98 and 2.30 MHz. The cavitation events were detected with a passive acoustic technique which relies upon the scattering of the irradiation field by the bubble clouds associated with the events. The results indicate that the threshold is independent of pulse duration and acoustic frequency for pulses longer than approximately 10 acoustic cycles. The threshold increases for shorter pulses. The cavitation events are likely to be associated with bubble clouds rather than single bubbles.     

Interaction between moving bubbles and an acoustic field: Controlling flows and nonlinear acoustic vision

Results from studying the interaction between gas bubbles and the field of a flow-through acoustic resonator, and the Raman scattering of acoustic waves by moving bubbles, are presented. The structure of the distribution of bubble concentration in the resonator is studied. It is shown that nonlinear scattering by moving bubbles can be used to image bubble objects.     

The interaction of shockwaves with a vapour bubble in boiling histotripsy: The shock scattering effect

Boiling histotripsy is a High Intensity Focused Ultrasound (HIFU) technique which uses a number of short pulses with high acoustic pressures at the HIFU focus to induce mechanical tissue fractionation. In boiling histotripsy, two different types of acoustic cavitation contribute towards mechanical tissue destruction: a boiling vapour bubble and cavitation clouds. An understanding of the mechanisms underpinning these phenomena and their dynamics is therefore paramount to predicting and controlling the overall size of a lesion produced for a given boiling histotripsy exposure condition. A number of studies have shown the effects of shockwave heating in generating a boiling bubble at the HIFU focus and have studied its dynamics under boiling histotripsy insonation. However, not much is known about the subsequent production of cavitation clouds that form between the HIFU transducer and the boiling bubble. The main objective of the present study is to examine what causes this bubble cluster formation after the generation of a boiling vapour bubble. A numerical simulation of 2D nonlinear wave propagation with the presence of a bubble at the focus of a HIFU field was performed using the k-Wave MATLAB toolbox for time domain ultrasound simulations, which numerically solves the generalised Westervelt equation. The numerical results clearly demonstrate the appearance of the constructive interference of a backscattered shockwave by a bubble with incoming incident shockwaves. This interaction (i.e., the reflected and inverted peak positive phase from the bubble with the incoming incident rarefactional phase) can eventually induce a greater peak negative pressure field compared to that without the bubble at the HIFU focus. In addition, the backscattered peak negative pressure magnitude gradually increased from 17.4 MPa to 31.6 MPa when increasing the bubble size from 0.2 mm to 1.5 mm. The latter value is above the intrinsic cavitation threshold of –28 MPa in soft tissue. Our results suggest that the formation of a cavitation cloud in boiling histotripsy is a threshold effect which primarily depends (a) the size and location of a boiling bubble, and (b) the sum of the incident field and that scattered by a bubble.     

Implementation of the Laser Diffraction Technique for Acoustic Cavitation Bubble Investigations

Knowledge of the acoustic cavitation cloud would be useful for improving ultrasound reactor design. Among the characterisation techniques, few are adapted to bubble investigations in an intense ultrasound field. Some problems raised by these measurements result from interactions between the acoustic pressure wave and the measuring light wave. This paper reports the implementation of the laser diffraction technique to determine the size and volume concentration of bubbles generated by a dipping horn operating at 20 kHz. Measurements were performed with a Malvern 2600 instrument. The size distribution, deduced from the diffraction pattern scattered by the bubble cloud crossed by a laser beam, is disturbed by the acoustic pressure wave involving deviation of a light beam at low diffusion angles (acousto‐optic effect). A bubble size correction procedure based on the subtraction of the light energy due to the ultrasound wave is described. The size measurements, and thus the correction procedure, were validated by a second laser technique based on a different measuring principle: phase Doppler interferometry. The measurement reliability was further confirmed by an original application of laser diffraction based on measurements performed just after sonication. These three methods lead to a mean bubble size (Sauter mean diameter) of about 10 μm at a high ultrasound power input. Concerning the void fraction, only measurements achieved after sonication and by laser diffraction predict a correct estimation of this parameter.     

气泡体积分数对沙质沉积物低频声学特性的影响

由于有机物质分解等原因,实际的海底沉积物中存在气泡,气泡的存在会显著影响沉积物低频段的声学特性,因此研究气泡对沉积物低频段声速的影响机理具有重要意义.考虑到外场环境的不可控性,在室内水池中搭建了大尺度含气非饱和沙质沉积物声学特性获取平台,在有界空间中应用多水听器反演方法首次获取了含气非饱和沙质沉积物300—3000 Hz频段内的声速数据(79—142 m/s),并同时利用双水听器法获取了同一频段的数据(112—121 m/s).在声波频率远低于沉积物中最大气泡的共振频率时,根据等效介质理论,将孔隙水和气泡等效为一种均匀流体,改进了水饱和等效密度流体近似模型.模型揭示了气泡对沉积物低频段声学特性的影响规律,理论上解释了沉积物中声速的降低.通过分析模型预报声速对模型参数的敏感性,根据测量得到的声速分布反演得到了沉积物不同区域的气泡体积分数,气泡体积分数从1.07%变化到2.81%.改进的模型为沉积物中气泡体积分数估计提供了一种新方法.     

Numerical simulation of the nonlinear ultrasonic pressure wave propagation in a cavitating bubbly liquid inside a sonochemical reactor

We investigate the acoustic wave propagation in bubbly liquid inside a pilot sonochemical reactor which aims to produce antibacterial medical textile fabrics by coating the textile with ZnO or CuO nanoparticles. Computational models on acoustic propagation are developed in order to aid the design procedures. The acoustic pressure wave propagation in the sonoreactor is simulated by solving the Helmholtz equation using a meshless numerical method. The paper implements both the state-of-the-art linear model and a nonlinear wave propagation model recently introduced by Louisnard (2012), and presents a novel iterative solution procedure for the nonlinear propagation model which can be implemented using any numerical method and/or programming tool. Comparative results regarding both the linear and the nonlinear wave propagation are shown. Effects of bubble size distribution and bubble volume fraction on the acoustic wave propagation are discussed in detail. The simulations demonstrate that the nonlinear model successfully captures the realistic spatial distribution of the cavitation zones and the associated acoustic pressure amplitudes.     

Optical nucleation of bubble clouds in a high pressure spherical resonator

An experimental setup for nucleating clouds of bubbles in a high-pressure spherical resonator is described. Using nanosecond laser pulses and multiple phase gratings, bubble clouds are optically nucleated in an acoustic field. Dynamics of the clouds are captured using a high-speed CCD camera. The images reveal cloud nucleation, growth, and collapse and the resulting emission of radially expanding shockwaves. These shockwaves are reflected at the interior surface of the resonator and then reconverge to the center of the resonator. As the shocks reconverge upon the center of the resonator, they renucleate and grow the bubble cloud. This process is repeated over many acoustic cycles and with each successive shock reconvergence, the bubble cloud becomes more organized and centralized so that subsequent collapses give rise to stronger, better defined shockwaves. After many acoustic cycles individual bubbles cannot be distinguished and the cloud is then referred to as a cluster. Sustainability of the process is ultimately limited by the detuning of the acoustic field inside the resonator. The nucleation parameter space is studied in terms of laser firing phase, laser energy, and acoustic power used.     

Multibubble sonoluminescence spectra of water which resemble single-bubble sonoluminescence

Multibubble sonoluminescence (MBSL) spectra of water from cavitation clouds were collected in the presence of different noble gases and at different acoustic intensities. Results show that at high acoustic intensity and with xenon as a dissolved gas the emission of the OH* radical becomes indiscernible from the continuum. These spectra resemble single-bubble sonoluminescence (SBSL) spectra. It is concluded that the source of emission in MBSL and SBSL can be the same, the difference in spectra is due to the higher temperature inside the bubble during SBSL.     

Bubble size distribution in acoustic droplet vaporization via dissolution using an ultrasound wide-beam method

Performance and efficiency of numerous cavitation enhanced applications in a wide range of areas depend on the cavitation bubble size distribution. Therefore, cavitation bubble size estimation would be beneficial for biological and industrial applications that rely on cavitation. In this study, an acoustic method using a wide beam with low pressure is proposed to acquire the time intensity curve of the dissolution process for the cavitation bubble population and then determine the bubble size distribution. Dissolution of the cavitation bubbles in saline and in phase-shift nanodroplet emulsion diluted with undegassed or degassed saline was obtained to quantify the effects of pulse duration (PD) and acoustic power (AP) or peak negative pressure (PNP) of focused ultrasound on the size distribution of induced cavitation bubbles. It was found that an increase of PD will induce large bubbles while AP had only a little effect on the mean bubble size in saline. It was also recognized that longer PD and higher PNP increases the proportions of large and small bubbles, respectively, in suspensions of phase-shift nanodroplet emulsions. Moreover, degassing of the suspension tended to bring about smaller mean bubble size than the undegassed suspension. In addition, condensation of cavitation bubble produced in diluted suspension of phase-shift nanodroplet emulsion was involved in the calculation to discuss the effect of bubble condensation in the bubble size estimation in acoustic droplet vaporization. It was shown that calculation without considering the condensation might underestimate the mean bubble size and the calculation with considering the condensation might have more influence over the size distribution of small bubbles, but less effect on that of large bubbles. Without or with considering bubble condensation, the accessible minimum bubble radius was 0.4 or 1.7 μm and the step size was 0.3 μm. This acoustic technique provides an approach to estimate the size distribution of cavitation bubble population in opaque media and might be a promising tool for applications where it is desirable to tune the ultrasound parameters to control the size distribution of cavitation bubbles.     

Spatial resolution limits for the reconstruction of acoustic source strength by inverse methods

One method for deducing the strength of an acoustic source distribution from measurement of the radiated field involves the inversion of the matrix of frequency response functions relating the field measurement points to the strengths of a number of point sources used to represent the source distribution. In practice, the frequency response function matrix to be inverted may very often be ill-conditioned. This ill-conditioning will also often result in an ill-posed problem and thus regularization algorithms are used to produce reasonable solutions. For this purpose, Tikhonov regularization has been applied, and generalized cross-validation (GCV) has been introduced as an effective method for determining the proper amount of regularization without prior knowledge of either the source distribution or the contaminating errors. In the present work, the emphasis is placed on the relationship between the spatial resolution of the reconstructed source distribution and the small singular values of the frequency response function matrix to be inverted. However, the use of Tikhonov regularization often suppresses the effect of small singular values and these are in turn often associated with high spatial frequencies of the source distribution. Thus, the process of regularization produces a useful estimate of the acoustic source strength distribution but with a limited spatial resolution. Furthermore, in the field of Fourier acoustics, the spatial resolution of the reconstructed source distribution is usually limited by the wavelength of the radiation. This paper expresses the relationship between estimation accuracy, spatial resolution, noise-level and source/sensor geometry, when a range of inverse sound radiation problems are regularised using Tikhonov regularization with GCV. The results presented form the basis of guidelines that enable the reconstruction of acoustic source strength with a resolution that is finer than the intrinsic half-wavelength limit.     

Sonoluminescence and dynamics of cavitation bubble populations in sulfuric acid

The detailed link of liquid phase sonochemical reactions and bubble dynamics is still not sufficiently known. To further clarify this issue, we image sonoluminescence and bubble oscillations, translations, and shapes in an acoustic cavitation setup at 23 kHz in sulfuric acid with dissolved sodium sulfate and xenon gas saturation. The colour of sonoluminescence varies in a way that emissions from excited non-volatile sodium atoms are prominently observed far from the acoustic horn emitter (“red region”), while such emissions are nearly absent close to the horn tip (“blue region”). High-speed images reveal the dynamics of distinct bubble populations that can partly be linked to the different emission regions. In particular, we see smaller strongly collapsing spherical bubbles within the blue region, while larger bubbles with a liquid jet during collapse dominate the red region. The jetting is induced by the fast bubble translation, which is a consequence of acoustic (Bjerknes) forces in the ultrasonic field. Numerical simulations with a spherical single bubble model reproduce quantitatively the volume oscillations and fast translation of the sodium emitting bubbles. Additionally, their intermittent stopping is explained by multistability in a hysteretic parameter range. The findings confirm the assumption that bubble deformations are responsible for pronounced sodium sonoluminescence. Notably the observed translation induced jetting appears to serve as efficient mixing mechanism of liquid into the heated gas phase of collapsing bubbles, thus potentially promoting liquid phase sonochemistry in general.     

丛聚的含气泡水对线性声传播的影响

含气泡水内气泡的空间分布会对线性声传播产生影响,导致实验结论与理论预测存在较大偏差.为解决这一问题,将准晶体近似引入到自洽方法中,导出了考虑空间分布时多分散含气泡水的等效声波波数.考虑到含气泡水内,气泡间存在小范围的聚集趋势(简称丛聚现象),在此基础上引入Neyman-Scott点过程描述了含气泡水内气泡的丛聚现象.分析发现,丛聚时,声速、声衰减的峰值将受到抑制,并向低频偏移,且抑制和频偏现象会随丛聚加剧而变强;随频率远离峰值段,丛聚对声传播的影响逐渐减弱.此外,考虑到空间分布的统计信息提取对相关研究的精确与否起到重要作用,引入了一种比例无偏估计,通过该方法获得了仿真环境下丛聚含气泡水模型的相速度及衰减系数,该建模及统计方法也可为相关实验工作提供理论基础.     

Radial oscillation and translational motion of a gas bubble in a micro-cavity

According to classical nucleation theory, a gas nucleus can grow into a cavitation bubble when the ambient pressure is negative. Here, the growth process of a gas nucleus in a micro-cavity was simplified to two “events”, and the full confinement effect of the surrounding medium of the cavity was considered by including the bulk modulus in the equation of state. The Rayleigh–Plesset-like equation of the cavitation bubble in the cavity was derived to model the radial oscillation and translational motion of the cavitation bubble in the local acoustic field. The numerical results show that the nucleation time of the cavitation bubble is sensitive to the initial position of the gas nucleus. The cavity size affects the duration of the radial oscillation of the cavitation bubble, where the duration is shorter for smaller cavities. The equilibrium radius of a cavitation bubble grown from a gas nucleus increases with increasing size of the cavity. There are two possible types of translational motion: reciprocal motion around the center of the cavity and motion toward the cavity wall. The growth process of gas nuclei into cavitation bubbles is also dependent on the compressibility of the surrounding medium and the magnitude of the negative pressure. Therefore, gas nuclei in a liquid cavity can be excited by acoustic waves to form cavitation bubbles, and the translational motion of the cavitation bubbles can be easily observed owing to the confining influence of the medium outside the cavity.