鄒正平,劉火星,唐海龍,萬敏,王洪偉,陳小龍,陳懋章,*

1.北京航空航天大學(xué) 能源與動力工程學(xué)院 航空發(fā)動機(jī)氣動熱力國家級重點實驗室,北京 100191 2.先進(jìn)航空發(fā)動機(jī)協(xié)同創(chuàng)新中心,北京 100191 3.北京航空航天大學(xué) 機(jī)械工程及自動化學(xué)院,北京 100191

伴隨著科技的進(jìn)步,人類對探索和利用空間越來越重視,因此對高超聲速飛行器的需求凸顯出來。高超聲速飛行器具備即時、高精度打擊能力,比現(xiàn)有各種隱形技術(shù)有更好的突防能力,且既具威懾性又具實用性,因而被稱為“常規(guī)的戰(zhàn)略武器”,將在未來的空天戰(zhàn)中成為主角。另一方面,現(xiàn)有的各種航天器多以多級火箭為動力,分級推動航天器進(jìn)入最終的空間軌道。這種方式的缺點是不能重復(fù)使用,酬載成本很高。降低有效載荷成本、重復(fù)使用一直是航天器動力追求的目標(biāo)。鑒于高超聲速飛行器所具備的巨大潛在戰(zhàn)略、戰(zhàn)術(shù)及經(jīng)濟(jì)價值,美、俄、歐、日等航空強(qiáng)國均已將其作為研究熱點,并給予高度重視和持續(xù)的資金支持,以期優(yōu)先占領(lǐng)此未來技術(shù)制高點。美國早已圍繞高超聲速飛行器及其動力裝置,開展了大量研究,如 Hy Fly(Hypersonic Flight Demonstrator)計劃、Hy Tech(Hypersonic Technology)計劃、Hyper-X(eXperimental Hypersonic)計劃、NGLT(Next Generation Launch Technology)計劃和高超聲速轟炸機(jī)計劃等,并陸續(xù)取得關(guān)鍵性的進(jìn)展[1]。如X-43A已成功進(jìn)行了7倍聲速和10倍聲速飛行試驗,X37B進(jìn)行了入軌飛行等。2013年,美國公布了其正在研發(fā)的SR-72高超聲速飛行器,飛行速度可達(dá)Ma=6,兼具偵察與攻擊雙重職能,計劃于2025年首飛,2030年裝備部隊。俄羅斯也在全力研制自己的高超聲速武器并給予優(yōu)先發(fā)展,以期盡快趕上美國。

動力是高超聲速飛行器最為關(guān)鍵的決定性技術(shù)。美國等航空強(qiáng)國針對高超聲速飛行器動力進(jìn)行了大量研究工作,但由于技術(shù)復(fù)雜,技術(shù)成熟度普遍較低,可用于水平起降重復(fù)使用高超聲速飛行器的動力技術(shù)更是如此,短時間還難以見到工程可用的方案。因此,動力已成為高超聲速飛行器的瓶頸所在。

若要降低航天器發(fā)射成本或?qū)崿F(xiàn)高超聲速飛行器水平起降,首先要考慮航空渦輪發(fā)動機(jī),在低馬赫數(shù)(0～3)飛行時,它是理想的動力,技術(shù)成熟、可靠性很高,其比沖是各類發(fā)動機(jī)中最高的,但在高馬赫數(shù)(＞3)時航空渦輪發(fā)動機(jī)性能急劇惡化且可靠性降低。而沖壓發(fā)動機(jī)在高馬赫數(shù)下性能優(yōu)越。因而,兼具航空渦輪發(fā)動機(jī)和沖壓發(fā)動機(jī)優(yōu)點的渦輪基沖壓組合發(fā)動機(jī)(Turbine-Based Combined Cycle,TBCC)成為了高超聲速飛行器動力的主要發(fā)展方向之一。但時至今日,TBCC仍難以解決模式轉(zhuǎn)換工程技術(shù)困難和“推力縫隙”、“背死重”以及高超聲速飛行時熱防護(hù)困難等問題,嚴(yán)重制約了其發(fā)展和工程應(yīng)用。

在高超聲速動力的探索過程中,預(yù)冷成為了主攻方向之一。所謂預(yù)冷是指將進(jìn)入發(fā)動機(jī)的高溫空氣進(jìn)行冷卻,使溫度降低到航空渦輪發(fā)動機(jī)能正常工作的溫度。各航空強(qiáng)國都提出過自己的預(yù)冷方案,并進(jìn)行了大量研究,如美國提出的射流預(yù)冷卻技術(shù)方案,日本和俄羅斯提出的預(yù)冷器技術(shù)方案,英國提出強(qiáng)預(yù)冷器與閉式循環(huán)組合的技術(shù)方案。在這些方案中最值得關(guān)注的是英國的方案。

英國噴氣發(fā)動機(jī)公司(Reaction Enging Limited,REL)在20世紀(jì)80年代提出了強(qiáng)預(yù)冷技術(shù)方案,30多年的沉寂之后于2012年11月宣布,該方案地面試驗取得成功,實現(xiàn)了強(qiáng)預(yù)冷,可將以Ma=5飛行的發(fā)動機(jī)進(jìn)口氣流滯止溫度降低到發(fā)動機(jī)能正常工作的進(jìn)氣溫度。以這種強(qiáng)預(yù)冷方案為基礎(chǔ),REL公司先后研發(fā)了兩種發(fā)動機(jī),即彎刀(Scimitar)和軍刀(Sabre,或譯為佩刀)。彎刀發(fā)動機(jī)采用常規(guī)的渦扇發(fā)動機(jī)技術(shù),加上強(qiáng)預(yù)冷系統(tǒng),用于Ma=5的高超聲速飛機(jī)。軍刀發(fā)動機(jī)從原理方案看,則是彎刀發(fā)動機(jī)與火箭發(fā)動機(jī)的組合,用于“云霄塔”航天飛行器Skylon,采用吸氣式發(fā)動機(jī)將飛行器加速到Ma=5.5,再用火箭進(jìn)一步推動飛行器,實現(xiàn)單級入軌。可將1 kg有效載荷送入地球低軌道的成本從目前的15 000英鎊/千克降至650英鎊/千克,與常規(guī)入軌方案相比,具有水平起降、重復(fù)使用、單級入軌、低酬載成本、兩次運行間隔時間短等一系列技術(shù)優(yōu)勢。

國際上對于此強(qiáng)預(yù)冷航空發(fā)動機(jī)技術(shù)給予了極高評價,認(rèn)為是“噴氣推進(jìn)發(fā)明以來的第二次革命”。英國大學(xué)與科學(xué)部長 Willetts發(fā)表評價:“這項技術(shù)會徹底改變未來的空中和太空旅行”,英國業(yè)內(nèi)人士預(yù)計2025年可最終在工程上實現(xiàn)這項技術(shù)。并認(rèn)為,強(qiáng)預(yù)冷航空發(fā)動機(jī)技術(shù)可使整個航空發(fā)動機(jī)面貌發(fā)生根本性變化,具有非常巨大的潛在的技術(shù)優(yōu)勢,可以把它看成是航空發(fā)動機(jī)領(lǐng)域的顛覆性技術(shù)。

本文將綜述國內(nèi)外高超聲速強(qiáng)預(yù)冷航空發(fā)動機(jī)相關(guān)技術(shù)的研究與進(jìn)展情況。

1 國外航空發(fā)動機(jī)強(qiáng)預(yù)冷技術(shù)研究

