王成鵬， 楊永陽， 程克明

（南京航空航天大學(xué) 空氣動(dòng)力學(xué)系，江蘇 南京210016）

0 引 言

近年來各種高超聲速飛行實(shí)驗(yàn)的成功［1－2］，表明了高空條件下超聲速燃燒的可實(shí)現(xiàn)性，也驗(yàn)證了采用超燃沖壓發(fā)動(dòng)機(jī)方案推進(jìn)吸氣式高超聲速飛行器的可行性；盡管還存在如推力富裕度很小等很多關(guān)鍵問題有待解決［3］，但吸氣式高超聲速飛行器相對(duì)火箭在長(zhǎng)航時(shí)飛行時(shí)所表現(xiàn)出的經(jīng)濟(jì)性仍吸引了更多的關(guān)注。這種高超聲速飛行器機(jī)體／推進(jìn)系統(tǒng)高度一體化和飛行速域?qū)挼奶攸c(diǎn)要求在概念設(shè)計(jì)階段各學(xué)科的模型分析精度就要達(dá)到較高的水平，而總體選型和性能分析需要進(jìn)行大量的方案評(píng)估，需要快速地得到航程、有效載荷等各項(xiàng)性能的優(yōu)化參數(shù)，所以需要各部件建立具有一定精度的快速分析流動(dòng)模型［4］；其中針對(duì)燃燒室內(nèi)超聲速燃燒流場(chǎng)的分析和計(jì)算最為關(guān)鍵，超聲速燃燒室流場(chǎng)結(jié)構(gòu)復(fù)雜，存在著激波系、膨脹波系、激波附面層相互干擾結(jié)構(gòu)等，化學(xué)非平衡效應(yīng)強(qiáng)烈；針對(duì)這種含化學(xué)反應(yīng)的多組分高溫高速三維流動(dòng)，詳細(xì)計(jì)算需要采用CFD技術(shù)，需要考慮大量的中間化學(xué)反應(yīng)過程，一般采用簡(jiǎn)化的化學(xué)動(dòng)力學(xué)模型進(jìn)行并行模擬，二維和三維流場(chǎng)計(jì)算量均相當(dāng)大［5］，而且這種計(jì)算中包含的先進(jìn)湍流燃燒模型尚不成熟。尋找一種快捷簡(jiǎn)單而又具有一定精度的超聲速流動(dòng)燃燒一維模型為吸氣式高超飛行器的初步設(shè)計(jì)和性能估算服務(wù)非常必要。

