張波濤，楊寶娥，楊岸龍 ，李 平,2

(1.西安航天動(dòng)力研究所 液體火箭發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710100;2.航天推進(jìn)技術(shù)研究院，陜西 西安 710100)

0 引言

針?biāo)▏娮⑵魇且环N可以連續(xù)調(diào)節(jié)噴注面積的噴注器，具有獨(dú)特的幾何特性和流場(chǎng)特征。針?biāo)▏娮⑵饔泻芏鄡?yōu)點(diǎn)，如結(jié)構(gòu)簡(jiǎn)單、深度調(diào)節(jié)能力、燃燒效率高(96%～99%)和固有的燃燒穩(wěn)定性等。

由于針?biāo)▏娮⑵黛F場(chǎng)濃密，試驗(yàn)采用的光學(xué)設(shè)備很難獲得霧場(chǎng)細(xì)節(jié)特征，因此數(shù)值仿真成為研究針?biāo)▏娮⑵鲊婌F過(guò)程的重要方法。目前關(guān)于針?biāo)▏娮⑵黛F化特性的數(shù)值仿真研究成果很少。王凱等基于Gerris開源軟件采用自適應(yīng)網(wǎng)格加密技術(shù)和分段線性(volume of fluid，VOF)幾何重構(gòu)方法對(duì)液液針?biāo)▏娮卧F化過(guò)程進(jìn)行數(shù)值仿真，分析了液膜撞擊液束的變形過(guò)程，指出液膜撞擊液束噴注單元構(gòu)型優(yōu)于液膜撞擊液膜。與之前采用傳統(tǒng)網(wǎng)格模型和雷諾平均方法求解霧化過(guò)程相比，Gerris數(shù)值方法已有很大的進(jìn)步，可以在減小計(jì)算量的同時(shí)精細(xì)地捕捉到霧化過(guò)程的細(xì)節(jié)特征。文獻(xiàn)[13-16]采用Gerris分別對(duì)內(nèi)混式直流氣液噴嘴、互擊式噴嘴和離心式噴嘴的外流場(chǎng)霧化過(guò)程進(jìn)行了仿真分析，得到了非常好的結(jié)果。Cheng等和Son等分別對(duì)簡(jiǎn)化的二維軸對(duì)稱液液針?biāo)▏娮⑵骱蜌庖横標(biāo)▏娮⑵黛F化角進(jìn)行數(shù)值仿真，給出了霧化角預(yù)測(cè)公式。Radhakrishnan等采用拉格朗日方法中的WAVE模型分析了二維軸對(duì)稱針?biāo)▏娮⑵鞯幕旌咸匦?，并通過(guò)獲得軸向氣體速度和液滴速度來(lái)分析氣體和液滴的相互作用。采用拉格朗日方法的優(yōu)點(diǎn)是通過(guò)對(duì)顆粒追蹤獲得霧場(chǎng)，顆粒對(duì)網(wǎng)格尺寸沒(méi)有要求，可以計(jì)算顆粒從破碎到蒸發(fā)的過(guò)程。由于拉格朗日方法沒(méi)有捕捉氣液兩相界面而無(wú)法細(xì)致分析液相破碎過(guò)程，只適合分析霧場(chǎng)特性。

以前研究者對(duì)液液針?biāo)▏娮⑵骱洼S向縫/徑向縫型氣液針?biāo)▏娮⑵髯隽艘欢ǖ难芯抗ぷ?，但目前鮮有關(guān)于軸向縫/徑向孔型氣液針?biāo)▏娮⑵髌扑檫^(guò)程的研究，鑒于對(duì)液束撞擊氣膜的基礎(chǔ)理論認(rèn)識(shí)不足，本文以軸向縫/徑向孔型的氣液針?biāo)▏娮⑵鳛檠芯繉?duì)象，采用網(wǎng)格自適應(yīng)加密技術(shù)、CLSVOF方法和SBES湍流方法對(duì)氣液針?biāo)▏娮⑵饕菏矒魵饽さ钠扑檫^(guò)程進(jìn)行細(xì)致仿真，并通過(guò)高速攝影試驗(yàn)結(jié)果對(duì)計(jì)算模型進(jìn)行驗(yàn)證。對(duì)氣液針?biāo)▏娮⑵黛F化模式、液束撞擊氣膜變形過(guò)程、流動(dòng)特性和霧化角進(jìn)行系統(tǒng)研究，從而認(rèn)識(shí)液束氣膜相互作用機(jī)制，為氣液針?biāo)▏娮⑵骼碚撗芯亢凸こ淘O(shè)計(jì)提供參考。

