環(huán)境壓力降低對(duì)底排二次燃燒影響的數(shù)值模擬

余文杰，余永剛

（南京理工大學(xué)能源與動(dòng)力工程學(xué)院，江蘇南京210094）

環(huán)境壓力降低對(duì)底排二次燃燒影響的數(shù)值模擬

余文杰，余永剛

（南京理工大學(xué)能源與動(dòng)力工程學(xué)院，江蘇南京210094）

為了揭示高空低壓環(huán)境下底排減阻率減小的機(jī)理,建立底排裝置尾部流場(chǎng)的化學(xué)非平衡流數(shù)學(xué)物理模型。其中二次燃燒模型采用10組分25步反應(yīng)的H2-CO燃燒模型,運(yùn)用統(tǒng)一算法的思路編程求解二維軸對(duì)稱方程組,對(duì)底排尾部流場(chǎng)進(jìn)行數(shù)值模擬。模擬結(jié)果和實(shí)驗(yàn)進(jìn)行對(duì)比驗(yàn)證,基本吻合。在此基礎(chǔ)上,對(duì)底排尾部流場(chǎng)以及燃燒特性進(jìn)行數(shù)值預(yù)測(cè),研究環(huán)境壓力降低對(duì)底排尾部二次燃燒的影響。結(jié)果表明：二次燃燒對(duì)底部加能的貢獻(xiàn)是熱排氣的6.4倍,是底排加能減阻的關(guān)鍵；隨著環(huán)境壓力的降低,模型尾部的環(huán)狀回流區(qū)內(nèi)H2的燃燒效率逐漸降低,中間產(chǎn)物H逐漸增多,燃燒逐漸變得不充分,導(dǎo)致底排減阻率明顯下降。

兵器科學(xué)與技術(shù)；底部排氣；二次燃燒；數(shù)值模擬；尾部流場(chǎng)；化學(xué)非平衡流

0　引言

彈丸在超聲速飛行時(shí)，會(huì)遇到較大的底部阻力。底排藥劑燃燒產(chǎn)生負(fù)氧型高溫氣體，排入彈丸尾部后與空氣中的氧氣接觸發(fā)生二次燃燒，燃燒釋放的能量使底部下游較長(zhǎng)的距離上保持高溫，使底部壓力顯著提高，底部阻力大大減?。?］。由于底排為亞聲速排氣，環(huán)境壓力的降低將影響底排裝置內(nèi)部藥劑的燃燒特性，進(jìn)而影響底排裝置外部的二次燃燒，并最終使底排減阻率減小。

文獻(xiàn)［2-3］對(duì)底部冷排氣的圓柱體模型的尾部流場(chǎng)進(jìn)行實(shí)驗(yàn)研究，并對(duì)底排尾部流場(chǎng)進(jìn)行了可視化。Bowman等［4］對(duì)底部熱排氣的圓柱體模型的底部壓力進(jìn)行實(shí)驗(yàn)研究，發(fā)現(xiàn)排氣溫度即使高達(dá)5 070 K，最佳底排減阻率也只有25%左右。丁則勝等［5］利用部分預(yù)混合可燃?xì)怏w為介質(zhì)，對(duì)底部燃燒進(jìn)行實(shí)驗(yàn)研究，能夠反映尾跡區(qū)中的二次燃燒效應(yīng)。發(fā)現(xiàn)底部燃燒時(shí)底排減阻率要比冷排氣時(shí)高出一個(gè)量級(jí)，最佳底排減阻率達(dá)到70%～80%.可見(jiàn)尾部的二次燃燒是底排增壓減阻的關(guān)鍵。丁則勝等［6］實(shí)驗(yàn)研究了環(huán)境壓力對(duì)底排性能的影響，發(fā)現(xiàn)隨著環(huán)境壓力下降，底排火焰溫度顯著降低。

