先進旋渦燃燒室鈍體結構參數選擇的數值分析*

孫海俊，曾卓雄，徐義華

(南昌航空大學飛行器工程學院，南昌 330063)

0 引言

穩定的燃燒火焰和較高的燃燒效率是燃燒室結構性能優化的主要目標，也是沖壓發動機研究所需要解決的關鍵問題。美國的Ramgen公司在駐渦燃燒室[1](TVC，trapped vortex combustor)的研究基礎上，提出了先進旋渦燃燒室[2](AVC，advanced vortex combustor)的概念。AVC是在燃燒室內布置前后兩個獨立的鈍體，當來流經過前鈍體時，在前后鈍體之間的凹腔處形成穩定的旋渦結構，為燃燒室充當穩定的點火源;另一方面，穩定旋渦的存在不僅可以有效地降低燃燒室的總壓損失和污染物的排放，還可以極大的提高燃燒效率[3]。然而旋渦的穩定和燃燒流動特性的優劣主要取決于前后鈍體結構參數的匹配關系。文獻[4－5]等依據 Mair[6]得出的 TVC 穩定駐渦條件的結構關系式對AVC內鈍體結構參數和布置進行了研究。雖然AVC和TVC在結構上有相似之處，但仍需要區別對待。

因此，對不同鈍體結構參數影響下的燃燒室流場進行數值模擬研究，以確定前后鈍體各參數間的最佳匹配關系，為進一步深入研究提供理論依據。

1 物理模型與結構參數

AVC物理模型如圖1所示，鈍體布置方式不同于文獻[7－8]之處在于將鈍體前壁面直接與通道齊平，將原來的大通道進氣口變成上下兩個小通道進氣口。燃燒室通道尺寸為360mm×100mm×100mm，前后鈍體寬度為100mm，鈍體的高和長以及兩鈍體間距S之間的匹配關系是文中研究的重點。

數值模擬采用不可壓N-S方程，湍流模型為realizable k-ε模型，近壁面采用壁面函數法。計算時入口區域延長，燃燒室進出口邊界條件分別為速度入口與壓力出口。

圖1 AVC物理模型

研究參數:

1)進口速度V=100m/s。

2)前鈍體高 H1=50mm，60mm，70mm。

3)后鈍體高為 H2，H2/H1=0.5，0.6，0.7。

4)凹腔長為 S，S/H1=0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9，1.0，1.1，1.2。

5)前鈍體長L1=40mm，50mm，60mm;后鈍體長L2，L2/L1=0.3，0.4，0.5，0.6，0.7，0.8，0.9，1.0。

研究H1、H2與S之間的關系時，L1=40mm，L2=20mm。

2 研究結果分析

2.1 兩種鈍體布置方式的比較分析

總壓損失系數定義為:

模型Ⅰ(見圖1)為所要研究的鈍體布置方式，模型Ⅱ為文獻[7－8]中的布置方式。在相同工況下對兩種模型進行數值模擬，其靜壓分布如圖2～圖3所示。

圖2 模型Ⅰ靜壓分布圖

總壓損失計算結果表明，模型Ⅰ的總壓損失系數為1.62%，模型Ⅱ總壓損失系數為24.3%。可見，模型Ⅰ的總壓損失較模型Ⅱ小很多。模型Ⅱ的總壓損失主要來源于兩方面，一是凹腔處形成的低壓槽;二是當主氣流流經前鈍體時，由于前鈍體阻擋，在前鈍體上下面上發生了擾流，從而形成局部低壓槽(見圖3)，局部低壓槽的存在不可避免增大了總壓損失。模型Ⅰ相對模型Ⅱ具有總壓損失小的明顯優勢，并且在燃燒室總長一定時，模型Ⅰ有較長的流動發展區，有利于尾跡的充分發展。

圖3 模型Ⅱ靜壓分布圖

下文對鈍體結構參數與流動性能關系的研究均基于模型Ⅰ的鈍體布置方式。

2.2 鈍體高度、間距對總壓損失的影響

前鈍體高度H1為50mm、60mm、70mm時，H2/H1、S/H1與總壓損失系數之間的關系如圖4～圖6所示。

隨H2/H1增大，圖中曲線呈現出兩條規律:1)曲線先降后升，即每一個H2/H1都存在一個S/H1使得總壓損失最小;2)曲線最小值點左移，即最小總壓損失對應的S/H1隨之減小。

圖4 H1=50mm時，S/H1與δ*間的關系曲線

圖5 H1=60mm時，S/H1與δ*間的關系曲線

圖6 H1=70mm時，S/H1與δ*間的關系曲線

每條曲線的最小值點所對應的結構參數當量比如表1所示。分析表1數據可知，隨前鈍體高度增加，相同的H2/H1、S/H1無量綱條件下，曲線最小值也隨之增大(總壓損失系數增大)。當H1=60mm和70mm時，最小值點所對應的結構參數關系均為H2/H1=0.7、S/H1=0.6，H1=50mm 時對應的結構參數為H2/H1=0.6、S/H1=0.6，并不滿足TVC 穩定駐渦形成結構關系式[6]，故AVC結構參數的選擇不能完全照搬TVC模式。