國(guó)內(nèi)外在描述超聲速燃燒主要物理過程的模型建立方面做了很多工作［6－14］，文獻(xiàn)［6］對(duì)一類超燃燃燒室一維流動(dòng)分析模型進(jìn)行了總結(jié)并給出了改進(jìn)模型，模型主要以 William H.Heiser和 David T.Pratt［7］提出的三段物理過程（絕熱壓縮，等壓釋熱和膨脹釋熱）為框架，基于實(shí)驗(yàn)中獲得的實(shí)驗(yàn)數(shù)據(jù)分布，人為確定三種物理過程的分界點(diǎn)，忽略摩擦、熱交換等次要因素，先求馬赫數(shù)分布，再求其他變量；這類以實(shí)驗(yàn)燃燒室靜壓數(shù)據(jù)為基礎(chǔ)的超燃燃燒室流場(chǎng)分析模型可以給出超燃沖壓發(fā)動(dòng)機(jī)燃燒室內(nèi)的超聲速燃燒過程的合理描述，在估算燃燒室性能方面有重要的參考價(jià)值［8］。在概念性飛行器設(shè)計(jì)的初步階段，可能無法獲得相應(yīng)的燃燒室壁面壓力數(shù)據(jù)供預(yù)測(cè)飛行器性能，此時(shí)研究者們發(fā)展了一類準(zhǔn)一維或稱為廣義一維數(shù)值分析方法，結(jié)合有限速率化學(xué)反應(yīng)模型［9－10］或無限反應(yīng)速率一步模型［11－14］，在短時(shí)間內(nèi)獲得燃燒室壓力、馬赫數(shù)等參數(shù)沿程分布，為飛行器設(shè)計(jì)提供性能分析或評(píng)價(jià)。這其中Timothy，Starkey等發(fā)展的一種“準(zhǔn)一維超燃流場(chǎng)模型”［9－10］較有特色，該模型在考慮面積變化、壁面摩擦、燃料噴射以及向燃燒室內(nèi)壁傳熱影響的基礎(chǔ)上，在一維管流控制方程中耦合了有限速率化學(xué)反應(yīng)方程，加入了燃料混合模型，相對(duì)無限反應(yīng)速率模型在一定程度上提高了精度，該模型已經(jīng)被用于吸氣式高超聲速飛行器的初步設(shè)計(jì)之中［9］，他們建立起來的準(zhǔn)一維流動(dòng)模型能夠預(yù)測(cè)出著火點(diǎn)，在著火點(diǎn)之后能夠準(zhǔn)確地計(jì)算出壁面沿程靜壓分布；但是由于此模型沒有考慮噴射點(diǎn)之前燃燒高壓前傳形成的激波串結(jié)構(gòu)的影響，因此在燃燒室著火點(diǎn)之前，計(jì)算結(jié)果與實(shí)驗(yàn)值差別很大；結(jié)合超聲速燃燒物理過程目前的理解程度，分析這種模型計(jì)算結(jié)果與實(shí)驗(yàn)值的差別還需從實(shí)驗(yàn)觀察入手，潘余等［15］采用光學(xué)觀測(cè)的方式對(duì)不同著火方式下燃燒室內(nèi)著火過程進(jìn)行了研究，給出了著火過程中燃燒室內(nèi)的流動(dòng)現(xiàn)象，他們發(fā)現(xiàn)在一些射流自燃工況中，氫氣首先在射流噴孔下游某個(gè)區(qū)域自燃，初始火焰形成后，燃燒區(qū)放熱使局部壓力升高，在燃燒區(qū)域前形成若干激波系，最終實(shí)現(xiàn)穩(wěn)定燃燒。Schmidt V 等［16－17］應(yīng)用高速攝影和高速陰影對(duì)著火過程進(jìn)行了觀測(cè)，拍攝了火焰、流場(chǎng)陰影圖像等，結(jié)果表明火焰首先起始于燃燒室下游，然后逆流傳播至噴射點(diǎn)上游，最后建立穩(wěn)定燃燒的狀態(tài)。對(duì)這些實(shí)驗(yàn)現(xiàn)象的觀察也許能夠解釋Timothy，Starkey等［9－10］人的“準(zhǔn)一維超燃流場(chǎng)模型”在激波串區(qū)域失效的原因：實(shí)驗(yàn)中的超聲速燃燒過程可以由兩個(gè)階段描述，第一個(gè)階段，燃料從噴射口噴出，與主流摻混，在著火點(diǎn)自燃，開始形成初始火焰，壓力升高；第二個(gè)階段，燃燒造成的高壓在附面層內(nèi)亞聲速區(qū)域前傳，產(chǎn)生的逆壓力梯度使壁面附面層發(fā)生分離，進(jìn)而導(dǎo)致超聲速芯流面積收縮，壓縮主流，形成激波串結(jié)構(gòu)和穩(wěn)定的燃燒過程。由于 “準(zhǔn)一維超燃流場(chǎng)模型”是一種一維模型，從燃燒室進(jìn)口沿程向出口計(jì)算，模擬了燃燒的第一個(gè)階段，卻無法模擬第二個(gè)階段，造成了在著火點(diǎn)之前的壓力計(jì)算不準(zhǔn)。

從較早的 William H.Heiser和David T.Pratt［7］的“一維流動(dòng)分析模型”到最近的 Timothy，Starkey［9－10］的“準(zhǔn)一維超燃流場(chǎng)模型”，隨著對(duì)超燃現(xiàn)象研究的深入，模型開始更細(xì)致地考慮基元化學(xué)反應(yīng)的影響，同時(shí)也考慮了摩擦、熱交換、燃料噴射等一些在過去模型中忽略的因素，有研究［18］表明這些被忽略的因素很可能會(huì)對(duì)參數(shù)計(jì)算造成較大的誤差。即使這樣，Timothy，Starkey［9－10］的“準(zhǔn)一維超燃流場(chǎng)模型”仍然不能在燃料噴口之前的預(yù)燃激波串區(qū)域準(zhǔn)確地計(jì)算出沿程壓力分布，究其原因，是因?yàn)樵诮⒛Ｐ偷倪^程中沒有考慮激波串結(jié)構(gòu)的影響；而實(shí)驗(yàn)中超燃現(xiàn)象表明在燃料著火之后，高壓在附面層亞聲速流中前傳形成激波串，這一現(xiàn)象在“準(zhǔn)一維超燃流場(chǎng)模型”中沒有體現(xiàn)?；谶@種分析，本文在超聲速燃燒實(shí)驗(yàn)現(xiàn)象觀察的基礎(chǔ)上，在“準(zhǔn)一維超燃流場(chǎng)模型”中耦合進(jìn)激波串結(jié)構(gòu)模型［19－20］，提高其壁面壓力的預(yù)測(cè)能力，而壁面壓力計(jì)算結(jié)果精度的提高將對(duì)超燃燃燒室壓力阻／推力的計(jì)算乃至整個(gè)飛行器的推阻性能計(jì)算有重要意義。