數(shù)值模擬方面，Sahu等［7］最先運(yùn)用冷排氣的方法對(duì)底排尾部流場(chǎng)進(jìn)行數(shù)值研究，發(fā)現(xiàn)隨著排氣參數(shù)的變化底壓先增加再減小。Gibeling等［8］針對(duì)底排燃燒問(wèn)題提出一種12步化學(xué)反應(yīng)的H2-CO燃燒模型，對(duì)小排氣參數(shù)I=0.002 2時(shí)的底排尾部流場(chǎng)進(jìn)行了數(shù)值研究。Choir等［9］根據(jù)端羥基聚丁二烯（HTPB）和高氯酸銨（AP）的燃燒特性，建立氫氣和一氧化碳的燃燒模型，對(duì)復(fù)合型增程彈的全流場(chǎng)進(jìn)行了數(shù)值研究。陸中兵等［10］對(duì)底部排氣彈三維湍流流場(chǎng)進(jìn)行了數(shù)值模擬。陳新虹等［11］采用熱排氣的方法數(shù)值研究了排氣能量對(duì)底排彈氣動(dòng)特性的影響。Shin等分別采用大渦模擬［12］和直接模擬［13］的方法對(duì)底排尾部流場(chǎng)進(jìn)行了數(shù)值模擬，研究尾部回流區(qū)的大小和形狀。關(guān)于環(huán)境壓力對(duì)尾部二次燃燒影響的數(shù)值研究未見(jiàn)文獻(xiàn)報(bào)道。

以往研究表明：亞大氣壓下底排減阻率降低的本質(zhì)原因是向底部區(qū)域加能的減?。?］。所加的能量來(lái)自于高溫排氣的熱能和尾部的二次燃燒所釋放的熱量。為了揭示環(huán)境壓力對(duì)底排增壓減阻影響的機(jī)理，本文在和文獻(xiàn)［5］中的實(shí)驗(yàn)對(duì)比的基礎(chǔ)上，對(duì)底排裝置尾部的化學(xué)非平衡流進(jìn)行編程模擬，研究二次燃燒對(duì)加能的貢獻(xiàn)以及環(huán)境壓力對(duì)二次燃燒的影響。為進(jìn)一步研究如何減小環(huán)境壓力所帶來(lái)的負(fù)面影響提供參考依據(jù)。

1　控制方程

假設(shè)模型尾部流場(chǎng)軸對(duì)稱，微分形式的二維軸對(duì)稱控制方程組如下：

式中：U為守恒變量；F、G為對(duì)流項(xiàng)通量；Fv、Gv為擴(kuò)散項(xiàng)通量；W、Q、S分別為湍流源項(xiàng)、軸對(duì)稱源項(xiàng)、化學(xué)反應(yīng)源項(xiàng)；x、y分別表示軸向和徑向；u、v分別為軸向速度和徑向速度；ρ為密度；p為壓力；k為湍動(dòng)能；ω為湍動(dòng)能耗散率；τ為粘性應(yīng)力；ρi、Yi、Di、hi分別為i組分的密度、質(zhì)量分?jǐn)?shù)、擴(kuò)散系數(shù)、焓；λ、μl、μt分別為熱傳導(dǎo)系數(shù)、層流粘性系數(shù)、湍流粘性系數(shù)；qx、qy為導(dǎo)熱熱流，形式為；e為單位體積總能，形式為為i組分的化學(xué)反應(yīng)質(zhì)量產(chǎn)生率。湍流模型為SST模型，系數(shù)σk、σω和湍流源項(xiàng)Sk、Sω的給定參見(jiàn)文獻(xiàn)［14］。

控制方程組中的軸對(duì)稱源項(xiàng)Q為對(duì)流項(xiàng)的軸對(duì)稱源項(xiàng)。以往的研究［15-16］中發(fā)現(xiàn)對(duì)流項(xiàng)的軸對(duì)稱源項(xiàng)對(duì)控制方程的影響起主導(dǎo)作用，故忽略了粘性項(xiàng)的軸對(duì)稱源項(xiàng)。

2　模擬方法

2.1物理模型

如圖1所示，模型為圓柱體結(jié)構(gòu)（有船尾），x、y分別表示軸向和徑向，R為模型最大半徑，船尾長(zhǎng)與模型最大半徑值相同，船尾角為5°，Ma∞、p∞、T∞分別為來(lái)流馬赫數(shù)、來(lái)流靜壓、來(lái)流靜溫，Tj、rj分別為排氣溫度、排氣口半徑。I為排氣參數(shù)，形式為I=.其中分別為質(zhì)量流率、最大截面積、來(lái)流密度、來(lái)流速度。模擬的具體參數(shù)值在表1中給出。

圖1　底排模型示意圖Fig.1 Schematic diagram of base bleed model

表1　模擬參數(shù)Tab.1 Simulation parameters

