基于一維聲柵共振場的大規(guī)模微粒并行排列的實(shí)驗(yàn)研究*

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

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

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

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

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

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

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

1 引言

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

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

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

2 理論模擬

2.1 聲柵共振聲場特征分析

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

圖1 聲柵樣品示意圖與歸一化透射譜,以及一維聲柵歸一化聲壓分布圖 (a) 聲柵樣品示意圖;(b) 平面波透過聲柵樣品的透射譜;(c) 共振頻率為39.6 kHz 時聲場分布圖;(d) 非共振頻率為38 kHz 時聲場分布圖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)生了周期分布局域聲場,該局域聲場可以用于操控微粒的運(yùn)動.在本實(shí)驗(yàn)中,擬操控泡沫微球半徑約r=0.5mm,工作頻率位于約39.6 kHz,可以計(jì)算得出kr≈0.37.雖然kr ≈0.37接近引起Mie 散射的頻率范圍,但由于微球的聲阻抗遠(yuǎn)大于空氣的聲阻抗,在該頻率處沒有引起微球共振導(dǎo)致的反常聲輻射力產(chǎn)生.因此,可以用Gor’kov 理論近似模擬顆粒在聲場中受到的聲輻射力[18,19].該聲輻射力的表達(dá)式如下:

其中U為空間聲場能量勢,可以表示為

為了進(jìn)一步了解泡沫球在共振頻率下聲柵表面所受的聲輻射力特征,展示了計(jì)算獲得的聲柵表面泡沫球在共振頻率39.6 kHz 時所受歸一化聲輻射力空間分布,如圖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方向一個周期內(nèi)有三個零值點(diǎn),坐標(biāo)分別為–0.25d,0,0.25d,在三個零值點(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é)果可以看出,微粒能夠在一個柵格周期內(nèi)穩(wěn)定吸附在–0.25d和0.25d兩個位置處.

圖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í)線)與Fz (紅色實(shí)線)曲線圖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)觀察

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

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

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

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