数值计算法实现同振式矢量水听器驻波管中灵敏度修正EI北大核心CSCD

大尺寸的矢量水听器在驻波管中进行灵敏度校准时会引起校准声场分布不均匀,使得灵敏度校准结果出现较大误差。在无法获取校准声场准确的解析表达式情况下,以实际驻波管校准装置为原型,建立声场仿真分析模型,通过数值计算得到矢量水听器和参考标量水听器的声压和加速度量值,并结合灵敏度校准公式和参考值导出灵敏度修正因子。利用修正因子,对不同尺寸矢量水听器灵敏度测试结果修正以后,与参考值最大绝对误差减小到2.0 dB下。结果表明数值计算法能有效减小非均匀声场引起的大尺寸矢量水听器灵敏度校准误差,扩展了矢量水听器驻波管校准装置的测试对象。     

驻波声场中悬浮临界密度及稳定性研究

本文以声场中物体为研究对象,理论上得到行波和驻波场中的声辐射压力方程.在驻波声场中引入临界悬浮密度概念,可作为物体能否在非线性声场中悬浮的判据,同时给出谐振腔移动速度的最大范围.更进一步,以实验参数作为数值计算的输入来指导实验,并结合实验结果讨论了驻波声场中样品密度和大小、发射面和反射面形状以及两者之间的距离、反射面的尺寸等因素对物体悬浮稳定性的影响,发现当物体尺寸和密度确定时,调控好谐振腔的长度,增加波腹处的声压是提升声悬浮稳定性的有效手段.     

热声发动机的CFD数值模拟

采用计算流体动力学(CFD)方法,对高频驻波热声发动机和热声斯特林发动机实验系统分别进行了二维和三维数值模拟。计算模型具有与实验系统相同的几何结构、尺寸和运行工况。对计算模型的有效性进行了研究,表明实现有限换热条件的板叠实物模型适合驻波热声发动机的模拟,而实现局域热平衡的多孔介质模型适合热声斯特林发动机的模拟。计算结果成功观测到了非线性的自激振荡演化过程,捕捉到了两种发动机的不同非线性现象。计算结果分别给出了两种热声发动机内部的声场分布特性和复杂流场。计算结果与实验结果的对比验证了CFD方法对高频驻波热声发动机和热声斯特林发动机模拟的有效性。     

格子玻尔兹曼方法对不同张角聚焦声束的建模

高强度聚焦超声(HIFU)是一种新型的无创治疗肿瘤新技术,其中换能器声场数值计算能够为HIFU治疗提供重要的依据。传统非线性KZK和SBE模型广泛应用于换能器声场数值计算,但依然存在某些不足。我们采用一种介观尺度的新型流体力学方法,即格子Boltzmann方法(LBM),基于2维9离散速度(D2Q9)格子构建了轴对称多弛豫参数LBM模型,并通过调节弛豫参数分析其对模型的影响;利用该模型对两个具有不同张角的球面聚焦换能器的声场进行数值模拟,并与KZK和SBE模型的计算结果进行比较。结果表明LBM模型能够很好地描述超声波的激发和传播机制,从流体力学的角度描述聚焦声场的分布,具有清晰的物理意义,且计算过程不受换能器张角的限制,在换能器声场的理论分析和模拟计算及其在HIFU治疗中的应用有着积极的意义。      

在直角坐标系用NGFM模拟复杂计算域水下高压气泡膨胀问题

采用NGFM(New version of Ghost Fluid Method)处理复杂计算域的固壁边界,用RGFM(Real Ghost Fluid Method)求解气-水界面附近网格节点的状态参数,从而在直角坐标系下对复杂计算域的水下高压气泡膨胀问题进行数值模拟。流场控制方程选用Euler方程,用五阶WENO格式离散空间导数项,二阶Runge-Kutta法离散时间导数项;气-水界面追踪使用Level Set方法,对Level Set方程,用五阶HJ-WENO(Hamilton-Jacobi WENO)和三阶Runge-Kutta法求解。将计算结果与任意坐标系下的结果进行对比,验证了NGFM在笛卡尔网格中处理复杂形状固壁边界的可行性。得到了水下流场压力等值线图、高压气泡的演变过程以及特定点处的压力-时间曲线。计算结果表明,高压气泡在固壁反射激波的作用下,膨胀过程受到抑制;强激波在固壁的反射会导致固壁附近出现大范围的空化流动。     

气泡引起的皇冠型水冢实验与数值研究

