低溫流體閃蒸噴霧研究

羅 盟，王鐵巖，張青松，王 楠，陳 慧，李文斌

(北京宇航系統(tǒng)工程研究所，北京 100076)

0 引言

液體在真空噴射過(guò)程中，由于壓力驟減(低于液體溫度對(duì)應(yīng)的飽和蒸汽壓)，液體處于過(guò)熱狀態(tài)(熱力學(xué)亞穩(wěn)定狀態(tài))，在外界小擾動(dòng)下，其迅速回歸到穩(wěn)定狀態(tài)，并伴隨復(fù)雜的物理相變，產(chǎn)生“爆炸”式的霧化、汽化現(xiàn)象稱(chēng)為閃蒸。對(duì)于閃蒸研究可以追溯到1960年，Brown等[1]開(kāi)展了水和Freon-11介質(zhì)的閃蒸試驗(yàn)研究，并指出閃蒸存在臨界過(guò)熱條件。之后，Shepherd等[2]研究了過(guò)熱極限下液滴表面閃蒸的泰勒不穩(wěn)定。Kita-mura等[3]、Park等[4]、Cleary等[5]研究提出了閃蒸的不同模式。Adachi等[6]在試驗(yàn)基礎(chǔ)上提出了過(guò)熱蒸發(fā)速率的半經(jīng)驗(yàn)計(jì)算模型。Lamanna等[7]、Luo等[8]試驗(yàn)研究了低溫流體的閃蒸特性，提出了判定完全閃蒸的無(wú)量綱參數(shù)。國(guó)內(nèi)嚴(yán)俊杰等[9]、周致福等[10]、季璨等[11]、Guo等[12]在閃蒸換熱特性、噴霧場(chǎng)特性等方面開(kāi)展了大量的工作，對(duì)閃蒸過(guò)程特性的認(rèn)識(shí)具有重要的意義。

由于閃蒸過(guò)程涉及復(fù)雜相變過(guò)程，其理論并不成熟，目前仍缺少合適的數(shù)值計(jì)算模型。并且，在航天工程應(yīng)用方面閃蒸現(xiàn)象具有較大的破壞作用，如歐洲阿里安火箭Aestus發(fā)動(dòng)機(jī)曾出現(xiàn)過(guò)推進(jìn)劑閃蒸“結(jié)冰”導(dǎo)致的點(diǎn)火燃燒問(wèn)題[13]。鑒于此，本文在Zuo等[14]、Schmehl等[15]、Ramcke等[16]研究基礎(chǔ)上，建立了低溫流體閃蒸計(jì)算模型，開(kāi)展了低溫流體的閃蒸過(guò)程氣-液兩相流仿真研究，并用試驗(yàn)結(jié)果對(duì)其進(jìn)行了驗(yàn)證，為航天工程應(yīng)用提供參考。

1 閃蒸噴霧數(shù)學(xué)物理模型

噴霧閃蒸可認(rèn)為是由大量的單個(gè)液滴閃蒸組成。液滴閃蒸是一個(gè)涉及相變、傳熱傳質(zhì)的復(fù)雜瞬態(tài)問(wèn)題。閃蒸發(fā)生時(shí)，液滴的傳熱傳質(zhì)機(jī)制主要包括內(nèi)部過(guò)熱引起的核沸騰傳熱和外部的對(duì)流、熱傳導(dǎo)、輻射等傳熱。在較低的過(guò)熱條件下，外部傳熱具有重要作用，但在較高的過(guò)熱度下，內(nèi)部過(guò)熱占據(jù)主導(dǎo)地位。因此，完整的閃蒸模型應(yīng)同時(shí)考慮液滴內(nèi)部過(guò)熱與外部傳熱。

1.1 過(guò)熱汽化換熱

研究推進(jìn)劑噴射閃蒸過(guò)程，考慮液滴完整真實(shí)的傳熱傳質(zhì)過(guò)程是十分復(fù)雜并且耗費(fèi)資源的。因此，本模型中將液滴內(nèi)部核沸騰作簡(jiǎn)化處理，采用Adachi等[6]提出的半經(jīng)驗(yàn)參數(shù)模型。該模型假設(shè)過(guò)熱引起的換熱具有對(duì)流換熱的形式，并通過(guò)修正有效傳熱系數(shù)來(lái)考慮液滴內(nèi)部過(guò)熱影響。其模型如下

