尹 婷，吳高楊，聶萬(wàn)勝

（裝備學(xué)院，北京101416）

0 引言

在液體火箭發(fā)動(dòng)機(jī)啟動(dòng)及穩(wěn)態(tài)工作過(guò)程中，時(shí)常伴隨有不穩(wěn)定燃燒現(xiàn)象發(fā)生，表現(xiàn)為燃燒室內(nèi)壓力的周期性振蕩[1]。不穩(wěn)定燃燒不僅影響推進(jìn)劑燃燒過(guò)程，而且激發(fā)機(jī)理復(fù)雜，燃燒過(guò)程中的霧化、蒸發(fā)、燃燒各子環(huán)節(jié)與燃燒室聲學(xué)結(jié)構(gòu)的耦合都會(huì)引起不穩(wěn)定燃燒[2]。因而，研究不穩(wěn)定燃燒問(wèn)題具有極大挑戰(zhàn)性，探索液滴蒸發(fā)燃燒過(guò)程對(duì)壓力振蕩的動(dòng)態(tài)響應(yīng)特性及其規(guī)律，是揭示不穩(wěn)定燃燒激勵(lì)機(jī)理、解決液體火箭發(fā)動(dòng)機(jī)不穩(wěn)定燃燒問(wèn)題的前提。有研究指出，液滴蒸發(fā)過(guò)程對(duì)壓力振蕩的動(dòng)態(tài)響應(yīng)特性極有可能是導(dǎo)致不穩(wěn)定燃燒的重要原因[3]。因此，解決不穩(wěn)定燃燒問(wèn)題、探索其激發(fā)機(jī)理，首先要從蒸發(fā)過(guò)程開(kāi)始研究。

國(guó)外率先開(kāi)展了壓力振蕩環(huán)境下液滴蒸發(fā)特性規(guī)律的研究，在理論分析與數(shù)值仿真方面[4-6]，主要考查了壓力振蕩引發(fā)的液滴蒸發(fā)速率，液滴表面溫度和物性參數(shù)的變化；而在實(shí)驗(yàn)研究方面[7-9]，主要采用揚(yáng)聲器和麥克風(fēng)的相互配合在密閉聲波導(dǎo)管內(nèi)形成微小的聲波振蕩，從而研究振蕩環(huán)境下液滴蒸發(fā)過(guò)程，對(duì)實(shí)驗(yàn)設(shè)備和技術(shù)要求較高。相較于國(guó)外數(shù)十年的研究，國(guó)內(nèi)對(duì)壓力振蕩環(huán)境下液滴蒸發(fā)特性的研究則顯得相對(duì)落后。尤其在實(shí)驗(yàn)研究方面，由于技術(shù)水平和實(shí)驗(yàn)條件的限制，相關(guān)研究還處于起步階段。

本文針對(duì)偏二甲肼液滴，設(shè)計(jì)了壓力振蕩環(huán)境實(shí)驗(yàn)系統(tǒng)，開(kāi)展了不同振蕩環(huán)境下液滴蒸發(fā)動(dòng)態(tài)響應(yīng)特性研究，建立了壓力振蕩環(huán)境液滴非穩(wěn)態(tài)蒸發(fā)模型，考察了振蕩頻率、幅值、燃燒室壓力及溫度等因素對(duì)液滴蒸發(fā)特性的影響規(guī)律。

1 實(shí)驗(yàn)系統(tǒng)設(shè)計(jì)

壓力振蕩環(huán)境液滴蒸發(fā)實(shí)驗(yàn)系統(tǒng)的工作原理是高壓氣源通過(guò)旋轉(zhuǎn)換氣閥產(chǎn)生振蕩氣流，不斷向燃燒室內(nèi)進(jìn)行充、放氣過(guò)程，從而持續(xù)形成壓力振蕩環(huán)境，通過(guò)高速攝像機(jī)記錄該環(huán)境下液滴蒸發(fā)過(guò)程，分析其特性規(guī)律。該實(shí)驗(yàn)系統(tǒng)由壓力振蕩產(chǎn)生子系統(tǒng)、燃燒室子系統(tǒng)、高速攝影子系統(tǒng)、燃料液滴生成子系統(tǒng)及數(shù)據(jù)采集控制子系統(tǒng)組成。

