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

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

一维非线性声波传播特性

针对一维非线性声波的传播问题进行了有限元仿真和实验研究. 首先推导了一维非线性声波方程的有限元形式, 含有高阶矩阵的非线性项导致声波具有波形畸变、谐波滋生、基频信号能量向高次谐波传递等非线性特性. 编制有限元程序对一维非线性声波进行了计算并对仿真得到的畸变非线性声波信号进行处理, 分析其传播性质和物理意义. 为验证有限元计算结果, 开展了水中的非线性声波传播的实验研究, 得到了不同输入信号幅度激励下和不同传播距离的畸变非线性声波信号. 然后对基波和二次谐波的传播性质进行详细讨论, 分析了二次谐波幅度与传播距离和输入信号幅度的变化关系及其意义, 拟合出二次谐波幅度随传播距离变化的方程并阐述了拟合方程的物理意义. 结果表明, 数值仿真信号及其频谱均与实验结果有较好的一致性, 证实计算方法和结果的正确性, 并提出了具有一定物理意义的二次谐波随传播距离变化的简单数学关系. 最后还对固体中的非线性声波传播性质进行了初步探讨. 本研究工作可为流体介质中的非线性声传播问题提供理论和实验依据.     

Contrast harmonic imaging

The behavior of ultrasound contrast agents depends highly on the acoustic pressure of the insonified ultrasound wave. For low pressure the expansion and compression is linear to the pressure, for medium acoustic pressure nonlinear behavior starts to occur and for high pressures, but still in the diagnostic range transient scattering can be noticed, resulting in an enhanced scattering followed by a disappearance of the bubble. The nonlinear and transient regime can be utilized for imaging of the contrast agent in or nearby tissue. The magnitude of the nonlinear signal from the contrast has to compete with the nonlinear component of the ultrasound wave, which is generated during propagation. It is shown that contrast is superior to tissue when using low frequencies and imaging the third or fourth harmonic of the transmitted frequency.     

Transfer functions of US transducers for harmonic imaging and bubble responses

Current medical diagnostic echo systems are mostly using harmonic imaging. This means that a fundamental frequency (e.g., 2 MHz) is transmitted and the reflected and scattered higher harmonics (e.g., 4 and 6 MHz), produced by nonlinear propagation, are recorded. The signal level of these harmonics is usually low and a well-defined transfer function of the receiving transducer is required. Studying the acoustic response of a single contrast bubble, which has an amplitude in the order of a few Pascal, is another area where an optimal receive transfer function is important.

We have developed three methods to determine the absolute transfer function of a transducer. The first is based on a well-defined wave generated by a calibrated source in the far field. The receiving transducer receives the calibrated wave and from this the transfer functions can be calculated. The second and third methods are based on the reciprocity of the transducer. The second utilizes a calibrated hydrophone to measure the transmitted field. In the third method, a pulse is transmitted by the transducer, which impinges on a reflector and is received again by the same transducer. In both methods, the response combined with the transducer impedance and beam profiles enables the calculation of the transfer function.

The proposed methods are useful to select the optimal piezoelectric material (PZT, single crystal) for transducers used in reception only, such as in certain 3D scanning designs and superharmonic imaging, and for selected experiments like single bubble behavior.

We tested and compared these methods on two unfocused single element transducers, one commercially available (radius 6.35 mm, centre frequency 2.25 MHz) the other custom built (radius 0.75 mm, centre frequency 4.3 MHz). The methods were accurate to within 15%.     

A simple model of ultrasound propagation in a cavitating liquid. Part I: Theory, nonlinear attenuation and traveling wave generation

The bubbles involved in sonochemistry and other applications of cavitation oscillate inertially. A correct estimation of the wave attenuation in such bubbly media requires a realistic estimation of the power dissipated by the oscillation of each bubble, by thermal diffusion in the gas and viscous friction in the liquid. Both quantities and calculated numerically for a single inertial bubble driven at 20 kHz, and are found to be several orders of magnitude larger than the linear prediction. Viscous dissipation is found to be the predominant cause of energy loss for bubbles small enough. Then, the classical nonlinear Caflish equations describing the propagation of acoustic waves in a bubbly liquid are recast and simplified conveniently. The main harmonic part of the sound field is found to fulfill a nonlinear Helmholtz equation, where the imaginary part of the squared wave number is directly correlated with the energy lost by a single bubble. For low acoustic driving, linear theory is recovered, but for larger drivings, namely above the Blake threshold, the attenuation coefficient is found to be more than 3 orders of magnitude larger then the linear prediction. A huge attenuation of the wave is thus expected in regions where inertial bubbles are present, which is confirmed by numerical simulations of the nonlinear Helmholtz equation in a 1D standing wave configuration. The expected strong attenuation is not only observed but furthermore, the examination of the phase between the pressure field and its gradient clearly demonstrates that a traveling wave appears in the medium.     

