齊紹富 蔡飛燕 田振 黃先玉5) 周娟 王金萍 李文成 鄭海榮 鄧科

1) (吉首大學(xué)信息科學(xué)與工程學(xué)院,吉首 416000)

2) (吉首大學(xué)物理與機(jī)電工程學(xué)院,吉首 416000)

3) (中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院,勞特伯生物醫(yī)學(xué)成像研究中心,深圳 518055)

4) (湖南醫(yī)藥學(xué)院醫(yī)學(xué)院,懷化 418000)

5) (綏寧縣第一中學(xué),邵陽(yáng) 422699)

聲操控微粒技術(shù)可以非接觸無(wú)損傷地控制聲場(chǎng)中的物體運(yùn)動(dòng),其在精密制造、材料工程、體外診斷等領(lǐng)域具有廣闊的應(yīng)用前景.傳統(tǒng)聲操控微粒技術(shù)一般采用自由聲場(chǎng),如利用單個(gè)換能器或陣列換能器產(chǎn)生的聚焦聲場(chǎng)、行波場(chǎng)或駐波場(chǎng)等.然而,一般單個(gè)換能器產(chǎn)生的聲場(chǎng)僅能操控單個(gè)微粒;而陣列換能器的驅(qū)動(dòng)系統(tǒng)復(fù)雜,導(dǎo)致操控器件成本高昂且難以微型化;因此,亟需研究新的聲場(chǎng)形態(tài)實(shí)現(xiàn)多樣性微粒操控.本工作中,采用單個(gè)換能器產(chǎn)生的平面波激發(fā)一維聲柵的共振聲場(chǎng),實(shí)驗(yàn)實(shí)現(xiàn)了大規(guī)模泡沫微球的周期排列操控.其操控機(jī)制是由于聲柵狹縫中法布里-珀羅諧振聲場(chǎng)與聲柵表面周期衍射場(chǎng)共振耦合,在聲柵表面形成周期分布的局域梯度聲場(chǎng),導(dǎo)致微粒在平行于聲柵表面受到聲捕獲力,在垂直于聲柵表面受到指向表面的聲吸引力,實(shí)現(xiàn)了微粒周期排列在聲柵表面上.該工作為利用超聲在空氣中大規(guī)模排列微粒提供了理論基礎(chǔ)和技術(shù)支持.

1 引言

聲操控微粒是利用聲輻射力捕獲、操縱和移動(dòng)微小顆粒的一項(xiàng)技術(shù).其中聲輻射力是由于處于聲場(chǎng)中的微粒對(duì)聲波的散射、吸收等導(dǎo)致聲場(chǎng)攜帶的動(dòng)量在聲場(chǎng)與微粒之間交換產(chǎn)生[1,2].聲操控微粒技術(shù)具有非接觸、無(wú)損傷、可穿透非透明介質(zhì)等優(yōu)點(diǎn),在精密制造、材料工程、體外診斷等領(lǐng)域具有廣闊應(yīng)用前景[3?8].

傳統(tǒng)聲操控微粒技術(shù)一般是基于單個(gè)換能器或者換能器陣列產(chǎn)生的自由聲場(chǎng)進(jìn)行微粒操控.其中單個(gè)換能器與全反射面形成的駐波可以實(shí)現(xiàn)微粒的停駐,但難以個(gè)性化操控微粒;面陣換能器可以實(shí)現(xiàn)微粒的任意運(yùn)動(dòng),但其驅(qū)動(dòng)電路系統(tǒng)復(fù)雜,導(dǎo)致操控器件成本高昂且難以微型化.近期,聲人工結(jié)構(gòu)可以構(gòu)建特定的調(diào)制聲場(chǎng)實(shí)現(xiàn)微粒的個(gè)性化操控,引起了人們的極大興趣,如: Li 等[9,10]利用聲子晶體板的周期結(jié)構(gòu)共振激勵(lì)板子零階非對(duì)稱(chēng)Lamb 波模式,獲得周期梯度局域強(qiáng)場(chǎng),實(shí)現(xiàn)水中大規(guī)模微粒的捕獲、篩選和移動(dòng);Melde 等[11]設(shè)計(jì)全息板結(jié)構(gòu)實(shí)現(xiàn)水表面輕軟微粒的任意形態(tài)排列;Memoli 等[12]利用超材料陣列調(diào)控相位空間分布獲得瓶形聲場(chǎng)分布并實(shí)現(xiàn)空氣中微粒的捕獲.由于調(diào)制聲場(chǎng)可以通過(guò)單陣元換能器與人工結(jié)構(gòu)結(jié)合產(chǎn)生,該裝置簡(jiǎn)單易集成,展現(xiàn)了優(yōu)越的實(shí)用性.因此,需要研究更多可用于微粒操控的人工聲場(chǎng),實(shí)現(xiàn)微粒的個(gè)性化操控.