(1)

(2)

式中，αf為內(nèi)部有效傳熱系數(shù)，Ap為液滴表面積，Tp為液滴溫度，Tb為液滴飽和溫度，L(Tb)為飽和溫度下的汽化潛熱,ΔT=Tp-Tb。

1.2 外部導(dǎo)熱和對(duì)流換熱

在較低的過(guò)熱條件下，閃蒸過(guò)程中液滴外部氣體導(dǎo)熱和對(duì)流換熱占據(jù)重要的作用。在液滴蒸發(fā)過(guò)程準(zhǔn)靜態(tài)、蒸發(fā)過(guò)程中始終保持理想球形，和液滴內(nèi)部溫度均勻分布的假設(shè)下，Zuo等[14]給出了液滴導(dǎo)熱和對(duì)流換熱產(chǎn)生的蒸發(fā)質(zhì)量流率

(3)

式中，r0為液滴初始半徑，λ為導(dǎo)熱系數(shù)，h為氣體焓值，cp為液滴比熱，Nu*為考慮沸騰的修正努塞爾數(shù)，其表達(dá)式為

(4)

Respe和Prspe為特定溫度Tspe下的雷諾數(shù)和普朗特?cái)?shù)。BT為Spading數(shù)，其表達(dá)式為

(5)

氣-液混合熱力學(xué)參數(shù)和輸運(yùn)特性參數(shù)均用特定溫度計(jì)算，該特定溫度采用Sparrow和Gregg提出的1/3規(guī)則[17]計(jì)算

(6)

式中，T∞為環(huán)境溫度，Tb為飽和溫度。

1.3 輻射換熱

因?yàn)榈蜏匾后w溫度較低，本文考慮了外界環(huán)境輻射對(duì)液滴表面蒸發(fā)的影響。由輻射換熱引起的液體蒸發(fā)質(zhì)量流率為

(7)

式中，ε為液滴表面發(fā)射系數(shù)，本計(jì)算取0.96，σ為斯蒂芬-波爾茲曼常數(shù)，Ap為液滴表面積，Tamb為環(huán)境溫度。

1.4 液滴能量方程

閃蒸過(guò)程中，液滴內(nèi)部溫度變化由過(guò)熱蒸發(fā)決定。液滴的能量方程為

(8)

1.5 液滴運(yùn)動(dòng)模型

運(yùn)動(dòng)顆粒軌跡的描述通過(guò)求解牛頓運(yùn)動(dòng)方程。在實(shí)際噴霧閃蒸霧化過(guò)程中，受到復(fù)雜的外力作用，如曳力、重力、壓力梯度引起的體積力、浮力等。本文忽略液滴之間的碰撞和融合，并認(rèn)為氣動(dòng)阻力是其受力主導(dǎo)因素。阻力模型采用修正的Schiller-Naumann模型[18]。液滴運(yùn)動(dòng)方程為

(9)

CD=

(10)

式中，vp和v分別為液滴和氣體速度矢量，BM為Spalding傳質(zhì)數(shù)，CD為阻力系數(shù)，μ為流體的黏度，ρp為液滴密度，dp為液滴直徑，Re為雷諾數(shù)。

1.6 氣-液兩相相互作用

低溫推進(jìn)劑真空噴射屬于兩相流動(dòng)。本文在歐拉-拉格朗日框架下采用CFD-DPM方法對(duì)其進(jìn)行仿真計(jì)算。通過(guò)求解三維非穩(wěn)態(tài)可壓縮N-S方程來(lái)描述連續(xù)相(氣體)，通過(guò)求解牛頓運(yùn)動(dòng)方程追蹤離散相運(yùn)動(dòng)軌跡(液滴噴霧)，氣-液兩相之間的相互作用通過(guò)在連續(xù)相方程中添加相應(yīng)的源項(xiàng)來(lái)考慮。三維非穩(wěn)態(tài)N-S方程可參考文獻(xiàn)[19]。兩相之間的質(zhì)量、動(dòng)量、能量源項(xiàng)分別如式(11)～(13)

