張?zhí)杖?莫潤(rùn)陽(yáng) 胡靜 陳時(shí) 王成會(huì) 郭建中

(陜西師范大學(xué)超聲學(xué)重點(diǎn)實(shí)驗(yàn)室, 西安 710062)

為探究外層彈性介質(zhì)對(duì)液體腔內(nèi)氣泡和粒子相互作用的影響, 從速度勢(shì)分布理論出發(fā), 結(jié)合拉格朗日方程得到了腔內(nèi)氣泡和粒子的運(yùn)動(dòng)方程, 分析了氣泡的共振頻率及聲波作用下粒子和氣泡間的相互作用對(duì)其平動(dòng)行為的影響. 結(jié)果表明, 介質(zhì)彈性和密度等特性均可改變腔內(nèi)氣泡的共振頻率, 表現(xiàn)為氣泡的共振頻率隨著球腔的半徑增大有先減小后增大的變化趨勢(shì), 且逐漸趨于無(wú)界液體中單泡的共振頻率. 氣泡和粒子在球形液體腔中的平動(dòng)受到聲場(chǎng)參數(shù)、外層彈性介質(zhì)特性、氣泡和粒子本身特征等因素的影響, 總體特征表現(xiàn)為粒子有向腔壁運(yùn)動(dòng)的趨勢(shì), 氣泡的平動(dòng)與氣泡和粒子間相互作用特征密切相關(guān).

1 引 言

彈性介質(zhì)包圍的液體腔內(nèi)部氣泡動(dòng)力學(xué)研究在地質(zhì)[1]、生物學(xué)[2,3]等領(lǐng)域有比較廣的應(yīng)用. 氣泡在液體中的非線性振動(dòng)可引起空化效應(yīng), 形成局部高溫高壓和沖擊波等極端現(xiàn)象, 在一定程度上可促進(jìn)化學(xué)反應(yīng)和物質(zhì)交換等. 近年來(lái), 許多學(xué)者對(duì)超聲波作用下氣泡的空化現(xiàn)象進(jìn)行了廣泛而深入的研究, 發(fā)展了氣泡處于無(wú)界液體中[4,5]或液體內(nèi)的界面附近[6]等情形下的氣泡動(dòng)力學(xué). 為簡(jiǎn)化分析,在研究螺旋槳面[7,8]、血管壁[9,10]附近的氣泡振動(dòng)時(shí)近似將界面簡(jiǎn)化為無(wú)限大硬壁面后, 借助于鏡像理論等研究氣泡與界面間的相互作用. 在生物學(xué)中, 液體中微泡運(yùn)動(dòng)通常受到血管或軟介質(zhì)的約束, 因此, 血管內(nèi)氣泡的動(dòng)力學(xué)行為的研究對(duì)分析空化氣泡的生物學(xué)影響尤為重要. 在特定的環(huán)境中, 液體可能存在于封閉腔體內(nèi), 當(dāng)液體空化時(shí),空化氣泡的振動(dòng)行為將呈現(xiàn)出不一樣的特征.Wang[11]和Doinikov[12]等的理論研究表明彈性介質(zhì)中一定大小的球形液體腔內(nèi)氣泡的共振頻率隨著氣泡半徑的增大而增大, 與在無(wú)界液體中氣泡振動(dòng)的自然頻率的變化趨勢(shì)完全相反, Vincent 等[13]通過(guò)實(shí)驗(yàn)得出彈性介質(zhì)中氣泡由于約束作用呈現(xiàn)快速振蕩的運(yùn)動(dòng)趨勢(shì). 這預(yù)示著球形液體腔周圍介質(zhì)的特性將對(duì)氣泡的動(dòng)力學(xué)行為形成顯著影響.