Improved contrast to tissue ratio at higher harmonics

The challenge in ultrasound contrast imaging is a better discrimination between the perfused tissue and the contrast bubbles, which is usually expressed by contrast to tissue ratio (CTR). Imaging based on the second harmonic frequency showed a higher CTR than imaging at the fundamental frequency. However, because of nonlinear propagation of ultrasound waves, harmonic frequencies are generated. These harmonic frequencies will be linearly reflected by the tissue and therefore limit the CTR at the second harmonic frequency. In order to reduce the scattering of tissue at harmonic frequencies and by that increase the CTR, nonlinear distortion has to be reduced. We demonstrate in this study that the CTR increases with the harmonic number. The increase is substantial when transmitting at lower frequencies. To take advantage of the higher harmonics (third, fourth, fifth and the ultraharmonics and termed here super harmonics), we have developed a new phased array transducer with a wide frequency band. In-vitro measurements using the new probe show an increase of 40 dB of the CTR for super harmonic components over the conventional second harmonic system. The increase in CTR is in agreement with the calculations using existing models for the response of encapsulated bubbles and known theory of nonlinear propagation.     

磁绝缘线振荡器中空间电荷的调制

 首次研究了磁绝缘线振荡器中射频场（包括辐射波场和空间电荷波场）对空间电荷的调制,得到了在饱和时谐波电流一阶分量的公式。这个公式表明：考虑空间电荷波时,在小信号和辐射波与空间电荷波同相的条件下,它将增强谐波电流。在大信号时,情况不确定;谐波电流一阶分量将随二极管上电压的增加而增加,随运行频率和饱和长度的增加而减少。     

Bubble nonlinear dynamics and stimulated scattering process

A complete understanding of the bubble dynamics is deemed necessary in order to achieve their full potential applications in industry and medicine. For this purpose it is first needed to expand our knowledge of a single bubble behavior under different possible conditions including the frequency and pressure variations of the sound field. In addition, stimulated scattering of sound on a bubble is a special effect in sound field, and its characteristics are associated with bubble oscillation mode. A bubble in liquid can be considered as a representative example of nonlinear dynamical system theory with its resonance, and its dynamics characteristics can be described by the Keller–Miksis equation. The nonlinear dynamics of an acoustically excited gas bubble in water is investigated by using theoretical and numerical analysis methods. Our results show its strongly nonlinear behavior with respect to the pressure amplitude and excitation frequency as the control parameters, and give an intuitive insight into stimulated sound scattering on a bubble. It is seen that the stimulated sound scattering is different from common dynamical behaviors, such as bifurcation and chaos, which is the result of the nonlinear resonance of a bubble under the excitation of a high amplitude acoustic sound wave essentially. The numerical analysis results show that the threshold of stimulated sound scattering is smaller than those of bifurcation and chaos in the common condition.     

光学微腔中倍频光场演化和光谱特性

在考虑光学微腔中二阶和三阶非线性效应的情况下,引入了可同时描述腔内基频和倍频光场的演化过程的Lugiato-Lefeve方程,分析了SiN微腔中二次谐波的产生,并讨论了各参数对腔内基频和倍频光场的影响.理论分析结果表明,失谐参量为0时,稳定后的基频光场为平顶脉冲的形式,而倍频光场呈正弦分布;失谐参量增加,将导致腔内基频和倍频光功率在演化过程中出现振荡,且最终稳定的光功率变弱,稳定后的光场分布为周期性变化;失谐参量的值过大,会使得微腔光场处于混沌状态.抽运光强较弱时,腔内可形成稳定的光场分布;抽运光强较强时,会导致腔内色散以及非线性效应过强,最终稳定的光场仍然呈周期性变化,且抽运光功率越强,光功率的演化曲线振荡越强.此外,选取特定的微腔尺寸,微腔可工作于"图灵环"状态.理论分析结果对研究光学微腔中二次谐波的产生有重要意义.     

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

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