1 物理模型與分析

物理模型的分析仍然借助于William H.Heiser和David T.Pratt［7］提出的超聲速燃燒模型，燃燒室構(gòu)型為常見的“等直段＋擴(kuò)張段”（圖1），隔離段進(jìn)口為經(jīng)過高超聲速進(jìn)氣道壓縮后的超聲速來流，燃料在3處噴入超聲速流場(chǎng)中，燃料與來流混合后在噴射點(diǎn)下游某點(diǎn)開始點(diǎn)燃，燃燒放熱形成的高壓壓縮上游超聲速來流，引起附面層分離，使得火焰區(qū)域得以逆流傳播，在d點(diǎn)和s點(diǎn)之間燃燒釋熱與超聲速有效流動(dòng)面積變化共同作用形成等壓區(qū)域；在s點(diǎn)之后面積擴(kuò)張的影響占主導(dǎo)，壓力沿程下降；在d點(diǎn)之前，一段激波串區(qū)域跨越燃料噴射點(diǎn)形成絕熱壓縮段銜接上游低壓與燃燒區(qū)高壓。下文將詳細(xì)敘述各個(gè)物理化學(xué)過程的?；幚矸椒ā?/p>

圖1 超聲速燃燒流動(dòng)物理模型Fig.1 Flow analysis model for supersonic combustor

2 超聲速燃燒流場(chǎng)控制方程

將燃燒室內(nèi)流動(dòng)介質(zhì)視為化學(xué)反應(yīng)完全氣體混合物，每一組元視為熱完全氣體，各項(xiàng)組元比熱與比熱比為溫度的函數(shù)；考慮流動(dòng)變量在x方向上的變化；多組分準(zhǔn)一維流剛性常微分控制方程組由連續(xù)方程（1）、動(dòng)量方程（2）、狀態(tài)方程（3）、組分質(zhì)量分?jǐn)?shù)方程（4）、能量方程（5）組成［9］：

公式（1）～（5）中A為燃燒室管道截面積（m2）；Cf為壁面摩擦系數(shù)；cp為比定壓熱容（J／（kg·K））；為以cp為自變量的函數(shù)，詳見文獻(xiàn)［9］；D為圓形截面管道直徑或非圓形截面管道當(dāng)量直徑（m）；h為單位質(zhì)量的焓（J／kg）；Ma為馬赫數(shù)；MW為分子量（kg／kmol為組分混合物分子量（kg／kmol）；˙m為質(zhì)量流量（kg／s）；p為壓力（Pa）；Pr為普朗特?cái)?shù)；T為溫度（K）；U為矢量速度（m／s）；x為軸向坐標(biāo)（m）；Y為質(zhì)量分?jǐn)?shù)；ε為燃料射流速度在x向上分量與主流速度之比；γ為比熱比；ρ為密度（kg／m3）；˙ωi為每一組元在多步反應(yīng)機(jī)理中的凈生成率（kmol／（s·m3））：

其中、對(duì)應(yīng)i組分在j個(gè)反應(yīng)中方程兩邊反應(yīng)物與生成物的化學(xué)當(dāng)量系數(shù)，kf，jkr，j為正向和逆向反應(yīng)速率系數(shù)，L表示基元反應(yīng)步數(shù)，N為組分個(gè)數(shù)，這種表達(dá)式非常便于化學(xué)動(dòng)力學(xué)過程的求解。式（4、5）中下標(biāo)added為燃料噴射添加的組元，aw為絕熱壁條件，ω為壁面條件，0為滯止?fàn)顟B(tài)。式（6）中的化學(xué)反應(yīng)機(jī)理在本文中采用了Jachimowski提出的13組分、32反應(yīng)的氫氣／空氣反應(yīng)模型［21］，該模型是專門針對(duì)高速燃燒問題設(shè)計(jì)的反應(yīng)模型，原模型中包括了氮?dú)獾幕瘜W(xué)反應(yīng)，考慮燃燒室溫度小于氮?dú)夥纸鉁囟?000K，所以僅考慮氫氧反應(yīng)過程。化學(xué)反應(yīng)生成項(xiàng)及其涉及的反應(yīng)特征時(shí)間尺度差異很大，使得上述控制方程組是剛性常微分方程組，采用剛性系統(tǒng)算法器VODE［22］求解，方程組（1～5）的書寫格式也是為了采用 VODE 求解。CHEMKIN－II［23］作為一種子程序庫(kù)為流動(dòng)燃燒耦合過程中的化學(xué)動(dòng)力學(xué)問題提供了計(jì)算工具，計(jì)算之前首先將氫氣／空氣反應(yīng)化學(xué)機(jī)理讀入CHEMKIN－II的解釋器（Interpreter），然后從熱力學(xué)數(shù)據(jù)庫(kù)中提出與機(jī)理相關(guān)的熱力學(xué)信息，、hi、MW、cp等各項(xiàng)參數(shù)由 CHEMKIN－II的 子程序庫(kù)計(jì)算得到，輸出的計(jì)算結(jié)果應(yīng)用于控制方程組（1～5）的相關(guān)項(xiàng)中。

