孫慧娟，張海濱，白博峰

(西安交通大學(xué)動(dòng)力工程多相流國家重點(diǎn)實(shí)驗(yàn)室，西安 710049)

0 引言

水霧通常用于消防、加濕、冷卻、食品加工行業(yè)等，近年來，水霧應(yīng)用到水下動(dòng)力系統(tǒng)中，為動(dòng)力裝置提供氧化劑、做功工質(zhì)和冷卻劑。以水沖壓發(fā)動(dòng)機(jī)為例，從燃燒室排出的混合燃?xì)膺M(jìn)入摻混室，水霧沿著垂直于摻混室圓管軸線方向進(jìn)入，與混合燃?xì)鈾M向交叉混合，完成水霧蒸發(fā)和兩相摻混過程，實(shí)現(xiàn)混合氣體降溫和增加噴管做功工質(zhì)的作用。較高的溫降效率和較好的摻混效果，對于提高發(fā)動(dòng)機(jī)的推力至關(guān)重要。國內(nèi)外對水沖壓發(fā)動(dòng)機(jī)的研究多側(cè)重于金屬燃料［2－8］、理論性能［9－13］、結(jié)構(gòu)組織方式［14－16］、熱力計(jì)算［17－19］及流場預(yù)示［20－22］等方面，缺乏對摻混機(jī)制的探索。因此，有必要研究不同因素影響下的兩相摻混機(jī)制，為發(fā)動(dòng)機(jī)的設(shè)計(jì)提供理論依據(jù)。

在噴霧射流與橫流氣體的摻混中，隨著射流對橫流的卷吸，在橫流中擴(kuò)散的液滴被夾帶，氣相渦結(jié)構(gòu)不斷長大，在有限空間內(nèi)發(fā)展，渦結(jié)構(gòu)對兩相摻混的影響很大。Salewski和Fuchs［23］采用大渦模擬證明反旋渦對(CVP)是一個(gè)主要的渦結(jié)構(gòu)，液滴的擴(kuò)散依賴CVP的強(qiáng)度。Bai等［24］通過 PIV觀察到在受限空間內(nèi)，CVP是影響摻混的主要渦結(jié)構(gòu)，且形態(tài)受入射角度影響較大。Choi等［25］用米氏散射技術(shù)測量橫流中的燃料霧化結(jié)構(gòu)，發(fā)現(xiàn)在噴霧底部形成兩個(gè)獨(dú)立的旋渦，認(rèn)為旋渦是由燃料與噴嘴壁面間的摩擦產(chǎn)生的剪切力造成的，旋渦卷吸橫流氣體到噴霧中，有利于強(qiáng)化摻混。Pan?o和 Moreira［26］揭示了橫流條件下，壁面附近的 3維渦是影響噴霧和壁面作用的主要因素。除了渦結(jié)構(gòu)，噴霧在橫流中的研究還涉及液滴速度分布、直徑分布、破碎、蒸發(fā)及控制摻混過程機(jī)理的探索［27－30］。通過大量研究，學(xué)者們普遍認(rèn)識(shí)到霧化壓力［31］、液體流量［32］、液滴/橫流速度比［28］和密度比［33］等對近噴霧區(qū)和遠(yuǎn)場的液滴分布和噴霧結(jié)構(gòu)均有不同程度的影響，噴霧和氣流的相互作用同樣也會(huì)影響液滴蒸發(fā)、液滴破碎和噴霧軌跡。

前期工作［34－35］從冷態(tài)摻混角度分析了不同水氣比、噴嘴霧化角度、液滴直徑、噴嘴入射角度對液滴群和橫流氣體的摻混過程影響。對冷態(tài)而言，液滴分布是衡量摻混的因素;而對熱態(tài)來說，溫度分布則是一個(gè)重要指標(biāo)。對渦結(jié)構(gòu)的影響因素很多，本文研究噴霧霧化錐角對蒸發(fā)的噴霧液滴在橫流中摻混的影響，從摻混度和溫降效果兩方面衡量霧化錐角對有相變的兩相摻混的影響程度，分析溫度分布規(guī)律和流場規(guī)律，獲得對摻混過程的認(rèn)知，為水沖壓發(fā)動(dòng)機(jī)概念原型的優(yōu)化設(shè)計(jì)和參數(shù)組織方式提供理論依據(jù)。