其中，壓力振蕩產(chǎn)生子系統(tǒng)是實(shí)驗(yàn)成功的關(guān)鍵。它包含了氣瓶、供氣室、交流伺服電機(jī)及其驅(qū)動(dòng)器、聯(lián)軸器、旋轉(zhuǎn)換氣閥以及穩(wěn)流艙。接通電源后，交流伺服驅(qū)動(dòng)器驅(qū)動(dòng)電機(jī)旋轉(zhuǎn)，帶動(dòng)了與電機(jī)同軸的旋轉(zhuǎn)換氣閥轉(zhuǎn)動(dòng)。換氣閥閥盤(pán)上等間距的分布了26個(gè)通氣孔，經(jīng)過(guò)供氣室過(guò)濾減壓作用的氣源通入旋轉(zhuǎn)換氣閥后，形成以一定頻率間歇性排出的氣流，經(jīng)過(guò)穩(wěn)流艙的穩(wěn)流作用后，通入燃燒室形成壓力振蕩環(huán)境。

燃料液滴生成子系統(tǒng)用于偏二甲肼液滴的生成，由燃料貯箱、噴注器、懸掛針頭、電磁閥和流量計(jì)組成。實(shí)驗(yàn)時(shí)通過(guò)N2提供擠壓動(dòng)力，將燃料貯箱內(nèi)的燃料擠出，經(jīng)噴注器在懸掛針頭的末端形成液滴，通過(guò)流量計(jì)確定燃料的流量，通過(guò)計(jì)算機(jī)控制電磁閥的通斷及時(shí)間間隔，從而控制生成液滴的大小和頻率。

燃燒室子系統(tǒng)是一個(gè)由10 mm厚的不銹鋼板焊接而成的長(zhǎng)方體密閉容器，其內(nèi)部體積為150 mm×150 mm×280 mm，分為頂蓋、身部和底蓋。頂蓋上留有通氣孔及測(cè)壓孔，與單液滴生成子系統(tǒng)的噴注器直接相連；身部側(cè)壁上留有觀察窗，方便進(jìn)行觀察和拍攝；底部裝有加熱器，能夠?qū)⑷紵已杆偌訜嶂翆?shí)驗(yàn)所需溫度。

高速攝影子系統(tǒng)用于拍攝記錄液滴蒸發(fā)過(guò)程，方便實(shí)驗(yàn)數(shù)據(jù)后期處理，通過(guò)計(jì)算機(jī)控制高速攝像機(jī)的拍攝范圍和頻率。

數(shù)據(jù)采集控制子系統(tǒng)由溫度、壓力傳感器和控制線路等組成，通過(guò)傳感器的測(cè)量，可以實(shí)時(shí)測(cè)量并記錄燃燒室內(nèi)振蕩壓力及溫度的變化。同時(shí)，操作計(jì)算機(jī)從而控制各管路閥門(mén)開(kāi)關(guān)與電磁閥的通斷，保證實(shí)驗(yàn)時(shí)各子系統(tǒng)間的緊密銜接與同步。

2 非穩(wěn)態(tài)蒸發(fā)模型

2.1 控制方程

對(duì)于壓力振蕩環(huán)境下液滴蒸發(fā)過(guò)程，由于壓力振蕩使得物性參數(shù)始終處于非穩(wěn)態(tài)變化中，因而需要建立適用于壓力振蕩環(huán)境的非穩(wěn)態(tài)蒸發(fā)模型。

連續(xù)性方程

式中：ρg和vg分別為氣體密度和速度；下標(biāo) g表示氣相。

能量守恒方程

式中：Cpg為氣相比熱容；kg為導(dǎo)熱系數(shù)。

組分守恒方程

式中：D為擴(kuò)散系數(shù)；Y為燃料蒸氣質(zhì)量分?jǐn)?shù)；W為組分源項(xiàng)。

式（1）～式（3）即為氣相流動(dòng)控制方程。液相控制方程具有類(lèi)似的表達(dá)形式，連續(xù)性方程

能量守恒方程

式中：下標(biāo)l表示液相；QT表示能量源項(xiàng)，無(wú)化學(xué)反應(yīng)時(shí)該項(xiàng)為零。

組分守恒方程

氣體滿(mǎn)足理想狀態(tài)方程

2.2 邊界條件

各區(qū)域邊界條件如下：

氣相遠(yuǎn)場(chǎng)處

液滴中心處的邊界條件為

氣液界面處于平衡狀態(tài)，滿(mǎn)足溫度分布的連續(xù)性，即

假設(shè)氣體為理想氣體，則有，