為解決高馬赫數(shù)飛行時,發(fā)動機(jī)進(jìn)口的高滯止溫度給發(fā)動機(jī)帶來的種種不利影響,對進(jìn)入發(fā)動機(jī)的空氣降溫的預(yù)冷技術(shù)日益受到重視,并得到了快速的發(fā)展。該技術(shù)即在常規(guī)發(fā)動機(jī)前增添預(yù)冷裝置,預(yù)先冷卻壓氣機(jī)進(jìn)口空氣。預(yù)冷可能帶來諸多好處:在高超聲速飛行時,降低進(jìn)氣溫度可擴(kuò)展高馬赫數(shù)飛行范圍,改善各部件工作條件,緩解發(fā)動機(jī)機(jī)體高溫防護(hù)問題；降低進(jìn)氣溫度可提高進(jìn)氣密度以增大進(jìn)氣質(zhì)量流量,進(jìn)而增大推力；高速滯止熱的適當(dāng)利用可提高循環(huán)熱效率。實現(xiàn)預(yù)冷的途徑主要有兩種:一是在壓氣機(jī)進(jìn)口噴入冷卻介質(zhì),其中的典型代表如美國MSE技術(shù)應(yīng)用公司提出的射流預(yù)冷卻TBCC發(fā)動機(jī)方案(Mass Injection Pre-compressor Cooling Turbine-Based Combined Cycle,MIPCC-TBCC)；二是利用預(yù)冷器,如日本的吸氣式渦輪沖壓膨脹循環(huán)發(fā)動機(jī)(Expander Air Turbo Ramjet Engine,ATREX)、俄羅斯的深冷空氣渦輪發(fā)動機(jī)(Deeply Cooled Air Turborocket,ATRDC)和英國REL公司等方案。

在射流預(yù)冷卻MIPCC-TBCC方案中,通過加裝在傳統(tǒng)渦輪發(fā)動機(jī)壓氣機(jī)前部的液體噴射裝置在進(jìn)氣道內(nèi)噴射冷卻介質(zhì),達(dá)到有效降低壓氣機(jī)進(jìn)口空氣溫度的目的,如圖1所示。美國在F100發(fā)動機(jī)基礎(chǔ)上發(fā)展MIPCC并進(jìn)行了試驗驗證,結(jié)果表明海平面推力增加了約1倍,在24.7 km的高度上最大飛行馬赫數(shù)可達(dá)3.5,可以有效拓寬原有飛行包線[2]。但是冷卻介質(zhì)的噴入會使得高空高速飛行時,壓氣機(jī)進(jìn)口來流氧含量下降,通常需要在壓氣機(jī)后注入氧化劑。同時該方案還存在比沖降低、發(fā)動機(jī)長度增加、來流總壓損失、進(jìn)氣總壓/總溫畸變以及所需冷卻介質(zhì)量大等問題。實際上為解決高溫大氣條件下起飛推力不足的問題,一個已得到實際應(yīng)用的辦法就是在進(jìn)氣道前方噴水。

圖1 MIPCC-TBCC發(fā)動機(jī)方案Fig.1 MIPCC-TBCC engine scheme

使用預(yù)冷器預(yù)冷的技術(shù)效率相對高,但設(shè)計難度增大,其中最為關(guān)鍵的是輕質(zhì)、高效的緊湊型預(yù)冷器的設(shè)計。日本提出的ATREX發(fā)動機(jī)方案即采用該技術(shù),通過安裝在進(jìn)氣道后的空氣/氫預(yù)冷器將進(jìn)口空氣冷卻[3-4](見圖2),并在預(yù)冷器前噴灑少量的酒精以有效阻止預(yù)冷器結(jié)冰[5]。地面縮尺模型試驗結(jié)果表明:當(dāng)風(fēng)扇進(jìn)口溫度低于160 K,且預(yù)冷器總壓恢復(fù)系數(shù)為0.9時,ATREX循環(huán)發(fā)動機(jī)的推力和比沖比無預(yù)冷時分別增加了2倍和1.5倍?；贏TREX發(fā)動機(jī)方案,日本宇航研究開發(fā)機(jī)構(gòu)(Japan Aerospace eX-ploration Agency,JAXA)進(jìn)一步研制了高速渦輪發(fā)動機(jī)(S-發(fā)動機(jī)),并于2014年2月成功進(jìn)行了在Ma=4條件下的地面試驗[6]。但是該方案發(fā)動機(jī)推重比較低,且無法保證整個飛行馬赫數(shù)范圍內(nèi)都擁有較高的比沖(其比沖在Ma=1,5時分別為3 500 s和2 500 s)[7],而且用氫直接冷卻,一旦預(yù)冷器發(fā)生氫脆問題,氫燃料泄漏到主流道,容易發(fā)生爆炸等安全問題；此外,它還存在冷卻所需的氫超過燃燒所需的氫,造成氫的浪費、換熱器單位質(zhì)量換熱量小以及結(jié)霜等問題。雖然存在諸多問題,S-發(fā)動機(jī)的實質(zhì)性進(jìn)展證明了預(yù)冷卻技術(shù)在提升渦輪發(fā)動機(jī)工作馬赫數(shù)方面具有巨大的潛力。

圖2 ATREX發(fā)動機(jī)方案Fig.2 ATREX engine scheme

俄羅斯在深冷空氣渦輪發(fā)動機(jī)(ATRDC)研究中同樣對預(yù)冷器冷卻技術(shù)進(jìn)行了探索。ATRDC也是采用燃料氫對進(jìn)口空氣進(jìn)行深度預(yù)冷并驅(qū)動渦輪,如圖3所示。該發(fā)動機(jī)擁有比ATREX更高的增壓比,壓氣機(jī)進(jìn)口溫度為98～112 K,壓比為40時,其在Ma=0～6范圍內(nèi)的平均比沖可達(dá)2 500 s,推重比高達(dá)18～22,但是預(yù)冷器很重,約占整個發(fā)動機(jī)質(zhì)量(不含進(jìn)氣道)的40%,它同樣存在預(yù)冷所需氫的量大于燃燒所需氫的量的問題,導(dǎo)致約1/2的氫在渦輪中膨脹后直接被排掉,造成了嚴(yán)重的資源浪費[8]。

圖3 ATRDC發(fā)動機(jī)方案Fig.3 ATRDC engine scheme

針對傳統(tǒng)TBCC和上述預(yù)冷動力方案的不足,在歐洲航天局和英國政府等的支持下,英國和比利時等聯(lián)合提出了一種適用于Ma=5的高超聲速航空發(fā)動機(jī)即“彎刀”(Scimitar)方案[9],這種發(fā)動機(jī)主要采用強(qiáng)預(yù)冷和閉式循環(huán)等新技術(shù),與TBCC相比,具有明顯優(yōu)勢:無論是低馬赫數(shù)還是高馬赫數(shù)飛行,耗油率皆可下降18%～23%；所有部件在全部飛行馬赫數(shù)范圍內(nèi)一直工作,不存在TBCC渦輪-沖壓模態(tài)之間的轉(zhuǎn)換和再啟動問題,以及渦輪或沖壓不工作時的“死重”問題；沒有TBCC模態(tài)切換時推力下降問題；由于強(qiáng)預(yù)冷,使高馬赫數(shù)時進(jìn)口溫度大幅下降,大大減輕了材料的高溫防護(hù)問題等,且與已有預(yù)冷動力方案相比不會大幅度增加發(fā)動機(jī)質(zhì)量,技術(shù)優(yōu)勢非常明顯。2012年11月,研制方噴氣發(fā)動機(jī)公司(REL)宣布地面試驗達(dá)到了強(qiáng)預(yù)冷的目標(biāo)(見圖4),即所發(fā)展的預(yù)冷技術(shù)具備在0.01 s瞬間將100 kg/s空氣的溫度由1 000℃降到-150℃的能力。這是一項意義非凡的重大技術(shù)突破,將整個改變航空發(fā)動機(jī)的面貌,有望成為未來最適用的高超聲速動力技術(shù)。

圖4 緊湊快速強(qiáng)換熱器試驗裝置Fig.4 Device diagram for compact precooler test

