水下等离子体声源的冲击波负压特性

基于修正的Rayleigh气泡脉动方程对水下等离子体声源放电产生的 强声冲击波的传播过程进行了分析; 利用Euler方程作为控制方程组, 建立了水下等离子体声源的聚束声场模型, 通过仿真计算获得的传播云图对冲击波负压的形成机理进行了直观的理论分析. 结果表明: 经过聚能反射罩反射汇聚得到的聚束波在反射稀疏波和水的惯性作用下, 聚束波周围水域产生了拉伸, 形成负压区, 如果拉伸力大于水的抗拉上限, 就会使得水中形成不连续现象, 即出现空化气泡; 此外聚能罩边缘处产生的衍射波进一步加剧了负压的产生, 边缘衍射波最终与拉伸波叠加, 使冲击波负压达到最大值; 通过对比仿真波形和实验波形, 从而验证和进一步揭示了冲击波负压的形成原因. 研究结果对认识水下冲击波的传播规律和进一步改进等离子体声源的设计具有指导意义. 关键词： 等离子体声源 冲击波负压 聚束声场模型 气泡     

液体薄层中环链状空化泡云结构稳定性分析

为分析超声空化的薄层液体中稳定的环状气泡链结构,本文考虑气泡间次级声辐射影响,得到了表征气泡间相互作用的气泡基本动力学方程以及次Bjerknes力的表达式,数值分析了气泡平衡半径、声波频率和声压对纯液体区可能出现的单气泡所受的次Bjerknes力,发现环形泡链能够吸引液体区内的新生的半径小于2μm的气泡,这可能是一定条件下环形气泡链能够稳定存在的原因.随着驱动声波压力增加,气泡数密度增加,气泡间的耦合作用增强,液体区内的环形泡链结构可能被液体区内出现的大气泡或者气泡团破坏,进而导致环形结构演变成柱状、雾状乃至整个液体区均充满空化泡的情况发生.通过高速摄影机观察了强声场作用下换能器辐射面外侧液体薄层内空化初生至形成空化云团簇的整个过程,在空化云团簇中发现了局部同步崩溃并形成类纯液体薄层的现象,该液体薄层边界随时间振荡持续约4个声周期后被空化云团簇吞没,局部类纯液体区出现的位置具有随机性.实验观察结果和理论预测具有很好的一致性.     

声场中气泡运动的混沌特性

以水为工作介质,考虑了液体的可压缩性,研究了声场中气泡的运动特性,模拟了声波频率、声压幅值、气泡初始半径以及液体的表面张力和黏滞系数的变化对气泡运动状态的影响. 分析了空化处理效果与气泡运动状态之间关系. 结果表明:气泡运动处于混沌状态,是提高声空化降解有机污染物能力的最重要因素. 关键词： 声空化 混沌 相图 功率谱图     

深海大深度声传播特性及直达声区水下声源距离估计

深海大深度声传播特性对在深海近海底进行水声目标探测和定位具有重要意义。利用一次南海中南部深海不完全声道中的脉冲声传播实验数据,分析了海底附近大深度声传播损失及脉冲多途传播特性,并根据直达波和海底-海面反射波的时延差与收发距离的关系,提出一种利用深海直达声区脉冲多途到达时间进行水下声源距离估计的方法。结果表明:当接收器深度位于南海深海海底附近而声源深度较浅时,直达声区水平宽度可达30 km,传播损失相对影区来说较小,有利于水下声源探测;直达声区的直达波与海底-海面反射波的到达时延差随着收发距离的增大单调减小,可被用于水下声源距离估计。得到水下声源的距离估计结果与实验GPS测量结果较为一致,距离估计均方误差为0.28 km。      

低频超声空化场中柱状泡群内气泡的声响应