Rayleigh 為解釋螺旋槳的空化腐蝕現(xiàn)象, 最早建立了不可壓縮液體內(nèi)氣泡動(dòng)力學(xué)的理論模型, 定性地分析了空化產(chǎn)生的機(jī)理, 為以后的空化研究奠定了基礎(chǔ). 在此基礎(chǔ)上考慮液體壓縮性、黏性及表面張力后對(duì)其修正得到R-P 方程, Keller 等[14]考慮聲輻射、黏性、表明張力和聲波的影響對(duì)氣泡動(dòng)力學(xué)進(jìn)行了研究. 由于液體中的氣泡大多以多氣泡的形式存在, 人們利用牛頓力學(xué)和分析力學(xué)處理動(dòng)力學(xué)問(wèn)題的手段提出了兩個(gè)和多個(gè)氣泡間相互作用的理論模型, 如Doinikov[15]構(gòu)建了兩氣泡體系的拉格朗日函數(shù)并結(jié)合拉格朗日方程得到了氣泡徑向和平動(dòng)的動(dòng)力學(xué)方程, Ilinskii 等[16]研究了考慮氣泡間相互作用后的耦合聲輻射, 以及通過(guò)構(gòu)建系統(tǒng)哈密頓函數(shù)結(jié)合正則方程得到氣泡徑向振動(dòng)和平動(dòng)的動(dòng)力學(xué)方程. 王成會(huì)等[17]考慮了氣泡間次級(jí)聲輻射的影響, 探討了氣泡的徑向振動(dòng)和平動(dòng)規(guī)律. 為了解釋沖擊波碎石術(shù)中腎結(jié)石碎片和空化團(tuán)簇的作用, 發(fā)展了無(wú)界液體中氣泡和粒子相互作用的模型. Hay 等[18]在模型中假設(shè)氣泡和粒子都可以自由移動(dòng)并精確到五階來(lái)描述氣泡-粒子的相互作用. Li 等[19]利用邊界積分法和實(shí)驗(yàn)研究了氣泡和懸浮球之間的相互作用. 超聲波作用下液體中的氣泡和懸浮顆粒物等均可能與附近的空泡產(chǎn)生相互作用并形成空化影響, 因此, 人們發(fā)展了無(wú)界液體中粒子與氣泡的運(yùn)動(dòng)模型. 然而, 在超聲萃取或者干燥等過(guò)程中, 可能存在生物組織(彈性體)包圍的液體腔內(nèi)的空化行為以及空化泡與粒子間的相互作用. 因此, 有必要發(fā)展受限條件下空化氣泡和粒子間的相互作用的動(dòng)力學(xué)模型, 為解釋超聲促進(jìn)物質(zhì)交換等奠定理論基礎(chǔ). 現(xiàn)有的理論和實(shí)驗(yàn)研究表明[11,13], 負(fù)壓或者聲波作用下, 彈性介質(zhì)包圍的封閉液體腔內(nèi)的氣泡可生長(zhǎng)并振蕩, 其動(dòng)力學(xué)行為將影響周圍介質(zhì)的力學(xué)性能等. 事實(shí)上, 包裹液體腔的介質(zhì)也會(huì)影響腔內(nèi)氣泡的平動(dòng), 而密閉液體腔內(nèi)的多泡體系或者粒子的存在也會(huì)影響氣泡徑向振動(dòng)和平動(dòng)行為, 為探索密閉腔內(nèi)液體中氣泡和粒子間的相互作用, 本文將發(fā)展氣泡-粒子系統(tǒng)的動(dòng)力學(xué)方程, 并分析氣泡和粒子的平動(dòng)特征.

2 理論推導(dǎo)

為探究彈性介質(zhì)包圍的液體腔內(nèi)氣泡和粒子的動(dòng)力學(xué)行為, 從氣泡、粒子及液體腔三者的速度勢(shì)分布函數(shù)出發(fā), 借助坐標(biāo)變換得到速度分布及系統(tǒng)的動(dòng)能和勢(shì)能表達(dá)式, 進(jìn)一步基于拉格朗日方程導(dǎo)出氣泡和粒子的動(dòng)力學(xué)方程.

2.1 速度勢(shì)函數(shù)