鑒于英國強(qiáng)預(yù)冷方案的巨大技術(shù)優(yōu)勢和已取得的重大突破,該方案已作為高超聲速飛機(jī)推進(jìn)技術(shù)的研究重點之一列入歐盟遠(yuǎn)期先進(jìn)推進(jìn)概念和技術(shù)計劃(Long-Term Advanced Propulsion Concepts and Technologies,LAPCAT),獲得了歐盟第六框架和第七框架計劃的經(jīng)費支持[10]。該項目預(yù)計將在2019年進(jìn)行整機(jī)測試,英國政府于2013年7月為該項目2014—2016年的研究投資6 000萬英鎊,為實現(xiàn)2025年發(fā)動機(jī)投入使用的目標(biāo),預(yù)計總投資80億歐元。該項目目前已吸引歐盟多個國家的關(guān)注和直接參與。同樣美國軍方對“軍刀”和“彎刀”發(fā)動機(jī)技術(shù)高度重視。2014年1月,美國空軍實驗室(Air Force Research Laboratory,AFRL)與英國噴氣發(fā)動機(jī)公司達(dá)成了合作研究和發(fā)展協(xié)議,要求與英方共享關(guān)鍵性的緊湊快速強(qiáng)換熱技術(shù),以支持美國的高超聲速飛機(jī)研發(fā)。同時,AFRL也開展了針對“軍刀”發(fā)動機(jī)的建模和性能評估工作,并于2015年初給出初步研究結(jié)果,認(rèn)為該新型發(fā)動機(jī)方案原理先進(jìn)可行,并且在工程實現(xiàn)上不存在重大技術(shù)障礙,他們將與REL繼續(xù)保持緊密合作進(jìn)一步深入研究相關(guān)關(guān)鍵技術(shù)。日本在2014年S-發(fā)動機(jī)試驗后,也改變原有研究思路,計劃在下一階段的預(yù)冷發(fā)動機(jī)研究過程中,參照英國“彎刀”方案思路,采用閉式換熱介質(zhì)循環(huán)和惰性冷卻介質(zhì)等相結(jié)合的改進(jìn)措施來進(jìn)一步發(fā)展高超聲速航空發(fā)動機(jī)技術(shù)。

2 緊湊快速強(qiáng)換熱機(jī)理及換熱器設(shè)計技術(shù)研究

強(qiáng)預(yù)冷航空發(fā)動機(jī)優(yōu)異性能的核心支撐是緊湊快速強(qiáng)換熱技術(shù),該創(chuàng)新性技術(shù)的核心是超臨界換熱介質(zhì)和微小尺度換熱單元的結(jié)合。處于超臨界狀態(tài)的流體,其黏性系數(shù)和擴(kuò)散系數(shù)接近氣態(tài)而傳熱系數(shù)和密度等接近液態(tài)。因此,用超臨界狀態(tài)流體作為換熱介質(zhì)可得到流動損失低、換熱量大的換熱效果。再結(jié)合幾何尺度為毫米或亞毫米級的基本換熱單元,可實現(xiàn)換熱器單位質(zhì)量換熱能力遠(yuǎn)超現(xiàn)有換熱技術(shù)的強(qiáng)換熱器。但是將超臨界與微小尺度二者相結(jié)合導(dǎo)致?lián)Q熱器內(nèi)部流動換熱機(jī)制更為復(fù)雜,必須對相關(guān)機(jī)理進(jìn)行研究,弄清楚超高熱流密度情況下微小尺度流動與換熱耦合機(jī)制及參數(shù)之間的關(guān)聯(lián)、超大的微通道長徑比對超臨界流動換熱的影響機(jī)制、超強(qiáng)溫度梯度引起的流動結(jié)構(gòu)變化及其與換熱的耦合機(jī)制等問題。

首要工作是機(jī)理等研究手段的發(fā)展,包括發(fā)展流熱耦合數(shù)值模擬程序,使之具備詳細(xì)模擬及分析緊湊快速強(qiáng)換熱器內(nèi)部復(fù)雜流動和換熱的能力；發(fā)展超臨界介質(zhì)微小尺度換熱流動實驗方法、流動測試方法并建立專用試驗平臺,解決緊湊快速強(qiáng)換熱器流動換熱機(jī)理試驗及性能試驗?zāi)芰栴}；探討該緊湊換熱技術(shù)在航空發(fā)動機(jī)設(shè)計中的應(yīng)用途徑,解決強(qiáng)預(yù)冷高超聲速航空發(fā)動機(jī)總體性能設(shè)計問題等。限于篇幅,有關(guān)試驗和發(fā)動機(jī)總體性能等方面的內(nèi)容,本文不作介紹。

2.1 高精度流熱耦合數(shù)值模擬方法及校驗

該數(shù)值模擬程序應(yīng)具備的主要特征是精度高、通用性好、功能豐富、對湍流模型和轉(zhuǎn)挨模型的調(diào)試改進(jìn)具有可擴(kuò)展性[11]。

1)采用基于預(yù)處理的時間推進(jìn)法,可求解由多孔材料、純流體以及純固體組成的耦合區(qū)域復(fù)雜流動和換熱問題。程序分別采用基于原始變量以及守恒變量的預(yù)處理方法來求解純流體域流動問題,而求解多孔區(qū)域流動與換熱問題只采用基于原始變量的預(yù)處理方法。

對于純流體域定常流場的預(yù)處理方法,其控制方程組可以寫為

式中:Ω為空間有界閉合區(qū)域；S為封閉曲面；W=[ρ ρvxρvθr ρvrρE]T為守恒變量,ρ 為密度,r為半徑,vx、vθ、vr分別為軸向、周向、徑向的速度,E為單位質(zhì)量流體的總能量；Fc、Fv和Q分別為對流項、擴(kuò)散項和源項。

經(jīng)過預(yù)處理后的Navier-Stokes方程,用原始變量Wp表示為

式中:Wp=[p vxvθvrT]T,T為溫度,p為壓力；ˉΓ為預(yù)處理矩陣[12]。

多孔介質(zhì)區(qū)域宏觀控制方程與純流體Navier-Stokes方程具有相似的方程形式,采用相似的預(yù)處理方法來解決時間推進(jìn)法在求解多孔域宏觀控制方程所遇到的方程剛性問題。根據(jù)不同的實際問題,多孔域能量方程采用局部非熱平衡(Local Thermal Equilibrium,LTNE)模型或者局部熱平衡(Local Thermal Nonequilibrium,LTE)模型。

2)添加基于局部變量的γ-Reθt轉(zhuǎn)挨模型,提高了有轉(zhuǎn)挨條件下的流熱耦合數(shù)值模擬精度。

Menter&Lantry提出的γ-Reθt轉(zhuǎn)挨模型是基于k-ω剪切應(yīng)力輸運(SST)湍流模型的四方程模型,這里只給出間歇因子γ輸運方程以及Reθt輸運方程。

守恒形式的間歇因子輸運方程為

式中:Tu為當(dāng)?shù)赝牧鞫龋籙為當(dāng)?shù)厮俣?；θ為動量厚度；s為流線弧長；F(λθ)代表了當(dāng)?shù)貕毫μ荻纫约巴牧鞫鹊挠绊?；γ為間歇因子；ui為速度分量；μT為渦黏系數(shù),μL和μ均為動力黏性系數(shù)；Pγ和Eγ分別為過渡源項和再層流化源項；λθ為流向壓力梯度。

3)添加保證通量守恒的非匹配流固交界面邊界條件,增強(qiáng)了程序處理復(fù)雜流熱耦合問題的能力。分區(qū)流熱耦合算法要求流固交界面上溫度和熱流通量連續(xù),程序中采用的是流體域向固體域傳遞邊界熱流量,固體域向流體域傳遞邊界溫度的方法。其中交界面上的溫度以及熱流密度采用了有限單元形函數(shù)插值方法進(jìn)行插值[13],在交界面上將固體域的溫度傳遞給流體域,同理將流體域的熱流密度傳遞給固體域,直到最終收斂。