本文在气泡群振动模型的基础上,考虑气泡间耦合振动的影响,得到了均匀柱状泡群内振动气泡的动力学方程,以此为基础分析了低频超声空化场中柱形气泡聚集区内气泡的非线性声响应特征.气泡间的耦合振动增加了系统对每个气泡的约束,降低了气泡的自然频率,增强了气泡的非线性声响应.随着气泡数密度的增加,气泡的自然共振频率降低,受迫振动气泡受到的抑制增强.数值分析结果表明:1)驱动声波频率越低,气泡的初始半径越小,气泡数密度变化对气泡最大半径变化幅度的影响越大;2)气泡振动幅值响应存在不稳定区,不稳定区域分布与气泡初始半径、驱动声波压力幅值、驱动声波频率等因素有关.在低频超声波作用下,对初始半径处在1—10μm之间的空化气泡而言,气泡初始半径越小,气泡最大半径不稳定区分布范围越大,表明小气泡具有更强的非线性特征.因此,气泡初始半径越小,声环境变化对空化泡声响应稳定性影响越显著.     

深海海底斜坡环境下的声传播

海底地形变化对声传播具有很大影响,在南海深海区域海底斜坡环境下进行了一次声传播实验,实验显示倾斜海底环境下声传播损失出现了一些不同于平坦海底环境下的现象,分析并解释了海底地形变化对产生声传播差异的原因.结果表明,海底斜坡对声波的反射增强作用可使斜坡上方的声传播损失减少约5 d B.当声波第一次入射到达的海底位置有较小幅度的山丘(凸起高度小于1/10海深)时,海底小山丘即可对声波有反射遮挡作用,导致在其反射区特定传播距离和深度上出现倒三角声影区,比平坦海底环境下相同影区位置处的传播损失增大约8 d B,影响深度可达海面以下1500 m.而海底斜坡对声波的反射阻挡作用使得从海面反射及水体向下折射的会聚区结构消失,只剩下从水体向上折射的会聚结构.因此,海底地形对深海声传播影响较大,在水下目标探测和性能评估等应用中应予以重视.     

不同水底介质对有限域中装药沉底爆炸特性的影响

通过开展装药沉底爆炸原理性实验研究,得到了不同水底条件对装药沉底爆炸气泡运动和冲击波压力的影响规律:沉底装药水下爆炸气泡通常呈半球形,依附在水底并同时急剧膨胀达到最大半径,随后气泡做收缩运动并连带水底介质颗粒迅速上浮,同时形状发生显著变化,如在石底和泥底时气泡在水底射流的作用下呈蘑菇状上浮,砂底时气泡则呈现出柱状上浮、未出现明显脉动运动即已坍塌破裂;沉底爆炸与水底介质发生强烈耦合作用而形成不同程度的反射波,通常水下石底的反射冲击波较泥底和砂底更强烈,且一般迭加后的冲击波峰值压力高于入射波阵面压力;不同水底介质条件下,装药与水底发生的强烈耦合作用均对水底造成强冲击破坏。     

硫酸中多气泡声致发光光谱

非线性声波方程与气泡脉动方程联立, 可以描述声空化云中的声场以及任何一个气泡的脉动过程,为数值计算空化场问题提供了理论框架.计算的声压分布变化可以用来计算单气泡动力学,了解任何位置处气泡发光过程以及气泡内气体温度和压强变化等. 对浓硫酸中氙气泡空化云的计算定性符合实验观测, 只有钠原子线谱的计算结果相比实验观测有些出入.     

含混合气泡液体中声波共振传播的抑制效应

当声波在含气泡的液体中传播时会出现共振传播现象,即在气泡的共振频率附近声衰减和声速会显著地增大,这是声空化领域的一个重要现象.以往的研究一般假设液体中只存在单一种类的气泡,因此忽略了声波共振传播的某些重要信息.本文研究了含混合气泡液体中声波的共振传播,混合气泡是指液体中包含多种静态半径不同的气泡.结果显示:在这种系统中存在声波共振传播的抑制效应,即与含单一种类气泡的系统相比,在含混合气泡的系统中声波的共振衰减和共振声速会明显变小.对于两种气泡混合、多种气泡混合以及气泡满足某种连续分布的系统,研究了抑制效应的本质和主要特征,此外还探究了黏性和空化率等对抑制效应的影响.本文的研究结果是对该领域现有知识的必要补充.     

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

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

Acoustic characteristics of pulse detonation engine with ellipsoidal reflector

