楊　祺，張文瑞，于錕錕

（蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室，蘭州　730000）

基于FLUENT的紅外窗口主動(dòng)冷卻裝置數(shù)值模擬分析

楊祺，張文瑞，于錕錕

（蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室，蘭州730000）

處于氣動(dòng)加熱環(huán)境下的紅外窗口，由于強(qiáng)熱流影響會(huì)導(dǎo)致窗口材料外形和光學(xué)性質(zhì)發(fā)生變化，嚴(yán)重影響光學(xué)窗口的正常工作。以紅外成像在高超聲速飛行器上的應(yīng)用為需求背景，針對(duì)氣動(dòng)加熱條件下紅外窗口的傳熱特性，建立紅外窗口三維非穩(wěn)態(tài)傳熱計(jì)算模型，利用FLUENT軟件對(duì)紅外窗口采用主動(dòng)冷卻與被動(dòng)熱防護(hù)效果進(jìn)行了數(shù)值模擬，驗(yàn)證了主動(dòng)冷卻技術(shù)在紅外窗口溫度控制的可行性，并對(duì)所需的冷量進(jìn)行了計(jì)算，給出了典型狀態(tài)下的計(jì)算結(jié)果。

熱防護(hù)；氣動(dòng)加熱；紅外窗口；數(shù)值模擬

0　引言

紅外窗口是紅外自動(dòng)尋的導(dǎo)彈等高超聲速飛行器的重要部件，飛行期間要經(jīng)受較嚴(yán)重的氣動(dòng)加熱［1］。在飛行高度為10 km，飛行速度4馬赫的情況下，彈體頭部駐點(diǎn)的溫度將達(dá)到930 K左右。氣動(dòng)加熱產(chǎn)生的高溫和窗口附近復(fù)雜的大氣環(huán)境會(huì)對(duì)紅外成像導(dǎo)引頭的末制導(dǎo)系統(tǒng)產(chǎn)生嚴(yán)重的熱干擾［2］，造成導(dǎo)引頭成像模糊、抖動(dòng)，甚至?xí)?dǎo)致目標(biāo)信號(hào)被紅外窗口的熱輻射所淹沒(méi)，為保證光學(xué)窗口正常工作并提高其工作能力，需采用一種能實(shí)現(xiàn)光學(xué)窗口冷卻的高效穩(wěn)定的冷卻技術(shù)。

熱防護(hù)系統(tǒng)（TPS）［3］是保證高超聲速飛行器安全、快捷飛行的關(guān)鍵技術(shù)之一。熱防護(hù)系統(tǒng)一般位于飛行器外部，是保護(hù)內(nèi)部結(jié)構(gòu)或者其它子系統(tǒng)溫度在使用的溫度范圍，同時(shí)承載一定力學(xué)載荷的結(jié)構(gòu)系統(tǒng)。熱防護(hù)系統(tǒng)的功能主要是保護(hù)飛行器內(nèi)部結(jié)構(gòu)和有效載荷的溫度不要超過(guò)設(shè)計(jì)使用溫度。

高超聲速飛行器的熱防護(hù)系統(tǒng)主要分為主動(dòng)熱防護(hù)系統(tǒng)、半主動(dòng)熱防護(hù)系統(tǒng)以及被動(dòng)熱防護(hù)系統(tǒng)三大類。在被動(dòng)熱防護(hù)方案中，熱量由飛行器表面輻射出去或被吸收，不需要工作流體（工質(zhì)）來(lái)排除。半主動(dòng)熱防護(hù)系統(tǒng)介于被動(dòng)和主動(dòng)冷卻方案之間，大部分的熱量靠工作流體或空氣流帶走。主動(dòng)冷卻熱防護(hù)系統(tǒng)采用特殊防熱結(jié)構(gòu)和材料，利用冷卻工質(zhì)（固態(tài)、液態(tài)、氣態(tài)）阻止或帶走氣動(dòng)熱，盡可能降低氣動(dòng)加熱對(duì)飛行器的影響，保護(hù)內(nèi)部結(jié)構(gòu)和其他系統(tǒng)的溫度不超過(guò)設(shè)定溫度。

1　主動(dòng)冷卻與被動(dòng)熱防護(hù)技術(shù)原理

1.1工作原理

主動(dòng)冷卻與被動(dòng)熱防護(hù)技術(shù)原理如圖1所示，使用隔熱結(jié)構(gòu)阻隔來(lái)自于外部高速氣流的熱量，在隔熱結(jié)構(gòu)升溫過(guò)程中，熱防護(hù)系統(tǒng)采用相變制冷原理，使存貯在飛行器中的制冷工質(zhì)經(jīng)換熱器相變后作用于被冷卻窗口，冷卻紅外窗口及附近結(jié)構(gòu)。

圖1　熱防護(hù)系統(tǒng)原理圖