網(wǎng)格采用弧長(zhǎng)法生成［17］。圖2為模型尾部區(qū)域網(wǎng)格圖，網(wǎng)格總數(shù)為29 000.遠(yuǎn)場(chǎng)采用無(wú)反射邊界條件，固壁采用無(wú)滑移邊界條件，中心軸線上采用對(duì)稱邊界條件，底排噴口的邊界條件直接給定。

根據(jù)Jachimowski［18］的H2和空氣的燃燒機(jī)理以及Gardiner的CO燃燒機(jī)理［19］，建立10組分25步反應(yīng)的H2-CO燃燒模型，作為尾部的二次燃燒模型。底排裝置排氣口組分參數(shù)參照文獻(xiàn)［8］。

圖2　底排模型尾部區(qū)域網(wǎng)格Fig.2 Grid of base region of base bleed model

2.2數(shù)值計(jì)算格式

采用有限體積法編程求解方程組。根據(jù)劉君等［20］的統(tǒng)一算法的思路，對(duì)方程組中的總能e進(jìn)行變換，然后將得到的新方程組分成3個(gè)部分分別進(jìn)行離散求解。其中，對(duì)流項(xiàng)通過(guò)改進(jìn)的AUSM+格式加入Van Leer限制器離散［21］，具有3階精度。擴(kuò)散項(xiàng)采用2階中心差分格式離散。時(shí)間項(xiàng)采用LUSGS隱式時(shí)間推進(jìn)方法［22］求解，具有1階精度。時(shí)間離散中粘性項(xiàng)采用近似隱式處理，湍動(dòng)能生成項(xiàng)顯示處理，組分方程中的擴(kuò)散項(xiàng)顯示處理?；瘜W(xué)反應(yīng)源項(xiàng)采用2階隱式梯形公式求解，以消除非平衡流動(dòng)的剛性問(wèn)題［20］。最后求出守恒變量后，通過(guò)牛頓迭代法求解溫度［23］，再通過(guò)分壓定理求壓強(qiáng)。和劉君等［20］的統(tǒng)一算法有些區(qū)別，這里的時(shí)間項(xiàng)通過(guò)LU-SGS隱式算法求解，并且湍流和Navier-Stokes方程之間采用全耦合方法求解。這使得計(jì)算的收斂速度加快，也更有利于程序由計(jì)算湍流改進(jìn)到計(jì)算湍流化學(xué)非平衡流的編制。同時(shí)仍不失劉君等［20］的統(tǒng)一算法在計(jì)算化學(xué)非平衡流時(shí)的優(yōu)點(diǎn)：對(duì)計(jì)算機(jī)內(nèi)存要求低，程序編制較簡(jiǎn)單。

3　數(shù)值模擬結(jié)果與分析

3.1數(shù)值模擬和實(shí)驗(yàn)結(jié)果的對(duì)比

選擇丁則勝等［5］的底部燃燒實(shí)驗(yàn)進(jìn)行對(duì)比。具體實(shí)驗(yàn)參數(shù)［4-5］如下：實(shí)驗(yàn)?zāi)Ｐ蜑閳A柱體，直徑為34.2 mm，排氣溫度為1 750 K，來(lái)流壓力為一個(gè)標(biāo)準(zhǔn)大氣壓，來(lái)流馬赫數(shù)為2.0，排氣口直徑為0.375倍圓柱體直徑，底排氣體為氫氣、氧化劑和稀釋劑的預(yù)混氣體。針對(duì)實(shí)驗(yàn)?zāi)Ｐ图皸l件進(jìn)行數(shù)值模擬，和實(shí)驗(yàn)進(jìn)行底壓的對(duì)比。

數(shù)值模擬迭代20 000步左右收斂，在普通PC機(jī)上計(jì)算10 h左右。圖3為底部平均壓強(qiáng)隨排氣參數(shù)變化的對(duì)比圖。從圖3可以看出，在底壓變化趨勢(shì)上模擬結(jié)果和實(shí)驗(yàn)結(jié)果較吻合。在小排氣參數(shù)時(shí)底壓速增，在I=0.009附近底壓達(dá)到峰值，然后隨排氣參數(shù)的增加底壓緩慢下降。在小排氣參數(shù)時(shí)，模擬結(jié)果與實(shí)驗(yàn)結(jié)果較吻合，隨排氣參數(shù)的增大，模擬結(jié)果略低于實(shí)驗(yàn)值，最大誤差在5%以內(nèi)。