Cumulative solutions of nonlinear longitudinal vibration in isotropic solid bars

Based on the strain invariant relationship and taking the high-order elastic energy into account, a nonlinear wave equation is derived, in which the excitation, linear damping, and the other nonlinear terms are regarded as the first-order correction to the linear wave equation. To solve the equation, the biggest challenge is that the secular terms exist not only in the fundamental wave equation but also in the harmonic wave equation （unlike the Duffing oscillator, where they exist only in the fundamental wave equation）. In order to overcome this difficulty and to obtain a steady periodic solution by the perturbation technique, the following procedures are taken： （i） for the fundamental wave equation, the secular term is eliminated and therefore a frequency response equation is obtained; （ii） for the harmonics, the cumulative solutions are sought by the Lagrange variation parameter method. It is shown by the results obtained that the second- and higher-order harmonic waves exist in a vibrating bar, of which the amplitude increases linearly with the distance from the source when its length is much more than the wavelength; the shift of the resonant peak and the amplitudes of the harmonic waves depend closely on nonlinear coefficients; there are similarities to a certain extent among the amplitudes of the odd- （or even-） order harmonics, based on which the nonlinear coefficients can be determined by varying the strain and measuring the amplitudes of the harmonic waves in different locations.     

两种气泡混合的声空化

将非线性声波方程和改进的Rayleigh-Plesset方程联立可以描述空化环境中的声场及相应的气泡动力学特征. 用时域有限差分方法模拟了圆柱形容器内两种气泡相互混合时的空化情况. 在烧杯内的稳态背景声场形成过程中, 瓶壁耗散吸收扮演了重要的角色. 在稳态背景声场的基础上, 分析了混合气泡与声场的相互作用、气泡之间的相互作用、混合情况下的频谱特性. 结果表明: 两种气泡平衡半径都不太大时, 气泡与声场的相互作用不强, 声场及气泡的行为也比较规律; 相反, 当其中一种气泡平衡半径相对比较大时, 声场与气泡具有较强的非线性相互作用, 声场及气泡的行为表现出复杂的特性.     

声波在水-含气沉积物界面的反射

含气泡海洋沉积物的声学特性是海底探测的重要问题。为了研究气泡存在对水-含气沉积物界面声反射系数的影响,本文基于气泡振动修正的Biot波动方程推导了气泡存在修正的Biot弹性模量,并结合水-沉积物界面的“开孔”边界条件推导了声波从水入射到水-含气沉积物界面的反射系数。数值分析表明气泡的振动导致反射系数呈现显著的频率特性。在气泡共振频率附近,由于气泡的共振引发的强散射和强衰减,使得反射系数很大,无论以何种角度入射,声波都很难进入含气泡的沉积物。本文研究结果表明,气泡半径、含量、声波频率以及入射角度都是影响水-含气沉积物界面反射系数的主要因素。      

"Compression-only" behavior: a second-order nonlinear response of ultrasound contrast agent microbubbles

Oscillating phospholipid-coated ultrasound contrast agent microbubbles display a so-called "compression-only" behavior, where it is observed that the bubbles compress efficiently while their expansion is suppressed. Here, a theoretical understanding of the source of this nonlinear behavior is provided through a weakly nonlinear analysis of the shell buckling model proposed by Marmottant et al. [J. Acoust. Soc. Am. 118, 3499-3505 (2005)]. It is shown that the radial dynamics of the bubble can be considered as a superposition of a linear response at the fundamental driving frequency and a second-order nonlinear low-frequency response that describes the negative offset of the mean bubble radius. The analytical solution deduced from the weakly nonlinear analysis shows that the compression-only behavior results from a rapid change of the shell elasticity with bubble radius. In addition, the radial dynamics of single phospholipid-coated microbubbles was recorded as a function of both the amplitude and the frequency of the driving pressure pulse. The comparison between the experimental data and the theory shows that the magnitude of compression-only behavior is mainly determined by the initial phospholipids concentration on the bubble surface, which slightly varies from bubble to bubble.     

Transient analysis of nonlinear Euler–Bernoulli micro-beam with thermoelastic damping,via nonlinear normal modes

