馬艷1)2) 林書玉1)? 徐潔1) 唐一璠1)

1)(陜西師范大學陜西省超聲學重點實驗室,西安 710062)2)(寧夏師范學院物理與電子信息工程學院,納米結構及功能材料工程技術研究中心,固原 756000)(2016年7月5日收到;2016年10月10日收到修改稿)

非球形效應對強聲場中次Bjerknes力的影響?

馬艷1)2) 林書玉1)? 徐潔1) 唐一璠1)

1)(陜西師范大學陜西省超聲學重點實驗室,西安 710062)2)(寧夏師范學院物理與電子信息工程學院,納米結構及功能材料工程技術研究中心,固原 756000)(2016年7月5日收到;2016年10月10日收到修改稿)

考慮了非球形氣泡在聲場中的形狀振動,推導了非球形氣泡和球形氣泡之間的次Bjerknes力方程,數值模擬了聲場中非球形氣泡和球形氣泡之間的次Bjerknes力和兩個球形氣泡之間的次Bjerknes力,并對非球形氣泡和球形氣泡之間的次Bjerknes力的影響因素進行了分析討論.研究結果表明:當驅動聲壓振幅大于非球形氣泡的Black閾值且又能使得非球形氣泡穩(wěn)定振動時,在第一個聲驅動周期內,非球形氣泡和球形氣泡之間的次Bjerknes力和兩個球形氣泡的次Bjerknes力方向差異較大,在大小上是兩個球形氣泡次Bjerkens力的數倍,且有著更長的作用距離.非球形氣泡和球形氣泡之間的次Bjerknes力取決于非球形氣泡的形狀模態(tài)、兩個氣泡初始半徑的比值、驅動聲壓振幅、氣泡間距和兩個氣泡的相對位置.

非球形氣泡,次Bjerknes力,形狀模態(tài)

1引 言

聲場中氣泡之間的相互作用是一個著名的聲學現(xiàn)象,會對高壓強場中氣泡的聚合產生重要影響,例如流光的形成和多泡聲致發(fā)光(MBSL)等[1?5],在含氣泡液體中,氣泡受到聲場的作用力主要涉及Bjerknes力.Bjerknes力的典型類型有兩種:一種是單個氣泡在聲場中聲壓波腹處受到的吸引力或排斥力,稱為主Bjerknes力,是聲場中氣泡所受到的主要聲場力;第二種是由于驅動聲場對氣泡作用引起氣泡振動,而使得相鄰氣泡之間產生的相互吸引或排斥力,稱為次Bjerknes力或氣泡之間的相互作用力.自從Bjerknes力被發(fā)現(xiàn)之后,許多作者已經對其進行了實驗和理論研究[6?14].Arston等用視頻記錄了超聲駐波引起的氣泡的聚集,Crum[10]給出了兩個球形氣泡之間的次Bjerknes力的簡單模型,提出了兩個球形氣泡之間的次Bjerknes力的方向僅取決于兩個球形氣泡是否在相位內振動,如果驅動頻率位于兩個氣泡的共振頻率之間,兩個氣泡互相排斥,反之則吸引的經典理論.在一段時間內經典次Bjerknes力理論曾是人們解釋氣泡之間相互作用的主要理論,但近年來發(fā)現(xiàn),氣泡之間的作用力是一個復雜的動力學現(xiàn)象.Yoshida和Takaak[11]用高速攝影機連接長焦距的顯微鏡觀測了兩個氣泡在次Bjerknes力作用下的運動,特別是次Bjerknes力方向的反轉,并與Crum經典次Bjerknes力理論做了比較,發(fā)現(xiàn)實驗現(xiàn)象并不能用經典的次Bjerknes力理論解釋.