如圖1 所示, 氣泡和粒子位于由彈性介質(zhì)包裹的球形液腔內(nèi). 為簡(jiǎn)化分析, 假定二者在運(yùn)動(dòng)過(guò)程中始終保持球形, 且二者球心均位于直徑z軸上,相對(duì)位置坐標(biāo)為zj(j= 1, 2); 以球腔中心為坐標(biāo)原點(diǎn), 建立球坐標(biāo)系(r,θ,φ) ; 分別以氣泡和粒子中心為坐標(biāo)原點(diǎn)建立局部球坐標(biāo)系(rj,θj,φj) . 若腔內(nèi)充滿理想液體, 則其運(yùn)動(dòng)將滿足拉普拉斯方程Δφ=0 , 其中φ為標(biāo)量速度勢(shì).

圖1 充滿液體的空腔中的氣泡和粒子Fig. 1. The bubble and particle in a liquid-filled cavity.

液體速度勢(shì)可表示為

式中,φ1和φ2分別為與氣泡和粒子運(yùn)動(dòng)相關(guān)的液體速度勢(shì),φc為腔內(nèi)壁面運(yùn)動(dòng)相關(guān)的液體速度勢(shì);求解拉普拉斯方程, 有

式中,Rj是為瞬時(shí)氣泡(粒子)半徑,Rc0是靜止時(shí)空腔的半徑,Pn是n階勒讓德多項(xiàng)式,b(nj)(t) 是與時(shí)間相關(guān)的函數(shù),b(nj)可表示為[20]

式中ξj=Rj/d,d(t)=z2(t)-z1(t) 是氣泡 與粒子中心的距離 .

為在氣泡、粒子及球腔表面應(yīng)用邊界條件, 需將速度勢(shì)φj及φc的表達(dá)式進(jìn)行坐標(biāo)變換. 在球腔中心為坐標(biāo)原點(diǎn)的球坐標(biāo)系中, 速度勢(shì)φj為[20]

式中

在以氣泡和粒子中心為坐標(biāo)原點(diǎn)的局部坐標(biāo)系中, 速度勢(shì)φc表示為[20]

式中

氣泡和粒子表面法向速度對(duì)應(yīng)的邊界條件為

得

由于腔體外側(cè)彈性介質(zhì)層僅發(fā)生微小形變, 故腔內(nèi)表面法向速度關(guān)系近似為

式中Rc(t) 是空腔的瞬時(shí)半徑. 將(3)式和(7)式代入(13)式, 得

將(8)式代入(14)式中, 得

由于液體具有不可壓縮性, 故液體體積守恒, 有

基于速度勢(shì)分布以及邊界條件可得到氣泡周圍介質(zhì)中速度分布, 進(jìn)而可以求氣泡徑向振動(dòng)和平動(dòng)過(guò)程中引起的系統(tǒng)動(dòng)能變化, 其推導(dǎo)過(guò)程參見附錄.

2.2 氣泡和粒子的運(yùn)動(dòng)方程

基于拉格朗日方程可以求得徑向振動(dòng)氣泡和粒子平動(dòng)動(dòng)力學(xué)方程. 系統(tǒng)拉格朗日函數(shù)L=T -U,其中T和U是系統(tǒng)的動(dòng)能和勢(shì)能. 系統(tǒng)動(dòng)能包括:腔內(nèi)液體動(dòng)能Tl、剛性粒子平動(dòng)動(dòng)能Tp和彈性介質(zhì)層動(dòng)能Ts, 即T=Tl+Tp+Ts. 由附錄知腔體內(nèi)液體動(dòng)能可近似表示為

粒子的動(dòng)能為[21]

式中,ρl為液體的密度,Vl為液體所占體積,ρp為粒子的密度.

空腔外分布有球殼狀彈性介質(zhì)層, 其殼層內(nèi)外半徑分布為Rc0和Rs. 由于液體不可壓縮, 空腔內(nèi)氣泡的徑向體積脈動(dòng)必然引起彈性殼層的運(yùn)動(dòng), 其動(dòng)能可以表示為[20]

