復雜熱環境下藍寶石共形頭罩氣動熱效應研究*

肖昊蘇，陳澄，姚睿，魏宇飛，張麗琴

(北京電子工程總體研究所，北京 100854)

復雜熱環境下藍寶石共形頭罩氣動熱效應研究*

肖昊蘇，陳澄，姚睿，魏宇飛，張麗琴

(北京電子工程總體研究所，北京 100854)

分析了藍寶石共形頭罩外流場特性，獲得了頭罩溫度、應力、應變和形變場；根據熱光效應和彈光效應計算了頭罩折射率分布。研究發現復雜熱環境下共形頭罩駐點處和前方區域氣動加熱較為嚴重，熱光效應對頭罩折射率場改變的影響遠比彈光效應大。研究成果對認識共形頭罩氣動加熱的機理以及頭罩防熱設計具有一定指導意義。

氣動光學窗口；熱輻射；光線追跡；建模研究

0 引言

隨著高速飛行器速度的不斷提高，飛行器窗口的氣動特性越來越受到重視。因為共形頭罩具有較好的氣動特性，所以共形頭罩將在高速飛行器上的應用越來越廣泛。

美國早在20世紀90年代就開始了共形頭罩的研究。美國國防部高級研究計劃局自1996年資助并發展共形光學項目[1]。美國Raytheon公司在2000年就實現了共形光學系統成像[2]。2001年，美國成功地將共形頭罩應用于導彈系統中，使得導彈的空氣動力降低了25%[2]。國內對共形光學的研究尚處于起步階段，但是已經認識到了共形光學系統的重要性。中航工業洛陽電光設備研究所張鵬等[3]，丁全心等[4]，中科院長春光機所李東熙[5-6]、姜振海[7]、魏群[ 8]、王超[9]等，兵器工業209所黃秋等[10]，北京理工大學李林[11]、常軍等[12]，也開展了共形光學系統設計研究。

目前，共形頭罩光學系統設計[1-12]方面的研究已經很多，而氣動熱[13-14]對共形頭罩的作用，及其對其折射率場的影響，卻鮮有報道。本文分析了氣動熱環境下共形頭罩的溫度場、應變場和形變場的分布規律，根據熱光效應和彈光效應的理論計算了共形頭罩的折射率場分布，得到的結果可以用于分析氣動熱環境下共形頭罩的光傳輸特性。

1 共形頭罩外流場特性仿真研究

采用Spalart-Allmaras湍流模型分析共形頭罩的外流場特性：

(1)

式中：Yν為組分ν的質量分數；S為化學反應引起的放熱和吸熱；u為來流速度。

Spalart-Allmaras模型是相對簡單的單方程模型，主要用于求解有關渦粘性的運輸方程。利用雷諾類推概念，建立湍流傳質和傳熱模型：

ui(τij)eff]+Sh，

(2)

式中：(τij)eff為粘性發熱。

高速飛行器共形頭罩頂點的曲率半徑為50 mm，長徑比為1.0。根據高速飛行器幾何形狀參數，建立了高速飛行器外流場CFD(computed fluid dynamics)計算網格模型，根據高速飛行器的飛行軌跡設置了外流場計算的邊界條件，最終獲得了不同飛行狀態下高速飛行器共形頭罩窗口外流場密度場以及共形頭罩窗口外表面壁面靜溫、靜壓和熱流密度分布，分別如圖1～4所示。高速飛行器共形頭罩外表面靜溫場，靜壓場和熱流密度場從前端點沿對稱軸方向逐漸減小，駐點或者頭罩前方的氣動加熱現象最為嚴重。

圖1 共形頭罩外流場平均密度分布(kg/m3)Fig.1 Average density distribution of the external flow field around the conformal dome (kg/m3)

2 共性頭罩氣動熱響應仿真研究

