劉作平, 周 健, 裴澤光, b, 白 濤, 陳廣鋒
(東華大學(xué)a.機(jī)械工程學(xué)院;b.紡織裝備教育部工程研究中心,上海201620)
基于微滴噴射技術(shù)的噴墨印刷按照微滴形成的方式主要分為連續(xù)式與按需式[1]。其中壓電驅(qū)動(dòng)的按需式噴墨印刷技術(shù)利用壓電材料的逆壓電效應(yīng)產(chǎn)生驅(qū)動(dòng)力從而噴射出微滴,其由于具有響應(yīng)速度快、效率高、可實(shí)現(xiàn)按需噴射等優(yōu)點(diǎn)而被廣泛應(yīng)用[2]。壓電式噴墨過程包含復(fù)雜的流體形態(tài)變化,因此通過數(shù)值模擬研究微滴的成形過程,有助于明晰微滴噴射的內(nèi)在機(jī)制,這對(duì)基于壓電式噴墨印刷系統(tǒng)的技術(shù)改進(jìn)具有重要意義。
近年來(lái),為滿足噴墨印刷精度與速度的要求,微滴噴射技術(shù)不斷向精密化和高速高效化發(fā)展。薛光懷等[3]設(shè)計(jì)并加工一種壓電式微滴噴頭,并在微滴噴射平臺(tái)上進(jìn)行微滴噴射試驗(yàn),通過改變高電平時(shí)間和增加壓電片的自由振蕩衰減阻值可提高微滴噴射的最高頻率;劉趙淼等[4]設(shè)計(jì)一種用于生成均勻微滴的壓電驅(qū)動(dòng)式微滴噴射裝置,采用數(shù)值模擬和試驗(yàn)相結(jié)合的方法,研究不同控制參數(shù)下的膜片振幅及其對(duì)生成的微滴的尺寸和均勻性的影響;周健等[5]設(shè)計(jì)并制造一種利用同軸氣流對(duì)微滴飛行過程進(jìn)行驅(qū)動(dòng)的同軸氣流輔助式壓電微滴噴頭。目前國(guó)內(nèi)外研究者對(duì)于同軸氣流霧化[6-7]、氣-液流動(dòng)聚焦[8-9]以及熔噴[10]等液體射流在同軸氣流作用下的形成與流動(dòng)特性已做了較多的研究與報(bào)道,但是對(duì)液滴在壓電膜片與同軸氣流共同驅(qū)動(dòng)下的形成規(guī)律研究得相對(duì)較少。本文利用多物理場(chǎng)計(jì)算軟件COMSOL Multiphysics對(duì)同軸氣流作用下壓電式微滴噴射過程和流體流動(dòng)規(guī)律進(jìn)行二維數(shù)值模擬,以明晰同軸氣流作用下微滴的形成與運(yùn)動(dòng)特征。
COMSOL Multiphysics是一個(gè)專業(yè)的基于偏微分方程對(duì)多物理場(chǎng)模型進(jìn)行建模和仿真的科學(xué)仿真軟件,其通過有限元方法模擬求解在科研和工程中能用偏微分方程(partial differential equation, PDE)描述的各種問題[11-14]。圖1為同軸氣流輔助式壓電微滴噴頭的結(jié)構(gòu)示意圖。為降低模型求解的計(jì)算量,考慮到進(jìn)液通道與進(jìn)氣通道建模的特殊性,將設(shè)有同軸氣流噴射槽的壓電噴頭中的微滴噴射過程由一個(gè)三維立體問題簡(jiǎn)化為二維平面模型。模擬計(jì)算時(shí)不考慮支撐體底座與上蓋以及密封圈等部件的影響,整個(gè)噴頭計(jì)算區(qū)域如圖2所示,其中儲(chǔ)液腔錐角為120°。噴頭關(guān)鍵結(jié)構(gòu)參數(shù)如表1所示。
圖1 設(shè)有同軸氣流噴射槽的壓電式微滴噴頭結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the airflow-assisted piezoelectric printhead
圖2 噴頭計(jì)算區(qū)域示意圖Fig.2 Computational domain of the printhead
表1 噴頭主要結(jié)構(gòu)參數(shù)Table 1 Main structural parameters of the printhead mm
模擬使用的液體為水(流體1),密度為1 000 kg/m3,動(dòng)力黏度為0.001 Pa·s,表面張力為72.75×10-3N/m,液體為不可壓縮流;噴嘴外部區(qū)域及空氣域中的流體為空氣(流體2),密度為1.225 kg/m3,動(dòng)力黏度為1.789×10-5Pa·s。在設(shè)有同軸氣流噴射槽的壓電式噴頭中壓電致動(dòng)器采用圓片薄膜式結(jié)構(gòu),由位于下方的薄銅片和上方的壓電陶瓷片黏結(jié)而成。其中:薄銅片的密度為8 960 kg/m3,彈性模量為1.1×1011Pa,泊松比為0.35;所用的壓電陶瓷材料為COMSOL材料庫(kù)中的PZT-5H,與實(shí)際壓電陶瓷材料相符,密度為7 500 kg/m3。其他主要參數(shù)包括剛度矩陣CE、介電常數(shù)矩陣λ以及壓電耦合矩陣CB[15],具體數(shù)值如式(1)~(3)所示。
(1)
(2)
(3)
壓電材料的壓電方程為
(4)
式中:S1、T1分別為水平方向的應(yīng)變和應(yīng)力;E3、D3分別為豎直方向的電場(chǎng)和電應(yīng)力;λ為外應(yīng)力為0或定值時(shí)的介電矩陣;d31為壓電常數(shù);CE為電場(chǎng)為0或恒定值時(shí)測(cè)得的壓電陶瓷剛度矩陣。
采用層流兩相流水平集的方法,將流體視為不可壓縮流,采用Navier-Stokes方程[16]對(duì)噴腔體內(nèi)與噴嘴內(nèi)的流體進(jìn)行模擬分析,如式(5)~(6)所示。
Fg+Fst+Fext+F
(5)
(6)
式中:p為壓力;I為單位張量;Fg為噴腔體內(nèi)流體的重力;Fst為兩種流體界面的表面張力;Fext為外界自由能引起的體積力;F為作用在壓電陶瓷片上的外力;ρ為密度,其是水平集方程的函數(shù),如式(7)所示;μ為動(dòng)力黏度,如式(8)所示。
ρ=ρ2+(ρ1-ρ2)φ
(7)
μ=μ2+(μ1-μ2)φ
(8)
式中:ρ1、ρ2分別為流體1與流體2的密度;μ1、μ2分別為流體1與流體2的動(dòng)力黏度;φ為水平集函數(shù);流體1對(duì)應(yīng)于φ>0.5的區(qū)域,流體2對(duì)應(yīng)于φ<0.5的區(qū)域。
兩相流交界面的表面張力Fst=σkδn,其中,σ為表面張力系數(shù),k為曲率,δ為集中于交界面的狄拉克函數(shù),n為交界面單位向量。k取決于水平集函數(shù)φ的二次派生,這會(huì)引起表面張力精度的降低,因此根據(jù)變形公式(9)來(lái)計(jì)算Fst。
(9)
Fext在使用相場(chǎng)交界面時(shí)需要被計(jì)算到方程中,其表達(dá)式如式(10)所示。
(10)
式中:?f/?φ為外界自由能源。
采用COMSOL有限元仿真軟件建立設(shè)有同軸氣流噴射槽的壓電噴頭的二維多物理場(chǎng)耦合模型,并添加兩相流水平集、流固耦合與固體力學(xué)物理場(chǎng)。其中:空氣域以及水域添加兩相流水平集多物理場(chǎng),并設(shè)置流體流動(dòng)屬性為層流。將薄銅片與噴腔體的接觸面設(shè)置為流固耦合界面,耦合類型選擇全耦合。固體域選擇為壓電陶瓷片與薄銅片,并將模型中的壓電陶瓷片設(shè)置為靜電物理場(chǎng),用以模擬壓電驅(qū)動(dòng)過程,施加在壓電片上的驅(qū)動(dòng)電壓波形如圖3所示,電壓波形參數(shù)如表2所示。此外,將壓電陶瓷片下邊界設(shè)置為接地屬性以進(jìn)行靜電場(chǎng)仿真。
圖3 驅(qū)動(dòng)電壓波形Fig.3 Waveform of driving voltage
表2 電壓波形參數(shù)Table 2 Parameters of voltage waveform
帶同軸氣流的壓電式微滴噴頭網(wǎng)格劃分圖如圖4所示。其中:壓電陶瓷片與銅片的網(wǎng)格類型為自由剖分四邊形網(wǎng)格,最大單元大小設(shè)置為0.05 mm;其他部分為自由三角形網(wǎng)格,設(shè)置網(wǎng)格大小為極細(xì)化。使用COMSOL軟件中的mesh refine功能對(duì)模型的關(guān)鍵計(jì)算區(qū)域壓電片及噴嘴部分進(jìn)行局部網(wǎng)格加密,其他區(qū)域保持相對(duì)較大的網(wǎng)格密度。每組仿真計(jì)算時(shí)都需要多次調(diào)整網(wǎng)格的數(shù)量以確保求解精度,最終調(diào)整得到每組仿真的網(wǎng)格數(shù)量為30萬(wàn)~40萬(wàn)。噴腔體進(jìn)液口初始條件為邊界無(wú)壓力并且設(shè)置為自動(dòng)進(jìn)液且抑制回流;進(jìn)氣口邊界條件為作用于壓電式噴頭的同軸氣流的法向速度值;兩相流交界面為噴嘴與空氣接觸面;微滴噴射出口為外界環(huán)境邊界,邊界條件設(shè)置為無(wú)壓力。將儲(chǔ)液腔內(nèi)壁、進(jìn)液通道內(nèi)壁面、噴嘴內(nèi)外壁面與端面、進(jìn)氣通道內(nèi)壁面設(shè)置為固壁,分析模型為瞬態(tài)分析,模型的仿真總時(shí)長(zhǎng)為15 000 μs,仿真時(shí)間步長(zhǎng)為250 μs。
圖4 噴頭網(wǎng)格劃分情況Fig.4 Grid generation of the printhead
圖5為同軸氣流流速v=0 m/s時(shí)一個(gè)噴射周期內(nèi)的速度分布圖,其中位于噴頭左側(cè)的圖例代表整個(gè)區(qū)域內(nèi)流場(chǎng)的速度。從圖5可以看出:當(dāng)同軸氣流流速為0 m/s時(shí),由于液體重力的作用,噴腔體的液體均流向下方并進(jìn)入噴嘴,此時(shí)噴嘴內(nèi)液體的流速最大,但由于此時(shí)噴腔體內(nèi)部還未充滿進(jìn)液,故噴嘴內(nèi)液體的速度約為0 m/s。隨著進(jìn)液口液體速度的增大,噴嘴內(nèi)的液體速度因壓電片向上運(yùn)動(dòng)而減小,隨后由于進(jìn)液口不斷進(jìn)液以及壓電片向下運(yùn)動(dòng)擠壓噴腔體內(nèi)的液體,進(jìn)液口與噴嘴內(nèi)液體的速度持續(xù)增大,約為0.12 m/s,如圖5(c)所示。由于沒有引入同軸氣流,噴嘴周圍的氣體速度變化極其微小。當(dāng)t=4 500 μs時(shí),位于噴嘴正下方的空氣域向左右兩側(cè)擴(kuò)散最終形成兩個(gè)非常小的渦旋,此時(shí)進(jìn)液口速度最大,約為0.23 m/s,如圖5(d)所示。
(a) t=0 μs
(b) t=1 500 μs
(c) t=3 000 μs
(d) t=4 500 μs
同軸氣流流速v=10 m/s時(shí)一個(gè)噴射周期內(nèi)的流場(chǎng)速度分布如圖6所示。從圖6可以看出:t=0 μs時(shí),進(jìn)氣口開始流入氣體,從進(jìn)氣口流入的氣體首先填充至進(jìn)氣口以及氣流噴射槽,進(jìn)氣口處的氣流流速約為0 m/s,此時(shí)從進(jìn)氣口流入的氣體對(duì)噴嘴下方空氣流動(dòng)的影響可以忽略,隨著氣流不斷地從進(jìn)氣口流入,氣體渦流區(qū)開始形成,在空氣域內(nèi)噴嘴出口兩側(cè)形成一對(duì)渦旋[17],此時(shí)進(jìn)氣口處氣流流速約為13.5 m/s。隨著進(jìn)氣口氣流流速的不斷增大以及從氣流噴射槽噴出的氣體體積不斷增大,渦流區(qū)的尺寸逐漸增大并且向噴嘴下游移動(dòng),t=3 000 μs時(shí)進(jìn)氣口氣流流速約為15.1 m/s,空氣域內(nèi)形成兩對(duì)渦旋,分別位于噴嘴出口兩側(cè)以及噴嘴出口正下方兩側(cè)??諝庥騼?nèi)充分進(jìn)氣后,氣體流速趨于穩(wěn)定,此時(shí)進(jìn)氣口氣流流速約為14.1 m/s,這時(shí)空氣域內(nèi)形成三對(duì)渦旋且逐漸增強(qiáng)擴(kuò)大并繼續(xù)向噴嘴下游移動(dòng),如圖6(d)所示。
(b) t=1 500 μs
(c) t=3 000 μs
(d) t=4 500 μs
噴頭內(nèi)部液體壓力場(chǎng)的變化對(duì)微滴形成過程具有重要的影響。圖7為同軸氣流流速v=0時(shí)一個(gè)噴射周期內(nèi)的流場(chǎng)壓強(qiáng)分布,其中,左側(cè)圖例表示除壓電陶瓷片外其他區(qū)域相對(duì)于標(biāo)準(zhǔn)大氣壓的壓強(qiáng),右側(cè)圖例表示壓電陶瓷片的運(yùn)動(dòng)位移(單位:m),下同。
(a) t=0 μs
(b) t=1 500 μs
(c) t=3 000 μs
(d) t=4 500 μs
由圖7(a)可以看出:受噴射液體自身重力的作用,儲(chǔ)液腔內(nèi)液體與硅膠管交界處B的壓強(qiáng)大于儲(chǔ)液腔內(nèi)液體與壓電陶瓷片交界處A以及玻璃噴嘴內(nèi)液體與空氣域交界處C的壓強(qiáng),此時(shí)最大壓強(qiáng)為61.1 Pa; 壓電陶瓷片受自身重力作用時(shí),壓電片位移約為0 μm。由圖7(b)可以看出:當(dāng)壓電片受壓電信號(hào)驅(qū)動(dòng)而向上彎曲時(shí),噴腔體內(nèi)液體壓強(qiáng)增大而噴嘴內(nèi)液體壓強(qiáng)保持不變,此時(shí)壓電片中心點(diǎn)位移即壓電片振幅為5.54 μm。 當(dāng)壓電片反向彎曲向下運(yùn)動(dòng)時(shí),噴腔體內(nèi)液體壓強(qiáng)減小,玻璃噴嘴內(nèi)部自進(jìn)口端至出口端液體壓強(qiáng)不斷增大,并在液體與空氣域交界的出口端面達(dá)到最大壓強(qiáng),如圖7(c)所示。當(dāng)壓電片自下方最大位移處向上運(yùn)動(dòng)時(shí),液體壓強(qiáng)持續(xù)減小,如圖7(d)所示。
同軸氣流流速v=10 m/s時(shí)一個(gè)噴射周期內(nèi)的流場(chǎng)壓強(qiáng)分布如圖8所示。
(a) t=0 μs
(b) t=1 500 μs
(c) t=3 000 μs
由圖8(a)可以看出,受噴射液體自身重力的作用,儲(chǔ)液腔內(nèi)液體與硅膠管交界處壓強(qiáng)最大,約為58.2 Pa,而壓電陶瓷片受自身重力的影響,其位移約為0 μm。當(dāng)壓電片受到壓電信號(hào)驅(qū)動(dòng)向上彎曲時(shí),噴腔體內(nèi)液體壓強(qiáng)增大而噴嘴內(nèi)液體及空氣域內(nèi)氣體的壓強(qiáng)保持不變,壓電片振幅為6.84 μm,與未施加同軸氣流的情況相比,振幅有所增加,如圖8(b)所示。當(dāng)壓電片受到壓電信號(hào)驅(qū)動(dòng)反向彎曲向下運(yùn)動(dòng)時(shí),噴腔體內(nèi)液體壓強(qiáng)不斷減小,玻璃噴嘴內(nèi)液體與空氣域交界面與進(jìn)液口處壓強(qiáng)達(dá)到最大值,噴頭內(nèi)的壓強(qiáng)變化情況與未施加同軸氣流時(shí)一致,空氣域的壓強(qiáng)因同軸氣流的引入而顯著減小,并且在噴嘴正下方兩側(cè)氣流擴(kuò)散而成的渦流區(qū)壓力值減小更為顯著,如圖8(c)所示。當(dāng)壓電片受壓電信號(hào)驅(qū)動(dòng)逐漸恢復(fù)初始狀態(tài)時(shí),噴腔體內(nèi)液體壓強(qiáng)繼續(xù)減小,噴嘴出口與空氣域交界處以及進(jìn)液入口處的壓強(qiáng)處于較大值,空氣域內(nèi)壓強(qiáng)較未施加同軸氣流時(shí)大幅減小,兩側(cè)氣流擴(kuò)散而成的渦流區(qū)壓強(qiáng)更低,如圖8(d)所示。
圖9為通過COMSOL數(shù)值仿真模擬得到的同軸氣流流速分別為0和5 m/s時(shí)的微滴噴射過程。
由圖9可知,微滴噴射過程可分為以下3個(gè)階段:
(1) 微滴的噴射與延伸階段。初始階段由于壓電片未發(fā)生變形,噴腔體內(nèi)壓強(qiáng)保持初始值,噴腔體內(nèi)部的流體受表面張力和分子間作用力的影響而處于相對(duì)穩(wěn)定狀態(tài)。當(dāng)t=1 500 μs時(shí)壓電片反向彎曲并在噴嘴處形成內(nèi)凹的彎月面。接著壓電片受脈沖信號(hào)作用向下彎曲擠壓噴腔體內(nèi)的液體,壓力波通過噴腔體中的液體傳播至噴嘴處開始形成液面凸起(t= 2 000~2 500 μs),隨后凸起的液滴迅速向下延伸為前端近似半球形的液柱(t=3 000~4 000 μs)。
(2) 微滴斷裂時(shí)刻。從圖9(a)可以看出,未施加同軸氣流時(shí),液滴的斷裂時(shí)間為t=4 500 μs。從圖9(b)可以看出,當(dāng)同軸氣流流速v=5 m/s時(shí),微滴斷裂時(shí)間為t=5 000 μs。由此可知同軸氣流的引入使微滴的斷裂時(shí)刻延后,這與周健等[5]的試驗(yàn)結(jié)果是相吻合的。
(3) 微滴飛行階段。微滴發(fā)生斷裂后因小曲率半徑而被剪斷的尾部末端的壓強(qiáng)較高,故尾部的液體向半球形頭端流動(dòng),出現(xiàn)無(wú)約束長(zhǎng)液帶退縮現(xiàn)象,如圖9(a)中t=4 500~6 500 μs與圖9(b)中t=5 000~ 7 000 μs時(shí)的微滴噴射圖所示。但由于微滴頭部速度仍大于尾部,此時(shí)微滴形成的長(zhǎng)液帶發(fā)生斷尾或多重碎裂,分裂后的微滴尾部逐漸由圓形收縮為小球狀而成為衛(wèi)星液滴。主液滴與衛(wèi)星液滴分離后均失去了噴嘴內(nèi)液體的推力作用并受到一定的空氣阻力,進(jìn)而噴射速度相比噴嘴中新噴射出來(lái)的液柱有所降低。于是從噴嘴擠壓出來(lái)的第二股液柱開始融合第一個(gè)周期結(jié)束后形成的衛(wèi)星液滴與主液滴,如圖9(a)中t= 7 000~8 500 μs與圖9(b)中t= 7 000~9 000 μs 時(shí)的微滴噴射圖所示。融合后進(jìn)一步分裂出主液滴與衛(wèi)星液滴,如圖9(a)中t=8 500~ 9 000 μs以及圖9(b)中t=9 000~9 500 μs時(shí)的微滴噴射圖所示。同軸氣流流速v=0 m/s時(shí)新形成的主液滴由于液滴合并時(shí)其動(dòng)能與表面能的相互轉(zhuǎn)換將以橢球形向下方飛行,如圖9(a)中t=9 000~10 500 μs時(shí)的微滴噴射圖所示,而同軸氣流流速v=5 m/s時(shí)新形成的主液滴則直徑變大,形狀更接近球形,如圖9(b)中t= 9 500~10 500 μs時(shí)的微滴噴射圖所示。
不同速度的同軸氣流作用下微滴即將脫離噴嘴時(shí)的形態(tài)如圖10所示。由圖10可知,在表面張力作用下,液柱前端逐漸縮聚成半球形,成為微滴頭部,此時(shí)液柱上存在兩個(gè)曲率半徑較小的位置,即A、 B處。其中,B處的曲率半徑比A處小,且位于噴嘴出口附近,A處位于圓柱狀液柱與微滴球狀頭部連接處。
(b) v=5 m/s
(c) v=10 m/s
(d) v=15 m/s
圖11為同軸氣流流速對(duì)微滴斷裂前A處微滴直徑D的影響及其隨時(shí)間的變化情況。由圖11可知,隨著同軸氣流流速的增大,A處液體增多,進(jìn)而微滴直徑變大。這是由于在液滴與氣流交匯區(qū)域的左右兩側(cè)形成軸對(duì)稱的若干股渦流,使得噴嘴正下方的壓強(qiáng)減小,因此在同軸氣流的作用下接近A處但位于其上游的表面液體沿液柱向下游加速流動(dòng),并且同軸氣流的流速越快,A處液體增加越多。
圖11 微滴斷裂前直徑隨時(shí)間變化Fig.11 Variation of droplet diameter before fracture over time
在微滴形成過程中,微滴的總長(zhǎng)度L被定義為微滴前端與噴嘴出口端面的垂直距離,其從微滴開始自噴嘴處凸起時(shí)逐漸增加,受到同軸氣流作用時(shí),A處上游的表面液體沿著液柱向下方加速流動(dòng)使得整個(gè)圓柱狀液柱段不斷變長(zhǎng)。微滴斷裂前微滴的長(zhǎng)度隨時(shí)間的變化情況如圖12所示。
圖12 微滴斷裂前長(zhǎng)度隨時(shí)間變化Fig.12 Variation of droplet length before fracture over time
由圖12可知,隨著同軸氣流流速的增加,微滴長(zhǎng)度變長(zhǎng)且微滴斷裂時(shí)刻延后。當(dāng)v=15 m/s時(shí),微滴長(zhǎng)度隨時(shí)間延長(zhǎng)而增加最為明顯;當(dāng)v=5與10 m/s時(shí),微滴長(zhǎng)度在t=3 500 μs之前變化并不明顯。這是由于較小的氣流造成的負(fù)壓對(duì)于微滴表面液體克服其內(nèi)部黏性力從而加速向下流動(dòng)的作用不明顯,因此不會(huì)使微滴長(zhǎng)度產(chǎn)生明顯變化。
利用COMSOL Multiphysics仿真模擬軟件,建立設(shè)有同軸氣流噴射槽的壓電式噴頭的微滴噴射流場(chǎng)模型,并通過數(shù)值模擬研究同軸氣流的引入對(duì)于壓電式噴頭在噴射過程中射流速度場(chǎng)與壓強(qiáng)場(chǎng)以及微滴噴射過程的影響。主要得出以下結(jié)論:
(1) 同軸氣流的引入使得空氣域內(nèi)形成軸對(duì)稱的成對(duì)渦旋,隨著氣流不斷從進(jìn)氣口進(jìn)入,渦旋數(shù)量增多且不斷加強(qiáng)擴(kuò)大并向噴嘴下游移動(dòng)。
(2) 壓電片的形變以及射流壓強(qiáng)場(chǎng)呈周期性變化,同軸氣流的引入使得壓電片振幅變大,空氣域的壓強(qiáng)變小,且位于噴嘴正下方兩側(cè)由同軸氣流擴(kuò)散形成的渦流區(qū)域的壓強(qiáng)變化更明顯。
(3) 隨著同軸氣流流速的增大,微滴斷裂時(shí)刻延后,微滴形狀由橢球形趨于球形,微滴延伸長(zhǎng)度增加,微滴直徑增大。