為探究結構參數對流場的影響機理，以表1中第二組結構參數為例，對凹腔內旋渦流場進行分析，其速度矢量分布如圖7～圖10所示。由圖可見，凹腔內均能形成正向分布的對稱旋渦。圖7因凹腔長度與后鈍體高度之比較大，致使旋渦被拉長，渦心偏離凹腔并超出后鈍體邊界，有效旋渦區變窄;另在后鈍體前壁處形成了一對小尺度反向旋渦，與正向渦在邊界層有剪切相互作用，從而削弱了大尺度渦的穩定發展勢力，并形成低壓槽，增大總壓損失。隨后鈍體高度的增加和凹腔長度的縮短，反向渦逐漸消失，大尺度渦中心隨之向凹腔內遷移，旋渦趨于穩定(見圖9);而當凹腔長度與后鈍體高度之比較小時(見圖10)，旋渦被完全擠壓在凹腔內，得不到充分的發展，無益于穩定旋渦和降低總壓損失。

表1 最小總壓損失系數對應的結構參數

圖7 H2/H1=0.5(L/H1=0.8)

圖8 H2/H1=0.6(S/H1=0.7)

圖9 H2/H1=0.7(S/H1=0.6)

圖10 H2/H1=0.8(S/H1=0.4)

后鈍體高度以及凹腔長度對旋渦穩定和總壓損失的影響相對較大，由上述分析可知，H1=60mm時H2/H1=0.7、S/H1=0.6 的結構關系對穩定旋渦和減小總壓損失具有較好的效果。

H1=50mm時，該組曲線的最小值點所對應的結構關系式為 H2/H1=0.6、S/H1=0.5，且總壓損失為三組數據中最小。三組數據中各最小總壓損失所對應的結構參數下的流場靜壓分布如圖11～圖13所示。由圖可見，凹腔處均形成了對稱分布的低壓槽，且發展較為充分，有助于旋渦的穩定。相比之下，結構參數為H1=50mm，H2/H1=0.6，S/H1=0.5 時(見圖11)，低壓槽中心脫離后鈍體阻擋區域，使旋渦容易受主氣流影響，不能穩定在凹腔范圍內，而其余兩組結構參數下的低壓槽中心均在鈍體阻擋范圍內，具有很好的穩定旋渦作用。

圖11 H1=50mm靜壓圖

圖12 H1=60mm靜壓圖

圖13 H1=70mm靜壓圖

由此可見，雖然第一組結構參數(H1=50mm時)下的總壓損失最小，但由于其低壓槽中心向邊偏移，在穩定旋渦方面并不十分理想;二、三組在穩定旋渦方面并無太大差異，但第二組參數下的總壓損失較第三組小，即第二組結構下的燃燒室綜合了旋渦穩定和總壓損失較小等優點，其對應的結構關系式H2/H1=0.7、S/H1=0.6應為理想的參數匹配關系式。

2.3 鈍體長度對壓力損失的影響

圖14 L2/L1與δ*間的關系曲線

圖14 所示。從圖可知，隨前后鈍體長度的增加，總壓損失系數均隨之增大，但幅度變化不大。這在于鈍體長度的增加，縮短了鈍體后的流動發展區，使湍流流動發展不充分，從而增大總壓損失;保證湍流流動的充分發展，有助于降低總壓損失。而鈍體長度的選擇并非越小越好，由于駐渦燃燒室內燃料注入通常采用壁面射流方式，所以還應考慮噴孔布置等因素。

綜合來看，相比鈍體高度和凹腔長度，鈍體長度對總壓損失的影響要小，在確定鈍體長度時，應結合總壓損失、湍流流動發展以及燃料射流孔布置等因素共同考慮。

基于相對最優結構參數關系(H1=60mm、H2/H1=0.7、S/H1=0.6)，取前鈍體長度分別為 30mm、40mm、50mm時，總壓損失系數與鈍體長度的關系如

3 結論

對鈍體布置方式和鈍體結構參數影響下的先進旋渦燃燒室流場的數值模擬結果表明:

1)燃燒室模型Ⅰ中鈍體布置方式具有較好的流動特性，且前后鈍體結構參數當量為H1/H2=0.7、S/H1=0.6時，凹腔內旋渦最為穩定對稱，且燃燒室總壓損失也較小;

2)AVC結構參數的選擇不能完全照搬TVC穩定駐渦形成結構關系式，需要區別對待;

3)增大鈍體長度會相應增大燃燒室總壓損失，確定鈍體長度時，應綜合考慮總壓損失、湍流流動發展以及燃料噴孔布置等因素。

[1]Hsu K Y，Goss L P，Trump D D，et al.Performance of a trapped-vortex combustor[C]//AIAA 33rd Aerospace Sciences Meeting and Exhibit，1995.

[2]KENDRICK D W，CHENEVERT B C，TRUEBLOOD B，et al.Combustion system development for the ramgen engine[J].Journal of Engineering for Gas Turbines and Power，2003(7):885－894.

[3]Edmonds R G，Steele R C，Williams J T，et al.Ultra-low NOxadvanced vortex combustor[C]//Proceedings of the ASME Turbo Expo 2006，Power for Land，Sea and Air，2006:255－262.

[4]鄧洋波，劉世青，鐘兢軍.AVC中鈍體布置與燃燒室流動特性研究[J].工程熱物理學報，2008，29(8):1416－1418.

[5]鄧洋波，劉世青，鐘兢軍.先進旋渦燃燒室燃燒特性數值模擬[J].大連海事大學學報，2008，34(3):21－24.

[6]Roquemore W，Shouse D，Burrus D，et al.Trapped vortex combustor concept for gas turbine engines[C]//39th AIAA Aerospace Sciences Meeting and Exhibit，Reno，Nevada，2001，0483.

[7]鄧洋波，劉世青，鐘兢軍.先進旋渦燃燒室流動與燃燒特性分析[J].航空動力學報，2009，24(3):488－493.

[8]鐘兢軍，劉世青.后駐體噴孔位置對駐渦腔流動冷態數值的影響[J].上海海事大學學報，2011，32(1):44－48.