氣液界面處滿(mǎn)足能量守恒，則

式中：Gs為液滴表面蒸發(fā)流率；L為蒸發(fā)潛熱；Ts為液滴表面溫度；ΔTs為Δt時(shí)刻前液滴表面溫度。

液滴表面蒸發(fā)流率由燃料和介質(zhì)兩部分組成

式中：Go為液滴表面介質(zhì)蒸發(fā)流率；Gf為燃料蒸發(fā)流率；Gs則為液滴表面蒸發(fā)流率。

液滴燃料蒸發(fā)流率表達(dá)式為

由于介質(zhì)不聚集在液滴表面，因而，聯(lián)立式（14）和（15），可得到液滴表面蒸發(fā)流率

3 實(shí)驗(yàn)結(jié)果及分析

壓力振蕩環(huán)境下液滴蒸發(fā)實(shí)驗(yàn)考察了5 Hz，13 Hz，26 Hz和130 Hz四個(gè)不同振蕩頻率下的液滴蒸發(fā)過(guò)程，同時(shí)改變振幅、燃燒室平均壓力、溫度等工況條件開(kāi)展了不同工況下的實(shí)驗(yàn)研究，獲得了豐富的實(shí)驗(yàn)數(shù)據(jù)，下面選取其中的26組實(shí)驗(yàn)工況進(jìn)行分析。

表1給出了26組壓力振蕩環(huán)境下液滴蒸發(fā)實(shí)驗(yàn)工況參數(shù)。其中，“最大振幅”指壓力振蕩波峰處到波節(jié)處的最大值，按照最大振幅占燃燒室平均壓力的百分比來(lái)考察；采樣頻率指壓力傳感器采集燃燒室內(nèi)壓力變化的頻率，為了保證測(cè)量精度，表中各工況的采樣頻率均為1 000 Hz，遠(yuǎn)大于燃燒室壓力振蕩頻率，因此能夠很好地反映燃燒室內(nèi)振蕩壓力的變化。

表1 壓力振蕩環(huán)境液滴蒸發(fā)實(shí)驗(yàn)工況參數(shù)Tab.1 Parameters of droplet evaporation in pressure oscillation environment

圖1為壓力振蕩環(huán)境下不同液滴蒸發(fā)過(guò)程。由圖可知，液滴尺寸隨著蒸發(fā)過(guò)程的進(jìn)行而逐步縮小。同時(shí)，部分工況下液滴會(huì)在噴注器末端左右擺動(dòng)，發(fā)生形變，如圖1（a）。這說(shuō)明液滴在蒸發(fā)過(guò)程中受到氣流影響，氣流會(huì)加速液滴蒸發(fā)過(guò)程，因而無(wú)法判斷液滴在蒸發(fā)過(guò)程中的變化是來(lái)自于氣流，還是壓力振蕩環(huán)境的影響；而圖1（b）和1（c）中的液滴則能夠非常穩(wěn)定地懸掛在噴注器末端，既沒(méi)有發(fā)生形變，也沒(méi)有出現(xiàn)位置的偏移，可以認(rèn)為此時(shí)液滴并沒(méi)有受到氣流影響。因此，在壓力振蕩環(huán)境下液滴蒸發(fā)的實(shí)驗(yàn)中，為了獲得更加準(zhǔn)確的實(shí)驗(yàn)數(shù)據(jù)，選擇蒸發(fā)過(guò)程中沒(méi)有出現(xiàn)擺動(dòng)和形變的液滴進(jìn)行分析。

圖1 壓力振蕩環(huán)境下液滴蒸發(fā)過(guò)程Fig.1 Droplet evaporation process in pressure oscillation environment

3.1 液滴生存時(shí)間

為了與壓力振蕩環(huán)境下液滴蒸發(fā)過(guò)程進(jìn)行對(duì)比，本文選取了3個(gè)典型實(shí)驗(yàn)工況，利用上述實(shí)驗(yàn)系統(tǒng)，關(guān)閉旋轉(zhuǎn)換氣閥，使燃燒室內(nèi)壓力維持恒定，開(kāi)展了壓力恒定環(huán)境下液滴蒸發(fā)實(shí)驗(yàn)研究，實(shí)驗(yàn)結(jié)果如表2所示。由表可知，基本相同的工況條件下，壓力振蕩環(huán)境中液滴生存時(shí)間明顯小于壓力恒定環(huán)境中液滴生存時(shí)間，盡管實(shí)驗(yàn)形成的壓力振蕩幅度有限，但液滴生存時(shí)間普遍減少了15 s以上。因此，壓力振蕩環(huán)境明顯促進(jìn)了液滴蒸發(fā)，大幅縮短了液滴生存時(shí)間。