4)添加了非定常模塊以及并行功能,實現(xiàn)了多孔/流體/固體非定常流熱耦合問題的快速求解。

通過一系列具有代表性的算例,分別對程序的對流項離散格式、各種湍流模型、轉(zhuǎn)挨模型、預(yù)處理方法、熱傳導(dǎo)求解器和流熱耦合方法等進(jìn)行了詳細(xì)的校驗和分析,結(jié)果表明各個模塊功能良好,具有較高的模擬精度和可靠性。

第一個校驗算例是對固體肋強(qiáng)化換熱進(jìn)行校驗,校驗算例來自文獻(xiàn)[14],平行通道上下壁面以交錯的方式各設(shè)置2個多孔/固體肋片,不可壓流體以恒定的速度以及溫度進(jìn)入平行通道,以通道高度和主流進(jìn)口來流條件定義的Re為100,通道內(nèi)流動為層流流動。

圖5(a)和圖5(b)分別給出了平行通道采用固體肋時其上下壁面的局部努塞爾數(shù)Nu分布,并與H.Y.Li采用Simpler方法得到的結(jié)果進(jìn)行了對比。圖中橫坐標(biāo)為x方向相對位置,縱坐標(biāo)為Nu,其中x為軸向當(dāng)?shù)匚恢?H為通道軸向長度。從中可以看出,由于固體肋的阻塞作用,使得流體以較高的速度流過肋片與通道之間的區(qū)域,故在下壁面設(shè)置肋片可以增強(qiáng)上壁面相應(yīng)位置處的換熱,Nu存在著峰值。本文結(jié)果與文獻(xiàn)結(jié)果符合得比較好,說明本程序能夠準(zhǔn)確地處理較為復(fù)雜的不可壓流熱耦合問題。

圖5 采用固體肋時通道上下壁面局部Nu分布(Re=100)Fig.5 Local Nu profiles at upper and lower walls of channel with solid baffles(Re=100)

圖6 可滲透壁面平行通道內(nèi)層流流動Fig.6 Laminar flow in a plate channel with permeable walls

圖7 不同注入率時軸向速度分布Fig.7 Axial velocity distributions at different coolant injection rates

第二個校驗算例是可滲透壁面層流平行通道流動。圖6給出了可滲透壁面平行通道示意圖,平行通道上下壁面均為可滲透壁面,流體以恒定的縱向速度Vw從下壁面流入通道,又以Vw的速度流出上壁面,通道內(nèi)流動為層流流動,且流動是不可壓的,當(dāng)通道內(nèi)的流動處于充分發(fā)展?fàn)顟B(tài)時,通道內(nèi)流體速度分布為橫向速度U(y)與縱向速度Vw的疊加[15]:

圖7給出了不同流體注入率時通道橫截面上橫向速度分布,并與解析解進(jìn)行了對比。圖中縱坐標(biāo)為相對位置,橫坐標(biāo)為橫向相對速度大小,其中Um為通道截面平均速度。從中可以看出,隨著流體注入率F的增大,通道內(nèi)的流體橫向速度分布越來越不對稱,橫向速度最大值變大,最大值位置也向上壁面移動。在不同的注入率下數(shù)值計算結(jié)果(離散點)與解析解(曲線)均符合得非常好,說明本程序計算的可靠性。

此外,還通過Rayleigh-Bénard流動、方腔周期振蕩頂蓋驅(qū)動流、多孔介質(zhì)栓塞流、多孔介質(zhì)中的Bénard流動等多個典型算例對該程序進(jìn)行了校驗和分析,驗證結(jié)果表明,該程序?qū)ξ⑿〕叨鹊亩嗫?流體/固體多區(qū)域流熱耦合問題具有較高的模擬精度,可以滿足緊湊快速強(qiáng)換熱器數(shù)值模擬要求[16]。

在發(fā)展流熱耦合數(shù)值模擬程序的同時,還針對緊湊快速強(qiáng)換熱器的特點,對商用軟件的預(yù)測精度進(jìn)行了校驗,分析湍流模型、邊界條件、浮升力以及流動加速等因素對緊湊快速強(qiáng)換熱器內(nèi)流動和換熱的影響。圖8給出了用商用軟件計算緊湊快速強(qiáng)換熱器內(nèi)流動得到的結(jié)果。圖中橫坐標(biāo)為沿流向的幾何相對位置,縱坐標(biāo)為溫度,其中d為換熱單元內(nèi)徑。黑色方點為實驗數(shù)據(jù)[17],各曲線為使用不同湍流模型得到的數(shù)值模擬結(jié)果。數(shù)值模擬結(jié)果與實驗結(jié)果的對比表明,合理地選取湍流模型能夠得到更為準(zhǔn)確的計算結(jié)果。

圖8 緊湊快速強(qiáng)換熱器內(nèi)流動和換熱計算結(jié)果Fig.8 Computing results of flow and heat transfer in compact precooler

圖9 多孔介質(zhì)強(qiáng)化換熱計算模型Fig.9 Computational models for heat transfer enhancement in porous media

2.2 緊湊快速強(qiáng)換熱機(jī)理

利用2.1節(jié)發(fā)展和校驗的數(shù)值模擬手段,本文對緊湊快速強(qiáng)換熱器內(nèi)部流動換熱機(jī)理進(jìn)行了研究,分別從換熱器內(nèi)、外換熱兩個方面入手,目的在于以最小的流動損失代價得到最大的強(qiáng)化換熱效果。

在內(nèi)換熱方面,首先對多孔介質(zhì)微小尺度強(qiáng)化換熱機(jī)理及Re、孔隙率、達(dá)西數(shù)、多孔介質(zhì)層厚度等參數(shù)的影響規(guī)律進(jìn)行了研究。采用多孔材料是實現(xiàn)強(qiáng)化換熱的有效手段,已有研究表明,由于多孔材料可以擁有微米甚至納米級的微觀結(jié)構(gòu),流體在流過多孔介質(zhì)時產(chǎn)生強(qiáng)烈混合,從而能夠大大改善換熱效果。限于篇幅,本文只給出了Re、達(dá)西數(shù)等參數(shù)的影響規(guī)律,分析表明高滲透率、高熱傳導(dǎo)系數(shù)以及高比表面積的多孔材料能更有效地增強(qiáng)換熱性能。雖然多孔材料能夠有效強(qiáng)化換熱,但是由于目前多孔材料在加工過程中的不確定性較大,難以實現(xiàn)關(guān)鍵參數(shù)的精確控制,距離實際應(yīng)用還有較大差距。本文進(jìn)一步對微小尺度圓管作為換熱基本單元情況下,超臨界介質(zhì)流動換熱、微小尺度流動換熱機(jī)理以及二者的耦合作用機(jī)制進(jìn)行了探索,初步得到了微小尺度圓管內(nèi)部超臨界介質(zhì)流動換熱的關(guān)鍵特征參數(shù)影響規(guī)律。研究結(jié)果表明以超臨界和微小尺度為特征的內(nèi)換熱可以得到很高的對流換熱系數(shù),完全可以滿足緊湊快速強(qiáng)換熱的要求。

在外換熱方面,主要考慮如何增大單位介質(zhì)換熱面積和提高外表面對流換熱系數(shù)。緊湊快速強(qiáng)換熱器工作在航空發(fā)動機(jī)環(huán)境中,外換熱介質(zhì)為空氣或燃?xì)?提高外換熱能力最為直接的方法是減小基本換熱單元的徑向尺寸,實現(xiàn)微小尺度換熱。而提高對流換熱系數(shù)則要從流動結(jié)構(gòu)入手,通過控制流動結(jié)構(gòu)的發(fā)展,實現(xiàn)以最小的流動損失為代價提高外表面對流換熱系數(shù)的目的。

2.2.1 多孔介質(zhì)微小尺度強(qiáng)化換熱機(jī)理