1 數(shù)學(xué)模型

1.1 連續(xù)相模型

對于噴霧橫流兩相流動(dòng)，采用離散顆粒群軌跡模型，連續(xù)相的數(shù)學(xué)描述采用歐拉方法，離散相采用拉格朗日方法描述，在商業(yè)軟件FLUENT的基礎(chǔ)上計(jì)算。根據(jù)質(zhì)量、動(dòng)量、能量和組分守恒，連續(xù)相采用Realizable k－ε 湍流模型，統(tǒng)一為

流場的數(shù)值求解采用SIMPLE算法，控制方程的離散格式采用二階迎風(fēng)格式。

1.2 離散相模型

不考慮噴霧液膜的破碎、液滴的形成過程細(xì)節(jié)，只關(guān)注噴霧初始形態(tài)和產(chǎn)生液滴的統(tǒng)計(jì)直徑。采用壓力－旋流霧化模型，即Schmidt等［36］提出的線性非穩(wěn)定液膜霧化模型，把液滴加入橫流氣體。采用TAB模型［37］考慮噴霧液滴的破碎，采用O＇Rourke統(tǒng)計(jì)顆粒模型［38］考慮顆粒的碰撞和聚并，主要考慮拖曳力和重力，液滴受力平衡方程為

其中:

式中 ug和up分別為氣相和液滴速度;μg為氣體動(dòng)力粘度;CD為拖曳力系數(shù);Rep為顆粒雷諾數(shù);τp為顆粒速度弛豫時(shí)間;Dp和ρp分別為液滴直徑和密度。

湍流的影響通過顆粒隨機(jī)步進(jìn)模型考慮，計(jì)算液滴運(yùn)動(dòng)時(shí)，氣相瞬時(shí)速度由平均速度和脈動(dòng)速度構(gòu)成:

兩相間動(dòng)量交換為

假設(shè)傳熱和傳質(zhì)速率相等，液滴蒸發(fā)速率為

通過對流傳熱和潛熱間的平衡得到液滴溫度:

兩相間質(zhì)量交換為

能量交換為

式中 λg為氣體熱導(dǎo)率;cp，g、cp，p分別為氣體和液滴的比定壓熱容;BM為傳質(zhì)系數(shù);L為液體的潛熱。

動(dòng)量、質(zhì)量和能量交換作為氣相方程求解的源項(xiàng)來實(shí)現(xiàn)雙向耦合。

液滴和壁面作用采用反彈模型［39］。反彈后的法向van和切向速度vat為

其中，e為恢復(fù)系數(shù)，表示成液滴入射方向和壁面夾角θi的函數(shù):

2 幾何結(jié)構(gòu)及模擬條件

計(jì)算采用的幾何模型為長500 mm，內(nèi)徑95 mm的圓柱腔，4個(gè)旋流霧化噴嘴垂直圓柱且沿周向均布在距離圓管入口100 mm處，見圖1。水由噴嘴霧化成索太爾直徑一定的小液滴，以不同的霧化錐角進(jìn)入圓柱腔，在橫流氣體作用下進(jìn)行蒸發(fā)和摻混。

圖1 幾何結(jié)構(gòu)Fig.1 Geometrical structure

氣相(空氣)入口速度30 m/s，溫度600 K，摻混室壓強(qiáng)101 325 Pa。單噴嘴流量 0.02 kg/s，液滴 D32為150～160 μm，初始溫度300 K，噴嘴垂直入射(圖1中α =90°，β =0°)，霧化錐角(圖1 中 φ)為40°、60°、80°、90°和 100°。

3 結(jié)果與分析

3.1 實(shí)驗(yàn)驗(yàn)證