隨后,多位學者對經典次Bjerknes力理論做了修正,使得氣泡之間的次Bjerknes力更接近液體中氣泡的真實情況.例如:Zabolotskaya[12]和Ida[13]通過線性模型和兩個球形氣泡的耦合振動理論證明次Bjerknes力的方向取決于氣泡之間的距離.盡管如此,實際聲場中,氣泡的振動和氣泡之間的作用力是非常復雜的聲學現(xiàn)象,有很多影響因素,而上述理論中都采用的是線性理論和球形氣泡模型.隨后Mettin等[14]利用考慮了液體可壓縮性的非線性Keller-Miksis模型方程研究了強聲場中球對稱小振蕩空化泡(R＜10μm)間的作用力,結果表明聲場強度和方向都會引起次Bjerknes力的變化,他們的研究結果在一定程度上解釋了一些經典次Bjerknes力理論所不能解釋的實驗現(xiàn)象,是次Bjerknes力研究中的一項重要成果.但是Mettin等的研究仍然建立在球形氣泡的基礎上,而聲致發(fā)光實驗表明:氣泡在振蕩時,已經不再是標準的球形[15?20],而當兩個氣泡之間的距離較近時,氣泡已經不能完全近似為球形氣泡,非球形效應對次Bjerknes力的影響將變得重要[14].強聲場中(Pdr＞1.0×105Pa),較近氣泡間距時氣泡的形變對氣泡之間的次Bjerkens力產生什么樣的影響,以及次Bjerkens力如何影響氣泡的振動和運動,是文本研究的重點內容.由于兩個非球形氣泡之間的次Bjerknes力非常復雜,兩個非球形氣泡的形狀模態(tài)的耦合與否均會對次Bjerknes力產生影響,因此本文只研究了一種特殊情況:一個氣泡為非球形氣泡,另一個氣泡近似為球形氣泡,以此來修正球形氣泡之間的次Bjerknes力理論,并研究強聲場中(Pdr＞1.0×105Pa)較近氣泡間距下,非球形效應對次Bjerknes力的影響以及試圖用研究結果來解釋強聲場中泡群的聚合現(xiàn)象例如“streamer formation”和MBSL等現(xiàn)象.

2非球形氣泡和球形氣泡的次Bjerknes力

對于聲場中兩個球形氣泡,由于氣泡在聲場中做周期性振動,而使得相鄰氣泡之間產生相互作用力,兩個球形氣泡之間的次Bjerknes力可表示為如下形式[14]:

(1)式中FB為兩個球形氣泡之間的次Bjerknes力,為液體密度;d為兩個氣泡之間的距離;為兩個球形氣泡的體積對時間的一階微分;〈〉表示在一個驅動周期內取平均值,令cB=則兩個球形氣泡之間的次Bjerknes力的方向取決于〈cB〉的正負,其中〈cB〉大于0為吸引力,〈cB〉小于0為排斥力.

當聲場中有兩個氣泡,假設氣泡1為非球形氣泡,氣泡2近似認為是球形氣泡(見圖1),非球形氣泡壁距氣泡中心的距離可表示為[20]

其中R1(t)是非球形氣泡1無擾動時的半徑,Yn(cosθ)為n階Legendre函數,an為非球形氣泡1的n階形狀模態(tài)振幅,且認為非球形氣泡為單模態(tài),不存在不同形狀模態(tài)的耦合.根據Eller類似的推導,非球形氣泡1對球形氣泡2的作用力為

(3)式中V2(t)為氣泡2的體積,?P1為非球形氣泡1的聲壓梯度,根據液體運動方程:

u1(t)為由氣泡1振動而引起的液體的振動速度,可表示為

其中φ1為由于非球形氣泡1振動而在液體中某一位置處產生的速度勢,可表示為[20]

聯(lián)立(5)和(6)式并對時間微分可得

將(4)和(7)式代入(3)式,然后在一個周期內對其取平均可得非球形氣泡和球形氣泡之間的次Bjerknes力為非球形氣泡和球形氣泡之間的次Bjerkne力的方向取決于〈cBNS〉的正負,其中,〈cBNS〉大于0為吸引力,〈cBNS〉小于0為排斥力.當氣泡2為球形氣泡時,an=0,˙an=0,¨an=0.(8)式經微分變形,便回到兩個球形氣泡之間的次Bjerkne力公式,即(1)式.

