弹性微管内气泡的非线性受迫振动

将弹性管壁视为膜弹性结构, 探索在外部声场作用下弹性微管内液柱-气泡-管壁构成耦合振动系统的非线性特征. 利用逐级近似法对系统非线性共振频率、基频和三倍频振动幅值响应、 分频激励共振机理等进行了理论分析. 基频和三倍频振动的幅-频响应数值结果表明: 气泡的轴向共振和管壁共振不能同时出现; 两垂直方向的振动均表现出幅值响应多值性, 进而可能引起系统的不稳定声响应; 三倍频振动在低频区的声响应强于高频区. 关键词： 弹性微管 受迫振动 非线性振动 气泡声响应     

A two-frequency acoustic technique for bubble resonant oscillation studies

A two-frequency acoustic apparatus has been developed to study the dynamics of a single gas or vapor bubble in water. An advantage of the apparatus is its capability of trapping a bubble by an ultrasonic standing wave while independently driving it into oscillations by a second lower frequency acoustic wave. For a preliminary application, the apparatus is used to study resonant oscillations. First, near-resonant coupling between the volume and the n = 3 shape oscillation modes of air bubbles at room temperature is studied, where n is the mode number. The stability boundary, amplitude versus frequency, of the volume oscillation forms a wedge centered at the resonant frequency, which qualitatively agrees with a theoretical prediction based on a phase-space analysis. Next, the resonant volume oscillations of vapor bubbles are studied. The resonant radius of vapor bubbles at 80 degrees C driven at 1682 Hz is determined to be 0.7 mm, in agreement with a prediction obtained by numerical simulation.     

微管内气泡的受迫振动

在气泡-液柱一维耦合振动模型的基础上对刚性微管两侧声压不相等时管内柱状气泡的轴向一维受迫振动进行了理论探索. 声压不均匀分布不影响气泡线性振动时的共振频率, 但振动幅度受到有效声压幅值的影响. 利用逐级近似法分析了管内非线性振动气泡的基频、三倍频和三分之一分频振动的幅-频响应关系, 结果表明当驱动声压超过0.1 MPa时, 气泡振动处于非线性状态. 非线性声响应特征主要表现为:基频和分频振动幅值响应的多值性; 三倍频振动在低频区响应强于高频区; 三分频振动在大于共振频率的频域内出现的概率更大.     

Nonlinear non-spherical oscillations of gas bubble under a periodic variation of ambient liquid pressure

A mathematical model is constructed for the bubble dynamics, in which the interphase surface variation is presented in the form of a series in spherical harmonics, and the equations are written with the accuracy up to the squared amplitude of the distortion of the spherical shape of the bubble. In the oscillation regimes close to periodic sonoluminescence of a single bubble in a standing acoustic wave, the character of air bubble oscillations in water was studied depending on the bubble initial radius and the amplitude of the liquid pressure variation. It was found that non-spherical oscillations of bounded amplitude can take place outside the region of linearly stable spherical oscillations. Both the oscillations with a period equal to one or several periods of the liquid pressure variation and aperiodic oscillations are observed. It is shown that neglecting the distortions in the form of spherical harmonics with large numbers (i > 3) may lead to a change of oscillation regimes. The influence of distortions on the bubble surface shape for the harmonics with i > 8 is insignificant.     

声场中空化气泡的耦合振动及形状不稳定性的研究

计算了两个具有非球形扰动的气泡所组成系统的能量,并基于Lagrange方程得到了有声相互作用的非球形气泡的动力学方程和形状稳定性方程,研究了声场中非球形气泡间相互作用力对非球形气泡的形状不稳定性和气泡形状模态振幅的影响.研究结果表明声场中具有非球形扰动的气泡之间的耦合方式有两种:形状耦合模式和径向耦合模式,气泡之间的耦合方式取决于气泡形状扰动模态.由形状耦合及径向耦合产生的气泡之间的相互作用力能够改变单个气泡的形状不稳定及形状模态振幅,具体影响因素取决于声场驱动条件、气泡形状模态、相邻气泡的初始半径.     

球状泡群内气泡的耦合振动