為實現(xiàn)強(qiáng)化內(nèi)換熱的目的,利用上述發(fā)展和校驗的多孔/流體/固體多區(qū)域三維流熱耦合數(shù)值模擬程序,對交錯肋片平行通道[18]、受限層流沖擊射流[19]以及層流平板發(fā)散冷卻[20]這3個多孔介質(zhì)強(qiáng)化換熱典型應(yīng)用中的流動換熱機(jī)理進(jìn)行了研究,其模型如圖9所示。圖9(a)中Lt為通道總長,Lc為通道中帶肋部分總長,U為氣流速度,L為肋間距,Li為第一個肋片與進(jìn)口間長度,t為肋厚,h為肋高,H為通道高度；圖9(b)中L為多孔介質(zhì)層長度,h為多孔介質(zhì)層厚度,B為槽寬,H為計算域厚度,vin為氣流進(jìn)口速度,Tin為氣流進(jìn)口溫度,Tw為底面溫度；圖9(c)中H為多孔介質(zhì)層厚度,Vc為冷卻流體速度,Tc為冷卻流體溫度,Uin為主流速度,Tin為主流溫度,δv為邊界層位移厚度。

圖10 不同Re、Da下壁面的局部Nu分布Fig.10 Local Nu profiles at lower wall with different Reynolds numbers and Darcy numbers

圖10 (a)給出了設(shè)置有交錯多孔/固體肋片的平行通道內(nèi)不同雷諾數(shù)時下壁面的局部Nu分布,橫坐標(biāo)為相對位置,縱坐標(biāo)為局部Nu(Nu=h L/λ,h為對流換熱系數(shù),L為特征長度,λ為導(dǎo)熱系數(shù))。從圖中可以看出,下壁面處的局部Nu分布呈現(xiàn)出周期性分布,由于多孔肋對于流體的阻礙作用,大部分流體會從多孔肋與下壁面的通道通過,會增強(qiáng)該處的對流換熱效果；另一方面,水平流動流體會逐漸改變方向,從而對下壁面有沖擊作用,也會增強(qiáng)該處的對流換熱作用；之后在兩個連續(xù)肋片之間,Nu出現(xiàn)極小值。另外,隨著Re的增大,Nu分布曲線整體上移,換熱能力增強(qiáng)。圖10(b)給出了不同Da數(shù)(Da=K/L2,K為多孔介質(zhì)滲透率,L為特征長度)時通道上下壁面的局部Nu分布,橫坐標(biāo)為相對位置,縱坐標(biāo)為局部Nu。顯然,Da對通道壁面局部Nu分布總體趨勢影響不大,但對Nu的峰值大小有很大的影響。隨著Da的減小,采用多孔肋的換熱特性與采用固體肋片的相類似。與固體肋片相比,采用多孔肋會使得流體傾斜沖擊壁面的能力減弱,故Nu極大值會減小,同時流體能夠穿透多孔肋會改善肋片之后的對流換熱,Nu極小值峰值也會有所削弱,這說明采用多孔肋會在整體上削弱換熱,但能使得通道壁面的換熱更加均勻。隨著Da的減小,多孔肋對流體的阻礙作用越來越小,Nu峰值進(jìn)一步減小。

圖11(a)給出了受沖擊壁面上鋪設(shè)有多孔介質(zhì)層的受限層流沖擊射流在不同Re時通道下壁面的Nu分布,圖中橫坐標(biāo)為相對位置,縱坐標(biāo)為Nu,ε為孔隙率,B為槽寬,Bi為畢渥數(shù),λs為固體導(dǎo)熱系數(shù),λf為流體導(dǎo)熱系數(shù),Cf為阻力系數(shù)。從中可以看出隨著Re的增大,滯止區(qū)的換熱性能逐漸增大,之后均沿流向急劇減小。另外,由于多孔介質(zhì)的存在使得光通道中下壁面區(qū)域存在的二次回流區(qū)消失,故相應(yīng)的二次峰值現(xiàn)象也消失,相應(yīng)Nu分布更加平滑。圖11(b)給出了相應(yīng)條件下不同Da時下壁面Nu的流向分布,橫坐標(biāo)為相對位置,縱坐標(biāo)為Nu。從中可以看出,隨著Da的增大,穿透多孔層沖擊到下壁面的冷卻流體增多,滯止區(qū)的Nu隨之增大；當(dāng)Da小到通道內(nèi)重新出現(xiàn)二次回流區(qū)時,下壁面Nu分布相應(yīng)的會出現(xiàn)二次峰值現(xiàn)象,但不明顯。

圖12(a)給出了典型平板發(fā)散冷卻中不同冷氣注入率時的多孔域/流體域交界面無量綱溫度分布,圖中橫坐標(biāo)為相對位置,縱坐標(biāo)為無量綱溫度(θ=(Ti-Tc)/(Tin-Tc),其中Ti為交界面溫度,Tc為冷卻介質(zhì)進(jìn)口溫度,Tin為主流進(jìn)口溫度)。從圖中可以看出,隨著冷氣注入率F的增加,多孔域/流體域交界面的溫度顯著降低,且隨著冷卻流體的不斷積累,越往下游邊界層厚度越大,相應(yīng)的溫度下降幅度越高。在本文的計算條件下,冷氣注入率不到0.2%時便能使得相對溫度下降到0.5,說明發(fā)散冷卻的冷卻效果非常好。圖12(b)給出了其他條件不變,不同主流Re時多孔域/流體域交界面的無量綱溫度分布,圖中橫坐標(biāo)為相對位置,縱坐標(biāo)為無量綱溫度。從中可以看出,Re的改變對于層流發(fā)散冷卻的冷卻效果還是有比較明顯的影響；在本文的計算條件下,隨著Re的增加,多孔域/流體域交界面上的溫度不斷降低,其冷卻效果不斷減弱。對于層流發(fā)散邊界層來說,Re主要改變的是邊界層特性,且隨著Re的增加,邊界層內(nèi)黏性力的作用減弱,邊界層厚度不斷降低,相應(yīng)的交界面上溫度和對流換熱系數(shù)升高,冷卻效果變差。另外,在本文計算條件下,不同的數(shù)量級時Re對冷卻效果的影響也不相同,在Re較小時,層流發(fā)散冷卻的冷卻效果對Re的變化更加敏感。

圖11 不同Re、Da時下壁面Nu分布Fig.11 Nu profiles at lower wall with different Reynolds numbers and Darcy numbers

圖12 不同主流Re、冷氣注入率多孔域/流體域交界面無量綱溫度分布Fig.12 Dimensionless temperature profiles at porous/fluid interface at different Reynolds numbers and coolant injection rates

2.2.2 微小尺度圓管內(nèi)部超臨界介質(zhì)流動換熱規(guī)律

提高壁面對流換熱系數(shù)是提高換熱能力的主要手段之一,研究氣動熱力幾何參數(shù)對換熱規(guī)律的影響,可為緊湊快速強(qiáng)換熱器一維設(shè)計參數(shù)的選取提供依據(jù)。限于篇幅,本文只給出圓管直徑D、進(jìn)口Re、進(jìn)口質(zhì)量流量密度G等參數(shù)對壁面換熱系數(shù)影響規(guī)律的分析。

圖13給出了在保證進(jìn)口Re相同和進(jìn)口質(zhì)量流量密度G相同的情況下,圓管直徑D對圓管壁面對流換熱系數(shù)的影響,其中橫坐標(biāo)為距入口的無量綱管徑D/D0(其中D0為研究中選取的特征直徑),縱坐標(biāo)為壁面對流換熱系數(shù)h,在這兩種情況下均保證對不同管徑D的流體單位質(zhì)量加熱量相同。從圖中可以看出,兩種情況下圓管壁面的對流換熱系數(shù)均隨D的減小而增大,并且存在一臨界直徑。

圖13 無量綱管徑對圓管壁面對流換熱系數(shù)的影響Fig.13 Heat transfer coefficient profiles at wall of tube for different dimensionless tube diameters