將共形頭罩外表面靜壓、靜溫和熱流密度分布作為頭罩熱響應計算的邊界條件，采用熱-結構直接耦合分析的方法，利用有限元軟件ANSYS獲得了藍寶石晶體頭罩在末制導最后時刻的溫度場、等效應變場和等效形變場分布，分別如圖5～8所示。頭罩具有最高溫度分布的區域與外表面具有最大熱流密度分布的區域一致，且頭罩溫度最高可以升至701.41 K。為了使探測器不發生飽和， 頭罩的溫度必須控制在450 K以下[15]。必須采用致冷措施來降低頭罩窗口的溫度。隨著頭罩溫度不斷增加，頭罩等效形變也不斷的增加，且頭罩最大等效形變為372 μm。隨著頭罩溫度和等效形變的不斷增加，頭罩等效應變也不斷增加，且頭罩最大等效應變為4.34×10-3。

圖2 共形頭罩外表面壁面靜溫分布(K)Fig.2 Static wall temperature distribution on the outside surface of the conformal dome (K)

圖3 共形頭罩外表面壁面靜壓分布(Pa)Fig.3 Static pressure distribution on the outside surface of the conformal dome (Pa)

圖4 共形頭罩外表面壁面熱流密度分布(W/m2)Fig.4 Heat flux distribution on the outside surface of the conformal dome (W/m2)

圖5 末制導最后時刻共形頭罩溫度場(K)Fig.5 Temperature filed of the conformal dome at the last moment of terminal guidance (K)

圖6 末制導最后時刻共形頭罩等效形變場(m)Fig.6 Sum deformation field of the conformal dome at the last moment of terminal guidance (m)

圖7 末制導最后時刻共形頭罩等效應變場Fig.7 Equivalent von mises strain field of the conformal dome at the last moment of terminal guidance

3 氣動熱環境下共形頭罩折射率場仿真研究

將頭罩熱響應分析得到的溫度場、熱應力場和熱應變場作為輸入條件，根據藍寶石晶體材料的熱光系數和彈光系數，可計算藍寶石晶體頭罩在氣動熱環境下折射率變化。

3.1 熱光效應

熱光效應是指光學介質的光學性質隨溫度變化發生變化。折射率隨溫度變化的關系為[16]

(3)

3.2 彈光效應

由于藍寶石單軸晶體屬于三方晶系，晶軸方向為Ox1,Ox2,Ox3，折射率橢球為一旋轉橢球體：

(4)

(5)

藍寶石晶體屬于三方晶系，三方晶系的彈光系數矩陣可以寫成[16]：

(6)

由式(6)可以得出：

(7)

綜合式(4)～(7)可知，當藍寶石單晶受到熱應變作用時，折射率橢球不再是旋轉橢球，且主軸方向發生改變，晶體不再是單軸晶體，而是變為雙軸晶體，為了簡化，只考慮3個主軸方向的折射率的變化，變化后的3個主軸方向的折射率n1′,n2′,n3′分別為

(8)

藍寶石晶體為單軸晶體，因此它具有o光折射率no和e光折射率ne。藍寶石單晶在296 K溫度下，4 μm波段o光和e光的折射率分別為no=1.677 0，ne=1.667 0，o光熱光系數為1.139×10-5K-1，e光熱光系數為1.282×10-5K-1。藍寶石晶體的彈光系數為P11=0.2，P12=0.078，P13=0，P14=0.03，P31=0，P33=0.252,P41=0.02和P44=0.09。

3.3 藍寶石晶體頭罩折射率場分析

受熱光效應的影響，末制導最后時刻藍寶石頭罩外內表面中心點的折射率改變量如表1所示。受彈光效應的影響，末制導最后時刻藍寶石頭罩外內表面中心點折射率改變量如表2所示。由表1和表2可知：藍寶石頭罩外表面中心點由彈光效應引起的折射率改變量最大值(2.7×10-4)僅為熱光效應引起的折射率改變量最小值(4.3×10-3)的6.3%；藍寶石頭罩內表面中心點由彈光效應引起的折射率改變量最大值(-6.5×10-5)僅為熱光效應引起的折射率改變量最小值(2.8×10-3)的6.1%。因此，熱光效應對藍寶石共形頭罩窗口折射率場分布的影響遠比彈光效應的影響大，熱光效應對藍寶石共形頭罩窗口光傳輸的影響也將遠比彈光效應的影響大。

