宋天，單智超，丁林

（1.重慶大學(xué)低品位能源利用技術(shù)及系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400044；2.中國(guó)航發(fā)四川燃?xì)鉁u輪研究院，四川 綿陽 621000)

燃?xì)馊臃治黾夹g(shù)是航空發(fā)動(dòng)機(jī)燃燒室效率測(cè)量、溫度測(cè)量及排放指標(biāo)評(píng)估最重要的技術(shù)手段。然而，在當(dāng)前的取樣技術(shù)中，對(duì)于燃燒室全包線范圍的出口氣流狀態(tài)，取樣探針始終無法將樣氣保持在合適的取樣溫度范圍，這給燃燒效率測(cè)量和燃燒室出口溫度計(jì)算帶來了較大的系統(tǒng)誤差。為了解決該問題，可通過換熱數(shù)值分析設(shè)計(jì)合理的取樣換熱結(jié)構(gòu)，進(jìn)而采用單獨(dú)設(shè)置取樣循環(huán)保溫系統(tǒng)對(duì)樣氣進(jìn)行控溫。其中涉及的關(guān)鍵問題之一是高溫高壓氣體管內(nèi)流動(dòng)熱流固耦合及散熱控制。國(guó)外在高超聲速流場(chǎng)、熱和結(jié)構(gòu)之間多場(chǎng)耦合問題方面進(jìn)行了長(zhǎng)期系統(tǒng)的研究，無論是基礎(chǔ)研究還是實(shí)際工程應(yīng)用研究都取得了許多有益的進(jìn)展，多場(chǎng)耦合分析的思想和應(yīng)用已經(jīng)逐漸深入工程實(shí)際。國(guó)內(nèi)研究主要集中于流場(chǎng)的氣動(dòng)加熱與固體結(jié)構(gòu)傳熱之間的耦合問題，對(duì)于高超聲速流場(chǎng)、熱和結(jié)構(gòu)之間多物理場(chǎng)耦合特征和規(guī)律的認(rèn)識(shí)還不夠深入，缺乏長(zhǎng)期的、系統(tǒng)的基礎(chǔ)研究，同時(shí)也缺乏實(shí)際的工程應(yīng)用研究，其相關(guān)多場(chǎng)耦合建模與分析方法和研究手段還有待進(jìn)一步發(fā)展和完善。

本文針對(duì)高溫燃?xì)馊舆^程中的高溫高壓氣體在取樣管內(nèi)的流動(dòng)換熱問題展開研究，詳細(xì)闡述了取樣管的物理模型、計(jì)算區(qū)域及網(wǎng)格劃分等，并對(duì)燃?xì)馊庸軆?nèi)的流體流動(dòng)Ma 數(shù)、壓力、換熱和出口溫度控制等進(jìn)行分析。研究結(jié)果可為燃?xì)馊臃治黾夹g(shù)的發(fā)展提供理論支撐。

1 物理模型

圖1 所示為高溫高壓燃?xì)馊庸芪锢砟Ｐ?。取樣管材料?04 不銹鋼。燃?xì)馔ㄟ^樣氣入口進(jìn)入取樣管，與管外冷卻介質(zhì)完成熱量交換后由樣氣出口流出。樣氣入口為變截面圓管，后續(xù)管道為等截面圓管。為了降低外部邊界條件對(duì)計(jì)算結(jié)果的影響，考慮取樣點(diǎn)外部燃?xì)猸h(huán)境。樣氣入口與外部流場(chǎng)形成耦合流動(dòng)邊界。沿取樣管軸向，樣氣與取樣管固體域界面考慮為熱-流-固耦合壁面邊界；取樣管外壁與冷卻水形成耦合壁面。表1 所示為不同取樣管入口條件及其詳細(xì)參數(shù)。

表1 不同入口邊界條件及詳細(xì)參數(shù)

圖1 物理模型

2 數(shù)值計(jì)算方法

2.1 控制方程

高溫高壓燃?xì)庠谌庸軆?nèi)流動(dòng)及換熱的三維效應(yīng)不可忽略。取樣過程是三維瞬態(tài)熱-流-固多場(chǎng)耦合過程。流動(dòng)與換熱過程中的守恒關(guān)系包括質(zhì)量守恒、動(dòng)量守恒和能量守恒，控制方程包括：

本文數(shù)值計(jì)算基于有限體積法，對(duì)流項(xiàng)采用二階迎風(fēng)格式，梯度項(xiàng)采用least squares cell-based 方法，界面上壓力求解采用二階精度插值，瞬態(tài)項(xiàng)采用隱式格式，采用COUPLED 算法實(shí)現(xiàn)離散方程的求解。由于樣氣流體處于高溫、高速、壓力梯度較大的流動(dòng)狀態(tài)，且存在近壁面黏性流動(dòng)，因此本文計(jì)算采用k-ω SST 湍流模型，適用于取樣過程的數(shù)值求解。

2.2 網(wǎng)格劃分

計(jì)算區(qū)域網(wǎng)格采用結(jié)構(gòu)化和非結(jié)構(gòu)化混合的網(wǎng)格策略，對(duì)每個(gè)分區(qū)進(jìn)行網(wǎng)格尺寸控制然后劃分網(wǎng)格。取樣管流體及固體域采用結(jié)構(gòu)化網(wǎng)格生成方案，如圖2 所示。彎管處采用多個(gè)拉伸塊保證貼體性，截面變化處切分塊保證網(wǎng)格映射。經(jīng)過網(wǎng)格無關(guān)性驗(yàn)證后，計(jì)算采用的取樣管網(wǎng)格數(shù)量為3181824。