Acoustic characteristics of a pulse detonation engine(PDE) with and without an ellipsoidal reflector are numerically and experimentally investigated. A two-dimensional(2 D) non-splitting unstructured triangular mesh Euler solver based on the space-time conservation element and solution element(CE/SE) method is employed to simulate the flow field of a PDE.The numerical results clearly demonstrate the external flow field of the PDE. The effect of an ellipsoidal reflector on the flow field characteristic near the PDE exit is investigated. The formation process of reflected shock wave and reflected jet shock are reported in detail. An acoustic measurement system is established for the PDE acoustic testing. The experimental results show that the ellipsoidal reflector changes the sound waveform and directivity of PDE sound. The reflected shock wave and reflected jet shock result in two more positive pressure peaks in the sound waveform. The ellipsoidal reflector changes the directivity of PDE sound from 20 to 0. It is found that the peak sound pressure level(PSPL) and overall sound pressure level(OASPL) each obtain an increment when the PDE is installed with a reflector. The maximum relative increase ratio of PSPL and OASPL are obtained at the focus point F2, whose values are 6.1% and 6.84% respectively. The results of the duration of the PDE sound indicate that the reflecting and focusing wave generated by the reflector result in the increment of A duration and B duration before and near focus point F2. Results show that the ellipsoidal reflector has a great influence on the acoustic characteristic of PDE sound. The research is helpful for understanding the influence of an ellipsoidal reflector on the formation and propagation process of PDE sound.     

Modeling of an electrohydraulic lithotripter with the KZK equation.

The acoustic pressure field of an electrohydraulic extracorporeal shock wave lithotripter is modeled with a nonlinear parabolic wave equation (the KZK equation). The model accounts for diffraction, nonlinearity, and thermoviscous absorption. A numerical algorithm for solving the KZK equation in the time domain is used to model sound propagation from the mouth of the ellipsoidal reflector of the lithotripter. Propagation within the reflector is modeled with geometrical acoustics. It is shown that nonlinear distortion within the ellipsoidal reflector can play an important role for certain parameters. Calculated waveforms are compared with waveforms measured in a clinical lithotripter and good agreement is found. It is shown that the spatial location of the maximum negative pressure occurs pre-focally which suggests that the strongest cavitation activity will also be in front of the focus. Propagation of shock waves from a lithotripter with a pressure release reflector is considered and because of nonlinear propagation the focal waveform is not the inverse of the rigid reflector. Results from propagation through tissue are presented; waveforms are similar to those predicted in water except that the higher absorption in the tissue decreases the peak amplitude and lengthens the rise time of the shock.     

Comparison of electrohydraulic lithotripters with rigid and pressure-release ellipsoidal reflectors. II. Cavitation fields.

Dramatically different cavitation was produced by two separate acoustic pulses that had different shapes but similar duration, frequency content, and peak positive and negative pressure. Both pulses were produced by a Dornier HM-3 style lithotripter: one pulse when the ellipsoidal reflector was rigid, the other when the reflector was pressure release. The cavitation, or bubble action, generated by the conventional rigid-reflector pulse was nearly 50 times longer lived and 3-13 times stronger than that produced by the pressure-release-reflector pulse. Cavitation durations measured by passive acoustic detection and high-speed video agreed with calculations based on the Gilmore equation. Cavitation intensity, or destructive potential, was judged (1) experimentally by the size of pits in aluminum foil detectors and (2) numerically by the calculated amplitude of the shock wave emitted by a collapsing bubble. The results indicate that the trailing positive spike in the pressure-release-reflector waveform stifles bubble growth and mitigates the collapse, whereas the trough after the positive spike in the rigid-reflector waveform triggers inertially driven growth and collapse. The two reflectors therefore provide a tool to compare effects in weakly and strongly cavitating fields and thereby help assess cavitation's role in lithotripsy.     

The effect of reflector geometry on the acoustic field and bubble dynamics produced by an electrohydraulic shock wave lithotripter