1.2相變制冷原理［4-5］

利用制冷劑本身的飽和蒸氣壓作為系統(tǒng)工作動(dòng)力，不需要配置額外的壓力源。經(jīng)過(guò)節(jié)流組件節(jié)流后的制冷劑進(jìn)入窗口周圍換熱器，在換熱器中相變蒸發(fā)吸熱，產(chǎn)生冷量實(shí)現(xiàn)紅外窗口降溫控溫。原理圖如圖2所示。

圖2　開(kāi)式節(jié)流制冷原理圖

2　光學(xué)窗口傳熱模型的建立

2.1計(jì)算模型

針對(duì)紅外窗口在某典型飛行器彈道條件下，氣動(dòng)加熱時(shí)的受熱情況，對(duì)采取基于相變制冷的主動(dòng)冷卻及被動(dòng)隔熱防護(hù)措施時(shí)所需的冷量進(jìn)行了計(jì)算。采用CFD/NHT商用軟件ANSYS13.0-FLUENT進(jìn)行計(jì)算。非穩(wěn)態(tài)項(xiàng)采用二階隱式格式，時(shí)間步長(zhǎng)為0.1 s，每時(shí)間步長(zhǎng)迭代次數(shù)為20次。計(jì)算使用的三維物理模型如圖3所示。

圖3　計(jì)算模型結(jié)構(gòu)剖面圖

控制方程為：

其中：keff為有效導(dǎo)熱率（k+keff，kf為湍流引致的導(dǎo)熱率，由模型中使用的湍流模型確定）；為組分j的擴(kuò)散通量。式（1）右邊的前三項(xiàng)分別表示由于熱傳導(dǎo)、組分?jǐn)U散、粘性耗散而引起的能量轉(zhuǎn)移。Sh包含化學(xué)反應(yīng)放（吸）熱以及任何其他的由用戶定義的體積熱源式（2）：

其中：顯焓h的定義（對(duì)理想氣體）為式（3）：

對(duì)不可壓流體：

式（2）、（3）中Yj為組分 j的質(zhì)量分?jǐn)?shù)：

其中：Tref為298.15 K。

2.2網(wǎng)格劃分

采用gambit軟件對(duì)網(wǎng)格進(jìn)行了劃分。為了生成高質(zhì)量計(jì)算網(wǎng)格，模型網(wǎng)格生成采用分塊結(jié)構(gòu)化網(wǎng)格生成技術(shù)，為了準(zhǔn)確計(jì)算溫度場(chǎng)，在計(jì)算域中的網(wǎng)格進(jìn)行了局部加密，以滿足其要求。網(wǎng)格數(shù)約為500 W，通過(guò)與多種網(wǎng)格數(shù)的計(jì)算結(jié)果對(duì)比，已驗(yàn)證計(jì)算結(jié)果是網(wǎng)格無(wú)關(guān)解。

2.3邊界條件設(shè)定

建立三維數(shù)學(xué)模型，模擬紅外窗口氣動(dòng)加熱條件下保證窗口溫度在一定范圍內(nèi)所需冷量。計(jì)算中選用SiO為窗口材料［6］，隔熱結(jié)構(gòu)選用氣凝膠，加熱面采用某典型氣動(dòng)加熱邊界條件，非加熱面采用對(duì)流冷卻邊界條件，艙內(nèi)流體溫度為333 K，對(duì)流換熱系數(shù)取為30 W/（m2·K），光學(xué)窗口初始溫度293 K，不考慮環(huán)境輻射影響，計(jì)算中用到的材料物性如表1所列。

表1　計(jì)算所需物性參數(shù)

3　計(jì)算結(jié)果及分析

3.1無(wú)冷卻條件下傳熱計(jì)算

在無(wú)冷卻條件下對(duì)紅外窗口1 000 s內(nèi)的傳熱進(jìn)行了計(jì)算。無(wú)冷卻時(shí)，模型溫度分布云圖及紅外窗口溫度如圖4、圖5所示，從圖中可以看出，隨著時(shí)間的推移，窗口溫度迅速升高，在1 000 s后，窗口溫度達(dá)到430 K左右。由圖4可以看出，由于結(jié)構(gòu)原因，窗口溫度分布不均勻。

圖4　無(wú)冷卻條件下窗口模型溫度分布云圖

圖5　無(wú)冷卻條件下紅外玻璃窗口溫度曲線圖

3.2冷量為250 W條件下的傳熱計(jì)算

通過(guò)主動(dòng)冷卻方式對(duì)窗口施加250 W的冷量后，根據(jù)模型進(jìn)行計(jì)算，得到非穩(wěn)態(tài)熱流傳導(dǎo)特性的結(jié)果，如圖6、圖7所示。

圖7　250 W冷量冷卻下紅外窗口溫度曲線圖