圖2 網(wǎng)格劃分

3 結(jié)果與討論

圖3a 為燃?xì)馊庸軆?nèi)Ma 數(shù)隨取樣管長(zhǎng)度的變化曲線。以無量綱數(shù)L*表示取樣管橫截面距離入口的相對(duì)位置，L*=0 代表取樣管入口位置，L*=1 代表取樣管出口位置。不同入口條件下取樣管內(nèi)Ma 數(shù)整體變化規(guī)律是相似的，由于取樣管入口段截面幾何特點(diǎn)為先變小后增大的縮放結(jié)構(gòu)，樣氣進(jìn)入取樣管后Ma 數(shù)首先急劇上升，然后大幅下降。條件1 和2 樣氣入口狀態(tài)相同，隨著出口壓力降低，管內(nèi)樣氣最大Ma 數(shù)增加，而且在L*>0.1后流動(dòng)也存在更大的波動(dòng)。由此可見，取樣管出口壓力的變化對(duì)管內(nèi)樣氣Ma 數(shù)存在明顯影響，直接影響實(shí)際工程取樣中樣氣流速和流量，將進(jìn)一步影響樣氣與冷卻液之間的熱量交換。為了進(jìn)一步分析燃?xì)馊舆^程中取樣管內(nèi)樣氣特性變化，圖3b 給出了2 個(gè)入口條件取樣管內(nèi)樣氣靜壓隨管長(zhǎng)的變化曲線。從圖中可以看出，當(dāng)L*>0.1 后，樣氣壓力隨取樣管長(zhǎng)度增加緩慢降低，當(dāng)取樣管進(jìn)出口壓差較大時(shí)，壓力變化存在較大的波動(dòng)。

圖3 不同出口壓力下(a)Ma 數(shù)和(b)靜壓沿L＊變化

圖4 為2 個(gè)燃?xì)馊庸r所對(duì)應(yīng)的管內(nèi)樣氣溫度隨管長(zhǎng)的變化曲線。為了將樣氣在取樣管出口處的溫度控制在取樣溫度要求（433K），通過液態(tài)水實(shí)現(xiàn)樣氣冷卻。條件1 與條件3 出入口邊界條件均取自實(shí)驗(yàn)實(shí)際工作條件，取樣管出口樣氣溫度與目標(biāo)值具有一定差距，通過改變出口壓力，例如條件2，將出口壓力由101.33kPa降至51.38kPa，可有效控制出口樣氣溫度在目標(biāo)溫度433K。

圖4 不同出口壓力下溫度沿L＊變化

圖5 給出了入口條件2 取樣管入口中心剖面的壓力、Ma 數(shù)和溫度分布云圖。取樣管入口為縮放幾何結(jié)構(gòu)，高溫高壓燃?xì)庠诹鹘?jīng)入口段時(shí)具有較大的變化梯度，結(jié)合圖3 ～5 的定量分析結(jié)果可以看出，樣氣Ma 數(shù)、壓力和溫度在入口段都存在劇烈的波動(dòng)。

圖5 條件2 取樣管入口壓力、Ma、溫度云圖

保持其他參數(shù)不變，通過改變?nèi)犹结槼隹诒硥?，調(diào)節(jié)樣氣在取樣管內(nèi)的流速和流量，進(jìn)而實(shí)現(xiàn)樣氣與冷卻水之間的換熱量調(diào)控，從而可實(shí)現(xiàn)取樣管出口樣氣溫度控制。為了對(duì)比分析不同入口條件下改變?nèi)庸艹隹诒硥簳r(shí)所得到的樣氣出口狀態(tài)，圖6 給出了取樣管出口樣氣平均溫度隨壓比Pout/Pamb的變化曲線。對(duì)于入口條件1 和入口條件2，在所測(cè)試的壓比變化范圍（0.30 ≤Pout/Pamb≤0.90），出口樣氣平均溫度可以實(shí)現(xiàn)390 ～430K 范圍內(nèi)的調(diào)節(jié)，出口溫度達(dá)到預(yù)期出口樣氣調(diào)節(jié)目標(biāo)。

圖6 出口溫度隨壓比變化及擬合曲線

4 結(jié)語

（1）由于取樣管入口段截面為先變小后增大的縮放結(jié)構(gòu)，樣氣壓力在取樣管入口處先驟降后上升。相應(yīng)地，樣氣Ma 數(shù)在入口段首先急劇上升至Ma>1，然后大幅下降。在L*>0.1 后，樣氣壓力和Ma 數(shù)隨取樣管長(zhǎng)度增加緩慢降低。（2）樣氣在入口段出現(xiàn)溫度驟降。在L*>0.1 后，樣氣溫度隨取樣管長(zhǎng)度從1800K 左右逐漸降低至433K 左右。（3）通過調(diào)節(jié)出口壓比可調(diào)節(jié)出口樣氣溫度，出口樣氣溫度隨出口壓比的增大而減小。在所測(cè)試的壓比變化范圍內(nèi)（0.30 ≤Pout/Pamb≤0.90），出口樣氣平均溫度可在390 ～430K 的范圍內(nèi)進(jìn)行調(diào)節(jié)。