圖1 聲場中的非球形氣泡和球形氣泡Fig.1.A nonspherical bubble and a spherical bubble in an acoustic fi eld.

當氣泡1和2在聲場中振動時,氣泡振動方程為考慮了氣泡之間相互作用的Keller-Miksis方程[9]:

(9a)和(10a)式中R2為氣泡2的半徑,點代表對時間求導數;c為液體中的聲速.其中:

(9b)和(10b)中的Pb1和Pb2分別表示氣泡1和2的內部壓強,ν為液體的黏滯系數,σ為表面張力系數,非球形氣泡1的n階形狀模態(tài)振幅an可表示為[21]

其中:

式中δ為氣泡1壁周圍黏滯引起的渦流層厚度,可表示為

3非球形氣泡和球形氣泡之間的次Bjerknes力的影響因素

當氣泡的初始形狀為非球形,在聲場的作用下,有兩種運動趨勢:要么隨時間的變化趨于球形,要么隨時間變化越來越偏離球形,最終破裂[16?20].對于第二種情況下的次Bjerknes力不屬于本文討論的范圍,在這里我們只討論第一種情況下的次Bjerknes力及其影響因素,即第一個振動周期內非球形氣泡的時均次Bjerknes力及其影響因素.經過1個周期或者幾個周期后,非球形氣泡的形變隨時間而趨于零,非球形氣泡和球形氣泡之間的相互作用力將回到兩個球形氣泡之間的相互作用力.為了研究這種由非球形效應引起的次Bjerknes力的變化,我們對球形氣泡和非球形氣泡的次Bjerknes力進行數值模擬,并試圖從模擬結果來探討非球形效應對次Bjerknes力的影響機理.

在以下的數值模擬中,模擬條件為:驅動頻率f=20kHz,液體靜壓強P0=1.01×105Pa,表面張力系數σ=7.25×10?2N/m,水中的聲速為:c=1490m/s,驅動聲壓振幅為Pdr=1.34×105Pa.

圖2是兩個氣泡相對位置不同時,非球形氣泡1(2μm)和球形氣泡2(4μm)在氣泡間距為200μm時,兩個氣泡之間次Bjerknes力的變化關系曲線.數值模擬結果表明:對于非球形氣泡和球形氣泡之間的次Bjerknes力而言,其大小與兩個氣泡的相對位置有關.在上述的模擬條件下,當非球形氣泡的長軸和兩個氣泡中心連線的角度不同,兩個氣泡之間的次Bjerknes力也不相同,但是相差不大,不同位置處的最大次Bjeknes力為?0.08041μN,最小次Bjerknes力為?0.07993μN,均為吸引力.同樣條件下兩個球形氣泡間的次Bjerknes力為?0.0296μN,也是吸引力.因此,在上述模擬條件下,氣泡的非球形效應使得氣泡之間的次Bjerknes力較球形氣泡間的次Bjerknes力的大小增加了2.7003—2.7166倍.以下的討論中,我們只討論非球形氣泡長軸平行兩個氣泡中心連線的情況.

圖2 非球形氣泡和球形氣泡之間次Bjerknes力隨兩個氣泡相對位置變化關系曲線(R10=2μm,R20=4μm,d=200μm)Fig.2.The secondary Bjerknes force between a nonspherical bubble and a spherical bubble vs relative position of two bubbles(R10=2 μm,R20=4 μm and d=200μm).

3.1 非球形氣泡形狀模態(tài)對次Bjerknes力的影響