從圖7中可以看出，對(duì)冷卻結(jié)構(gòu)施加250 W冷量后，在前200 s內(nèi)，窗口溫度上升平緩，甚至有下降趨勢(shì)，到200 s后，窗口溫度逐漸升高，1 000 s后，窗口平均溫度約為370 K，可見(jiàn)施加冷量后，窗口的溫度得到了明顯的控制。從圖6中可以看出，被冷卻后的窗口溫度分布均勻性得到了很大提升。

3.3冷量為500 W條件下的傳熱計(jì)算

圖8、圖9為對(duì)冷卻結(jié)構(gòu)施加500 W冷量時(shí)，紅外窗口的溫度分布云圖及升溫曲線。從圖中可以看出，紅外窗口的溫度出現(xiàn)兩個(gè)拐點(diǎn)，在約220 s后溫度緩慢上升，最終1 000 s后紅外窗口的溫度被控制在335 K左右。窗口的溫度分布均勻性比250 W冷量時(shí)得到進(jìn)一步改善。

圖8　500 W冷量冷卻下窗口模型溫度分布云圖

圖9　500 W冷量冷卻下紅外窗口溫度曲線圖

3.4熱模擬分析結(jié)果

在無(wú)冷卻、加冷量分別為250 W、500 W的條件下，1 000 s內(nèi)紅外窗口平均溫度曲線對(duì)比圖如10所示。

圖10　三種條件下1 000 s內(nèi)紅外窗口平均溫度對(duì)比圖

從圖10中可以看出，通過(guò)主動(dòng)冷卻裝置給紅外窗口施加冷源后，窗口的溫度得到了有效的控制，相較于無(wú)冷源和250 W冷源的情況，在1 000 s內(nèi)，施加500 W冷源時(shí)窗口溫度上升趨勢(shì)較為平緩，窗口平均溫度滿足要求，窗口溫度均勻性好。

4　結(jié)論

通過(guò)主動(dòng)冷卻與被動(dòng)熱防護(hù)相結(jié)合，在1 000 s內(nèi)，對(duì)窗口區(qū)域施加500 W的冷量可將紅外窗口的溫度控制在要求的范圍內(nèi)。

［1］耿湘人，桂業(yè)偉，賀立新，等，紅外窗口不同冷卻方式下的結(jié)構(gòu)傳熱和熱應(yīng)力特性計(jì)算研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2008.26（3）：329-333.

［2］范志剛，張亞萍，裴揚(yáng)威.氣動(dòng)熱環(huán)境下高速飛行器光學(xué)窗口光傳輸數(shù)值仿真研究［J］.紅外與毫米波學(xué)報(bào)，2007，26（5）：396-400.

［3］劉雙.高超聲速飛行器熱防護(hù)系統(tǒng)主動(dòng)冷卻機(jī)制與效能評(píng)估［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2010.

［4］陳正剛，王小軍，朱恩寶，等.開(kāi)放式制冷系統(tǒng)與應(yīng)用前景分析［J］.真空與低溫，2011，17（z1）：622-625.

［5］唐小偉.開(kāi)式節(jié)流制冷技術(shù)在局部環(huán)境溫度控制中的應(yīng)用研究［J］.真空與低溫，2010，16（4）：223-226.

［6］李東輝，夏新林，艾青.氣動(dòng)加熱下光學(xué)窗口傳熱特性的數(shù)值模擬［J］.宇航材料工藝，2012（2）：30-33.

NUMERICAL SIMULATION ON INITIATIVE COOLING OF OPTICAL WINDOWS BASED ON FLUENT

YANG Qi，ZHANG Wen-rui，YU Kun-kun
（Science and Technology on Vacuum Technology and Physics Laboratory，Lanzhou Institute of Space Technology and Physics，LanZhou730000，China）

In the aero-dynamic environment，optical windows shape and optical property will be changed because of high heat flow，then the employment of the transmission of optical windows will be effected.According to the requirement of IR imaging technique applied in hypersonic aircraft，the heat transfer characteristic in optical windows under aerodynamic heating were studied and the three dimensions unsteady state heat conduction equations in unstructured were established.The effect of initiative cooling and passive heat insulation thermal protection technology were calculated with the engineering prediction software-Fluent.The feasibility of initiative cooling technology on temperature control of optical windows were validated.The capacity were caculate and some computer simulation results at classical state were presented.

thermal protection；aerodynamic heating；Optical Windows；Numerical Simulation

TB657

A

1006-7086（2015）02-0092-04

10.3969/j.issn.1006-7086.2015.02.007

2014-08-26

真空低溫技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室基金項(xiàng)目（9140C550211120C5501）

楊祺（1983-），女，甘肅省蘭州人，碩士，從事低溫制冷方面的研究。E-mail：qiqi94520@126.com。