表2 壓力振蕩環(huán)境及壓力恒定環(huán)境液滴蒸發(fā)典型實(shí)驗(yàn)工況參數(shù)對(duì)比Tab.2 Typical experimental parameters of droplet evaporation inpressure oscillation environment and constant environment

3.2 液滴蒸發(fā)速率

壓力振蕩環(huán)境不但減小了液滴生存時(shí)間，對(duì)液滴蒸發(fā)速率也具有不可忽視的作用。根據(jù)高速攝像機(jī)記錄的液滴蒸發(fā)過(guò)程，測(cè)量到液滴直徑隨時(shí)間的變化，從而計(jì)算獲得液滴蒸發(fā)速率實(shí)驗(yàn)值。圖2分別給出了工況Test 9（377 K，0.144 MPa）和Test 6（393 K，0.378 MPa）液滴蒸發(fā)速率與振蕩壓力波形。由圖可知，兩種工況下液滴蒸發(fā)速率分別可達(dá)到35 mg/s和27 mg/s，而壓力恒定環(huán)境下相同工況的液滴蒸發(fā)速率不超過(guò)0.1 mg/s，遠(yuǎn)小于壓力振蕩環(huán)境中的液滴蒸發(fā)速率。如圖3所示，為壓力恒定環(huán)境下液滴蒸發(fā)速率，工況分別為 377 K，0.100 MPa和 393 K，0.380 MPa，與圖3中的實(shí)驗(yàn)工況十分接近，其最大蒸發(fā)速率分別為0.08 mg/s和0.042 mg/s。相比于振蕩環(huán)境下液滴蒸發(fā)速率，兩者相差兩個(gè)數(shù)量級(jí)。

此外，圖2還反映出，液滴蒸發(fā)速率隨著壓力振蕩的變化而呈現(xiàn)出相應(yīng)的變化趨勢(shì)，且蒸發(fā)速率波形相位滯后壓力振蕩波形接近180°，即燃燒室內(nèi)振蕩壓力處于波谷時(shí)，液滴蒸發(fā)速率接近波峰處，而當(dāng)振蕩壓力處于波峰時(shí)，液滴蒸發(fā)速率則接近波谷處。

圖2 壓力振蕩環(huán)境下液滴蒸發(fā)速率Fig.2 Droplet evaporation rate in pressure oscillation environment

圖3 壓力恒定環(huán)境下液滴蒸發(fā)速率Fig.3 Droplet evaporation rate in constant pressure environment

4 仿真結(jié)果及分析

通過(guò)實(shí)驗(yàn)雖然能夠?qū)毫φ袷幁h(huán)境下液滴蒸發(fā)過(guò)程進(jìn)行直觀地觀察，但由于實(shí)驗(yàn)條件的限制，液滴蒸發(fā)速率較低、環(huán)境因素影響較多，只能定性地得到結(jié)論，而無(wú)法定量獲得因壓力振蕩而造成的蒸發(fā)特征參數(shù)的變化。因此，基于前文壓力振蕩環(huán)境下液滴非穩(wěn)態(tài)蒸發(fā)模型，通過(guò)仿真獲得了不同振蕩頻率、振幅、燃燒室平均壓力及溫度條件下液滴蒸發(fā)過(guò)程特征參數(shù)，彌補(bǔ)了實(shí)驗(yàn)的不足，有助于進(jìn)一步了解燃燒不穩(wěn)定問(wèn)題。

4.1 振蕩頻率對(duì)液滴蒸發(fā)過(guò)程影響