將數(shù)值模擬結(jié)果與自行設(shè)計(jì)實(shí)驗(yàn)的結(jié)果進(jìn)行對比。實(shí)驗(yàn)段為截面95 mm×95 mm的矩形腔，入口空氣溫度673.15 K，速度30 m/s。實(shí)驗(yàn)系統(tǒng)可參見文獻(xiàn)［24］。單噴嘴丙酮流量0.001 kg/s，初始溫度284 K，初始液滴D32=48 μm，霧化錐角80°。圖2比較了噴嘴后不同橫截面豎直中心線的溫度降(Dc－n表示當(dāng)前截面與噴嘴所在截面的距離)，圖2顯示模擬和實(shí)驗(yàn)結(jié)果吻合較好。

圖2 實(shí)驗(yàn)與模擬的橫截面豎直中心線的溫度降比較Fig.2 Comparsions of temperature drops at the center lines of the cross sections between experimental and simulation results

3.2 不同霧化錐角的摻混結(jié)果

根據(jù)冷態(tài)氣液摻混的研究結(jié)果［34］發(fā)現(xiàn)，霧化角度120°時(shí)，液滴分布明顯不如90°均勻，隨著霧化角度的逐漸增大，液滴分布更不均勻?；谠摽紤]，霧化角度在大于90°時(shí)只取100°。研究了5種霧化錐角40°、60°、80°、90°和 100°下的噴霧在橫流中的蒸發(fā)摻混過程。以摻混室出口的溫降水平和摻混效果作為判別霧化錐角對摻混過程利弊的準(zhǔn)則。溫降效果為

式中 Tin為摻混室入口橫流氣體溫度;Tave，out為摻混室出口混合氣體平均溫度。

摻混度體現(xiàn)了摻混室橫截面溫度或組分分布的均勻程度。本文根據(jù)灰關(guān)聯(lián)分析的思想，提出了摻混度的定義。該方法是分析灰色系統(tǒng)中各因素間關(guān)聯(lián)程度的一種量化方法，根據(jù)序列曲線幾何形狀的相似程度來判斷灰色過程發(fā)展態(tài)勢的關(guān)聯(lián)程度。

式中 X為所關(guān)注的參量，如溫度、組分濃度等;Xi為截面上區(qū)域i對應(yīng)的參量X;ˉX表示截面參量平均值;ωi為權(quán)重因子，表示不同區(qū)域的參量對于整個(gè)截面的影響程度，本文取截面各處權(quán)重因子相等。

摻混度在0～1之間，值越大說明摻混效果越好。本文采用溫度作為關(guān)注參量。

圖3給出了5種霧化錐角下的溫降效果和摻混度。隨霧化錐角增大，溫降效果提高，但霧化錐角大于90°后，溫降效果增加速率減緩。摻混度先增大，霧化錐角大于90°后，摻混度略有降低。

圖3 不同霧化錐角影響下的摻混和溫降效果Fig.3 Mixing degree and cooling effect under different spray angles

3.3 溫度分布

圖4給出了5種霧化錐角下?lián)交旖孛娴臏囟鹊母怕屎瘮?shù)分布(PDF)。

式中 ΔT為溫度小量;Δs為(T－1/2ΔT，T+1/2ΔT)的溫度范圍對應(yīng)的橫截面的小區(qū)域之和;S為覆蓋截面所有溫度范圍的Δs總和。

圖4 不同霧化角下不同摻混截面溫度的PDF分布Fig.4 Temperature distributions(PDF)at different cross sections under different spray angles

某溫度對應(yīng)的PDF值越高，說明橫截面上該溫度分布范圍大。

圖4表明，不同霧化錐角下，隨摻混距離增加，截面溫度分布越趨向更低的溫度范圍。初始截面，溫度分布在橫流初溫600 K附近，摻混距離越長，截面上低溫區(qū)域所占的比例越大。霧化錐角40°和60°時(shí)，雖然隨著摻混距離變長，摻混截面溫度降低，但低溫區(qū)對應(yīng)的PDF值偏低，即低溫區(qū)域所占比例較小，溫降效果較差。霧化錐角80°、90°和100°時(shí)，初始摻混截面的平均溫度就略低于較小的霧化錐角對應(yīng)的初始截面平均溫度。隨霧化錐角增加，較低溫度區(qū)域的溫度跨度增大，低溫度的PDF值有增加趨勢。霧化錐角100°時(shí)，出口截面平均溫度最低。所以，該霧化錐角帶來的溫降效果最好。但由于初始噴霧的中空形態(tài)，霧化錐角越大，液滴越趨向近壁面區(qū)域，液滴分布的不均勻性，影響了摻混截面的溫度分布。霧化錐角小于90°時(shí)，摻混室中心向壁面處溫度逐漸升高，而霧化錐角100°時(shí)，摻混室近壁面處溫度低于管中心溫度，霧化錐角的增大，造成了溫度變化趨勢的反轉(zhuǎn)。

