超聲激勵(lì)薄液膜Faraday波形成機(jī)理

高國富，李康，李瑜，向道輝，趙波

（河南理工大學(xué) 機(jī)械與動(dòng)力工程學(xué)院，焦作454000）

隨著金屬粉末注射成型技術(shù)、球柵電子封裝技術(shù)以及金屬3D打印制造技術(shù)等先進(jìn)粉末冶金技術(shù)的發(fā)展，高性能金屬粉末的制備受到了廣泛關(guān)注［1］。當(dāng)前金屬粉末的主要制備方法有等離子旋轉(zhuǎn)電極法、等離子霧化法、均勻顆粒成型法和氣霧化法等［1-2］。其中均勻顆粒成型法和氣霧化法應(yīng)用最為廣泛。氣霧化法存在噴嘴結(jié)構(gòu)復(fù)雜、收得率低、衛(wèi)星粉形成機(jī)率大等缺點(diǎn)［3］。均勻顆粒成型法因受到噴嘴直徑的影響，較難制備直徑小于100μm的金屬顆粒［4］。超聲輔助靜電噴涂概念的提出為微納米顆粒的獲取提供了新的解決方案。超聲輔助靜電噴涂的原理是薄液膜受到高頻超聲振動(dòng)，在液膜表面形成大量微細(xì)駐波。同時(shí)，在靜電場的作用下，從波峰處分離出直徑均勻的帶電液滴，實(shí)現(xiàn)獲得均勻顆粒目的［5］。

液體受到周期性垂直激勵(lì)形成的振動(dòng)頻率為超聲激勵(lì)頻率1／2的液體表面波統(tǒng)稱為Faraday波［6］。Benjamin和Ursell［7］推導(dǎo)了圓柱形容器中自由表面波運(yùn)動(dòng)的振幅方程，并用穩(wěn)定性理論得出不同條件下表面波振動(dòng)頻率可為超聲激勵(lì)頻率的1／2、1或者3／2倍。Eisenmenger［8］基于線性流變振蕩理論對(duì)Faraday波的形成進(jìn)行理論分析和實(shí)驗(yàn)驗(yàn)證，測量了不同超聲激勵(lì)頻率下毛細(xì)波波長及形成毛細(xì)波的初始激勵(lì)振幅，給出激勵(lì)振幅與毛細(xì)波振幅的關(guān)系。Peskin和Raco［9］對(duì)薄液膜進(jìn)行線性穩(wěn)定性分析，導(dǎo)出了超聲霧化的霧粒直徑與超聲激勵(lì)頻率、換能器振幅和薄液膜厚度之間的關(guān)系。Sindayihebura等［10-11］分析了超聲激勵(lì)頻率、振幅、平均膜厚、液體黏度等對(duì)駐波的影響。菅永軍等［12-14］忽略表面張力的影響，用兩變量時(shí)間展開法得到具有立方項(xiàng)以及底部驅(qū)動(dòng)項(xiàng)影響的非線性振幅方程。當(dāng)受迫振動(dòng)頻率較低時(shí)，表面張力可忽略，但在高頻條件下表面張力不可忽略。Higginbotham 等［15］通過實(shí)驗(yàn)觀察到低頻振蕩器產(chǎn)生與亞諧波Faraday波相關(guān)的液滴，高頻振蕩器產(chǎn)生與諧波響應(yīng)一致的液滴；中頻振蕩器可根據(jù)驅(qū)動(dòng)振幅激發(fā)任何一種響應(yīng)的液滴。在高頻振蕩器中，可以通過增大驅(qū)動(dòng)振幅改變霧化液滴的尺寸。Li和Umemura［16］基于不可壓縮歐拉方程對(duì)Faraday波不穩(wěn)定性進(jìn)行了兩相流仿真，從物理上探討了液體韌帶是如何從振動(dòng)的液體層脫離。劉財(cái)興等［17］從實(shí)驗(yàn)與理論兩方面研究了垂直激勵(lì)低黏度硅油在低頻范圍內(nèi)Faraday波的特性，發(fā)現(xiàn)隨著超聲激勵(lì)頻率增大，駐波圖案的波峰數(shù)增加，模態(tài)復(fù)雜化，對(duì)應(yīng)的驅(qū)動(dòng)振幅減小。Sheldrake等［18］系統(tǒng)地探索了振幅、頻率、體積（或深度）、溫度和大氣壓力對(duì)垂直壁圓形容器中Faraday波形的影響。趙文定等［19］由Faraday波振幅方程導(dǎo)出了穩(wěn)定條件下的色散關(guān)系，利用參量共振方程得到了在亞簡諧條件下本征波矢的取值范圍。Tadrist等［20］導(dǎo)出了由垂直振動(dòng)液體表面產(chǎn)生的Faraday表面波模型。