圖3　底部平均壓強(qiáng)隨排氣參數(shù)變化曲線Fig.3 Area-averaged base pressure as function of exhaust parameter

3.2二次燃燒對(duì)底排加能的影響

對(duì)底部熱排氣和底部二次燃燒模型的尾部溫度場(chǎng)以及尾部單位體積內(nèi)能分布進(jìn)行比較，研究二次燃燒對(duì)底排加能的影響。

圖4為底部熱排氣和底部二次燃燒時(shí)模型尾部溫度分布圖。圖中可以看出：底部熱排氣時(shí)，排氣出噴口后溫度迅速下降，在下游x/R=1.0處溫度已經(jīng)降到1 000 K左右。底部二次燃燒時(shí)，整個(gè)尾部區(qū)域充滿高溫，在下游x/R=10.0處，中軸線上的溫度才開(kāi)始低于排氣溫度?？偟膩?lái)說(shuō)，二次燃燒使尾部的高溫區(qū)域大大增加。

圖4 模型尾部區(qū)域溫度分布（p∞=101 325 Pa，Tj=1 550 K，I=0.010 7）Fig.4 Temperature contour of base region（p∞=101 325 Pa，Tj=1 550 K，I=0.010 7）

圖5為模型底部冷排氣、熱排氣、二次燃燒時(shí)尾部單位體積內(nèi)能分布，其中內(nèi)能為無(wú)量綱值。無(wú)量綱的單位體積內(nèi)能的參考值為無(wú)窮遠(yuǎn)處的來(lái)流值，可表示為

式中：e∞、ρi∞、hi∞、p∞分別為無(wú)窮遠(yuǎn)處的單位體積內(nèi)能、組分i的密度、組分i的焓、壓力。從圖5可以看出：冷排氣時(shí)，底部有一個(gè)內(nèi)能較低的區(qū)域；排氣溫度升高后，尾部?jī)?nèi)能較低的區(qū)域有所減??；底部存在二次燃燒時(shí)尾部出現(xiàn)一個(gè)內(nèi)能較高的區(qū)域。對(duì)計(jì)算區(qū)域進(jìn)行積分得到單位體積內(nèi)能的平均值。底部冷排氣、熱排氣、二次燃燒時(shí)的內(nèi)能平均值分別為1.003 89、1.005 07、1.012 59.熱排氣使內(nèi)能平均值增加0.001 18，在此基礎(chǔ)上二次燃燒使內(nèi)能平均值又增加0.007 52，可見(jiàn)二次燃燒的加能效果是熱排氣的6.4倍。

圖5　尾部單位體積內(nèi)能分布（p∞=101 325 Pa，I=0.010 7，冷排氣Tj=293 K，熱排氣、二次燃燒Tj=1 550 K）Fig.5 Internal energy contours per unit volume of base region（p∞=101 325 Pa，I=0.010 7，Tj=293 K for cold base bleed，Tj=1 550 K for hot base bleed and secondary combustion）

3.3環(huán)境壓力降低對(duì)底排尾部流場(chǎng)的影響

圖6為不考慮二次燃燒的情況下，不同來(lái)流壓力時(shí)模型尾部區(qū)域無(wú)量綱密度場(chǎng)以及流線圖。圖6中密度的無(wú)量綱化參考值為各自的來(lái)流密度。由圖6可知，在相同的排氣參數(shù)下，亞大氣壓下的模型尾部無(wú)量綱化的密度場(chǎng)分布以及流線圖和常壓下基本相同。因?yàn)橐环矫?，在?lái)流溫度T∞以及來(lái)流組分不變的情況下，來(lái)流密度ρ∞隨來(lái)流壓力p∞的減小呈正比減小。根據(jù)公式，當(dāng)來(lái)流速度v∞和排氣參數(shù)I不變時(shí)，來(lái)流密度ρ∞和排氣的質(zhì)量流率呈正比。所以排氣的質(zhì)量流率應(yīng)隨來(lái)流壓力p∞的減小呈正比減小。另一方面，由于底排為亞聲速排氣，來(lái)流壓力的減小會(huì)影響到底排裝置內(nèi)部，使得排氣的壓力相應(yīng)減小，因?yàn)榕艢鉁囟炔蛔?，所以排氣的密度也相?yīng)減小。在圖6中，不同來(lái)流壓力下，底排模型尾部的ρ/ρ∞場(chǎng)基本相同，可見(jiàn)底排模型尾部區(qū)域的密度基本隨來(lái)流密度的減小呈正比減小，即排氣的密度基本隨來(lái)流壓力的減小呈正比減小。由于底排的質(zhì)量流率和排氣密度基本都隨來(lái)流壓力的減小呈正比減小，使得排氣的初速度隨來(lái)流壓力的減小基本不變。所以不同來(lái)流壓力下，工作在相同排氣參數(shù)下的模型尾部流線圖也基本一樣。