1 物理模型和計(jì)算方法

1.1 物理模型

研究對(duì)象為液體中心配置的軸向縫/徑向孔型氣液針?biāo)▏娮⑵?，液體推進(jìn)劑和氣體推進(jìn)劑分別從徑向孔和軸向環(huán)縫噴出，軸向氣膜噴出后到與液束相撞的距離定義為跳躍距離。為了分析單孔液束與氣膜相撞的破碎過(guò)程，在液路噴嘴設(shè)置兩個(gè)對(duì)稱的徑向液束孔。圖1給出了針?biāo)▏娮⑵魇疽鈭D，結(jié)構(gòu)參數(shù)如表1所示。為了充分研究無(wú)量綱參數(shù)局部動(dòng)量比對(duì)液束撞擊氣膜破碎過(guò)程的影響，工況設(shè)計(jì)選取局部動(dòng)量比為0.5～5.5，覆蓋較大的變化范圍，同時(shí)局部動(dòng)量比的變化包含了液束和氣膜速度的變化，具體工況參數(shù)如表2所示。其中局部動(dòng)量比可以定義為徑向液束動(dòng)量與徑向液束孔等寬的軸向氣膜動(dòng)量之比，即

(1)

式中：為徑向矩形孔長(zhǎng)度；為徑向矩形孔寬度；為軸向氣膜厚度。

圖1 氣液針?biāo)▏娮⑵魇疽鈭DFig.1 Schematic diagram of the injector

表1 針?biāo)▏娮⑵鹘Y(jié)構(gòu)參數(shù)

表2 工況條件

1.2 計(jì)算模型

由于物理模型具有對(duì)稱性特征，數(shù)值計(jì)算時(shí)只對(duì)物理模型的1/4區(qū)域進(jìn)行計(jì)算，將對(duì)稱面設(shè)置為對(duì)稱面邊界條件，計(jì)算域如圖2所示。噴嘴結(jié)構(gòu)參數(shù)與物理模型結(jié)構(gòu)參數(shù)一致，外流場(chǎng)計(jì)算域半徑為27 mm，在軸向方向針?biāo)^至下游出口距離為3 mm。采用水作為液體模擬介質(zhì)，空氣作為氣體模擬介質(zhì)。液體和氣體入口均為速度入口邊界，計(jì)算域外邊界為壓力出口邊界，背壓為大氣壓環(huán)境。壓力速度耦合方法采用SIMPLEC(semi-implicit method for pressure linked equation consistent)算法，壓力離散法為PRESTO(pressure staggering option)方法。計(jì)算域網(wǎng)格采用結(jié)構(gòu)網(wǎng)格，初始網(wǎng)格數(shù)量為6×10，使用八叉樹網(wǎng)格自適應(yīng)加密算法對(duì)氣液兩相界面進(jìn)行局部加密，液路噴嘴出口到計(jì)算域下游出口區(qū)域初始網(wǎng)格為80 μm，設(shè)置網(wǎng)格加密/粗化等級(jí)為2級(jí)，加密后最小網(wǎng)格約為20 μm。網(wǎng)格自適應(yīng)函數(shù)設(shè)置為相界面曲率，即網(wǎng)格會(huì)根據(jù)氣液兩相界面曲率值進(jìn)行粗化或加密。

圖2 計(jì)算域Fig.2 Computation domain

1.3 控制方程

由于氣液針?biāo)▏娮⑵鳉饽ち魉傩∮?.3，因此可以看作不可壓縮流體。不考慮液相蒸發(fā)過(guò)程，只需求解連續(xù)方程和動(dòng)量方程。

連續(xù)方程為

(2)

動(dòng)量方程為

(3)

(4)

1.4 CLSVOF方法