圖4給出了不同振蕩頻率下液滴蒸發(fā)過(guò)程特征參數(shù)，其中，各算例除了振蕩頻率不同外，其余參數(shù)均保持一致。由圖4（a）可知，液滴蒸發(fā)速率會(huì)隨著燃燒室內(nèi)壓力的振蕩而產(chǎn)生相應(yīng)振蕩，這與實(shí)驗(yàn)結(jié)論相同。此外，振蕩頻率的增加加劇了蒸發(fā)速率的振蕩，但對(duì)其平均值并無(wú)影響，因而蒸發(fā)過(guò)程中液滴直徑的變化趨勢(shì)基本相同，如圖4（c）所示。分析認(rèn)為，振蕩壓力通過(guò)改變環(huán)境中氣體的密度來(lái)影響液滴蒸發(fā)速率，而較高的頻率使得液滴表面周?chē)鷼怏w密度變化加快，增加了燃料蒸氣進(jìn)入環(huán)境的速度，從而加劇了蒸發(fā)速率振蕩，但由于液滴處于相同的溫度、壓力環(huán)境中，從周?chē)h(huán)境吸收熱量的能力基本相同，因而蒸發(fā)速率的平均值相同。由圖4（b）可知，不同頻率下液滴表面溫度的變化趨勢(shì)基本相同。相較于100 Hz及以上的振蕩頻率，頻率較低時(shí)（5 Hz，13 Hz），液滴表面溫度振蕩更為明顯。這是因?yàn)榈皖l振蕩導(dǎo)致液滴吸收的熱量與蒸發(fā)實(shí)際需要的熱量不平衡，從而使得液滴表面溫度產(chǎn)生了較大波動(dòng)。綜合分析認(rèn)為，100 Hz以上的振蕩會(huì)影響液滴蒸發(fā)速率，而較低頻率振蕩（＜100 Hz）則會(huì)對(duì)液滴表面溫度及液滴直徑的變化趨勢(shì)產(chǎn)生影響。因此，研究壓力振蕩環(huán)境下的蒸發(fā)問(wèn)題要綜合考慮不同頻率振蕩所帶來(lái)的影響。

圖4 不同振蕩頻率下液滴蒸發(fā)特征參數(shù)Fig.4 Droplet evaporation characteristic parameters at different oscillation frequency

4.2 振蕩幅值對(duì)液滴蒸發(fā)過(guò)程影響

圖5給出了相同環(huán)境下不同振蕩幅值的液滴蒸發(fā)過(guò)程特征參數(shù)，算例中燃燒室平均壓力為0.1 MPa，振幅分別為燃燒室壓力的5%，10%和20%。由圖可知，壓力振蕩幅值越大，液滴蒸發(fā)速率振蕩越劇烈，液滴表面溫度振蕩越劇烈。分析認(rèn)為，較大振幅引發(fā)了液滴表面附近環(huán)境壓力的變化，從而對(duì)液滴蒸發(fā)過(guò)程產(chǎn)生了較大影響，使得蒸發(fā)速率和液滴表面溫度產(chǎn)生明顯振蕩；而較小振幅不足以使液滴表面環(huán)境壓力產(chǎn)生明顯變化，因而對(duì)液滴蒸發(fā)過(guò)程影響較小。

圖5 不同振蕩幅值下液滴蒸發(fā)特征參數(shù)Fig.5 Droplet evaporation characteristic parameters at different oscillation amplitude

4.3 燃燒室平均壓力對(duì)液滴蒸發(fā)過(guò)程影響

圖6給出了不同燃燒室平均壓力下液滴蒸發(fā)過(guò)程特征參數(shù)。由圖可知，增加燃燒室平均壓力大幅增加了液滴穩(wěn)態(tài)蒸發(fā)時(shí)的平衡溫度，液滴表面達(dá)到平衡溫度的時(shí)間延長(zhǎng)，如圖6（b），液滴直徑減小更為緩慢，如圖6（c）。此外，由圖6（a）可知，液滴表面未達(dá)到平衡溫度前，壓力越大，液滴蒸發(fā)速率越?。划?dāng)液滴表面達(dá)到平衡溫度后，三個(gè)工況下液滴蒸發(fā)速率基本相同。

圖6 不同燃燒室平均壓力下液滴蒸發(fā)特征參數(shù)Fig.6 Droplet evaporation characteristics parameters at different chamber average pressure

同時(shí)，增加燃燒室平均壓力，抑制了蒸發(fā)速率、液滴表面溫度等參數(shù)的振蕩。這是因?yàn)?，盡管算例增加了燃燒室平均壓力，但壓力振蕩幅值并沒(méi)有相應(yīng)增加，即振幅所占燃燒室平均壓力百分比減小，從而導(dǎo)致蒸發(fā)過(guò)程對(duì)振蕩的響應(yīng)程度減弱。因此，不僅燃燒室平均壓力影響液滴蒸發(fā)過(guò)程，振幅所占燃燒室平均壓力百分比對(duì)蒸發(fā)過(guò)程同樣具有重要作用。由此得出結(jié)論：在研究壓力振蕩環(huán)境下的蒸發(fā)問(wèn)題時(shí)，不能單純地研究燃燒室平均壓力或振蕩幅值，而是要將兩者結(jié)合起來(lái)共同分析，把振幅占燃燒室平均壓力的百分比作為影響因素之一來(lái)進(jìn)行考察。