圖14 圓管壁面對流換熱系數(shù)隨進(jìn)口Re的變化關(guān)系Fig.14 Heat transfer coefficient profiles at wall of tube with different Reynolds numbers of inlet

圖14 給出了管徑D相同時,圓管壁面對流換熱系數(shù)隨進(jìn)口Re的變化關(guān)系,其中橫坐標(biāo)為進(jìn)口Re,縱坐標(biāo)為壁面對流換熱系數(shù)h,保證進(jìn)口Re不同時流體單位質(zhì)量加熱量相同。從圖中可以看出,圓管壁面換熱系數(shù)隨進(jìn)口Re的增大而增大。這是由于隨著進(jìn)口Re的增大,速度邊界層厚度減小,而加熱量相同意味著在相同流向位置處的物性參數(shù)相同,即普朗特數(shù)Pr相同,此時的溫度邊界層厚度也隨著隨度邊界層的減小而減小,因此壁面對流換熱系數(shù)增加。上述研究結(jié)果表明,減小管徑D和提高進(jìn)口Re均可以有效提高圓管壁面換熱系數(shù)h。

2.2.3 微小尺度圓管內(nèi)超臨界介質(zhì)流動換熱熵產(chǎn)分析

良好的換熱能力無疑是緊湊快速強(qiáng)換熱器重要的性能指標(biāo)之一,但若從氣動角度考慮,由于管內(nèi)被加熱的超臨界氦及管外預(yù)冷空氣均將最終進(jìn)入渦輪內(nèi)部進(jìn)行做功,因此有效控制換熱過程中造成的管內(nèi)外工質(zhì)做功能力的下降,對于航空發(fā)動機(jī)用換熱器的設(shè)計而言是極為重要的,即換熱器的設(shè)計追求換熱效率與流動損失等指標(biāo)的綜合性能最優(yōu)。而流動和換熱等因素引起的熵產(chǎn)之間關(guān)系的分析有助于理解換熱過程中熵產(chǎn)生的原因,可為設(shè)計提供參考。本文將熵產(chǎn)Sgen分為兩大部分,一部分為黏性耗散引起的耗散熵產(chǎn)SgenF,另一部分是由于導(dǎo)熱過程存在溫差所引起的加熱熵產(chǎn)SgenH,其表達(dá)式為

式中:u為速度；ε為湍流耗散系數(shù)；λ為導(dǎo)熱系數(shù)；α為熱擴(kuò)散系數(shù)；αT為湍流熱擴(kuò)散系數(shù)；δij為克羅內(nèi)克符號。式(9)等號右邊兩項分別代表統(tǒng)計平均流場內(nèi)的黏性耗散熵產(chǎn)與由湍流脈動導(dǎo)致的黏性耗散熵產(chǎn),式(10)等號右邊兩項分別代表統(tǒng)計平均流場內(nèi)的換熱熵產(chǎn)與由湍流脈動導(dǎo)致的換熱黏性耗散熵產(chǎn)。

本文結(jié)合對不同管徑換熱管內(nèi)超臨界氦的湍流流動數(shù)值模擬,對流場內(nèi)總熵產(chǎn)及其兩個組成部分的情況進(jìn)行分析。圖15為相同來流雷諾數(shù)下熵產(chǎn)隨管徑的變化趨勢,圖16為相同質(zhì)量流量密度下熵產(chǎn)隨管徑的變化趨勢。

圖15 相同來流雷諾數(shù)下熵產(chǎn)隨無量綱管徑的變化Fig.15 Entropy generation profiles for different dimensionless tube diameters with the same freestream Reynolds number

圖16 相同質(zhì)量流量密度下熵產(chǎn)隨無量綱管徑的變化Fig.16 Entropy generation profiles for different dimensionless tube diameters with the same mass flux

首先從圖15(b)中可以看出,隨管徑的增加,耗散熵產(chǎn)呈現(xiàn)大幅度減小,造成這一現(xiàn)象的原因在于,在相同的來流雷諾數(shù)下,管徑越小,來流速度越大,這兩點共同造成流場徑向速度梯度的增加,從而增強(qiáng)了平均流場內(nèi)的剪切耗散作用,最終使流動耗散隨管徑的減小而增強(qiáng)。而加熱熵產(chǎn)隨管徑的減小則是由于小管徑內(nèi)溫度梯度較大所致。雖然在同樣的熱流密度下小管徑內(nèi)由于換熱系數(shù)增強(qiáng)使得平均壁面-流體溫差減小,但由于該減小程度不如管徑的減小程度,而造成溫度梯度隨管徑的增加。導(dǎo)致相同來流雷諾數(shù)下熵產(chǎn)隨管徑的減小而增大。對于圖16所示情況,由于來流速度相同,同樣會出現(xiàn)在管徑減小下的速度梯度增加,造成耗散熵產(chǎn)的加劇。而導(dǎo)熱溫差導(dǎo)致的加熱熵產(chǎn)也同樣是管徑減小導(dǎo)致的溫度梯度增大所致。從熵產(chǎn)分析結(jié)果可以初步看出,減小管徑D雖然可以增強(qiáng)換熱性能,但相應(yīng)地也會造成流場各處總熵產(chǎn)的升高,導(dǎo)致流動系統(tǒng)能量品質(zhì)的降低更為嚴(yán)重。所以,設(shè)計時管徑的選擇應(yīng)綜合考慮流動損失與換熱效果。

圖17 單排微小尺度圓柱管束總壓恢復(fù)系數(shù)和對流換熱系數(shù)隨進(jìn)口Re變化Fig.17 Total pressure recovery and heat transfer coefficient profiles with different Reynolds numbers of inlet at a single row of micro-size tubes

圖18 多排微小尺度圓柱管束總壓恢復(fù)系數(shù)和對流換熱系數(shù)隨進(jìn)口Re變化Fig.18 Total pressure recovery and heat transfer coefficient profiles with different Reynolds numbers of inlet at micro-size tube matrix

2.2.4 微小尺度圓柱擾流換熱機(jī)理

采用緊湊排列圓柱管束來實現(xiàn)強(qiáng)化外換熱的同時會增大流動損失。因此,需要通過對緊湊排列圓柱管束中流動換熱參數(shù)影響規(guī)律的研究,得到綜合考慮流動損失和換熱性能后的最佳參數(shù)范圍,以便為緊湊快速強(qiáng)換熱器設(shè)計提供參考。圖17給出了單排微小尺度圓柱管束中的流動換熱規(guī)律,其中橫坐標(biāo)為進(jìn)口Re,圖17(a)和圖17(b)的縱坐標(biāo)分別為總壓恢復(fù)系數(shù)ζ和對流換熱系數(shù)h。從圖中可以看出,隨著進(jìn)口Re的增加,單排和多排管束的總壓恢復(fù)系數(shù)以及壁面平均對流換熱系數(shù)總體變化趨勢相同,即總壓恢復(fù)系數(shù)逐漸降低,對流換熱系數(shù)逐漸升高。在單排管束時,總壓恢復(fù)系數(shù)的變化趨勢基本呈線性規(guī)律,而在多排管束情況(見圖18)下,隨著Re的增加,總壓恢復(fù)系數(shù)的下降趨勢逐漸加劇。在單排管束和多排管束情況下,壁面平均對流換熱系數(shù)隨Re增加而增大,但增大趨勢逐漸變緩。這就意味著在選取進(jìn)口Re時,應(yīng)根據(jù)設(shè)計要求的流動和換熱性能選擇適中的范圍。