式中ρs為固體的密度. 系統(tǒng)勢(shì)能包括氣泡體積勢(shì)能和腔壁外層介質(zhì)彈性勢(shì)能, 即U=U1+U2, 其中U1可表示為[21]

式中P1(t) 是氣泡表面的壓力,FDj是氣泡和粒子所受的黏性阻力, 在液體勢(shì)能表達(dá)式中加入氣泡和粒子在平動(dòng)時(shí)受到黏性阻力相關(guān)的項(xiàng); 聲壓P1和黏滯力可以表示為[21]:

式中Pg1=P0+2σl/R10是氣泡內(nèi)部的平衡蒸氣壓;R10是氣泡1 的平衡半徑;cl是液體中的聲速;γ是液體的比熱;σl是表面張力;ηl是液體的黏度;P0是液體中的靜態(tài)壓力;Pac是驅(qū)動(dòng)聲壓;f是聲波頻率;N是循環(huán)數(shù), 其取值決定高斯波的作用時(shí)間;υj表示氣泡(粒子)所在位置與彼此相互作用有關(guān)的液體體元振速, 其表達(dá)式為[15]

需要特別說(shuō)明的是, 對(duì)粒子(j= 2), 在 (23)式中有=0 .

空腔外彈性介質(zhì)層的勢(shì)能Us可表示為[20]

式中μ為固體的剪切模量.

將動(dòng)能表達(dá)式(17)式—(19)式和勢(shì)能表達(dá)式(20)式和(24)式代入拉格朗日方程

可得密閉球形液體腔內(nèi)氣泡和粒子的動(dòng)力學(xué)方程分別為:

式中

動(dòng)力學(xué)方程(27)中出現(xiàn)了腔壁瞬時(shí)半徑Rc對(duì)時(shí)間的高階倒數(shù), 根據(jù)液體體積不變, 結(jié)合(16)式和(27)式, 約去高階小量, 近似有

2.3 球腔內(nèi)氣泡的共振頻率

為分析密閉腔內(nèi)氣泡相對(duì)位置對(duì)氣泡共振頻率的影響, 設(shè)

且有|x1|≤R10,|xc|≤Rc0, 若腔外側(cè)介質(zhì)厚度趨于無(wú)窮, 即Rs?Rc0, 忽略(26)式展開式中的高階小量, 可得到氣泡線性振動(dòng)方程, 即

其中

由此可以看出, 球腔內(nèi)氣泡的共振頻率與腔外介質(zhì)密度、彈性、液體表面張力系數(shù)、環(huán)境壓力及氣泡在球腔內(nèi)的位置等因素有關(guān).

3 數(shù)值分析

基于密閉腔體中液體內(nèi)氣泡和粒子間的相互作用相關(guān)的動(dòng)力學(xué)方程(27)式—(29)式可分析氣泡的徑向振動(dòng)和平動(dòng)以及粒子的平動(dòng), 采用的參數(shù)有:ρl=1000 kg/m3,cl=1500 m/s ,ηl=0.001 Pa·s ,σl=0.072 N/m,ρs=1233 kg/m3,cs=2111 m/s ,γ=1.4 ,μ=6 MPa ,f=0.7 MHz ,P0=101.3 kPa ,Pa=2 atm (1 atm = 101325 Pa),R1=R2=5 μm ,ρp=1.5ρ,Rc0=500 μm,Rs=Rc0+500 μm,d=10 μm,N= 10. 空腔外彈性介質(zhì)層設(shè)為水凝膠[20].