圖3 氣泡的徑向振動和氣泡1的2,3階形狀模態(tài)振幅隨時間的變化(a)氣泡1,2的徑向振動隨時間的變化;(b)氣泡1的2階形狀模態(tài)隨時間的變化;(c)氣泡1的3階形狀模態(tài)隨時間的變化Fig.3.Radial oscillations of two bubbles and the amplitudes of di ff erent shape modes of the fi rst bubblevs time during onedrivingperiod:(a)The radial oscillations of two bubbles;(b)the amplitude of the second-order shape mode of the fi rst bubble;(c)the amplitude of the third-order shape mode of the fi rst bubble.

當非球形氣泡的初始半徑為2μm,球形氣泡的初始半徑為4μm,氣泡間距為100μm時非球形氣泡的2和3階形狀模態(tài)振幅如圖3(b)和圖3(c)所示,可以看出,在驅動聲壓振幅為1.34×105Pa,驅動頻率為20kHz的驅動聲場下,非球形氣泡的2和3階形狀模態(tài)振幅只在第一個振動周期內隨時間變化,當非球形氣泡的振動時間超過一個周期后,其2和3階形狀模態(tài)振幅趨于零,其振動恢復到球形氣泡的振動狀態(tài).

圖4(a)是上述驅動聲場下,相距100μm的兩個球形氣泡在一個驅動周期內的cB隨時間的變化關系,圖4(b)是同樣驅動條件下相距100μm的非球形氣泡和球形氣泡一個驅動周期內的cBNS隨時間的變化關系(n=2),cB和cBNS反映了兩個氣泡之間的次Bjerknes力的方向和大小在一個周期內隨時間的變化關系.對比圖4(a)和圖4(b)可以看出:在一個驅動周期內,兩個球形氣泡之間的相互作用力隨時間的變化關系與同樣條件下非球形氣泡和球形氣泡之間相互作用力隨時間的變化關系差別較大.通過計算可以得到:兩個球形氣泡之間的〈cB〉=?9.1498×10?20m6/s2,為排斥力;而非球形氣泡對球形氣泡的〈cBNS〉=1.2940×10?19m6/s2,為吸引力.也就是說在上述的驅動條件和氣泡間距下,非球形氣泡和球形氣泡之間的次Bjerknes力的大小是同樣條件兩個球形氣泡之間的次Bjerknes力的1.414倍,且方向發(fā)生改變.

圖4 兩個氣泡之間cB隨時間的變化 (a)球形氣泡之間cB隨時間的變化;(b)非球形氣泡和球形氣泡之間的cBNS隨時間的變化(n=2);(c)非球形氣泡和球形氣泡之間的cBNS隨時間的變化(n=3)Fig.4.The curve of cBbetween two bubbles vs time during one driving period:(a)cBbetween two spherical bubbles;(b)cBNSbetween a nonspherical bubble and a spherical bubble(n=2);(c)cBNSbetween a nonspherical bubble and a spherical bubble(n=3).

同樣的驅動聲場下,當非球形氣泡的形狀模態(tài)階數為n=3時,數值模擬結果表明(圖4(c)):非球形氣泡和球形氣泡的〈cBNS〉=3.2678×10?19m6/s2,為吸引力.非球形氣泡和球形氣泡之間的次Bjerknes力的大小是同樣驅動聲場中兩個球形氣泡之間的次Bjerknes力的3.571倍,且方向發(fā)生改變.

綜上所述:氣泡的非球形效應會對氣泡的次Bjerknes力產生較大影響,能夠在一定的條件下使得球形氣泡間的次Bjerknes力改變方向,并且增大數倍.具體原因可能是:非球形氣泡在聲場中振動,除了徑向振動外還有形狀振動,對于我們的研究情況而言,雖然形狀振動時間非常短,往往在一個或者幾個周期內迅速衰減為零,但這種徑向振動和形狀振動的耦合及形變結束時形變能量轉移為振動能量的結果使得一個周期內非球形氣泡的輻射聲場迅速增加,使之不同于球形氣泡徑向振動所形成的輻射聲場,這種變化使得非球形氣泡和球形氣泡之間的次Bjerknes力迅速增大,在一定條件下氣泡的非球形效應也可能使得氣泡之間的次Bjerknes力方向發(fā)生改變.而且在上述研究條件下非球形氣泡形狀模態(tài)階數越大,其形狀振動與徑向振動耦合產生的輻射聲壓梯度越大,對周圍氣泡產生更強的作用力.