綜上所述，F(xiàn)araday波多在較低的激勵(lì)頻率下研究，對(duì)超聲激勵(lì)下Faraday波的理論和實(shí)驗(yàn)研究較少。Faraday波的波長與激勵(lì)頻率有關(guān)，通過增大激勵(lì)頻率可得到微納尺度波長的表面波，進(jìn)而有利于獲得微納顆粒。本文針對(duì)超聲激勵(lì)頻率為35 kHz條件下液膜表面形成Faraday波進(jìn)行有限元仿真和實(shí)驗(yàn)，對(duì)Faraday波的形成機(jī)理進(jìn)行探究。

1 CFD有限元仿真分析

對(duì)薄液膜施加垂直其表面的超聲振動(dòng)，在液膜表面激發(fā)排列整齊的正方形表面波陣列，即Faraday波。通過CFD有限元仿真對(duì)Faraday波的形成機(jī)理進(jìn)行探究，并對(duì)Faraday波的形成過程及相關(guān)條件進(jìn)行討論。

1.1 計(jì)算模型

表面波的幾何形狀與液體體積有關(guān)，而較少的液體體積不影響仿真結(jié)果。三維計(jì)算模型雖然直觀，但考慮到計(jì)算時(shí)長和硬件限制，二維計(jì)算模型時(shí)間成本較低，且不影響仿真結(jié)果，故選擇二維計(jì)算模型。

該計(jì)算模型采用超聲激勵(lì)頻率為35 kHz，所用液體為蒸餾水（室溫為20℃）。蒸餾水的表面張力為0.072 N／m，動(dòng)力黏度為0.001 N·s，密度為1 000 kg／m3［21］。

表面波波長由開爾文公式［22］計(jì)算：

式中：l為表面波波長；σ為液體表面張力；ρ為液體密度；f為超聲激勵(lì)頻率。

將相關(guān)數(shù)據(jù)代入式（1）可知，35 kHz超聲激勵(lì)水薄膜形成表面波的波長約為113.87μm，因此模型長度設(shè)置為1mm，約為理論波長的8.7倍，液膜厚度為0.05mm。為研究薄液膜表面波的運(yùn)動(dòng)情況，設(shè)置液體域上方的空氣域高度為液膜厚度的2倍，即為0.1mm。因此，CFD有限元仿真二維計(jì)算模型由2個(gè)矩形組成，分別為上方空氣域矩形及下方流體域矩形?？諝庥蚓匦伍L度為1mm，寬度為0.1 mm，流體域矩形長度為1 mm，寬度為0.05mm［23］。

在該二維計(jì)算模型中，以非結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算域進(jìn)行網(wǎng)格劃分，考慮到表面波波長和液膜底部側(cè)壁及空氣域上邊界對(duì)表面波影響較小，液膜底部和側(cè)壁一個(gè)網(wǎng)格的長度設(shè)置為10μm，空氣域上邊界一個(gè)網(wǎng)格長度設(shè)置為10μm，空氣域兩側(cè)邊界和流體域上邊界對(duì)表面波的形成影響較大，因此對(duì)這兩部分網(wǎng)格進(jìn)行加密處理，設(shè)置一個(gè)網(wǎng)格長度為5μm，使每個(gè)表面波波長方向上劃分有20個(gè)網(wǎng)格，以保證計(jì)算的有效性［24］。計(jì)算模型如圖1所示。

圖1 計(jì)算模型Fig.1 Calculation model