圖6　底排模型尾部區(qū)域無(wú)量綱密度場(chǎng)以及流線圖（Tj=1 550 K，I=0.010 7）Fig.6 Density contour and streamline chart of base region（Tj=1 550 K，I=0.010 7）

表2為不考慮二次燃燒的情況下，不同來(lái)流壓力p∞下的底排減阻率RCDB.由表2可見(jiàn)，在相同的排氣參數(shù)下，亞大氣壓下的底排減阻率和常壓下相差很小。來(lái)流壓力由常壓降到31 500 Pa，底排減阻率只下降了2.8%。說(shuō)明相同排氣參數(shù)下，不考慮二次燃燒時(shí)，來(lái)流壓力的減小對(duì)底排增壓減阻的影響很小，即只要保持排氣參數(shù)相同，排氣的質(zhì)量流率不同所帶來(lái)的影響很小。

表2　不考慮二次燃燒時(shí)的底排減阻率（Tj=1 550 K,I=0.010 7）Tab.2 Drag reduction rate of base bleed without secondary combustion（Tj=1 550 K,I=0.010 7）

圖7 模型尾部區(qū)域溫度場(chǎng)以及流線圖（Tj=1 550 K，I=0.010 7）Fig.7 Temperature contours and streamline charts of the base region（Tj=1 550 K，I=0.010 7）

圖7為不同來(lái)流壓力p∞時(shí)，模型尾部區(qū)域溫度場(chǎng)以及流線圖。圖7中上半部分為考慮二次燃燒的情況，下半部分為不考慮二次燃燒的情況。由圖7可見(jiàn)：不考慮二次燃燒時(shí)，來(lái)流壓力的變化對(duì)模型尾部的溫度場(chǎng)基本沒(méi)有影響。模型尾部的溫度分布規(guī)律為：排氣口附近溫度最高，下游溫度逐漸降低，在環(huán)狀回流區(qū)內(nèi)溫度降到1 100 K左右，在主回流區(qū)內(nèi)溫度降到750 K左右。考慮二次燃燒后，來(lái)流壓力的變化對(duì)模型尾部的溫度場(chǎng)以及回流區(qū)大小和位置影響很大。來(lái)流壓力為101 325 Pa時(shí)，兩個(gè)回流區(qū)內(nèi)以及下游部分區(qū)域的溫度都升高到2 000 K以上，主回流區(qū)變得很小。當(dāng)來(lái)流壓力降到68 700 Pa時(shí)，尾部區(qū)域的溫度明顯降低，2 000 K左右的高溫區(qū)域只集中在兩個(gè)回流區(qū)內(nèi)，且主回流區(qū)逐漸變大。當(dāng)來(lái)流壓力降到31 500 Pa時(shí)，只在底部固壁附近溫度略微升高，下游溫度緩慢降低，主回流區(qū)內(nèi)溫度只有1 350 K左右，且主回流區(qū)和無(wú)二次燃燒時(shí)相比基本一樣大，只是略微前移。總的來(lái)說(shuō)，二次燃燒釋放的熱量使尾部溫度升高。但是隨著來(lái)流壓力的減小，尾部區(qū)域溫度的升高程度逐漸減小。

表3為考慮二次燃燒的情況下，不同來(lái)流壓力p∞下的底排減阻率RCDB.對(duì)比表3和表2可看出：來(lái)流壓力分別為101325 Pa、68700 Pa、31500 Pa時(shí)，二次燃燒盡管使底排減阻率分別升高了54.2%、46.1%、30.4%，但是隨著來(lái)流壓力的降低，二次燃燒對(duì)底排增壓減阻的貢獻(xiàn)卻逐漸變小。

表3　考慮二次燃燒時(shí)的底排減阻率（Tj=1 550 K, I=0.010 7）Tab.3 Drag reduction rate of base bleed with secondary combustion（Tj=1 550 K,I=0.010 7）

