化學激光器燃燒室傳熱過程分析

王杰，尹燕華，郭建增，顏飛雪

（1．中國船舶重工集團公司第七一八研究所，河北邯鄲056027; 2．南京理工大學瞬態物理國家重點實驗室，江蘇南京 210094）

0 引言

燃燒室是HF/DF化學激光器中工作溫度最高的部件，其用途在于為光腔泵浦反應提供足夠的F原子。HF/DF激光器燃燒室溫度必須足夠高（大于1 500 K），以保證過量含氟氧化劑（NF3或SF6等）充分解離分解成自由F原子，且不會復合成為F2分子。燃燒室內壁在高溫燃氣作用下，可能會發生變形、裂紋、皺曲等故障。為了降低燃燒室壁面溫度和溫度梯度，往往對燃燒室壁面采取有效的冷卻措施［1－4］。冷卻措施的加入，必然造成相當熱量的損失，化學激光器燃燒室設計要求在熱量損失最小的前提下保證激光器正常高效運轉。Robert［2］等的研究表明，降低化學激光器燃燒室熱量損失能提高激光器效率。他們對減小水冷HYLTE化學激光器模塊熱量損失方法進行了探討，設計了新型低熱損失燃燒室。然而，化學激光器燃燒室熱量損失的最小化與激光器燃燒室基本結構、燃燒室熱量的產生及傳遞規律有關。本文在實驗測量的基礎上，通過理論分析對化學激光器燃燒室傳熱過程和傳熱模型進行了分析和討論。

1 燃燒室熱量的產生

圖1 化學激光器燃燒室結構示意圖Fig.1Schematic of chemical laser combustor

燃燒室是HF/DF化學激光器的重要組成部分，其結構如圖1所示。燃燒室上游與噴注器連接，下游與噴管列陣相連。噴注器結構復雜，其作用包含2個方面:一是用來注入燃燒室燃燒反應的反應物（氧化劑、燃料和稀釋劑）;二是實現燃燒室的點火過程，觸發燃燒反應的進行。燃燒段構成化學激光器燃燒室側壁。噴管列陣主要是將燃燒室產生的含氟原子氧化劑流和參與光腔反應的燃料氣快速充分混合，使燃燒室產生的含氟原子氧化劑氣流冷卻至低溫，使氟原子濃度有效“凍結”［5］。

化學激光器燃燒室熱量的產生來自于燃燒反應釋放的熱能。對于F2－D2－He體系的HF化學激光而言，燃燒室燃燒反應是

對NF3－C2H4－He體系的DF化學激光而言，燃燒室燃燒反應是

式中，Q為燃燒反應釋放的熱量。2個反應中氧化劑（F2和NF3）是過量的，燃燒反應釋放的熱量取決于燃料氣的摩爾流率。燃燒反應產生的熱量，一部分用來加熱燃氣使之達到穩定的滯止溫度，一部分傳遞給燃燒室壁損失了，其余部分都用來分解解離過量含氟氧化劑，產生自由氟原子。燃燒室的滯止溫度可以通過加入的稀釋劑和過量的含氟氧化劑的量來調節控制。

2 傳熱模型及分析

圖2是化學激光器燃燒室壁面傳熱示意圖。燃燒室的內壁面與高溫燃氣有對流換熱，其熱流通量為qhc。燃燒室的內壁面與高溫燃氣有輻射換熱，其熱流通量為qhr。燃燒室壁面內有軸向、徑向和周向的導熱，但在實際中經常忽略軸向和周向的導熱，只考慮沿壁厚（徑向）方向的導熱，其熱量通量為qr。燃燒室的外壁面與外界環境接觸，有自然對流，其熱流通量用qcc表示。燃燒室外壁面向外界環境的輻射換熱熱流通量為qcr。在熱平衡（穩態）條件下，燃燒室壁面的熱平衡方程可表示為

式中各項均與燃燒室壁面溫度有關。

圖2 燃燒室壁面傳熱示意圖Fig.2Schematic of combustor heat transfer