由此可以看出,當聲場中的驅動聲壓振幅大于氣泡的Black空化閾值,非球形氣泡長軸與氣泡間距平行時,在較近的氣泡間距下,非球形氣泡和球形氣泡之間的次Bjerknes力和同樣初始尺寸的兩個球形氣泡之間的次Bjerknes力的差異較大,與Mettin理論中預測的強聲場中(＞1×105Pa)兩個球形氣泡之間次Bjerknes力的變化有一定差異,而這些差異也可能是強聲場中形成的“stream formation”和MBSL的原因.

3.2 氣泡間距對球形氣泡和非球形氣泡的次Bjerknes力的影響

兩個氣泡之間的距離是影響氣泡次Bjerknes力的重要因素,兩個球形氣泡的次Bjerknes力會隨著氣泡間距的增大而逐漸減小,當兩個氣泡之間的距離大于一定值時,兩個球形氣泡就會成為兩個無耦合的自由氣泡,此時兩個球形氣泡之間的次Bjerknes力可忽略不計.為了研究非球形效應對氣泡之間的次Bjerknes力的影響因素,我們研究了不同氣泡間距下非球形氣泡和球形氣泡之間的次Bjerknes力,并與同樣條件下的兩個球形氣泡之間的次Bjerknes力做了對比.在以下的討論中,我們只考慮非球形氣泡做2階形狀模態(tài)振動的情形.

圖5是不同氣泡間距下的兩個氣泡的次Bjerknes力的對數隨氣泡1的初始半徑變化曲線,圖5(a)—(d)中黑色虛線代表非球形氣泡和球形氣泡的次Bjerknes力對數的變化,紅色實線代表了同樣條件下的兩個球形氣泡的次Bjerknes力對數的變化,對比可以發(fā)現(xiàn),當非球形氣泡的初始半徑保持不變(R20=2μm),球形氣泡2的初始半徑從0.5μm變化到10μm的過程中,當兩個氣泡的初始半徑為一定比例時,非球形氣泡和球形氣泡之間的次Bjerknes力是同樣條件下兩個球形氣泡之間的次Bjerknes力的幾倍到十幾倍.而除此之外,非球形氣泡和球形氣泡之間的次Bjerknes力近似等于同樣條件下的兩個球形氣泡之間的次Bjerknes力.當兩個氣泡之間的距離增大到500μm和1000μm時,通過數值模擬可以看出,兩個球形氣泡之間的次Bjerknes力非常小,趨于0,而同樣尺寸的非球形氣泡和球形氣泡之間的次Bjerknes力也隨著氣泡間距的增加出現(xiàn)了衰減,但相較同樣條件下的兩個球形氣泡之間的次Bjerknes力,仍然保持一個較大的數量級,如圖5(c)—(d).這說明:非球形氣泡和球形氣泡之間的次Bjerknes力有更大的作用距離,能夠使得相鄰的氣泡在較遠的距離處仍然受到相當大小的次Bjerknes力作用,該結果與文獻[11]中的實驗結果相一致.

數值模擬結果表明:非球形氣泡和球形氣泡之間的次Bjerknes力,隨氣泡間距的變化關系不同于同樣尺寸下兩個球形氣泡之間的次Bjerknes力隨氣泡間距的變化關系,這是由于非球形氣泡的形狀振動,使得不同的氣泡間距下,不同初始半徑的球形氣泡和非球形氣泡對的耦合模式完全不同,雖然會出現(xiàn)次Bjerkens力隨氣泡間距減小的整體趨勢,但仍然會出現(xiàn)一些特殊比例的氣泡對的次Bjerknes力在小氣泡間距下較小或者在大氣泡間距下較大(圖5(a)—(d)).除此之外,數值模擬結果表明非球形氣泡和球形氣泡之間的次Bjerknes力的作用距離更長,當氣泡距離大于1000μm以上才逐漸衰減為0.

