• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    預(yù)冷型組合循環(huán)發(fā)動機技術(shù)

    2022-03-16 05:31:20羅佳茂楊順華母忠強張彎洲
    空氣動力學(xué)學(xué)報 2022年1期
    關(guān)鍵詞:冷卻劑馬赫數(shù)預(yù)冷

    羅佳茂,楊順華,母忠強,張彎洲,游 進(jìn)

    (北京流體動力科學(xué)研究中心,北京 100011)

    0 引 言

    隨著航空技術(shù)的不斷進(jìn)步,人類開始追求越來越高的飛行馬赫數(shù)。作為高速飛行器的核心部件,全速域發(fā)動機在飛行器實現(xiàn)高速飛行中發(fā)揮著決定性作用[1]。雖然單一類型發(fā)動機(渦輪發(fā)動機、沖壓發(fā)動機和火箭發(fā)動機等)技術(shù)相對成熟,但均無法滿足全速域且可重復(fù)使用的要求,如:常規(guī)渦輪發(fā)動機最高飛行速度不超過馬赫數(shù)2.5,飛行高度30 km以下;超燃沖壓發(fā)動機需在馬赫數(shù)3.0以上才能起動,且機動性較差;火箭發(fā)動機能在全速域內(nèi)工作,亦不受高度限制,但比沖低,也難以重復(fù)使用。

    為解決全速域飛行器動力問題,各研究機構(gòu)結(jié)合單一發(fā)動機工作特點提出了多種類型的組合循環(huán)發(fā)動機,常見類型包括:渦輪基組合循環(huán)發(fā)動機(Turbine Based Combined Cycle,TBCC)、火箭基組合循環(huán)發(fā)動機(Rocket Based Combined Cycle,RBCC)和空氣渦輪火箭發(fā)動機(Air Turbo-Rocket,ATR)等[2-3]。組合循環(huán)發(fā)動機在結(jié)構(gòu)布局和工作區(qū)間上對單一類型發(fā)動機進(jìn)行了有機融合,具有工作速域?qū)?、比沖性能較好、可重復(fù)使用和經(jīng)濟性好等優(yōu)點,在未來單級入軌(Single Stage to Orbit,SSTO)和兩級入軌(Two Stage to Orbit,TSTO)可重復(fù)使用飛行器中有廣闊的應(yīng)用前景。

    若低速段采用渦輪類發(fā)動機作為主動力,而高速段采用超燃沖壓發(fā)動機提供主動力,則組合循環(huán)發(fā)動機會面臨工作速域和推力的銜接問題,需采取相關(guān)措施擴展渦輪類發(fā)動機或超燃沖壓發(fā)動機工作包線。

    預(yù)冷型組合循環(huán)發(fā)動機是利用低溫介質(zhì)對來流高溫空氣進(jìn)行預(yù)冷或液化的一種新型組合循環(huán)發(fā)動機,具有工作速域?qū)?、比沖高和推重比大等優(yōu)點,是組合發(fā)動機的重要發(fā)展方向,在未來空天領(lǐng)域應(yīng)用潛力大。美國、日本、英國和俄羅斯等國已對預(yù)冷型組合循環(huán)發(fā)動機開展了諸多研究工作,取得了大量理論和試驗結(jié)果[4-5]。目前預(yù)冷型組合循環(huán)發(fā)動機面臨的主要問題集中在輕質(zhì)高效的預(yù)冷器設(shè)計與加工制造技術(shù)、開式/閉式循環(huán)系統(tǒng)設(shè)計技術(shù)、全速域循環(huán)參數(shù)匹配優(yōu)化技術(shù)、高效增壓技術(shù)和大功率渦輪功提取技術(shù)。

    本文首先對國內(nèi)外幾類典型的預(yù)冷型組合循環(huán)發(fā)動機研究進(jìn)展進(jìn)行了詳細(xì)介紹,然后通過數(shù)值建模計算分析了預(yù)冷型組合循環(huán)發(fā)動機熱力循環(huán)中的幾個關(guān)鍵難點,為預(yù)冷型組合循環(huán)發(fā)動機的熱力循環(huán)設(shè)計和部件設(shè)計提供參考。

    1 預(yù)冷型組合循環(huán)發(fā)動機發(fā)展歷程

    預(yù)冷型組合循環(huán)發(fā)動機發(fā)展大致經(jīng)歷了三個階段(圖1):20世紀(jì)六七十年代,主要以LACE循環(huán)為代表的液化空氣循環(huán)發(fā)動機;20世紀(jì)八九十年代,以RB545為代表的深度預(yù)冷循環(huán)發(fā)動機;20世紀(jì)90年代中期到現(xiàn)在,以SABRE、ATREX為代表的適度預(yù)冷循環(huán)發(fā)動機。

    圖1 預(yù)冷變循環(huán)組合發(fā)動機研制歷程Fig. 1 Research progress of pre-cooled variable-cycle combined engine

    1.1 液化空氣循環(huán)——LACE與ACES發(fā)動機

    20世紀(jì)60年代,Marquardt公司提出了液化空氣循環(huán)發(fā)動機(Liquid Air Cycle Engines,LACE)[6]。該發(fā)動機融合了吸氣式發(fā)動機和火箭發(fā)動機工作模式,在液氫液氧火箭發(fā)動機基礎(chǔ)上研發(fā)而來,最高能工作至馬赫數(shù)6~7。基本原理如圖2所示,在吸氣模態(tài),發(fā)動機利用低溫高熱沉的液氫作為冷卻劑將捕獲的空氣液化,再將液態(tài)空氣增壓后噴入傳統(tǒng)火箭發(fā)動機燃燒室。工作模態(tài)轉(zhuǎn)換為火箭模態(tài)后,利用攜帶的液氧代替液化空氣作為氧化劑[7]。

    圖2 LACE發(fā)動機原理圖[7]Fig. 2 Principle diagram of the LACE cycle[7]

    LACE發(fā)動機的主要優(yōu)點在于吸氣模態(tài)利用空氣中的氧氣作為氧化劑,減少了氧化劑的攜帶,同時吸氣模態(tài)和火箭模態(tài)共用噴管,不僅減輕了重量,也讓空間利用更緊湊。不足之處在于與其他吸氣式發(fā)動機相比,由于冷凝空氣需要消耗大量冷卻劑氫,導(dǎo)致比沖較低,僅為800 s左右。地面狀態(tài),LACE發(fā)動機所需的冷卻劑氫/空氣質(zhì)量比約0.2~0.35,數(shù)倍于燃燒化學(xué)恰當(dāng)比,來流空氣溫度升高后,所需冷卻劑氫量更大,這制約了吸氣模態(tài)發(fā)動機比沖性能[6]。

    1993至1996年,俄羅斯Keldysh研究中心對LACE發(fā)動機開展了一系列計算、優(yōu)化設(shè)計和試驗,并提出了地面推力500 kN的發(fā)動機方案(圖3)[8]和推力10 kN的燃燒室(圖4)[9],結(jié)果證明:熱交換器效率能達(dá)到2000~3000 W/(m2·K);獲得了5~8 K的超冷液態(tài)空氣;壓氣機工作溫度100~120 K;液態(tài)空氣/液氫和液氧/液氫雙模態(tài)燃燒室能實現(xiàn)點火工作[8-10]。

    圖3 500 kN地面推力LACE發(fā)動機布局[8]Fig. 3 Layout of LACE engine with a sea level thrust of 500 kN[8]

    圖4 10 kN推力LACE發(fā)動機燃燒室結(jié)構(gòu)[9]Fig. 4 Chamber structure of the LACE engine with a thrust of 10 kN[9]

    LACE發(fā)動機難點在于液化空氣中的氧氣和氮氣所需的溫度在100 K以下,此溫度已低于二氧化碳和水的凝固點,固態(tài)二氧化碳與冰會堵塞冷凝器,且熱交換器表面霜凍問題無良好的解決辦法。攜帶大量氫也將占用飛行器較大空間。另外,能否實現(xiàn)穩(wěn)定的模態(tài)轉(zhuǎn)換暫無試驗結(jié)果支撐[6,8]。

    ACES( Air Collection and Enrichment System) 發(fā)動機是在LACE發(fā)動機基礎(chǔ)上增加了液氧分離器而發(fā)展出來的衍生型。從20世紀(jì)60年代至21世紀(jì)初,多個機構(gòu)對其進(jìn)行了改進(jìn)研制,比較有代表性的為安德魯斯航天公司提出的Alchemist ACES發(fā)動機,其基本工作原理如圖5所示[11]。ACES發(fā)動機同樣需要用低溫冷卻劑氫將來流空氣液化,但液化后還需將液氧與液氮進(jìn)行分離,部分液氧用于供應(yīng)吸氣模態(tài)燃燒室,剩余部分注入液氧罐為之后的火箭模態(tài)備用。

    圖5 ACES發(fā)動機原理圖[11]Fig. 5 Principle diagram of the ACES engine[11]

    相比于LACE發(fā)動機,ACES發(fā)動機由于加入了液氧分離器而變得更復(fù)雜。該部件的設(shè)計旨在減少飛行器的起飛重量,因為起飛狀態(tài)液氧罐處于空瓶狀態(tài)。實際上,液化空氣需消耗大量液氫,所需液氫與空氣質(zhì)量比接近0.2,因此起飛狀態(tài)需攜帶更多液氫,并未減少飛行器起飛重量[6]。而且液氫密度低,會占用飛行器較大體積,給飛行器空間布局的設(shè)計也造成了困難[11-12]。

    ACES發(fā)動機的關(guān)鍵技術(shù)有兩點[13]:一是液氧分離器,因為液氧與液氮物理屬性相近,分離器的設(shè)計將變得更復(fù)雜困難;二是輕質(zhì)高效的氫預(yù)冷器。Hendrick等[14]通過試驗驗證了氫預(yù)冷器的可行性,并證明了預(yù)冷器幾乎無泄漏(僅泄漏0.1 ug/s的氫),綜合性能較高,之后還發(fā)展了多種類型的管殼式鋁材質(zhì)預(yù)冷器[14]。

    空氣液化循環(huán)發(fā)動機對冷卻劑的消耗量極大,發(fā)動機系統(tǒng)復(fù)雜,技術(shù)難點多,經(jīng)濟性和實用性均較差。

    1.2 空氣深冷循環(huán)——RB545和ATRDC發(fā)動機

    20世紀(jì)80年代,英國開始發(fā)展可重復(fù)使用航天飛機技術(shù)。英國BAE系統(tǒng)公司提出了水平起降、單級入軌空天飛機(HOTOL)。HOTOL飛機采用羅?羅公司Alan Band等提出的預(yù)冷吸氣式火箭發(fā)動機方案——RB545,其基本原理如圖6所示[15]。

    圖6 RB545發(fā)動機基本原理[15]Fig. 6 Principle diagram of the RB545 engine[15]

    與LACE發(fā)動機類似,RB545發(fā)動機有兩種工作模式—吸氣模式和純火箭模式。不同之處在于,RB545在吸氣模式僅對空氣進(jìn)行深度預(yù)冷,但為了避免空氣相變過程的能量轉(zhuǎn)換,最低只將空氣冷卻至約80 K,并不使空氣液化,這種改進(jìn)節(jié)省了冷卻劑用量。吸氣模 式 空 氣 流 量240 kg/s,燃 燒 室 溫 度2 700 K,推 力340 kN。火箭模式采用自身攜帶的液氫液氧工作,推力室室壓15 MPa,氫流量23.3 kg/s,燃燒室溫度3 400 K,推力735 kN[4,15]。

    為簡化發(fā)動機結(jié)構(gòu),兩種工作模式共用推力室,因此吸氣模式下也需達(dá)到較高的推力室壓力,約15 MPa,這就要求壓氣機有較高的壓比。為滿足壓氣機高壓比的功率需求,大約2/3的高焓高壓氫將用于驅(qū)動渦輪。膨脹后的氫已處于低壓狀態(tài),不滿足推力室壓力要求,只能排放掉,少部分膨脹后的氫進(jìn)入沖壓燃燒室燃燒掉。大量氫的排放也導(dǎo)致RB545比沖僅有1 508 s,混合比10.4[3]。

    蘇聯(lián)在RB545的基礎(chǔ)上發(fā)展出了深冷空氣渦輪火箭發(fā)動機(Deeply Cooled Air Turborocket, ATRDC),其主要部件包括進(jìn)氣道、預(yù)冷器、壓氣機、燃燒室和氫加熱器,其中燃燒室分為內(nèi)外兩管路,內(nèi)管以吸氣模式運行,外管以火箭模式運行,原理如圖7所示,工作速域馬赫數(shù)0~6。ATRDC發(fā)動機同樣采用液氫為冷卻劑,首先對來流空氣進(jìn)行深度冷卻,再經(jīng)氫加熱器進(jìn)一步吸熱后達(dá)到1 000~1 200 K,之后高焓高壓狀態(tài)氫驅(qū)動渦輪機工作,過剩的氫將由單獨噴管排出[1]。

    圖7 ATRDC發(fā)動機原理圖[16]Fig. 7 Principle diagram of the ATRDC engine[16]

    由于預(yù)冷后空氣溫度低至98~112 K,ATRDC能擁有最高40的壓比,工作速域內(nèi)比沖3 300~1 500 s,推重比預(yù)估高達(dá)18~22[16]。該型發(fā)動機的缺點在于預(yù)冷器重量過大,約占整個發(fā)動機(不含進(jìn)氣道)的40%[16],且預(yù)冷對冷卻劑消耗量大,一半的氫在膨脹后將被直接排放掉。若將ATRDC發(fā)動機與沖壓發(fā)動機進(jìn)行組合,將過剩的氫注入沖壓發(fā)動機燃燒,將有效提升發(fā)動機推力和比沖性能。由圖8可以看出,充分集成后的ATRDC與沖壓組合發(fā)動機最高比沖能提升至4 000 s以上[16]。

    圖8 ATRDC及其組合發(fā)動機比沖分布[16]Fig. 8 Specific impulse distribution of ATRDC and its combined engines[16]

    空氣深冷循環(huán)發(fā)動機性能高,但仍然面臨冷卻劑過剩問題,在經(jīng)濟性和實用性上優(yōu)勢并不算突出。

    1.3 適度預(yù)冷閉式膨脹循環(huán)——SABRE和Scimitar發(fā)動機

    1988年,由于RB545發(fā)動機技術(shù)風(fēng)險高、所需資金也過大,羅?羅公司停止了對該型發(fā)動機的投資,HOTOL項目被迫終止。但Alan Band并未灰心,1989年與HOTOL項目其余數(shù)名核心成員(Richard Varvil、John Scott等)創(chuàng)立了反應(yīng)引擎公司(Reaction Engines Limited,REL),繼續(xù)研發(fā)空天飛行器及其動力系統(tǒng)。經(jīng)過5年努力,于1994年提出了“云霄塔”(SKYLON)單級入軌空天飛機概念,同時在LACE和RB545的基礎(chǔ)上通過增設(shè)氦閉式循環(huán)而形成新的發(fā)動機方案——協(xié)同吸氣式火箭發(fā)動機(Synergistic Air-Breathing Rocket Engine,SABRE),即“佩刀”發(fā)動機[17],基本結(jié)構(gòu)如圖9所示。通過引入第三介質(zhì)(氦)閉式循環(huán),一是避免了“氫脆”問題,二是氦閉式布萊頓循環(huán)將部分熱量用于驅(qū)動空氣壓氣機和氦渦輪,分流了冷卻劑(氫)需承受的熱量,可有效減少冷卻劑用量[18]。

    圖9 SABRE發(fā)動機結(jié)構(gòu)示意圖[17]Fig. 9 Baseline of SABRE engine[17]

    截至2013年,Alan Bond團(tuán)隊已掌握了SABRE-3發(fā)動機多項關(guān)鍵技術(shù),尤其是輕質(zhì)高效預(yù)冷器的設(shè)計和加工技術(shù)[19],2003年與布里斯托大學(xué)合作突破了預(yù)冷器霜控技術(shù),開展了大量預(yù)冷器試驗并證明了設(shè)計的預(yù)冷器能滿足SABRE飛行需求[20],為之后的進(jìn)一步改進(jìn)提供了技術(shù)支撐[21]。表1給出了SABRE-3部分性能參數(shù)。熱力循環(huán)原理如圖10所示。

    圖10 SABRE-3發(fā)動機熱力循環(huán)原理圖[18]Fig. 10 Principle diagram of the SABRE-3 engine thermal cycle[18]

    表1 SABRE-3發(fā)動機關(guān)鍵運行參數(shù)[21]Table 1 Key operating parameters for SABRE-3 engine[21]

    與RB545類似,早期的SABRE-3發(fā)動機吸氣模態(tài)和火箭模態(tài)共用燃燒室和噴管,這就要求吸氣模態(tài)空氣需達(dá)到較高壓力。為提升壓比,空氣將被深度冷卻至140 K左右,壓比更是達(dá)到140,為滿足壓氣機高壓比的需求,需消耗大量液氫用于預(yù)冷和驅(qū)動壓氣機工作,總需求量達(dá)到2.8倍當(dāng)量比,提高比沖成了下一步工作的關(guān)鍵。2015年,REL公司公開了SABRE-4發(fā)動機系統(tǒng),為減少液氫用量,SABRE-4采用雙喉道噴管解耦了吸氣模態(tài)和火箭模態(tài)推力室,盡管結(jié)構(gòu)重量有所增加,但降低了冷卻劑用量,當(dāng)量比降低至1.2左右。熱力循環(huán)原理如圖11所示。冷卻劑用量的降低,也使壓氣機功率降低,壓比降低至14~30,因此允許更高的預(yù)冷器出口溫度,空氣預(yù)冷后溫度380 K,冷卻深度的減弱也削減了霜凍的影響[22]。

    圖11 SABRE-4發(fā)動機熱力循環(huán)原理圖[3]Fig. 11 Principle diagram of the SABRE-4 engine thermal cycle[3]

    由于大量使用輕質(zhì)合金材料,SABRE-4發(fā)動機馬赫數(shù)5.0狀態(tài)推重比高達(dá)6~14,與噴氣發(fā)動機和沖壓發(fā)動機比均有明顯優(yōu)勢。圖12給出了SABRE發(fā)動機與其他類型發(fā)動機推重比對比結(jié)果[21,23]。吸氣模態(tài)能大幅提升發(fā)動機比沖性能,SABRE-4吸氣模態(tài)最大比沖能接近5200 s,火箭模態(tài)比沖460 s[24]。

    圖12 SABRE-4與其余類型發(fā)動機性能比較[23]Fig. 12 Performance comparison of SABRE-4 with other forms of engines[23]

    2009年,REL公司啟動了“渦輪基組合吸氣式火箭發(fā)動機關(guān)鍵技術(shù)試驗研究”項目(Experimental Investigation of Key Technologies for a Turbine Based Combined Rocket Engine),該項目獲得了歐洲空間局(ESA)的支持,REL公司計劃通過三個階段完成“佩刀”發(fā)動機關(guān)鍵技術(shù)的攻關(guān)[25]。第1階段(2009~2011年),REL聯(lián)合英國BAE系統(tǒng)公司開展三方面工作—氧化劑冷卻燃燒室技術(shù)、補償噴管技術(shù)和熱交換器模塊研制。2009年2月,ESA投資100萬歐元在該階段項目上[26]。第2階段(2011~2013年),主要目標(biāo)為驗證全尺寸熱交換器,同時也開始了新一輪進(jìn)氣道氣動問題研究。該階段開展了一系列風(fēng)洞試驗,2012年通過試驗將空氣冷卻至130 K,2013年完成全尺寸預(yù)冷器演示驗證試驗。空氣動力有限公司(Gas Dynamics Limited)的風(fēng)洞試驗也檢驗了變幾何進(jìn)氣道理論模型。第3階段(2013~2016年),REL公司主要開展SABRE地面靜態(tài)演示驗證發(fā)動機研究,檢驗SABRE發(fā)動機熱力循環(huán)。為支撐2016年SKYLON驗證機項目,研發(fā)了艙體驗證飛行器(Nacelle Test Vehicle,NTV),驗證進(jìn)排氣系統(tǒng)和飛發(fā)一體性能。該階段經(jīng)費開銷約2.25億英鎊[22]。

    預(yù)冷器作為SABRE發(fā)動機關(guān)鍵部件,技術(shù)進(jìn)展備受關(guān)注。按照原三階段計劃,REL公司在2009年前完成了全尺寸的預(yù)冷器模塊研制,采用鉻鎳鐵合金718制造了毛細(xì)管管徑1 mm、壁厚30 μm、耐壓達(dá)60 MPa的預(yù)冷器組件,見圖13(a)[20,25]。第2階段對預(yù)冷器的研究工作正常開展,試驗驗證了集成的預(yù)冷器(圖13(b))能在低于373 K的溫度下持續(xù)工作超過5 min。第3階段REL公司搭建了空氣來流溫度最高1250 K的預(yù)冷器熱考核試驗臺,見圖13(c)[27]。2019年4月該公司完成了SABRE-4發(fā)動機全尺寸預(yù)冷器馬赫數(shù)3.3條件高溫?zé)峥己嗽囼?,預(yù)冷器產(chǎn)生熱交換功率約1.5 MW,將700 K的空氣在0.05 s內(nèi)冷卻至373 K。同年10月,實施了預(yù)冷器馬赫數(shù)5.0熱考核試驗,能在0.05 s時間內(nèi)將來流高溫空氣冷卻1000 K以上。

    圖13 SABRE發(fā)動機預(yù)冷器及試驗臺Fig. 13 Pre-cooler of the SABRE engine and its test facility

    近幾年在突破預(yù)冷器技術(shù)后,REL公司獲得了英國和美國多家機構(gòu)的資金支持,建造了多個試驗平臺,完成了多項SABRE發(fā)動機關(guān)鍵技術(shù)試驗驗證[27-29]。

    歐洲空間局于2005年提出了“長期先進(jìn)推進(jìn)概念與技術(shù)”(Long-Term Advanced Propulsion Concepts and Technologies, LAPCAT)計劃。在LAPCAP項目牽引下,REL公司提出了400噸級高聲速民航機A2概念,并配套設(shè)計了馬赫數(shù)5預(yù)冷型發(fā)動機—Scimitar,如圖14所示[30],同樣采用氫為燃料,使用氦進(jìn)行閉式循環(huán)。

    圖14 Scimitar發(fā)動機示意圖[30]Fig. 14 Schematic of the Scimitar engine[30]

    Scimitar發(fā)動機屬于SABRE發(fā)動機的衍生型,基本原理和關(guān)鍵技術(shù)與SABRE發(fā)動機十分相似,但設(shè)計目標(biāo)不同。SABRE發(fā)動機目標(biāo)瞄準(zhǔn)單級/兩級入軌,Scimitar發(fā)動機瞄準(zhǔn)高聲速巡航,對比沖和使用壽命有更高要求,設(shè)計壽命15000 h。Scimitar發(fā)動機采用三激波進(jìn)氣道,動能效率0.9,安裝了6段預(yù)冷器,每段70模塊,使壓氣機進(jìn)氣總溫控制在635 K以內(nèi);壓比進(jìn)一步降低,采用兩軸對轉(zhuǎn)壓氣機,總壓比4.07。相對于SABRE發(fā)動機,Scimitar更貼近民用,在現(xiàn)有技術(shù)水平上可實現(xiàn)性更高[30]。

    Scimitar發(fā)動機有兩種工作模式,由渦扇發(fā)動機提供主推力的B模式,以及由ATR和旁路沖壓提供主推力的P模式,在馬赫數(shù)2.5進(jìn)行模態(tài)轉(zhuǎn)換。地面推力372 kN,燃料當(dāng)量比0.8,模態(tài)轉(zhuǎn)換推力272.8~313.1 kN。該型發(fā)動機可在亞聲速和超聲速進(jìn)行循環(huán),亞聲速可在馬赫數(shù)0.9以9688 s比沖巡航,超聲速保持馬赫數(shù)5以4195 s比沖巡航,經(jīng)濟性較好[30]。圖15給出Scimitar發(fā)動機兩種巡航模式結(jié)構(gòu)布局[31]。近幾年,由于沒有后續(xù)項目支持,Scimitar發(fā)動機進(jìn)展緩慢,大多數(shù)關(guān)鍵技術(shù)停留在方案階段。

    圖15 Scimitar發(fā)動機巡航狀態(tài)工作模式[31]Fig. 15 Operation modes of the Scimitar engine at cruise conditions[31]

    1.4 適度預(yù)冷開式膨脹循環(huán)——ATREX發(fā)動機

    1986年,為滿足低成本、可重復(fù)使用航天運載系統(tǒng)動力需求,日本宇航研究院(Institute of Space and Astronautical Science,ISAS)聯(lián)合多家研究機構(gòu)提出了膨脹循環(huán)空氣渦輪火箭發(fā)動機(Air-Turboramjet with Expander Cycle,ATREX)[32]。該型發(fā)動機采用低溫液氫燃料為冷卻劑,在進(jìn)氣道內(nèi)安裝預(yù)冷器對來流空氣進(jìn)行預(yù)冷,充分吸熱后的燃料達(dá)到高焓高壓狀態(tài),驅(qū)動渦輪帶動壓氣機工作,再噴入燃燒室燃燒。為防止冷卻管束表面結(jié)冰,在進(jìn)氣錐上加裝了介質(zhì)噴射裝置,通過向空氣內(nèi)噴射液氮/液氧/醇類來抑制結(jié)冰[33]。

    早期的ATREX發(fā)動機采用葉尖渦輪驅(qū)動壓氣機工作,如圖16(a),而葉尖渦輪的C/C復(fù)合材料氧化問題始終未有效解決,因此,2000年后又采用傳統(tǒng)后置渦輪布局代替原有的葉尖渦輪布局,如圖16(b),同時將渦輪進(jìn)氣溫度降至1200 K,發(fā)動機綜合性能與葉尖渦輪方案相當(dāng)[34]。

    圖16 ATREX發(fā)動機原理示意圖[34]Fig. 16 Schematic of ATREX engine[34]

    1992年之后日本開始對ATREX開展試驗研究,最早利用縮比模型對ATREX進(jìn)行了約30次地面試車,該模型在未安裝預(yù)冷器條件下測得1400 s的比沖和4500 N的推力[35]。預(yù)冷器作為ATREX發(fā)動機的核心部件之一,是ISAS研究的重點,1992~1994年,先后對三種布局(分別為Type-I、Type-II、Type-III)的預(yù)冷器開展了風(fēng)洞試驗,1995年后進(jìn)行了地面點火試驗[33]。預(yù)冷器的試驗結(jié)果證明采用液氮為冷卻劑最低能將來流空氣從293 K冷卻至209 K,若采用氫為冷卻劑則能降低至171 K,總壓恢復(fù)系數(shù)最高0.9以上[32]。圖17展示了ATREX發(fā)動機地面點火試驗臺。根據(jù)對三種布局預(yù)冷器地面試驗結(jié)果來看,Type-Ⅰ和Type-Ⅱ有較大的管徑和管束間距,總壓損失較小,但是熱交換率無法達(dá)到預(yù)計值。Type-Ⅲ減小管徑尺寸并縮小管束間距后,管壁容易迅速結(jié)冰,壓力損失陡增[36]。向進(jìn)氣道內(nèi)噴入乙醇能有效抑制預(yù)冷器表面結(jié)冰,結(jié)合地面試驗結(jié)果評估可知整個飛行過程(地面狀態(tài)至30 km)所需的乙醇總量不超過載油量的3%[37]。

    圖17 ATREX地面點火試驗臺[38]Fig. 17 Ground firing test of the ATREX engine[38]

    后期地面熱試驗(圖18、圖19)采用液態(tài)甲醇或者氣態(tài)甲醇來抑制結(jié)冰,能降低80%的壓力損失。以液氫作為冷卻劑和燃料。ATREX-12、ATREX-13和ATREX-14的地面點火試驗測得2400~3540 N的推力,空氣流量6.05~7.28 kg/s,比沖10200~16300 N·s/kg[33]。日本航空航天研究局(JAXA)也提出了基于ATREX發(fā)動機的飛行試驗樣機和天地往返飛行試驗,并開展商業(yè)應(yīng)用探索[38],但截止目前尚無商業(yè)應(yīng)用報告。

    圖18 液態(tài)甲醇ATREX-500發(fā)動機試驗臺[33]Fig. 18 Ground test of the ATREX-500 engine with liquid methanol injection[33]

    圖19 氣態(tài)甲醇ATREX-500發(fā)動機試驗臺[33]Fig. 19 Ground test of the ATREX-500 engine with gaseous methanol injection[33]

    1.5 其他預(yù)冷型組合循環(huán)發(fā)動機

    20世紀(jì)90年代末,Vladimir Balepin提出了將液體火箭發(fā)動機(Liquid Rocket Engine,LRE)和深度預(yù)冷渦噴發(fā)動機(Deeply Cooled Turbojet,DCTJ)有機結(jié)合而成的新型預(yù)冷型組合循環(huán)發(fā)動機——KLIN發(fā)動機[39](圖20),該型發(fā)動機的設(shè)計旨在為單級入軌(Single Stage to Orbit,SSTO)飛行器提供動力。發(fā)動機的火箭單元全速域參與工作。

    圖20 KLIN發(fā)動機原理圖[40]Fig. 20 Schematic of KLIN engine[40]

    經(jīng)預(yù)冷作用后的渦噴發(fā)動機最高可工作至馬赫數(shù)6.0~6.5。地面狀態(tài)渦噴發(fā)動機來流空氣被冷卻劑氫預(yù)冷至110 K,而馬赫數(shù)6.0狀態(tài)能被冷卻至200~250 K,期間空氣不會被液化。由于進(jìn)氣溫度低,渦噴發(fā)動機得以使用較高的壓比,提高了發(fā)動機熱循環(huán)效率。但在深度預(yù)冷的同時,又帶來了預(yù)冷器結(jié)冰的問題,KLIN發(fā)動機在預(yù)冷器前加裝了液氧噴注裝置來控制結(jié)冰,液氧噴注的另一個作用也是為了冷卻空氣。與空氣液化循環(huán)發(fā)動機相比,KLIN發(fā)動機簡化了裝置,減輕了重量,發(fā)動機得以獲取更高的推重比[40]。

    KLIN發(fā)動機可垂直起飛,亞聲速段火箭發(fā)動機也會提供較大比例的推力,跨聲速后火箭發(fā)動機開始節(jié)流,只提供較小的推力。馬赫數(shù)6.0左右,預(yù)冷渦噴發(fā)動機關(guān)閉,火箭發(fā)動機重新起動,工作模態(tài)轉(zhuǎn)換為火箭發(fā)動機單獨工作。低速段采用吸氣模式獲取氧化劑,減輕了所需攜帶液氧的重量,比沖約為液氫/液氧火箭發(fā)動機的兩倍。另外,垂直起飛可省去起落架及其附件的重量,進(jìn)一步減輕了飛行器起飛重量[41]。

    KLIN發(fā)動機的主要難點在于低溫高壓比壓氣機,馬赫數(shù)0~6.0變化時,壓氣機進(jìn)氣流量和總壓變化較大,需重新設(shè)計工作范圍更寬的壓氣機[16]。

    2003年日本JAXA在ATREX研究基礎(chǔ)上面向民用航空提出了預(yù)冷渦噴發(fā)動機(Pre-Cooled Turbojet,PCTJ)方案,如圖21所示,并加工了縮比尺度的預(yù)冷渦噴發(fā)動機S-PCTJ,用于驗證發(fā)動機循環(huán)可行性[42-44]。目前已完成進(jìn)氣道馬赫數(shù)1.5~5.0風(fēng)洞試驗、馬赫數(shù)2飛行試驗和馬赫數(shù)4風(fēng)洞試驗,結(jié)合CFD評估驗證了PCTJ發(fā)動機能完成馬赫數(shù)0~5連續(xù)工作[45-46]。

    圖21 S-PCTJ發(fā)動機結(jié)構(gòu)[43]Fig. 21 Schematic of S-PCTJ engine[43]

    西安航天動力研究所張蒙正等針對可重復(fù)使用航天運輸或臨近空間快速投放平臺動力需求,提出了一種預(yù)冷空氣渦輪火箭(Pre-Cooling Air Turbo Rocket,PATR)組合動力系統(tǒng)方案[47],如圖22所示,主體布局為中間的火箭發(fā)動機和環(huán)繞火箭發(fā)動機的旁路沖壓發(fā)動機,外加第三介質(zhì)循環(huán)系統(tǒng)。采用氦為第三介質(zhì),閉式循環(huán)方式對空氣進(jìn)行預(yù)冷,并驅(qū)動空氣壓氣機和氫泵工作。以液氫作為最終熱沉,對高溫氦進(jìn)行冷卻,再驅(qū)動氦壓氣機-渦輪,膨脹后的氫一部分進(jìn)入火箭發(fā)動機燃燒室燃燒,另一部分導(dǎo)入沖壓燃燒室燃燒。PATR發(fā)動機工作速域馬赫數(shù)0~6,火箭發(fā)動機在全速域內(nèi)工作,也是推力的主要來源,旁路沖壓工作速域初步設(shè)計為馬赫數(shù)1.8~6.0。

    圖22 PATR發(fā)動機原理[48]Fig. 22 Schematic diagram of PATR engine[48]

    對PATR發(fā)動機熱力循環(huán)系統(tǒng)(圖23)和部件關(guān)鍵技術(shù)的分析表明[49],該型發(fā)動機的技術(shù)難點在于熱力參數(shù)的匹配、輕質(zhì)高效預(yù)冷器設(shè)計及制造技術(shù)、燃燒室熱防護(hù)技術(shù)、高效氦渦輪技術(shù)和低溫空氣壓縮技術(shù)等,未來團(tuán)隊將持續(xù)對系統(tǒng)方案進(jìn)行優(yōu)化。

    圖23 PATR發(fā)動機熱力循環(huán)系統(tǒng)[49]Fig. 23 Thermodynamic cycle of the PATR engine[49]

    2 預(yù)冷型發(fā)動機熱力循環(huán)問題分析

    預(yù)冷技術(shù)是預(yù)冷型組合循環(huán)發(fā)動機的典型特征,但預(yù)冷技術(shù)的引入也為發(fā)動機的設(shè)計帶來了諸多問題。以下從熱力循環(huán)角度來分析預(yù)冷型組合循環(huán)發(fā)動機的主要技術(shù)難點。

    2.1 預(yù)冷與燃燒對冷卻劑流量的需求矛盾問題

    無論是空氣液化循環(huán)、空氣深冷循環(huán)還是空氣適度預(yù)冷循環(huán),都是以低溫燃料作為冷卻劑,消耗燃料熱沉對來流空氣進(jìn)行預(yù)冷和對燃燒室壁面進(jìn)行熱防護(hù)。然而,冷卻對燃料的需求和燃燒對燃料的需求在大部分工作狀態(tài)都無法保持相同,上述的數(shù)型發(fā)動機幾乎都面臨冷卻對燃料需求量大于燃燒需求量的問題,不僅降低了發(fā)動機比沖,也增加了飛行器載荷負(fù)擔(dān)。

    以典型預(yù)冷型發(fā)動機為例,建立預(yù)冷器熱交換計算模型和燃燒室釋熱模型,對預(yù)冷和燃燒對燃料需求進(jìn)行對比分析[50]。

    式中: ηPC為 預(yù)冷器熱交換效率;q和qmax為空氣與冷卻劑在預(yù)冷器內(nèi)的熱流及其最大值;Cc為燃料(冷卻劑)熱容量流率,即燃料質(zhì)量流率與定壓熱容的乘積;Cmin為冷卻劑與空氣熱容量流率的較小者;cp_c、cp_a、cp_g和cp_mix分別為燃料、空氣、燃?xì)夂突旌蠚獾谋葻?;Tc_1和Tc_2分別為冷卻劑在預(yù)冷器入口和出口的總溫;Ta_1和Ta_2分別為空氣在預(yù)冷器入口和出口的總溫;T6和T7為 混合氣和燃?xì)鉁囟?;?與 φb為燃料當(dāng)量比及燃燒室實際燃燒掉的燃料當(dāng)量比,γ為燃料混合空氣燃燒的化學(xué)恰當(dāng)比;H為燃料反應(yīng)焓;ηb為釋熱效率。式(1)和式(2)用于評估不同冷卻劑流量對空氣的冷卻效果,式(3)用于評估燃燒室燃燒釋熱后燃?xì)鉁囟取?/p>

    分別選擇氫和甲烷為燃料,計算不同飛行速域下不同當(dāng)量比燃料對空氣的預(yù)冷效果,計算結(jié)果如圖24所示。可以看出,相同冷卻劑當(dāng)量比條件下,氫對空氣的預(yù)冷效果高于甲烷。1.0~2.0倍當(dāng)量比的氫在馬赫數(shù)0~4.5速域內(nèi)能將空氣冷卻51~476 K,而相同流量的甲烷在馬赫數(shù)0~4.0范圍僅能將空氣冷卻24~182 K。供應(yīng)1~2倍當(dāng)量比的燃料流量,能將空氣進(jìn)行適度預(yù)冷,基本滿足能將渦輪型核心機工作速域擴展至馬赫數(shù)4.0以上的需求。

    圖24 氫與甲烷對空氣的預(yù)冷效果Fig. 24 Pre-cooling effect of hydrogen and methane on air

    經(jīng)換熱后的燃料將被噴入燃燒室進(jìn)行燃燒,但燃燒室能消耗的燃料不僅受氧化劑流量的限制,還受燃燒室材料耐高溫極限和燃燒室壁面熱防護(hù)效果限制。以加力燃燒室為例,假設(shè)燃燒室最高可承受溫度為2 473 K,計算燃料與空氣混合燃燒達(dá)到最高溫所需的燃燒的燃料當(dāng)量比,結(jié)果如圖25所示,圖中 φb為可燃燒掉的燃料當(dāng)量比??梢钥闯鰺o論是甲烷燃料還是氫燃料,供應(yīng)冷卻劑當(dāng)量比較大時燃燒室均能燃燒更多燃料,原因為更高的冷卻劑當(dāng)量比降低了來流總溫,燃燒室入口溫度也得以降低,因此在燃燒室溫度極限內(nèi)可燃燒更多燃料。

    圖25 燃燒室溫度限制內(nèi)燃燒所需燃料Fig. 25 Fuel requirement within the temperature limit of the combustor

    對甲烷燃料,馬赫數(shù)1.5以下,燃燒室基本能燃燒1倍當(dāng)量比的燃料。馬赫數(shù)1.5~4.0能燃燒0.8~1.0倍當(dāng)量比的燃料,氫的熱值更高,因此馬赫數(shù)0~4.5內(nèi)在燃燒室溫度極限內(nèi)僅能燃燒掉0.6~0.8倍當(dāng)量比。很明顯,燃燒室無法消耗完預(yù)冷需求的燃料量,若將過剩的燃料直接排放,會造成資源的浪費,并降低發(fā)動機比沖。

    為充分利用燃料熱沉和熱值,提高發(fā)動機綜合性能,需平衡發(fā)動機冷卻用燃料和燃燒用燃料。有以下幾項建議措施:(1)優(yōu)化預(yù)冷器布局,提高預(yù)冷器熱交換效率。預(yù)冷器熱交換效率的提高可使更少燃料達(dá)到相同預(yù)冷效果。(2)優(yōu)化預(yù)冷器結(jié)構(gòu),降低預(yù)冷器對空氣造成的壓損。預(yù)冷器會對空氣造成總壓損失,導(dǎo)致發(fā)動機推力相對降低。(3)提高燃燒室材料耐高溫極限。由以上評估計算可知,2473 K的燃燒室總溫上限無法滿足氫燃料燃燒需求,若能提高燃燒室壁面材料耐高溫極限,燃燒室便可承受更高當(dāng)量比的燃料燃燒。(4)采用更高效的熱防護(hù)措施,包括主動冷卻和被動冷卻。燃燒室壁面得到更好的熱防護(hù)措施后,能容許燃燒室有更高燃燒溫度。

    2.2 預(yù)冷器與壓氣機參數(shù)匹配問題

    預(yù)冷器的研制是預(yù)冷型組合循環(huán)發(fā)動機最關(guān)鍵的技術(shù)之一,其冷卻管道布局方式一般為管束式微小通道叉排分布,屬于管殼式換熱器的一種[51]。其中,日本ATREX發(fā)動機的預(yù)冷器為軸向?qū)ΨQ型(圖26),SABRE發(fā)動機的為周向?qū)ΨQ型(圖27)[52]。預(yù)冷器在發(fā)動機內(nèi)屬于上游部件,其性能對下游部件甚至發(fā)動機總體性能影響較大,相關(guān)性最高的無疑是壓氣機。一方面,預(yù)冷器降低了來流空氣總溫,讓壓氣機能在更高馬赫數(shù)下工作,壓氣機進(jìn)氣溫度的降低也容許壓氣機在材料耐高溫極限內(nèi)有更高的壓比;另一方面,預(yù)冷器的加裝給空氣帶來了總壓損失,需通過壓氣機來彌補。

    圖26 軸向?qū)ΨQ型預(yù)冷器結(jié)構(gòu)Fig. 26 Structure layout of the axial symmetry pre-cooler

    圖27 周向?qū)ΨQ型預(yù)冷器結(jié)構(gòu)Fig. 27 Baseline of circumferential symmetry pre-cooler

    將預(yù)冷器與壓氣機視作一體化部件,如圖28所示,引入熵函數(shù)對預(yù)冷器/壓氣機一體化部件參數(shù)匹配問題進(jìn)行設(shè)計。

    圖28 預(yù)冷器與壓氣機一體化結(jié)構(gòu)Fig. 28 Configuration of combined pre-cooler and compressor

    式中,Smc為預(yù)冷器/壓氣機一體化構(gòu)型的熵增函數(shù);和為 來流空氣和冷卻劑 的平均比熱;為 空氣流量;ka和kc分 別為空氣和冷卻劑的比熱比;Ta_3為壓氣機出口空氣溫度;πc和 τc分別為壓氣機壓比和溫比;ηc為壓氣機絕熱效率;Kcpm為冷卻劑與空氣的比熱流之比;σc為冷卻劑流經(jīng)預(yù)冷器的總壓恢復(fù)系數(shù)。式(4)中右側(cè)第一項意義為空氣在預(yù)冷器和壓氣機內(nèi)的溫度變化帶來的熵增,第二項為空氣在預(yù)冷器和壓氣機中由于壓力變化帶來的熵增,第三項為預(yù)冷器內(nèi)冷卻劑溫度和壓力變化帶來的熵增。

    引用的熵函數(shù)Smc反映了預(yù)冷器/壓氣機一體部件的熵增情況,是非目的的能量轉(zhuǎn)換,用于表達(dá)預(yù)冷器/壓氣機一體結(jié)構(gòu)在熱勢差效應(yīng)和功熱轉(zhuǎn)換過程的能量耗散總和,在較小的熵函數(shù)Smc條件下設(shè)計的預(yù)冷器與壓氣機匹配參數(shù)能獲得更高的循環(huán)效率。這種方法用一項指標(biāo)代替了原來由總壓恢復(fù)系數(shù)、預(yù)冷器熱交換效率、壓氣機絕熱效率等多個參數(shù)組成的評價體系,簡化了在多個性能參數(shù)間進(jìn)行折中選取的過程。

    圖29給出了預(yù)冷器消耗甲烷燃料當(dāng)量比 φ與壓氣機壓比 πc在 熵函數(shù)Smc圖上的分布,可以看出,相同壓比條件下,燃料當(dāng)量比越高則熵增越大,意味著燃料當(dāng)量比越高,發(fā)動機循環(huán)效率會越低;而在相同燃料當(dāng)量比條件下,壓比的提升也會一定程度上提高熵增,這是因為壓比的升高會帶來溫比的升高,溫比的升高帶來的熵增效果將超過壓比。圖30則給出了與圖29對應(yīng)工況下壓氣機出口總溫與燃料當(dāng)量比在熵函數(shù)圖上的分布。

    圖29 燃料當(dāng)量比與壓氣機壓比在熵函數(shù)圖上的關(guān)系Fig. 29 Relation between fuel equivalent ratio and compressor pressure ratio on the entropy function map

    圖30 燃料當(dāng)量比與壓氣機出口溫度在熵函數(shù)圖上的關(guān)系Fig. 30 Relation between fuel equivalent ratio and compressor exit temperature on the entropy function map

    根據(jù)發(fā)動機性能需求,可設(shè)計不同的當(dāng)量比-壓比(φ-πc)協(xié)同工作線,其中壓氣機出口溫度上限為900 K。若要求壓氣機壓比在馬赫數(shù)0~3各速度點均能達(dá)到最大值,則φ-πc工作線為圖29和圖30中工作線1.結(jié)合圖30可知馬赫數(shù)2以上限制壓比的主要因素為壓氣機耐高溫極限,此時馬赫數(shù)2狀態(tài)壓比能達(dá)到24.1,馬赫數(shù)3狀態(tài)能達(dá)到14.2。此工作線雖然有較高壓比,但熵增較高,能量耗散大,燃料消耗大,發(fā)動機熱循環(huán)效率也較低。

    為充分利用燃料熱沉和化學(xué)能,控制燃料當(dāng)量比在最高在1.0左右,保證燃料不會過剩,設(shè)計的φ-πc工作線則為圖29和圖30中的工作線2。此工作線顯示了一倍當(dāng)量比的燃料能在壓氣機耐高溫極限內(nèi)能使壓氣機在馬赫數(shù)2條件下獲得23的壓比,馬赫數(shù)3狀態(tài)仍有9的壓比。此工作線熵增不大,說明能量耗散不高,整個飛行速域壓比較高,燃料也沒有過剩情況,是比較理想的當(dāng)量比-壓比協(xié)同工作線。

    若要求發(fā)動機熱效率高,能量利用率高,則需保持熵增極小,可得到圖29和圖30中工作線3,此線始終保持較低熵增,消耗的燃料當(dāng)量比也低,熱循環(huán)效率高。但缺點也很明顯,即壓比無法提升,對發(fā)動機比推力和比沖性能帶來負(fù)面影響,另外,冷卻用燃料也不足以供應(yīng)燃燒室燃燒。

    2.3 壓氣機與渦輪共同工作問題

    壓氣機與渦輪共同工作問題也是預(yù)冷型組合循環(huán)發(fā)動機的關(guān)鍵技術(shù),與傳統(tǒng)渦輪發(fā)動機相比,最大的區(qū)別在于驅(qū)動渦輪工作的工質(zhì)由燃?xì)飧某闪烁哽矢邏籂顟B(tài)的可壓縮氣體,開式膨脹循環(huán)發(fā)動機一般為燃料(氫、甲烷等),閉式膨脹循環(huán)一般為氦。而對于大多數(shù)工況,渦輪工質(zhì)的質(zhì)量流量遠(yuǎn)低于來流空氣質(zhì)量流量,導(dǎo)致渦輪能為壓氣機提供的功率較小,限制了壓氣機壓比的提升。

    以典型開式膨脹循環(huán)預(yù)冷發(fā)動機核心機為例,建立了壓氣機壓比與渦輪落壓比數(shù)值計算模型,聯(lián)立流量平衡方程和功率平衡方程研究壓氣機壓比與渦輪落壓比之間的關(guān)系,以及壓氣機和渦輪單位功率變化規(guī)律[53],控制方程如下:

    式中, πT和 ηT分 別為渦輪落壓比和效率;pc_4和pc_5分別為燃料在渦輪入口和出口的壓力;Tc_4和Tc_5為燃料在渦輪入口和出口溫度;m˙c為 燃料質(zhì)量流;ηm為核心機的軸機械傳動效率。

    圖31展示了以甲烷為燃料的壓氣機壓比和渦輪落壓比之間的關(guān)系。假設(shè)渦輪落壓比最高為15,渦輪前溫度不超952 K,渦輪效率0.88~0.90,壓氣機效率0.77~0.82??梢钥闯?,地面狀態(tài)1倍當(dāng)量比燃料在渦輪最大落壓比范圍內(nèi)最高能使壓氣機壓比達(dá)到1.8,而壓氣機壓比要達(dá)到3以上則需要約2倍當(dāng)量比的燃料。馬赫數(shù)2.0狀態(tài),2倍當(dāng)量比的燃料在最高落壓比限制內(nèi)最高只能獲得2.8左右的壓比。到馬赫數(shù)4.0狀態(tài),來流空氣總溫和總壓明顯升高,壓縮更難,1~2倍當(dāng)量比的燃料僅能滿足1.3~1.7的壓比需求。

    圖31 壓氣機壓比與渦輪落壓比關(guān)系Fig. 31 Relation between compressor pressure ratio and turbine pressure ratio

    若在馬赫數(shù)0~4速域內(nèi)按照最大壓比設(shè)計πc?πT工作線,如圖31中工作線I所示,2倍當(dāng)量比燃料能在馬赫數(shù)0~4速域內(nèi)獲得3.35~1.6的壓比,此工作線消耗燃料較多,比沖低。若僅供應(yīng)1倍當(dāng)量比燃料,按照最高落壓比設(shè)計,則 πc?πT工作線如圖中工作線Ⅱ所示,此工作線有較高比沖,但最高壓比在1.9以下,壓比過低??紤]到較高的落壓比會給甲烷工質(zhì)的渦輪設(shè)計帶來較大設(shè)計難度,可適當(dāng)提高燃料當(dāng)量比,并降低渦輪落壓比,得到圖中工作線Ⅲ,此工作線當(dāng)量比1.5左右,最高落壓比不超過10,壓比范圍2.35~1.5,與氫工質(zhì)或燃?xì)怛?qū)動核心機渦輪相比,壓比仍較低。采用小流量甲烷作為驅(qū)動渦輪的工質(zhì),壓氣機壓比較小。

    對應(yīng)圖31壓氣機壓比變化規(guī)律在馬赫數(shù)0~4.0范圍內(nèi),考察壓氣機單位功率/渦輪單位功率與壓比/燃料當(dāng)量比的關(guān)系。圖32為不同速域下壓氣機單位功率隨壓比和當(dāng)量比的變化關(guān)系。圖33為不同速域下渦輪單位功率隨壓比和當(dāng)量比變化關(guān)系。由圖32可以看出,壓比相同情況下,燃料當(dāng)量比對壓氣機單位功率影響較弱,單位功率主要由壓比控制,飛行速度越高單位功率隨壓比增加越快。而無論是地面狀態(tài)還是高速狀態(tài),壓氣機單位功率上限均相近,約120 kw·s/kg。因此,可根據(jù)壓氣機最大單位功率來設(shè)計不同飛行速度下最大壓比。

    圖32 壓氣機單位功率變化Fig. 32 Variation of the compressor specific power

    與壓氣機相比,燃料當(dāng)量比對渦輪單位功率(相對于燃料流量)影響更明顯,當(dāng)量比越高,渦輪單位功率越低,如圖33所示。同樣可以看出,不同速域下以甲烷為工質(zhì)的渦輪最大單位功率達(dá)1200 kw·s/kg。高速狀態(tài)單位功率隨壓比上升迅速,達(dá)到最大值后渦輪無法滿足壓氣機更大壓比對功率的需求。因此,提高渦輪單位功率才能提高壓氣機單位功率和壓比。

    圖33 渦輪單位功率隨壓氣機壓比變化Fig. 33 Variation of the turbine specific power with respect to the compressor pressure ratio

    渦輪單位功率及驅(qū)動渦輪的工質(zhì)和渦輪落壓比相關(guān)。圖34顯示了甲烷工質(zhì)渦輪單位功率與渦輪落壓比關(guān)系??梢?,低速狀態(tài)較高單位功率要求落壓比達(dá)15以上,高速狀態(tài)需達(dá)8以上。與傳統(tǒng)核心機相比,預(yù)冷膨脹循環(huán)發(fā)動機的渦輪工質(zhì)流量與空氣流量相差甚大,導(dǎo)致渦輪的單位功率(相對于燃料流量)需達(dá)極高水平才能滿足壓氣機功率需求。更大的落壓比可以提高渦輪單位功率,但要求渦輪有更多級數(shù),且工質(zhì)有更高膨脹速度,再匹配壓氣機轉(zhuǎn)速相等,給渦輪的設(shè)計帶來了極大難度。預(yù)冷膨脹循環(huán)發(fā)動機壓氣機與渦輪的設(shè)計有待進(jìn)一步研究。

    圖34 渦輪單位功率隨渦輪落壓比變化Fig. 34 Variation of the turbine specific power with respect to the turbine pressure ratio

    3 結(jié)束語

    本文介紹了各型預(yù)冷組合循環(huán)發(fā)動機的發(fā)展歷程,并分析了各型發(fā)動機的優(yōu)缺點和技術(shù)難點。預(yù)冷型組合循環(huán)發(fā)動機從概念提出至如今,熱力循環(huán)方案經(jīng)歷了多次迭代改進(jìn),既有基本原理的繼承,也有新技術(shù)的加持。在總結(jié)了各型預(yù)冷組合循環(huán)發(fā)動機熱力循環(huán)特點后,對發(fā)動機熱力循環(huán)中的三個主要難點進(jìn)行了數(shù)值建模分析,結(jié)合計算結(jié)果對預(yù)冷型組合循環(huán)發(fā)動機典型關(guān)鍵技術(shù)概括如下:

    1)采用第三介質(zhì)閉式膨脹循環(huán)的預(yù)冷型發(fā)動機(SABRE和PATR等)是預(yù)冷型發(fā)動機未來重要的發(fā)展方向。該型發(fā)動機部件多,熱力循環(huán)系統(tǒng)復(fù)雜,存在諸多大溫差和大壓差節(jié)點,對部件性能和材料要求極高。除預(yù)冷器外還包括其余轉(zhuǎn)子部件,需要突破這些部件的輕質(zhì)、高效、小體積和高強度設(shè)計及制造技術(shù)。

    2)全速域范圍內(nèi)的熱力循環(huán)參數(shù)匹配是預(yù)冷型組合循環(huán)發(fā)動機的難點和關(guān)鍵。高速狀態(tài)需消耗大量冷卻劑,發(fā)動機預(yù)冷深度和比沖難以兼顧,根據(jù)發(fā)動機目標(biāo)性能要求調(diào)節(jié)冷卻劑流量才能獲取折中的性能。預(yù)冷型發(fā)動機預(yù)冷深度的優(yōu)化改進(jìn)正是反映了設(shè)計師們在尋求冷卻劑消耗量與發(fā)動機總體性能之間的平衡點。

    3)采用膨脹循環(huán)模式的預(yù)冷型組合發(fā)動機驅(qū)動渦輪工作的工質(zhì)流量低、焓差大,渦輪單位功率高、落壓比大,給壓氣機和渦輪的構(gòu)型設(shè)計以及協(xié)同工作控制策略設(shè)計造成極大困難,寬速域范圍內(nèi)該類型發(fā)動機的壓氣機與渦輪一體化設(shè)計與控制研究還較少,其關(guān)鍵技術(shù)有待突破。

    得益于設(shè)計和加工制造技術(shù)的快速發(fā)展,隨著預(yù)冷型組合循環(huán)發(fā)動機的進(jìn)一步研究,高性能的預(yù)冷組合循環(huán)發(fā)動機有望逐步進(jìn)入實際應(yīng)用。我國預(yù)冷型組合循環(huán)發(fā)動機研究目前還停留在方案論證和部件關(guān)鍵技術(shù)驗證上,發(fā)動機整體研制進(jìn)程與國外存在較大差距,在理論分析和機理試驗研究上還有大量工作亟需開展。

    猜你喜歡
    冷卻劑馬赫數(shù)預(yù)冷
    核電站主冷卻劑泵可取出部件一體化吊裝檢修工藝探索
    基于ANSYS的LNG接收站典型管道預(yù)冷分析
    煤氣與熱力(2022年4期)2022-05-23 12:44:58
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    小型LNG氣化站預(yù)冷方法對比分析
    煤氣與熱力(2021年7期)2021-08-23 01:11:10
    載荷分布對可控擴散葉型性能的影響
    反應(yīng)堆冷卻劑pH對核電廠安全運行影響研究
    中國核電(2017年1期)2017-05-17 06:10:13
    冷卻劑泄漏監(jiān)測系統(tǒng)在核電廠的應(yīng)用
    中國核電(2017年1期)2017-05-17 06:10:05
    不同預(yù)冷方式對紅提葡萄的預(yù)冷效果
    中國果菜(2015年8期)2015-12-10 06:37:50
    冷卻劑管道取樣管焊縫裂紋分析
    焊接(2015年8期)2015-07-18 10:59:14
    NF-6連續(xù)式跨聲速風(fēng)洞馬赫數(shù)控制方式比較與研究
    大香蕉久久网| 成人国产一区最新在线观看| 亚洲专区字幕在线| 一级毛片电影观看| 中文字幕人妻熟女乱码| 免费观看av网站的网址| 欧美激情 高清一区二区三区| 精品人妻在线不人妻| 91精品三级在线观看| 国产成人精品久久二区二区91| 美女主播在线视频| av片东京热男人的天堂| 丝袜美腿诱惑在线| 国产亚洲精品第一综合不卡| 国产亚洲一区二区精品| 久久久久久免费高清国产稀缺| 中文字幕精品免费在线观看视频| 99久久综合免费| 97精品久久久久久久久久精品| 久久国产精品男人的天堂亚洲| 老司机影院毛片| av福利片在线| 日本av手机在线免费观看| 国产精品一区二区精品视频观看| 欧美性长视频在线观看| 不卡av一区二区三区| 一级毛片精品| 欧美日韩亚洲国产一区二区在线观看 | 热re99久久国产66热| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 电影成人av| 操美女的视频在线观看| 一本大道久久a久久精品| 午夜免费成人在线视频| 肉色欧美久久久久久久蜜桃| 日本撒尿小便嘘嘘汇集6| 亚洲欧洲日产国产| 大片电影免费在线观看免费| 69精品国产乱码久久久| 一个人免费在线观看的高清视频 | 国产又爽黄色视频| 亚洲全国av大片| 国产一区有黄有色的免费视频| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 国产成人精品无人区| 91麻豆av在线| 下体分泌物呈黄色| 天天躁夜夜躁狠狠躁躁| 日韩制服丝袜自拍偷拍| 午夜激情av网站| 中文字幕最新亚洲高清| 国产极品粉嫩免费观看在线| 成人黄色视频免费在线看| 免费高清在线观看视频在线观看| 欧美精品av麻豆av| 欧美乱码精品一区二区三区| 老司机深夜福利视频在线观看 | 狠狠精品人妻久久久久久综合| 精品福利观看| 欧美精品啪啪一区二区三区 | 极品少妇高潮喷水抽搐| 国产免费福利视频在线观看| 国产精品 欧美亚洲| 欧美日韩av久久| 欧美黑人精品巨大| 久久性视频一级片| 美女福利国产在线| 狠狠精品人妻久久久久久综合| 美女中出高潮动态图| 啦啦啦啦在线视频资源| 一本综合久久免费| 婷婷丁香在线五月| 午夜成年电影在线免费观看| 国产片内射在线| 久久久久久久久免费视频了| 好男人电影高清在线观看| 亚洲久久久国产精品| 69av精品久久久久久 | 日韩欧美一区视频在线观看| 老熟妇仑乱视频hdxx| 老汉色∧v一级毛片| 韩国高清视频一区二区三区| 欧美日韩亚洲高清精品| 精品久久久久久久毛片微露脸 | 亚洲国产欧美在线一区| 夫妻午夜视频| 欧美另类亚洲清纯唯美| 少妇人妻久久综合中文| 国产黄色免费在线视频| 亚洲七黄色美女视频| 久久精品国产亚洲av高清一级| 这个男人来自地球电影免费观看| 日韩大码丰满熟妇| av片东京热男人的天堂| 精品一区在线观看国产| 亚洲午夜精品一区,二区,三区| 激情视频va一区二区三区| 精品福利观看| 欧美激情久久久久久爽电影 | 久久天堂一区二区三区四区| 1024香蕉在线观看| 99香蕉大伊视频| 国产成人免费无遮挡视频| 亚洲全国av大片| 精品国产一区二区三区久久久樱花| 各种免费的搞黄视频| 超碰成人久久| 丁香六月天网| 国产欧美亚洲国产| www.999成人在线观看| 午夜福利,免费看| 国产免费av片在线观看野外av| 国产在线观看jvid| 日韩精品免费视频一区二区三区| 久热这里只有精品99| 黄色视频,在线免费观看| 久久久久久亚洲精品国产蜜桃av| 青青草视频在线视频观看| 丰满少妇做爰视频| 久久热在线av| 高清av免费在线| 欧美黄色淫秽网站| 国产亚洲午夜精品一区二区久久| 久久国产精品人妻蜜桃| 日韩免费高清中文字幕av| 久久人人97超碰香蕉20202| 搡老岳熟女国产| 国产在线观看jvid| 丝瓜视频免费看黄片| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩一区二区三区在线| 精品国产超薄肉色丝袜足j| 亚洲成人免费电影在线观看| 精品视频人人做人人爽| 国产一区有黄有色的免费视频| a级片在线免费高清观看视频| 亚洲欧美精品自产自拍| √禁漫天堂资源中文www| 亚洲精品第二区| 精品一区二区三卡| 亚洲精品美女久久av网站| 蜜桃在线观看..| 另类亚洲欧美激情| 18禁国产床啪视频网站| 免费观看人在逋| 精品熟女少妇八av免费久了| 性少妇av在线| 国产日韩一区二区三区精品不卡| www.自偷自拍.com| 午夜福利在线观看吧| 免费在线观看日本一区| 国产xxxxx性猛交| 桃红色精品国产亚洲av| 日日爽夜夜爽网站| 欧美 日韩 精品 国产| 欧美午夜高清在线| 欧美黑人欧美精品刺激| 国产精品久久久人人做人人爽| 午夜影院在线不卡| av视频免费观看在线观看| 99re6热这里在线精品视频| 欧美在线一区亚洲| 丁香六月欧美| 欧美变态另类bdsm刘玥| 99国产精品一区二区蜜桃av | 午夜91福利影院| 亚洲午夜精品一区,二区,三区| 激情视频va一区二区三区| 91字幕亚洲| 男女无遮挡免费网站观看| 欧美日本中文国产一区发布| 成人三级做爰电影| 国产高清视频在线播放一区 | 777久久人妻少妇嫩草av网站| kizo精华| 最新的欧美精品一区二区| 美女扒开内裤让男人捅视频| 欧美另类亚洲清纯唯美| 亚洲成国产人片在线观看| 中文字幕色久视频| 亚洲 欧美一区二区三区| 水蜜桃什么品种好| 一二三四在线观看免费中文在| 丁香六月天网| 中文字幕人妻丝袜制服| 91大片在线观看| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜制服| 一级毛片精品| 18禁国产床啪视频网站| av又黄又爽大尺度在线免费看| 五月天丁香电影| 日本五十路高清| 亚洲精品美女久久av网站| 亚洲精品粉嫩美女一区| 免费高清在线观看日韩| 香蕉国产在线看| 首页视频小说图片口味搜索| 天堂8中文在线网| 久久久久精品国产欧美久久久 | 天堂8中文在线网| 又黄又粗又硬又大视频| 国产不卡av网站在线观看| 国产成人av教育| 岛国在线观看网站| 国产精品自产拍在线观看55亚洲 | 深夜精品福利| 大香蕉久久网| 成年美女黄网站色视频大全免费| 国产男女超爽视频在线观看| 国产成人精品无人区| 免费在线观看黄色视频的| 亚洲综合色网址| 天天躁日日躁夜夜躁夜夜| 极品少妇高潮喷水抽搐| 麻豆国产av国片精品| www日本在线高清视频| 高清视频免费观看一区二区| 波多野结衣一区麻豆| 久久综合国产亚洲精品| 老司机在亚洲福利影院| 国产精品一二三区在线看| 丝瓜视频免费看黄片| 午夜视频精品福利| 少妇被粗大的猛进出69影院| a级毛片在线看网站| 精品一区二区三区四区五区乱码| 国产成人av激情在线播放| 久久久久久免费高清国产稀缺| 肉色欧美久久久久久久蜜桃| 亚洲中文日韩欧美视频| a级片在线免费高清观看视频| 一本色道久久久久久精品综合| 婷婷色av中文字幕| 久久天躁狠狠躁夜夜2o2o| 亚洲成人手机| 国产精品久久久人人做人人爽| av网站免费在线观看视频| 视频区图区小说| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧洲精品一区二区精品久久久| 久久人人爽av亚洲精品天堂| 一边摸一边抽搐一进一出视频| 纯流量卡能插随身wifi吗| 18禁国产床啪视频网站| 一进一出抽搐动态| 亚洲精品在线美女| 国产免费av片在线观看野外av| 少妇人妻久久综合中文| 999久久久国产精品视频| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 999精品在线视频| 精品人妻熟女毛片av久久网站| 99精品久久久久人妻精品| 国产精品免费大片| 国产亚洲午夜精品一区二区久久| 精品国产国语对白av| 亚洲精品在线美女| 高清视频免费观看一区二区| 成人手机av| 亚洲精品国产一区二区精华液| 男人爽女人下面视频在线观看| 国产在线观看jvid| 性少妇av在线| 国产深夜福利视频在线观看| 最黄视频免费看| 又紧又爽又黄一区二区| 亚洲av欧美aⅴ国产| 91九色精品人成在线观看| 亚洲精品美女久久久久99蜜臀| 自拍欧美九色日韩亚洲蝌蚪91| 12—13女人毛片做爰片一| av有码第一页| 青草久久国产| e午夜精品久久久久久久| 熟女少妇亚洲综合色aaa.| 精品国内亚洲2022精品成人 | 国产主播在线观看一区二区| 男女免费视频国产| 国产男女内射视频| 国产日韩欧美在线精品| 女人精品久久久久毛片| 久久精品国产a三级三级三级| 国产色视频综合| 国产成+人综合+亚洲专区| 少妇精品久久久久久久| 久久久久精品人妻al黑| 母亲3免费完整高清在线观看| 69精品国产乱码久久久| 亚洲精品国产一区二区精华液| 又黄又粗又硬又大视频| 在线观看一区二区三区激情| 老司机影院毛片| 国产成人免费观看mmmm| 高清欧美精品videossex| 亚洲国产日韩一区二区| 一级毛片女人18水好多| 亚洲成av片中文字幕在线观看| 多毛熟女@视频| 国产一区有黄有色的免费视频| 久久久精品国产亚洲av高清涩受| 成年av动漫网址| 波多野结衣一区麻豆| 一边摸一边做爽爽视频免费| 交换朋友夫妻互换小说| 久久久久久久久免费视频了| 狠狠精品人妻久久久久久综合| 香蕉丝袜av| 国产av又大| 国产区一区二久久| 大码成人一级视频| av电影中文网址| 久久人人97超碰香蕉20202| 成在线人永久免费视频| 精品久久久久久电影网| 黑人巨大精品欧美一区二区mp4| 大型av网站在线播放| 一本大道久久a久久精品| 91精品三级在线观看| 日韩一卡2卡3卡4卡2021年| 97人妻天天添夜夜摸| 久久亚洲精品不卡| 天天添夜夜摸| 色精品久久人妻99蜜桃| 成年女人毛片免费观看观看9 | 九色亚洲精品在线播放| 成年动漫av网址| 久久精品亚洲av国产电影网| 黄网站色视频无遮挡免费观看| 一级片免费观看大全| 欧美激情久久久久久爽电影 | 国产亚洲av片在线观看秒播厂| 日本黄色日本黄色录像| a级毛片在线看网站| 免费久久久久久久精品成人欧美视频| 99久久99久久久精品蜜桃| 999久久久国产精品视频| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 一边摸一边抽搐一进一出视频| 国产精品自产拍在线观看55亚洲 | 蜜桃在线观看..| 美女大奶头黄色视频| 69精品国产乱码久久久| 国产人伦9x9x在线观看| 精品少妇一区二区三区视频日本电影| 久久狼人影院| 最近最新免费中文字幕在线| 亚洲精品久久午夜乱码| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 男人操女人黄网站| 免费一级毛片在线播放高清视频 | 国产成人影院久久av| 欧美日韩亚洲国产一区二区在线观看 | 美女午夜性视频免费| 国产精品成人在线| 国产1区2区3区精品| av线在线观看网站| 老汉色av国产亚洲站长工具| 男女免费视频国产| 国产成人一区二区三区免费视频网站| 精品高清国产在线一区| 午夜视频精品福利| 视频区欧美日本亚洲| 国产片内射在线| 他把我摸到了高潮在线观看 | 91大片在线观看| 日韩欧美国产一区二区入口| 大陆偷拍与自拍| 老司机福利观看| 侵犯人妻中文字幕一二三四区| 欧美黄色片欧美黄色片| 99国产综合亚洲精品| 欧美日韩一级在线毛片| 丰满少妇做爰视频| 一二三四在线观看免费中文在| 欧美激情高清一区二区三区| 免费一级毛片在线播放高清视频 | 久久这里只有精品19| 欧美少妇被猛烈插入视频| 一本久久精品| www.自偷自拍.com| 人成视频在线观看免费观看| 淫妇啪啪啪对白视频 | 高清视频免费观看一区二区| 亚洲国产av影院在线观看| 伊人久久大香线蕉亚洲五| 亚洲欧美成人综合另类久久久| 自线自在国产av| 极品人妻少妇av视频| 精品亚洲成a人片在线观看| 亚洲免费av在线视频| 亚洲成人免费av在线播放| 一区在线观看完整版| 97人妻天天添夜夜摸| 美女高潮到喷水免费观看| 99久久国产精品久久久| 啦啦啦视频在线资源免费观看| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 亚洲精品乱久久久久久| 久久久久网色| 成年人免费黄色播放视频| 大片电影免费在线观看免费| 青春草视频在线免费观看| 国产精品 欧美亚洲| 欧美变态另类bdsm刘玥| 久久人人爽av亚洲精品天堂| 亚洲中文av在线| 国产成人av激情在线播放| 亚洲精品久久成人aⅴ小说| 欧美日韩一级在线毛片| 一区二区三区乱码不卡18| 一个人免费看片子| 交换朋友夫妻互换小说| 色视频在线一区二区三区| 乱人伦中国视频| 免费观看av网站的网址| 丰满迷人的少妇在线观看| 天天影视国产精品| 亚洲激情五月婷婷啪啪| 久久 成人 亚洲| 久热这里只有精品99| 国产亚洲午夜精品一区二区久久| 免费在线观看完整版高清| 国产成人av激情在线播放| 国产精品久久久人人做人人爽| 精品一区二区三区四区五区乱码| 午夜两性在线视频| 一本一本久久a久久精品综合妖精| 人人妻人人澡人人看| 日韩免费高清中文字幕av| 国产熟女午夜一区二区三区| 大片免费播放器 马上看| 久久人妻福利社区极品人妻图片| 夜夜骑夜夜射夜夜干| 午夜精品国产一区二区电影| 亚洲久久久国产精品| 老司机午夜福利在线观看视频 | 最新在线观看一区二区三区| 成年人午夜在线观看视频| 国产免费现黄频在线看| 青春草亚洲视频在线观看| 国产亚洲av片在线观看秒播厂| 麻豆av在线久日| 亚洲欧美精品综合一区二区三区| 亚洲人成电影观看| www.自偷自拍.com| 欧美日韩亚洲国产一区二区在线观看 | 久久久国产一区二区| av天堂在线播放| 午夜精品久久久久久毛片777| 97人妻天天添夜夜摸| 叶爱在线成人免费视频播放| 国产高清国产精品国产三级| 日韩欧美国产一区二区入口| 伊人亚洲综合成人网| 精品高清国产在线一区| 国产主播在线观看一区二区| bbb黄色大片| 麻豆乱淫一区二区| 99国产精品一区二区蜜桃av | 国产人伦9x9x在线观看| 亚洲精品久久成人aⅴ小说| av又黄又爽大尺度在线免费看| 18禁观看日本| 国产黄色免费在线视频| 亚洲成人手机| 久久精品熟女亚洲av麻豆精品| 国产熟女午夜一区二区三区| 色老头精品视频在线观看| 色94色欧美一区二区| 久久久久久久久免费视频了| 美女中出高潮动态图| 午夜老司机福利片| 日韩中文字幕欧美一区二区| 母亲3免费完整高清在线观看| 国产亚洲一区二区精品| 午夜老司机福利片| 亚洲精品久久成人aⅴ小说| 精品一区二区三区av网在线观看 | 最近中文字幕2019免费版| 精品福利观看| 中文字幕人妻熟女乱码| 久久国产精品男人的天堂亚洲| 免费久久久久久久精品成人欧美视频| 波多野结衣av一区二区av| 国产麻豆69| 国产日韩欧美视频二区| 国产伦人伦偷精品视频| 欧美另类一区| 精品人妻一区二区三区麻豆| 一区二区三区精品91| 曰老女人黄片| 日本a在线网址| 亚洲欧美一区二区三区久久| 亚洲精品粉嫩美女一区| 久久国产精品人妻蜜桃| 欧美激情极品国产一区二区三区| 欧美国产精品va在线观看不卡| 免费人妻精品一区二区三区视频| 韩国高清视频一区二区三区| a在线观看视频网站| 久久国产精品影院| 日韩熟女老妇一区二区性免费视频| 一级毛片女人18水好多| 精品亚洲成a人片在线观看| 三级毛片av免费| 国产无遮挡羞羞视频在线观看| 十八禁人妻一区二区| 色播在线永久视频| 啦啦啦 在线观看视频| 久久ye,这里只有精品| 国产主播在线观看一区二区| 午夜福利影视在线免费观看| 国产高清国产精品国产三级| 精品一区二区三卡| 欧美另类亚洲清纯唯美| 在线观看一区二区三区激情| 亚洲综合色网址| 久久精品熟女亚洲av麻豆精品| 亚洲国产中文字幕在线视频| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 亚洲国产欧美日韩在线播放| 久久天躁狠狠躁夜夜2o2o| 男人爽女人下面视频在线观看| 丝袜美腿诱惑在线| 91精品国产国语对白视频| www.av在线官网国产| 精品第一国产精品| 亚洲av美国av| 夫妻午夜视频| 国产av一区二区精品久久| 午夜成年电影在线免费观看| 两个人免费观看高清视频| 亚洲精品国产一区二区精华液| 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 在线观看www视频免费| 最新的欧美精品一区二区| 岛国在线观看网站| 狠狠婷婷综合久久久久久88av| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女 | 欧美激情久久久久久爽电影 | 国产精品1区2区在线观看. | 亚洲伊人久久精品综合| 国产精品久久久久久精品古装| 大香蕉久久成人网| 久久久久久人人人人人| 亚洲国产av新网站| 宅男免费午夜| 一级毛片女人18水好多| 国产成人精品在线电影| 99九九在线精品视频| 国产欧美日韩精品亚洲av| 黄色怎么调成土黄色| 少妇被粗大的猛进出69影院| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 爱豆传媒免费全集在线观看| 少妇精品久久久久久久| 男女国产视频网站| 国产精品99久久99久久久不卡| 亚洲精品久久久久久婷婷小说| 啦啦啦啦在线视频资源| 这个男人来自地球电影免费观看| 日韩一区二区三区影片| 18在线观看网站| 秋霞在线观看毛片| 在线观看免费日韩欧美大片| 婷婷色av中文字幕| 日韩有码中文字幕| www.自偷自拍.com| 欧美黄色淫秽网站| 亚洲国产毛片av蜜桃av| 自拍欧美九色日韩亚洲蝌蚪91| 久久狼人影院| 国产日韩一区二区三区精品不卡| 国产精品九九99| 97人妻天天添夜夜摸| 99热网站在线观看| 日韩精品免费视频一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 99国产精品免费福利视频| 99国产精品99久久久久| 99精国产麻豆久久婷婷| 岛国在线观看网站| 超色免费av| 久久久久久久大尺度免费视频| 一区二区av电影网| videos熟女内射| 各种免费的搞黄视频| 新久久久久国产一级毛片| tocl精华| 乱人伦中国视频| 午夜激情久久久久久久| 久久毛片免费看一区二区三区| 精品一区二区三区av网在线观看 | 免费在线观看视频国产中文字幕亚洲 | 国产欧美日韩一区二区精品| 999久久久国产精品视频| 欧美黄色淫秽网站| 久久精品国产亚洲av香蕉五月 | 精品一区二区三区四区五区乱码| 亚洲男人天堂网一区| 精品少妇内射三级|