界面捕捉方法主要有VOF方法和Level-Set方法，由于VOF方法是基于體積函數(shù)計(jì)算網(wǎng)格內(nèi)流體的體積分?jǐn)?shù)，并非對(duì)氣液兩相界面進(jìn)行捕捉，因此VOF方法具有界面不連續(xù)的缺點(diǎn)。Level-Set 方法是把相界面視為零等值面，通過(guò)求解流場(chǎng)中各個(gè)位置的符號(hào)距離函數(shù)獲得零等值面，缺點(diǎn)是通過(guò)多次迭代會(huì)造成質(zhì)量損失。為了解決VOF方法和Level-Set方法中的缺陷并利用兩種方法的優(yōu)勢(shì)，Sussman提出了將VOF方法和Level-Set方法結(jié)合起來(lái)的方法，即CLSVOF方法。CLSVOF方法兼具VOF方法保證最小質(zhì)量損失和Level-Set方法對(duì)氣液界面連續(xù)捕捉的優(yōu)點(diǎn)，可以計(jì)算連續(xù)液體射流的變形過(guò)程。

1.5 湍流模型

數(shù)值模擬方法有直接數(shù)值模擬、大渦模擬和雷諾平均方法，直接數(shù)值模擬的優(yōu)勢(shì)是不需要湍流模型而直接求解三維定常Navier-Stokes(N-S)方程，從而得到瞬時(shí)湍流流場(chǎng)，但直接數(shù)值模擬對(duì)網(wǎng)格量要求非常高，在實(shí)際工程應(yīng)用中還很少。大渦模擬的特點(diǎn)是用濾波方法把大尺度渦和小尺度渦分開，采用直接數(shù)值模擬求解大尺度渦，通過(guò)亞網(wǎng)格應(yīng)力尺度模型近似模擬小渦對(duì)大渦的影響。雖然大渦模擬相比直接數(shù)值模擬對(duì)網(wǎng)格要求小一些，但計(jì)算量依舊很大。雷諾平均方法計(jì)算量很小，平均化后準(zhǔn)確性較差。但雷諾平均的優(yōu)點(diǎn)是渦黏模型在平衡湍流(湍動(dòng)能生成等于湍動(dòng)能耗散)有很好的適用性，在這種流動(dòng)中沒(méi)必要采用大渦模擬。大渦模擬則適用于非平衡的復(fù)雜湍流。因此在近平衡湍流區(qū)和非平衡湍流區(qū)分別采用雷諾平均和大渦模擬的組合模型是計(jì)算復(fù)雜湍流運(yùn)動(dòng)比較經(jīng)濟(jì)的一種方案。

為細(xì)致捕捉液束與氣膜相互作用產(chǎn)生的湍流結(jié)構(gòu)，采用Menter等提出的SBES方法，在近壁面采用RANS模式，在湍流核心區(qū)采用LES模式，以保證大尺度湍流脈動(dòng)被直接求解。應(yīng)力混合形式為

(5)

(6)

(7)

(8)

上述表達(dá)式各項(xiàng)含義詳見文獻(xiàn)[24-25]。

2 試驗(yàn)方法

2.1 試驗(yàn)系統(tǒng)

試驗(yàn)系統(tǒng)由供應(yīng)系統(tǒng)、測(cè)量系統(tǒng)、試驗(yàn)件和拍攝系統(tǒng)組成，試驗(yàn)中采用過(guò)濾后的水和空氣作為推進(jìn)劑的模擬介質(zhì)，背壓為大氣壓環(huán)境。通過(guò)高壓空氣對(duì)水貯箱進(jìn)行增壓后向噴注器供應(yīng)水，噴注器氣體介質(zhì)直接由高壓氣源供應(yīng)。供應(yīng)管路上設(shè)置科氏流量計(jì)測(cè)量水和氣體的質(zhì)量流量，在噴注器前安裝壓力傳感器測(cè)量噴前壓力。氣液針?biāo)▏娮⑵髟囼?yàn)件結(jié)構(gòu)參數(shù)與物理模型結(jié)構(gòu)參數(shù)一致，為使氣體介質(zhì)在軸向環(huán)縫出口處分布均勻，氣體分兩路進(jìn)入試驗(yàn)件，在試驗(yàn)裝置中匯合后以軸向氣膜形式噴出。液體介質(zhì)從中心管路進(jìn)入試驗(yàn)件，從液路噴嘴設(shè)置的兩個(gè)對(duì)稱徑向液束孔噴出后與軸向氣膜發(fā)生撞擊霧化。使用LED光源照射噴霧場(chǎng)，采用高速攝影相機(jī)拍攝噴霧場(chǎng)，相機(jī)型號(hào)為Phantom V12.1，試驗(yàn)中拍攝頻率為3 kHz，圖像分辨率為640×480像素，曝光時(shí)間為10 μs。