針對(duì)燃料垂直地還是平行地與主流混合，將采取不同的混合模型。垂直噴注依靠燃料射流的能力獲得混合效果，文獻(xiàn)［9］中給出了一種燃料垂直噴射的混合模型，假定在一維范疇內(nèi)，燃料的添加過程在一定長(zhǎng)度內(nèi)發(fā)生，從燃料噴射點(diǎn)向下游定義一混合長(zhǎng)度Lmix，在混合長(zhǎng)度內(nèi)燃料僅僅發(fā)生了混合，燃料生產(chǎn)率˙ωi被設(shè)定為0，當(dāng)達(dá)到混合長(zhǎng)度后，反應(yīng)開始發(fā)生，直至燃料到了燃燒室的末端或反應(yīng)已達(dá)平衡。本文針對(duì)燃料垂直噴射工況仍然采用這種模型，不同于文獻(xiàn)［9］中的人為給定Lmix，基于文獻(xiàn)［24］中相似理論分析和實(shí)驗(yàn)數(shù)據(jù)結(jié)果，采用下式來確定燃料空氣混合長(zhǎng)度Lmix：

其中d為噴孔直徑（m），j、∞分別表示噴流與主流條件，ξc為噴流中心線上燃料濃度，這里用燃料摩爾分?jǐn)?shù)表示，文獻(xiàn)［24］采用式（7）來表達(dá)噴流中心線上燃料濃度在噴孔下游沿流向方向的變化趨勢(shì)，本文將可燃上限對(duì)應(yīng)的氫燃料摩爾分?jǐn)?shù)代入ξc，來確定Lmix。對(duì)于采用支板與主流平行噴射燃料的工況，可壓縮混合長(zhǎng)度由下式［13］來計(jì)算：

Mc為燃料與主流混合層對(duì)流馬赫數(shù)，c為聲速。冷流混合實(shí)驗(yàn)結(jié)果與多維計(jì)算結(jié)果表明公式（7、8）僅適用各自垂直和平行主流噴射燃料的工況，如果實(shí)際的燃燒室模型采用與主流成一定角度噴射燃料或采用凹腔等穩(wěn)焰結(jié)構(gòu)，均需根據(jù)實(shí)際情況建立混合模型。

3 激波串模型

如前所述，燃料噴入超聲速來流發(fā)生混合、燃燒，形成的高壓將通過附面層向上游傳播，在隔離段中形成激波串結(jié)構(gòu)（圖1中u－d段）。第2節(jié)中的常微分控制方程組在隔離段進(jìn)口賦初值求解，考慮燃料噴射、混合、燃燒模型，卻沒有考慮反壓前傳形成的激波串對(duì)沿程參數(shù)的影響，所以利用控制方程組求出沿程壓力分布后，還需建立激波串模型。激波串為隔離段進(jìn)口超聲速來流在壁面摩擦、燃燒反壓造成的逆壓力梯度共同作用下減速擴(kuò)壓形成的復(fù)雜流動(dòng)結(jié)構(gòu)，包含激波相交反射、激波／附面層的相互作用、氣流分離等，即使采用多維RANS（Reynolds－averaged Navier－Stokes）數(shù)值模擬手段刻畫激波串結(jié)構(gòu)目前也是一大難題，原因在于湍流模型在處理分離流、可壓縮流體、激波不穩(wěn)定性效應(yīng)等現(xiàn)象時(shí)精度過低，而目前研究者對(duì)這些在激波串流動(dòng)結(jié)構(gòu)中包含的流動(dòng)現(xiàn)象的物理機(jī)制了解的很不深入。自1958年Crocco［25］首次提出“無激波理論”激波串理論模型以來，至今基于理論分析和實(shí)驗(yàn)研究已發(fā)展了多個(gè)激波串（或偽激波）模型，最具代表性的包括 Matsuo［26］發(fā)展的“流量平均”理論模型和Waltrup［19］的基于實(shí)驗(yàn)數(shù)據(jù)的激波串公式，Matsuo的“流量平均”模型采用流量平均的方法處理激波串之前的包含附面層的截面平均參數(shù)，控制體包括整個(gè)偽激波結(jié)構(gòu)（激波串＋下游的混合段），下游的靜壓、總壓恢復(fù)等參數(shù)可以通過上游的流量平均參數(shù)求出，模型計(jì)算值經(jīng)大量實(shí)驗(yàn)值驗(yàn)證表明具有很高的精度，但該模型無法給出激波串區(qū)域的沿程參數(shù)變化曲線；激波串區(qū)域流動(dòng)分離導(dǎo)致的超聲速流動(dòng)有效面積無法精確確定，導(dǎo)致在一維范疇上處理激波串問題還需借助實(shí)驗(yàn)分析或經(jīng)驗(yàn)方法，本文采用Waltrup的激波串公式來處理激波串區(qū)域的壓力分布，該公式基于大量一維軸對(duì)稱等直圓管內(nèi)激波串結(jié)構(gòu)實(shí)驗(yàn)數(shù)據(jù)，考慮了馬赫數(shù)、進(jìn)口附面層動(dòng)量厚度以及基于進(jìn)口附面層動(dòng)量厚度的雷諾數(shù)等因素的影響，并根據(jù)影響的程度確定了影響因子，確定了管內(nèi)增壓比與激波串長(zhǎng)度的關(guān)系式，因其實(shí)用性而被后人在計(jì)算激波串區(qū)域參數(shù)時(shí)廣泛引用：