Wang[11]分析了介質(zhì)體積彈性對(duì)彈性介質(zhì)包圍的液體腔中心的氣泡的共振頻率的影響, 發(fā)現(xiàn)隨著液體腔半徑與氣泡半徑比值增大共振頻率逐漸減小并趨近于無(wú)界液體中單氣泡共振頻率. 為更清楚地認(rèn)識(shí)液體腔周圍介質(zhì)對(duì)氣泡共振頻率的影響,基于(37)式我們分析R10= 2 μm 的氣泡的相對(duì)共振頻率ω10/ωm隨半徑比Rc0/R10的變化 (ωm是無(wú)界液體中單氣泡固有頻率), 如圖2 所示, 其中氣泡相對(duì)位置z1/Rc0分別為0, 0.3, 0.5, 對(duì)比發(fā)現(xiàn)相對(duì)位置對(duì)氣泡的共振頻率的影響較小, 幾乎可以忽略不計(jì). 隨著半徑比Rc0/R10的增加, 氣泡的共振頻率有先減小后增大的變化趨勢(shì), 在半徑比7 附近存在極小值. 當(dāng)半徑比較小時(shí)(大約Rc0/R10< 4),氣泡相對(duì)頻率比大于1, 而后相對(duì)頻率比小于1 且隨著液體腔半徑的增大逐漸趨近于1, 與Wang[11]預(yù)測(cè)結(jié)果存在一定的差異, 主要原因在于其未計(jì)入外層介質(zhì)密度的影響. 從密閉腔內(nèi)液體中氣泡的共振頻率表達(dá)式可以看出, 腔外介質(zhì)的彈性的存在可使得共振頻率增加, 然后隨著半徑比的增加, 介質(zhì)彈性對(duì)共振頻率的影響將逐漸減弱, 外層介質(zhì)密度對(duì)共振頻率的影響逐漸占主導(dǎo)地位, 進(jìn)而導(dǎo)致共振頻率比小于1. 隨著腔體半徑的進(jìn)一步增大, 腔壁外側(cè)介質(zhì)特性對(duì)氣泡的影響逐漸減弱, 氣泡的共振頻率趨近于液體中自由泡的共振頻率. 因此, 液體腔外層介質(zhì)對(duì)氣泡的振動(dòng)行為的影響會(huì)隨著腔體尺寸的變化而不同.

圖2 相對(duì)頻率隨空腔半徑與氣泡半徑比值的變化Fig. 2. The relative frequency changes with the ratio of cavity radius to bubble radius.

氣泡的初始半徑不同, 聲場(chǎng)中氣泡的振動(dòng)行為也不同. 為分析氣泡半徑對(duì)粒子和氣泡系統(tǒng)平動(dòng)的影響, 選取初始半徑分別為2, 3, 4, 5 μm 氣泡做數(shù)值分析, 如圖3 所示. 在氣泡和粒子中心間距不變的條件下, 當(dāng)氣泡半徑為2 μm 時(shí), 氣泡和粒子均向腔壁平動(dòng); 隨著氣泡半徑的增大, 氣泡向液體腔中心跳躍式平動(dòng)而粒子向腔壁運(yùn)動(dòng)且其路徑有起伏跳躍特征, 且在氣泡半徑為4 μm 時(shí)相互作用力最強(qiáng). Wu 等[21]發(fā)現(xiàn)無(wú)界液體中氣泡和粒子間相互作用強(qiáng)弱與氣泡的初始半徑相關(guān)并表現(xiàn)為互相吸引. 由于液體腔外層的彈性介質(zhì)會(huì)對(duì)氣泡與粒子間相互作用產(chǎn)生影響, 從而使得液體腔中的氣泡和粒子與無(wú)界液體中的氣泡和粒子的運(yùn)動(dòng)趨勢(shì)不同,并在氣泡半徑不同時(shí)表現(xiàn)出相互吸引或排斥的特征.

圖3 氣泡初始半徑對(duì)氣泡和粒子的平動(dòng)的影響Fig. 3. The effect of initial radius of bubble on the translation of bubble and particle.

