王　婕,余永剛,劉　焜,2

(1.南京理工大學 能源與動力工程學院,南京 210094;2.上海齊耀動力有限公司,上海 201213)

?

對撞式噴嘴在大氣環(huán)境中噴霧實驗與數(shù)值模擬

王婕1,余永剛1,劉焜1,2

(1.南京理工大學 能源與動力工程學院,南京 210094;2.上海齊耀動力有限公司,上海 201213)

摘要:為了研究對撞式噴嘴在大氣環(huán)境中的噴霧特性,設計了對撞式噴嘴噴霧實驗裝置。采用相位多普勒粒子動態(tài)分析儀(PDA),測量了不同工況下含能液體模擬工質(zhì)霧化場液滴的特征參數(shù),獲得了不同噴射壓力下液滴索特爾平均直徑D32的分布規(guī)律以及噴霧場液滴軸向速度vz沿軸向和徑向的分布規(guī)律。在實驗基礎上,建立噴霧場二維非穩(wěn)態(tài)簡化模型,氣相采用k-ε湍流模型,液滴破碎采用泰勒類比破碎模型。結(jié)果表明:噴霧場從起始到統(tǒng)計穩(wěn)定的發(fā)展過程約需15 ms,液滴索特爾平均直徑D32的分布規(guī)律、軸向速度的計算值與實驗結(jié)果基本吻合。

關鍵詞:液體火箭發(fā)動機;對撞式噴嘴;霧化特性;數(shù)值模擬

液體火箭發(fā)動機由于其性能和工作可靠性等方面的優(yōu)勢,在航天領域具有廣泛的用途[1-3]。液體火箭發(fā)動機有一個關鍵部件噴嘴直接影響發(fā)動機的工作性能,而對撞式噴嘴通過兩股射流的對撞霧化,具有響應快、能迅速混合和燃燒、易于制造的優(yōu)點,是液體火箭發(fā)動機常用的燃料噴射器。

針對對撞式噴嘴的霧化特性,國內(nèi)外學者展開了一系列的研究。張蒙正等[4-5]用激光全息及圖像處理技術,研究兩股對撞射流的撞擊狀況及霧化性能,獲得了噴嘴的霧化性能參數(shù)。Anderson等人[6]對兩股射流的撞擊霧化特性進行了實驗研究和數(shù)值模擬。在不同的流動狀態(tài)和噴嘴結(jié)構(gòu)下,實驗測量了液滴粒徑、液膜破裂的長度和寬度、液膜表面波的結(jié)構(gòu),并利用線性穩(wěn)定性理論對液膜破碎進行建模,對破碎長度和液滴粒徑進行了模擬計算。孫紀國等[7]用激光散射粒子儀和探針對一種40°撞擊角的對撞式噴嘴做了實驗研究。結(jié)果表明,在一定的噴嘴壓降下,噴霧形狀為橢圓形,噴嘴壓降越大,離撞擊點越遠,霧化質(zhì)量越好。Bremond等[8]針對射流沖擊霧化,研究了從2個相同的圓柱形射流斜碰撞形成液膜及其破碎的過程特征,分析了液膜形狀、邊緣尺寸、流體速度場、液滴的形狀和尺寸大小等參數(shù)。Choo等[9]針對射流速度分布對2個沖擊射流形成液膜的速度和厚度做了研究,對于3種不同的射流速度,通過求解質(zhì)量、動量和能量守恒方程對液膜的厚度和速度分布進行預測,并與實驗結(jié)果進行對比。Syed等[10]實驗研究了對撞噴嘴霧化凝膠推進劑,重點分析了噴嘴幾何形狀對霧化效果的影響。Baek等[11]研究了含有納米粒子凝膠材料對撞噴嘴的霧化特性,重點研究了不同射流速度對工質(zhì)霧化特性的影響。Avulapati等[12]實驗研究了采用空氣輔助的方法對于對撞射流霧化特性的影響。把氣體射流引導到兩股液體射流撞擊點處以提升霧化效果,重點分析了液體粘度、表面張力、液體射流間的對撞角度、氣體噴射射流直徑對噴霧特性的影響。鄭剛等[13]研究了撞擊角度對于對撞式噴嘴霧化特性的影響,采用CLSVOF方法對對撞式噴嘴霧化過程進行數(shù)值模擬。結(jié)果表明,撞擊角增加,液膜破碎長度減小,液膜破碎程度和表面波幅值增加。

本文針對某含能液體模擬工質(zhì),采用相位多普勒粒子動態(tài)分析儀(PDA),采集了對撞式噴嘴在大氣環(huán)境的霧化特性參數(shù),側(cè)重分析壓力變化對噴霧場霧滴軸向速度和平均粒徑等參數(shù)的影響,并在實驗的基礎上建立簡化的二維霧化場模型,對霧化過程進行了數(shù)值模擬。