在本工作中,基于前期理論研究微粒在一維聲柵表面的受力特征[13],構(gòu)建了空氣中的一維聲柵系統(tǒng),并實(shí)驗(yàn)觀(guān)察到泡沫微球在聲柵表面周期排列.首先探討了聲柵表面的共振聲場(chǎng)特征;接著,研究了微粒在共振聲場(chǎng)中受到的聲輻射力特征,發(fā)現(xiàn)微粒在平行于聲柵表面受到捕獲力,在垂直于聲柵表面受到指向表面的吸引力;進(jìn)一步,制備聲柵樣品并搭建實(shí)驗(yàn)平臺(tái),實(shí)現(xiàn)了泡沫微球在聲柵表面的周期排列操控.該工作為空氣中大規(guī)模并行操控微粒提供了實(shí)驗(yàn)基礎(chǔ).

2 理論模擬

2.1 聲柵共振聲場(chǎng)特征分析

我們研究的聲柵是由鋼板刻蝕周期分布的狹縫構(gòu)成,該聲柵結(jié)構(gòu)放置在空氣中.為了便于觀(guān)察,選取了聲柵三個(gè)周期結(jié)構(gòu)作為計(jì)算單元,如圖1(a)所示,其中鋼板的厚度t=6.25 mm,狹縫寬度c=1 mm,狹縫周期為d=8.34 mm .被操控的泡沫球微粒放置在聲柵表面,其半徑約為r=0.5 mm .材料的聲學(xué)參數(shù) [14]如表1 所列.

表1 材料聲學(xué)參數(shù)Table 1.Material acoustic parameters.

采用 COMSOL Multiphysics?軟件 [15]模擬聲柵的聲場(chǎng)特征.計(jì)算模型采用壓力聲學(xué)和固體力學(xué)模塊,研究類(lèi)型采用頻域模塊;左右邊界采用周期性邊界條件,用于模擬模型周期結(jié)構(gòu);上下邊界采用平面波輻射邊界條件,用于模擬開(kāi)放空間;下邊界設(shè)置為平面波入射聲場(chǎng),入射方向是從底面向上激發(fā);泡沫球位于聲柵上表面,其位置坐標(biāo)如圖1(a)所示.圖1(b)展示了聲柵系統(tǒng)的聲透射譜,可以看到在38—41 kHz 之間,有一個(gè)透射增強(qiáng)峰,該頻率約為39.6 kHz;圖1(c)展示了系統(tǒng)在該共振頻率點(diǎn)的聲場(chǎng)分布,為了對(duì)比聲場(chǎng)特征;圖1(d)展示了系統(tǒng)在非共振頻率38 kHz 時(shí)的聲場(chǎng)分布,其中紅色和藍(lán)色分別代表聲壓的極大值和極小值,聲場(chǎng)強(qiáng)度用最大值歸一化.對(duì)比共振頻率和非共振頻率的聲場(chǎng)特征,可以看到在共振頻率處,聲柵表面聲場(chǎng)得到了極大增強(qiáng);而且在共振頻率處,聲場(chǎng)除了局域在狹縫里,在柵格表面也存在有梯度分布;在一個(gè)周期柵格表面有兩個(gè)零值聲壓[16].該共振聲場(chǎng)是由于聲柵表面激發(fā)的衍射波與狹縫內(nèi)部的法布里-珀羅諧振腔共振后兩者發(fā)生耦合產(chǎn)生[13,17].