3.4 摻混度變化規(guī)律

圖5為摻混度隨摻混距離的變化規(guī)律。摻混度隨摻混距離的增加，由初始摻混截面的較大摻混度降低，或降低后再略有升高。

圖5 不同霧化角錐下的摻混度隨摻混距離的變化Fig.5 Mixing degrees under different spray angles

從圖5可見，初始摻混截面溫度較高，由于蒸發(fā)距離很短，液滴蒸發(fā)帶來的降溫量小，溫度基本在600 K附近。所以，整體截面溫度偏差較小，摻混度反而較高。隨著摻混距離增大，液滴蒸發(fā)量增加，截面溫度逐漸降低，從摻混室中心到近壁面區(qū)域的溫度不均勻分布，造成了摻混度的降低。要獲得較高摻混度，要求溫度范圍區(qū)域小，且低溫對應(yīng)的PDF峰值大，該情況下?lián)交燧^均勻。霧化錐角40°和60°時(shí)，在截面溫度覆蓋范圍內(nèi)，PDF值都較小，且相差不大，即不同溫度的分布區(qū)域相當(dāng)，沒有較大范圍的溫度相近區(qū)，造成了摻混度較低。霧化錐角80°和90°時(shí)，最后達(dá)到的摻混度相差不大，90°時(shí)摻混度達(dá)到霧化錐角研究范圍內(nèi)的最大值，而霧化錐角100°時(shí)，PDF峰值降低較多，且溫度覆蓋范圍有所變寬，摻混度開始下降。

3.5 速度矢量分布

圖6給出了穩(wěn)定摻混區(qū)的代表性速度矢量圖。摻混過程可分為4個(gè)區(qū)域:噴霧區(qū)、過渡區(qū)、穩(wěn)定區(qū)和充分發(fā)展區(qū)［40］。在每種霧化錐角下，噴霧區(qū)的流場較混亂，穩(wěn)定區(qū)形成穩(wěn)定的渦結(jié)構(gòu)影響摻混發(fā)展［24］，充分發(fā)展區(qū)渦強(qiáng)度和尺度逐漸減弱。噴霧霧化形成的空心錐造成初始液滴分布不均，噴霧對氣流的剪切卷吸形成多個(gè)渦，促進(jìn)局部液滴向整個(gè)摻混室的擴(kuò)散，在一定程度上，有利于改善溫降效果和摻混效果。霧化錐角40°和60°時(shí)，噴霧霧錐外緣形成規(guī)則的多渦對，促進(jìn)中心的液滴向其他區(qū)域擴(kuò)散。但由于渦尺度不夠大，不足以促使整個(gè)截面的溫度趨于均勻，只能作用于小區(qū)域，使溫度局部均勻，這也造成了較小霧化錐角時(shí)摻混度較低。霧化錐角增大后，相對規(guī)則的多渦結(jié)構(gòu)演化為更為混亂的大小不等的多尺度渦結(jié)構(gòu)，更有利于帶動(dòng)整個(gè)截面的摻混，摻混度較高。

圖6 不同霧化錐角下穩(wěn)定摻混區(qū)的截面速度矢量圖Fig.6 Velocity vectors in stable mixing regions under different spray angles

4 結(jié)論

(1)在霧化錐角40°～100°研究范圍內(nèi)，隨著霧化錐角增大，溫降效果逐漸提高，霧化角大于90°后，提高速率變慢;摻混度先增大，霧化錐角大于90°后，摻混度略有降低。

(2)在較小霧化錐角下，溫度的PDF值較小，且相差不大，截面溫度覆蓋范圍較寬，造成了摻混度較低;在中等霧化錐角下，溫度的PDF峰值較大，溫度相近區(qū)域較大，摻混度較高，90°時(shí)達(dá)到霧化錐角研究范圍內(nèi)的最大值;霧化錐角繼續(xù)增大，PDF峰值降低較多，且截面溫度覆蓋范圍變寬，摻混度下降。