式中D為圓管直徑，Ma1為激波串前主流馬赫數(shù)，θ1為激波串前附面層動(dòng)量損失厚度，Reθ為基于激波串前附面層動(dòng)量損失厚度的雷諾數(shù)，p1為激波串前壁面靜壓，p為對(duì)應(yīng)激波串中軸向位置x點(diǎn)的壁面靜壓。如果管道為矩形截面，則采用下式計(jì)算激波串區(qū)域壓力分布：

將公式（9）中的改為，圓管直徑D改為管道進(jìn)口高度H。如果實(shí)際工程應(yīng)用中需要考慮管道進(jìn)口（進(jìn)氣道出口）非對(duì)稱附面層來流工況，可以采用下式［20］計(jì)算：

式中Dθ為基于附面層動(dòng)量損失厚度的表征進(jìn)氣道非對(duì)稱出口條件的非對(duì)稱度因子，θmax、θmin分別為管道進(jìn)口最大與最小附面層動(dòng)量損失厚度，根據(jù)實(shí)驗(yàn)數(shù)據(jù)取為0.3。

由于燃料的噴射、摻混、著火延遲，高壓在燃料噴射點(diǎn)之后形成，圖1表明反壓形成的激波串結(jié)構(gòu)跨越噴射點(diǎn)被分成兩段（u－3段和3－d段），為了確定激波串的位置，文獻(xiàn)［27］發(fā)現(xiàn)3－d段軸向長(zhǎng)度x3d與上游來流附面層動(dòng)量損失厚度θ1、燃燒生成的高壓ps／p1有關(guān)：

在圖1的d點(diǎn)和s點(diǎn)之間，考慮超聲速有效流動(dòng)面積擴(kuò)張與燃燒加熱的共同作用，假定為等壓段。

總結(jié)整個(gè)流場(chǎng)分析模型的計(jì)算過程：1）求解第二節(jié)中的超聲速燃燒流場(chǎng)控制方程組，確定燃燒著火點(diǎn)和燃燒生成的高壓ps；2）由公式（12）確定激波串在燃料噴射口下游的長(zhǎng)度和圖1中d點(diǎn)的位置；3）由ps／p1根據(jù)管道形狀和來流工況確定采用公式（9～11）中的一個(gè)確定激波串的總長(zhǎng)度xud和ud段壁面壓力的分布；得到整個(gè)燃燒室沿程靜壓分布后，可以根據(jù)需要求解相應(yīng)控制方程得到其他流動(dòng)變量變化曲線。

4 模型算例驗(yàn)證

首先選擇了一種常見的“等直段＋擴(kuò)張段”超燃燃燒室基本構(gòu)型進(jìn)行流場(chǎng)計(jì)算，算例為Billg和Grenleski［28］研究的一個(gè)軸對(duì)稱的、壁面垂直噴射燃料的超燃燃燒室模型，如圖2所示，燃燒室的組成：一段長(zhǎng)0.3312m的截面積為0.0038m2的圓柱，之后緊接著的是長(zhǎng)0.6096m的發(fā)散型圓臺(tái)，圓臺(tái)進(jìn)出口面積比為2；在x＝0處沿周向均勻分布8個(gè)直徑為2.64mm的燃料噴射口。根據(jù)實(shí)驗(yàn)條件，確定計(jì)算模型與數(shù)值模擬的來流條件：Ma＝3.19，p＝54000Pa，T＝920K，壁溫Tω＝500K；假定來流空氣的體積組成是氮78%，氧21%，氬1%；氫燃料溫度705K，油氣當(dāng)量比為0.5，噴流馬赫數(shù)為1。