振动气泡形成辐射场影响其他气泡的运动, 故多气泡体系中气泡处于耦合振动状态. 本文在气泡群振动模型的基础上, 考虑气泡间耦合振动的影响, 得到了均匀球状泡群内振动气泡的动力学方程, 以此为基础分析了气泡的非线性声响应特征. 气泡间的耦合振动增加了系统对每个气泡的约束, 降低了气泡的自然共振频率, 增强了气泡的非线性声响应. 随着气泡数密度的增加, 振动气泡受到的抑制增强; 增加液体静压力同样可抑制泡群内气泡的振动, 且存在静压力敏感区(1–2 atm, 1 atm=1.01325×10⁵ Pa); 驱动声波对气泡振动影响很大, 随着声波频率的增加, 能够形成空化影响的气泡尺度范围变窄. 在同样的声条件、泡群尺寸以及气泡内外环境下, 初始半径小于5 μm 的气泡具有较强的声响应. 气泡耦合振动会削弱单个气泡的空化影响, 但可延长多气泡系统空化泡崩溃发生的时间间隔和增大作用范围, 整体空化效应增强.     

Nonlinear oscillation and acoustic scattering of bubbles

The scattered acoustic pressure and scattered cross section of bubbles is studied using the scattered theory of bubbles. The nonlinear oscillations of bubbles and the scattering acoustic fields of a spherical bubble cluster are numerically simulated based on the bubble dynamic and fluid dynamic. The influences of the interaction between bubbles on scattering acoustic field of bubbles are researched. The results of numerical simulation show that the oscillation phases of bubbles are delayed to a certain extent at different positions in the bubble cluster, but the radii of bubbles during oscillation do not differ too much at different positions. Furthermore, directivity of the acoustic scattering of bubbles is obvious. The scattered acoustic pressures of bubbles are different at the different positions inside and outside of the bubble cluster. The scattering acoustic fields of a spherical bubble cluster depend on the driving pressure amplitude, driving frequency, the equilibrium radii of bubbles, bubble number and the radius of the spherical bubble cluster. These theoretical predictions provide a further understanding of physics behind ultrasonic technique and should be useful for guiding ultrasonic application.     

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.     

Forced linear oscillations of microbubbles in blood capillaries

A theoretical investigation of the forced linear oscillations of a gas microbubble in a blood capillary, whose radius is comparable in size to the bubble radius is presented. The natural frequency of oscillation, the thermal and viscous damping coefficients, the amplitude resonance, the energy resonance, as well as the average energy absorbed by the system, bubble plus vessel, have been computed for different kinds of gas microbubbles, containing air, octafluropropane, and perflurobutane as a function of the bubble radius and applied frequency. It has been found that the bubble behavior is isothermal at low frequencies and for small bubbles and between isothermal and adiabatic for larger bubbles and higher frequencies, with the viscous damping dominating over the thermal damping. Furthermore, the width of the energy resonance is strongly dependent on the bubble size and the natural frequency of oscillation is affected by the presence of the vessel wall and position of the bubble in the vessel. Therefore, the presence of the blood vessel affects the way in which the bubble absorbs energy from the ultrasonic field. The motivation of this study lies in the possibility of using gas microbubbles as an aid to therapeutic focused ultrasound treatments.     

Numerical investigation of the effect of dual frequency sonication on stable bubble dynamics

A computational study aiming to simulate an oxygen single acoustic bubble oscillation under a dual-frequency sonication was presented in this paper. The non-linear response of the bubble to the superposition of two fields of ultrasonic waves was investigated through dynamics parameters, collapse ratios and average velocities. The main goal of this analyze is to link the properties of the wave resulting from the dual-frequency excitation to the dynamics behavior of the bubble. The obtained results prove that, in contrast with the mono-frequency, coupling a wave to lower frequencies enhances the collapse duration and raises the compression ratio in the case of 35 kHz, while associating any of the studied waves to a higher frequency elevates the number of bubble oscillations during a time interval as compared to mono-frequency. The total sonochemical production has been investigated in accordance with the dynamics results, as well as the proportions of the three predominant free radicals, that show a dependency on the value of the basic frequency.     

Influence of heat transfer on the heating of an oscillating bubble wall and on the transition from periodic oscillations to chaotic ones

A simple scheme is proposed for estimating both the heating of the liquid at the wall of a bubble performing radial oscillations and the influence of different factors on the oscillation amplitude and phase. The driving force amplitude corresponding to the transition from periodic bubble oscillations to chaotic ones proves to be independent of whether the heat conduction is taken into account or not.     

Radial and translational oscillations of an acoustically levitated bubble in aqueous ethanol solutions