(11)

(12)

(13)

式中，Vcell為計(jì)算網(wǎng)格的體積。

1.7 物性參數(shù)

考慮到低溫推進(jìn)劑閃蒸過(guò)程中溫度變化范圍較大(忽略壓力變化)，常物性參數(shù)計(jì)算帶來(lái)的影響較大，如圖1所示。因此本文基于NIST 9.1數(shù)據(jù)庫(kù)，采用多項(xiàng)式擬合來(lái)描述基于溫度變化的氣液物性參數(shù)(如密度、比焓，比熱、黏性、導(dǎo)熱率等)

圖1 模擬流體的物性Fig.1 Properties of the simulated fluids

φ(T)=a0+a1T+a2T2+a3T3+a4T4

(14)

通過(guò)以上分析建立了低溫流體閃蒸噴霧計(jì)算模型。本文將此模型植入FLUENT軟件中，當(dāng)液滴溫度高于其周?chē)h(huán)境壓力對(duì)應(yīng)的飽和溫度時(shí)，調(diào)用閃蒸模型；當(dāng)溫度低于飽和溫度時(shí)，應(yīng)用FLUENT的擴(kuò)散控制蒸發(fā)模型(D2蒸發(fā)規(guī)律)。

2 試驗(yàn)簡(jiǎn)介

圖2給出了試驗(yàn)裝置圖。試驗(yàn)研究的是低溫推進(jìn)劑經(jīng)過(guò)單孔直噴嘴噴入低壓環(huán)境的過(guò)程。方形真空室為三面開(kāi)玻璃窗結(jié)構(gòu)，上端設(shè)計(jì)一個(gè)液氮儲(chǔ)罐，噴嘴、控制閥等浸泡在里面。試驗(yàn)時(shí)，通過(guò)向小液氮罐加壓，控制其中液氮溫度。推進(jìn)劑管路穿過(guò)小液氮罐，通過(guò)與液氮的換熱保證合適的推進(jìn)劑噴射溫度。沿燃燒室中線(xiàn)軸線(xiàn)方向布置T型熱電偶(直徑1.0 mm)，每?jī)烧唛g隔20 mm，用以測(cè)量噴霧沿軸線(xiàn)的溫度分布。燃燒室上布置壓力傳感器，用以測(cè)量真空室內(nèi)壓力。噴嘴頭腔有熱電偶和壓力傳感器測(cè)口，用以測(cè)量噴射流體的溫度和壓力。噴霧形態(tài)采用高速紋影技術(shù)記錄。

圖2 試驗(yàn)裝置圖Fig.2 Experimental setup

3 數(shù)值計(jì)算

3.1 計(jì)算方法

本研究基于CFD-DPM方法，在歐拉-拉格朗日框架下模擬氣-液(液滴)兩相流動(dòng)。通過(guò)求解三維非穩(wěn)態(tài)可壓縮雷諾平均N-S方程(U-RANS)計(jì)算連續(xù)相流場(chǎng)，通過(guò)求解牛頓運(yùn)動(dòng)方程追蹤離散相運(yùn)動(dòng)軌跡。氣-液兩相之間的相互作用是通過(guò)在連續(xù)相方程中添加相應(yīng)的源項(xiàng)來(lái)考慮。采用兩方程標(biāo)準(zhǔn)k-ε湍流模型預(yù)測(cè)連續(xù)相湍流流動(dòng)，采用隨機(jī)軌道追蹤模型來(lái)考慮連續(xù)相湍流效應(yīng)對(duì)離散相的作用。仿真計(jì)算采用基于密度的隱式求解器，PISO壓力-速度耦合算法，動(dòng)量方程、湍動(dòng)能、湍流耗散率方程對(duì)流項(xiàng)均采用二階迎風(fēng)格式離散，時(shí)間項(xiàng)采用一階隱式。

3.2 計(jì)算模型及邊界條件

