Loading [MathJax]/jax/output/HTML-CSS/jax.js

搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

陈时 张迪 王成会 张引红

陈时, 张迪, 王成会, 张引红. 含混合气泡液体中声波共振传播的抑制效应. 物理学报, 2019, 68(7): 074301. doi: 10.7498/aps.68.20182299
引用本文: 陈时, 张迪, 王成会, 张引红. 含混合气泡液体中声波共振传播的抑制效应. 物理学报, 2019, 68(7): 074301. doi: 10.7498/aps.68.20182299
Chen Shi, Zhang Di, Wang Cheng-Hui, Zhang Yin-Hong. Restraining effect of resonant propagation of acousticwaves in liquids with mixed bubbles. Acta Phys. Sin., 2019, 68(7): 074301. doi: 10.7498/aps.68.20182299
Citation: Chen Shi, Zhang Di, Wang Cheng-Hui, Zhang Yin-Hong. Restraining effect of resonant propagation of acousticwaves in liquids with mixed bubbles. Acta Phys. Sin., 2019, 68(7): 074301. doi: 10.7498/aps.68.20182299

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

陈时, 张迪, 王成会, 张引红

Restraining effect of resonant propagation of acousticwaves in liquids with mixed bubbles

Chen Shi, Zhang Di, Wang Cheng-Hui, Zhang Yin-Hong
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 当声波在含气泡的液体中传播时会出现共振传播现象, 即在气泡的共振频率附近声衰减和声速会显著地增大,这是声空化领域的一个重要现象.以往的研究一般假设液体中只存在单一种类的气泡, 因此忽略了声波共振传播的某些重要信息. 本文研究了含混合气泡液体中声波的共振传播, 混合气泡是指液体中包含多种静态半径不同的气泡. 结果显示:在这种系统中存在声波共振传播的抑制效应, 即与含单一种类气泡的系统相比, 在含混合气泡的系统中声波的共振衰减和共振声速会明显变小. 对于两种气泡混合、多种气泡混合以及气泡满足某种连续分布的系统, 研究了抑制效应的本质和主要特征, 此外还探究了黏性和空化率等对抑制效应的影响. 本文的研究结果是对该领域现有知识的必要补充.
    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.
      PACS:
      43.35.+d(Ultrasonics, quantum acoustics, and physical effects of sound)
      47.35.Rs(Sound waves)
      47.55.dp(Cavitation and boiling)
      通信作者: 陈时, chenshi@snnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11474191)和中央高校基本科研业务费专项基金(批准号: GK201703012)资助的课题.
      Corresponding author: Chen Shi, chenshi@snnu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11474191), and the Fundamental Research Funds for the Central Universities, China (Grant No. GK201703012).

    声波在含气泡液体中的传播是声空化领域的重要问题, 是很多声波应用的理论基础, 如声化学、超声处理和超声成像等. 对这类问题的研究有助于深刻理解空化器件内声场的分布情况, 为超声处理方法的大型化和工业化[1-4]提供解决的思路.

    在含气泡液体中声波的一个重要特点是会出现共振传播现象, 即在气泡的共振频率附近声衰减和声速会显著地增大[5,6]. 这是一个普遍的现象, 是很多应用(或者潜在运用)的基础, 但是目前对它的认识并不全面. 关于含气泡液体中声波的传播问题, 很多文献忽略了共振传播现象, 如文献[710]; 文献[6, 1118]虽然考虑了共振传播现象, 但是假设液体中所有气泡的静态半径(没有波动时气泡的半径)是相同的, 即液体中只包含单一种类的气泡, 因此遗漏了共振传播方面很多重要的信息; 文献[5]虽然涉及含混合气泡液体中声波的共振传播, 但是只是简单的介绍, 并没有给出系统深入的分析, 因此很多重要的知识并没有被揭示.

    近年来张鹏利等[19]、苗博雅和安宇[20]、王德鑫和那仁满都拉[21]对双气泡(不同静态半径)系统的振动问题进行了调查, 揭示了一些有价值的信息, 但是这些文献并未考虑声波的传播问题(气泡振动和声波传播是两类不同的问题). 显然在包含混合气泡(液体中存在多种不同静态半径的气泡)的液体中, 不同静态半径气泡间的相互作用会对声传播产生重要的影响, 并完善声波共振传播的知识.

    本文运用有效介质理论系统地研究了含混合气泡液体中声波的共振传播性质. 研究结果显示在这些系统中存在声波共振传播的抑制效应, 即和含单一种类气泡的系统相比, 在含混合气泡的系统中声波的共振衰减和共振声速会明显地变小. 本文系统地研究了这种抑制效应的本质和主要特点, 此外还考虑了空化率和黏性等对抑制效应的影响, 研究结果是对含气泡液体中声波共振传播知识的必要补充.

    在含混合气泡的液体中, 声波的波动方程可以表示为[6]

    (1β0)c202P1t22P13μρ0c202P1t=ρ02β1t23μ(1β0)2β1t, (1)

    其中P是某个时空点的压强, ρ是纯液体的密度, β是单位体积的系统包含的气泡的体积(即气体的体积分数); 这些量可以展开为如下的形式:

    {P=P0+P1+P2+ρ=ρ0+ρ1+ρ2+β=β0+β1+β2+, (2)

    其中Fi是物理量F的第i阶小量; F0是静态时(即不存在声波时)系统的各物理量; P0, ρ0c0分别是静态时系统的压强、纯液体的密度和纯液体中的声速; μ是液体的黏滞系数(体积黏滞系数, 或称为黏度); β0是不存在声波时系统的体积分数, 称为空化率.

    假设液体中包含M种不同静态半径的气泡, 第j种气泡的半径为Rj, 数密度是nj(单位体积内包含第j种气泡的数量), 数量百分比是θj=nj/n(n是气泡总的数密度, n=Mj=1nj). 气体的体积分数可以表示为

    β=43πnMj=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)(PBjP), (6)

    其中j=1,2,,M.; ˙Rj¨Rj分别是第j种气泡的半径对时间的一阶和二阶导数; P是气泡所在位置处声波的压强, 它由方程(1)确定. PBj由下面的方程确定

    PBj=PGj2σRj4μ˙RjRj, (7)

    其中PGj是第j种气泡内气体的压强, 它随时间变化; σ是液体的表面张力系数.

    本文只讨论声波的线性传播, 所有的一阶小量都包含一个和时间有关的因子, 即exp(iωt), 其中i是单位虚数, ω是声波的角频率. 这里, 所有的零阶量都是和时间无关的常数. 如果考虑气泡和液体之间的热交换[6,22], 在小振幅声波驱动近似下, PGj可以表示为

    PGj=PGj0(1φjRj1Rj0), (8)

    其中PGj0=2σRj0+P0, Rj0PGj0分别是静态时第j种气泡的半径和泡内气体的压强. φj是一个和角频率ω有关的复数,

    φj=3γ13(γ1)iχj[gjcoth(gj)1], (9)

    其中χj=D/(ωR2j0), D是气体的热扩散系数, gj=(i/χj)1/2, γ是气体的等压和等容热容量之比, coth()是双曲余切函数.

    若要把气泡的振动方程(6)线性化, 则需将方程(2)中的P和方程(4)的R的展开式代入方程(6), 且只保留到一阶小量. 因为ωRj0c01, 所以可以忽略(ω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)的线性化形式, 它表示气泡半径对声压的响应.

    由方程(5)和(10)可得

    β1=EP1, (11)

    其中E=4πMj=1R2j0Rj1njfj.

    把方程(11)代入方程(1), 可得波动方程的最终形式

    2P1+b2P1t2+c2P1t=0, (12)

    其中b=ρ0E(1β0)c20, c=3μρ0c203μ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ωtikrx)exp(kix). 这个平面波的声速是v=ω/kr, 衰减系数是α=ki.

    本节用数值计算例子来显示含混合气泡液体中线性声波的传播特性. 在所有的计算中, 如果文中或者图中没有特别标明, 那么用到的参数是: γ=1.4, D=2.4×105m2/s, c0=1.5×103m/s, P0=1.01×105Pa, σ=102N/m, μ=2×103Pas, β0=0.0001. 对于含两种混合气泡的液体系统, 第一种气泡的静态半径是R1=20μm, 第二种气泡的静态半径是R2=60μm. 如果文中或者图中对某个参数有特别标明, 那么该参数就替换为标明的数据.

    图1显示了含两种混合气泡(M=2)的液体中声波的衰减谱. 虚线对应θ1=1的情况, 此时液体中只存在第一种气泡(静态半径为R1); 点划线对应θ1=0(或者θ2=1)的情况, 此时液体中只存在第二种气泡(静态半径为R2). 小气泡(静态半径小的气泡, 即R1)对应高频的共振衰减峰(在其他的文献[6]中已经阐明共振衰减峰是由气泡的共振引起的, 以后简称衰减峰), 大气泡(静态半径大的气泡, 即R2)对应低频的衰减峰. 点线和实线表示液体中包含两种混合气泡的情况, 此时衰减谱存在两个峰.

    图 1 含两种混合气泡的液体中声波的衰减谱\r\nFig. 1. Attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.
    图 1  含两种混合气泡的液体中声波的衰减谱
    Fig. 1.  Attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.

    与只存在一种气泡的情况相比(在空化率相同的条件下), 存在两种混合气泡时衰减峰的峰值会明显变小. 当θ1=1时(只存在小气泡), α的峰值是αmax1=1797.7m1(对应图1的第二个衰减峰); 此时衰减峰的带宽是Δω1, 它是α大于或等于峰值的一半所对应的角频率间隔. 当θ1=0时(只存在大气泡), α的峰值是αmax2=708.1m1, 衰减峰的带宽是Δω2; 当θ1=0.95时(液体中存在混合的小气泡和大气泡), 衰减谱的第二个峰的峰值是1095.9m1, 它小于αmax1, 此时第一个峰的峰值也小于αmax2; 当θ1=0.9时, 衰减谱的第一个峰的峰值是589.1m1, 它小于αmax2, 此时第二个峰的峰值也小于αmax1. 这是一个普遍的现象, 称为混合气泡对声波共振衰减的抑制效应. 众所周知, 气泡的共振会导致声衰减的极大增加, 气泡振动越强声衰减也会越大[6]. 在含混合气泡的液体中, 声波共振衰减抑制的根源是气泡振动被抑制了.

    图2显示了在含两种混合气泡的液体中, 声波衰减谱的相对峰值和相对带宽与θ1(或者θ2)的关系. 图2(a)对应着衰减谱中的第二个峰(它由小气泡的共振引起), 其相对峰值是衰减峰峰值除以αmax1, 相对带宽是带宽除以Δω1. 图2(b)对应着第一个峰(它由大气泡的共振引起), 其相对峰值是衰减峰峰值除以αmax2, 相对带宽是带宽除以Δω2.

    图 2 含两种混合气泡的液体中声波衰减谱的相对峰值和相对带宽与${\theta _1}$(或者${\theta _2}$)的关系\r\nFig. 2. Relationships between ${\theta _1}$(${\theta _2}$)and relative peak values or bandwidths of attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.
    图 2  含两种混合气泡的液体中声波衰减谱的相对峰值和相对带宽与θ1(或者θ2)的关系
    Fig. 2.  Relationships between θ1(θ2)and relative peak values or bandwidths of attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.

    图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.7ms1变成5561.5ms1, 其抑制率是89884ms1. 在第一个峰处(由大气泡共振引起), θ2由1.00变成0.80, 共振声速峰值由13207.4ms1变成12365.8ms1, 其抑制率是4208ms1. 小气泡共振声速峰值的抑制率是大气泡的21.36倍.

    图 3 在含两种混合气泡的液体中, 不同${\theta _1}$(或者${\theta _2}$)对应的声波的声速谱 图中${\theta _1}$分别是0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95, 1.00, 或者${\theta _2}$分别是0, 0.05, 0.10, 0.20, 0.40, 0.60, 0.80, 1.00. 左边的箭头表示${\theta _2}$增加的方向, 右边的箭头表示${\theta _1}$增加的方向\r\nFig. 3. Speed spectrums for different ${\theta _1}$(or ${\theta _2}$) in liquids with two kinds of bubbles. Here ${\theta _1}$ is 0, 0.20, 0.40, 0.60, 0.80, 0.90, 0.95 and 1.00, respectively. ${\theta _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 ${\theta _2}$, the gray arrow indicates the direction of increase of ${\theta _1}$.
    图 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.

    图4显示了在含两种混合气泡的液体中, 不同β0下声波的相对衰减谱. 当β0不变时, 衰减系数是角频率的函数α(β0,ω), 其最大值是αmax(β0). 相对衰减系数定义为αr=α(β0,ω)/αmax(β0). αmax(β0)β0的增加而增加, 当β0106, 104, 103102时, αmax(β0)分别是15.84, 649.11, 2217.107231.88m1. 这是比较容易理解的, 单个气泡对声波有一定的耗散作用, 当β0增加时, 单位体积内气泡的数量增多, 因此气泡对声波的衰减变大.

    图 4 在含两种混合气泡的液体中, ${\beta _0}$不同时声波的相对衰减谱 图中${\theta _1} = 0.8$, ${\beta _0}$分别是${10^{ - 6}}$, ${10^{ - 4}}$, ${10^{ - 3}}$和${10^{ - 2}}$, 箭头表示该角频率区域${\beta _0}$增加的方向\r\nFig. 4. Attenuation spectrums for different ${\beta _0}$ in liquids with two kinds of bubbles. Here, ${\theta _1} = 0.8$, ${\beta _0}$ is ${10^{ - 6}}$, ${10^{ - 4}}$, ${10^{ - 3}}$ and ${10^{ - 2}}$, respectively. The arrow indicates the direction of increase of ${\beta _0}$.
    图 4  在含两种混合气泡的液体中, β0不同时声波的相对衰减谱 图中θ1=0.8, β0分别是106, 104, 103102, 箭头表示该角频率区域β0增加的方向
    Fig. 4.  Attenuation spectrums for different β0 in liquids with two kinds of bubbles. Here, θ1=0.8, β0 is 106, 104, 103 and 102, respectively. The arrow indicates the direction of increase of β0.

    图4可见, 在相对衰减谱的两个峰之间的(角频率)区域内, αrβ0的增加而增加, 但是它们会逐渐趋于一个特定的曲线(如图中的实线). 当β0再继续增加时, 该区域的αr曲线和图中的实线区别很小, 无法在图中显示出来. 当β0增加时(单位体积内气泡的数量增多), 气泡间的平均距离减少, 气泡间的相互作用力变大, 相对衰减谱中两个峰的带宽变大[6], 以至于在两个峰之间的区域内声波的相对衰减增强. 但是这种增强效应具有饱和性, 即当β0增加到一定值时, 声波的相对衰减会趋于稳定.

    图5显示了在含两种混合气泡的液体中, 不同β0对应的声速谱. 在声速谱的两个峰处, 声速先随β0的增加而增加, 然后达到最大值, 再随β0的增大而减少. 可见β0对声波的共振衰减和声速的影响是不同的, 前者随β0的增加而单调增加, 后者却是先增加而后减少. 在两个峰之间的区域内, 声速随β0的增加而单调地减少. 由图4可知, 在这个区域内相对衰减随β0的增加而增大, 但是声速却正好相反, 这是一个比较反常的现象.

    图 5 在含两种混合气泡的液体中, ${\beta _0}$不同时声波的声速谱 图中${\theta _1} = 0.8$, ${\beta _0}$分别是${10^{ - 6}}$, ${10^{ - 4}}$, ${10^{ - 3}}$和${10^{ - 2}}$, 箭头表示该角频率区域${\beta _0}$增加的方向\r\nFig. 5. Speed spectrums for different ${\beta _0}$ in liquids with two kinds of bubbles. Here, ${\theta _1} = 0.8$, ${\beta _0}$ is ${10^{ - 6}}$, ${10^{ - 4}}$, ${10^{ - 3}}$ and ${10^{ - 2}}$, respectively. The arrow indicates the direction of increase of ${\beta _0}$.
    图 5  在含两种混合气泡的液体中, β0不同时声波的声速谱 图中θ1=0.8, β0分别是106, 104, 103102, 箭头表示该角频率区域β0增加的方向
    Fig. 5.  Speed spectrums for different β0 in liquids with two kinds of bubbles. Here, θ1=0.8, β0 is 106, 104, 103 and 102, respectively. The arrow indicates the direction of increase of β0.

    图6显示了在含两种混合气泡的液体中黏度不同时声波的衰减谱. 由图可见: 1)在两个衰减峰处, 衰减系数的峰值随黏度的增加反而减小. 在这两个区域内, 声衰减主要由气泡的共振引起, 而黏性对气泡的共振起抑制作用(即黏度越大共振越弱), 因此衰减系数的峰值随黏度的增加反而减少; 2)在两个衰减峰之间的区域内, 衰减系数随黏度的增加而增加. 这说明在这个区域内声衰减主要是由液体的黏性引起, 而不是由气泡的共振引起.

    图 6 在两种混合气泡的液体中黏度不同时声波的衰减谱\r\nFig. 6. Attenuation spectrums for different viscosity in liquids with two kinds of bubbles.
    图 6  在两种混合气泡的液体中黏度不同时声波的衰减谱
    Fig. 6.  Attenuation spectrums for different viscosity in liquids with two kinds of bubbles.

    图7 显示了在含多种混合气泡的液体中声波的衰减谱. 由图可见: 1)在含多种混合气泡的液体中声波共振传播的抑制效应依然存在, 且混合气泡的种类越多, 声波共振传播的抑制效应越明显, 即当M越大时, 衰减系数的峰值变得越小; 2)当M比较小时, 衰减谱上存在多个衰减峰(一般是M个峰), 当M增加到一定值时, 衰减谱变成一条光滑的曲线(如M=20时, 图7中的粗实线). 再继续增加M, 衰减谱曲线的变化很小, 无法在图上显示出来. 这说明对M而言, 声波共振传播的抑制效应存在饱和现象. 当M很大时, 抑制效果会趋于稳定, 而不是持续增加.

    图 7 在含多种混合气泡的液体中气泡的种类数不同时声波的衰减谱. 图中${\theta _j} = 1/M$, ${R_j} = 20 \!+ \!40\left( {j\! -\! 1} \right)/\left( {M \!-\! 1} \right)$$\left( {{\text{μ}}{\rm{m}}} \right)$($j = 1, 2, \cdots, M$).\r\nFig. 7. Attenuation spectrums for different $M$. Here, ${\theta _j} = 1/M$, ${R_j} = 20 + 40\left( {j - 1} \right)/\left( {M - 1} \right)\left( {{\text{μ}}{\rm{m}}} \right)$, $j = 1, 2, $$\cdots, M$.
    图 7  在含多种混合气泡的液体中气泡的种类数不同时声波的衰减谱. 图中θj=1/M, Rj=20+40(j1)/(M1)(μm)(j=1,2,,M).
    Fig. 7.  Attenuation spectrums for different M. Here, θj=1/M, Rj=20+40(j1)/(M1)(μm), j=1,2,,M.

    图8显示了气泡数满足不同分布时声波的衰减谱. 点线对应的分布是θMi+1=ai+b(其中b=0.1/M, a=(1bM)/Mi=1i), 这是一个气泡数偏向小气泡的分布, 衰减谱的带宽是6.032×105Hz; 粗实线对应的分布是θj=1/M, 这是一个均匀的分布, 衰减谱的带宽是4.083×105Hz; 虚线对应的分布是θi=ai+b, 这是一个气泡数偏向大气泡的分布, 衰减谱的带宽是3.106×105Hz. 可以得出如下的结论: 气泡数分布越偏向大气泡, 共振衰减峰的带宽就越小, 反之则越大.

    图 8 气泡满足不同分布时声波的衰减谱 细实线是参考曲线, 它对应着存在两种混合气泡的情况; 对其他曲线而言, 各个气泡的静态半径分别是R = 20 + 40(j – 1)/ $\left( {M - 1} \right)\left( {{\text{μ}}{\rm{m}}} \right)$, $j = 1, 2, \cdots, M$, $M = 21$\r\nFig. 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\left( {j - 1} \right)/\left( {M - 1} \right)\left( {{\text{μ}}{\rm{m}}} \right)$, $j = 1, 2, \cdots, M$.
    图 8  气泡满足不同分布时声波的衰减谱 细实线是参考曲线, 它对应着存在两种混合气泡的情况; 对其他曲线而言, 各个气泡的静态半径分别是R = 20 + 40(j – 1)/ (M1)(μ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(j1)/(M1)(μm), j=1,2,,M.

    这里试着解释上面的结论, 前面已知小气泡共振特征量的抑制效果远比大气泡的大. 当气泡数分布偏向大气泡时, 小气泡的数量就比较少, 小气泡对应的高频共振衰减被极大地抑制了, 因此衰减谱的带宽比较小. 当气泡数分布偏向小气泡时, 小气泡的数量比较大, 小气泡对应的高频共振衰减的抑制效果较小; 虽然大气泡的数量较小, 但是大气泡的特征量的抑制效果本来就弱, 因此大气泡对应的低频共振衰减的抑制效果也较小. 结果是衰减峰带宽比较大.

    本文运用有效介质理论研究了含混合气泡液体中声波的共振传播特性, 得出的主要结论概括如下.

    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 333Google Scholar

    Ying C F 2008 J. Appl. Acoust. 27 333Google 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 732Google 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 30Google Scholar

    Wang C H, Lin S Y 2008 J. Shaanxi Norm. Univ. 36 30Google Scholar

    [9]

    王勇, 林书玉, 张小丽 2014 物理学报 63 034301Google Scholar

    Wang Y, Lin S Y, Zhang X L 2014 Acta Phys. Sin. 63 034301Google Scholar

    [10]

    陈伟中 2018 应用声学 37 675Google Scholar

    Chen W Z 2018 J. Appl. Acoust. 37 675Google Scholar

    [11]

    朱哲民, 杜功焕 1995 声学学报 6 425

    Zhu Z M, Du H G 1995 Acta Acust. 6 425

    [12]

    An Y 2012 Phys. Rev. E 85 016305Google Scholar

    [13]

    Vanhille C, Campospozuelo C 2013 Ultrason. Sonochem. 20 963Google Scholar

    [14]

    Vanhille C, Campospozuelo C 2012 Ultrason. Sonochem. 19 217Google Scholar

    [15]

    Lebon G S B, Tzanakis I, Djambazov G, et al. 2017 Ultrason. Sonochem. 37 660Google Scholar

    [16]

    Zhang Y, Guo Z, Du X 2018 Appl. Therm. Eng. 133 483Google Scholar

    [17]

    Zhang Y, Du X 2015 Ultrason. Sonochem. 26 119Google Scholar

    [18]

    Trujillo F J 2018 Ultrason. Sonochem. 47 75Google Scholar

    [19]

    张鹏利, 林书玉, 张涛 2013 中国科学: 物理学 力学 天文学 43 249

    Zhang P L, Lin S Y, Zhang T 2013 Sci. Sin. Phys. Mech. Astron. 43 249

    [20]

    苗博雅, 安宇 2015 物理学报 64 225Google Scholar

    Miao B Y, An Y 2015 Acta Phys. Sin. 64 225Google Scholar

    [21]

    王德鑫, 那仁满都拉 2018 物理学报 67 037802Google Scholar

    Wang D X, Naranmandula 2018 Acta Phys. Sin. 67 037802Google 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)

  • 图 1  含两种混合气泡的液体中声波的衰减谱

    Fig. 1.  Attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.

    图 2  含两种混合气泡的液体中声波衰减谱的相对峰值和相对带宽与θ1(或者θ2)的关系

    Fig. 2.  Relationships between θ1(θ2)and relative peak values or bandwidths of attenuation spectrums of acoustic waves in liquids with two kinds of bubbles.

    图 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.

    图 4  在含两种混合气泡的液体中, β0不同时声波的相对衰减谱 图中θ1=0.8, β0分别是106, 104, 103102, 箭头表示该角频率区域β0增加的方向

    Fig. 4.  Attenuation spectrums for different β0 in liquids with two kinds of bubbles. Here, θ1=0.8, β0 is 106, 104, 103 and 102, respectively. The arrow indicates the direction of increase of β0.

    图 5  在含两种混合气泡的液体中, β0不同时声波的声速谱 图中θ1=0.8, β0分别是106, 104, 103102, 箭头表示该角频率区域β0增加的方向

    Fig. 5.  Speed spectrums for different β0 in liquids with two kinds of bubbles. Here, θ1=0.8, β0 is 106, 104, 103 and 102, respectively. The arrow indicates the direction of increase of β0.

    图 6  在两种混合气泡的液体中黏度不同时声波的衰减谱

    Fig. 6.  Attenuation spectrums for different viscosity in liquids with two kinds of bubbles.

    图 7  在含多种混合气泡的液体中气泡的种类数不同时声波的衰减谱. 图中θj=1/M, Rj=20+40(j1)/(M1)(μm)(j=1,2,,M).

    Fig. 7.  Attenuation spectrums for different M. Here, θj=1/M, Rj=20+40(j1)/(M1)(μm), j=1,2,,M.

    图 8  气泡满足不同分布时声波的衰减谱 细实线是参考曲线, 它对应着存在两种混合气泡的情况; 对其他曲线而言, 各个气泡的静态半径分别是R = 20 + 40(j – 1)/ (M1)(μ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(j1)/(M1)(μ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 333Google Scholar

    Ying C F 2008 J. Appl. Acoust. 27 333Google 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 732Google 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 30Google Scholar

    Wang C H, Lin S Y 2008 J. Shaanxi Norm. Univ. 36 30Google Scholar

    [9]

    王勇, 林书玉, 张小丽 2014 物理学报 63 034301Google Scholar

    Wang Y, Lin S Y, Zhang X L 2014 Acta Phys. Sin. 63 034301Google Scholar

    [10]

    陈伟中 2018 应用声学 37 675Google Scholar

    Chen W Z 2018 J. Appl. Acoust. 37 675Google Scholar

    [11]

    朱哲民, 杜功焕 1995 声学学报 6 425

    Zhu Z M, Du H G 1995 Acta Acust. 6 425

    [12]

    An Y 2012 Phys. Rev. E 85 016305Google Scholar

    [13]

    Vanhille C, Campospozuelo C 2013 Ultrason. Sonochem. 20 963Google Scholar

    [14]

    Vanhille C, Campospozuelo C 2012 Ultrason. Sonochem. 19 217Google Scholar

    [15]

    Lebon G S B, Tzanakis I, Djambazov G, et al. 2017 Ultrason. Sonochem. 37 660Google Scholar

    [16]

    Zhang Y, Guo Z, Du X 2018 Appl. Therm. Eng. 133 483Google Scholar

    [17]

    Zhang Y, Du X 2015 Ultrason. Sonochem. 26 119Google Scholar

    [18]

    Trujillo F J 2018 Ultrason. Sonochem. 47 75Google Scholar

    [19]

    张鹏利, 林书玉, 张涛 2013 中国科学: 物理学 力学 天文学 43 249

    Zhang P L, Lin S Y, Zhang T 2013 Sci. Sin. Phys. Mech. Astron. 43 249

    [20]

    苗博雅, 安宇 2015 物理学报 64 225Google Scholar

    Miao B Y, An Y 2015 Acta Phys. Sin. 64 225Google Scholar

    [21]

    王德鑫, 那仁满都拉 2018 物理学报 67 037802Google Scholar

    Wang D X, Naranmandula 2018 Acta Phys. Sin. 67 037802Google Scholar

    [22]

    Keller J B, Kolodner I I 1956 J. Appl. Phys. 2 71152

  • [1] 邱春扬, 沈壮志, 姚博文. 涡流场中气泡破碎对声空化效果的影响. 物理学报, 2025, 74(4): 044301. doi: 10.7498/aps.74.20241523
    [2] 郑雅欣, 那仁满都拉. 可压缩液体中气泡的声空化特性. 物理学报, 2022, 71(1): 014301. doi: 10.7498/aps.71.20211266
    [3] 郑雅欣, 那仁满都拉. 可压缩液体中气泡的声空化特性. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211266
    [4] 清河美, 那仁满都拉. 不同类型气泡组成的混合泡群声空化特性. 物理学报, 2020, 69(18): 184301. doi: 10.7498/aps.69.20200381
    [5] 清河美, 那仁满都拉. 空化多泡中大气泡对小气泡空化效应的影响. 物理学报, 2019, 68(23): 234302. doi: 10.7498/aps.68.20191198
    [6] 王汝佳, 吴士平, 陈伟. 热粘弹波在变温非均匀合金熔体中的传播. 物理学报, 2019, 68(4): 048101. doi: 10.7498/aps.68.20181923
    [7] 李佳蔚, 鹿力成, 郭圣明, 马力. warping变换提取单模态反演海底衰减系数. 物理学报, 2017, 66(20): 204301. doi: 10.7498/aps.66.204301
    [8] 张攀, 赵雪丹, 张国华, 张祺, 孙其诚, 侯志坚, 董军军. 垂直载荷下颗粒物质的声波探测和非线性响应. 物理学报, 2016, 65(2): 024501. doi: 10.7498/aps.65.024501
    [9] 史久林, 郭鹏峰, 黄育, 钱佳成, 王泓鹏, 刘娟, 何兴道. 温度、湿度及压强对激光在水中衰减特性的影响. 物理学报, 2015, 64(2): 024215. doi: 10.7498/aps.64.024215
    [10] 沈壮志. 声驻波场中空化泡的动力学特性. 物理学报, 2015, 64(12): 124702. doi: 10.7498/aps.64.124702
    [11] 苗博雅, 安宇. 两种气泡混合的声空化. 物理学报, 2015, 64(20): 204301. doi: 10.7498/aps.64.204301
    [12] 王勇, 林书玉, 张小丽. 含气泡液体中的非线性声传播. 物理学报, 2014, 63(3): 034301. doi: 10.7498/aps.63.034301
    [13] 王勇, 林书玉, 张小丽. 声波在含气泡液体中的线性传播. 物理学报, 2013, 62(6): 064304. doi: 10.7498/aps.62.064304
    [14] 张军, 曾新吾, 陈聃, 张振福. 水下强声波脉冲负压的产生和空化气泡运动. 物理学报, 2012, 61(18): 184302. doi: 10.7498/aps.61.184302
    [15] 沈壮志, 林书玉. 声场中气泡运动的混沌特性. 物理学报, 2011, 60(10): 104302. doi: 10.7498/aps.60.104302
    [16] 何兴道, 夏健, 史久林, 刘娟, 李淑静, 刘建安, 方伟. 水的衰减系数及有效增益长度对受激布里渊散射输出能量的影响. 物理学报, 2011, 60(5): 054207. doi: 10.7498/aps.60.054207
    [17] 刘 娟, 白建辉, 倪 恺, 景红梅, 何兴道, 刘大禾. 受激布里渊散射对激光在水中衰减特性的影响. 物理学报, 2008, 57(1): 260-264. doi: 10.7498/aps.57.260
    [18] 刘海军, 安宇. 空化单气泡外围压强分布. 物理学报, 2004, 53(5): 1406-1412. doi: 10.7498/aps.53.1406
    [19] 杜启振, 杨慧珠. 线性黏弹性各向异性介质速度频散和衰减特征研究. 物理学报, 2002, 51(9): 2101-2108. doi: 10.7498/aps.51.2101
    [20] 王克斌, 李士, 唐孝威. 应用共振吸收谱仪测量Al对γ射线的衰减系数. 物理学报, 1981, 30(9): 1279-1283. doi: 10.7498/aps.30.1279
  • 期刊类型引用(2)

    1. 赵宁宁,肖新宇,凡凤仙,苏明旭. 基于蒙特卡罗原理的混合颗粒三相体系声衰减计算模型研究. 物理学报. 2022(07): 217-224 . 百度学术
    2. 清河美,那仁满都拉. 空化多泡中大气泡对小气泡空化效应的影响. 物理学报. 2019(23): 167-175 . 百度学术

    其他类型引用(6)

计量
  • 文章访问数:  8043
  • PDF下载量:  73
  • 被引次数: 8
出版历程
  • 收稿日期:  2018-12-29
  • 修回日期:  2019-01-28
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-04-05

/

返回文章
返回