液體腔體內(nèi)氣泡和粒子動(dòng)力學(xué)模型可用于分析超聲植物提取、生物組織空化等, 應(yīng)用領(lǐng)域不同, 選用的聲波頻率范圍不同. 為探索頻率的影響, 利用(37)式計(jì)算了5 μm 氣泡的共振頻率約為0.72 MHz并選取頻率分別為0.7, 1.0, 1.5 MHz 的聲波分析液體腔內(nèi)氣泡和粒子的平動(dòng), 如圖4 所示. 當(dāng)頻率為0.7 MHz 接近氣泡的共振頻率時(shí), 氣泡和粒子間相互作用最強(qiáng), 且表現(xiàn)為氣泡向腔中心平動(dòng), 而粒子向腔壁運(yùn)動(dòng). 隨著頻率的增加, 氣泡的徑向振動(dòng)行為相對(duì)逐漸減弱, 二者之間的相互作用也減弱, 在驅(qū)動(dòng)聲波頻率為1.5 MHz 時(shí), 氣泡和粒子間僅有輕微的相互遠(yuǎn)離運(yùn)動(dòng)特征.

圖4 驅(qū)動(dòng)聲波頻率對(duì)氣泡和粒子的平動(dòng)的影響Fig. 4. The influence of driving sound frequency on the translation of bubble and particle.

聲波頻率一定時(shí), 驅(qū)動(dòng)聲壓大小決定氣泡振動(dòng)強(qiáng)弱, 進(jìn)而影響氣泡與粒子的平動(dòng)特征, 圖5 中曲線對(duì)應(yīng)聲壓幅值分別為0.5, 1.0, 1.5, 3.0 atm. 在同樣的聲場(chǎng)條件下, 粒子相較氣泡而言其平動(dòng)更顯著, 主要原因在于氣泡的膨脹推動(dòng)粒子向腔壁運(yùn)動(dòng). 隨著驅(qū)動(dòng)聲壓幅值的增加, 氣泡的徑向振動(dòng)越強(qiáng), 推動(dòng)粒子向腔壁運(yùn)動(dòng)的位移越大. 在聲壓幅值小于1.5 atm 時(shí)氣泡進(jìn)發(fā)生較小的平動(dòng)位移, 而當(dāng)聲壓幅值增大到3.0 atm 時(shí), 氣泡向液體腔中心運(yùn)動(dòng). Vincent 等[13]在實(shí)驗(yàn)中觀察到凝膠中的液體腔內(nèi)氣泡由于外層彈性介質(zhì)的影響將向中心位置運(yùn)動(dòng), 而在聲波作用下, 氣泡和粒子間的相互作用在一定程度上能夠平衡外層彈性介質(zhì)的影響, 但當(dāng)驅(qū)動(dòng)聲壓較高時(shí), 二者相互作用增強(qiáng), 將表現(xiàn)為更為復(fù)雜的相互吸引和排斥的特征.

液體腔的大小對(duì)密閉腔內(nèi)氣泡的運(yùn)動(dòng)狀態(tài)有著顯著影響, 為更進(jìn)一步分析影響特征, 改變內(nèi)腔半徑分別為100, 300, 500 μm 分析粒子和氣泡的平動(dòng), 如圖6 所示. 在不同尺度的液體腔內(nèi)的粒子和氣泡在聲波作用下的運(yùn)動(dòng)行為表現(xiàn)為: 粒子朝腔壁運(yùn)動(dòng)且隨著內(nèi)腔半徑的增大粒子的運(yùn)動(dòng)位移變大; 對(duì)于氣泡而言, 在內(nèi)腔半徑較小時(shí), 外層彈性介質(zhì)對(duì)氣泡的振動(dòng)行為的影響更為顯著, 因此, 在一定程度影響可調(diào)節(jié)粒子和氣泡相互吸引或排斥的狀態(tài), 表現(xiàn)為腔內(nèi)徑100 μm 時(shí)氣泡向著腔壁遠(yuǎn)離粒子運(yùn)動(dòng); 隨著內(nèi)腔半徑的增大, 粒子向腔壁運(yùn)動(dòng)的位移變大, 聲波作用下氣泡則朝向球腔中心有小的位移變化.

圖5 驅(qū)動(dòng)聲壓幅值對(duì)氣泡和粒子的平動(dòng)的影響Fig. 5. The influence of sound pressure amplitude on the translation of bubble and particle.

圖6 腔體大小對(duì)氣泡和粒子平動(dòng)的影響Fig. 6. The influence of cavity size on the translation of bubble and particle.