計(jì)算模型及邊界條件如圖3所示。為節(jié)省計(jì)算資源，取模型1/4為計(jì)算域，采用對(duì)稱(chēng)邊界條件。連續(xù)相采用氣相速度入口，壓力出口，絕熱無(wú)滑移壁面邊界條件。離散相壁面邊界采用彈性碰撞的反射邊界。入口離散相液滴大小以Rosin-Rammler分布規(guī)律給定。前期對(duì)計(jì)算空間離散和時(shí)間離散進(jìn)行了敏感性分析。計(jì)算網(wǎng)格約180萬(wàn)，連續(xù)相時(shí)間步長(zhǎng)1×10-5s，離散相追蹤時(shí)間步長(zhǎng)1×10-5s。

圖3 計(jì)算域、網(wǎng)格和邊界條件Fig.3 Computational domain, mesh and boundary conditions

4 結(jié)果與分析

本節(jié)通過(guò)液氮和液氧閃蒸的試驗(yàn)結(jié)果驗(yàn)證了仿真模型。由于液氧和液氮的閃蒸規(guī)律具有一定相似性，4.3～4.5節(jié)以液氧為例對(duì)仿真結(jié)果進(jìn)行了詳細(xì)的分析。

4.1 噴霧形狀分析

圖4和圖5分別給出了仿真計(jì)算的液氧、液氮閃蒸噴霧中心截面上形態(tài)與試驗(yàn)紋影對(duì)比圖。右半邊紋影圖中的黑色線(xiàn)條是熱電偶，直徑為1.0 mm。頂端噴嘴出口直徑為0.5 mm，可視化窗口寬度約為52 mm。試驗(yàn)條件如表1所示。

表1 計(jì)算工況條件

由圖4～圖5可以看出，仿真計(jì)算較好地預(yù)測(cè)了液氧、液氮的閃蒸噴霧形態(tài)。該噴霧形態(tài)明顯不同于典型的射流噴霧(小錐形噴霧角)。閃蒸噴霧在噴嘴出口附近表現(xiàn)出很大的噴霧角和徑向穿透距離，這是由于低溫推進(jìn)劑在噴入低壓環(huán)境的劇烈相變導(dǎo)致。

(a)仿真結(jié)果 (b)試驗(yàn)紋影圖圖4 液氧噴霧Fig.4 LOx spray contour

(a)仿真結(jié)果 (b)試驗(yàn)紋影圖圖5 液氮噴霧Fig.5 LN2 spray contour

以液氧工況為例，液氧噴入低壓環(huán)境后，在噴嘴出口，由于壓力的突然降低，液氧處于過(guò)熱狀態(tài)，在外界壓力擾動(dòng)下釋放過(guò)熱回到穩(wěn)定狀態(tài)。在回歸穩(wěn)定狀態(tài)過(guò)程中是氣-液傳熱傳質(zhì)的強(qiáng)烈非平衡過(guò)程，過(guò)熱液體中會(huì)產(chǎn)生大量初始?xì)馀莺?。初始?xì)馀莺丝梢允蔷|(zhì)結(jié)核也可以是異質(zhì)結(jié)核。由分子運(yùn)動(dòng)理論可知，液體中分子團(tuán)能量的不均勻性導(dǎo)致了流體密度起伏，局部低密度區(qū)形成初始的空化核心(均質(zhì)結(jié)核核心)，而異質(zhì)初始結(jié)核核心主要由固體接觸表面或者液體內(nèi)部雜質(zhì)提供。在過(guò)熱液體中形成一個(gè)半徑為r的氣泡核心后引起系統(tǒng)的自由焓的變化(等于形成該氣泡外界所需做出的功)為[20-21]

(15)

式中，σ為氣泡表面張力，Δμ為過(guò)熱液體與氣體的化學(xué)勢(shì)差，第一項(xiàng)恒為正，表示過(guò)熱液體生成氣泡需要做表面功。

隨著r值增大，系統(tǒng)的自由焓變隨氣泡大小變化曲線(xiàn)(紅線(xiàn)為均質(zhì)結(jié)核情況，綠線(xiàn)為異質(zhì)結(jié)核情況)如圖6所示。可以看出，該曲線(xiàn)存在最大值，其對(duì)應(yīng)的氣泡半徑為臨界半徑。當(dāng)氣泡半徑大于臨界半徑時(shí)，生成氣泡引起的系統(tǒng)自由焓變是減小的(生成氣泡所需的功減小)，這有利于氣泡的生長(zhǎng)過(guò)程。在此情況下，大量氣泡生成，產(chǎn)生劇烈的核沸騰，氣泡的進(jìn)一步生長(zhǎng)使噴射流體破碎、霧化，伴隨著劇烈汽化，從而形成了如圖6所示的“爆炸式”噴霧形態(tài)。