In this paper an Euler–Bernoulli model has been used for vibration analysis of micro-beams with large transverse deflection. Thermoelastic damping is considered to be the dominant damping mechanism and introduced as imaginary stiffness into the equation of motion by evaluating temperature profile as a function of lateral displacement. The obtained equation of motion is analyzed in the case of pure single mode motion by two methods; nonlinear normal mode theory and the Galerkin procedure. In contrast with the Galerkin procedure, nonlinear normal mode analysis introduces a nonconventional nonlinear damping term in modal oscillator which results in strong damping in case of large amplitude vibrations. Evaluated modal oscillators are solved using harmonic balance method and tackling damping terms introduced as an imaginary stiffness is discussed. It has been shown also that nonlinear modal analysis of micro-beam with thermoelastic damping predicts parameters such as inverse quality factor, and frequency shift, to have an extrema point at certain amplitude during transient response due to the mentioned nonlinear damping term; and the effect of system?s characteristics on this critical amplitude has also been discussed.     

Nonlinear wave propagation along periodic-loaded transmission line

A high frequency transmission line periodically loaded with varactor diodes is presented to study nonlinear wave propagation. The nonlinearity and dispersion characteristics are experimentally and theoretically analysed. Experimental results on shock wave formation and harmonic frequency generation are found to be in good agreement with theoretical predictions. The influence of dispersion on second-harmonic frequency generation is particularly discussed yielding large values of the efficiency.     

高功率高效率三次谐波二维数值模拟研究

 对大口径、高功率情况下的三次谐波转换,选择KDP晶体TypeⅠ/TypeⅡ匹配角度失谐的三倍频方案,并从非线性三波耦合方程组出发,采用离散傅立叶变换和四阶龙格库塔(R-K)积分方法,编制了二维模拟三次谐波转换的计算程序。计算了各种晶体厚度和不同失谐角条件下,二、三次谐波的转换特性和效率。并对入射基频光为六阶超高斯分布时的三倍频器进行了优化,三倍频转换效率达80%,此时具有较大的动态范围和较高的三倍频转换效率。     

Effect of an elastic wall on the dynamics of an encapsulated microbubble: A simulation study

The purpose of the present simulation study is to reveal how confining surfaces with different mechanical properties affect the acoustic response of a contrast agent microbubble. To this end, numerical simulations are carried out for three types of walls: a plastic (OptiCell) wall, an aluminium wall, and a biological tissue. For each wall, the behaviour of contrast microbubbles of three sizes is investigated. The spectral characteristics of the scattered pressure produced by the microbubbles are compared for two cases: the bubble oscillates far away from the wall and the same bubble oscillates in the immediate vicinity of the wall. The results of the simulations allow one to make the following main conclusions. The effect of the OptiCell wall on the acoustic bubble response is stronger than that of the aluminium and tissue walls. Changes in the bubble response near the wall are stronger when bubbles are excited above their resonance frequency. Considering changes in the fundamental and the 2nd harmonic with respect to the peak values of these components at different bubble radii, it is found that the changes are stronger for smaller bubbles and that the changes in the 2nd harmonic are stronger than those in the fundamental. These results allow one to gain an insight into conditions under which the effect of an elastic wall on the acoustic response of a contrast agent microbubble is easier to be detected.     

非线性声波方程二次谐波高阶近似解及实验验证

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差。利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解。在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对幅度随传播距离和输入信号均呈现先增加后减小的趋势,高阶近似解即使在传播距离不太短和输入信号不太弱的宽泛条件下也能更好地吻合实验结果。得到的高阶近似解不仅弥补了现有二次谐波幅度理论公式的不足,进行非线性声学实验时还可扩展测量范围,充分利用实验数据也能提高非线性系数的测量精度。      

Optimization of acoustic scattering from dual-frequency driven microbubbles at the difference frequency

The second harmonic radiation of acoustically driven bubbles is a useful discriminant for their presence in clinical ultrasound applications. It is useful because the scatter from a bubble at a frequency different from the driving can have a contrast-to-tissue ratio better than at the drive frequency. In this work a technique is developed to optimize the scattering from a microbubble at a frequency different from the driving. This is accomplished by adjusting the relative phase and amplitudes of the components of a dual-frequency incident ultrasound wave form. The investigation is focused primarily on the example of dual-mode driving at frequencies of 1 MHz and 3 MHz, with the scattering optimized at 2 MHz. Bubble radii of primary interest are 0.5 to 2 microm and driving amplitudes to 0.5 atm. Bubbles in this size range are sensitive to modulation of driving. It is shown that an optimal forcing scheme can increase the target response eightfold or more. This suggests new applications in imaging and in bubble detection.