圖2 Billig超燃實(shí)驗(yàn)?zāi)Ｐ虵ig.2 Billig's experimental combustor configuration

圖3中的虛線為求解第二節(jié)中的超聲速燃燒流場(chǎng)控制方程組所得壓力沿程分布，計(jì)算過程采用了公式（7）所確定的氫垂直噴射混合模型，燃料在x＝0處開始噴射入主流，經(jīng)過混合燃燒過程在x＝0.078m處達(dá)到最大壓力點(diǎn)，而后在擴(kuò)張段壓力沿程下降，可見這段計(jì)算過程可以確定燃料的著火點(diǎn)和很好的預(yù)測(cè)擴(kuò)張段壓力分布［9］，圖4給出了對(duì)應(yīng)圖3中虛線壓力數(shù)據(jù)的溫度和組分沿程分布，表明了對(duì)應(yīng)過程中H2O的生成與OH基的分布。考慮燃燒高壓在附面層內(nèi)亞聲速區(qū)域的前傳過程，采用第三節(jié)中的激波串模型來處理氫燃料著火點(diǎn)之前的壁面壓力分布，其中激波串長(zhǎng)度采用公式（9）計(jì)算，結(jié)果如圖3實(shí)線所示，可見考慮激波串的模型很好的計(jì)算了壁面壓力分布，考慮該模型的較高精度和快速計(jì)算特征，可以用于燃燒室的概念設(shè)計(jì)和性能優(yōu)化中。為了更好地理解燃燒室內(nèi)流動(dòng)特征，圖5給出了采用文獻(xiàn)［5］中的三維數(shù)值模擬方法計(jì)算得到的燃燒室中心面上的等馬赫數(shù)線圖譜，可見噴口附近形成的激波串結(jié)構(gòu)。

圖3 不同模型壁面壓力計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比Fig.3 Comparison of theoretical and experimental wall pressure distributions

圖4 未考慮激波串模型時(shí)沿程溫度與組分計(jì)算結(jié)果Fig.4 Temperature and mass fraction of H2，H2O，and OH distributions calculated using the model without shock train

圖5 數(shù)值模擬Billig燃燒室等馬赫數(shù)線圖譜Fig.5 Mach number contours in Billig′s experimental combustor

在各種形式的超燃燃燒室數(shù)值和實(shí)驗(yàn)研究中，帶后向臺(tái)階的垂直噴氫方式也非常常見，第二個(gè)算例選擇了Billg和Grenleski［28］的一個(gè)軸對(duì)稱的、帶后向臺(tái)階的超燃燃燒室模型，燃燒室?guī)缀螚l件見圖6，燃燒室同樣由“等直段＋擴(kuò)張段”組成，與第一個(gè)算例不同的是，等直段由兩段不同直徑的部分組成，擴(kuò)張段也由不同擴(kuò)張角的兩段發(fā)散型圓臺(tái)連接而成，氫氣在后向臺(tái)階附近上游垂直噴射。第二個(gè)算例與第一個(gè)算例采用同一個(gè)超聲速噴管，計(jì)算模型來流條件同第一個(gè)算例，所不同的是氫當(dāng)量比提高至0.93。圖7給出模型計(jì)算得出的壁面壓力分布與實(shí)驗(yàn)數(shù)據(jù)的對(duì)比，考慮激波串的模型計(jì)算所得激波串長(zhǎng)度為0.0542m，沿主流方向于噴口之后的延伸距離為0.0051m，盡管存在臺(tái)階，模型在x＝0附近仍很好地計(jì)算了激波串區(qū)域的壓力分布，數(shù)值模擬結(jié)果也表明臺(tái)階的存在對(duì)激波串結(jié)構(gòu)靠近壁面的分離區(qū)流場(chǎng)結(jié)構(gòu)影響較?。辉趚＝0.1～0.2m處計(jì)算結(jié)果的偏低表明等壓假設(shè)的不完善，但作為一種高效快捷的模型計(jì)算結(jié)果可以接受。

圖6 帶后向臺(tái)階的超燃燃燒室模型Fig.6 A backward－facing step combustor configuration

圖7 后向臺(tái)階模型壁面壓力計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比Fig.7 Comparison of theoretical and experimental wall pressure distributions is a step－cylinder－cone combustor