2.2 圖像處理方法

由于液束與氣膜撞擊破碎過(guò)程是瞬態(tài)過(guò)程，為了采用統(tǒng)一的標(biāo)準(zhǔn)分析液束與氣膜撞擊形成的霧化角，對(duì)每個(gè)工況下拍攝的1 000張圖像的平均圖像測(cè)量霧化角，圖像處理過(guò)程如圖3所示。圖3(a)為高速攝影相機(jī)拍攝的原始圖像。對(duì)拍攝的圖像進(jìn)行增強(qiáng)處理，如圖3(b)所示，然后對(duì)1 000張?jiān)鰪?qiáng)處理的圖像求平均圖像，如圖3(c)所示。最后將平均圖像轉(zhuǎn)化為二值圖，如圖3(d)所示。由于液束剛噴出后與氣膜充分相互作用，因此定義霧化角為液束噴出后迎風(fēng)面在徑向2.5 mm處的點(diǎn)和噴嘴出口上沿連線與軸向方向形成的夾角，如圖4所示。為了避免不對(duì)稱性引起的測(cè)量誤差，最終的霧化角取和的平均值。

圖3 圖像處理過(guò)程Fig.3 Image processing

圖4 霧化角Fig.4 Spray angle

2.3 算例驗(yàn)證

選取相同區(qū)域內(nèi)軸向氣膜速度80 m/s時(shí)動(dòng)量比為3.35和5.45的試驗(yàn)結(jié)果驗(yàn)證數(shù)值計(jì)算結(jié)果，如圖5所示?？梢远ㄐ缘乜闯鰯?shù)值仿真獲得的液束宏觀形態(tài)和試驗(yàn)結(jié)果很相似。為了定量分析液束撞擊氣膜形成的霧化角，可以看出數(shù)值計(jì)算霧化角與試驗(yàn)霧化角吻合非常好。

圖5 試驗(yàn)結(jié)果與數(shù)值仿真對(duì)比Fig.5 Comparison between experimental and numerical results

3 結(jié)果與討論

3.1 霧化模式

圖6為氣液針?biāo)▏娮⑵鬏S向氣膜速度為80 m/s時(shí)不同動(dòng)量比下的霧場(chǎng)瞬態(tài)圖像，從圖中可以看出液束與氣膜撞擊后破碎為液滴要經(jīng)過(guò)復(fù)雜的過(guò)程。當(dāng)動(dòng)量比較小時(shí)液束剛噴出即斷裂為液團(tuán)，隨后液團(tuán)破碎為液塊和液滴。隨著動(dòng)量比增大，液束破碎長(zhǎng)度和霧化角均增大，液束迎風(fēng)面在氣動(dòng)力的作用下有小尺度波動(dòng)產(chǎn)生，隨后小尺度波動(dòng)在氣動(dòng)力的作用下沿著液束表面發(fā)展為大尺度波動(dòng)，最終因?yàn)闅鈩?dòng)力大于表面張力和黏性力而導(dǎo)致液束斷裂為液塊，發(fā)生一次霧化，接著大液塊在氣動(dòng)力的作用下發(fā)生二次霧化，破碎為眾多小液滴。同時(shí)，液束迎風(fēng)面在氣動(dòng)力的作用下向液束兩側(cè)運(yùn)動(dòng)，有小液滴從液束兩側(cè)剝離，這一現(xiàn)象始終貫穿于液束發(fā)展的全過(guò)程。

圖6 不同動(dòng)量比下的霧場(chǎng)圖像Fig.6 Spray images with various momentum ratios

當(dāng)動(dòng)量比增大到5以上時(shí)，液束在噴嘴出口近場(chǎng)范圍內(nèi)變形很小，液束迎風(fēng)面在氣動(dòng)力的作用下向兩側(cè)運(yùn)動(dòng)，形成液膜包裹液束的形狀，迎風(fēng)面有表面波形成，氣動(dòng)力促使液束上邊緣和下邊緣相撞，兩側(cè)形成很多液絲，隨后液絲破碎為液滴。

液束撞擊氣膜時(shí)的霧化模式受到局部動(dòng)量比的影響，可以根據(jù)試驗(yàn)結(jié)果分為3種典型的模式：振蕩破碎、剪切破碎和液束主導(dǎo)的波動(dòng)破碎。低動(dòng)量比對(duì)應(yīng)著振蕩破碎，中等動(dòng)量比對(duì)應(yīng)著剪切破碎，高動(dòng)量比對(duì)應(yīng)著液束主導(dǎo)的波動(dòng)破碎。