气泡在近自由液面运动过程中与自由液面发生强烈的耦合作用, 自由液面会出现极其复杂的物理现象, 形成的水冢类型繁多. 本文将针对皇冠型水冢, 基于势流理论, 建立气泡与自由液面耦合作用的数值模型, 并在气泡完成射流向下运动时, 忽略气泡对自由液面的影响, 继续模拟自由液面的运动过程. 同时利用高速摄影对近自由面的电火花气泡进行实验研究, 数值结果与实验符合良好, 相对误差在10%以内. 通过数值计算, 发现了围裙卷缩和主峰珠化等特殊的物理现象, 研究了气泡初始条件与韦伯数对皇冠型水冢动态特性的影响, 旨为皇冠型水冢的研究提供参考. 关键词： 气泡 自由液面 皇冠型水冢 韦伯数     

超声场中气泡的耦合运动

超声场中气泡除径向振动外,还可能会平动并且相互影响.本文以考虑邻近气泡次级声辐射影响后的球形气泡径向振动模型为基础,结合气泡在声场中受到的力,利用数值方法研究了平面波声场中不同尺寸的两气泡径向振动和平动规律.发现尺寸较大气泡的径向振动具有一定的本征性特征且具有较大平动位移.利用高速摄影系统定性地观察了气泡运动和泡群分布...     

用阴影法观察水声模型的声场

我们建立了一套用于显示模拟海底反射声场的光学系统。实验水声模型是由水-玻璃粉-玻璃-铝组成的三层模拟海底结构,声源是一个管状换能器,其轴线平行于分界面。实验用这套系统记录了从水下三个界面反射的六组波的声场图像。并且分辨出了它们各自的传播路径。从声场图像也可以得出水-玻璃粉界面、玻璃粉-玻璃界面的面波声速,玻璃粉和玻璃中的体波声速以及它们的厚度,结果与预先测出的参数符合。     

含气泡液体中气泡振动的研究

研究了含气泡液体中单个气泡在驱动声场一定情况下的振动过程. 让每次驱动声场作用的时间特别短, 使气泡半径发生微小变化后再将其变化反馈到气泡群对驱动声场的散射作用中去, 从而可以得到某单个气泡周围受气泡散射影响后的声场, 接着再让气泡在该声场作用下做短时振动, 如此反复. 通过这样的方法, 研究了液体中单个气泡的振动情况并对其半径变化进行了数值模拟, 结果发现, 在液体中含有大量气泡的情况下, 某单个气泡的振动过程明显区别于液体中只有一个气泡的情况. 由于大量气泡和驱动声场的相互作用, 使气泡半径的变化存在多种不同的振动情况, 在不同的气泡大小和含量的情况下, 半径变化过程分别表现为: 在平衡位置附近振荡的过程; 周期性的空化过程; 一次空化过程后保持某一大小振荡的过程; 增长后维持某一大小振荡的过程等. 所以, 对于含气泡液体中气泡振动的研究, 在驱动声场一定的情况下, 必须考虑气泡含量的因素. 关键词： 含气泡液体 超声空化 散射 数值模拟     

球形腔聚焦换能器的非线性声场建模

球形腔聚焦换能器是一种特殊形式的聚焦换能器。为理论证实球形腔聚焦换能器能突破传统超声聚焦在聚焦精度和聚焦增益上的限制,采用Westervelt非线性方程并结合时域有限差分法,建立了球形腔聚焦换能器的非线性声场的数值模型。数值计算了直径为120 mm的0.6 MHz球形腔聚焦换能器的非线性声场,并与传统球壳形聚焦换能器进行了对比。当激励声压为100 kPa时,球形腔聚焦换能器与同尺寸壳形聚焦换能器相比,焦点正声压增益提高约8.5倍,且焦域精度更高,-6 dB聚焦区域在z方向减小约20倍,达到次波长尺度。研究表明球形腔聚焦换能器在高强度聚焦超声精细治疗上具有潜在的应用前景。      

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

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

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.     

两种气泡混合的声空化

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

Characterization of an acoustic cavitation bubble structure at 230 kHz

A generic bubble structure in a 230 kHz ultrasonic field is observed in a partly developed standing wave field in water. It is characterized by high-speed imaging, sonoluminescence recordings, and surface cleaning tests. The structure has two distinct bubble populations. Bigger bubbles (much larger than linear resonance size) group on rings in planes parallel to the transducer surface, apparently in locations of driving pressure minima. They slowly rise in a jittering, but synchronous way, and they can have smaller satellite bubbles, thus resembling the arrays of bubbles observed by Miller [D. Miller, Stable arrays of resonant bubbles in a 1-MHz standing-wave acoustic field, J. Acoust. Soc. Am. 62 (1977) 12]. Smaller bubbles (below and near linear resonance size) show a fast "streamer" motion perpendicular to and away from the transducer surface. While the bigger bubbles do not emit light, the smaller bubbles in the streamers show sonoluminescence when they pass the planes of high driving pressure. Both bubble populations exhibit cleaning potential with respect to micro-particles attached to a glass substrate. The respective mechanisms of particle removal, though, might be different.     