圖6 自由焓與氣泡大小關(guān)系Fig.6 Free enthalpy versus sbubble size

4.2 噴霧溫度分析

閃蒸現(xiàn)象的一個(gè)重要研究對(duì)象是溫度分布特性。圖7和圖8分別給出了液氧、液氮噴霧軸線(xiàn)溫度分布仿真計(jì)算和試驗(yàn)結(jié)果。圖中紅色虛線(xiàn)是噴霧軸線(xiàn)上離散相液滴溫度，紅色實(shí)線(xiàn)是軸線(xiàn)上連續(xù)相氣體的溫度分布。由于試驗(yàn)中熱電偶測(cè)量的是熱電偶端頭局部區(qū)域內(nèi)噴霧的溫度，即氣-液“混合團(tuán)”的溫度。因此，本文對(duì)數(shù)值計(jì)算的液滴和氣體溫度采用1/3規(guī)則進(jìn)行了平均，以此代替計(jì)算的局部氣-液“混合團(tuán)”溫度，如圖7和圖8 中星號(hào)曲線(xiàn)所示。可以看出，仿真計(jì)算結(jié)果與試驗(yàn)吻合較好，這也驗(yàn)證了上述建立的閃蒸模型的正確性。由圖看出，低溫液體噴霧在噴嘴出口附近(Z<20 mm或Z/d<40)溫度由111 K急劇下降至85 K，降幅達(dá)26%，隨后在流場(chǎng)下游區(qū)域溫度分布平緩，降幅僅約5.8%。這間接預(yù)示了低溫液體在噴入真空燃燒室后，閃蒸現(xiàn)象隨即發(fā)生，并在噴嘴出口附近區(qū)域基本完成了過(guò)熱液體的閃蒸過(guò)程。

圖7 液氧噴霧溫度沿軸向分布Fig.7 LOx temperature distribution along the spray centerline

圖8 液氮噴霧溫度沿軸向分布Fig.8 LN2 temperature distribution along the spray centerline

4.3 噴霧粒徑變化分析

圖9和圖10分別給出了氧液滴粒徑D32(索泰爾直徑)沿噴霧軸向和徑向分布規(guī)律。圖9中的散點(diǎn)是在55 ms時(shí)刻軸線(xiàn)附近追蹤的顆粒大小分布。紅色曲線(xiàn)是0～55 ms時(shí)間段內(nèi)的索泰爾直徑沿軸線(xiàn)分布的統(tǒng)計(jì)平均結(jié)果。可以看出，在噴嘴出口附近，氧液滴尺寸呈現(xiàn)出非常大的梯度分布規(guī)律，隨后液滴尺寸變化逐漸平緩。在0～20 mm范圍內(nèi)，液滴尺寸從大約9.6 μm減小到8.5 μm，其減小幅度占整個(gè)計(jì)算域內(nèi)減小幅值的70%。粒徑變化與溫度變化規(guī)律一致，這驗(yàn)證了過(guò)熱閃蒸主要發(fā)生在噴嘴出口附近。圖中綠色點(diǎn)劃線(xiàn)是對(duì)散點(diǎn)圖進(jìn)行了數(shù)據(jù)擬合，結(jié)果表明閃蒸過(guò)程中液滴大小較符合指數(shù)變化規(guī)律。若假設(shè)液滴在徑向各個(gè)截面的分布規(guī)律相似，則通過(guò)粒徑的軸向分布規(guī)律可以大致估算整個(gè)計(jì)算域內(nèi)的蒸發(fā)量約占噴注量的34%。

圖9 液氧噴霧粒徑沿軸向分布Fig.9 LOx droplet size distribution along the spray centerline

圖10 液氧噴霧粒徑沿徑向分布Fig.10 LOx droplet size distribution in the spray radial direction