3.1.1 振蕩破碎

當(dāng)動(dòng)量比較小時(shí)，由于液束動(dòng)量小于氣膜動(dòng)量，當(dāng)氣膜運(yùn)動(dòng)到液束出口時(shí)會(huì)間斷性地堵住液束出口，使得液束振蕩表現(xiàn)為液團(tuán)間歇性脫落特征。圖7給出了動(dòng)量比為0.73時(shí)液束在氣膜中的破碎過(guò)程時(shí)序圖，圖片時(shí)間間隔Δ為2/3 ms，紅色虛線圈內(nèi)代表液束局部變形過(guò)程。當(dāng)氣膜堵住液束出口時(shí)，液束噴注壓降減小，當(dāng)液束噴注通道內(nèi)壓力增大到一定程度后有一股液束噴出，隨后液束噴注通道內(nèi)壓力減小，這時(shí)氣膜又會(huì)把液束出口堵住。這個(gè)過(guò)程表現(xiàn)為液團(tuán)間歇性脫落，脫落后的大塊液團(tuán)在氣動(dòng)力的作用下破碎為小液塊和液滴。當(dāng)運(yùn)動(dòng)到針?biāo)^下方時(shí)會(huì)形成回流區(qū)，氣流向針?biāo)^部下方的軸向中心線靠攏，霧化后的液滴較大且空間分布不均勻，形成收口形霧場(chǎng)。

圖7 振蕩破碎在不同時(shí)刻的霧場(chǎng)Fig.7 Spray field of oscillation breaking at different times

3.1.2 剪切破碎

液束在氣膜作用下的剪切破碎類似于液束在橫向氣流中的剪切破碎過(guò)程，圖8給出了動(dòng)量比為3.35的破碎過(guò)程時(shí)序圖，圖片時(shí)間間隔Δ為2/3 ms，從時(shí)序圖中可以看出剪切破碎的霧場(chǎng)波動(dòng)很小。液束噴出后在氣膜作用下會(huì)發(fā)生彎曲變形并斷裂為液塊，液塊在氣動(dòng)力的剪切作用下破碎為小液滴，同時(shí)液束兩側(cè)受到氣膜剪切作用有液滴剝離，剪切破碎霧化模式動(dòng)量比范圍為1～5，對(duì)應(yīng)的霧場(chǎng)形態(tài)為扇形霧場(chǎng)。

圖8 剪切破碎在不同時(shí)刻的霧場(chǎng)Fig.8 Spray field of shearing breaking at different times

3.1.3 液束主導(dǎo)的波動(dòng)破碎

當(dāng)動(dòng)量比大于5時(shí)會(huì)發(fā)生液束穿透氣膜的現(xiàn)象，圖9給出了動(dòng)量比為5.45的破碎過(guò)程時(shí)序圖，圖片時(shí)間間隔Δ為2/3 ms。液束噴出后在氣膜的作用下變形很小，液束剛噴出后從液束表面剝離的液滴很少。隨著液束的發(fā)展，液束迎風(fēng)面在氣膜的作用下向兩側(cè)移動(dòng)，液束上邊緣形成彎曲的膜狀，與液膜包裹液束的現(xiàn)象十分類似。隨著液束發(fā)展，液束上邊緣的液膜和下邊緣發(fā)生自身相撞導(dǎo)致液束波動(dòng)，在液束迎風(fēng)面處形成明顯的波峰和波谷。液束兩側(cè)的液膜在氣動(dòng)力和波動(dòng)的作用下形成凸起的小液柱和液絲，隨后脫落產(chǎn)生液滴。這種霧化模式以液束自身相撞引起的波動(dòng)破碎為主，氣動(dòng)破碎只起到輔助作用。

圖9 液束主導(dǎo)的波動(dòng)破碎在不同時(shí)刻的霧場(chǎng)Fig.9 Spray field of wave breaking dominated by liquid jet at different times

通過(guò)分析動(dòng)量比對(duì)霧化模式和霧場(chǎng)形態(tài)的影響，認(rèn)為選取中等動(dòng)量比時(shí)的剪切霧化模式及扇形霧場(chǎng)可以在獲得較大霧化角的同時(shí)具有液滴粒徑均勻且液滴空間分布均勻的特性。

