張留歡，杜 泉，逯婉若，張蒙正，李光熙

(西安航天動(dòng)力研究所，陜西 西安710100)

0 引言

空氣渦輪火箭發(fā)動(dòng)機(jī)(Air Turbo Rocket, ATR)將火箭發(fā)動(dòng)機(jī)與渦輪噴氣發(fā)動(dòng)機(jī)組合[1-2]，使渦輪入口氣體參數(shù)與飛行條件解耦，拓寬了發(fā)動(dòng)機(jī)工作范圍，其可作為臨近空間飛行器的新型動(dòng)力。

從熱力循環(huán)方式來(lái)說(shuō)，ATR發(fā)動(dòng)機(jī)主要分為燃?xì)獍l(fā)生器循環(huán)和膨脹循環(huán)兩種方案[3-7]。燃?xì)獍l(fā)生器循環(huán)ATR發(fā)動(dòng)機(jī)多使用雙組元(液氧/液氫、液氧/甲烷及H2O2/煤油等)作為推進(jìn)劑[8-9]，這種方式發(fā)動(dòng)機(jī)自帶氧化劑，且推進(jìn)劑低溫或熱穩(wěn)定性差；膨脹循環(huán)ATR發(fā)動(dòng)機(jī)(Air-Turbo-Ram Expander Cycle，ATREX)使用液氫作為推進(jìn)劑[10-11]，其做功能力強(qiáng)，燃燒熱值高，但密度低，儲(chǔ)存容積大且對(duì)材料性能要求高，操作性較差[12]。

甲醇(CH3OH)作為燃料在汽車發(fā)動(dòng)機(jī)等內(nèi)燃機(jī)方面的應(yīng)用研究較多[13-15]，其成分單一，特性穩(wěn)定，屬于清潔燃料。一般利用發(fā)動(dòng)機(jī)排氣余熱將常溫液體甲醇加熱、蒸發(fā)[16]，并在催化劑作用下，將甲醇裂解[17-19]為H2和CO等，進(jìn)入發(fā)動(dòng)機(jī)內(nèi)燃燒。甲醇裂解氣中氫氣含量最高可達(dá)67%，可充分利用氫氣的燃燒特性，但不需要解決儲(chǔ)氫的問(wèn)題[20]。同時(shí)，裂解氣中碳?xì)浠衔锖繕O少，結(jié)焦問(wèn)題不明顯。此外，張琴等[21]研究了甲醇作為吸熱燃料應(yīng)用于超燃沖壓發(fā)動(dòng)機(jī)主動(dòng)冷卻技術(shù)的可行性，給出了不同涂層催化劑對(duì)甲醇分解的影響結(jié)果。

本文提出了一種基于甲醇裂解工作的ATR發(fā)動(dòng)機(jī)(下文簡(jiǎn)稱甲醇ATR發(fā)動(dòng)機(jī))方案，詳細(xì)闡述了該發(fā)動(dòng)機(jī)工作原理，并開(kāi)展了發(fā)動(dòng)機(jī)性能計(jì)算，研究了甲醇裂解度對(duì)發(fā)動(dòng)機(jī)性能的影響規(guī)律，可為后續(xù)吸氣式組合發(fā)動(dòng)機(jī)方案研究及推進(jìn)劑選取提供新的思路。

1 工作原理及性能

甲醇ATR發(fā)動(dòng)機(jī)主要包括進(jìn)氣道、壓氣機(jī)、氣體發(fā)生裝置、渦輪、燃燒室、尾噴管及換熱裝置等組件。甲醇ATR發(fā)動(dòng)機(jī)的特點(diǎn)主要是其氣體發(fā)生裝置利用甲醇受熱、催化裂解產(chǎn)生H2和CO等氣體的特性。甲醇裂解的反應(yīng)方程為：CH3OH+90kJ→2H2+CO，為吸熱反應(yīng)。圖1為甲醇ATR發(fā)動(dòng)機(jī)工作原理示意圖。在發(fā)動(dòng)機(jī)工作過(guò)程中，首先采用單獨(dú)設(shè)計(jì)安裝的燃?xì)獍l(fā)生器(固體或液體)產(chǎn)生高溫高壓富燃燃?xì)怛?qū)動(dòng)發(fā)動(dòng)機(jī)短暫工作，之后一定壓力、流量的甲醇進(jìn)入發(fā)動(dòng)機(jī)，通過(guò)換熱裝置實(shí)現(xiàn)升溫裂解，產(chǎn)生的高溫高壓氣相物進(jìn)入氣體發(fā)生裝置，直接驅(qū)動(dòng)渦輪帶動(dòng)壓氣機(jī)工作，大氣中的空氣經(jīng)壓氣機(jī)增壓后直接進(jìn)入渦輪后的燃燒室，在燃燒室內(nèi)和經(jīng)過(guò)渦輪做功后的可燃?xì)怏w進(jìn)行燃燒，生成的高溫燃?xì)馔ㄟ^(guò)噴管膨脹產(chǎn)生推力。