作為一種重要的燃料噴注方式，支板噴射可以將燃料噴入超聲速來流的主流當(dāng)中并獲得良好的混合效果。第三個(gè)算例選擇了Anderson和Gooderum［29］的支板噴射氫燃料超燃試驗(yàn)，如圖8所示，燃燒室由“等直段＋擴(kuò)張段”段構(gòu)成，氫燃料由支板上噴口垂直噴入主流，燃料噴口位于圖8中x＝0處，燃燒室截面為矩形，沿程保持等寬度0.17m，其他尺寸見圖中標(biāo)識(shí)。噴管出口條件為馬赫數(shù)Ma＝2.7，靜壓P＝75777Pa，靜溫T＝1031K，油氣當(dāng)量比為0.619，壁面溫度假定為Tw＝500K，來流空氣組分按照質(zhì)量百分?jǐn)?shù)定為：O2：27.69%，N2：50.52%，H2O：21.79%。圖9中考慮激波串的模型計(jì)算所得激波串長(zhǎng)度為0.1562m，在燃燒室膨脹段計(jì)算壓力分布與實(shí)驗(yàn)數(shù)據(jù)吻合很好，在x＝0～0.2m區(qū)間模型計(jì)算結(jié)果沒有反映出實(shí)驗(yàn)數(shù)據(jù)的跳躍趨勢(shì)，是由于模型沒有考慮支板與壁面間的激波反射相交現(xiàn)象。

圖8 Anderson超燃實(shí)驗(yàn)?zāi)Ｐ虵ig.8 Anderson′s experimental combustor configuration

圖9 Anderson實(shí)驗(yàn)?zāi)Ｐ捅诿鎵毫τ?jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比Fig.9 Comparison of theoretical and experimental wall presuure distributions in Anderson′s experimental combustor

通過對(duì)上述3例的比較分析，發(fā)現(xiàn)本文所建立模型能很好地反映超聲速燃燒室流動(dòng)中的物理化學(xué)過程，并和實(shí)驗(yàn)數(shù)據(jù)相符。從而認(rèn)為該模型是有效的。

5 結(jié) 論

對(duì)考慮了激波串結(jié)構(gòu)模型的超聲速燃燒流場(chǎng)中物理化學(xué)過程的?；幚矸椒ㄟM(jìn)行了研究，建立了意義明確、精度較高的流場(chǎng)分析模型，計(jì)算了“等直段＋擴(kuò)張段”、后向臺(tái)階和支板噴射三種構(gòu)型的燃燒室內(nèi)氫燃料的超聲速燃燒問題，在每一算例中，均能夠計(jì)算出與實(shí)驗(yàn)數(shù)據(jù)分布趨勢(shì)一致的壁面沿程壓力數(shù)據(jù)，該模型適用于采用超燃沖壓發(fā)動(dòng)機(jī)推進(jìn)的高超聲速飛行器概念設(shè)計(jì)階段的流場(chǎng)性能預(yù)測(cè)分析和參數(shù)優(yōu)化，也可用于地面、飛行實(shí)驗(yàn)數(shù)據(jù)分析，加深對(duì)燃燒室內(nèi)氣動(dòng)熱力過程的理解。論文給出了處理超燃燃燒室內(nèi)流動(dòng)問題的框架，所述方法具有一定的參考價(jià)值，研究人員可以根據(jù)具體的燃燒室構(gòu)型（如凹腔結(jié)構(gòu)、多點(diǎn)噴射等）在本文的基礎(chǔ)上構(gòu)建適合的混合、燃燒、激波串模型對(duì)流場(chǎng)進(jìn)行分析。

［1］PEEBLES C.Road to Mach 10：lessons learned from the X－43Aflight research program［R］.AIAA Inc.，Reston，VA，2007.

［2］SMART M K，HASS N E，PAULL A.Flight data analysis of the HyShot 2scramjet flight experiment［J］.AIAAJournal，2006，44（10）：2375－2366.

［3］姜宗林.關(guān)于吸氣式高超聲速推進(jìn)技術(shù)研究的思考［J］.力學(xué)進(jìn)展，2006，39（4）：398－405.

［4］羅世彬.高超聲速飛行器機(jī)體／發(fā)動(dòng)機(jī)一體化及總體多學(xué)科設(shè)計(jì)優(yōu)化方法研究［D］.［博士論文］.長(zhǎng)沙：國(guó)防科技大學(xué)，2004.

［5］王成鵬，楊永陽，劉晨，等.超燃燃燒室流場(chǎng)計(jì)算方法比較分析［J］.航空動(dòng)力學(xué)報(bào)，2009，24（5）：963－969.

［6］張鵬，俞剛.超燃燃燒室一維流場(chǎng)分析模型的研究［J］.流體力學(xué)實(shí)驗(yàn)與測(cè)量，2003，17（1）：88－92.