圖19和圖20分別給出了單排圓管在4個不同時刻的展向渦量與耗散熵產(chǎn)云圖及時均云圖,圖中Vorticity Z表示展向渦量,其中Z向為沿管長方向。由圖可知,圓柱后渦對的強(qiáng)剪切引起了很強(qiáng)的耗散熵產(chǎn),并遷移到下游較遠(yuǎn)的區(qū)域。相對而言,圓管外流動中換熱引起的熵產(chǎn)主要集中在圓管壁面附近(見圖21),其影響范圍相對集中,因此從減小熵產(chǎn)的角度來看,應(yīng)主要從減小耗散引起的熵產(chǎn)入手。在本文研究的情況下,管束中耗散引起的熵產(chǎn)主要包括兩個方面,即脈動耗散引起的熵產(chǎn)和平均流耗散引起的熵產(chǎn)。

圖19 單排圓管時均和瞬時展向渦量分布云圖Fig.19 Time-averaged and transient vorticity Z distributions of a single row of tubes

圖20 單排圓管時均和瞬時耗散熵產(chǎn)分布云圖Fig.20 Time-averaged and transient entropy generation due to energy dissipation for a single row of tubes

圖21 單排圓管時均和瞬時換熱熵產(chǎn)分布云圖Fig.21 Time-averaged and transient entropy generation due to heat transfer for a single row of tubes

圖22和圖23分別給出了多排圓管情況下展向渦量和耗散熵產(chǎn)的非定常時均云圖。與單排管束的時均解對比可以看出,在多排管束中,由于前一排管束尾跡對下游管束的作用,使得下游管束前緣處的耗散熵產(chǎn)明顯增大。同時,由于管束排列比較緊湊,管束尾跡受下游圓管固壁的限制,被擠壓在較窄的范圍內(nèi),并在此范圍內(nèi)產(chǎn)生了較大的熵產(chǎn)。

圖22 多排圓管情況下時均展向渦量分布Fig.22 Time-averaged vorticity Z distributions for tube matrix

圖23 多排圓管情況下時均耗散熵產(chǎn)分布Fig.23 Time-averaged entropy generation due to energy dissipation distributions for tube matrix

2.3 緊湊快速強(qiáng)換熱器設(shè)計方法發(fā)展

由于傳統(tǒng)的換熱器設(shè)計方法已不再適用于緊湊快速強(qiáng)換熱器設(shè)計,本文利用前述低維模型,發(fā)展了緊湊快速強(qiáng)換熱器的確定性設(shè)計方法,并利用該方法完成了用于高超聲速強(qiáng)預(yù)冷航空發(fā)動機(jī)不同位置的多種不同結(jié)構(gòu)形式的緊湊快速強(qiáng)換熱器確定性設(shè)計方案。其中,毛細(xì)管式換熱器能將高溫來流在流經(jīng)換熱器的瞬間冷卻至發(fā)動機(jī)可以正常工作的溫度,同時換熱器單位重量的換熱量超過100 k W/kg。圖24為毛細(xì)管式換熱器示意圖。

為保證設(shè)計方案滿足設(shè)計要求,采用數(shù)值模擬手段對設(shè)計結(jié)果進(jìn)行了檢驗。利用毛細(xì)管換熱器的簡化模型進(jìn)行數(shù)值模擬,模型中空氣流經(jīng)10排毛細(xì)管。為減小計算網(wǎng)格量,每排以一根毛細(xì)管進(jìn)行模擬,保證兩種工質(zhì)進(jìn)口流速與設(shè)計值相同,網(wǎng)格數(shù)總計780萬,網(wǎng)格示意圖如圖25所示。

圖24 毛細(xì)管式緊湊換熱器Fig.24 Micro-size tube compact precooler

圖25 毛細(xì)管換熱器數(shù)值模擬網(wǎng)格示意圖Fig.25 Meshes used in simulations of micro-size tube compact precooler

冷卻后空氣出口溫度的數(shù)值模擬結(jié)果與設(shè)計值對比如圖26所示,圖中橫坐標(biāo)為沿毛細(xì)管高度方向的相對位置,縱坐標(biāo)為空氣出口溫度。結(jié)果證明換熱器的設(shè)計方案初步滿足設(shè)計要求。

圖26 冷卻后空氣出口溫度數(shù)值模擬結(jié)果與設(shè)計值對比Fig.26 Comparison of cooled air temperature at outlet between computing result and design value

3 緊湊快速強(qiáng)換熱器先進(jìn)加工及檢測技術(shù)

圖27 不同溫度下GH4169晶粒尺寸效應(yīng)Fig.27 Grain size effect of GH4169 at different temperature

3.1 薄壁毛細(xì)管制造及檢測技術(shù)

薄壁毛細(xì)管的制造是緊湊快速換熱器的關(guān)鍵。其壁厚和均勻度嚴(yán)重影響到后續(xù)的換熱能力、換熱效率和結(jié)構(gòu)可靠性。此外,薄壁毛細(xì)管后續(xù)的定位裝配過程要求兩端端口具有很高的同心度,以保證釬料均勻分布在待焊接頭處。這些都對薄壁毛細(xì)管的制造提出了非?？量痰囊?為了成功制造薄壁毛細(xì)管,必須深入掌握其流動變形機(jī)理,包括微尺度變形中的晶粒尺寸效應(yīng)、幾何尺寸效應(yīng)關(guān)系,以及溫度對其尺寸效應(yīng)的影響等問題。

圖27分別給出了不同溫度下應(yīng)力隨晶粒尺寸的變化趨勢,以分析溫度對于GH4169晶粒尺寸效應(yīng)的影響,圖中橫坐標(biāo)為晶粒尺寸的倒數(shù),縱坐標(biāo)為應(yīng)力。圖例中,25、200、400分別代表變形溫度,0.002等數(shù)字代表真實應(yīng)變大小。

從圖中可以看出,在不同溫度、不同應(yīng)變下的晶粒尺寸效應(yīng)曲線中,都存在分段現(xiàn)象,且分段點隨著應(yīng)變的增加而變大；常溫下存在兩個分段點,第一個分段點隨著應(yīng)變增加變化趨勢明顯,第二個分段點和應(yīng)變沒有太大關(guān)系,這表明當(dāng)晶粒小到一定程度后,晶粒對應(yīng)力的影響可以忽略,應(yīng)力主要受應(yīng)變的影響,即晶粒尺寸效應(yīng)現(xiàn)象不明顯；當(dāng)在200℃ 的時候,分段點呈微微下凹的趨勢,在同一應(yīng)變下的線性趨勢最好,應(yīng)力水平主要受不同應(yīng)變的影響較大；在400℃下,分段點隨應(yīng)變增加而增大現(xiàn)象明顯。

將不同溫度下屈服時的晶粒尺寸效應(yīng)擬合線斜率放在一起,如圖28所示,圖中橫坐標(biāo)為晶粒尺寸,縱坐標(biāo)為屈服應(yīng)力?？梢赃M(jìn)一步看出明顯的分段現(xiàn)象,另外,存在一臨界晶粒尺寸,晶粒大于該尺寸時,斜率隨溫度下降明顯,且溫度越高下降速度越慢；晶粒小于該尺寸時,斜率隨溫度上升而下降,但下降速度基本不變。這就要求在選擇成形工藝路線時,需要考慮成形溫度、臨界晶粒尺寸等因素對于成形力的影響。

圖29給出了試制的毛細(xì)管樣品掃描電子顯微鏡的照片,通過SEM照片,可以看到其內(nèi)壁粗糙度滿足要求,而對于使用中影響最大的成形管壁內(nèi)部孔洞缺陷,在隨機(jī)選擇的多個樣品的SEM照片中均未發(fā)現(xiàn)成形孔洞缺陷,表明薄壁毛細(xì)管在拉拔過程中不存在成形孔洞缺陷問題。

表1則給出了3種不同熱處理狀態(tài)(軟態(tài)、半硬態(tài)和硬態(tài))的薄壁毛細(xì)管的單向拉伸力學(xué)性能檢測結(jié)果,結(jié)果表明不同狀態(tài)的抗拉強(qiáng)度不同,半硬態(tài)與硬態(tài)的薄壁毛細(xì)管抗拉強(qiáng)度相似,而軟態(tài)薄壁毛細(xì)管抗拉強(qiáng)度相對較低。圖30給出了半硬態(tài)成形下的單拉曲線圖,圖中橫坐標(biāo)為位移,縱坐標(biāo)為拉力。