1噴霧實驗裝置

為了研究對撞式噴嘴在大氣環(huán)境中的噴霧特性,設計了噴霧實驗裝置,如圖1所示,由壓氣機、儲液箱、流量計、壓力表、PDA、噴嘴等組成。實驗開始先將模擬工質(zhì)注入儲液箱,調(diào)節(jié)壓氣機到指定實驗壓力,再對PDA試驗系統(tǒng)進行手動調(diào)焦。待噴霧穩(wěn)定后,用PDA系統(tǒng)進行測量,通過計算機輸出數(shù)據(jù)。

實驗采用的對撞式噴嘴結(jié)構(gòu)如圖2所示,兩股射流間撞擊角2θ=30°,噴孔直徑d0=0.23 mm,中心距l(xiāng)0=3 mm。本次實驗采用PDA實驗系統(tǒng),該系統(tǒng)可以同時測量球形顆粒尺寸和速度,對霧化場無接觸、無干擾,粒徑測量范圍為0.5～200 μm,測量精度在1%以內(nèi);速度測量范圍為-300～1 000 m/s,測量精度在0.1%以內(nèi)。

圖2　對撞式噴嘴示意圖

坐標軸確定方法如圖3所示,兩噴嘴連線中點為坐標原點O,噴嘴出口方向為z軸,與z軸垂直的平面是Oxy平面,在該平面內(nèi)與z軸相交且平行于噴嘴連線方向為x軸方向,垂直于噴嘴連線方向為y方向。r為測量點在Oxy平面投影到原點O的距離。

圖3　坐標系示意圖

2實驗結(jié)果與分析

霧化工質(zhì)采用含能液體推進劑模擬工質(zhì),粘度為1×10-3Pa·s,噴射壓力分別取2.2 MPa和2.6 MPa。噴射壓力為2.2 MPa與2.6 MPa下各截面液滴平均直徑D32的變化如圖4所示。索特爾平均直徑體現(xiàn)了液滴表面積和體積的相對大小,其值大小反映了液滴與周圍介質(zhì)進行熱質(zhì)交換的強弱。當p=2.2 MPa時,液滴平均直徑D32沿軸向先減小后逐步增加。當p=2.6 MPa時,液滴平均直徑D32逐步增加,當軸向距離大于49.5 mm后,液滴平均直徑D32增加速率減小。這表明,隨著軸向距離的增加,液滴受到沿程阻力,速度減小,液滴發(fā)生聚合與二次碰撞的概率增加,使液滴的粒徑增加。p為2.6 MPa情況下的D32均小于p為2.2 MPa情況下的D32,壓力增加使液滴間相互碰撞作用力增大,液滴粒徑減小,這表明增大噴射壓力可以減小D32。

圖4　霧化液滴平均直徑D32沿軸向的變化規(guī)律

噴射壓力分別為2.2 MPa和2.6 MPa,不同液滴索特爾平均直徑D32的數(shù)目分布P如圖5所示。可以看出,2種壓力下大部分液滴直徑集中在52～68 μm范圍內(nèi),隨著壓力的增大,液滴索特爾平均直徑D32在44 μm以下數(shù)目減少。在p=2.2 MPa時,液滴D32在68～100 μm之間的數(shù)目與在52～68 μm間的數(shù)目相當。p=2.6 MPa時,索特爾平均直徑大于68 μm的數(shù)目是液滴索特爾平均直徑在52～68 μm之間的60%。隨著壓力的增大,在44～68 μm間的液滴數(shù)目增加,2.6 MPa情況比2.2 MPa情況下液滴D32數(shù)目整體增加10.7%,這說明隨著噴射壓力的增加,射流動量也增大,大氣環(huán)境中空氣擾動加劇,液滴在下游發(fā)生二次碰撞的概率增加,液滴索特爾平均直徑往減小的方向發(fā)展,霧化液滴尺寸分布均勻。

圖6給出了噴射壓力為2.6 MPa下霧化場液滴軸向速度沿軸向和徑向的分布特性。液滴軸向速度u沿軸向先增加后減小,在79.7 mm和102 mm截面處,隨著徑向距離的增加,液滴軸向速度整體變化趨勢相同,都是逐漸減小,這是由于液滴往下游發(fā)展,有2股射流碰撞產(chǎn)生能量損失以及沿程阻力的存在,導致液滴速度逐步衰減。在r<17.5 mm之處,79.7 mm橫截面上的液滴軸向速度大于102 mm橫截面處的液滴軸向速度;r>17.5 mm時,102 mm橫截面處的液滴軸向速度大于79.7 mm橫截面處液滴的軸向速度。

圖5　大氣中各測點的液滴平均直徑D32數(shù)目分布

3簡化模型及數(shù)值模擬

3.1物理模型