如圖2所示，燃燒室燃燒反應產生的高溫燃氣絕熱溫度為Thg，燃燒室的內壁面溫度為Thw，外壁面溫度為Tcw，燃燒室的壁面初始溫度與外界環境溫度為T0。沿燃燒室壁面厚度方向的導熱熱流通量計算可近似用平壁公式，即

式中:λw為燃燒室壁的導熱系數;δ為燃燒室壁厚度。

燃燒室的內壁面對流換熱與燃氣溫度、壓力、燃氣流動速度及燃氣成分等參數有關，可根據下式計算:

式中，h為燃燒室內壁面的對流換熱系數。由于激光器燃燒室燃氣溫度較高，流動速度相對較大，燃氣的流動表現為湍流流動的特點，其對流傳熱為強迫對流傳熱。傳熱系數計算可用下面相似準則［6］計算:

按努塞爾數、雷諾數和普朗特數的定義展開，有

整理，即得強迫對流傳熱系數

式中:cp為定壓比熱容;λ為導熱系數;μ為粘性系數;m˙為燃氣質量流率;A為燃氣流通截面積;d為特征尺寸。

由于構成激光器燃燒室的噴注器、燃燒段和噴管列陣的結構各不相同，燃氣與之接觸作用的方式也各不相同。因此在計算對流換熱系數時，往往需要采用不同的修正方法對換熱系數的計算進行修正，如噴管內壁面的對流換熱系數經常采用巴茲公式進行計算［2］。因此，采用式（4）求解得到的對流換熱系數只是局部對流換熱系數。若令

式中:heff為有效換熱系數;hi為各局部換熱系數。可將式（3）改寫為

燃燒室的外壁面與空氣的對流傳熱為空氣自由對流傳熱，其傳熱系數可以采用格拉曉夫準則或瑞利準則進行計算［6］。高溫燃氣與燃燒室的內壁面的輻射傳熱可用下式計算:

式中，hhr為輻射換熱系數，滿足

其中:σ為斯忒藩－波耳茲曼常數;εhg為高溫燃氣的發射率;αhg為高溫燃氣對內壁面輻射吸收率;ε'w為灰體壁面的有效發射率。

燃燒室的外壁面與外界環境的輻射換熱也可以用相同的方法求解。對于式（5），實際上是采用對流換熱的形式來表達輻射換熱。將式（2）～式（5）代入式（1），即可對激光器燃燒室傳熱過程及相關參數進行求解。化學激光器的內壁面溫度通常需要進行迭代求解，先假設1個內壁面溫度，迭代計算直至結果收斂。

然而，化學激光器在短時間運行中，壁面的導熱過程來不及達到穩定，激光器就結束工作了。圖3是化學激光器運行5 s時燃燒室外壁面溫度隨時間變化曲線圖，圖中的2條曲線分別表示2次不同實驗測量結果。由圖可知，在激光器運行的5 s內，燃燒室外壁面溫度只有微弱的升高，即使在激光器運行結束后的5 s時，外壁面溫度仍然沒有達到穩定。整個傳熱過程屬于非穩態傳熱過程，燃燒室外壁面溫度梯度逐漸增大。

激光器燃燒室在非穩態傳熱情況下，不能再采用式（1）所示的穩態傳熱模型進行傳熱及相關參數計算。由圖3可知，在激光器運行的5 s內，可近似認為無熱量自燃燒室外壁面散出，即qcc和qcr等于0。對激光器燃燒室而言，處于絕熱狀態;對燃燒室內的高溫燃氣而言，熱量損失是由燃燒室壁面吸熱造成。在這種情況下，激光器燃燒室壁可當作無限大平板的一半進行處理。采用經典的分離變量法求解無內熱源的一維不穩定導熱微分方程，即可得到燃燒室壁內的溫度分布:

圖3 燃燒室外壁面溫度變化曲線Fig.3Temperature curve of combustor ektexine

式中:θi為初始過余溫度;Fo為傅立葉數，滿足

其中:λ為燃燒室壁面材料導熱系數;ρ為材料密度; c為比熱容;t為燃燒室運行時間;βn為超越方程Bi =βntanβn的根，其中Bi為畢渥數。