(3)小霧化錐角摻混穩(wěn)定時(shí)，產(chǎn)生對稱的多渦對結(jié)構(gòu);而較大霧化錐角時(shí)，產(chǎn)生混亂的大小不等的多尺度渦結(jié)構(gòu)，有利于促進(jìn)整個(gè)摻混截面的溫度均勻性分布。

［1］繆萬波，夏智勛，郭健，等.金屬/水反應(yīng)沖壓發(fā)動(dòng)機(jī)理論性能計(jì)算與分析［J］.推進(jìn)技術(shù)，2005，12(6):563－566.

［2］Foote J P，Thompson B R，Lineberry J T.Combustion of aluminum with steam for underwater propulsion［C］//Roy G.Advances in chemical propulsion.Boca Raton，F(xiàn)L，CRC Press，2002:133－146.

［3］Risha G A，Huang Y，Yetter R A，et al.Combustion of aluminum particles with steam and liquid water［R］.AIAA 2006－1154.

［4］Risha G A，Sabourin J L，Yang V，et al.Combustion and conversion efficiency of nano aluminum－water mixtures［J］.Combust.Sci.Technol.，2008，180(12):2127－2142.

［5］Huang Y，Risha G A，Yang V，et al.Effect of particle size on combustion of aluminum particle dust in air［J］.Combust.Flame，2009，156(1):5－13.

［6］李芳，張為華，張煒，等.水反應(yīng)金屬燃料能量特性分析［J］.固體火箭技術(shù)，2005，28(4):256－259.

［7］劉冠鵬，郭效德，段紅珍，等.鎂水反應(yīng)的實(shí)驗(yàn)研究及機(jī)理初探［J］.固體火箭技術(shù)，2008，31(5):497－500.

［8］周星，張煒，李是良.鎂粉的高溫水反應(yīng)特性研究［J］.固體火箭技術(shù)，2009，32(3):302－305.

［9］Miller T F，Herr J D.Green rocket propulsion by reaction of Al and Mg powders and water［R］.AIAA 2004－4037.

［10］Miller T F，Garza A B.Finite rate calculations of magnesium combustion in vitiated oxygen and steam atmospheres［R］.AIAA 2006－4066.

［11］羅凱，黨建軍，王育才，等.金屬/水反應(yīng)水沖壓發(fā)動(dòng)機(jī)系統(tǒng)性能估算［J］.推進(jìn)技術(shù)，2004，25(6):495－498.

［12］朱千穩(wěn)，劉冰，夏志勛，等.噴嘴特性對水沖壓發(fā)動(dòng)機(jī)性能的影響［J］.導(dǎo)彈與航天運(yùn)載技術(shù)，2010，6(310):11－13.

［13］Huang L Y，Xia Z X，Hu J X，et al.Performance study of a water ramjet engine［J］.Sci.China，Ser.E，2011，54(4):877－882.

［14］霍東興，何國強(qiáng)，陳林泉，等.進(jìn)水方式對水沖壓發(fā)動(dòng)機(jī)性能的影響［J］.固體火箭技術(shù)，2010，33(1):31－33.

［15］黃利亞，夏志勛，胡建新，等.鎂基水沖壓發(fā)動(dòng)機(jī)燃燒室長度及進(jìn)水距離設(shè)計(jì)方法［J］.航空學(xué)報(bào)，2010，32(1):67－74.

［16］黃利亞，夏志勛，張為華，等.水沖壓發(fā)動(dòng)機(jī)試驗(yàn)水燃比選擇方法［J］.航空學(xué)報(bào)，2010，31(9):1740－1745.

［17］田維平，蔡體敏，陸賀建，等.水沖壓發(fā)動(dòng)機(jī)熱力計(jì)算［J］.固體火箭技術(shù)，2006，29(2):95－98.

［18］安慶友.金屬/水反應(yīng)燃料沖壓發(fā)動(dòng)機(jī)熱力計(jì)算與仿真［D］.西安:西北工業(yè)大學(xué)，2007.