針對對撞式噴嘴的特點,對噴霧過程做出如下簡化假設:

①將對撞射流噴霧場近似作為二維氣液兩相流場處理;

②噴霧過程中液滴為離散相,不考慮重力的作用;

③氣相采用k-ε湍流模型;

④液滴破碎采用泰勒類比破碎模型,液滴曳力選球形,計算過程中不考慮液滴變形。

3.2數(shù)學模型

3.2.1氣相控制方程

質(zhì)量守恒方程為

(1)

動量守恒方程為

(2)

湍流模型采用k-ε模型,k方程為

(3)

ε方程為

(4)

式中:μt為湍流粘性系數(shù),Gk是由平均速度梯度引起的湍動能k的生成項;Gb是由浮力產(chǎn)生的湍動能k的生成項;YM為在可壓縮湍流中脈動擴張的貢獻;C1ε,C2ε,C3ε為經(jīng)驗常數(shù);σk和σε分別是與湍動能k和耗散率ε對應的Prandtl數(shù)。

3.2.2離散相控制方程

液滴的作用力平衡方程在笛卡爾坐標系下的形式(i方向)為

(5)

式中:Fd(v-vd)為液滴的單位質(zhì)量曳力;v為氣相速度;vd為液滴速度;Fi為i方向上的其他作用力。

3.3幾何模型與網(wǎng)格劃分

利用Gambit對對撞式噴嘴在大氣環(huán)境中的噴射霧化進行二維建模,噴嘴尺寸如圖2所示,大氣環(huán)境中幾何模型如圖7所示,計算區(qū)域為200 mm×150 mm,總網(wǎng)格數(shù)約120 000個,在噴嘴附近處氣液相互作用比較激烈,對這一區(qū)域網(wǎng)格做了加密處理。

圖7　大氣環(huán)境中幾何模型

3.4數(shù)值模擬結(jié)果與分析

對噴射壓力為2.6 MPa時的對撞式噴嘴在大氣環(huán)境中的霧化特性進行數(shù)值模擬,霧化場液滴粒徑分布隨時間變化如圖8所示。

圖8　霧化場液滴粒徑分布序列圖

由圖可見,兩股射流在撞擊點相撞后,繼續(xù)沿軸向和徑向發(fā)展,形成以撞擊點為頂點的錐形霧化場,錐形外部有飛濺的細小液滴,當t>15 ms時霧化場參數(shù)達到統(tǒng)計穩(wěn)態(tài)。

噴嘴噴射壓力為2.6 MPa時對撞式噴嘴霧化場液滴沿軸向截面的液滴索特爾平均直徑D32的模擬結(jié)果與實驗結(jié)果如圖9所示。由圖可見,液滴D32數(shù)值模擬結(jié)果趨勢與實驗結(jié)果相同,平均誤差為6.53%。

圖10給出了噴嘴噴射壓力為2.6 MPa時對撞式噴嘴霧化場液滴u沿軸向和在102 mm截面處沿徑向的數(shù)值模擬結(jié)果與實驗結(jié)果。液滴沿軸向的軸向速度的數(shù)值模擬結(jié)果比實驗結(jié)果要小,平均偏小5.9%。在102 mm截面處,液滴軸向速度數(shù)值模擬結(jié)果發(fā)展趨勢與實驗結(jié)果基本一致,液滴軸向速度沿徑向逐漸減小,是射流能量損失造成的。數(shù)值模擬與實驗結(jié)果有偏差,一部分是因為噴霧過程是一個三維過程,把它簡化為二維過程帶來的偏差;另一部分是因液滴碰撞破碎與聚合模型以及初始條件設定產(chǎn)生的誤差。

圖9　液滴D32數(shù)值模擬結(jié)果與實驗結(jié)果軸向分布比較

圖10　液滴軸向速度數(shù)值模擬結(jié)果與實驗結(jié)果沿軸向和徑向分布對比

4結(jié)論

綜合本文實驗與數(shù)值模擬結(jié)果,得出如下結(jié)論:

①隨著噴射壓力的增加,液滴間相互碰撞作用力增大,液滴索特爾平均直徑D32減小,噴射壓力在2.6 MPa時,液滴索特爾平均直徑D32沿軸向的增加速度小于噴射壓力為2.2 MPa時液滴索特爾平均直徑D32的增加速度。

②在一定的截面上,測量點離中心軸越遠,液滴軸向速度越小,由于液滴往下游發(fā)展,有兩股射流碰撞產(chǎn)生能量損失以及沿程阻力的存在,導致液滴速度逐步衰減。

③數(shù)值模擬結(jié)果與實驗結(jié)果吻合較好,平均誤差在5.9%～6.5%之間,說明本文采用的理論模型和數(shù)值模擬方法是可行的。

參考文獻