3.2 變形過(guò)程

液束在氣膜中發(fā)生彎曲變形的同時(shí)橫截面也會(huì)變形，由于試驗(yàn)中高速攝影拍攝的照片只能獲得液束宏觀變形情況，為了深入分析液束橫截面變形，采用數(shù)值仿真方法獲得液束在不同徑向位置的橫截面。由上述分析可知不同動(dòng)量比下的液束變形過(guò)程不同，為了細(xì)致分析液束橫截面變形過(guò)程，選取低、中、高這3個(gè)動(dòng)量比下的液束橫截面變形過(guò)程進(jìn)行分析。液束噴嘴出口至氣膜變形前的外邊緣距離為4.5 mm，圖10給出了不同動(dòng)量比下液束在噴出后0.2 ～4.5 mm的橫截面圖，圖中Δ為截面距噴嘴出口的距離。當(dāng)動(dòng)量比為1.02時(shí)，液束剛噴入氣膜中即發(fā)生明顯的變形，液束兩側(cè)在氣膜的剪切作用下有液體脫落，隨著液束前緣進(jìn)一步向液束兩側(cè)運(yùn)動(dòng)，液束先發(fā)展為“T”形，隨后被展向拉伸為薄膜狀，當(dāng)液束運(yùn)動(dòng)到4.5 mm時(shí)已完全破碎。隨著動(dòng)量比的增加，液束隨著運(yùn)動(dòng)方向變形減慢。當(dāng)動(dòng)量比為3.35時(shí)，液束噴出噴嘴還沒(méi)有明顯變形，當(dāng)運(yùn)動(dòng)到4.5 mm位置處液束已經(jīng)有大于2/3的液束橫截面發(fā)生變形。當(dāng)動(dòng)量比增大到5.45時(shí)，液束運(yùn)動(dòng)到1.5 mm幾乎沒(méi)有變形，液束前緣向兩側(cè)運(yùn)動(dòng)也明顯減少，運(yùn)動(dòng)到4.5 mm處液束截面類似梯形。

圖10 液束橫截面變形過(guò)程Fig.10 Cross-section deformation process of liquid jet

沿軸向氣膜流動(dòng)方向看液束形態(tài)，液束與氣膜撞擊后形成一個(gè)扇形噴霧區(qū)，其展向?qū)挾入S著垂直高度增加而增大，形成一定的擴(kuò)張角。定義擴(kuò)張角為液束破碎前液束主體兩側(cè)邊緣形成的夾角，如圖11所示。擴(kuò)張角隨著動(dòng)量比增大而減小，這是由于當(dāng)動(dòng)量比較小時(shí)，液束動(dòng)量小使其快速向兩側(cè)擴(kuò)展后破碎為液團(tuán)，霧扇向兩側(cè)拓展很多。隨著動(dòng)量比增加，液束相對(duì)動(dòng)量增加，液束在氣膜作用下不易變形，迎風(fēng)面兩側(cè)剝離的液絲也減少，因此形成的擴(kuò)張角也較小。

圖11 軸向方向的噴霧扇Fig.11 Spray fan structure in axial view

3.3 流動(dòng)特性

為了分析液束與氣膜撞擊時(shí)的流動(dòng)特性，圖12為不同垂直高度處橫截面的體積分?jǐn)?shù)圖和速度矢量圖，在液束根部處氣膜繞矩形液束的流動(dòng)與圓柱繞流有相似之處，都會(huì)在液束后緣兩側(cè)形成一系列對(duì)轉(zhuǎn)的尾跡渦。但也有很多不同之處，Δ表示截面距噴嘴出口的徑向距離，在Δ=0.5 mm處液束前緣只發(fā)生了微小變形，這個(gè)微小變形使得氣流在前緣變形后面形成渦，隨后氣流繞過(guò)液束前緣貼近液束兩側(cè)壁面流動(dòng)，由于速度差對(duì)液束兩側(cè)形成剪切作用。隨著液束迎風(fēng)面展向?qū)挾冗M(jìn)一步增加，在Δ=1.5 mm處液束橫截面近似“T”形，液束前緣明顯的變形引起氣流在液束兩側(cè)形成一系列渦結(jié)構(gòu)，這些復(fù)雜的渦結(jié)構(gòu)增強(qiáng)了氣液兩相的相互作用。