圖10給出離噴嘴出口不同軸向距離截面處，液滴大小沿徑向分布規(guī)律?？梢钥闯觯诳拷鼑婌F中心區(qū)域，液滴在各個(gè)不同截面的大小基本均勻。在同一截面上，噴霧外圍液滴尺寸沿徑向緩慢減小然后急劇增大，即較大的液滴聚集在噴霧外圍，這是由于大直徑液滴的慣性較大，受周?chē)B續(xù)相流場(chǎng)的影響較小，因此其更多地流向噴霧外圍。Munnannur等[22]以及Lasheras等[23]的兩相流噴霧試驗(yàn)和計(jì)算也得到了類(lèi)似的粒徑大小徑向分布規(guī)律。

4.4 噴霧脈動(dòng)速度分析

在追蹤液滴過(guò)程中，液滴的運(yùn)動(dòng)由液滴和連續(xù)相氣體的相互作用決定。因此，液滴的脈動(dòng)速度場(chǎng)也從一方面反映了整個(gè)噴霧流場(chǎng)的湍流程度。式(16)給出了追蹤液滴的均方根脈動(dòng)速度(RSM)

(16)

圖11給出了沿噴霧軸線(xiàn)不同截面處，液氧噴霧脈動(dòng)速度徑向分布曲線(xiàn)?？梢钥闯?，氧液滴在噴嘴附近(如Z=10 mm)表現(xiàn)出較大的速度脈動(dòng)(～4%)，并且沿徑向分布波動(dòng)較大。隨著軸向距離增加，速度脈動(dòng)大幅減小，在噴嘴下游較遠(yuǎn)處(如Z>100 mm)，脈動(dòng)速度小于1%，并且在徑向外圍分布均勻。這是由于噴嘴出口附近的劇烈閃蒸引起的連續(xù)相流場(chǎng)的較大脈動(dòng)，進(jìn)而作用于液滴，引起液滴較大的速度脈動(dòng)；而隨著軸向距離增加，閃蒸急劇減弱，整個(gè)噴霧流場(chǎng)趨于平穩(wěn)。

圖11 液氧噴霧脈動(dòng)速度沿徑向分布Fig.11 LOx velocity fluctuation along the spray radial direction

4.5 傳熱傳質(zhì)分析

由以上的分析可知，低溫流體在噴入低壓環(huán)境后，在噴嘴出口附近區(qū)域便發(fā)生了劇烈的閃蒸現(xiàn)象，導(dǎo)致大的噴霧霧化角和徑向擴(kuò)張，溫度急劇降低，液滴粒徑迅速減?。辉趪娮煜掠?，隨著距離增加，以上參數(shù)變化平緩。閃蒸的這些特性是由其傳熱傳質(zhì)決定的，本節(jié)對(duì)閃蒸過(guò)程中的傳熱傳質(zhì)影響進(jìn)行了分析。

圖12給出了液氧噴霧沿軸線(xiàn)方向的蒸發(fā)速率變化規(guī)律。由圖可知，在液氧噴入低壓環(huán)境后，強(qiáng)烈的閃蒸發(fā)生，產(chǎn)生的蒸發(fā)流量在10-9kg/s量級(jí)。之后，在較小的距離內(nèi)(Z<20 mm)，閃蒸蒸發(fā)流量迅速降至10-11kg/s量級(jí)。在此區(qū)域，對(duì)流和導(dǎo)熱引起的蒸發(fā)流量較低，在10-12～10-11kg/s量級(jí)。因此，在靠近噴嘴出口附近，閃蒸占據(jù)絕對(duì)主導(dǎo)地位，高于其他換熱至少2個(gè)數(shù)量級(jí)。隨著噴霧流向噴嘴下游，由于液滴的過(guò)熱度大幅減小，閃蒸引起的蒸發(fā)流量快速減小，而對(duì)流和導(dǎo)熱逐漸占據(jù)主導(dǎo)作用。由于試驗(yàn)前的預(yù)冷過(guò)程，真空燃燒室周?chē)跏辑h(huán)境溫度較低(約220 K)。因此，整個(gè)過(guò)程中的輻射換熱較弱。從計(jì)算結(jié)果可以看出，其產(chǎn)生的蒸發(fā)流量較低，基本維持在10-13kg/s量級(jí)。