圖5 (網刊彩色)不同氣泡間距下兩個氣泡之間的次Bjerknes力的對數隨氣泡2初始半徑的變化曲線(Pdr=1.34×105Pa,R10=2μm) (a)d=100μm;(b)d=200μm;(c)d=500μm;(d)d=1000μmFig.5.(color online)The curve of the logarithm of the secondary Bjerknes force between two bubbles vs the radius of the second bubble(Pdr=1.34×105Pa,R10=2μm):(a)d=100μm;(b)d=200μm;(c)d=500μm;(d)d=1000μm.

3.3 驅動聲壓振幅對球形氣泡和非球形氣泡之間的次Bjerknes力的影響

圖6是不同驅動聲壓振幅下的非球形氣泡和球形氣泡之間的次Bjerknes力的對比圖,其中圖6(a)是兩個氣泡之間距離為200μm時,兩個球形氣泡的次Bjerknes力隨驅動壓強振幅的變化關系.圖6(b)是兩個氣泡之間距離為200μm時,非球形氣泡和球形氣泡之間的次Bjerknes力隨驅動壓強的變化關系,對比數據可以發(fā)現(xiàn),兩個球形氣泡之間的次Bjerknes力和非球形氣泡和球形氣泡之間的次Bjerknes力的大小都隨著驅動聲壓振幅的增加而增加.當驅動聲壓振幅為1.30×105Pa和1.32×105Pa時,兩個球形氣泡之間的次Bjerknes力和同樣氣泡間距、同樣初始尺寸的非球形和球形氣泡之間的次Bjerknes力的大小相差不多.當驅動聲壓振幅增大到1.34×105Pa時,數值模擬結果表明,非球形氣泡和球形氣泡之間的次Bjerknes力是同樣條件下兩個球形氣泡之間的次Bjerknes力的幾倍甚至十幾倍以上.

綜上所述：當驅動聲壓振幅使得非球形氣泡能夠穩(wěn)定振動而不破裂時,驅動聲壓振幅越大,非球形效應對次Bjerknes力的影響越大,特別是當驅動聲壓振幅超過兩個氣泡的Black閾值后,使得非球形氣泡和球形氣泡之間的次Bjerknes力的大小迅速增大,遠大于兩個球形氣泡之間的Bjerknes力.非球形效應隨聲壓振幅增加的原因是,驅動聲壓越強,非球形氣泡的徑向振動幅度和形狀振動幅度也越大,特別是當驅動聲壓振幅增加到非球形氣泡的空化閾值之上,其在一個驅動周期內徑向振動和形狀振動耦合形成的平均輻射聲場梯度將大幅度增加,此時,非球形氣泡和球形氣泡之間的次Bjerknes力也迅速增大.因此,驅動聲壓振幅越大,非球形效應對次Bjerknes力的影響越顯著.

圖6 (網刊彩色)不同驅動聲壓振幅下兩個氣泡之間的次Bjerknes力隨氣泡2初始半徑的變化(R10=2μm,d=200μm)(a)兩個球形氣泡之間的次Bjerknes力;(b)非球形氣泡和球形氣泡的次Bjerknes力Fig.6.(color online)The curve of the secondary Bjerknes force between two bubblesvsthe radius of the second bubble(R10=2 μm,d=200 μm):(a)The secondary Bjerknes force between two spherical bubblest;(b)the secondary Bjerknes force between anonspherical bubble and a pherical bubble.

4結 論