圖8為不同來(lái)流壓力p∞時(shí)，模型尾部區(qū)域H2質(zhì)量分?jǐn)?shù)分布圖。圖8中上半部分為考慮二次燃燒的情況，下半部分為不考慮二次燃燒的情況。由圖8可見(jiàn)：不考慮二次燃燒時(shí)，隨著來(lái)流壓力的減小，模型尾部區(qū)域的H2質(zhì)量分?jǐn)?shù)分布變化很小?？紤]二次燃燒后，在底部固壁附近以及下游主回流區(qū)內(nèi)，H2質(zhì)量分?jǐn)?shù)明顯減小，且隨著來(lái)流壓力的減小，底部固壁附近的H2質(zhì)量分?jǐn)?shù)的減小程度逐漸被削弱。說(shuō)明在主回流區(qū)和環(huán)狀回流區(qū)內(nèi)均存在H2的燃燒，但是隨著來(lái)流壓力的減小，環(huán)狀回流區(qū)內(nèi)H2質(zhì)量消耗率降低，其燃燒效率逐漸降低。

圖9為不同來(lái)流壓力p∞時(shí)，模型尾部區(qū)域H質(zhì)量分?jǐn)?shù)分布圖。由圖9可見(jiàn)：隨著來(lái)流壓力的減小，模型尾部的H質(zhì)量分?jǐn)?shù)峰值逐漸增大，且逐漸向環(huán)狀回流區(qū)靠近。說(shuō)明隨著來(lái)流壓力的減小，混合氣體在環(huán)狀回流區(qū)附近燃燒不充分，導(dǎo)致中間產(chǎn)物增多，這與圖8結(jié)果一致。并且從圖5還可看出，常壓下，底排氣體排出后先沿著主回流區(qū)和環(huán)狀回流區(qū)之間的夾縫流入剪切層，H2和O2在剪切層中開(kāi)始反應(yīng)，反應(yīng)中的混合氣體在流入兩個(gè)回流區(qū)以及下游的過(guò)程中繼續(xù)燃燒直至燃完。所以在剪切層中中間產(chǎn)物最多，在兩個(gè)回流區(qū)以及下游區(qū)域逐漸變小。

圖8　模型尾部區(qū)域H2質(zhì)量分?jǐn)?shù)分布圖（Tj=1 550 K，I=0.010 7）Fig.8 Mass fraction of hydrogen in base region（Tj=1 550 K，I=0.010 7）

綜上所述，底排燃?xì)馀懦龊笱刂鴥蓚€(gè)回流區(qū)之間的夾縫流入剪切層，在剪切層中遇到O2開(kāi)始反應(yīng)，反應(yīng)中的混合氣體在流入兩個(gè)回流區(qū)以及下游的過(guò)程中繼續(xù)燃燒。隨著來(lái)流壓力的減小，環(huán)狀回流區(qū)內(nèi)H2的燃燒效率逐漸降低，中間產(chǎn)物H逐漸增多，燃燒逐漸變得不充分，導(dǎo)致底排減阻率明顯降低。

圖9　模型尾部區(qū)域H質(zhì)量分?jǐn)?shù)分布圖（Tj=1 550 K，I=0.010 7）Fig.9 Mass fraction of hydrogen atom in base region（Tj=1 550 K，I=0.010 7）

4　結(jié)論

1）二次燃燒使模型尾部的高溫區(qū)域劇增，對(duì)尾部區(qū)域加能的貢獻(xiàn)是熱排氣的6.4倍，是底排加能減阻的關(guān)鍵。

2）在相同排氣參數(shù)下，不考慮二次燃燒時(shí)，來(lái)流壓力的減小對(duì)模型尾部流場(chǎng)以及底排減阻率的影響很小。考慮二次燃燒后，模型尾部的溫度升高，但是隨著來(lái)流壓力的減小，尾部區(qū)域溫度的升高程度逐漸減小，且隨著來(lái)流壓力的減小，二次燃燒對(duì)底排增壓減阻的貢獻(xiàn)也逐漸變小。