圖12 液氧噴霧蒸發(fā)速率沿軸向分布Fig.12 LOx sprays evaporation mass flow rates along the spray centerline

4.6 閃蒸固態(tài)相變

值得一提的是，在非常大的過(guò)熱條件下，試驗(yàn)中觀察到了液氮噴霧復(fù)雜的氣-液-固相變。固態(tài)氮更多地聚集在第一個(gè)熱電偶上(距離噴嘴出口20 mm)，并且靠近噴霧中心附近，如圖13所示。這也從一方面驗(yàn)證了本文數(shù)值模擬結(jié)果的合理性，即過(guò)熱蒸發(fā)在靠近噴嘴區(qū)域占據(jù)主導(dǎo)地位，從而導(dǎo)致該區(qū)域換熱量非常大，產(chǎn)生了復(fù)雜的相變過(guò)程。在液氧噴射試驗(yàn)過(guò)程中，由于試驗(yàn)設(shè)備可提供的真空條件受限(高于40 mbar)，而氧氣三相點(diǎn)壓力為1.46 mbar，遠(yuǎn)低于試驗(yàn)真空條件，因此并未觀察到固態(tài)氧。相比而言，氮?dú)馊帱c(diǎn)壓力約為125 mbar，較容易觀察其凝固相變過(guò)程。

圖13 閃蒸過(guò)程中氮凝固現(xiàn)象Fig.13 Nitrogen solidification during the flashing

閃蒸的強(qiáng)烈換熱引起的固態(tài)相變會(huì)給火箭動(dòng)力系統(tǒng)帶來(lái)潛在的危害，在實(shí)際工程設(shè)計(jì)中應(yīng)予以重視。如發(fā)動(dòng)機(jī)高空二次點(diǎn)火過(guò)程中，推進(jìn)劑噴射閃蒸可能會(huì)產(chǎn)生固態(tài)推進(jìn)劑塊，導(dǎo)致點(diǎn)火條件(點(diǎn)火處局部氧/燃比)偏離設(shè)計(jì)工況，帶來(lái)點(diǎn)火可靠性問(wèn)題；或大量的固態(tài)推進(jìn)劑團(tuán)聚集在發(fā)動(dòng)機(jī)壁面附近，推力室點(diǎn)火時(shí)可能產(chǎn)生爆燃，引起燃燒不穩(wěn)定，甚至造成燃燒室結(jié)構(gòu)件破壞。

5 結(jié)論

本文建立了低溫流體噴射閃蒸計(jì)算模型，在歐拉-拉格朗日框架下采用CFD-DPM方法開(kāi)展了低溫液體真空噴射兩相流仿真計(jì)算，并進(jìn)行了試驗(yàn)驗(yàn)證。結(jié)果表明該模型可以很好地預(yù)測(cè)閃蒸噴霧特性，包括閃蒸噴霧形態(tài)和溫度特性。研究結(jié)果表明：

1)低溫流體噴入低壓環(huán)境后，在噴嘴出口處便發(fā)生閃蒸現(xiàn)象，呈現(xiàn)出很大的噴霧角。

2)噴霧溫度、液滴大小在噴嘴附近呈現(xiàn)非常大的梯度變化。在出噴口較短的距離，噴霧溫度便降至接近環(huán)境壓力對(duì)應(yīng)的飽和溫度，隨后溫度變化緩慢。

3)過(guò)熱閃蒸主要發(fā)生在噴嘴出口附近區(qū)域(x/d～40)，在此區(qū)域過(guò)熱閃蒸占據(jù)絕對(duì)主導(dǎo)地位(高出其他換熱2個(gè)數(shù)量級(jí))；噴嘴下游，對(duì)流換熱與熱傳導(dǎo)影響占據(jù)主導(dǎo)地位；由于環(huán)境溫度較低，整個(gè)過(guò)程中輻射換熱量較小(低于其他換熱1～2個(gè)數(shù)量級(jí))。

4)閃蒸帶來(lái)的強(qiáng)烈換熱可能導(dǎo)致噴嘴附近發(fā)生復(fù)雜的氣-液-固相變，會(huì)對(duì)發(fā)動(dòng)機(jī)可靠點(diǎn)火帶來(lái)潛在風(fēng)險(xiǎn)。