［19］Yang Y J，He M G.A theoretical investigation of thermodynamic performance for a ramjet based on a magnesium－water reaction［J］.Proc.IMechE，Part M:J.Engineering for the Maritime Environment，2010，224(1):61－72.

［20］薛英.固體發(fā)動(dòng)機(jī)補(bǔ)燃室內(nèi)燃?xì)鈸交炖碚摷皩?shí)驗(yàn)研究［D］.西安:西北工業(yè)大學(xué)，2003.

［21］田維平，蔡體敏，孫展鵬，等.水沖壓發(fā)動(dòng)機(jī)摻混燃燒數(shù)值分析［J］.固體火箭技術(shù)，2007，30(1):9－11.

［22］繆萬波，夏智勛，胡建新，等.金屬/水反應(yīng)沖壓發(fā)動(dòng)機(jī)內(nèi)流場數(shù)值模擬［J］.推進(jìn)技術(shù)，2007，28(2):186－189.

［23］Salewski M，F(xiàn)uchs L.Dispersion of circular，non－circular，and swirling spray jets in crossflow［C］//Lamballais E，F(xiàn)riedrich R，Geurts B J，et al.Direct and large－eddy simulation VI.Dordrecht，Springer，2006:319－328.

［24］Bai B F，Zhang H B，Liu L，et al.Experimental study on turbulent mixing of spray droplets in crossflow［J］.Exp.Therm.Fluid Sci.，2009，33(6):1012－1020.

［25］Choi J，Lee S，Bae C.Spray and flow－field interaction of gasoline direct injection［J］.Atomization Spray，2004，14:159－174.

［26］Pan?o M R O，Moreira A L N.Experimental characterization of an intermittent gasoline spray impinging under crossflow conditions［J］.Atomization Spray，2005，15:201－222.

［27］Uthuppan J，Aggarwal S K，Grinstein F F，et al.Particle dispersion in a transitional axisymmetric jet:a numerical simulation［J］.AIAA J.，1994，32(10):2004－2014.

［28］Ghosh S，Hunt J C R.Spray jets in a cross－flow［J］.J.Fluid Mech.，1998，365:109－136.

［29］Ariyapadi S，Berruti F，Balachandar R.Effect of cross flow on the spray characteristics of an industrial feed nozzle［J］.Chem.Eng.Comm.，2003，190(12):1681－1704.

［30］Shi H，Kleinstreuer C.Simulation and analysis of highspeed droplet spray dynamics［J］.J.Fluids Eng.，2007，129(5):621－633.

［31］Nouri J M，Whitelaw J H.Gasoline sprays in uniform crossflow［J］.Atomization Spray，2007，17(7):621－640.

［32］Lin J Z，Qian L J，Xiong H B，et al.Effects of operating conditions on droplet deposition onto surface of atomization impinging spray［J］.Surf.Coat.Technol.，2009，203(12):1733－1740.

［33］Herrmann M.The influence of density ratio on the primary atomization of a turbulent liquid jet in crossflow［J］.Proc.Combust.Inst.，2011，33(2):2079－2088.

［34］孫慧娟，張海濱，劉利，等.氣體－顆粒摻混的數(shù)值模擬研究［J］.工程熱物理學(xué)報(bào)，2009，30(2):241－245.

［35］孫慧娟，張海濱，劉利，等.噴嘴入射方式對氣體－液滴群摻混的影響［J］.固體火箭技術(shù)，2009，32(5):492－495，510.

［36］Schmidt D P，Nouar I，Senecal P K，et al.Pressure－swirl atomization in the near field［R］.SAE－1999－01－0496，1999.

［37］O＇Rourke P J，Amsden A A.The Tab method for numerical calculation of spray droplet breakup［R］.SAE－872089，1987.

［38］O＇Rourke P J.Statistical properties and numerical implementation of a model for droplet dispersion in turbulent gas［J］.J.Comput.Phys.，1989，83(2):345－360.

［39］Bai C，Gosman A D.Development of methodology for spray impingement simulation［R］.SAE－950283，1995.

［40］劉利.橫流中噴霧摻混特性［D］.西安:西安交通大學(xué)，2010.