3）排氣排出后沿著兩個(gè)回流區(qū)之間的夾縫流入剪切層，在剪切層中遇到O2開(kāi)始反應(yīng)，反應(yīng)中的混合氣體在流入兩個(gè)回流區(qū)以及下游的過(guò)程中繼續(xù)燃燒。隨著來(lái)流壓力的減小，環(huán)狀回流區(qū)內(nèi)H2的燃燒逐漸變得不充分，導(dǎo)致底排減阻率明顯降低。

（

）

［1］丁則勝，邱光純，劉亞飛，等.固體燃料底部排氣空氣動(dòng)力研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，1991，9（3）：300-307. DING Ze-sheng，QIU Guang-chun，LIU Ya-fei，et al.An aerodynamic investigation of base bleed by solid fuel［J］.Acta Aerodynamica Sinica，1991，9（3）：300-307.（in Chinese）

［2］Mathur T，Dutton J C.Velocity and turbulence measurements in a supersonic base flow with mass bleed［J］.AIAA Journal，1996，34（6）：1153-1159.

［3］Bourdon C J，Dutton J C.Visualization of a central bleed jet in an axisymmetric compressible base flow［J］.Physics of Fluids，2003，15（2）：499-510.

［4］Bowman J E，Clayden W A.Cylindrical afterbodies at M=2 with hot gas ejection［J］.AIAA Journal，1968，6（12）：2429-2431.

［5］丁則勝，羅榮，陳少松，等.底部燃燒減阻性能的若干參數(shù)影響研究［J］.彈道學(xué)報(bào)，1996，8（4）：79-83. DING Ze-sheng，LUO Rong，CHEN Shao-song，et al.A study of some parameters influence on performance of drag reduction by base burning［J］.Journal of Ballistics，1996，8（4）：79-83.（in Chinese）

［6］丁則勝，陳少松，劉亞飛，等.底排性能的環(huán)境壓力效應(yīng)［J］.彈道學(xué)報(bào)，2002，14（1）：88-92. DING Ze-sheng，CHEN Shao-song，LIU Ya-fei，et al.Influence of ambient pressure on base bleed［J］.Journal of Ballistics，2002，14（1）：88-92.（in Chinese）

［7］Sahu J，Nietubicz C J，Steger J L.Navier-Stokes computations of projectile base flow with and without base injection［J］.AIAA Journal，1985，23（9）：1348-1355.

［8］Gibeling H J，Buggeln R C.Projectile base bleed technology part 1：analysis and results，AD-A258459［R］.Glastonbury，CT：Scientific Research Associates，1992.

［9］Jachimowski C J.An analytical study of the hydrogen-air reaction mechanism with application to scramjet combustion，NASA-TP-2791［R］.Hampton，VA：Langley Research Center，1988.

［10］陸中兵，丁玨，周彥煌，等.超聲速飛行底部排氣彈三維湍流流場(chǎng)數(shù)值模擬［J］.南京理工大學(xué)學(xué)報(bào)，2007，31（1）：27-30. LU Zhong-bing，DING Jue，ZHOU Yan-huang，et al.Numerical investigation of 3-D turbulent flow fields for supersonic flying base bleed projectile［J］.Journal of Nanjing University of Science and Technology，2007，31（1）：27-30.（in Chinese）

［11］陳新虹，黃華，周志超，等.排氣能量對(duì)底部排氣彈氣動(dòng)特性影響的數(shù)值模擬［J］.兵工學(xué)報(bào)，2010，31（4）：447-452. CHEN Xin-hong，HUANG Hua，ZHOU Zhi-chao，et al.Numerical simulation of base bleed energy affecting aerodynamic performance of base bleed projectiles［J］.Acta Armamentarii，2010，31（4）：447-452.（in Chinese）

［12］Shin J R，Cho D R，Won S H，et al.Hybrid RANS/LES study of base-bleed flows in supersonic mainstream［C］∥15th AIAAInternational Space Planes and Hypersonic Systems and Technologies Conference.Dayton，Ohio：AIAA，2008.

［13］Shin J R，Choi J Y.DES study of base and base-bleed flows with dynamic formulation of DES constant［C］∥49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition.Orlando，F(xiàn)lorida：AIAA，2011.

［14］Menter F R.Two-equation eddy-viscosity turbulence models for engineering application［J］.AIAA Journal，1994，32（8）：1598-1605.

［15］余文杰，余永剛，倪彬.底部排氣圓柱體模型尾部流場(chǎng)的數(shù)值模擬［J］.彈道學(xué)報(bào)，2014，26（1）：7-12. YU Wen-jie，YU Yong-gang，NI Bin.Numerical simulation of base flow field over a cylindrical model with base bleed［J］.Journal of Ballistics，2014，26（1）：7-12.（in Chinese）