采用多相流模型（VOF）對(duì)氣液交界面的運(yùn)動(dòng)情況進(jìn)行追蹤顯示。在計(jì)算模型中，除空氣域的上邊界條件為壓力出口（pressure-out）之外，其余的邊界條件均為壁（wall），液體域中的側(cè)壁無剪切滑移作用［25］，具體如圖1所示。通過動(dòng)網(wǎng)格及用戶自定義函數(shù)（UDF）功能，在流體域下邊界施加超聲振動(dòng)［26］。設(shè)施加超聲振動(dòng)位移為

式中：ω為超聲激勵(lì)角頻率，；A為超聲激勵(lì)振幅；t為時(shí)間。

流體域底部作簡諧運(yùn)動(dòng)，其軸向速度為

通過相圖和流線圖，對(duì)超聲激勵(lì)下液體表面波的形成過程進(jìn)行分析。

1.2 仿真結(jié)果與分析

利用動(dòng)網(wǎng)格和 UDF功能對(duì)流體域施加35 kHz超聲振動(dòng)。圖2所示為超聲激勵(lì)頻率35 kHz、超聲激勵(lì)振幅8μm條件下不同時(shí)刻薄液膜的氣液相圖。

圖2 超聲激勵(lì)頻率35 kHz、超聲激勵(lì)振幅8μm條件下不同時(shí)刻的Faraday波相圖Fig.2 Phase diagram of Faraday wave at differentmoments under ultrasonic excitation frequency 35 kHz and ultrasonic excitation amplitude 8μm

從圖2中可知，當(dāng)施加超聲激勵(lì)后，薄液膜表面從靜止?fàn)顟B(tài)開始振動(dòng)，經(jīng)歷一段時(shí)間（7.2×10－5s），在液膜表面中心位置存在振幅較小的表面波動(dòng)；隨著時(shí)間推移（1.3×10－4s），液膜表面波動(dòng)較明顯，振幅不斷增大；當(dāng)經(jīng)歷多個(gè)周期超聲激勵(lì)（1.9×10－4s）后，液膜表面波動(dòng)現(xiàn)象逐漸趨于穩(wěn)定，振幅不再增加，波峰波谷隨著時(shí)間作交替變換，波節(jié)位置基本保持不變。

圖3為超聲激勵(lì)頻率35 kHz時(shí)，不同超聲激勵(lì)振幅的氣液相圖。如圖3所示，當(dāng)施加超聲激勵(lì)振幅為2、4μm時(shí)，表面波動(dòng)不明顯；當(dāng)施加超聲激勵(lì)振幅增加為6、8、10μm時(shí)，薄液膜表面出現(xiàn)較明顯的均勻表面波，且中間位置表面波一致性好；當(dāng)施加超聲激勵(lì)振幅為12μm時(shí)，表面波振幅較大，形成脖子，表面波波峰脖子以上液體有脫離液膜主體的趨勢(shì)，即超聲霧化即將發(fā)生。

由圖3可知，當(dāng)超聲激勵(lì)振幅超過某一臨界值α?xí)r，在液膜表面產(chǎn)生表面波；超聲振幅增加，表面波振幅隨之增加；當(dāng)超聲激勵(lì)振幅超過臨界值β時(shí)，薄液膜表面的表面波波峰處形成脖子，脖子以上液體不再恢復(fù)到薄液膜主體，即發(fā)生超聲霧化現(xiàn)象。在不同超聲激勵(lì)振幅下，表面波波長基本無變化，即某一種液體在同一超聲激勵(lì)頻率下，改變超聲激勵(lì)振幅僅改變表面波振幅，對(duì)表面波波長無影響。

圖4所示為施加超聲激勵(lì)頻率35 kHz，超聲激勵(lì)振幅8μm條件下，1.9×10－4s時(shí)刻的氣液相圖，表1為圖4中各個(gè)波峰位置處的坐標(biāo)值及形成的表面波波長。表面波波長為后一波峰坐標(biāo)位置與前一波峰位置坐標(biāo)的差值，最終得出該時(shí)刻表面波平均波長。

圖3 超聲激勵(lì)頻率35 kHz條件下不同超聲激勵(lì)振幅的Faraday波相圖Fig.3 Phase diagram of Faraday wave under different ultrasonic excitation amplitudes at ultrasonic excitation frequency 35 kHz