The radial and translational oscillations of a single cavitation bubble in a standing ultrasound wave were investigated experimentally at various driving acoustic pressures for aqueous ethanol solutions with different bulk molar fractions of ethanol range of 0-1.3 × 10(-3). The results show that both the lower and upper stability thresholds of the acoustic driving pressure decreased as the concentration of ethanol was increased. At a given driving pressure the ambient and maximum bubble sizes increased with increasing ethanol concentration. In addition, as the ethanol was increased, the sonoluminescence intensity decreased while the bubble dynamics remained largely unchanged. The translational oscillation of the levitated bubble, however, became increasingly violent with increasing ethanol concentration. The displacement of the bubble reached 0.7 mm at the highest concentration studied (1.3 × 10(-3)) and the maximum bubble size was found to change as the bubble jumped up and down. This bubble translation may be responsible for the decrease of the acoustic driving pressure threshold and suggests that repetitive injection of ethanol molecules into the bubble takes place. These results may account for the different sensitivities of single bubble and multi-bubble sonoluminescence to the presence of volatile additives.     

Enhancement of oscillation amplitude of cavitation bubble due to acoustic wake effect in multibubble environment

Acoustic cavitation occurs in ultrasonic treatment causing various phenomena such as chemical synthesis, chemical decomposition, and emulsification. Nonlinear oscillations of cavitation bubbles are assumed to be responsible for these phenomena, and the neighboring bubbles may interact each other. In the present study, we numerically investigated the dynamic behavior of cavitation bubbles in multi-bubble systems. The results reveal that the oscillation amplitude of a cavitation bubble surrounded by other bubbles in a multi-bubble system becomes larger compared with that in the single-bubble case. It is found that this is caused by an acoustic wake effect, which reduces the pressure near a bubble surrounded by other bubbles and increases the time delay between the bubble contraction/expansion cycles and sound pressure oscillations. A new parameter, called “cover ratio” is introduced to quantitatively evaluate the variation in the bubble oscillation amplitude, the time delay, and the maximum bubble radius.     

限位液滴瞬时失重自激振荡

为探索重力瞬变引起的约束液滴自激振荡机理, 本文利用落塔装置模拟短时微重力环境并借助高速CCD记录圆形限位基片上液滴整个过程的运动情况. 自激振荡是微重力下液滴形态的重整恢复过程, 边界的限位作用使得液滴在整个运动过程接触线钉扎不变, 具体可分为两个阶段: 首先是振荡的高低点位置高度渐进上升的液滴形态变化阶段, 与重力环境渐进变化有关; 而后是平衡位置附近的阻尼衰减振荡阶段, 此时振荡的频率恒定, 振幅衰减类似孤立黏性液滴的指数衰减过程. 对于第二阶段, 在高低点等位置处存在高度不变过程, 高度起伏变化时液滴振荡模式类似自由液滴二阶振荡, 高度不变时振荡模式类似自由液滴三阶振荡. 此外, 对于本实验体系的恒定接触面积的钉扎约束, 液滴的体积量不同时, 内驱振荡的阶段和模式不变, 但具体的振荡过程有所不同. 对于大体积液滴, 会在初始振荡的中间位置出现高度不变现象, 并且随振荡逐渐消失; 而小液滴中间位置则不存在此现象, 波形较一致; 第二阶段小体积液滴振幅衰减的阻尼率更大, 无量纲频率也更高.     

一种水中爆炸气泡脉动实验研究方法

 在2 m×2 m×2 m的实验水箱中开展小当量PETN炸药水中爆炸气泡脉动实验时,利用弹性波从声阻抗高的物质传入声阻抗低的物质时的“减震缓冲”原理,采用在水箱壁贴低声阻抗材料的方法,有效降低了水箱壁反射冲击波对气泡脉动过程的影响,获得了清晰的气泡脉动过程图像、气泡水射流形成过程图像和气泡脉动压力曲线。将水箱实验结果与8 kg TNT当量爆炸水池实验结果对比,得到的最大相对误差仅为5.1%,验证了水箱实验方法的准确性,为炸药水中爆炸气泡脉动现象研究提供了一种简便有效的实验方法。     

Shape stability of a gas bubble in a soft solid

Predicting the onset of non-spherical oscillations of bubbles in soft matter is a fundamental cavitation problem with implications to sonoprocessing, polymeric materials synthesis, and biomedical ultrasound applications. The shape stability of a bubble in a Kelvin-Voigt viscoelastic medium with nonlinear elasticity, the simplest constitutive model for soft solids, is analytically investigated and compared to experiments. Using perturbation methods, we develop a model reducing the equations of motion to two sets of evolution equations: a Rayleigh-Plesset-type equation for the mean (volume-equivalent) bubble radius and an equation for the non-spherical mode amplitudes. Parametric instability is predicted by examining the natural frequency and the Mathieu equation for the non-spherical modes, which are obtained from our model. Our theoretical results show good agreement with published experiments of the shape oscillations of a bubble in a gelatin gel. We further examine the impact of viscoelasticity on the time evolution of non-spherical mode amplitudes. In particular, we find that viscosity increases the damping rate, thus suppressing the shape instability, while shear modulus increases the natural frequency, which changes the unstable mode. We also explain the contributions of rotational and irrotational fields to the viscoelastic stresses in the surroundings and at the bubble surface, as these contributions affect the damping rate and the unstable mode. Our analysis on the role of viscoelasticity is potentially useful to measure viscoelastic properties of soft materials by experimentally observing the shape oscillations of a bubble.     