由圖6可知，增加燃燒室平均壓力大幅增加了液滴穩(wěn)態(tài)蒸發(fā)時(shí)的平衡溫度，液滴表面達(dá)到平衡溫度的時(shí)間延長(zhǎng)，如圖6（b），液滴直徑減小更為緩慢，如圖6（c）。此外，由圖6（a）可知，液滴表面未達(dá)到平衡溫度前，壓力越大，液滴蒸發(fā)速率越小；當(dāng)液滴表面達(dá)到平衡溫度后，三個(gè)工況下液滴蒸發(fā)速率基本相同。同時(shí)，增加燃燒室平均壓力，抑制了蒸發(fā)速率、液滴表面溫度等參數(shù)的振蕩。這是因?yàn)?，盡管算例增加了燃燒室平均壓力，但壓力振蕩幅值并沒(méi)有相應(yīng)增加，即振幅所占燃燒室平均壓力百分比減小，從而導(dǎo)致蒸發(fā)過(guò)程對(duì)振蕩的響應(yīng)程度減弱。因此，不僅燃燒室平均壓力影響液滴蒸發(fā)過(guò)程，振幅所占燃燒室平均壓力百分比對(duì)蒸發(fā)過(guò)程同樣具有重要作用。由此得出結(jié)論：在研究壓力振蕩環(huán)境下的蒸發(fā)問(wèn)題時(shí)，不能單純地研究燃燒室平均壓力或振蕩幅值，而是要將兩者結(jié)合起來(lái)共同分析，把振幅占燃燒室平均壓力的百分比作為影響因素之一來(lái)進(jìn)行考察。

4.4 燃燒室溫度對(duì)液滴蒸發(fā)過(guò)程影響

圖7給出了不同燃燒室溫度下液滴蒸發(fā)過(guò)程特征參數(shù)。由圖7可知，高溫環(huán)境在加快了液滴蒸發(fā)速率的同時(shí)，也加劇了蒸發(fā)速率的振蕩。為了便于觀察，圖7（a）顯示了蒸發(fā)時(shí)間0～3 s內(nèi)蒸發(fā)速率的變化。當(dāng)燃燒室溫度較大時(shí)（＞1 000 K），液滴穩(wěn)態(tài)蒸發(fā)時(shí)的平衡溫度較高，促進(jìn)了蒸發(fā)速率大幅增加，且液滴表面在較短時(shí)間內(nèi)迅速達(dá)到平衡溫度，使得液滴蒸發(fā)速率迅速上升至最大值；由于快速蒸發(fā)加快了燃料蒸氣的產(chǎn)生，使得液滴表面燃料蒸氣質(zhì)量分?jǐn)?shù)迅速增加，從而抑制了液滴蒸發(fā)，造成蒸發(fā)速率的減小，因而高溫下液滴蒸發(fā)速率呈現(xiàn)出先增加后減小的趨勢(shì)；而燃燒室溫度較低時(shí)（＜1 000 K），增加溫度對(duì)液滴蒸發(fā)速率的促進(jìn)作用較小，因而液滴蒸發(fā)速率并未呈現(xiàn)出先增加后減小的趨勢(shì)。同時(shí)，較高的溫度環(huán)境增大了液滴表面與周?chē)h(huán)境的溫度梯度，使得環(huán)境的傳熱能力增強(qiáng)，因而液滴表面能夠更快地達(dá)到平衡溫度，如圖7（b）所示；液滴蒸發(fā)更為迅速，直徑減小得更快，如圖7（c）所示。分析認(rèn)為，蒸發(fā)的加快使得液滴表面燃料蒸氣增多，而環(huán)境壓力的振蕩造成燃料蒸氣不斷進(jìn)入到氣相環(huán)境中，液滴蒸發(fā)又不斷地產(chǎn)生新的燃料蒸氣。溫度越高，燃料蒸氣進(jìn)入環(huán)境的速度就越快，蒸發(fā)產(chǎn)生燃料蒸氣的速度也越快，液滴表面燃料蒸氣質(zhì)量分?jǐn)?shù)處于不斷變化中，從而使得液滴蒸發(fā)速率產(chǎn)生大幅振蕩。