圖1 聲柵樣品示意圖與歸一化透射譜,以及一維聲柵歸一化聲壓分布圖 (a) 聲柵樣品示意圖;(b) 平面波透過(guò)聲柵樣品的透射譜;(c) 共振頻率為39.6 kHz 時(shí)聲場(chǎng)分布圖;(d) 非共振頻率為38 kHz 時(shí)聲場(chǎng)分布圖Fig.1.Schematic diagram of the grating sample with normalized transmission spectrum,and normalized sound pressure distribution of the 1D grating: (a) Schematic diagram of the grating sample;(b) transmission spectrum of the plane wave through the grating sample;(c) the pressure field at resonant frequency of 39.6 kHz;(d) the pressure field at off resonant frequency of 38 kHz.

2.2 微球在聲柵表面所受聲輻射力特征分析

聲柵結(jié)構(gòu)的共振耦合導(dǎo)致其表面產(chǎn)生了周期分布局域聲場(chǎng),該局域聲場(chǎng)可以用于操控微粒的運(yùn)動(dòng).在本實(shí)驗(yàn)中,擬操控泡沫微球半徑約r=0.5mm,工作頻率位于約39.6 kHz,可以計(jì)算得出kr≈0.37.雖然kr ≈0.37接近引起Mie 散射的頻率范圍,但由于微球的聲阻抗遠(yuǎn)大于空氣的聲阻抗,在該頻率處沒(méi)有引起微球共振導(dǎo)致的反常聲輻射力產(chǎn)生.因此,可以用Gor’kov 理論近似模擬顆粒在聲場(chǎng)中受到的聲輻射力[18,19].該聲輻射力的表達(dá)式如下:

其中U為空間聲場(chǎng)能量勢(shì),可以表示為

為了進(jìn)一步了解泡沫球在共振頻率下聲柵表面所受的聲輻射力特征,展示了計(jì)算獲得的聲柵表面泡沫球在共振頻率39.6 kHz 時(shí)所受歸一化聲輻射力空間分布,如圖2 所示;其中圖2(a)顏色表示聲輻射力大小,箭頭表示聲輻射力方向;圖2(b)是泡沫球離聲柵Z方向0.24d處,X方向和Z方向的聲輻射力Fx,Fz與 ?X的關(guān)系,其中 ?X的變化范圍是–0.5d到0.5d.從圖2(a)可以看出,泡沫球在聲柵表面受到的聲輻射力呈周期變化,當(dāng)泡沫球逐漸靠近聲柵表面,所受聲輻射力逐漸增大.從圖2(b)可以看出,在聲柵表面處,微粒在水平X方向一個(gè)周期內(nèi)有三個(gè)零值點(diǎn),坐標(biāo)分別為–0.25d,0,0.25d,在三個(gè)零值點(diǎn)中,–0.25d和0.25d兩點(diǎn)微粒受到的是捕獲力,微?？梢苑€(wěn)定停駐,0 點(diǎn)是非穩(wěn)態(tài)停住點(diǎn)[20].從圖2 可以進(jìn)一步得出,在–0.25d和0.25d兩點(diǎn)處,微粒受到垂直Z方向的力是負(fù)的[13,21];即與波傳播方向相反,指向板面方向.因此,從理論分析結(jié)果可以看出,微粒能夠在一個(gè)柵格周期內(nèi)穩(wěn)定吸附在–0.25d和0.25d兩個(gè)位置處.