本文在Crum和Mettin等關于球形氣泡的次Bjerknes力的研究基礎上,考慮了聲場中非球形氣泡的形狀振動對次Bjerknes力的影響,推導了非球形氣泡和球形氣泡之間的次Bjerknes力方程,數值模擬了非球形氣泡和球形氣泡之間的次Bjerknes力并與兩個球形氣泡之間的次Bjerknes力進行了對比.

研究發(fā)現(xiàn):當聲場驅動聲壓大于非球形氣泡的Black閾值時且使得非球形氣泡能夠穩(wěn)定振動時,非球形氣泡和球形氣泡之間的次Bjerknes力在大小和方向上都不同于兩個球形氣泡之間的次Bjerknes力,作用力在數量級上大于同樣條件下的兩個球形氣泡之間的次Bjerknes力;非球形效應對次Bjerknes力的影響取決于非球形氣泡的形狀模態(tài)、兩個氣泡初始半徑的比值、驅動聲壓振幅、氣泡間距和兩個氣泡的相對位置.

在以往的研究中,線性條件下兩個球形氣泡之間的次Bjerknes力非常小,人們往往忽略其影響,對于氣泡聚合等實驗現(xiàn)象不能給出合理的解釋.我們的研究發(fā)現(xiàn),在非線性條件下,單個氣泡的非球形效應使得次Bjerknes力較Mettin非線性理論中球形氣泡間的次Bjerknes力在方向和數量級上有一定的差異,能夠使一定尺寸的氣泡對的次Bjerknes力增加幾倍到十幾倍.研究結果雖然只反映了第一個驅動周期次Bjerknes力的變化,但由于在第一個振動周期內非球形效應的影響,已經使得氣泡的運動狀態(tài)發(fā)生了較大變化,進而影響了氣泡隨后的運動.非球形效應使得兩個氣泡之間具有相當大小的作用力,使得氣泡聚合或者分離現(xiàn)象能夠發(fā)生,產生一些穩(wěn)定的氣泡結構,也能夠改變快速運動的小氣泡的運動趨勢,這可能也是“stream formation”的主要原因;對于強聲場中的氣泡群,較大的次Bjerknes力能夠使氣泡的空化效應更容易發(fā)生,這可能也是實驗中多泡聲致發(fā)光更易發(fā)生的一個重要原因.兩個氣泡都發(fā)生形變及氣泡形變耦合對次Bjerknes力的影響將是我們下一步的工作.

[1]Anthony H,Kaper T 2001 J.Fluid Mech.445 377

[2]Thomas J M,Sean M C 1997 J.Acoust.Soc.Am.102 1522

[3]Rossello J M,Dellavale D,Bonetto F J 2015 Ultrason.Sonochem.22 59

[4]Yuan L,Joseph K 2013 Phys.Fluids 25 073301

[5]Eller A 1968 J.Acoust.Soc.Am.43 107

[6]Alexander A D 1997 J.Acoust.Soc.Am 102 747

[7]David R,Pierre T B 2011 Phys.Fluids 23 042003

[8]Mohammad A A 2011 J.Acoust.Soc.Am.130 3321

[9]Rasoul S B,Nastaran R 2010 Phys.Rev.E 82 016316

[10]Crum L A 1975 J.Acoust.Soc.Am.57 1363

[11]Yoshida K J,Takaaki F 2011 J.Acoust.Soc.Am 130 135

[12]Zabolotskaya E A 1984 Sov.Phys.Scoust 30 365

[13]Ida M 2007 Phys.Rev.E 76 04309

[14]Mettin R,Akhatov I,Parlitz U 1997 Phys.Rev.E 56 2924

[15]Shao W H,Chen W Z 2013 J.Acoust.Soc.Am.133 119

[16]Prosperetti A 1977 Q.Appl.Math 34 339

[17]Bogoyavlenskiy V A 2000 Phy.Rev.E 62 2158

[18]Pelekasis N A,Tsamopouslos J A 1990 Phys.Fluids A 2 1328