表1 藍寶石頭罩外內表面中心點由熱光效應引起的折射率改變量Table 1 Refractive index variations of the central points on the outside and inside surfaces of sapphire dome caused by the thermo-optical effect

表2 藍寶石頭罩外內表面中心點由彈光效應引起的折射率改變量

4 結束語

采用Spalart-Allmaras湍流模型分析了高速飛行器藍寶石共形頭罩外流場特性；利用熱-結構直接耦合分析的方法獲得了高速飛行器藍寶石共形頭罩溫度場、應變場和形變場，根據熱光效應和彈光效應的原理計算了共形頭罩外內表面中心點折射率變化情況。通過分析可知，氣動熱環境下共形頭罩窗口溫度場、應變場和形變場分布是不均勻的，受熱光效應和彈光效應的影響，共形頭罩窗口的折射率場分布也是非均勻的，且熱光效應對藍寶石共形頭罩窗口折射率場分布的影響遠比彈光效應的影響大。本文的研究成果對認識高速飛行器共形頭罩氣動加熱的機理和防熱設計具有一定的指導意義。

[1] KNAPPA J, MILLS, P, HEGG G, et al. Conformal Optics Risk Reduction Demonstration [C]∥SPIE 2001, 4375: 146-153.

[2] TROTTA A. Precision Conformal Optics Technology Program [C]∥SPIE 2001, 4375: 96-107.

[3] 張鵬, 趙春竹, 崔慶豐. 航空共形光學窗口的設計方法[J]. 光學學報, 2013,33(6):1-6. ZHANG Peng, ZHAO Chun-zhu, CUI Qing-feng. Design Method of an Aerial Conformal Optical Window[J]. Acta Optica Sinica, 2013，33(6)：1-6.

[4] 丁全心, 熊鐘秀. 應用于高速飛行器的共形光學系統設計研究[J]. 紅外與激光工程, 2010, 39(S): 1-4. DING Quan-xin, XIONG Zhong-xiu.Study on Conformal Optical System Design for High Speed Aircraft[J]. Infrared and Laser Engineering, 2010, 39(S): 1-4.

[5] 李東熙, 盧振武, 孫強, 等. 基于實際光線追跡的共形光學系統設計[J]. 紅外與激光工程, 2008, 37(5): 834-838. LI Dong-xi, LU Zhen-wu, SUN Qiang, et al. Design of Conformal Optical System based on Real Ray Tracing[J]. Infrared and Laser Engineering, 2008, 37(5): 834-838.

[6] 李東熙, 盧振武, 孫強, 等. 利用Wassermann-Wolf原理設計共形光學系統[J]. 光子學報, 2008, 37(4): 776-779. LI Dong-xi, LU Zhen-wu, SUN Qiang, et al. Research on Conformal Optical System Design Using Wassermann-Wolf Principle[J]. Acta Photonica Sinica, 2008, 37(4): 776-779.

[7] 姜振海, 王超, 魏群, 等. 超音速共形整流罩風洞試驗及其光機特性[J]. 光學精密工程, 2012, 20(9): 1999-2005. JIANG Zhen-hai, WANG Chao, WEI Qun, et al. Wind Tunnel Experiment of Supersonic Conformal Dome and Its Optical and Structure Characteristics[J]. Optics and Precision Engineering, 2012, 20(9): 1999-2005.

[8] 魏群, 王超, 姜湖海, 等. 中波紅外彈載共形光學系統研制[J]. 紅外與激光工程, 2013, 42(5): 1298-1301. WEI Qun, WANG Chao, JIANG Hu-hai, et al. Manufacture of Mid-Wave Infrared Missile Used Conformal Optical System[J]. Infrared and Laser Engineering, 2013, 42(5): 1298-1301.

[9] 黃秋, 陳亦慶, 高志峰. 共形光學[J]. 光學技術, 2009, 35(5): 729-735. HUANG Qiu，CHEN Yi-qing,GAO Zhi-feng. Coformal Optics Technology[J]. Optical Technique, 2009, 35(5): 729-735.