Stable multibubble sonoluminescence bubble patterns

Multibubble standing wave patterns can be generated from a flat piezoceramic transducer element radiating into water. By adding a second transducer positioned at 90 degrees from the transducer generating the standing wave, a 3-dimensional volume of stable single bubbles can be established. Further, the addition of the second transducer stabilizes the bubble pattern so that individual bubbles may be studied. The size of the bubbles and the separation of the standing waves depend on the frequency of operation. Two transducers, operating at frequencies above 500 kHz, provided the most graphic results for the configuration used in this study. At these frequencies stable bubbles exhibit a bright sonoluminescence pattern. Whereas stable SBSL is well-known, stable MBSL has not been previously reported. This paper includes discussions of the acoustic responses, standing wave patterns, and pictorial results of the separation of individual bubble sonoluminescence in a multibubble sonoluminescence environment.     

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.     

Physical investigation of acoustic waves induced by the oscillation and collapse of the single bubble

The objective of this paper is to apply both experimental and numerical methods to investigate acoustic waves induced by the oscillation and collapse of a single bubble. In the experiments, the schlieren technique is used to capture the temporal evolution of the bubble shapes, and the corresponding acoustic waves. The results are presented for the single bubble generated by a low-voltage bubble generator in the free field of water. During the numerical simulations, a three-dimensional (3D) weakly compressible model is introduced to investigate the single bubble dynamics, including the generation and propagation of acoustic waves. The results show that (1) Compression wave, rarefaction wave and shock wave are generated during expansion stage, collapse stage and rebound stage of the bubble respectively. (2) Compression waves are induced by the rapid expansion of the bubble and eventually steepen into one shock wave propagating outward in the liquid, then another strong shock wave is emitted at the final collapse stage. The velocity and pressure of the liquid field increases after the shock wave. (3) Rarefaction waves are generated during the collapse stage due to the contraction of the bubble. The rarefaction wave reduces the liquid pressure and its spatial distribution is dispersive. The pressure of these acoustic waves and their effect on the liquid velocity attenuate with the increase of propagation distance.     

Acoustic characteristics of bubble bursting at the surface of a high-viscosity liquid

An acoustic pressure model of bubble bursting is proposed.An experiment studying the acoustic characteristics of the bursting bubble at the surface of a high-viscosity liquid is reported.It is found that the sudden bursting of a bubble at the high-viscosity liquid surface generates N-shape wave at first,then it transforms into a jet wave.The fundamental frequency of the acoustic signal caused by the bursting bubble decreases linearly as the bubble size increases.The results of the investigation can be used to understand the acoustic characteristics of bubble bursting.     

Nonlinear dynamics of a cavitation bubble pair near a rigid boundary in a standing ultrasonic wave field

The dynamics of a micrometer-sized bubble pair in water near a rigid boundary under standing ultrasonic wave excitation is investigated in this study. The viscous effect in the boundary layer at the air-water interface is considered following the viscous correction model. The evolution of the bubble surface at the collapsing stage of the bubble pair is presented for different parameter sets. The field pressure near the rigid boundary, which is induced by the oscillating bubble pair, and the high-speed water jet at the collapse stage, form the main focus of the analysis. This reveals that a horizontal configuration of the bubble pair retards the strength of the bubble jet towards the boundary, whilst a vertical configuration, especially with differently-sized bubbles, can enhance the bubble collapse. This study may help to understand the interaction of multiple bubbles in an acoustic field and its application to surface cleaning.     

Collective bubble dynamics near a surface in a weak acoustic standing wave field

The transport of bubbles to a neighboring surface is very important in surface chemistry, bioengineering, and ultrasonic cleaning, etc. This paper proposes a multi-bubble transport method by using an acoustic standing wave field and establishes a model that explains the multi-bubble translation by expressing the balance between Bjerknes forces and hydrodynamic forces on a bubble in a liquid medium. Results indicated that the influence of primary Bjerknes force, secondary Bjerknes force, and buoyancy force on the bubble translation depends on the position of the target bubble in the acoustic field. Moreover, it was found that increasing the size of a bubble or pressure amplitude can accelerate the bubble motion and enhance the bubble-bubble interaction. The secondary Bjerknes force between two bubbles can switch from an attractive one when they oscillate in phase to a repulsive one when the bubble oscillations are out of phase. These findings provide an insight into the multi-bubble translation near a surface and can be applied to future bubble motion control studies, especially in drug delivery, sonoporation, and ultrasonic cleaning.