Thermal radiation characteristics of conformal dome in aero-dynamic environment
-
摘要: 高温整流罩产生强烈的红外辐射,在探测器接收面上形成背景噪声,严重影响成像质量。为了评估气动热环境下高速飞行器共形整流罩热辐射对探测器性能的影响,建立了气动热环境下整流罩热流固耦合计算模型、整流罩热辐射发射以及传输模型,计算了整流罩非均匀温度场、形变场、应力以及应变场,仿真得到了共形整流罩热辐射对高速飞行器光学系统成像质量的影响。研究结果表明:随着时间的不断增长,探测器接收面的最大辐照度也逐渐增大。在第10 s时,0°攻角下共形整流罩热干扰辐照度最小值为0.094 W/m2,最大值为0.108 W/m2。和相同工况下的共形整流罩相比,球形整流罩10 s时的最大辐照度分别是椭球面和抛物面形整流罩最大辐照度的近12倍和7倍,即共形整流罩产生的干扰对探测系统的影响较小,不会对探测器的探测性能产生致命性的影响。Abstract: The high-temperature dome produces intense infrared radiation, which forms the background noise on the detector receiving surface and seriously affects the image quality. To assess the performance degradation as a result of aerodynamic heating of the conformal dome flying at high speed, this paper established the multi-physics coupling model of conformal domes under aerodynamic heating environment, the dome thermal radiation emission and transmission model. Based on these models, the temperature field, the stress field, the strain field and the deformation field of conformal domes were obtained, and the imaging quality of high-speed aircraft affected by heat radiation of conformal dome was analyzed. Results show that as time goes on, the maximum irradiance of the detector is increasing. At 10 s, the minimum value of the thermal interference irradiance of the conformal shroud at 0° attack angle is 0.094 W/m2, and the maximum value is 0.108 W/m2. Compared with conformal domes under the same working condition, the maximum irradiance of the spherical dome at 10 s is about 12 times and 7 times of the maximum irradiance of ellipsoid and parabolic dome, respectively. Conditions thermal radiation interference of the conformal dome has less effect on the detection system comparing with spherical dome, and does not have a fatal effect on the detection performance of the detector.
-
Keywords:
- aero-optics /
- conformal dome /
- heat radiation /
- irradiance
-
引言
光干涉测量方法具有非接触、高效率、高灵敏度等特点,也就必然成为新世纪精密测量技术的重要发展方向[1]。干涉图相位提取是干涉测量中的关键一步,相位提取精度的优劣直接影响最终检测的精度[2]。快速傅里叶变换法(fast fourier transform,FFT)具有处理瞬变波面相位的能力,可有效克服机械振动和气流的影响[3-5],因此成为大口径光学元件测量的常用方法之一。
快速傅里叶变换法(FFT)于20世纪80年代被提出并应用于条纹图的相位提取中[6],对干涉条纹图进行FFT相位提取法的目的是提取出频谱中包含被测波面信息的正一级频谱,也称为滤波过程[7-8],滤波过程的好坏直接影响相位提取精度的高低。目前国内外文献针对FFT相位提取技术主要围绕算法中几个关键步骤进行研究,如文献[9]对干涉图延拓进行分析;文献[10]提出了一种基于FFT时移特性的叠栅条纹细分方法;文献[11]较为综合地分析了影响FFT相位提取法的边缘误差、窗函数、滤波器设计、干涉图延拓及载波条纹数等因素;文献[12]将傅里叶变换相位提取法应用于仿真及实验采集电子散斑干涉(electronic speckle interference,ESPI)信息处理中。为得到更高的相位提取精度,对二维FFT相位提取法的进一步研究仍是国内外研究热点。
针对FFT相位提取技术的滤波过程进行研究,通过选取不同滤波窗对计算机仿真干涉图进行相位提取,对比不同类型的滤波窗滤出的正一级频谱分布形状及其对相位提取精度的影响,最后利用实验采集干涉图对算法的可靠性进行验证。
1 基本原理
1.1 快速傅里叶变换(FFT)基本原理
对于一般的实验干涉仪,引入空间载频
${f_x}$ 、${f_y}$ 后的干涉条纹的强度分布可以表示为$$\begin{split} f(x,y) = & {I_d}(x,y) + c(x,y) \cdot {{\rm e}^{{\rm i}[{\rm{2}}{\text{π}} {f_x}x + {\rm{2}}{\text{π}}{f_y}y]}} + {c^ * }(x,y) \cdot \\ & {{\rm e}^{ - {\rm i}[{\rm{2}}{\text{π}} {f_x}x + {\rm{2}}{\text{π}} {f_y}y]}} \\ \end{split} $$ (1) 其中
$$c(x,y) = {1 / 2} \cdot b(x,y){{\rm e}^{{\rm i}[\varphi (x,y)]}} $$ (2) 式中:
${I_d}(x,y)$ 是干涉条纹的背景光强分布;$b(x,y)$ 是干涉条纹的调制度分布;$\varphi (x,y)$ 为含有待测波面相位信息的相位分布函数;$ * $ 表示复共轭[13]。对(1)式进行傅里叶变换,将空域中的干涉条纹信息转移到频域内,得到的频谱分布函数表示为
$$\begin{split} {\rm{F(}}{f_{\rm{1}}}{\rm{,}}{f_{\rm{2}}}{\rm{)}} = & {\rm{A(}}{f_{\rm{1}}}{\rm{,}}{f_{\rm{2}}}{\rm{)}} + C({f_{\rm{1}}} - {f_x},{f_{\rm{2}}} - {f_y}) + \\ & C*({f_{\rm{1}}} + {f_x},{f_{\rm{2}}} + {f_y}) \end{split} $$ (3) 式中
$A({f_1},{f_2})$ 、$C({f_1} - {f_x},{f_2} - {f_y})$ 、$C*({f_1} + {f_x},{f_2} + {f_y})$ 分别表示零级、正一级、负一级频谱分布函数。选取一个中心频率为$({f_x},{f_y})$ 的滤波器将正一级频谱$C({f_1} - {f_x},{f_2} - {f_y})$ 提取出来并平移至原点,得到$C({f_x},{f_y})$ ,对其进行二维傅里叶逆变换,即可得到$C(x,y)$ 。由(2)式可得所求相位$$\varphi (x,y) = \arctan \frac{{{\rm{lm}}[c(x,y)]}}{{{\rm Re} [c(x,y)]}} $$ (4) 式中
${\rm Re} [c(x,y)]$ 及${\rm{lm}}[c(x,y)]$ 分别为$C(x,y)$ 的实部和虚部。1.2 滤波窗函数
如图1(a)为仿真干涉图的频谱分布,要准确无误地从该频谱图中提取正一级频谱,就要选取适当的滤波窗,但若选取滤波窗不当,就会造成有效频谱信息的缺失[14]。如图1(b)所示为选取滤波窗为圆域形窗口时滤出的正一级频谱,其频谱信息出现明显截断现象,为了减少截断误差,有必要对几种通用的滤波器进行分析与比较。目前二维FFT进行干涉测试数据处理时,使用较多的滤波器有海明(Hamming)窗、汉宁(Hanning)窗、高斯(Gaussian)窗、布莱克曼(Blackman)窗。
根据干涉图频谱分布的这一特点,选取与其分布形状相似的几种滤波窗函数,表1为几种滤波窗函数的一维表达式。
表 1 几种滤波窗函数Table 1. Several filter window functionsWindow type Function Parameter Blackman $\omega (n) = \left\{ \begin{array}{l} 0.42 - 0.5\cos (\frac{{2{\text{π}} n}}{{N - 1}}) + \\ 0.08\cos (\frac{{4{\text{π}} n}}{{N - 1}})\quad \quad \quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $ main lobe width ${\rm{12}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 56 dB Gaussian $\begin{array}{l} \omega (n) = \;\exp \left[ { - \frac{1}{2}{{\left( {\frac{{n - N/2}}{{\sigma N/2}}} \right)}^2}} \right]\quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \sigma \leqslant 0.5 \\ \end{array} $ main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 44 dB Hanning $\omega (n) = \left\{ \begin{array}{l} 0.5\left[ {{\rm{1}} - \cos (\frac{{2{\text{π}} n}}{{N - 1}})} \right]\quad \quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $ main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 31 dB Hamming $\omega (n) = \left\{ \begin{array}{l} 0.5{\rm{4 - 0}}{\rm{.46}}\cos (\frac{{2{\text{π}} n}}{{N - 1}})\quad 0 \leqslant n \leqslant N - 1 \\ \quad \quad \quad \quad 0\quad \quad \quad \quad \quad {\rm otherwise} \\ \end{array} \right.\quad $ main lobe width ${\rm{8}}{\text{π}} {\rm{/}}N$ Side lobe peak attenuation 41 dB 根据表1中各类窗函数的分布函数绘制时频振幅响应分布,选取N=30作为采样点,为直观起见,绘制其一维空/频分布如图2,其对应二维空/频分布,可由一维空/频分布经旋转得到。
根据图2多种滤波窗函数的时域/频域响应及表1中参数可知,不同窗函数主要差别在于主瓣宽度及旁瓣衰减速率。其中主瓣宽度主要影响信号的能量分布,旁瓣峰值衰减速率影响能量的泄露程度[15],旁瓣越高,能量泄露越严重,衰减速率就越慢[16]。对于FFT滤波过程,要保证滤出正一级频谱的同时滤掉其他频谱成分,且要最大程度地保留正一级频谱沿x轴和y轴延伸区域内的频谱成分,需综合考虑滤波窗的频谱分布形状及对应参数。
2 仿真分析
2.1 滤波处理
针对第一节提出的频谱截断现象,选取Hamming、Hanning、Gaussian、Blackman四种与正一级频谱分布形状相近的滤波窗对仿真干涉图的频谱图(图1(a))进行滤波处理,其中仿真干涉图尺寸为300×300像素,零频频谱中心频率位于(151,151),正一级频谱中心频率位于(164,164),即滤波窗中心频率确定。根据距离公式求得二者像素间距为18.384 7像素,则滤出正一级频谱对应窗口函数的滤波半径可以此为参考。如图3分别为其对应频谱响应(滤波半径取值为18)及选取该窗口滤出的正一级频谱分布图。
分析图3(a)~(d)可得:经Hamming窗和Hanning窗提取的正一级频谱较完整地保留了原始频谱成分,没有其他频谱成分的混入。分析其造成的原因主要是由于主瓣宽度覆盖了正一级频谱带宽,且第一旁瓣高度较低,未出现频谱泄露现象,即未混入零频及负一级频谱;而经Gaussian、Blackman窗提取的正一级频谱分布虽较完整地保留了原始频谱成分,但有其他频谱成分的混入。分析对应窗口函数及频谱响应分布得到,其主瓣宽度覆盖正一级频谱宽度,但由于其第一旁瓣有一定幅值,造成了在此窗函数x及y方向上引入较多频谱泄露,即混入的零频及正一级频谱。
2.2 相位误差分析
将不同滤波窗函数提取出的正一级频谱进行移中并做逆傅里叶变换,利用反正切函数求得包裹相位,再利用离散余弦(discrete-cosin-transform,DCT)解包裹法提取连续相位,最后进行波面拟合,得到波面分布。计算不同滤波窗口选取下的FFT相位提取法得到波面的波面峰谷值(peak-to-valley,PV)和波面均方根值(root-mean-square,RMS),分析相位提取精度并与原始波面进行残差计算。其中原始波面分布评价参数PVa为0.246 4λ,RMSa为0.056 4λ,经FFT法得到的波面分布评价参数为PV1、RMS1,记ΔPV=PV1−PVa,ΔRMS=RMS1−RMSa。其值如表2所示。
表 2 不同滤波窗下FFT算法结果比较 (λ=632.8 nm)Table 2. Comparison of FFT algorithm results in different filtering windows (λ=632.8 nm)Window Type PV1/λ ∆PV/λ RMS1/λ ∆RMS/λ Hamming 0.257 9 0.008 5 0.056 3 −0.000 1 Hanning 0.258 6 0.012 2 0.056 3 −0.000 1 Gaussian 0.262 0 0.015 6 0.056 3 −0.000 1 Blackman 0.248 3 0.001 9 0.056 2 −0.000 2 对比上述ΔPV及ΔRMS可知,选取Blackman窗及Hamming窗处理得到波面的评价参数在数值上较其他窗口函数更逼近原始波面的评价参数,结果优于Gaussian窗及Hanning窗。其中最优窗口类型为Hamming窗,其波面峰谷值残差为0.008 5λ,波面均方根值残差为0.000 1λ。
3 实验验证
3.1 实验装置
为验证仿真结果的可行性和精度,分别采用二维傅里叶变换法(2D-FFT)和移相干涉法(PSI)对实验采集的干涉图进行处理并比较其测试结果。实验阶段所用的菲索干涉仪是美国Zygo公司生产的Verifire PE激光干涉仪,该干涉仪是共光路面形计量干涉仪,光源为低功率的632.8 nm的氦氖激光,采用PZT移相干涉原理可以对多种光学元件的面形进行检测,其干涉仪装置如图4所示。
3.2 算法处理
使用Zygo干涉仪采集实验干涉图如图5(a)所示,其尺寸为300×300像素。经消除白噪声后的条纹图如图5(b),对处理后的条纹图进行二维傅里叶变换,得到频谱分布如图5(c)。选取Hamming窗将正一级频谱取出并滤掉其他频谱,并将正一级频谱平移到整个频谱中心,如图5(d)~(e);对移中后的正一级频谱做二维傅里叶逆变换,利用反正切函数求取包裹相位如图5(f),采用离散余弦解包裹法进行处理,其结果如图5(g)。
选取其他3种滤波窗(Hanning、Blackman、Gaussian)对实验采集条纹图进行上述相位提取的算法处理,图6(a)~(d)所示为上述4种滤波窗提取出的实验干涉图正一级频谱分布。
对图6进行分析可得,实验干涉图在不同滤波窗处理下的正一级频谱分布与仿真干涉图的结果呈现相同的趋势。即Hamming窗和Hanning窗较完整地保留正一级频谱且无其他频谱成分混入,而Blackman窗及Gaussian窗仍有其他频谱成分混入。
3.3 结果分析
表3为用FFT相位提取法(不同滤波窗口)得到的波面及PSI法测得波面的评价参数汇总。
表 3 FFT(不同滤波窗)&PSI算法的结果比较 (λ=632.8 nm)Table 3. Comparsion of results in FFT(different filtering windows) and PSI (λ=632.8 nm)Method PV/λ RMS/λ PSI 0.445 4 0.002 7 FFT (Hamming Window) 0.436 1 0.002 2 FFT (Hanning Window) 0.417 7 0.002 8 FFT (Gaussian Window) 0.326 2 0.002 2 FFT (Blackman Window) 0.527 3 0.002 9 由表3可知,选取Hanning窗及Hamming窗处理得到波面的评价参数在数值上较其他窗口函数更逼近PSI法所得波面评价参数,结果优于Blackman窗及Gaussian窗。其中最优窗口类型为Hamming窗,其波面峰谷值残差为0.009 3λ,波面均方根值残差为0.000 5λ。
4 结论
针对二维傅里叶变换相位提取算法滤波过程中滤波窗的选取进行研究。首先分析干涉条纹图的频谱分布及各类滤波窗口函数的参数及其频谱响应;其次进行仿真实验,针对干涉图的频谱分布特点选取合适的窗口进行滤波处理并进行相位提取,分析不同窗口选取下对相位提取精度的影响;最后利用立式菲索干涉仪进行实验验证并与移相干涉测量法所测的结果进行对比。实验及仿真结果表明:选取Hamming窗提取的正一级频谱较完整地保留了原始频谱成分,相位提取精度优于0.01λ,可进一步应用于大口径光学元件的测量中。
-
表 1 300 K蓝宝石晶体的物理特性
Table 1 Physical properties of sapphire crystal at 300 K
密度
/(kg/m3)熔点/K 膨胀系数
/(10-6K-1)比热容
/(J/(kg·K))弹性模量
/GPa导热系数
/(W/(m·K))屈服强度
/(MPa)泊松比 3 980 2 050 5.3 750 344 36 300 0.27 -
[1] 殷兴良.现代光学新分支学科——气动光学[J].中国工程科学, 2006, 7(12):1-6. doi: 10.3969/j.issn.1009-1742.2006.12.001 Yin Xinliang.A new subdiscipline of contem porary optics-aero-optics[J].Enfineering Science, 2006, 7(12):1-6. doi: 10.3969/j.issn.1009-1742.2006.12.001
[2] 马毅飞, 赵文平.窗口辐射对红外成像探测影响的研究[J].系统工程与电子技术, 2005, 27(3):427-430. doi: 10.3321/j.issn:1001-506X.2005.03.013 Ma Yifei, Zhao Wenping. Effects of window radiation on infrared imaging detection[J].Systems Engineeting and Electronics, 2005, 27(3):427-430. doi: 10.3321/j.issn:1001-506X.2005.03.013
[3] Burns S P, Howell J R, Klein D E. Assessment of the swartz-wendroff finite element formulation for combined mode heat transfer, radiative heat transfer current research[J]. ASME HTD, 2002, 276:1-7.
[4] Brandon S, Derby J J. A finite element method for internal radiative heat transfer and its application to analysis of the growth of semitransparent crystals[J].ASME HTD, Fundamentals of Radiation Heat Transfer, 1992, 106:1-6.
[5] 范志刚, 于翠萍, 肖昊苏, 等.气动光学头罩热辐射效应数值仿真研究[J].应用光学, 2012, 33(4):649-653. http://d.old.wanfangdata.com.cn/Periodical/yygx201204001 Fan Zhigamg, Yu Cuiping, Xiao Haosu, et al.Numerical simulation of termal radiation effect about optical dome in aerodynamic thermal environment[J].Journal of Applied Optics, 2012, 33(4):649-653. http://d.old.wanfangdata.com.cn/Periodical/yygx201204001
[6] 刘珂, 陈宝国, 李丽娟.空空导弹红外导引头技术发展趋势及关键技术[J].激光与红外, 2011, 41(10):1117-1121. doi: 10.3969/j.issn.1001-5078.2011.10.012 Liu Ke, Chen Baoguo, Li Lijuan.Development tendency and key technology of IR seeker for air-to-air missile[J].Laser & Infrared, 2011, 41(10):1117-1121. doi: 10.3969/j.issn.1001-5078.2011.10.012
[7] Frase B S, Hemingway A. High performance faceted domes for tactical and strategic missiles[J]. SPIE, 1994, 2286:486-492. http://cn.bing.com/academic/profile?id=ed22b60a788f39ad134eb870b4c73d05&encoded=0&v=paper_preview&mkt=zh-cn
[8] Knapp D J, Mills J P, Hegg R G, et al. Conformal optics risk reduction demonstration[J]. SPIE, 2001, 4375:147. http://cn.bing.com/academic/profile?id=7de2efb40080c09dba9797a53e90a2b6&encoded=0&v=paper_preview&mkt=zh-cn
[9] Shorey A B, Kordonski W, Tracy J, et al. Developments in the finishing of domes and conformal optics[J]. SPIE, 2007, 6545:65450Q1-65450Q9. http://cn.bing.com/academic/profile?id=abbf2cdb9a1e169e8d03e44586dc4971&encoded=0&v=paper_preview&mkt=zh-cn
[10] Shannon R R. Overview of conformal optics[J]. SPIE, 1999, 3705:180-188. doi: 10.1117-12.354622/
[11] 王勖成.有限单元法[M].北京:清华大学出版社, 2003:441-464. Wang Maocheng.Finite element method[M].Beijing:Tsinghua University Press, 2003:441-464.
[12] Wakaki M, Kudo K, Shibuya T.Physical properties and data of optical materials[C]. Boca Raton: CRC Press, 2007: 383-535.
[13] Visknata R, Anderson E E. Heat transfer semitransparent solids[J].Advances in Heat Transfer, 1975, 11:317-441. doi: 10.1016/S0065-2717(08)70077-7
[14] 张琳.梯度折射率介质内辐射传递方程数值模拟的有限元法[D].哈尔滨: 哈尔滨工业大学学位论文, 2009. Zhang Lin.Finite element method for solving radiative transfer equation in graned index mediium[D].Haerbin: Harbin Institute of Technology, 2009: 32-63.
-
期刊类型引用(7)
1. 张瑜,张峰,郭芮,苏瑛,张云龙,许增奇,王阜超. 磁流变抛光二次非球面工件位置对刀校正. 应用光学. 2022(03): 518-523 . 本站查看
2. 刘春阳,李洛克. 基于轮廓测量的集料表面粗糙度表征方法试验研究. 工业安全与环保. 2021(03): 15-19+25 . 百度学术
3. 夏良静. 工艺参数对光学元件加工精度影响的数学模型研究. 激光杂志. 2020(08): 201-204 . 百度学术
4. 张晋烨,张海云,赵玉刚,张军,王显达. 磨料振动对曲轴抛光的影响. 组合机床与自动化加工技术. 2019(01): 46-48 . 百度学术
5. 高丰,李琦,向阳,宁旭. 高清折反式视频显微镜. 光学技术. 2019(02): 136-141 . 百度学术
6. 刘红丽. 散粒磨料振动光饰技术专利分析. 科技风. 2019(14): 141 . 百度学术
7. 马旭德,张东梅,尚春民. 固着磨料球面高速研磨加工技术研究. 新型工业化. 2018(04): 34-42 . 百度学术
其他类型引用(4)