[10] 李林, 劉家國, 李卓, 等. 共形光學設計研究[J]. 光學技術, 2006, 32(S): 509-515. LI Lin, LIU Jia-guo, LI Zhuo, et al. Design and Study of Conformal Optics[J]. Optical Technique, 2006, 32(S): 509-515.

[11] 常軍, 馮萍, 何伍斌. 共形光學系統設計研究[J]. 航空科學技術, 2013(3):76-78. CHANG Jun, FENG Ping, HE Wu-bin. The Research of Conformal Optical System Design[J]. Aeronautical Science & Technology, 2013(3):76-78.

[12] 王超, 姜湖海, 朱瑞飛, 等. 共形光學導引頭瞄視誤差分析與修正[J]. 光學學報, 2013, 33(9):1-7. WANG Chao, JIANG Hu-hai, ZHU Rui-fei, et al. Boresight Error Analysis and Correction of Conformal Infrared Seeker[J]. Acta Optica Sinica, 2013, 33(9):1-7.

[13] 張勝濤, 張慶兵, 王友進, 等. 高超聲速氣動熱化學非平衡效應數值分析研究[J]. 現代防御技術, 2012, 40(5): 167-173. ZHANG Sheng-tao, ZHANG Qing-bing, WANG You-jin, et al. Numerical Study of Chemical Non-Equilibrium Effect on Hypersonic Aerodynamic Heating[J]. Modern Denfence Technology, 2012, 40(5): 167-173.

[14] 趙旭, 丁海河, 張慶兵. 高超聲速過渡區氣動熱的DSMC數值模擬研究[J]. 現代防御技術，2010, 38(6): 42-48. ZHAO Xu,DING Hai-he,ZHANG Qing-bing.DSMC Numerical Simulation Research on Aerodynamic Heating at Hypersonic Transition Regime[J].Mordern Defence Technology,2010,38(6):42-48.

[15] 馬毅飛, 趙文平. 窗口輻射對紅外成像探測影響的研究[J]. 系統工程與電子技術, 2005, 27(3):429-430. MA Yi-fei, ZHAO Wen-ping. Effects of Window Radiation on Infrared Imaging Detection[J]. System Engineering and Electronics, 2005, 27(3): 429-430.

[16] 范志剛，肖昊蘇，李輝. 氣動光學頭罩光傳輸數值仿真[J]. 應用光學，2011，32(2)：189-194. FAN Zhi-gang, XIAO Hao-su, LI Hui. Numerical Simulation for Optical Transmission of Aero-Optical Dome [J]. Journal of Applied Optics, 2011, 32(2): 189-194.

Research on Aerodynamic Heating of Conformal Sapphire Dome in Complicated Thermal Environment

XIAO Hao-su, CHEN Cheng, YAO Rui, WEI Yu-fei,ZHANG Li-qin

(Beijing Institute of Electronic System Engineering,Beijing 100854, China)

The external flow field property of the conformal sapphire dome is analyzed. The temperature, stress and strain fields of the dome are obtained. The refractive index field of the dome is computed based on the theory of thermo-optical and elastic-optical effects. It is found that the aerodynamic heating is more serious at stagnation point and in the front of the dome and that the thermo-optical effect has a much greater influence on the refractive index variation of the dome than the elastic-optical effect. This research might play a guiding role in recognizing the mechanism of aerodynamic heating and the heat protection design of the conformal dome.

conformal dome; aerodynamic heating; refractive index field; thermo-optical effect; elastic-optical effect; heat protection design

2014-10-10；

2014-12-23

青年創新基金(LSQC1201-D)

肖昊蘇(1985-)，男，湖南衡陽人。工程師，博士，主要從事光電探測技術研究。

通信地址：100854 北京142信箱30分箱 E-mail：xiaohaosu1985@163.com

10.3969/j.issn.1009-086x.2015.04.003

TJ760

A

1009-086X(2015)-04-0013-06