驻波声场中一对气泡的振动特征及其相互作用力,声学论文.docx
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《驻波声场中一对气泡的振动特征及其相互作用力,声学论文.docx》由会员分享,可在线阅读,更多相关《驻波声场中一对气泡的振动特征及其相互作用力,声学论文.docx(18页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、驻波声场中一对气泡的振动特征及其相互作用力,声学论文摘 要: 该研究从驻波声场特征出发,研究一对气泡在驻波声场中的振动特征及气泡之间的互相作用力。研究结果表示清楚:驻波声场中波腹位置附近的气泡具有强烈的能量转换能力,能够在小振幅声压下做非线性振动,发生剧烈的空化现象,波节位置附近的气泡只做微小的线性振荡;驻波声场中处于波腹位置处的两个气泡之间的互相作用最大,波节位置处两个气泡之间的互相作用力趋于0;驻波声场中两个气泡之间互相作用力的大小和方向取决于两个气泡在驻波声场中的位置、驱动声压振幅、驱动频率和两个气泡的初始半径等因素。 本文关键词语 : 驻波声场;气泡;能量转换;互相作用力; Abstr
2、act: The characteristics of bubbles and the interaction between two bubbles were studied based on the characteristics of the acoustic standing wave field. The results show that in the standing wave sound field, the bubbles which are near the antinode position have a strong energy conversion ability.
3、 A bubble can vibrate nonlinearly under a small amplitude sound pressure and the cavitation can happen in the acoustic standing wave field. The interaction between two bubbles is the largest at the antinode position, and the interaction between two bubbles at the node position tends to 0. The magnit
4、ude and direction of the interaction between two bubbles depend on the amplitude of the driving sound pressure, the driving frequency and the initial radius of the two bubbles in the acoustic standing wave field. Keyword: acoustic standing wave field; bubble; energy conversion; interaction between b
5、ubbles; 超声空化试验发现,声场中的气泡并不是均匀存在,而是构成一些稳定的气泡构造1,这些气泡构造甚至能够存在数秒,关于声场中气泡怎样构成这样的稳定构造,当前理论上不清楚,但是,气泡之间的互相作用力是其构成的主要原因2,3,4,声场中气泡之间的互相作用力主要是指次Bjerknes力,是由两个相邻气泡之中的某一个气泡产生的“第二声场 使两个气泡之间产生的作用力。自次Bjerknes力发现之后,Crum5、Oguz等6和Mettin等7分别进行了有意义的研究工作,得到了气泡线性振动和非线性振动下,气泡之间次Bjerknes力的相关规律,为声场中气泡之间的互相作用提供了理论基础8,9,10,但
6、是他们的研究都是建立在平面声场中11,12,而超声反响器中,由于反响器尺寸一定,所以不可避免的在超声反响器中会构成驻波声场13,驻波声场中气泡之间的互相作用力具有什么样的特征,对空化气泡的振动特性及气泡构造产生什么样的影响,这些问题是研究超声反响器中气泡空化效应的重要因素,也是研究超声反响器应用的重要途径。 本文基于Mettin等关于次Bjerknes力的基本理论,研究了超声反响器驻波声场中两个气泡的振动特征和两个气泡之间的次Bjerknes力的特征及影响因素,为驻波声场中泡群构成的气泡构造和空化效应提供了理论根据,也为超声反响器的应用提供了相关理论。 1 、驻波声场中双气泡动力学方程 超声反
7、响器中,入射声波可表示为 Pi=PAej(t-kx) (1) 式中: PA为声压幅值; 为声波角频率; k为声波波数; x为声场中的位置。假设入射声波在超声反响器中完全反射,根据声波的叠加原理,合成声场为驻波声场,合成声压可表示为 Pdr=2PAcos kxejt (2) 假如声场中的两个气泡在振动的经过中仍然保持球形,则气泡振动的动力学方程能够用Keller-Miksis方程14来表示,详细形式为 式中: R1和R2分别为任意时刻气泡1、气泡2的半径; d为两个气泡之间的距离; 为液体的密度; c为液体中的声速; Ps1可表示为 Ps1=(P0+2R10)(R10R1)3?2R1?4R?1R
8、1?Pdr?P0(5) 式中: P0为液体的静压力; 为液体的外表张力; 为液体的黏滞系数; 为气体绝热压缩系数。声场中,两个气泡之间的互相作用力通常称为次Bjerknes力,可表示为 FB=?F12?=?4d2?V?1V?2?er (6) 式中: 为时间平均; F12为气泡1遭到气泡2的作用力; V1和V2分别为气泡1、气泡2的体积,为了方便描绘叙述次Bjerknes力的方向,通常引入次Bjerknes力系数,其定义为 fB=4?V?1V?2? (7) fB的正负代表了气泡之间的力为吸引力或者排挤力,假如fB 0,则两个气泡互相吸引,假如fB 0,则两个气泡互相排挤。 2、 驻波声场中气泡的
9、非线性振动与气泡之间的互相作用 本文的主要研究结果是基于双空化气泡的耦合振动方程式(3)和式(4)得到的,式(3)、式(4)分别是基于气泡1、气泡2半径的二阶非线性常微分方程,可通过数值迭代法得到它们的数值解15,16。本文采用的是软件MATLAB中基于四阶五级Runge-Kutta变步长算法,这种算法能够采用误差作为检测指标的闭环控制,既保证有较高的运算速度和计算精度,也能保证预期的数值稳定性17,是求解气泡非线性常微分方程数值解的重要方式方法。数值模拟中,气泡1、气泡2半径的二阶常微分方程的初始条件为R1(t=0)=R10, R2(t=0)=R20, dR1/dt(t=0)=0, dR2/
10、dt(t=0)=0。 在下面的数值模拟中,我们假设两个气泡之间的距离保持不变,只考虑入射波实部对气泡的影响,在这里基础上我们分别研究了入射波声压振幅为0.9105 Pa, 驱动频率20 kHz的平面波声场和驻波声场中初始半径为2 m的气泡的最大半径随声场位置的变化关系,如此图1所示。图1中,曲线、曲线分别为驻波声场和平面声场中气泡最大半径随声场位置的变化曲线,数值模拟结果表示清楚在驻波声场中,由于入射声波和反射声波的叠加,在不同的位置处,会出现驱动声压振幅的不同,不同位置处气泡的收缩振动和气泡的初始位置有关,处于驻波声场的声压波腹x=0,x=0.037,x=0.074(0,1/2声波长,1个声
11、波长)附近的气泡能够做剧烈的收缩振动,气泡在声压波腹处的最大振动半径能够到达初始半径的42.72倍。而同样初始半径的气泡在平面声场中,由于驱动声压振幅与声场中的位置无关,因此气泡在整个收缩扩张经过中如曲线所示,不同位置处的气泡在收缩振荡经过中,气泡的最大扩张半径均一样,在驱动声压振幅为0.9105 Pa时,整个平面声场中一个驱动周期内,气泡的最大半径只要初始半径的1.41倍左右,讲明小振幅驱动的平面声场中,气泡只在初始半径周围做微弱的收缩振动,不发生剧烈的收缩溃灭现象。上述研究结果表示清楚:驻波声场中,由于声波的叠加,使声场中出现了声压振幅的变化,造成了声场的不均匀性,但也是由于声压的叠加,使
12、小振幅驱动下在声压振幅极大值位置处的气泡能够剧烈振动,产生剧烈的空化效应。 图1 气泡1的径向振动与气泡位置的关系曲线 Fig.1 Curve of radial oscillation of the first bubble and the bubble position 当驱动频率为20 kHz 时,驻波声场和平面声场中,x=1/50 声波长,初始半径为2 m, 相距200 m 的两个气泡组成的气泡对中气泡1的Rmax/Rmin随声压振幅的变化曲线, 如此图2所示。 Rmax/Rmin反映了气泡的能量转换能力18,比照能够看出,同样的入射波声压振幅下,驻波声场中的气泡在声压波腹附近具有更好
13、的能量转换能力,能够在小振幅入射声压振幅下做剧烈的非线性振动,例如:当入射声压振幅只要0.85 Pa时,驻波声场中气泡的Rmax/Rmin即可到达89.09,而同样的入射声压下,平面声场中气泡的Rmax/Rmin只要1.467,此时的气泡只做微弱的收缩振动。驻波声场中,气泡能量转换能力的加强,源于驻波声场中,由于反响器的反射,使反响器部分区域的驱动声压振幅加倍,进而使气泡的非线性振动加强,但是数值模拟结果表示清楚:在驻波声场中,气泡固然能够在小振幅驱动声压下实现良好的能量转换能力,但仍然对驱动声压振幅有一定的要求。在上述条件下,我们发现当驱动声压振幅低于0.65105 Pa下面,气泡在驻波声场
14、中波腹位置处不能做强非线性振动,气泡振动的最大半径与最小半径的比值在1.739左右,这表示清楚在上述条件下,当驻波声场的驱动声压振幅低于0.65105 Pa时,即便处在波腹位置处的气泡也不能出现剧烈的空化现象。 图2 气泡1膨胀收缩比随驱动声压振幅的变化曲线 Fig.2 Curve of contraction ratio and amplitude of driving pressure 声场中,两个气泡之间的互相作用力与气泡在声场中的位置有很大的关系,图3是驱动声压振幅为0.8105 Pa ,驱动频率为 20 kHz的驻波声场和平面波声场中两个初始半径为2 m, 相距200 m的气泡在一个
15、声驱动周期内次Bjerknes力随气泡位置(1个声波长)的变化关系曲线图,如此图3所示。图3中:曲线为驻波声场;曲线为平面声场;x=0.018 6(1/4声波长)和x=0.055 9(3/4声波长)为波节位置;x=0.037 2(1/2声波长)和x=0.074 5(1个声波长)为波腹位置。数值模拟结果表示清楚:在驻波声场中,处于波腹位置处的两个气泡之间的互相作用最大,而波节位置处两个气泡之间的互相作用力趋于0,可忽略不计。详细原因是:处于声压波腹处的气泡能够做剧烈的振动,其产生的声场能对周围气泡产生更大的作用力,而位于声压波节位置处的气泡只能做微弱的线性振动,其对周围气泡的互相作用力也非常微弱
16、,可忽略不计。数值模拟结果还表示清楚:驻波声场中,一个声波长范围内,气泡位于0和1个声波长位置处,两个气泡之间的互相作用力最大,气泡位于1/2个声波长位置处,两个气泡之间的互相作用力次之;而在平面波声场中,在上述驱动条件下,两个气泡之间的次Bjerknes力近似为0。以上结果讲明在小振幅驱动的声场中,声波的叠加使气泡的非线性振动加强,进而使气泡之间的次Bjerknes力增加,气泡之间的互相作用不能够忽略不计,这可能也是超声反响器中构成气泡构造的主要原因。 图3 气泡之间的次Bjerknes力随气泡位置的变化曲线 Fig.3 Curve of the secondary Bjerknes for
17、ce and the bubble position 驻波声场中气泡1对气泡2的作用力F12与两个气泡的振动状态有关。驱动声压振幅为0.8105 Pa, 驱动频率为20 kHz的驻波声场中,位于1/50个声波长位置处的两个初始半径为2 m和3 m的气泡的作用力F12和每个气泡的收缩振动速度与时间的关系,如此图4所示,主要数值模拟结果如表1所示。在0.72个声驱动周期时,气泡1处于收缩状态,气泡2处于扩张状态,此时气泡1对气泡2的作用力F12为排挤力;在0.82个声驱动周期时,气泡1处于迅速扩张状态,气泡2处于扩张状态,此时F12为吸引力,计算结果表示清楚:气泡的收缩振动速度越大,气泡1对气泡2
18、的作用力越大。在平面波声场中,同样的驱动条件下,两个初始半径为2 m和3 m的气泡之间的瞬时作用力也遵循同样的规律,如此图5所示。在0.01个驱动周期时,气泡1的振动速度为-1.67 m/s, 表示清楚气泡1处于收缩状态,气泡2的振动速度为1.31 m/s, 表示清楚气泡2处于扩张状态,此时气泡1对气泡2的作用力F12为1.7610-11 N,为排挤力,当两个气泡都处于扩张或者收缩状态时,F12为吸引力,而当两个气泡中一个处于扩张状态,另一个处于收缩状态时,气泡1对气泡2的作用力为排挤力。 图4 驻波声场气泡之间的作用力F12和振动速度随时间的变化曲线 Fig.4 Curve of the a
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 农业相关
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内