Energy levels of oscillations in a nonlinear transmission line filled with saturated ferrite

Stationary solutions in the form of traveling electromagnetic waves in a uniform transmission line filled with saturated ferrite in the absence of damping are studied. Expressions are obtained for the period of nonlinear oscillations as a function of their amplitude. The period of linear oscillations is shown to linearly increase with the current through the line. The group velocity of the stationary wave increases with the oscillation amplitude. A relationship is found between the electromagnetic oscillations and oscillations of the magnetization vector. The oscillations become increasingly anisochronous with increasing the external magnetic field.     

Resonance frequencies of lipid-shelled microbubbles in the regime of nonlinear oscillations

Knowledge of resonant frequencies of contrast microbubbles is important for the optimization of ultrasound contrast imaging and therapeutic techniques. To date, however, there are estimates of resonance frequencies of contrast microbubbles only for the regime of linear oscillation. The present paper proposes an approach for evaluating resonance frequencies of contrast agent microbubbles in the regime of nonlinear oscillation. The approach is based on the calculation of the time-averaged oscillation power of the radial bubble oscillation. The proposed procedure was verified for free bubbles in the frequency range 1-4 MHz and then applied to lipid-shelled microbubbles insonified with a single 20-cycle acoustic pulse at two values of the acoustic pressure amplitude, 100 kPa and 200 kPa, and at four frequencies: 1.5, 2.0, 2.5, and 3.0 MHz. It is shown that, as the acoustic pressure amplitude is increased, the resonance frequency of a lipid-shelled microbubble tends to decrease in comparison with its linear resonance frequency. Analysis of existing shell models reveals that models that treat the lipid shell as a linear viscoelastic solid appear may be challenged to provide the observed tendency in the behavior of the resonance frequency at increasing acoustic pressure. The conclusion is drawn that the further development of shell models could be improved by the consideration of nonlinear rheological laws.     

Theoretical investigation on resonance characteristics of a vapor bubble based on Laplace transform method

In the present paper, resonance characteristics of the vapor bubble oscillating in an acoustic field are investigated analytically. The analytical solution of the non-dimensional perturbation of the instantaneous bubble radius during the transient process in the initial oscillation stage is explicitly obtained and physically analyzed at the resonance situation based on the Laplace transform method. And the typical oscillation behaviors obtained from the analytical solution are thoroughly exhibited and analyzed in the time and frequency domains. In addition, the corresponding oscillation behaviors at the non-resonance situation are also investigated for the purpose of comparisons. Through our investigation, several essential conclusions can be drawn as follows: (1) The analytical solution of the non-dimensional perturbation of the instantaneous bubble radius can be divided into four terms according to the physical meaning. Among them, it is the term related to the acoustic field that causes the progressively violent bubble oscillation. (2) The vapor bubble with a smaller equilibrium radius could respond faster and more significantly to the acoustic field during the oscillation. (3) The bubble oscillation characteristics always exhibit significant differences at the resonance and non-resonance situations in both the time and frequency domains, even if the difference between the natural frequency of the oscillating vapor bubble and the angular frequency of the acoustic field is greatly small.     

Dependence of the characteristics of bubbles on types of sonochemical reactors

Computer simulations of bubble oscillations in liquid water irradiated by an ultrasonic wave have revealed that the characteristic of bubbles depends on types of sonochemical reactors: a horn-type reactor and a standing-wave type reactor. When the acoustic amplitude is large at 20 kHz, the bubble content is mostly water vapor even at the end of the bubble collapse and the temperature inside a bubble at the collapse is relatively low. On the other hand, when the acoustic amplitude is relatively low, the bubble content is mostly noncondensable gas at the end of the bubble collapse and the bubble temperature is relatively high. In a horn-type sonochemical reactor, the former type of bubbles are dominant because many bubbles exist near the horn-tip where the acoustic amplitude is large, while in a standing-wave type reactor the latter type of bubbles are dominant because the Bjerknes force gathers bubbles at a region where acoustic amplitude is relatively low.