[1]CASIANO M J,HULKA J R,YANG V.Liquid-propellant rocket engine throttling:a comprehensive review[J].Journal of Propulsion and Power,2010,26(5):897-923.

[2]SUTTON G P.History of liquid propellant rocket engines in the united states[J].Journal of Propulsion and Power,2003,19(6):978-1 007.

[3]沈赤兵,王克昌,陳啟智.國外小推力液體火箭發(fā)動機的最新進展[J].上海航天,1996,3:41-45.

SHEN Chi-bing,WANG Ke-chang,CHEN Qi-zhi.The latest progress in small thrust liquid rocket engines abroad[J].Aerospace Shanghai,1996,3:41-45.(in Chinese)

[4]張蒙正,傅永貴,張澤平,等.兩股互擊式噴嘴霧化研究及應用[J].推進技術,1999,20(2):74-77.

ZHANG Meng-zheng,FU Yong-gui,ZHANG Ze-ping,et al.Spray property research and application of unlike impinging injector[J].Journal of Propulsion Technology,1999,20(2):74-77.(in Chinese)

[5]張蒙正,張澤平,李鰲,等.兩股互擊式噴嘴霧化性能實驗研究[J].推進技術,2000,21(1):58-60.

ZHANG Meng-zheng,ZHANG Ze-ping,LI Ao,et al.Experimental research on spray properties of unlike impinging injectors[J].Journal of Propulsion Technology,2000,21(1):58-60.(in Chinese)

[6]ANDERSON W E,RYAN H M,PAL S,et al.Fundamental studies of impinging liquid jets,AIAA-1992-0458[R].1992.

[7]孫紀國,王玨,沈赤兵,等.一種40°撞擊角雙股自擊式噴嘴試驗[J].推進技術,2002,23(3):207-208,222.

SUN Ji-guo,WANG Jue,SHEN Chi-bing,et al.Experimental for a like-doublet injector with 40° impingement angle[J].Journal of Propulsion Technology,2002,23(3):207-208,222.(in Chinese)

[8]BREMOND N,VILLERMAUX E.Atomization by jet impact[J].Journal of Fluid Mechanics,2006,549:273-306.

[9]CHOO Y J,KANG B S.The effect of jet velocity profile on the characteristics of thickness and velocity of the liquid sheet formed by two impinging jets[J].Physics of Fluids,2007,19(11):1-7.

[10]FAKHRI S,LEE G J,YETTER R.Atomization and spray characteristics of gelled-propellant simulants formed by two impinging jets,AIAA 2009-5241[R].2009.

[11]BAEK G,KIM S,HAN J.Atomization characteristics of impinging jets of gel material containing nanoparticles[J].Journal of Non-Newtonian Fluid Mechanics,2011,166(21-22):1 272-1 285.

[12]AVULAPATI M M,VENKATA R R.Experimental studies on air-assisted impinging jet atomization[J].International Journal of Multiphase Flow,2013,57:88-101.

[13]鄭剛,聶萬勝,何博,等.撞擊角對撞擊式噴嘴霧化特性影響研究[J].推進技術,2015,36(4):608-613.

ZHENG Gang,NIE Wan-sheng,HE Bo,et al.Effects of impingement angle on atomization characteristics of impinging jets injector[J].Journal of Propulsion Technology,2015,36(4):608-613.(in Chinese)

Experiment and Simulation Studies on Impinging Nozzle in Atmospheric Environment

WANG Jie1,YU Yong-gang1,LIU Kun1,2

(1.School of Energy and Power Engineering,Nanjing University of Science and Technology,Nanjing 210094,China;2.Shanghai MicroPowers Co.Ltd,Shanghai 201213,China)

Abstract:To study the spray characteristics of impinging nozzle in atmospheric environment,the experiment device of impinging nozzle was designed.The distribution characteristics of the spray parameters in the spray field under different conditions were measured with PDA system.Sauter mean diameter(D32)and axial velocity(vz)along the axial and radial direction were gained under different injection pressures.Based on the experiment results,the two-dimensional unsteady simplified model was established.The gas phase was simulated by usingk-εturbulence model,and Taylor analogy breakup model was used in simulating droplets’breakup.The results show that the time of the spray field being statistically stable takes about 15 ms.The simulated results ofD32andvzagree well with the experimental data.

Key words:liquid-propellant rocket engine;impinging nozzle;spray characteristic;numerical simulation

收稿日期:2015-12-02

基金項目:教育部博士點基金項目(20113219110024)

作者簡介:王婕(1992- ),女,碩士研究生,研究方向為含能材料的霧化。E-mail:13770831835@163.com。 通訊作者:余永剛(1963- ),男,教授,研究方向為含能材料的霧化與燃燒。E-mail:yonggangy@21cn.com。

中圖分類號:V434

文獻標識碼:A

文章編號:1004-499X(2016)02-0047-06