-
当声波在含气泡的液体中传播时会出现共振传播现象, 即在气泡的共振频率附近声衰减和声速会显著地增大,这是声空化领域的一个重要现象.以往的研究一般假设液体中只存在单一种类的气泡, 因此忽略了声波共振传播的某些重要信息. 本文研究了含混合气泡液体中声波的共振传播, 混合气泡是指液体中包含多种静态半径不同的气泡. 结果显示:在这种系统中存在声波共振传播的抑制效应, 即与含单一种类气泡的系统相比, 在含混合气泡的系统中声波的共振衰减和共振声速会明显变小. 对于两种气泡混合、多种气泡混合以及气泡满足某种连续分布的系统, 研究了抑制效应的本质和主要特征, 此外还探究了黏性和空化率等对抑制效应的影响. 本文的研究结果是对该领域现有知识的必要补充.There is the resonant propagation phenomenon of acoustic wave in bubbly liquid, i.e., the attenuation coefficient and the velocity of acoustic wave in range of resonant frequencies of bubbles can become very large. In previous papers, generally adopted was a simplified assumption that there is a single type of bubble in a liquid. It restricts our understanding of the resonant propagation phenomenon. In this paper the resonant propagation of acoustic wave in a liquid with mixed bubbles is studied. Here, static radii of bubbles are different from each other. Research results show that there is a restraining effect of the resonant propagation of acoustic wave in liquid with mixed bubbles. The attenuation coefficient and the velocity of acoustic wave in the liquid with mixed bubbles are obviously less than those in the liquid containing bubbles with the same static radius. The nature of the restraining effect is that the resonant vibration of bubbles is weakened due to the interaction between bubbles with different static radii. Some important properties of the restraining effect are investigated for all kinds of liquid systems with mixed bubbles. Moreover, the effect of the viscosity and the rate of cavitation on the restraining effect are also studied. Research results are shown as follows. 1) Comparing with bigger bubbles, resonant characteristic quantities (such as the attenuation coefficient and the velocity of acoustic wave) caused by smaller bubbles can be reduced more obviously; 2) the efficiency of the restraining effect increases with the number of types of bubbles increasing, however, it will gradually reach to a stable value when the number of types of bubbles is large; 3) the bandwidth of the resonant absorption of acoustic wave is dramatically affected by the distribution function of the percentage of the number of bubbles. The bandwidth of the resonant absorption will become large as the percentage of the number of smaller bubbles increases.
-
Keywords:
- mixed bubbles /
- acoustic cavitation /
- attenuation coefficient /
- linear acoustic wave
1. 引 言
声波在含气泡液体中的传播是声空化领域的重要问题, 是很多声波应用的理论基础, 如声化学、超声处理和超声成像等. 对这类问题的研究有助于深刻理解空化器件内声场的分布情况, 为超声处理方法的大型化和工业化[1-4]提供解决的思路.
在含气泡液体中声波的一个重要特点是会出现共振传播现象, 即在气泡的共振频率附近声衰减和声速会显著地增大[5,6]. 这是一个普遍的现象, 是很多应用(或者潜在运用)的基础, 但是目前对它的认识并不全面. 关于含气泡液体中声波的传播问题, 很多文献忽略了共振传播现象, 如文献[7—10]; 文献[6, 11—18]虽然考虑了共振传播现象, 但是假设液体中所有气泡的静态半径(没有波动时气泡的半径)是相同的, 即液体中只包含单一种类的气泡, 因此遗漏了共振传播方面很多重要的信息; 文献[5]虽然涉及含混合气泡液体中声波的共振传播, 但是只是简单的介绍, 并没有给出系统深入的分析, 因此很多重要的知识并没有被揭示.
近年来张鹏利等[19]、苗博雅和安宇[20]、王德鑫和那仁满都拉[21]对双气泡(不同静态半径)系统的振动问题进行了调查, 揭示了一些有价值的信息, 但是这些文献并未考虑声波的传播问题(气泡振动和声波传播是两类不同的问题). 显然在包含混合气泡(液体中存在多种不同静态半径的气泡)的液体中, 不同静态半径气泡间的相互作用会对声传播产生重要的影响, 并完善声波共振传播的知识.
本文运用有效介质理论系统地研究了含混合气泡液体中声波的共振传播性质. 研究结果显示在这些系统中存在声波共振传播的抑制效应, 即和含单一种类气泡的系统相比, 在含混合气泡的系统中声波的共振衰减和共振声速会明显地变小. 本文系统地研究了这种抑制效应的本质和主要特点, 此外还考虑了空化率和黏性等对抑制效应的影响, 研究结果是对含气泡液体中声波共振传播知识的必要补充.
2. 含混合气泡液体中声波的有效介质理论
在含混合气泡的液体中, 声波的波动方程可以表示为[6]
(1−β0)c20∂2P1∂t2−∇2P1−3μρ0c20∇2∂P1∂t=ρ0∂2β1∂t2−3μ(1−β0)∇2∂β1∂t, (1) 其中
P 是某个时空点的压强,ρ 是纯液体的密度,β 是单位体积的系统包含的气泡的体积(即气体的体积分数); 这些量可以展开为如下的形式:{P=P0+P1+P2+⋯ρ=ρ0+ρ1+ρ2+⋯β=β0+β1+β2+⋯, (2) 其中
Fi 是物理量F 的第i 阶小量;F0 是静态时(即不存在声波时)系统的各物理量;P0 ,ρ0 和c0 分别是静态时系统的压强、纯液体的密度和纯液体中的声速;μ 是液体的黏滞系数(体积黏滞系数, 或称为黏度);β0 是不存在声波时系统的体积分数, 称为空化率.假设液体中包含
M 种不同静态半径的气泡, 第j 种气泡的半径为Rj , 数密度是nj (单位体积内包含第j 种气泡的数量), 数量百分比是θj=nj/n (n 是气泡总的数密度,n=∑Mj=1nj ). 气体的体积分数可以表示为β=43πn∑Mj=1R3jθj. (3) Rj 可以展开为如下的形式:Rj=Rj0+Rj1+Rj2+⋯, (4) 其中
Rj0 是静态时第j 种气泡的半径,Rj1 是第j 种气泡半径的1阶小量,Rj2 是第j 种气泡半径的2阶小量. 把方程(4)代入方程(3), 可得β0=43π∑Mj=1R3j0nj,β1=4π∑Mj=1R2j0Rj1nj. (5) 第
j 种气泡的振动方程可以用Keller和Kolodner方程来描写[5,6,22]:(1−˙Rjc0)Rj¨Rj+32(1−˙Rj3c0)˙R2j=1ρ0(1+˙Rjc0+Rjc0ddt)(PBj−P), (6) 其中
j=1,2,⋯,M .;˙Rj 和¨Rj 分别是第j 种气泡的半径对时间的一阶和二阶导数;P 是气泡所在位置处声波的压强, 它由方程(1)确定.PBj 由下面的方程确定PBj=PGj−2σRj−4μ˙RjRj, (7) 其中
PGj 是第j 种气泡内气体的压强, 它随时间变化;σ 是液体的表面张力系数.本文只讨论声波的线性传播, 所有的一阶小量都包含一个和时间有关的因子, 即
exp(iωt) , 其中i 是单位虚数,ω 是声波的角频率. 这里, 所有的零阶量都是和时间无关的常数. 如果考虑气泡和液体之间的热交换[6,22], 在小振幅声波驱动近似下,PGj 可以表示为PGj=PGj0(1−φjRj1Rj0), (8) 其中
PGj0=2σRj0+P0 ,Rj0 和PGj0 分别是静态时第j 种气泡的半径和泡内气体的压强.φj 是一个和角频率ω 有关的复数,φj=3γ1−3(γ−1)iχj[gjcoth(gj)−1], (9) 其中
χj=D/(ωR2j0) ,D 是气体的热扩散系数,gj=(i/χj)1/2 ,γ 是气体的等压和等容热容量之比,coth() 是双曲余切函数.若要把气泡的振动方程(6)线性化, 则需将方程(2)中的
P 和方程(4)的R 的展开式代入方程(6), 且只保留到一阶小量. 因为ωRj0c0≪1 , 所以可以忽略(ωRj0c0)2 有关的项. 最后可以得到[6]Rj1=fjP1, (10) 其中
fj=−(ω2j0−ω2+2ibjω)−11ρ0Rj0 ,ω2j0=PGj0ρ0R2j0[Re(φj)−2σRj0PGj0] ,bj=2μρ0R2j0+PGj02ωρ0R2j0Im(φj)+ω2Rj02c0 ,Re(φj) 和Im(φj) 分别表示对φj 取实部和虚部. 方程(10)是方程(6)的线性化形式, 它表示气泡半径对声压的响应.β1=EP1, (11) 其中
E=4π∑Mj=1R2j0Rj1njfj .∇2P1+b∂2P1∂t2+c∇2∂P1∂t=0, (12) 其中
b=ρ0E−(1−β0)c20 ,c=3μρ0c20−3μE1−β0 .假设平面声波向
x 方向传播, 那么声压可以表示为P1=ˉP1exp(iωt)exp(−ikx), (13) 方程(13)是用复数形式来表示压强
P1 , 真实的压强应该是P1 的实部. 把方程(13)代入方程(12)可得k 的表达式k2=−bω2(1+icω), (14) 从方程(14)可见
k 是一个复数. 可以把它表示为k=kr+iki , 其中kr>0 ,ki<0 . 那么声压可以表示为P1=ˉP1exp(iωt−ikrx)exp(kix) . 这个平面波的声速是v=ω/kr , 衰减系数是α=−ki .3. 数值计算例子与讨论
本节用数值计算例子来显示含混合气泡液体中线性声波的传播特性. 在所有的计算中, 如果文中或者图中没有特别标明, 那么用到的参数是:
γ=1.4 ,D=2.4×10−5m2/s ,c0=1.5×103m/s ,P0=1.01×105Pa ,σ=10−2N/m ,μ=2×10−3Pa⋅s ,β0=0.0001 . 对于含两种混合气泡的液体系统, 第一种气泡的静态半径是R1=20μm , 第二种气泡的静态半径是R2=60μm . 如果文中或者图中对某个参数有特别标明, 那么该参数就替换为标明的数据.图1显示了含两种混合气泡(
M=2 )的液体中声波的衰减谱. 虚线对应θ1=1 的情况, 此时液体中只存在第一种气泡(静态半径为R1 ); 点划线对应θ1=0 (或者θ2=1 )的情况, 此时液体中只存在第二种气泡(静态半径为R2 ). 小气泡(静态半径小的气泡, 即R1 )对应高频的共振衰减峰(在其他的文献[6]中已经阐明共振衰减峰是由气泡的共振引起的, 以后简称衰减峰), 大气泡(静态半径大的气泡, 即R2 )对应低频的衰减峰. 点线和实线表示液体中包含两种混合气泡的情况, 此时衰减谱存在两个峰.与只存在一种气泡的情况相比(在空化率相同的条件下), 存在两种混合气泡时衰减峰的峰值会明显变小. 当
θ1=1 时(只存在小气泡),α 的峰值是αmax1=1797.7m−1 (对应图1的第二个衰减峰); 此时衰减峰的带宽是Δω1 , 它是α 大于或等于峰值的一半所对应的角频率间隔. 当θ1=0 时(只存在大气泡),α 的峰值是αmax2=708.1m−1 , 衰减峰的带宽是Δω2 ; 当θ1=0.95 时(液体中存在混合的小气泡和大气泡), 衰减谱的第二个峰的峰值是1095.9m−1 , 它小于αmax1 , 此时第一个峰的峰值也小于αmax2 ; 当θ1=0.9 时, 衰减谱的第一个峰的峰值是589.1m−1 , 它小于αmax2 , 此时第二个峰的峰值也小于αmax1 . 这是一个普遍的现象, 称为混合气泡对声波共振衰减的抑制效应. 众所周知, 气泡的共振会导致声衰减的极大增加, 气泡振动越强声衰减也会越大[6]. 在含混合气泡的液体中, 声波共振衰减抑制的根源是气泡振动被抑制了.图2显示了在含两种混合气泡的液体中, 声波衰减谱的相对峰值和相对带宽与
θ1 (或者θ2 )的关系. 图2(a)对应着衰减谱中的第二个峰(它由小气泡的共振引起), 其相对峰值是衰减峰峰值除以αmax1 , 相对带宽是带宽除以Δω1 . 图2(b)对应着第一个峰(它由大气泡的共振引起), 其相对峰值是衰减峰峰值除以αmax2 , 相对带宽是带宽除以Δω2 .由图2(a)可见, 当气泡1(小气泡)的百分比(
θ1 )由1逐渐减少时, 第二个衰减峰的相对峰值和相对带宽都会持续减少. 由图2(b)可见, 当气泡2(大气泡)的百分比(θ2 )由1逐渐减少时, 第一个衰减峰的相对峰值和相对带宽也都会持续减少. 这意味着混合气泡对声波共振衰减的抑制效应表现在两个方面: 其一是衰减峰的峰值降低, 其二是衰减峰的带宽减少.从图2还可以看出, 混合气泡对小气泡的共振特征量(如气泡共振引起的声衰减系数和声速等)的抑制效果要比对大气泡的大得多. 对第二个衰减峰而言(图2(a)), 当小气泡的百分比
θ1 由1减少到0.95时, 其相对峰值由1减少到0.61, 相对峰值的抑制率是7.8(某物理量的抑制率定义为物理量的减少量除以气泡百分比的减少量, 抑制率越大抑制效果越大); 相对带宽由1减少到0.688, 抑制率是6.24. 对第一个衰减峰而言(图2(b)), 当大气泡的百分比θ2 由1减少到0.20时, 相对峰值由1减少到0.917, 抑制率是0.104; 相对带宽由1减少到0.847, 抑制率是0.19. 小气泡共振衰减的相对峰值和相对带宽的抑制率分别是大气泡的75倍和32.84倍.图3显示了在含两种混合气泡的液体中声波的声速谱. 粗实线对应着
θ2=1 , 即只存在大气泡的情况. 点划线对应着θ1=1 , 即只存在小气泡的情况. 与只存在一种气泡的情况相比, 存在混合气泡时声速的峰值要明显减小, 称它为混合气泡对共振声速的抑制效应. 众所周知气泡的共振会导致声速的增加, 共振振动越强声速增加越大[6], 共振声速抑制意味着气泡振动被抑制了. 共振声速或者声衰减抑制效应的本质都是气泡振动的抑制, 以后把它们统称为声波共振传播的抑制效应. 从图3还可以看出混合气泡对小气泡的共振声速的抑制效果要远比对大气泡的大. 如在声速谱的第二个峰(由小气泡共振引起)处,θ1 由1.00变成0.95, 共振声速峰值由10055.7m⋅s−1 变成5561.5m⋅s−1 , 其抑制率是89884m⋅s−1 . 在第一个峰处(由大气泡共振引起),θ2 由1.00变成0.80, 共振声速峰值由13207.4m⋅s−1 变成12365.8m⋅s−1 , 其抑制率是4208m⋅s−1 . 小气泡共振声速峰值的抑制率是大气泡的21.36倍.图 3 在含两种混合气泡的液体中, 不同 (或者θ1 )对应的声波的声速谱 图中θ2 分别是0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95, 1.00, 或者θ1 分别是0, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80, 1.00. 左边的箭头表示θ2 增加的方向, 右边的箭头表示θ2 增加的方向θ1 Fig. 3. Speed spectrums for different (orθ1 ) in liquids with two kinds of bubbles. Hereθ2 is 0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95 and 1.00, respectively.θ1 is 0, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80 and 1.00, respectively. The white arrow indicates the direction of increase ofθ2 , the gray arrow indicates the direction of increase ofθ2 .θ1 图4显示了在含两种混合气泡的液体中, 不同
β0 下声波的相对衰减谱. 当β0 不变时, 衰减系数是角频率的函数α(β0,ω) , 其最大值是αmax(β0) . 相对衰减系数定义为αr=α(β0,ω)/αmax(β0) .αmax(β0) 随β0 的增加而增加, 当β0 为10−6 ,10−4 ,10−3 和10−2 时,αmax(β0) 分别是15.84 ,649.11 ,2217.10 和7231.88m−1 . 这是比较容易理解的, 单个气泡对声波有一定的耗散作用, 当β0 增加时, 单位体积内气泡的数量增多, 因此气泡对声波的衰减变大.由图4可见, 在相对衰减谱的两个峰之间的(角频率)区域内,
αr 随β0 的增加而增加, 但是它们会逐渐趋于一个特定的曲线(如图中的实线). 当β0 再继续增加时, 该区域的αr 曲线和图中的实线区别很小, 无法在图中显示出来. 当β0 增加时(单位体积内气泡的数量增多), 气泡间的平均距离减少, 气泡间的相互作用力变大, 相对衰减谱中两个峰的带宽变大[6], 以至于在两个峰之间的区域内声波的相对衰减增强. 但是这种增强效应具有饱和性, 即当β0 增加到一定值时, 声波的相对衰减会趋于稳定.图5显示了在含两种混合气泡的液体中, 不同
β0 对应的声速谱. 在声速谱的两个峰处, 声速先随β0 的增加而增加, 然后达到最大值, 再随β0 的增大而减少. 可见β0 对声波的共振衰减和声速的影响是不同的, 前者随β0 的增加而单调增加, 后者却是先增加而后减少. 在两个峰之间的区域内, 声速随β0 的增加而单调地减少. 由图4可知, 在这个区域内相对衰减随β0 的增加而增大, 但是声速却正好相反, 这是一个比较反常的现象.图6显示了在含两种混合气泡的液体中黏度不同时声波的衰减谱. 由图可见: 1)在两个衰减峰处, 衰减系数的峰值随黏度的增加反而减小. 在这两个区域内, 声衰减主要由气泡的共振引起, 而黏性对气泡的共振起抑制作用(即黏度越大共振越弱), 因此衰减系数的峰值随黏度的增加反而减少; 2)在两个衰减峰之间的区域内, 衰减系数随黏度的增加而增加. 这说明在这个区域内声衰减主要是由液体的黏性引起, 而不是由气泡的共振引起.
图7 显示了在含多种混合气泡的液体中声波的衰减谱. 由图可见: 1)在含多种混合气泡的液体中声波共振传播的抑制效应依然存在, 且混合气泡的种类越多, 声波共振传播的抑制效应越明显, 即当
M 越大时, 衰减系数的峰值变得越小; 2)当M 比较小时, 衰减谱上存在多个衰减峰(一般是M 个峰), 当M 增加到一定值时, 衰减谱变成一条光滑的曲线(如M=20 时, 图7中的粗实线). 再继续增加M , 衰减谱曲线的变化很小, 无法在图上显示出来. 这说明对M 而言, 声波共振传播的抑制效应存在饱和现象. 当M 很大时, 抑制效果会趋于稳定, 而不是持续增加.图8显示了气泡数满足不同分布时声波的衰减谱. 点线对应的分布是
θM−i+1=ai+b (其中b=0.1/M ,a=(1−bM)/∑Mi=1i ), 这是一个气泡数偏向小气泡的分布, 衰减谱的带宽是6.032×105Hz ; 粗实线对应的分布是θj=1/M , 这是一个均匀的分布, 衰减谱的带宽是4.083×105Hz ; 虚线对应的分布是θi=ai+b , 这是一个气泡数偏向大气泡的分布, 衰减谱的带宽是3.106×105Hz . 可以得出如下的结论: 气泡数分布越偏向大气泡, 共振衰减峰的带宽就越小, 反之则越大.图 8 气泡满足不同分布时声波的衰减谱 细实线是参考曲线, 它对应着存在两种混合气泡的情况; 对其他曲线而言, 各个气泡的静态半径分别是R = 20 + 40(j – 1)/ ,(M−1)(μm) ,j=1,2,⋯,M M=21 Fig. 8. Attenuation spectrums for different distribution functions of bubble number. The thin solid line indicates the systems with two kinds of bubbles. Here, R = ,20+40(j−1)/(M−1)(μm) .j=1,2,⋯,M 这里试着解释上面的结论, 前面已知小气泡共振特征量的抑制效果远比大气泡的大. 当气泡数分布偏向大气泡时, 小气泡的数量就比较少, 小气泡对应的高频共振衰减被极大地抑制了, 因此衰减谱的带宽比较小. 当气泡数分布偏向小气泡时, 小气泡的数量比较大, 小气泡对应的高频共振衰减的抑制效果较小; 虽然大气泡的数量较小, 但是大气泡的特征量的抑制效果本来就弱, 因此大气泡对应的低频共振衰减的抑制效果也较小. 结果是衰减峰带宽比较大.
4. 结 论
本文运用有效介质理论研究了含混合气泡液体中声波的共振传播特性, 得出的主要结论概括如下.
1)在含混合气泡的液体中存在声波共振传播的抑制效应. 这是一个普遍的现象, 具体的表述是: 和含单一种类气泡的液体相比, 在含混合气泡的液体中, 声波的共振衰减系数和声速都会明显地变小. 这种抑制效应的本质是液体中气泡的共振振动被抑制了.
2)在含混合气泡的液体中, 声波共振传播的抑制效应的主要特征有: 小气泡共振特征量的抑制效果远比大气泡的大; 当气泡数分布越偏向大气泡时, 衰减峰的带宽越小, 反之则越大; 抑制效果随
M 的增加而增加, 但是存在饱和效应, 即当M 很大时声波共振传播的抑制效果趋于稳定.3)在含混合气泡的液体中, 各种因素对声衰减谱和声速谱的主要影响有: 在声衰减谱的衰减峰处黏性对声衰减起抑制作用, 而在衰减峰之间的区域内黏性对声衰减起促进作用; 在衰减峰处空化率对声波的共振衰减和声速的影响是不同的, 前者随空化率的增加而单调增加, 后者却是先增加而后减少; 在衰减峰之间的区域内相对衰减随空化率的增加而增加, 但是这种增加具有饱和性, 即空化率达到较大值时, 相对衰减趋于稳定; 在衰减峰之间的区域内声速随空化率的增加而减少.
[1] 陈伟中 2014 声空化物理 (北京: 科学出版社)第371−421页
Chen W Z 2014 Acoustic Cavitation Physics (Beijing: Science Press) pp371−421 (in Chinese)
[2] Christopher E B 1995 Cavitation and Bubble Dynamics (Oxford: OxfordvUniversity Press) pp15−47
[3] 应崇福 2008 应用声学 27 333
Google Scholar
Ying C F 2008 J. Appl. Acoust. 27 333
Google Scholar
[4] 王寻, 陈伟中, 杨景 2015 声学技术 34 33
Wang X, Chen W Z, Yang J 2015 Tech. Acoust. 34 33
[5] Commander K W, Prosperetti A 1989 J. Acoust. Soc. Am. 85 732
Google Scholar
[6] 徐贞, 张迪, 陈时, 等 2018 中国科学: 物理学 力学 天文学 48 044301
Xu Z, Zhang D, Chen S, et al. 2018 Sci. Sin. Phys. Mech. Astron. 48 044301
[7] 王勇, 林书玉, 张小丽 2013 物理学报 62 064304
Wang Y, Lin S Y, Zhang X L 2013 Acta Phys. Sin. 62 064304
[8] 王成会, 林书玉 2008 陕西师范大学学报 36 30
Google Scholar
Wang C H, Lin S Y 2008 J. Shaanxi Norm. Univ. 36 30
Google Scholar
[9] 王勇, 林书玉, 张小丽 2014 物理学报 63 034301
Google Scholar
Wang Y, Lin S Y, Zhang X L 2014 Acta Phys. Sin. 63 034301
Google Scholar
[10] 陈伟中 2018 应用声学 37 675
Google Scholar
Chen W Z 2018 J. Appl. Acoust. 37 675
Google Scholar
[11] 朱哲民, 杜功焕 1995 声学学报 6 425
Zhu Z M, Du H G 1995 Acta Acust. 6 425
[12] An Y 2012 Phys. Rev. E 85 016305
Google Scholar
[13] Vanhille C, Campospozuelo C 2013 Ultrason. Sonochem. 20 963
Google Scholar
[14] Vanhille C, Campospozuelo C 2012 Ultrason. Sonochem. 19 217
Google Scholar
[15] Lebon G S B, Tzanakis I, Djambazov G, et al. 2017 Ultrason. Sonochem. 37 660
Google Scholar
[16] Zhang Y, Guo Z, Du X 2018 Appl. Therm. Eng. 133 483
Google Scholar
[17] Zhang Y, Du X 2015 Ultrason. Sonochem. 26 119
Google Scholar
[18] Trujillo F J 2018 Ultrason. Sonochem. 47 75
Google Scholar
[19] 张鹏利, 林书玉, 张涛 2013 中国科学: 物理学 力学 天文学 43 249
Zhang P L, Lin S Y, Zhang T 2013 Sci. Sin. Phys. Mech. Astron. 43 249
[20] 苗博雅, 安宇 2015 物理学报 64 225
Google Scholar
Miao B Y, An Y 2015 Acta Phys. Sin. 64 225
Google Scholar
[21] 王德鑫, 那仁满都拉 2018 物理学报 67 037802
Google Scholar
Wang D X, Naranmandula 2018 Acta Phys. Sin. 67 037802
Google Scholar
[22] Keller J B, Kolodner I I 1956 J. Appl. Phys. 2 71152
-
图 3 在含两种混合气泡的液体中, 不同
θ1 (或者θ2 )对应的声波的声速谱 图中θ1 分别是0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95, 1.00, 或者θ2 分别是0, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80, 1.00. 左边的箭头表示θ2 增加的方向, 右边的箭头表示θ1 增加的方向Fig. 3. Speed spectrums for different
θ1 (orθ2 ) in liquids with two kinds of bubbles. Hereθ1 is 0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95 and 1.00, respectively.θ2 is 0, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80 and 1.00, respectively. The white arrow indicates the direction of increase ofθ2 , the gray arrow indicates the direction of increase ofθ1 .图 8 气泡满足不同分布时声波的衰减谱 细实线是参考曲线, 它对应着存在两种混合气泡的情况; 对其他曲线而言, 各个气泡的静态半径分别是R = 20 + 40(j – 1)/
(M−1)(μm) ,j=1,2,⋯,M ,M=21 Fig. 8. Attenuation spectrums for different distribution functions of bubble number. The thin solid line indicates the systems with two kinds of bubbles. Here, R =
20+40(j−1)/(M−1)(μm) ,j=1,2,⋯,M . -
[1] 陈伟中 2014 声空化物理 (北京: 科学出版社)第371−421页
Chen W Z 2014 Acoustic Cavitation Physics (Beijing: Science Press) pp371−421 (in Chinese)
[2] Christopher E B 1995 Cavitation and Bubble Dynamics (Oxford: OxfordvUniversity Press) pp15−47
[3] 应崇福 2008 应用声学 27 333
Google Scholar
Ying C F 2008 J. Appl. Acoust. 27 333
Google Scholar
[4] 王寻, 陈伟中, 杨景 2015 声学技术 34 33
Wang X, Chen W Z, Yang J 2015 Tech. Acoust. 34 33
[5] Commander K W, Prosperetti A 1989 J. Acoust. Soc. Am. 85 732
Google Scholar
[6] 徐贞, 张迪, 陈时, 等 2018 中国科学: 物理学 力学 天文学 48 044301
Xu Z, Zhang D, Chen S, et al. 2018 Sci. Sin. Phys. Mech. Astron. 48 044301
[7] 王勇, 林书玉, 张小丽 2013 物理学报 62 064304
Wang Y, Lin S Y, Zhang X L 2013 Acta Phys. Sin. 62 064304
[8] 王成会, 林书玉 2008 陕西师范大学学报 36 30
Google Scholar
Wang C H, Lin S Y 2008 J. Shaanxi Norm. Univ. 36 30
Google Scholar
[9] 王勇, 林书玉, 张小丽 2014 物理学报 63 034301
Google Scholar
Wang Y, Lin S Y, Zhang X L 2014 Acta Phys. Sin. 63 034301
Google Scholar
[10] 陈伟中 2018 应用声学 37 675
Google Scholar
Chen W Z 2018 J. Appl. Acoust. 37 675
Google Scholar
[11] 朱哲民, 杜功焕 1995 声学学报 6 425
Zhu Z M, Du H G 1995 Acta Acust. 6 425
[12] An Y 2012 Phys. Rev. E 85 016305
Google Scholar
[13] Vanhille C, Campospozuelo C 2013 Ultrason. Sonochem. 20 963
Google Scholar
[14] Vanhille C, Campospozuelo C 2012 Ultrason. Sonochem. 19 217
Google Scholar
[15] Lebon G S B, Tzanakis I, Djambazov G, et al. 2017 Ultrason. Sonochem. 37 660
Google Scholar
[16] Zhang Y, Guo Z, Du X 2018 Appl. Therm. Eng. 133 483
Google Scholar
[17] Zhang Y, Du X 2015 Ultrason. Sonochem. 26 119
Google Scholar
[18] Trujillo F J 2018 Ultrason. Sonochem. 47 75
Google Scholar
[19] 张鹏利, 林书玉, 张涛 2013 中国科学: 物理学 力学 天文学 43 249
Zhang P L, Lin S Y, Zhang T 2013 Sci. Sin. Phys. Mech. Astron. 43 249
[20] 苗博雅, 安宇 2015 物理学报 64 225
Google Scholar
Miao B Y, An Y 2015 Acta Phys. Sin. 64 225
Google Scholar
[21] 王德鑫, 那仁满都拉 2018 物理学报 67 037802
Google Scholar
Wang D X, Naranmandula 2018 Acta Phys. Sin. 67 037802
Google Scholar
[22] Keller J B, Kolodner I I 1956 J. Appl. Phys. 2 71152
期刊类型引用(2)
1. 赵宁宁,肖新宇,凡凤仙,苏明旭. 基于蒙特卡罗原理的混合颗粒三相体系声衰减计算模型研究. 物理学报. 2022(07): 217-224 . 百度学术
2. 清河美,那仁满都拉. 空化多泡中大气泡对小气泡空化效应的影响. 物理学报. 2019(23): 167-175 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 8043
- PDF下载量: 73
- 被引次数: 8