［7］HEISER W H，PRATT D T.Hypersonic airbreathing propulsion［M］.washington DC.AIAA，1994.

［8］俞剛，李建國(guó).氫空氣超聲速燃燒研究［J］.流體力學(xué)實(shí)驗(yàn)與測(cè)量，1999，13（1）：1－12.

［9］O′BRIEN T F，STARKEY R P，LEWIS M J.Quasi－onedimensional high－speed engine model with finite－rate chemistry［J］.JournalofPropulsionandPower，2001，17（6）：1366－1374.

［10］STARKEY R P，LEWIS M J.Sensitivity of hydrocarbon combustion modeling for hypersonic missile Design［J］.JournalofPropulsionandPower，2003，19（1）：89－97.

［11］王蘭，邢建文，鄭忠華，等.超燃沖壓發(fā)動(dòng)機(jī)內(nèi)流性能的一維評(píng)估［J］.推進(jìn)技術(shù)，2008，29（6）：641－645.

［12］陳強(qiáng)，陳立紅，顧洪斌，等.釋熱分布對(duì)超燃沖壓發(fā)動(dòng)機(jī)性能的影響及優(yōu)化［J］.推進(jìn)技術(shù)，2009，30（2）：135－138.

［13］BIRZER C，DOOLAN C J.Quasi－one－dimensional model of hydrogen－fueled scramjet combustors［J］.Journalof PropulsionandPower，2009，25（6）：1220－1225.

［14］王振鋒，白菡塵，劉初平.燃燒效率一維評(píng)價(jià)的影響因素研究［J］.實(shí)驗(yàn)流體力學(xué)，2008，22（1）：11－16.

［15］潘余，劉衛(wèi)東，梁劍寒，等.模型超燃沖壓發(fā)動(dòng)機(jī)內(nèi)著火過程分析［J］.力學(xué)學(xué)報(bào)，2009，41（4）：455－462.

［16］SENDER J，SCHMIDT V，OSCHWALD M.Highly time resolved observation of the ignition process of a LOx／GH2coaxial injected spray［R］.AIAA 2001－3942，2001.

［17］SCHMIDT V，WEPLER U，HAIDN O J，et al.Characterization of the primary ignition process of a coaxial GH2／LOx spray［R］.AIAA 2004－1167，2004.

［18］ROGERS R C，SCHEXNAYDER C J.Chemical kinetics analysis of hydrogen－air ignition and reaction times［R］.NASA TP 1856，1981.

［19］WALTRUP P J，BILLIG F S.The structure of shock waves in cylindrical ducts［J］.AIAAJournal，1973，11（10）：1404－1408.

［20］王成鵬，張堃元，程克明.非對(duì)稱來流隔離段流動(dòng)特性研究［J］.推進(jìn)技術(shù)，2006，27（5）：436－440.

［21］JACHIMOWSKI C J.An analytical study of the hydrogen－air reaction with application to scramjet combustion［R］.NASA TP 2791，F(xiàn)eb.1988.

［22］BROWN P N，BYNE G D，HINDMARCH A C VODE.A variable－coefficient ordinary differential equation solver［J］.SIAMJournalonScientificandStatisticalComputing，1989，10（5）：1038－1051.

［23］KEE R J，RUPLEY F M，MILLER J A.Chemkin－Ⅱ：a fortran chemical kinetics package for the analysis of gas phase chemical kinetics［R］.SAND89－8009B，Sandia National Lab.，Albuquerque，NM，April 1989.

［24］HASSELBRINK E F，MUNGAL M G.Transverse jets and jet flames：Part 1.Scaling laws for strong transverse jets［J］.J.FluidMech.，2001，443：1－25.

［25］CROCCO L.One－dimensional treatment of steady gas dynamics［M］.in：emmons HW，editor.Fundamentals of gasdynamics，Princeton：Princeton University Press，1958：110－130.

［26］MATSUO K，MIYAZATO Y.Mass averaging pseudoshock model in a straight flow passage［J］.Proc.IMechE Part G：J.Aerospace Engineering，1999，213（6）：365－375.

［27］WALTRUP P J，BILLIG F S.Prediction of precombustion wall pressure distributions in scramjet engines［J］.J.Spacecraft，1973，10（9）：620－622.

［28］BILLIG F S，DUGGER G L，WALTRUP P J.Inletcombustor interface problems in scramjet engines［R］.The First International Symposium on Air Breathing Propulsion，Marseille，F(xiàn)rance，June 1972.

［29］ANDERSON G Y，GOODERUM P B.Exploratory tests of two strut fuel injectors for supersonic combustion［R］.NASA TN D－7581，1974.