由式（6）可以計算任一時刻燃燒室壁內的溫度分布。取x=0和x=δ，就可得到燃燒室內壁面和外壁面的過余溫度，分別為:

在工程計算中，對式（7）和式（8）中的無窮級數，取級數的第1項就足夠準確。當Fo≥0.3時，計算誤差不超過1%。因此，燃燒室內、外壁面過余溫度可表示為:

對圖3所示的激光器燃燒室非穩態問題，應用式（9）～（11）可計算出燃燒室內、外壁面溫度，燃燒室沿壁面平均溫度，即

燃燒室壁面吸熱量為:

通常材料的基本熱物性參數（如比熱容、導熱系數等）是隨溫度變化而變化的，但在工程計算中可以將這些物性參數看作是不隨溫度變化的常數。采用上面的傳熱模型及公式就可以對燃燒室傳熱過程與相關參數進行計算。在化學激光器中，構成燃燒室的各壁面形狀不是如圖2所示的“理想”壁面，其形狀結構相對比較復雜（如構成燃燒室各壁面內的復雜氣流通道等），傳熱模型中的一些相關參數也是隨時間、溫度和位置等變化的，因此要準確有效地計算化學激光器燃燒室的傳熱過程及相關參數，還需要大量的實驗結果對基本的傳熱模型進行修正和驗證。

3 結語

通過理論分析和試驗對比，建立了化學激光器燃燒室基本傳熱模型。激光器長時間運行過程中，燃燒室熱量傳遞達到穩定狀態時，高溫燃氣對燃燒室內壁的輻射傳熱和對流傳熱熱流通量之和與燃燒室壁面導熱熱流通量相等，與燃燒室外壁面的輻射傳熱和對流傳熱熱流通量之和相等。采用迭代方法可以對燃燒室傳熱過程及相關參數（燃燒室內壁面溫度等）進行計算。激光器在短時間運行過程中，燃燒室壁面的導熱過程來不及達到穩定，激光器就結束工作了，燃燒室傳熱過程為非穩態傳熱過程。對激光器燃燒室而言，處于絕熱狀態;對燃燒室內的高溫燃氣而言，熱量損失是由燃燒室壁面吸熱造成。通過將燃燒室壁轉化為無限大平板的一半的導熱問題，即能對激光器燃燒室非穩態傳熱過程進行計算。在實際的化學激光器中，構成燃燒室各壁面的形狀各不相同，高溫燃氣與各壁面作用過程也不相同，因此在接下來的工作中，將進一步通過實驗對化學激光器燃燒室穩態或非穩態傳熱過程及傳熱模型進行檢驗和修正。

［1］NAGAI C K，CARLSON L W，GIEDT R R，et al．Heattransfer characteristics of an advanced DF chemical laser［R］．AIAA Paper:74－684．

［2］WALDO R E，BETTS J A，GRAVES B R，et al．Approach for reduced heat loss in HYLTE laser module［R］．AIAA Paper，1997:1－17．

［3］唐力鐵，卓海濤，陸啟生．DF激光器燃燒室反應流場的一種新的計算求解方法［J］．強激光與粒子束，2005，17 （11）:1625－1629．

TANG Li-tie，ZHUO Hai-tao，LU Qi-sheng．Numerical solution to multi-component chemical equilibrium flows in DF chemical laser combustion chamber［J］．High Power Laser and Particle Beams，2005，17（11）:1625－1629．

［4］PATTERSON S，GRAVESB，BETTSJ，etal．HF fundamental and overtone testing using the HYLTE water cooled module［R］．AIAA Paper，1999－3472．

［5］格羅斯R W F，博特J F．化學激光手冊［M］．北京:科學出版社，1987．

GROSS R W F，BOTT J F．Handbook of chemical lasers［M］．Beijing:Science Press，1987．

［6］鄭亞，陳軍，鞠玉濤，等．固體火箭發動機傳熱學［M］．北京:北京航空航天大學出版社，2006．

ZHENG Ya，CHEN Jun，JU Yu-tao，et al．Heat transfer of solid rocket engine［M］．Beijing:Beijing University of Aeronautics and Astronautics Press，2006．