A theoretical model for the propagation of shock wave from an axisymmetric reflector was developed by modifying the initial conditions for the conventional solution of a nonlinear parabolic wave equation (i.e., the Khokhlov-Zabolotskaya-Kuznestsov equation). The ellipsoidal reflector of an HM-3 lithotripter is modeled equivalently as a self-focusing spherically distributed pressure source. The pressure wave form generated by the spark discharge of the HM-3 electrode was measured by a fiber optic probe hydrophone and used as source conditions in the numerical calculation. The simulated pressure wave forms, accounting for the effects of diffraction, nonlinearity, and thermoviscous absorption in wave propagation and focusing, were compared with the measured results and a reasonably good agreement was found. Furthermore, the primary characteristics in the pressure wave forms produced by different reflector geometries, such as that produced by a reflector insert, can also be predicted by this model. It is interesting to note that when the interpulse delay time calculated by linear geometric model is less than about 1.5 micros, two pulses from the reflector insert and the uncovered bottom of the original HM-3 reflector will merge together. Coupling the simulated pressure wave form with the Gilmore model was carried out to evaluate the effect of reflector geometry on resultant bubble dynamics in a lithotripter field. Altogether, the equivalent reflector model was found to provide a useful tool for the prediction of pressure wave form generated in a lithotripter field. This model may be used to guide the design optimization of reflector geometries for improving the performance and safety of clinical lithotripters.     

Focusing of shock waves induced by optical breakdown in water

The focusing of laser-generated shock waves by a truncated ellipsoidal reflector was experimentally and numerically investigated. Pressure waveform and distribution around the first (F(1)) and second foci (F(2)) of the ellipsoidal reflector were measured. A neodymium doped yttrium aluminum garnet laser of 1046 nm wavelength and 5 ns pulse duration was used to create an optical breakdown at F(1), which generates a spherically diverging shock wave with a peak pressure of 2.1-5.9 MPa at 1.1 mm stand-off distance and a pulse width at half maximum of 36-65 ns. Upon reflection, a converging shock wave is produced which, upon arriving at F(2), has a leading compressive wave with a peak pressure of 26 MPa and a zero-crossing pulse duration of 0.1 mus, followed by a trailing tensile wave of -3.3 MPa peak pressure and 0.2 mus pulse duration. The -6 dB beam size of the focused shock wave field is 1.6 x 0.2 mm(2) along and transverse to the shock wave propagation direction. Formation of elongated plasmas at high laser energy levels limits the increase in the peak pressure at F(2). General features in the waveform profile of the converging shock wave are in qualitative agreement with numerical simulations based on the Hamilton model.     

Shock wave-inertial microbubble interaction: methodology, physical characterization, and bioeffect study

A method of generating in situ shock wave-inertial microbubble interaction by a modified electrohydraulic shock wave lithotripter is proposed and tested in vitro. An annular brass ellipsoidal reflector (thickness = 28 mm) that can be mounted on the aperture rim of a Dornier XL-1 lithotripter was designed and fabricated. This ring reflector shares the same foci with the XL-1 reflector, but is 15 mm short in major axis. Thus, a small portion of the spherical shock wave, generated by a spark discharge at the first focus (F1) of the reflector, is reflected and diffracted by the ring reflector, producing a weak shock wave approximately 8.5 microseconds in front of the lithotripter pulse. Based on the configuration of the ring reflector (different combinations of six identical segments), the peak negative pressure of the preceding weak shock wave at the second focus (F2) can be adjusted from -0.96 to -1.91 MPa, at an output voltage of 25 kV. The preceding shock wave induces inertial microbubbles, most of which expand to a maximum size of 100-200 microns, with a few expanding up to 400 microns before being collapsed in situ by the ensuing lithotripter pulse. Physical characterizations utilizing polyvinylidene difluoride (PVDF) membrane hydrophone, high-speed shadowgraph imaging, and passive cavitation detection have shown strong secondary shock wave emission immediately following the propagating lithotripter shock front, and microjet formation along the wave propagation direction. Using the modified reflector, injury to mouse lymphoid cells is significantly increased at high exposure (up to 50% with shock number > 100). With optimal pulse combination, the maximum efficiency of shock wave-induced membrane permeabilization can be enhanced substantially (up to 91%), achieved at a low exposure of 50 shocks. These results suggest that shock wave-inertial microbubble interaction may be used selectively to either enhance the efficiency of shock wave-mediated macromolecule delivery at low exposure or tissue destruction at high exposure.     