液體腔內(nèi)由于粒子的存在影響氣泡的徑向振動(dòng)和平動(dòng). 為更進(jìn)一步分析粒子密度對(duì)密閉腔內(nèi)氣泡的運(yùn)動(dòng)狀態(tài)的影響, 比較粒子密度ρp分別為0.5ρ,ρ和1.5ρ時(shí)粒子和氣泡的平動(dòng)狀態(tài)變化, 如圖7 所示. 隨著粒子密度的增大, 粒子向腔壁運(yùn)動(dòng)的位移逐漸減小,原因在于重粒子平動(dòng)狀態(tài)發(fā)生變化需要更大的推力作用.對(duì)氣泡而言,粒子密度越大氣泡向著球腔中心運(yùn)動(dòng)位移也越小,因此氣泡和粒子間的相互作用同時(shí)影響二者的平動(dòng)行為.

圖7 粒子密度對(duì)氣泡和粒子平動(dòng)的影響Fig.7.The influence of particle density on the translation of bubble and particle.

4 結(jié) 論

本文構(gòu)建了彈性介質(zhì)包裹的液體球腔內(nèi)粒子和氣泡間的相互作用動(dòng)力學(xué)模型,基于拉格朗日方程得到了表征氣泡徑向振動(dòng)和平動(dòng)、粒子平動(dòng)動(dòng)力學(xué)方程.基于此方程,可分析聲場(chǎng)參數(shù)、球腔約束條件(如腔體大小)、氣泡和粒子本身特征等因素對(duì)系統(tǒng)運(yùn)動(dòng)的影響.球形液體腔外層介質(zhì)的彈性、密度、腔體結(jié)構(gòu)等均可影響氣泡的共振頻率,由于介質(zhì)彈性的影響,當(dāng)氣泡半徑與球腔半徑比較接近時(shí),氣泡共振頻率大于無(wú)界液體中的共振頻率;隨著球腔半徑的增加,氣泡共振頻率逐漸減小到低于無(wú)界液體中的共振頻率的情形.由于氣泡的徑向振動(dòng)形成次級(jí)聲輻射引發(fā)粒子和氣泡間的相互作用,進(jìn)而影響二者在液體腔內(nèi)的平動(dòng).氣泡的徑向振動(dòng)可推動(dòng)粒子向著腔壁運(yùn)動(dòng),且氣泡的振動(dòng)越強(qiáng),粒子在聲波持續(xù)時(shí)間內(nèi)的平動(dòng)位移越大,同樣,氣泡的平動(dòng)也會(huì)受到二者相互作用的調(diào)制,在不同的相互作用情形下存在向球腔中心運(yùn)動(dòng)和向腔壁運(yùn)動(dòng)兩種情形.

氣泡在球狀液體腔內(nèi)的運(yùn)動(dòng)可改變液體內(nèi)的速度場(chǎng)分布,同時(shí)影響周圍的物質(zhì)運(yùn)動(dòng)狀態(tài).氣泡和粒子間的相互作用在一定程度上可調(diào)控氣泡的運(yùn)動(dòng)行為,進(jìn)而增強(qiáng)液體腔內(nèi)外物質(zhì)交換,可為解釋超聲植物有效成分提取及生物組織內(nèi)的空化行為提供理論基礎(chǔ).

附 錄

為了得到空腔中液體的動(dòng)能表達(dá)式,需要利用前文求得的液體運(yùn)動(dòng)的解進(jìn)行推導(dǎo).腔內(nèi)液體動(dòng)能可表示為

式中,ρl為與液體的密度,Vl為液體所占體積.根據(jù)變換關(guān)系,(A 1)式可進(jìn)一步表示為

式中,Sc,S1,S2分別表示空腔、氣泡和粒子的表面;nc,n1,n2分別表示垂直于空腔、氣泡和粒子表面并指向液體的單位向量.結(jié)合(2)式—(15)式后液體速度勢(shì)(1)式代入(A2)式,得

將(11)式、(12)式代入(A 3)式,液體動(dòng)能的具體表達(dá)式為