圖2 泡沫球位置變化在X方向上從–1.5d到1.5d,在Z方向上從0.24d 到1.2d,半徑約為 r=0.5 mm 的泡沫球在一維聲柵上所受歸一化聲輻射力空間分布圖 (a) 黑色三角形箭頭表示聲輻射力所指方向,顏色深淺表示聲輻射力的大小,灰色矩形表示一維聲柵;(b) X 方向上泡沫球位置變化從–0.5d 到0.5d 的歸一化聲輻射力Fx (藍(lán)色實(shí)線(xiàn))與Fz (紅色實(shí)線(xiàn))曲線(xiàn)圖Fig.2.Spatial distribution of the normalized acoustic radiation force on a one-dimensional sound grid for a foam sphere with a radius of about r=0.5 mm varying from–1.5d to 1.5d in the X-direction and from 0.24d to 1.2d in the Z-direction: (a) The black triangle arrow indicates the direction of the acoustic radiation force,the color shade indicates the size of the acoustic radiation force,and the gray rectangle indicates the one-dimensional acoustic grid;(b) plot of normalized acoustic radiation force Fx (blue solid line)versus Fz (red solid line) for the change in position of the foam sphere from–0.5d to 0.5d in the X-direction.

3 實(shí)驗(yàn)觀(guān)察

我們搭建的聲操控微粒實(shí)驗(yàn)平臺(tái)如圖3 所示.中心頻率為40 kHz 的空氣聲換能器連接功率放大器和信號(hào)發(fā)生器;不銹鋼聲柵結(jié)構(gòu)通過(guò)模具加工制備,放置在換能器上方;泡沫微球任意排布在聲柵表面;攝像機(jī)放置在聲柵頂面用于觀(guān)察實(shí)驗(yàn)現(xiàn)象.

圖3 聲輻射力操控實(shí)驗(yàn)平臺(tái)示意圖Fig.3.Schematic diagram of acoustic radiation force manipulation experimental platform.

圖4 展示了超聲換能器關(guān)閉和開(kāi)啟時(shí)的排列實(shí)驗(yàn)效果圖.實(shí)驗(yàn)錄制的視頻在附件中,其中截取的瞬時(shí)圖片如圖4 所示.從圖4(a)可以看到,未開(kāi)啟超聲時(shí),聲柵結(jié)構(gòu)表面泡沫球微粒都是隨機(jī)放置;圖4(b)可以看到開(kāi)啟超聲后,且聲場(chǎng)頻率為39.6 kHz時(shí),泡沫球在每個(gè)聲柵表面整齊地排成兩列,其排列的位置是位于離狹縫–0.25d和0.25d處,獲得了與理論分析一致的結(jié)果.

圖4 泡沫球操控實(shí)驗(yàn)效果圖 (a) 超聲開(kāi)啟前泡沫球分布圖;(b) 超聲開(kāi)啟后泡沫球排列圖Fig.4.Experimental effect of foam ball manipulation:(a) Distribution of foam balls before ultrasonic opening;(b) arrangement of foam balls after ultrasonic opening.

4 結(jié)論

本文在實(shí)驗(yàn)上實(shí)現(xiàn)了基于一維聲柵共振場(chǎng)的大規(guī)模微粒周期排列操控.其主要的物理機(jī)制是由于平面聲場(chǎng)激勵(lì)聲柵產(chǎn)生共振,導(dǎo)致聲柵狹縫內(nèi)部的法布里-珀羅諧振聲場(chǎng)與表面的衍射聲場(chǎng)共振耦合,在聲柵表面產(chǎn)生周期分布的局域梯度聲場(chǎng),使得聲柵表面微粒受到平行于聲柵表面的捕獲力和垂直于聲柵表面的吸引力,從而獲得周期排列操控.我們的實(shí)驗(yàn)結(jié)果與有限元理論模擬結(jié)果一致.由于聲柵結(jié)構(gòu)簡(jiǎn)單易加工,無(wú)需復(fù)雜大規(guī)模電路裝置就可以實(shí)現(xiàn)周期分布的局域梯度聲場(chǎng),該工作為大規(guī)模操控微粒提供了一種簡(jiǎn)單易行的方法.