王浩然 藍(lán)君2) 陳佳惠 李義豐2)?
1) (南京工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院, 南京 211800)
2) (近代聲學(xué)教育部重點(diǎn)實(shí)驗(yàn)室, 南京大學(xué)聲學(xué)研究所, 南京 210093)
近年來, 在空氣聲中實(shí)現(xiàn)聲波分量的局域增強(qiáng)這項(xiàng)研究工作逐漸受到關(guān)注.若聲場(chǎng)能夠被大幅度增強(qiáng), 這將會(huì)發(fā)掘許多有價(jià)值的潛在應(yīng)用, 例如聲傳感器、聲通信等.在自然界中, 由于聲能具有很低的能量密度, 空氣中聲場(chǎng)局域增強(qiáng)效應(yīng)的實(shí)現(xiàn)通常比較困難.因此, 探索新的物理機(jī)制和方法去實(shí)現(xiàn)聲場(chǎng)增強(qiáng)具有重要的意義.
聲學(xué)超構(gòu)材料作為一種人工專門設(shè)計(jì)的復(fù)合結(jié)構(gòu), 可通過局域共振機(jī)理來調(diào)控彈性波的能帶結(jié)構(gòu), 從而獲得自然材料所不具備的物理特性, 如負(fù)等效體積模量[1,2]、負(fù)等效質(zhì)量密度[3-5]、零折射率[6]等.超構(gòu)材料在對(duì)聲波傳播方式的調(diào)控方面起重要作用, 目前, 基于亥姆霍茲共振器[7]和腔諧振器[8,9]的局域共振機(jī)理設(shè)計(jì)的超構(gòu)材料已經(jīng)實(shí)現(xiàn)了聲場(chǎng)增強(qiáng)效應(yīng).Yuan等[10,11]分別通過優(yōu)化亥姆霍茲共振器底部和頸部結(jié)構(gòu)來降低聲學(xué)粘性損耗, 實(shí)現(xiàn)了低頻聲波的局域增強(qiáng)效應(yīng).Song[12]等設(shè)計(jì)了一種新型聲腔型雙壁超構(gòu)材料, 其能夠在亞波長(zhǎng)的空氣隙內(nèi)發(fā)生強(qiáng)烈的局域聲場(chǎng)增強(qiáng)效應(yīng).此外, 近年來, 超構(gòu)材料的聲學(xué)Mie共振得到越來越多的關(guān)注, 相比局域共振, Mie共振具有更為豐富的共振模式, 如單極子、偶極子、四極子等更高的多極子共振模式.利用聲學(xué)Mie共振的不同共振模式可分別實(shí)現(xiàn)不同效果的局域聲場(chǎng)增強(qiáng).Zhu等[13]提出一種偶極子共振的亞波長(zhǎng)尺度下的指向性傳感結(jié)構(gòu), 結(jié)構(gòu)內(nèi)聲壓增幅與卷曲空間外殼的折疊程度相關(guān).Zhao等[14]設(shè)計(jì)了一種簡(jiǎn)并Mie共振的亞波長(zhǎng)殼結(jié)構(gòu), 在不改變聲源方向性的前提下, 可以實(shí)現(xiàn)單級(jí)源和多級(jí)源的聲波增強(qiáng)發(fā)射.上述方法都通過不同共振形式和結(jié)構(gòu)實(shí)現(xiàn)了聲場(chǎng)增強(qiáng)效應(yīng).
另外, 近年來利用聲學(xué)超構(gòu)材料設(shè)計(jì)實(shí)現(xiàn)的非對(duì)稱聲傳輸開始受到研究者的高度重視[15-17].研究人員已經(jīng)通過非線性效應(yīng)[18,19]、近零折射率材料特性[20,21]、超構(gòu)表面的表面消逝波等原理實(shí)現(xiàn)了非對(duì)稱聲傳輸.Jin等[18]和Liang等[19]利用非線性介質(zhì)和一維聲子晶體相結(jié)合的方法提出了一種聲整流模型.Li等[22]提出了一種由近零折射率超構(gòu)材料設(shè)計(jì)而成的聲學(xué)棱鏡模型, 其利用近零折射率材料對(duì)入射聲波角度的高度依賴性, 實(shí)現(xiàn)了聲波單向傳輸.Shen等[23]利用該原理, 通過采用梯度折射率超構(gòu)表面和近零折射率超構(gòu)表面組合而成的雙層復(fù)合結(jié)構(gòu), 產(chǎn)生了強(qiáng)非對(duì)稱聲傳輸效應(yīng).Song等[24]基于超構(gòu)表面的表面消逝波實(shí)現(xiàn)了非對(duì)稱聲傳輸效果.若能將聲場(chǎng)增強(qiáng)和非對(duì)稱聲傳輸這兩種功能在一個(gè)系統(tǒng)中同時(shí)實(shí)現(xiàn), 則有望實(shí)現(xiàn)具有非對(duì)稱性質(zhì)的聲場(chǎng)增強(qiáng)效應(yīng), 其將在隔聲、非對(duì)稱聲學(xué)器件等方面具有重要的應(yīng)用價(jià)值.
基于此, 本文構(gòu)建了一種多腔型基本單元結(jié)構(gòu), 利用其組成不同形狀的超構(gòu)材料, 并使每種超構(gòu)材料都含有豐富的Mie共振特性.聲場(chǎng)增強(qiáng)效應(yīng)是基于基本單元的聲腔和系統(tǒng)結(jié)構(gòu)之間相互耦合產(chǎn)生的單極子Mie共振來實(shí)現(xiàn).由于設(shè)計(jì)的超構(gòu)材料對(duì)聲波入射方向具有依賴性, 利用這一特性可以實(shí)現(xiàn)聲場(chǎng)的非對(duì)稱增強(qiáng).研究表明, 對(duì)稱型超構(gòu)材料的聲場(chǎng)增強(qiáng)效應(yīng)不受聲波入射方向的影響,而非對(duì)稱型超構(gòu)材料的聲場(chǎng)增強(qiáng)效應(yīng)具有較強(qiáng)的方向依賴性.
本文設(shè)計(jì)的多腔型基本單元的結(jié)構(gòu)如圖1(a)所示, 該基本單元的外部形狀是1個(gè)圓形, 其由1個(gè)圓形中心腔和12個(gè)諧振腔共同組成, 可以分為4個(gè)相同的部分, 其中12個(gè)諧振腔以3個(gè)為一組圍繞中心腔的四周均勻分布.中心腔和諧振腔內(nèi)填充的媒質(zhì)均為空氣, 腔壁選取的材料為環(huán)氧樹脂.腔壁的厚度為t, 中心腔的開口寬度均為a, 諧振腔的開口寬度均為w, 圓形中心腔的半徑為r1,基本單元的外部圓形的半徑為R1.所有的結(jié)構(gòu)參數(shù)分別為t = 3 mm, a = 5 mm, w = 4 mm, r1=2.7 cm和R1= 5 cm.環(huán)氧樹脂的材料參數(shù)分別為密度 ρ1=1050kg/m3、楊氏模量E = 5.08 GPa和泊松比v = 0.35; 背景媒質(zhì)空氣的密度和聲速分別為ρ0= 1.21 kg/m3和c0= 343 m/s.本文在r1和R1之間引入兩個(gè)具有相同場(chǎng)的虛擬層, 基于聲散射理論可以計(jì)算出多腔型基本單元的各階Mie共振所對(duì)應(yīng)的共振頻率.將入射聲源視為平面波pin=p0ei(krcosφ-ωt), 其中 p0為聲壓振幅, i表示虛部, k是空氣中的波數(shù), ω 是角頻率, φ 是相位, r是半徑.聲壓分布可以用Bessel (Jm)函數(shù)和Hankel(Hm)函數(shù)表示:
圖1 (a)多腔型基本單元的結(jié)構(gòu)示意圖; (b)特征頻率為672 Hz時(shí), 基本單元在單極子Mie共振模式下的聲強(qiáng)分布圖; (c)當(dāng)頻率為672 Hz的聲波從結(jié)構(gòu)左側(cè)垂直入射時(shí), 多腔型基本單元的聲強(qiáng)分布圖Fig.1.(a) Structural diagram of multiple-cavity unit; (b) sound intensity distribution of the unit in monopole Mie resonance mode at the characteristic frequency of 672 Hz; (c) sound intensity distribution of the unit when the sound wave with frequency of 672 Hz is normally incident from the left side.
其中 δ =4w/2πR1是四通道圓的占空比,k0=ω/c0是空氣中的波數(shù), k′是基本單元的等效波數(shù),Cm—Gm是各層m階的展開系數(shù).當(dāng) | Dm/Cm| 為1時(shí), 代表散射聲場(chǎng)達(dá)到最大, 此時(shí)可以求得第m階Mie共振模式的共振頻率.當(dāng)m = 0時(shí), 單極子Mie共振模式所對(duì)應(yīng)的共振頻率可以近似地表示為[25], 其中 ηr為基本單元的折射率, 根據(jù)反演法[26]求得基本單元的等效折射率 ηr約等于2.1, 則理論計(jì)算得到該基本單元的單極子Mie共振模式所對(duì)應(yīng)的共振頻率為666 Hz.當(dāng)m為其他正整數(shù)時(shí), 可以得到其他共振模式所對(duì)應(yīng)的共振頻率.
采用有限元仿真軟件COMSOL的聲固耦合模塊對(duì)多腔型基本單元的聲學(xué)特性進(jìn)行仿真驗(yàn)證.在Mie共振單元結(jié)構(gòu)中通過調(diào)整空氣通道內(nèi)聲波的相位可實(shí)現(xiàn)單極子共振、偶極子共振以及其他多重Mie共振.其中在單極子共振模式下, 聲能能夠聚集在結(jié)構(gòu)內(nèi)部, 且內(nèi)部中心處的聲強(qiáng)大小是入射聲波的數(shù)倍, 因此本文將采用多腔型基本單元的單極子Mie共振在結(jié)構(gòu)內(nèi)部實(shí)現(xiàn)聲場(chǎng)增強(qiáng)效應(yīng)[27].圖1(b)顯示了特征頻率為672 Hz時(shí)所對(duì)應(yīng)的單極子Mie共振模式下的聲強(qiáng)分布, 與上述利用Mie散射理論所得的單極子共振頻率(666 Hz)基本一致, 可以明顯觀察到多腔型基本單元結(jié)構(gòu)內(nèi)部出現(xiàn)聲場(chǎng)增強(qiáng)的現(xiàn)象, 且聲壓大小遠(yuǎn)比結(jié)構(gòu)外部的大,這也證實(shí)了該多腔型基本單元可用于實(shí)現(xiàn)聲場(chǎng)增強(qiáng)效應(yīng).圖1(c)是頻率為672 Hz的平面波從左側(cè)垂直入射到基本單元后的聲強(qiáng)分布圖, 其中入射平面波的聲壓振幅為1 Pa.觀察此圖可以發(fā)現(xiàn), 聲能主要集中在該基本單元的內(nèi)部, 且其內(nèi)部中心點(diǎn)處的聲場(chǎng)強(qiáng)度可高達(dá)30 Pa2, 但這個(gè)增強(qiáng)的聲場(chǎng)沿著圓形中心腔的徑向衰減, 單元外部的強(qiáng)度均接近于0.因此, 設(shè)計(jì)多腔型基本單元具有良好的聲場(chǎng)增強(qiáng)效應(yīng).
基于該多腔型基本單元, 設(shè)計(jì)了一系列由多個(gè)基本單元組成的復(fù)合結(jié)構(gòu), 用以實(shí)現(xiàn)復(fù)合結(jié)構(gòu)的單極子共振.如圖2(a)所示, 將8個(gè)相同的多腔型基本單元組成了一個(gè)外部形狀為正方形的對(duì)稱型超構(gòu)材料, 基本單元的結(jié)構(gòu)參數(shù)與圖1(a)的相同,相鄰兩個(gè)基本單元中心點(diǎn)之間的距離d為15 cm.圖2(b)給出了該對(duì)稱型超構(gòu)材料在特征頻率為669 Hz時(shí)所對(duì)應(yīng)的復(fù)合單極子Mie共振模式下的聲強(qiáng)分布.在該頻率點(diǎn)處, 設(shè)計(jì)的正方形對(duì)稱型超構(gòu)材料處于單極子共振模式, 可以將聲能集中在系統(tǒng)的中心區(qū)域, 實(shí)現(xiàn)聲場(chǎng)增強(qiáng)效應(yīng).
圖2 (a)由8個(gè)多腔型基本單元組成的正方形系統(tǒng)示意圖; (b)特征頻率為669 Hz時(shí), 正方形系統(tǒng)在復(fù)合單極子Mie共振模式下的聲強(qiáng)分布圖Fig.2.(a) Schematic diagram of a square system composed of eight multiple-cavity units; (b) sound intensity distribution of the square system in the compound monopole Mie resonance mode at the characteristic frequency of 669 Hz.
上述的研究通過使用特征頻率分析求得該正方形系統(tǒng)的單極子Mie共振模式所對(duì)應(yīng)的共振頻率為669 Hz.圖3(a)和圖3(b)分別是頻率為669 Hz的平面波從該正方形系統(tǒng)的左側(cè)和右側(cè)垂直入射時(shí)系統(tǒng)的聲強(qiáng)分布圖, 可以發(fā)現(xiàn), 該正方形系統(tǒng)的復(fù)合單極子共振模式在669 Hz頻率處被激發(fā), 聲場(chǎng)增強(qiáng)效應(yīng)在此頻率處實(shí)現(xiàn).同時(shí), 聲波從左側(cè)和右側(cè)垂直入射時(shí)聲強(qiáng)的分布均相同, 說明該系統(tǒng)具有良好的對(duì)稱性.為了與多腔型基本單元作對(duì)比,與其樣品同尺寸同分布的8個(gè)硬邊界圓形結(jié)構(gòu)被放置在相同位置, 相應(yīng)的聲強(qiáng)分布如圖3(c)所示.通過對(duì)比圖3(a)、圖3(b)與圖3(c), 可以發(fā)現(xiàn)在相同入射聲波的情況下, 正方形對(duì)稱型超構(gòu)材料樣品內(nèi)部中心點(diǎn)處的聲強(qiáng)明顯比硬邊界圓形結(jié)構(gòu)的要大, 這意味著超構(gòu)材料樣品具有更好的聲能聚焦效果, 有效地體現(xiàn)了該對(duì)稱型超構(gòu)材料在復(fù)合單極子共振模式下實(shí)現(xiàn)的聲場(chǎng)增強(qiáng)效應(yīng).
圖3 聲波從(a)左側(cè)和(b)右側(cè)垂直入射到8個(gè)多腔型基本單元組成的正方形系統(tǒng)時(shí)的聲強(qiáng)分布圖; (c)聲波從左側(cè)垂直入射到由8個(gè)硬邊界圓形結(jié)構(gòu)組成的系統(tǒng)時(shí)的聲強(qiáng)分布Fig.3.Sound intensity distributions of the square system composed of eight multiple-cavity units when the sound wave is normally incident from (a) the left side and (b) the right side, respectively; (c) sound intensity distribution of the system composed of eight hard boundary circular structures when the sound wave is normally incident from the left side.
除此之外, 還可以利用該多腔型基本單元組成了其他不同形狀的對(duì)稱型超構(gòu)材料, 驗(yàn)證不同的對(duì)稱結(jié)構(gòu)是否存在聲場(chǎng)增強(qiáng)效應(yīng).結(jié)構(gòu)設(shè)計(jì)分別如圖4(a)—(c)所示, 利用不同數(shù)量的多腔型基本單元設(shè)計(jì)了外部形狀分別為圓形、矩形、正六邊形的對(duì)稱型超構(gòu)材料.其中圓形對(duì)稱型超構(gòu)材料的半徑為d1= 18 cm, 矩形和正六邊形對(duì)稱型超構(gòu)材料中相鄰兩個(gè)基本單元中心點(diǎn)之間的距離分別為d2= 16 cm和d3= 20 cm, 并且這些多腔型基本單元的結(jié)構(gòu)參數(shù)與圖1(a)中的完全一致.對(duì)于這3個(gè)系統(tǒng), 圓形、矩形和正六邊形超構(gòu)材料對(duì)應(yīng)的復(fù)合單極子Mie共振模式所對(duì)應(yīng)的共振頻率分別為648, 629和621 Hz.在共振頻率處, 當(dāng)聲波以1 Pa的強(qiáng)度從左側(cè)垂直入射到這3個(gè)超構(gòu)材料系統(tǒng)時(shí)的聲強(qiáng)分布分別如圖4(d)—(f)所示, 聲場(chǎng)增強(qiáng)效應(yīng)均存在于這3個(gè)對(duì)稱型超構(gòu)材料的中心區(qū)域, 其內(nèi)部中心點(diǎn)處的聲場(chǎng)強(qiáng)度分別達(dá)到了15,20和20 Pa2.由此可以推斷, 由該多腔型基本單元構(gòu)成的不同形狀的對(duì)稱型超構(gòu)材料均可以實(shí)現(xiàn)聲場(chǎng)增強(qiáng)效應(yīng), 這種靈活可調(diào)的聲場(chǎng)增強(qiáng)系統(tǒng)可應(yīng)用于許多特定場(chǎng)景, 例如聲通信、聲源或噪聲定位等.
圖4 由多腔型基本單元組成的(a)圓形、(b)矩形和(c)正六邊形對(duì)稱型超構(gòu)材料的結(jié)構(gòu)示意圖; 在復(fù)合單極子Mie共振模式下, (d)圓形超構(gòu)材料、(e)矩形超構(gòu)材料和(f)正六邊形超構(gòu)材料的聲強(qiáng)分布圖, 對(duì)應(yīng)的工作頻率分別為648, 629, 621 HzFig.4.Schematic diagram of the (a) circular, (b) rectangular and (c) regular hexagonal symmetric metamaterials composed of multiple-cavity units.Sound intensity distributions of the (d) circular metamaterial, (e) rectangular metamaterial and (f) regular hexagonal metamaterial under the compound monopole Mie resonance mode with the working frequencies of 648, 629 and 621 Hz,respectively.
除了利用該多腔型基本單元設(shè)計(jì)對(duì)稱型超構(gòu)材料, 本文還利用該基本單元構(gòu)造了左右非對(duì)稱超構(gòu)材料.如圖5(a)所示, 設(shè)計(jì)了一個(gè)形狀為等邊三角形的左右非對(duì)稱型的超構(gòu)材料, 在此結(jié)構(gòu)中相鄰兩個(gè)基本單元中心點(diǎn)之間的距離均為d4= 28 cm.如圖5(b)所示, 頻率為637 Hz的平面波以1 Pa的強(qiáng)度從左側(cè)垂直入射, 該非對(duì)稱型超構(gòu)材料實(shí)現(xiàn)了聲場(chǎng)的局部增強(qiáng)效應(yīng), 聲能集中在系統(tǒng)內(nèi)部的中心區(qū)域, 中心點(diǎn)處聲強(qiáng)可高達(dá)100 Pa2, 實(shí)現(xiàn)了較強(qiáng)的聲場(chǎng)增強(qiáng)效應(yīng).為了探究該非對(duì)稱型超構(gòu)的聲場(chǎng)增強(qiáng)效應(yīng)是否受入射聲波角度的影響, 繪制了結(jié)構(gòu)內(nèi)部點(diǎn)A處的聲強(qiáng)隨入射聲波角度的變化情況.如圖5(c)所示, 點(diǎn)A處的聲強(qiáng)分布關(guān)于0°對(duì)稱,當(dāng)入射聲波的角度趨于—90°和90°時(shí), 點(diǎn)A處的聲強(qiáng)趨于0 Pa2.而當(dāng)聲波從左側(cè)垂直入射時(shí)( θ =0°),點(diǎn)A處的聲強(qiáng)達(dá)到最大值, 隨著角度的逐漸增大,點(diǎn)A處的聲強(qiáng)也隨之減小, 因此, 該非對(duì)稱型超構(gòu)材料內(nèi)部的聲場(chǎng)增強(qiáng)效應(yīng)與入射聲波的角度有密切關(guān)系.從圖5(a)的幾何圖形可以發(fā)現(xiàn), 當(dāng)平面波沿著該非對(duì)稱超構(gòu)材料左側(cè)以60°入射時(shí), 相當(dāng)于平面波從右側(cè)垂直入射.結(jié)合圖5(b)的聲強(qiáng)分布可知, 平面波從左側(cè)垂直入射的情況下, 結(jié)構(gòu)內(nèi)部中心點(diǎn)處的聲強(qiáng)約為100 Pa2, 而由圖5(d)可知,當(dāng)平面波從右側(cè)垂直入射的情況下, 結(jié)構(gòu)內(nèi)部中心點(diǎn)處的聲強(qiáng)僅為10 Pa2.因此, 設(shè)計(jì)的非對(duì)稱超構(gòu)材料對(duì)側(cè)向垂直入射聲波實(shí)現(xiàn)了高效的非對(duì)稱聲場(chǎng)增強(qiáng)效應(yīng), 且通過調(diào)節(jié)該等邊三角形的傾角可對(duì)聲場(chǎng)增強(qiáng)效應(yīng)進(jìn)行調(diào)控.
圖5 (a) 6個(gè)多腔型基本單元組成的等邊三角形非對(duì)稱型超構(gòu)材料的結(jié)構(gòu)示意圖; (b) 當(dāng)頻率為637 Hz的聲波從結(jié)構(gòu)左側(cè)垂直入射時(shí), 非對(duì)稱超構(gòu)材料的聲強(qiáng)分布圖; (c) 非對(duì)稱超構(gòu)材料內(nèi)部點(diǎn)A處的聲強(qiáng)隨入射聲波角度的變化情況; (d) 當(dāng)頻率為637 Hz的聲波從結(jié)構(gòu)右側(cè)垂直入射時(shí), 非對(duì)稱超構(gòu)材料的聲強(qiáng)分布圖Fig.5.(a) Structural diagram of the equilateral triangular asymmetric metamaterial composed of six multiple-cavity units;(b) sound intensity distribution of the asymmetric metamaterial when the sound wave with frequency of 637 Hz is normally incident from the left side; (c) acoustic intensity curve at the center point A of the asymmetric metamaterial with different angles of incident wave; (d) sound intensity distribution of the asymmetric metamaterial when the sound wave is normally incident from the right side with frequency of 637 Hz.
通過對(duì)該非對(duì)稱型超構(gòu)材料的研究, 可以發(fā)現(xiàn), 聲波入射方向不同, 聲場(chǎng)增強(qiáng)效應(yīng)也不同.對(duì)于本文設(shè)計(jì)的左右非對(duì)稱系統(tǒng)而言, 當(dāng)聲波從結(jié)構(gòu)左側(cè)垂直入射時(shí), 能夠在系統(tǒng)內(nèi)部實(shí)現(xiàn)聲場(chǎng)增強(qiáng)效應(yīng), 而當(dāng)聲波從結(jié)構(gòu)右側(cè)垂直入射時(shí), 其中心點(diǎn)處的聲強(qiáng)將大大減小, 聲能聚集于各個(gè)多腔型基本單元中.由此可知, 由該基本單元所設(shè)計(jì)的非對(duì)稱型超構(gòu)材料, 具有較強(qiáng)的方向依賴性, 利用該特性通過設(shè)計(jì)由多腔型基本單元組成的非對(duì)稱結(jié)構(gòu)可實(shí)現(xiàn)非對(duì)稱聲場(chǎng)增強(qiáng)效應(yīng), 進(jìn)而實(shí)現(xiàn)隔聲、聲傳感等方面的應(yīng)用.
本文利用多腔型基本單元構(gòu)造了不同結(jié)構(gòu)的超構(gòu)材料, 研究聲波入射方向、激發(fā)頻率、結(jié)構(gòu)對(duì)稱性等因素對(duì)聲場(chǎng)增強(qiáng)效應(yīng)的影響.結(jié)果表明, 不同的結(jié)構(gòu)在外部聲波激勵(lì)下, 有著不同效果的聲場(chǎng)增強(qiáng)效應(yīng).對(duì)稱型超構(gòu)材料的聲場(chǎng)增強(qiáng)效應(yīng)不受入射方向的影響, 這一特性使其在聲通信、聲源或噪聲定位等方面有重要的應(yīng)用價(jià)值.非對(duì)稱型超構(gòu)材料具有較強(qiáng)的方向依賴性, 聲波入射方向不同, 系統(tǒng)中心點(diǎn)處的聲強(qiáng)也不同.利用該多腔基本單元組成非對(duì)稱型超構(gòu)材料能夠?qū)崿F(xiàn)非對(duì)稱聲場(chǎng)增強(qiáng), 甚至單向聲場(chǎng)增強(qiáng), 這一特性使其在隔聲、聲傳感等方面具有潛在的應(yīng)用價(jià)值.