声驻波场中空化泡的动力学特性

以水为工作介质, 考虑了液体的可压缩性, 研究了驻波声场中空化泡的运动特性, 模拟了驻波场中各位置处空化泡的运动状态以及相关参数对各位置处空化泡在主Bjerknes力作用下运动方向的影响. 结果表明: 驻波声场中, 空化泡的运动状态分为三个区域, 即在声压波腹附近空化泡做稳态空化, 在偏离波腹处空化泡做瞬态空化, 在声压波节附近, 空化泡在主Bjerknes 力作用下, 一直向声压波节处移动, 显示不发生空化现象; 驻波场中声压幅值增加有利于空化的发生, 但声压幅值增加到一定上限时, 压力波腹区域将排斥空化泡, 并驱赶空化泡向压力波节移动, 不利于空化现象的发生; 当声频率小于初始空化泡的共振频率时, 声频率越高, 由于主Bjerknes 力的作用将有更多的空化泡向声压波节移动, 不利于空化的发生, 尤其是驻波场液面的高度不应是声波波长的1/4; 当声频率一定时, 空化泡初始半径越大越有利于空化现象的发生, 但当空化泡的初始半径超过声频率的共振半径时, 由于主Bjerknes力的作用将有更多的空化泡向声压波节移动, 不利于空化的发生.     

声波在含气泡液体中的线性传播

为了探讨含气泡液体对声波传播的影响, 研究了声波在含气泡液体中的线性传播. 在建立含气泡液体的声学模型时引入气泡含量的影响,建立气泡模型时引用 Keller的气泡振动模型并同时考虑气泡间的声相互作用,得到了经过修正的气泡振动方程. 通过对含气泡液体的声传播方程和气泡振动方程联立并线性化求解,在满足 (ω R₀)/c << 1 的前提下,得到了描述含气泡液体对声波传播的衰减系数和传播速度. 通过数值分析发现,在驱动声场频率一定的情况下,气泡含量的增加及气泡的变小均会导致衰减系数增加和声速减小;气泡的体积分数和大小一定时, 驱动声场频率在远小于气泡谐振频率的情况下,声速会随驱动频率的增加而减小; 气泡间的声相互作用对声波传播速度及含气泡液体衰减系数的影响不明显.最终认为气泡的大小、 数量和驱动声场频率是影响声波在含气泡液体中线性传播的主要因素. 关键词： 含气泡液体 线性声波 声衰减系数 声速     

Comparison of electrohydraulic lithotripters with rigid and pressure-release ellipsoidal reflectors. I. Acoustic fields.

The most common lithotripter, a Dornier HM-3, utilizes an underwater spark to generate an acoustic pulse and a rigid ellipsoidal reflector to focus the pulse on the kidney stone to be comminuted. The pulse measured in water with a PVDF membrane hydrophone at the external focus of the ellipsoid was a 1-microsecond positive-pressure spike followed by a 3-microsecond negative-pressure trough. When we replaced the rigid reflector in our experimental lithotripter with a pressure-release reflector, the pulse was a 1.6-microsecond trough followed by a 0.6-microsecond positive spike. The waveforms are nearly time inverses (i.e., their spikes and troughs are reversed). The frequency spectra, the maximum peak positive pressures P+ (42 MPa, rigid and 43 MPa, pressure-release), and the maximum peak negative pressures P- (-12 MPa and -14 MPa) are comparable. The maximum P- occurred 20 mm closer to the reflector than did the maximum P+, for both reflectors. However, the spatial maxima of the peak pressures (P+ and P-) produced by the pressure-release reflector were located 20 mm nearer to the reflector than those produced by the rigid reflector. Qualitative explanation of the waveforms and the location of pressure maxima as well as comparison to previous theoretical and experimental results is given. The alternate waveform produced by the pressure-release reflector may be a tool in determining the role of cavitation in lithotripsy because cavitation is highly sensitive to waveform.