［16］代淑蘭，許厚謙，王兵.含高速運(yùn)動(dòng)彈丸的膛口二次燃燒并行數(shù)值模擬［J］.彈道學(xué)報(bào)，2009，21（1）：83-86. DAI Shu-lan，XU Hou-qian，WANG Bing.Numerical simulation of secondary muzzle flash including high-speed projectile using parallel computation method［J］.Journal of Ballistics，2009，21（1）：83-86.（in Chinese）

［17］武頻，趙潤(rùn)祥，郭錫福.弧長(zhǎng)網(wǎng)格生成法及其應(yīng)用［J］.南京理工大學(xué)學(xué)報(bào)，2002，26（5）：482-485. WU Pin，ZHAO Run-xiang，GUO Xi-fu.Arc length method of grid generation and its application［J］.Journal of Nanjing University of Science and Technology，2002，26（5）：482-485.（in Chinese）

［18］Jachimowski C J.An analytical study of the hydrogen-air reaction mechanism with application to scramjet combustion，NASA-TP-2791［R］.Washington，DC：NASA，1988.

［19］Gardiner W C.Combustion chemistry［M］.New York：Springer-Verlag，1984.

［20］劉君，張涵信，高樹(shù)椿.一種新型的計(jì)算化學(xué)非平衡流動(dòng)的解耦方法［J］.國(guó)防科技大學(xué)學(xué)報(bào)，2000，22（5）：19-22. LIU Jun，ZHANG Han-xin，GAO Shu-chun.A new uncoupled method for numerical simulation of nonequilibrium flow［J］.Journal of National University of Defense Technology，2000，22（5）：19-22.（in Chinese）

［21］梁德旺，王可.AUSM+格式的改進(jìn)［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2004，22（4）：404-409. LIANG De-wang，WANG Ke.Improvement of AUSM+scheme［J］.Acta Aerodynamica Sinica，2004，22（4）：404-409.（in Chinese）

［22］Yoon S，Jameson A.Lower-upper symmetric Gauss-Seidel method for the Euler and Navier-Stokes equations［J］.AIAA Journal，1988，26（9）：1025-1026.

［23］劉晨.復(fù)雜燃燒流場(chǎng)數(shù)值模擬方法研究［D］.南京：南京航空航天大學(xué)，2009. LIU Chen.Numericalmethods for complex combustion flow fields［D］.Nanjing：Nanjing University of Aeronautics and Astronautics，2009.（in Chinese）

Numerical Simulation on the Effect of Reduction in Ambient Pressure on the Secondary Combustion of Base Bleed

YU Wen-jie，YU Yong-gang
（School of Energy and Power Engineering，Nanjing University of Science and Technology，Nanjing 210094，Jiangsu，China）

In order to investigate the mechanism of the decrease in drag reduction rate of base bleed in the subatmospheric pressure，a mathematical and physical model about the base flow with chemical nonequilibrium of base bleed is established.H2-CO combustion model which consists of 10 components and 25 reactions is used for secondary combustion.Two-dimensional axisymmetric equations are programmatically computed using a set of uniform numerical process methods.The base flow field of base bleed is simulated.Simulation results are validated with experiment data.On this basis，the base flow field and combustion characteristics are numerically predicted.The results show that the added energy released from secondary combustion is 6.4 times of the added energy of hot base bleed.Secondary combustion is a key for energy increase and drag reduction.With the decrease of the ambient pressure，the combustion efficiency of H2decreases gradually in the annular recirculation zone of the tail of the model，meanwhile the intermediate product of H increases gradually.These make the secondary combustion become more and more insufficiency，resulting in decreasing significantly the drag reduction rate of base bleed.

ordnance science and technology；base bleed；secondary combustion；numerical simula-tion；base flow field；chemical non-equilibrium flow

V211.3

A

1000-1093（2015）03-0443-08

10.3969/j.issn.1000-1093.2015.03.010

2014-06-04

國(guó)家自然科學(xué)基金項(xiàng)目（51176076）

余文杰（1986—），男，博士研究生。E-mail：spacecow@sina.com；余永剛（1963—），男，教授，博士生導(dǎo)師。E-mail：yyg801@njust.edu.cn