圖4 超聲激勵(lì)頻率35 kHz、超聲激勵(lì)振幅8μm條件下1.9×10－4 s時(shí)刻的Faraday波相圖Fig.4 Phase diagram of Faraday wave at 1.9×10－4 s under ultrasonic excitation frequency 35 kHz and ultrasonic excitation amplitude 8μm

從圖4可知，薄液膜表面形成了結(jié)構(gòu)一致的表面波。在薄液膜表面有9個(gè)波峰，第1個(gè)波峰和第9波峰的方向不是豎直向上，這可能是由于表面張力的存在使處于邊緣處的液膜有一定坡度，造成波峰方向有偏差，其余波峰方向豎直向上，且相鄰波峰之間長度基本一致。通過取點(diǎn)測量，各個(gè)波峰橫向坐標(biāo)如表1所示，9個(gè)波峰之間的平均距離為109.74μm，即表面波平均波長為109.74μm，與式（1）計(jì)算的理論值113.87μm的誤差為3.6%。

圖5為一個(gè)超聲激勵(lì)周期表面波的流線圖。其中，T為一個(gè)超聲激勵(lì)周期，ti為仿真經(jīng)歷時(shí)間。超聲激勵(lì)頻率為35 kHz，超聲激勵(lì)振幅為8μm。從流線圖可知：

1）遠(yuǎn)離薄液膜的空氣域速度方向始終和超聲激勵(lì)（液膜底面）的速度方向保持一致，薄液膜附近空氣域的速度方向有偏差，此部分空氣域速度變化較快，液體域速度存在滯后現(xiàn)象。

表1 圖4中各個(gè)波峰位置處的坐標(biāo)及波長Tab le 1 Coord inates and wavelength at each peak position in Fig.4

2）超聲激勵(lì)豎直向下的速度逐漸減小轉(zhuǎn)向?yàn)樨Q直向上并增大到最大值的過程中（圖5（a）、（b）和（f）），波峰左側(cè)空氣域和流體域的速度矢量組成一個(gè)順時(shí)針閉環(huán)，波峰右側(cè)空氣域和流體域的速度矢量組成一個(gè)逆時(shí)針閉環(huán)。表面波波峰和波谷先短暫增加然后衰減至液膜表面基本為平面，空氣抑制表面波增長。

3）超聲激勵(lì)豎直向上的速度由最大值逐漸減小轉(zhuǎn)向并增加到豎直向下最大值的過程中（圖5（c）～（e）），波峰左側(cè)空氣域和流體域的速度矢量組成一個(gè)逆時(shí)針閉環(huán)，波峰右側(cè)空氣域和流體域的速度矢量組成一個(gè)順時(shí)針閉環(huán)，空氣促進(jìn)波峰波谷增長。

4）圖5中x＝0.5mm處表面波在（a）中為波峰，在（f）中為波谷，即一個(gè)超聲激勵(lì)周期約為液膜表面波運(yùn)動(dòng)周期的1／2，超聲激勵(lì)與表面波運(yùn)動(dòng)存在不斷變化的相位差，一個(gè)相位差變化周期約為2個(gè)超聲激勵(lì)周期。圖6所示為圖5中x＝0.5mm處表面波與超聲激勵(lì)的相位變化，圖中η為位移量。

表面波是外加周期性垂直激勵(lì)后流體域和空氣域相互作用形成的。超聲激勵(lì)豎直向下的速度逐漸減小轉(zhuǎn)向?yàn)樨Q直向上且不斷增加到最大值的過程中，由于液體域慣性較大，液體域波峰處液體保持原來的運(yùn)動(dòng)趨勢(shì)，波峰處空氣域速度方向隨超聲激勵(lì)的變化而變化，因此在波峰左側(cè)形成順時(shí)針?biāo)俣仁噶块]環(huán)，在波峰右側(cè)形成逆時(shí)針?biāo)俣仁噶块]環(huán)，此過程中空氣域抑制波峰和波谷的增長；同理，當(dāng)超聲激勵(lì)豎直向上的速度由最大值逐漸減小轉(zhuǎn)向并增加到最大值的過程中，在波峰左側(cè)形成逆時(shí)針矢量閉環(huán)，在波峰右側(cè)形成順時(shí)針矢量閉環(huán)，空氣域促進(jìn)波峰波谷的增長。