為方便公式說(shuō)明，圖2給出了ATR發(fā)動(dòng)機(jī)組成及截面示意圖?？紤]發(fā)動(dòng)機(jī)遠(yuǎn)前端來(lái)流截面(序號(hào)0)及進(jìn)氣道入口截面(序號(hào)1)，給定發(fā)動(dòng)機(jī)入口(即壓氣機(jī)入口)截面序號(hào)為2，發(fā)動(dòng)機(jī)出口(即尾噴管出口)截面序號(hào)為8。

圖1 甲醇ATR發(fā)動(dòng)機(jī)工作原理示意圖Fig.1 Sketch of working principle for methanol ATR engine

圖2 ATR發(fā)動(dòng)機(jī)及其截面示意圖Fig.2 Sketch of ATR engine and its sections

根據(jù)甲醇ATR發(fā)動(dòng)機(jī)工作原理，建立了發(fā)動(dòng)機(jī)工作平衡方程。

1)渦輪與壓氣機(jī)物理轉(zhuǎn)速平衡:

nc=nt

(1)

式中：n為物理轉(zhuǎn)速；下標(biāo)c為壓氣機(jī)；下標(biāo)t為渦輪。

2)發(fā)動(dòng)機(jī)流道壓力平衡(不考慮氣體沿程壓力損失)：

(2)

式中：p為發(fā)動(dòng)機(jī)截面氣體總壓；π為壓比；下標(biāo)2，4和7分別為壓氣機(jī)入口截面、渦輪入口截面和噴管喉部截面。

3)氣體質(zhì)量流量平衡：

2006年5月20日，湖南省鳳凰縣、貴州省雷山縣由于其獨(dú)特的苗族銀飾鍛造技術(shù)，順利入選國(guó)務(wù)院第一批國(guó)家級(jí)非物質(zhì)文化遺產(chǎn)代表名錄。

m7=m2+m4

(3)

式中m為單位時(shí)間通過(guò)某截面氣體質(zhì)量流量。

4)渦輪與壓氣機(jī)功率平衡(不考慮機(jī)械功損失)：

Lc=Lt

(4)

式中L為功率。

基于上述基本平衡方程，建立了發(fā)動(dòng)機(jī)計(jì)算模型，發(fā)動(dòng)機(jī)設(shè)計(jì)點(diǎn)(100%轉(zhuǎn)速百分比)部分組件參數(shù)取值見(jiàn)表1。發(fā)動(dòng)機(jī)設(shè)計(jì)點(diǎn)空氣流量為10 kg/s，同時(shí)考慮到換熱裝置設(shè)計(jì)難度，給定甲醇裂解氣總溫(渦輪前氣體總溫)為1 000 K。

表1 甲醇ATR發(fā)動(dòng)機(jī)設(shè)計(jì)點(diǎn)組件參數(shù)(部分)