[19]Xie C G,An Y 2003 Acta Phys.Sin.52 102(in Chinese)[謝崇國,安宇2003物理學報52 102]

[20]Plesset M S 1954 J.Appl.Phys.25 96

[21]Brenner M P,Lohse D,Dupon T F 1995 Phys.Rev.Lett.75 954

PACS:43.25.+y,43.35.+d,47.55.ddDOI:10.7498/aps.66.014302

*Project supported by the National Natural Science Foundation of China(Grant Nos.11374200,11674206).

?Corresponding author.E-mail:sylin@snnu.edu.cn

In fl unece of nonspherical e ff ects on the secondary Bjerknes force in a strong acoustic fi eld?

Ma Yan1)2)Lin Shu-Yu1)?Xu Jie1)Tang Yi-Fan1)

1)(Shaanxi Key Laboratory of Ultrasonics,Shaanxi Normal University,Xi’an 710062,China)2)(College of Physics and Electronic Information Engineering,Engineering Research Center of Nanostructure and Functional Materials,Ningxia Normal University,Guyuan 756000,China)(Received 5 July 2016;revised manuscript received 10 October 2016)

The secondary Bjerknes force between bubbles in an acoustic fi eld is a well-known acoustic phenomenon.The theoretical researches of the secondary Bjerknes force mainly focus on the case of two spherical bubbles.The secondary Bjerknes force between two spherical bubbles,calculated based on the linear equations,is very small and negligible.Therefore these theoretical researches donot give a good explanation for the phenomenon,such as“streamer formation”and multi-bubble sonoluminescence(MBSL).Experiments of sonoluminescence show that the shapes of the bubbles in a sound fi eld are not entirely spherical.Nonspherical e ff ects have an important in fl uence on the secondary Bjerknes force when two bubbles come close to each other in a strong acoustic fi eld(＞ 1.0× 105Pa).How the shape distortion of a nonspherical bubble causes the secondary Bjerknes force between two bubbles to change,and how the secondary Bjerknes force a ff ects the oscillations and movements of bubbles are major problems which we are to solve in the present research.The expression of the secondary Bjerknes force between a nonspherical bubble and a spherical bubble is obtained by considering the shape oscillation of a nonspherical bubble.We numerical simulate the secondary Bjerknes force between a nonspherical bubble and a spherical bubble based on the nonlinear oscillation equations of two bubbles,and compare the secondary Bjerknes force between a nonspherical bubble and a spherical bubble with the secondary Bjerknes force between two spherical bubbles in the same condition.We discuss the in fl uence of nonspherical e ff ects on the secondary Bjerknes force between two bubbles.The results show that when the amplitude of driving pressure is greater than the Blake threshold of a nonspherical bubble and makes the bubble oscillate stably,the secondary Bjerknes force between this nonspherical bubble and a spherical bubble is di ff erent from the secondary Bjerknes force between two spherical bubbles in direction and magnitude.The secondary Bjerknes force between a nonspherical bubble and a spherical bubble is much bigger than that between two spherical bubbles.The interactional distance of the secondary Bjerknes force between a nonspherical bubble and a spherical bubble is longer than that between two spherical bubbles.The secondary Bjerknes force between a spherical bubble and a nonspherical bubble depends on the radii of two bubbles,distance between two bubbles,shape mode of the nonspherical bubble and the amplitude of driving pressure.Our research is closer to the actual bubbles in liquid.We also prove that big mutual interaction between bubbles is the main cause for froming a stable structure between bubbles.For bubbles,big mutual interaction causing the cavitation becomes easier.These results are important for explaining the phenomenon in an acoustic fi eld,such as“streamer formation”and MBSL.

nonspherical bubbles,secondary Bjerknes force,shape mode

10.7498/aps.66.014302

?國家自然科學基金(批準號:11374200,11674206)資助的課題.

?通信作者.E-mail:sylin@snnu.edu.cn