2 實(shí)驗(yàn)分析

為測試超聲激勵(lì)表面波的形成狀態(tài)，實(shí)驗(yàn)采用35 kHz超聲波發(fā)生器和圓錐形變幅桿組成超聲振動(dòng)系統(tǒng)，通過超景深顯微鏡（VHX-2000）觀察變幅桿發(fā)射端薄液膜表面的運(yùn)動(dòng)情況，實(shí)驗(yàn)裝置如圖7所示。

將工作液（蒸餾水）均勻涂覆在變幅桿發(fā)射端。通過超聲振動(dòng)系統(tǒng)對(duì)液膜施加縱向振動(dòng)，激勵(lì)液膜表面形成Faraday波。在超聲振動(dòng)系統(tǒng)上方搭建超景深顯微鏡，觀測液膜表面狀態(tài)。圖8為通過超景深顯微鏡觀察到的薄液膜表面超聲激勵(lì)產(chǎn)生的Faraday波陣列。

在圖8中，液膜表面存在均勻分布點(diǎn)陣，其中亮點(diǎn)為表面波波峰或波谷，彼此相間，表面波波節(jié)位置基本不變，2個(gè)亮點(diǎn)之間的距離為表面波波長的1／2，3個(gè)亮點(diǎn)之間的距離為一個(gè)表面波波長。通過測量，3個(gè)亮點(diǎn)之間距離的平均值約為117.6μm。因此，實(shí)驗(yàn)獲得表面波的振動(dòng)頻率約為超聲激勵(lì)頻率的1／2，理論誤差為3.3%。實(shí)驗(yàn)測得的表面波波長大于理論值，這可能是由開爾文公式未考慮重力波的影響和測量誤差的存在造成。

圖7 實(shí)驗(yàn)裝置平臺(tái)Fig.7 Experimental device platform

圖8 表面駐波陣列Fig.8 Surface standing wave array

3 結(jié) 論

通過實(shí)驗(yàn)和CFD有限元仿真分析，對(duì)薄液膜在超聲激勵(lì)下形成的表面波進(jìn)行了研究。

1）通過CFD有限元仿真分析，同一超聲激勵(lì)頻率下，當(dāng)超聲激勵(lì)振幅超過臨界值α?xí)r，經(jīng)歷數(shù)個(gè)超聲激勵(lì)周期后，薄液膜表面激發(fā)形成Faraday波。隨超聲激勵(lì)振幅增加，F(xiàn)araday波振幅增大，當(dāng)超聲激勵(lì)振幅超過臨界值β時(shí)，會(huì)產(chǎn)生超聲霧化現(xiàn)象。此外，表面波波長不隨超聲激勵(lì)振幅的改變而改變。

2）通過對(duì)一個(gè)超聲激勵(lì)周期表面波流線圖分析，在超聲激勵(lì)豎直向下的速度逐漸減小轉(zhuǎn)向?yàn)樨Q直向上且不斷增加到最大值的過程中，波峰處存在速度矢量閉環(huán)，空氣抑制表面波的增長；當(dāng)超聲激勵(lì)豎直向上的速度由最大值逐漸減小并轉(zhuǎn)向增大到豎直向下最大值的過程中，波峰處存在速度矢量閉環(huán)，空氣促進(jìn)表面波的生長。液膜表面的表面波是在外部周期性激勵(lì)下由空氣與液體相互作用形成的。表面波波峰或波谷與超聲激勵(lì)存在不斷變化的相位差，一個(gè)相位差變化周期約為2個(gè)超聲激勵(lì)周期。

3）超聲振動(dòng)激勵(lì)薄液膜表面形成Faraday波，實(shí)驗(yàn)測量的Faraday波波長與理論計(jì)算值相比較，得出35 kHz超聲激勵(lì)形成Faraday波的振動(dòng)頻率約為超聲激勵(lì)頻率的1／2，與CFD有限元仿真分析結(jié)果一致。