圖28 不同溫度下屈服時的晶粒尺寸效應(yīng)擬合線斜率Fig.28 Different slopes fitted at various temperatures for grain size effect of yield stress

圖29 毛細(xì)管壁厚電鏡照片F(xiàn)ig.29 SEM photograph for wall thickness of micro-size tube

表1 不同成型方法毛細(xì)管力學(xué)性能Table 1 Mechanical properties of micro-size tube in different forming processes

圖30 半硬態(tài)單拉試驗結(jié)果Fig.30 Results of uniaxial-tension test in middle-hard state

3.2 薄壁毛細(xì)管的組合焊接

薄壁毛細(xì)管的組合焊接,直接關(guān)系到最終產(chǎn)品的質(zhì)量。在一般釬料的制備方法上,進(jìn)行適合工況的釬料研發(fā),掌握釬料成分對于焊接強(qiáng)度的作用規(guī)律；找出釬焊接頭的組織、性能隨溫度、釬縫間隙等焊接參數(shù)的演變規(guī)律；探明不同焊接熱循環(huán)過程中接頭內(nèi)部殘余內(nèi)應(yīng)力的分布規(guī)律；組合釬焊還要求掌握密集毛細(xì)管的定位裝配方法,探究多個焊接接頭間熱力學(xué)影響規(guī)律；最后還需制定大量薄壁毛細(xì)管釬焊接頭焊接質(zhì)量表征方法,合理評價焊接質(zhì)量。

本文研究中對現(xiàn)有的BNi-2釬料進(jìn)行改進(jìn),并通過釬焊時間的合理控制以保障釬焊質(zhì)量,工藝試驗表明:通常情況下,釬縫越大,越容易在釬縫中央部位形成連續(xù)金屬間化合物的脆性相,導(dǎo)致釬縫的脆性增加；這是由于釬料向母材內(nèi)部擴(kuò)散不均勻?qū)е碌?本研究中通過試驗得到合理的釬焊間隙。為了改善釬縫組織和提高接頭韌性,焊后進(jìn)一步在合適溫度下保溫一定時間。最終通過多次試驗確定合理的工藝參數(shù)。圖31給出了兩種不同形式焊接試驗件的照片,圖32則給出了焊接處的SEM照片,表明焊接接頭處的釬料溶解均勻,沒有凹凸不平等表面缺陷。

圖31 毛細(xì)管焊接試驗件Fig.31 A test sample of welded micro-size tube

圖32 焊接部位的釬料填充情況Fig.32 Filling result of solder at the welded position

4 結(jié) 論

國內(nèi)外的研究表明,強(qiáng)預(yù)冷高超聲速航空發(fā)動機(jī)技術(shù)是一項具有非常巨大的潛在技術(shù)優(yōu)勢和前瞻性的共性技術(shù),特別是REL公司的地面試驗已實現(xiàn)了強(qiáng)預(yù)冷的技術(shù)指標(biāo),值得引起關(guān)注。目前國內(nèi)也已開展了強(qiáng)預(yù)冷航空發(fā)動機(jī)相關(guān)研究,并在微尺度流熱耦合換熱機(jī)理、流熱耦合數(shù)值模擬、緊湊強(qiáng)換熱器設(shè)計制造等方面均取得一定進(jìn)展,為下一步的研究打下了良好的理論和技術(shù)基礎(chǔ)。

[1] Tang M,Hamilton B A,Chase R L.The quest for hypersonic flight with air-breathing propulsion,AIAA-2008-2546[R].Reston:AIAA,2008.

[2] Balepin V,Engers R,Terry S.MIPCC technology development,AIAA-2003-6929[R].Reston:AIAA,2003.

[3] Miyagi H,Miyagawa T,Monji T,et al.Combined cycle engine research in Japanese HYPR project,AIAA-1995-2751[R].Reston:AIAA,1995.

[4] Tanatsugu N,Sato T,Balepin V.Development study on ATREX engine[J].Acta Astronautica,1997,41(2-8):851-862.

[5] Harada K,Tanatsugu N,Sato T.Development study of a precooler for the air-turboramjet expander-cycle engine[J].Journal of Propulsion and Power,2001,17(6):1233-1238.

[6] Zhang D B.Analysis of Japanese precooling high-speed turbine engines[J].Aeronautical Information,2014,1599(4):1-4(in Chinese).張東寶.日本預(yù)冷卻高速渦輪發(fā)動機(jī)研究取得重要進(jìn)展[J].航空情報,2014,1599(4):1-4.

[7] Sato T,Kobayashi H,Tanatsugu N.Development study of the precooler of ATREX engine,AIAA-2003-6985[R].Reston:AIAA,2003.

[8] Vladimir B.High speed propulsion cycles,RTO-AVTVKI LS CSP-07-5052[R].Rhode Saint Genese:VKI,2007.

[9] Steelant J.Sustained hypersonic flight in Europe:technology drivers for LAPCAT II,AIAA-2009-7240[R].Reston:AIAA,2009.

[10] Steelant J.Sustained hypersonic flight in Europe:first technology achievements within LAPCAT II,AIAA-2011-2243[R].Reston:AIAA,2011.

[11] Zhang H J,Zou Z P,Li Y,et al.Preconditioned densitybased algorithm for conjugate porous/fluid/solid domains[J].Numerical Heat Transfer Part A,2011,60(2):129-153.

[12] Weiss J M,Smith W A.Preconditioning applied to variable and constant density flows[J].AIAA Journal,1995,33(11):2050-2057.

[13] Zhang B,Han J L.Multi-field coupled computing platform and thermal transfer of hypersonic thermal protection structures[J].Acta Aeronautica et Astronautica Sinica,2011,32(3):400-409(in Chinese).張兵,韓景龍.多場耦合計算平臺與高超聲速熱防護(hù)結(jié)構(gòu)傳熱問題研究[J].航空學(xué)報,2011,32(3):400-409.

[14] Li H Y,Leong K C,Jin L W,et al.Analysis of fluid flow and heat transfer in a channel with staggered porous blocks[J].International Journal of Thermal Sciences,2010,49(6):950-962.

[15] Nikitin N V,Pavellev A A.Turbulent flow in a channel with permeable walls,direct numerical simulation and results of three-parameter model[J].Fluid Dynamics,1998,33(6):826-832.

[16] Zhang H J.Investigation of numerical conjugate heat transfer method and coupling mechanism for hybrid porous/fluid/solid domains[D].Beijing:Beihang University,2013(in Chinese).張紅軍.多孔/流體/固體多區(qū)域流熱耦合數(shù)值模擬方法以及耦合機(jī)制研究[D].北京:北京航空航天大學(xué),2013.

[17] Li Y.A 3-D conjugate heat transfer solver and methodology research[D].Beijing:Beihang University,2011(in Chinese).李宇.三維流/熱耦合數(shù)值模擬程序的發(fā)展及方法研究[D].北京:北京航空航天大學(xué),2011.

[18] Zhang H J,Zou Z P,Song S H,et al.Investigations of heat transfer enhancement in a channel with staggered porous ribs by the preconditioned density-based algorithm[J].Numerical Heat Transfer Part A,2015,67(2):1370-1385.

[19] Zhang H J,Zou Z P.Investigation of a confined laminar impinging jet on a plate with a porous layer using the preconditioned density-based algorithm[J].Numerical Heat Transfer Part A,2012,61(4):241-267.

[20] Zhang H J,Zou Z P.Numerical investigation of laminarplate transpiration cooling by preconditioned density-based algorithm[J].Numerical Heat Transfer Part A,2012,62(10):761-779.