圖7 不同燃燒室溫度下液滴蒸發(fā)特征參數(shù)Fig.7 Droplet evaporation characteristics parameters at different chamber temperature

5 結(jié)論

利用壓力振蕩環(huán)境實(shí)驗(yàn)系統(tǒng)，開(kāi)展了壓力振蕩環(huán)境下液滴蒸發(fā)實(shí)驗(yàn)研究，分析了壓力振蕩對(duì)液滴蒸發(fā)過(guò)程的影響，建立壓力振蕩環(huán)境液滴非穩(wěn)態(tài)蒸發(fā)模型，并基于該模型考察了振蕩頻率、振幅，燃燒室平均壓力及溫度等因素對(duì)液滴蒸發(fā)特性的影響規(guī)律，得出了以下主要結(jié)論：

1）與壓力恒定環(huán)境相比，壓力振蕩環(huán)境對(duì)液滴蒸發(fā)過(guò)程具有明顯的促進(jìn)作用，極大地增加了液滴蒸發(fā)速率，縮短了液滴生存時(shí)間；此外，液滴蒸發(fā)速率隨著壓力振蕩的變化呈現(xiàn)出相應(yīng)的變化趨勢(shì)，且蒸發(fā)速率波形相位滯后壓力振蕩波形接近 180°。

2）增加振蕩頻率會(huì)加劇蒸發(fā)速率的振蕩，但對(duì)蒸發(fā)速率平均值并無(wú)影響，不同頻率下的液滴直徑變化趨勢(shì)并無(wú)明顯區(qū)別。

3）在研究壓力振蕩環(huán)境下的蒸發(fā)問(wèn)題時(shí)，不能單純地考察燃燒室平均壓力或者振蕩幅值的影響作用，而要將兩者結(jié)合起來(lái)共同分析，把振幅占燃燒室平均壓力的百分比作為影響因素之一來(lái)進(jìn)行考察，因?yàn)榘俜直仍酱螅旱握舭l(fā)速率及表面溫度振蕩越劇烈。

4）燃燒室的高溫環(huán)境在增加液滴蒸發(fā)速率，促進(jìn)蒸發(fā)的同時(shí)，也會(huì)增強(qiáng)壓力振蕩對(duì)蒸發(fā)過(guò)程的影響效應(yīng)，加劇液滴蒸發(fā)速率振蕩。因此，液體火箭發(fā)動(dòng)機(jī)燃燒室中的超高溫環(huán)境會(huì)加劇不穩(wěn)定燃燒。

[1]程心.液體火箭發(fā)動(dòng)機(jī)燃燒不穩(wěn)定性[M].北京：國(guó)防工業(yè)出版社,1965.

[2]蘇凌宇,聶萬(wàn)勝,劉衛(wèi)東.低頻壓力振蕩環(huán)境下運(yùn)動(dòng)液滴蒸發(fā)過(guò)程的準(zhǔn)穩(wěn)態(tài)模型[J].推進(jìn)技術(shù),2011,32(5)：670-675.

[3]YANG V.Liquid propellant rocket engine combustion instability[M].USA：American Institute of Aeronautics&Astronautics,1995.

[4]STRAHLE W C.A theoretical study of unsteady droplet burning：transients and periodic solutions[D].Princeton,NJ：Department of Aeronautical Engineering,Princeton University,1963.

[5]SHEU W J,LIOU N C.Effect of temporal variation of pressureonvaporizationofliquiddroplets[J].International Journal of Heat and Mass Transfer.1999(42)：4043-4054.

[6]YANG Vigor,MICHAEL M M.Liquid rocket motor combustion stability using coaxial injectors and supercritical droplet combustion and dynamics,SA-SA GF49620-93-1-0126[R].[S.l.]：[s.n.],1995.

[7]OKAI K,MORIUE O,ARAKI M,et al.Combustion of single droplets and droplet pairs in a vibrating field under microgravity[J].Proceeding of the Combustion Institute,2000,28：977-983.

[8]TANABLE M,MORITA T,AOKI K,et al.Influence of standing sound waves on droplet combustion[J].Proceeding of the Combustion Institute,2000,28：1007-1013.

[9]TANABLE M,KUWAHARA T,SATOH K,et al.Droplet combustion in standing sound waves[J].Proceeding of the Combustion Institute,2005(30)：1957-1964.