基于上述發(fā)動(dòng)機(jī)參數(shù)，對(duì)地面狀態(tài)甲醇ATR發(fā)動(dòng)機(jī)性能進(jìn)行了計(jì)算。計(jì)算過(guò)程中，假設(shè)甲醇完全裂解為H2和CO，則進(jìn)入渦輪的甲醇裂解氣中H2和CO的質(zhì)量比為1∶7。圖3和圖4分別給出了完全裂解條件下甲醇ATR發(fā)動(dòng)機(jī)推力F、比沖Isp和燃燒室參數(shù)(總壓p、總溫T、余氣系數(shù)α)隨發(fā)動(dòng)機(jī)轉(zhuǎn)速n變化結(jié)果(橫坐標(biāo)為轉(zhuǎn)速百分比)。燃燒室參數(shù)在較大程度上影響著發(fā)動(dòng)機(jī)推力性能。圖3顯示，隨著發(fā)動(dòng)機(jī)轉(zhuǎn)速上升，發(fā)動(dòng)機(jī)推力逐漸上升，比沖先增大后減小。在90%轉(zhuǎn)速時(shí)，發(fā)動(dòng)機(jī)比沖最高約793 s(介于火箭發(fā)動(dòng)機(jī)與渦輪發(fā)動(dòng)機(jī)之間)。甲醇裂解氣H2和CO在渦輪中做功后進(jìn)入燃燒室與空氣摻混燃燒，隨著轉(zhuǎn)速上升，燃燒室余氣系數(shù)逐漸減小，燃燒總溫與總壓逐漸升高，具體如圖4所示。

圖3 推力、比沖隨轉(zhuǎn)速變化結(jié)果Fig.3 Results of thrust and specific impulse with rotating speed

圖4 燃燒室參數(shù)隨轉(zhuǎn)速變化結(jié)果Fig.4 Results of combustor parameters with rotating speed

2 裂解度影響

甲醇裂解存在不完全現(xiàn)象，此時(shí)裂解氣主要含有H2和CO以及未裂解的甲醇?xì)怏w[22]。采用特殊催化劑(如銅基/鈀基催化劑[21-22])可提高甲醇裂解度N(甲醇裂解度定義為已裂解甲醇與全部甲醇質(zhì)量比)。圖5和圖6分別給出了基于上述發(fā)動(dòng)機(jī)組件參數(shù)的甲醇ATR發(fā)動(dòng)機(jī)甲醇流量mf、推力及比沖等隨裂解度N(75%,80%，85%，90%，95%及100%)、轉(zhuǎn)速(60%，70%，80%，90%及100%)等的變化結(jié)果。其中，假設(shè)已裂解的甲醇全部反應(yīng)為H2和CO氣體。

圖5 不同裂解度、不同轉(zhuǎn)速下甲醇流量變化結(jié)果Fig.5 Mass flow rate of methanol at different cracking ratios and rotating speeds

圖6 不同裂解度、不同轉(zhuǎn)速下推力和比沖變化結(jié)果Fig.6 Thrust and specific impulse at different cracking ratios and rotating speeds

圖7給出了1 000 K溫度條件下甲醇裂解氣比熱容cp[23]隨裂解度變化的結(jié)果。圖7顯示，隨著裂解度逐漸升高，裂解氣比熱容逐漸升高。由于渦輪等熵膨脹功有LT=mf·cp·ΔT[24]，在相同的膨脹功條件下，渦輪前裂解氣的cp越高，需要供應(yīng)的甲醇流量mf越小，發(fā)動(dòng)機(jī)比沖越高，即同一轉(zhuǎn)速下，甲醇裂解度越高，裂解氣cp數(shù)值越大，甲醇流量越小，發(fā)動(dòng)機(jī)比沖越高。

圖7 不同裂解度時(shí)裂解氣cp變化結(jié)果Fig.7 Specific heat capacity of cracking gas at different cracking ratios

3 結(jié)論

提出了一種基于甲醇裂解工作的ATR發(fā)動(dòng)機(jī)方案，并對(duì)裂解度100%的甲醇ATR發(fā)動(dòng)機(jī)性能進(jìn)行計(jì)算，研究了裂解度對(duì)發(fā)動(dòng)機(jī)性能的影響規(guī)律，獲得以下結(jié)論：

1)隨著發(fā)動(dòng)機(jī)轉(zhuǎn)速上升，推力逐漸上升，比沖先增加后減小。在90%轉(zhuǎn)速百分比時(shí)，比沖最高。

2)在同一轉(zhuǎn)速下，隨著裂解度升高，比沖逐漸升高。在裂解度100%、轉(zhuǎn)速百分比90%條件下，發(fā)動(dòng)機(jī)比沖最高達(dá)到約793 s。

3)同一轉(zhuǎn)速下，甲醇裂解度越高，裂解氣比熱容數(shù)值越大。