圖13給出了液束與氣膜撞擊時(shí)的液束表面速度云圖和速度矢量圖，動(dòng)量比為3.35，軸向氣膜速度和徑向液束速度分別為80 m/s和4.69 m/s。從圖中可以看出液束迎風(fēng)面速度矢量成輻射狀，液束受氣膜影響速度分布梯度變化較大，液束主體速度接近4.69 m/s，迎風(fēng)面波峰處速度接近30 m/s。同時(shí)液束迎風(fēng)面向兩側(cè)運(yùn)動(dòng)形成膜狀，液膜在氣流的剪切作用下液膜邊緣處速度可達(dá)30 m/s。從液束兩側(cè)剝離的液滴粒徑越小，速度越大，剛剝離的小液滴速度可接近40 m/s，說(shuō)明小液滴在氣流中的跟隨性更好。

圖13 液束表面速度云圖和速度矢量圖Fig.13 Volume fraction and velocity contour of liquid jet surface and velocity vector diagram

經(jīng)以上分析可知液束與氣膜相撞后會(huì)在迎風(fēng)面和兩側(cè)受到氣膜的正應(yīng)力和剪切應(yīng)力，并產(chǎn)生一系列渦結(jié)構(gòu)，這些渦結(jié)構(gòu)增強(qiáng)了氣液兩相之間的相互作用，有助于液束霧化和氣液混合。

3.4 霧化角

圖14 不同動(dòng)量比下的霧化角Fig.14 Spray angle with various momentum ratios

4 結(jié)論

為了全面認(rèn)識(shí)氣液針?biāo)▏娮⑵饕菏矒魵饽さ钠扑檫^(guò)程，本文采用數(shù)值仿真結(jié)合試驗(yàn)的方法對(duì)氣液針?biāo)▏娮⑵饕菏矒魵饽て扑檫^(guò)程進(jìn)行分析研究，得到以下結(jié)論。

1)采用網(wǎng)格自適應(yīng)加密算法、CLSVOF方法并結(jié)合SBES湍流模型的方法，對(duì)氣液針?biāo)▏娮⑵饕菏矒魵饽さ钠扑檫^(guò)程進(jìn)行相界面追蹤，首次給出了液束撞擊氣膜破碎過(guò)程的物理圖畫，獲得了液束變形、斷裂及表面波等破碎特征。同時(shí)與高速攝影相機(jī)結(jié)果定性和定量對(duì)比都可以較好地吻合，故表明該數(shù)值方法在研究氣液噴嘴霧化過(guò)程方面具有較好的準(zhǔn)確性，可以為揭示氣液針?biāo)▏娮⑵髌扑闄C(jī)理提供一種高精度的計(jì)算方法。

2)氣液針?biāo)▏娮⑵鞯撵F化模式和霧場(chǎng)形態(tài)與傳統(tǒng)噴注器有顯著不同。根據(jù)動(dòng)量比，霧化模式可分為振蕩破碎、剪切破碎和液束主導(dǎo)的波動(dòng)破碎，霧場(chǎng)形態(tài)可分為收口形霧場(chǎng)和扇形霧場(chǎng)。

3)液束與氣膜撞擊破碎過(guò)程與液束在橫向氣流中破碎過(guò)程有相似之處，不同點(diǎn)是當(dāng)動(dòng)量比足夠大時(shí)液束會(huì)穿透氣膜。液束與氣膜撞擊后迎風(fēng)面會(huì)向兩側(cè)運(yùn)動(dòng)，橫截面從矩形先發(fā)展為“T”形，隨后被展向拉伸為薄膜。當(dāng)氣動(dòng)力大于表面張力和黏性力時(shí)，液束斷裂。

4)氣膜繞矩形液束流動(dòng)與圓柱繞流有相似之處，都會(huì)在液束后緣兩側(cè)形成一系列對(duì)轉(zhuǎn)的尾跡渦。不同之處在于液束橫截面發(fā)展為“T”形后液束前緣明顯變形引起氣流在液束兩側(cè)形成一系列渦結(jié)構(gòu)，增強(qiáng)了氣液相互作用。

6)經(jīng)綜合分析在寬范圍動(dòng)量比下的霧化角、霧化模式和霧場(chǎng)形態(tài)，認(rèn)為選取中等動(dòng)量比下的剪切霧化模式可以在有較大霧化角的同時(shí)液滴粒徑均勻且液滴空間分布均勻。