Numerical Simulation of Radiation Characteristics of Aircraft Exhaust Systems with Different Nozzles
-
摘要: 排气系统是飞行器最主要的红外辐射源,其喷管的形状类型对排气系统红外辐射强度的大小及分布有很大影响。本文建立了3种不同类型喷管的三维模型,在此基础上运用ANSYS软件模拟了各自排气系统的温度场分布,结合Curtis-Godson(C-G)谱带法对各类型喷管红外辐射特性进行了计算与对比研究。结果表明:在出口面积相同的条件下,二元矩形S弯喷管的尾焰核心区域面积最小,约为轴对称圆形喷管的60%;在矩形喷口的宽边探测面上,二元矩形S弯喷管的红外辐射强度最小。3类喷管中,二元矩形S弯喷管隐身性能最好,二元矩形喷管次之,轴对称圆形喷管最差。Abstract: The exhaust system is the most important infrared radiation source of an aircraft, and the shape of the nozzle contributes to the infrared radiation characteristics of the exhaust system. Three types of 3D nozzles were built, and the temperature field of the plume was simulated using ANSYS14.5. Then, the spectral infrared radiation characteristics of the plume were obtained using the single band Curtis-Godson (C-G) approximation method. The results show that under the same exit area, the core area of the S-shaped nozzle plume is minimum and is approximately 60% of the axisymmetric circular nozzle plume; in the rectangular nozzle wide edge detection surface, the infrared radiation of the dual rectangular nozzle is minimum, and among the three types of nozzles, the stealthy performance of the two-element rectangular S curved nozzle is the best, the two-element rectangular nozzle takes the second place, and the axial symmetrical circular nozzle is the worst.
-
Keywords:
- aircraft /
- nozzle /
- plume /
- temperature field /
- infrared radiation
-
0. 引言
编织复合材料是由碳纤维材料为增强体,以环氧树脂为基质,由不同方向纤维铺层叠加组成的新型复合材料,其具有高比强度、高比刚度、耐腐蚀、耐老化、耐热性的优点。编织复合材料由于制造成本低、生产效率高、机械性能突出、具有纤维互锁的网状结构[1]。编织复合材料在制造和适用过程中会碳纤维复合材料实现轻量化的同时,还能提供产品相当的安全性和可靠性,已被广泛应用于航空航天、汽车、轮船、运动器械等[2-4]。NASA喷气发动机风扇机匣、波音787飞机风扇机匣都采用的是编织复合材料。复合材料结构是由纤维铺层叠加而成,在面临碰撞和冲击时,纤维和纤维之间都是由基体传递载荷,冲击容易造成材料基体开裂、分层、脱粘等缺陷,编织层压板冲击损伤模式主要是层间分层,不会出现纤维劈裂等损伤模式[5]。在结构表面看不出来损伤缺陷,目视检测也无法检测出缺陷信息,但其内部可能因瞬时冲击能量的吸收,产生了严重的损伤,使得材料结构失效、材料性能降低,对结构承载能力造成大幅度下降[6-7]。因此,有效及时检测编织复合材料内部损伤或缺陷情况,对于保证复合材料结构件可靠性具有重要的意义。
超声红外热成像技术以超声作为热激励源,激励过程中仅界面类缺陷或闭合裂纹发热,正常区域不发热,属于选择性激励方式,更适合于不规则形状材料界面缺陷或闭合裂纹检测。因此目前超声红外热成像技术主要应用在对金属裂纹的检测[8-11]。韩梦[12]等采用超声红外热成像对铝金属的疲劳微裂纹进行检测,微裂纹最小宽度为20 μm;习小文[13]等采用超声红外热成像技术对航空发动机叶片裂纹进行检测,最小检测裂纹宽度为0.5 μm。G. Busse等[14]最早将超声红外锁相热成像检测技术引入到复合材料界面类缺陷检测。田干等[15]采用超声红外热成像技术对航空复合材料进行数值仿真研究,建立含裂纹缺陷复合材料的有限元模型;金国锋、张炜等[16-17]通过数值计算和试验研究了超声红外热成像技术对复合材料冲击损伤检测的适用性;李胤等[18]采用超声红外热像检测对复合材料层合板24 J和29 J冲击损伤进行检测研究,并与超声C扫进行对比,结果表明超声红外热像检测可检测出复合材料不同损伤形式,并且比超声C扫检测速度快、检测精度高。
目前对冲击能量损伤研究大部分基于1个或2个冲击能量对复合材料损伤的研究,为了解不同冲击能量对编织复合材料损伤及扩展情况,采用10 J、20 J、30 J、40 J、50 J冲击能量对不同编织复合材料进行冲击,并对冲击过程中的冲击能量-时间曲线进行了记录,采用超声红外热成像技术对不同冲击能量损伤试件进行检测与分析。
1. 基本原理
1.1 超声热成像技术的检测原理
超声红外热成像技术是将超声波激励和红外热成像两者结合的新型无损检测方法,属于主动式红外无损检测方法的一种。超声红外热成像技术的检测原理如图 1所示,是通过将高频超声波(20~40 kHz)耦合注入待测试件,超声波经过界面耦合在试件中传播,在超声波的作用下,试件缺陷由于塑性变形、摩擦作用、粘弹性效应等产生热量,使得裂纹、损伤区域温度升高,在热传导作用下,其裂纹或损伤位置对应的试件表面温度也相应的升高,利用红外热像仪采集试件表面的红外热序列图像,对红外序列图像进行处理便可得到其内部缺陷信息。相较于其他激励方式的红外热成像技术,例如:闪光灯激励、红外灯激励、激光激励等,这些激励方式都属于面激励,不论是缺陷区域还是无缺陷区域,试件表面被整体加热,表面吸收的能量通过热传导到试件内部,遇到隔热性缺陷时,热波反射回表面。超声红外热成像技术是对试件局部加热,其仅对缺陷的区域进行加热,无缺陷区域不进行加热,属于选择性加热方式,在超声激励下产生相对位移的缺陷区域温度都将升高,相当于是体加热方式,并且从缺陷区域通过热传导的方式传播到试件表面,仅传播了单次的距离。
1.2 超声热成像技术的检测原理
超声波在试件内界面缺陷处衰减产生热量,缺陷处则成为热源。由此可知,超声波激励产生的热流率q的传播满足热传导方程:
$$ {k_x}\frac{{{\partial ^2}T}}{{\partial {x^2}}} + {k_y}\frac{{{\partial ^2}T}}{{\partial {y^2}}} + {k_z}\frac{{{\partial ^2}T}}{{\partial {z^2}}}{\text{ + }}q{\text{ = }}\rho c\frac{{\partial T}}{{\partial t}} $$ (1) 式(1)中:T为温度,℃;kx、ky、kz分别为在三个维度上的导热率,W/(m·℃);ρ为密度,kg/m3,c为比热容,J/(kg·℃)。
热传导方程(1)的初始条件为:
$$ T(x, y, z, t)=T_{0}\qquad t=0$$ (2) 式(2)中:T0为环境温度,℃。
热传导方程支配热传导及其扩散过程,试件检测过程时间比较短,表面温度变化比较小,忽视试件表面对流和热辐射,则:
$$ \left\{ \begin{array}{l} \frac{{\partial T}}{{\partial t}} - \alpha {\nabla ^2} = q(t) \hfill \\ T(0) = {T_0} \hfill \\ \end{array} \right. $$ (3) 式(3)中:α为热扩散率,m2/s;∇2为拉普拉斯算子,
$$ {{\mathit{\nabla}} ^2}{\text{ = }}\frac{{{\partial ^2}}}{{\partial {x^2}}} + \frac{{{\partial ^2}}}{{\partial {y^2}}} + \frac{{{\partial ^2}}}{{\partial {z^2}}} $$ 分别对上式进行关于空间(x, y, z)的傅里叶变换得到:
$$ \left\{ \begin{array}{l}\frac{{\rm{d}}{T}^{f}}{{\rm{d}}t}-\alpha {ϵ}^{2}{T}^{f}={q}^{f}(ϵ, t)\\ {T}^{f}(ϵ, 0)={T}_{0}{}^{f}\end{array} \right. $$ (4) 式(4)中:qf为对生热率q作温度T的傅里叶变换;$ ϵ $为傅里叶变换过程中的中间变量,没有特殊的物理含义;Tf为对温度T的傅里叶变换。
常微分方程的解析解为:
$$ {T}^{f}(ϵ, t)={T}_{0}{}^{f}{{\rm{e}}}^{-\alpha {ϵ}^{2}t}+{\displaystyle {\int }_{0}^{t}{q}^{f}}(ϵ, t){{\rm{e}}}^{-\alpha {ϵ}^{2}(t-\tau )}{\rm{d}}\tau $$ (5) 对Tf(∈, t)进行傅里叶反变换得到:
$$ \begin{array}{l}T\left(X, t\right)={T}_{0}{{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K(X-ϵ, t){\rm{d}}ϵ+\\ {{\displaystyle \int }}_{0}^{t}{\rm{d}}\tau {{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K\left(X-ϵ, t-\tau \right){q}^{f}(ϵ, \tau ){\rm{d}}ϵ\end{array} $$ (6) 式(6)为超声波在试件内界面缺陷热源公式,式子表示,热源公式是一个复杂函数表达式。
2. 试验与结果分析
2.1 碳纤维复合材料冲击损伤
试验用的碳纤维复合材料是T700编织复合材料层压板,编织技术制造网格状堆叠而成,铺层方式为-45/90/45/0,经过浸渍固化直接形成复合材料结构,试件共5件,试件尺寸均为150 mm×100 mm×4 mm,对5件试件进行冲击试验,冲击能量分别是10 J、20 J、30 J、40 J、50 J,冲击能量-时间曲线如图 2所示,冲击后试件光学图像如图 3所示,图 3(a)为试件正面,为冲击面,图 3(b)为试件背面,光学图像上很难看到试件表面有损伤。
2.2 超声热成像检测系统
超声红外热成像系统如图 4所示,系统主要包括超声激励电源、超声枪、红外热像仪、夹具、图像处理和控制系统。超声枪包括换能器、调幅器、工具杆,最大功率为1200 W,激励频率为20 kHz,超声激励时间可设置为0.2~20 s。通过超声激励电源箱控制器调节输出功率大小,预紧力大小通过机械控制,预紧力大小最大可设置1000 N。红外热像仪型号为FLIR A655sc,属于非制冷型红外热像仪,工作波段为7.5~14 μm,室温下等效噪声温差为0.04 K,图像分辨率为640 pixel×480 pixel。
试验过程中超声枪不直接作用于试件表面,试件固定在夹具上,夹具固定在平台上,超声枪顶在夹具上面,在试验过程中,需要确保超声波有效注入到夹具中,夹具和试件要接触良好,超声枪在20 kHz的激励下,振幅约为几十微米,如果超声枪与夹具之间接触不好,则超声能量无法传递到夹具上,因此设置超声枪与夹具的预紧力为200 N,预紧力是通过旋转预紧力手柄迫使预紧力系统单元向前移动,内部弹簧受到挤压,从而对超声枪头产生压力,超声枪头压紧夹具。校正好红外热像仪,超声激励时间为1.5 s,采集频率为50 Hz,采集时间为10 s。
2.3 检测结果分析
试验中共有5种不同冲击能量的试件,选取30 J能量冲击能量的试件3#进行分析,图 5所示为超声红外热成像系统检测试件3#所示的红外序列热图,从图 5中可以看出,超声激励前,看不到任何热斑信息,随着超声持续激励,图像中热斑信息逐渐清晰,0.08 s时稍微可以看到零星散点,到0.5 s时逐渐形成了较为清晰的热斑,超声激励时间为1.5 s,红外序列中时间为1.6 s时热斑最为清晰,在超声激励结束后,冲击损伤区域不能持续发热,仅能靠自身热传导把热斑区域的能量传递到低温区域,逐渐达到温度平衡,从红外序列图中表现为图像越来越模糊,对比度越来越低,到8.66 s时已不能分辨出热斑轮廓信息,随着热传导的进行,将不能检测到任何热斑信息,所检测到的红外图像与0.02 s时图像一致。
图 6为不同位置处A点、B点、C点、D点、E点的温升变化曲线,由图 6曲线可以看出,超声激励后,损伤区域温度迅速上升,且各点温度上升不尽相同,A点为试件3#的温升最高点,最高温升可达到8℃,B点和C点为其他区域温升点,这3个点温升曲线有一个特点,温升曲线都是呈现先上升再下降的过程,这3个点代表着试件内部有持续发热的热源,这些热源就是在超声激励下,损伤处由于摩擦生热等原因产生的热源。D点温升曲线只有温升的过程,没有下降的过程,从图像上看,C点位热斑区中间位置,D点处于C点上方,D点处于热斑区域的边缘,边缘位置损伤比较轻微,温升没有那么高,但依然会在持续激励过程中温度不断上升,等热激励结束以后,随着热扩散的进行温度将下降,由于其处于热斑的边缘区域,热斑处的热量会由于热扩散向低温区域传播,C点和D点的温度由于热扩散会逐渐接近,因此导致D点温度下降曲线基本是平的。E点代表着无热斑区域,即无损伤区域,温升曲线基本是平坦的,E点对应位置的内部没有发热源,与上文分析一致。
为了更形象描述不同帧之间空间温度变化曲线,分别取图 5所示前6个时间点的对应图像的空间温度变化曲线,直线所对应的位置为图 7所示红色空间直线对应的位置,左边是起始点,对应的距离为0,右边对应的距离为141.5 mm。曲线横坐标为距离,纵坐标为温度,环境温度约为22℃,t=0.02 s对应的空间曲线是超声激励之前的空间温度曲线,温度基本相同,没有太明显变化,t=0.08 s红外图像上可以看到零星散点,红色空间直线正好位于散点上,在空间温度曲线上正好可以看到损伤处对应位置处有凸起,随着超声激励的加载,损伤区域空间温度曲线温度逐渐升高,并且其对应的宽度稍微有所增加,但温度升高的点仅限于损伤区域,非损伤区域对应的曲线基本是不变的,重叠在一起;损伤区域t=0.24 s对应的空间温度曲线就比t=0.08 s对应的空间温度曲线要高,同样的损伤区域t=0.50 s对应的空间温度曲线比t=0.24 s对应的空间温度曲线高,到超声激励结束时间附近t=1.6 s,空间温度曲线升高最高值,图 7对应的最高曲线。超声激励结束后,损伤区域不再主动发热,内部热源消失,温度开始下降,空间温度曲线也随之下降,t=2.28 s损伤区域对应的曲线的值就比t=1.6 s对应曲线值有所下降。图像中红色箭头指示的空间距离约60 mm时间为t=0.50 s和t=1.60 s损伤区域的边缘点,温度比其旁边两个尖峰都要低,当超声激励结束以后,随着热扩散的进行,t=2.28 s对应的温度反而要比t=1.60 s对应的温度要高,此时该区域对应的2个尖峰也变成了一个大而平滑的尖峰曲线,在红外序列图像中就表现为图像变模糊,与图 5所示图像相对应。
2.4 不同冲击能量试件结果与分析
图 8为不同冲击能量冲击试件的检测图像,从光学图像上都看不到表面的损伤信息;从检测图上看,冲击能量为10 J的1#试件检测图上有3个小的热斑区域,冲击能量为20 J的2#试件看出热斑区域连成了一个整体,冲击能量为30 J的3#试件可以看到热斑已延展开,冲击能量为40 J的4#试件检测图上看出形成的热斑延展成一条长的缝,冲击能量为50 J的5#试件检测热图上热斑几乎贯穿对角线,并且还有额外2条分支。试件是二维编织技术制作的二维编织复合材料,编织使得纤维之间连成了有机的整体,使得二维编织复合材料展现出了较强的整体性。在冲击能量的作用下,外表面几乎看不出损伤信息,当编织复合材料受到小的外部冲击后,不足以破坏有机的整体,损伤区域仅限于冲击点周围;当编织复合材料受到大的外部冲击后,其内部将会产生大的面积分层,损伤展现出延展性。
2.5 损伤区域识别
衡量超声红外热成像技术对编织复合材料检测能力关键是对损伤区域进行提取,采用双阈值分割,对图 5中时间为t=0.50 s、t=1.60 s、t=2.28 s图像进行分割,分割效果如图 9所示,无法达到理想的分割效果。底部区域是超声激励中夹具位置,夹具在超声激励作用下温度升高,通过热传导会把试件底部的温度升高,采用双阈值分割方法无法把该区域给剔除,对比3张图分割效果可以看出,选用不同时间的图像,分割的结果不尽相同。
由于红外序列图中单帧图像由于热传导一直在发生变化,采用单帧图像分割算法会由于选取的帧不同而导致分割区域的差异,为了准确提取损伤区域,本文采用曲线分类算法,利用每一个像素点的温升曲线进行分类,定义损伤区域强弱程度,温升曲线有上升沿和下降沿的点为热斑中心区域,该区域损伤比较严重,图 6所示对应的A、B、C三个点,该像素点置红色;温升曲线有上升沿但无下降沿的点为损伤比较轻微,图 6所示对应的D点,该像素点置绿色;温升曲线基本不变化的为无损伤区域,图 6所示的E点,该像素点置黑色。不同冲击能量试件分割效果如图 10所示。对比于双阈值分割,采用曲线分类算法分割效果比较好,且图像中有两种颜色表示,红色区域基本对应于检测图像中的亮热斑区域。
表 1为不同冲击能量试件损伤面积统计,从表格中可以看出,冲击能量越大,损伤面积越大,但是试件4#损伤面积反而比试件3#损伤面积小。注意到在冲击试验中,图 2所示冲击能量-时间曲线,试件4#峰值冲击能量是40 J,正常冲击结束处的能量是30 J左右,但冲击结束处的能量只有20 J,反而低于试件3#冲击结束处的能量24 J,所以造成这种对应关系不匹配的情况。
表 1 不同冲击能量试件损伤面积Table 1. Damage area of specimens with different impact energySpecimen Red area/mm2 Total area/mm2 1#-10J 630 926 2#-20J 1475 2525 3#-30J 2909 3929 4#-40J 2896 3383 5#-50J 4517 6291 由于试件4#冲击结束处能量偏低,因此把其他4个试件的冲击能量和损伤面积作为损伤面积关系图,如图 11所示,可以发现损伤面积和冲击能量成线性关系,与碳纤维层合板冲击损伤面积和冲击能量关系曲线相同[19]。
3. 结论
编织复合材料冲击损伤主要涉及复合材料的分层损伤,这些损伤无法采用目视检测,且损伤会使得结构承载能力严重降低,从而威胁复合材料构件的安全使用。制作了编织复合材料层压板试件共5件,并对不同冲击能量下冲击损伤开展了超声红外热成像技术研究。超声激励试件时,冲击损伤区域由于摩擦等因素产生热量,相当于内部有热源的存在;无损伤区域不生热,使用红外热像仪采集红外图像序列便可实现对损伤区域检测,本文以试件3#为例对超声激励过程的温升曲线、空间温度曲线进行了详细的分析。并对10~50 J冲击能量试件进行了对比分析,最后采用曲线分类算法对损伤区域进行了定量识别,发现编织复合材料损伤面积和冲击能量成线性关系。后续的工作是采用其他技术手段(如超声C扫描无损检测技术)确认试件损伤区域面积,以及对损伤区域的评估,将对损伤定量化及曲线分类算法得到更进一步提高。
-
表 1 喷管参数
Table 1 Parameters of three kind of nozzles
Nozzle type Nozzle length/m Nozzle height/m Nozzle width/m Nozzle inlet diameter/m Axisymmetric 1 0.225 0.225 0.6 Rectangle 1 0.2 0.8 0.6 S-pipeline 1 0.2 0.8 0.6 -
[1] Aguiar P, Brett D, Brandon N P. Solid oxide fuel cell/gas turbine hybrid system analysis for high-altitude long-endurance unmanned aerial vehicles[J]. International Journal of Hydrogen Energy, 2008, 33(23): 7214-7223. DOI: 10.1016/j.ijhydene.2008.09.012
[2] Suresh Patel, Christopher Blac. Statistical modeling of F/A-22 flight test Buffet data for probabilistic analysis[C]// Aiaa/asme/asce/ahs/asc Structures, Structural Dynamics and Materials Conference, 2005: DOI: 10.2514/6.2005-2289.
[3] Roblin A, Baudoux P E, Chervet P. UV missile-plume signature model[C]//Proc of SPIE, 2002, 4718: 344-355.
[4] Bakker E J, Fair M L, Schleijpen H M A. Modeling multispectral imagery data with NIRATAM v3.1 and NPLUME v1.6[C]// Proceedings of SPIE, 1999, 3699: 80-91.
[5] Watkins W R. Thermal background modeling and its use in science and technology[C]//Proceedings of SPIE, 1992, 1687: 84-106.
[6] 张海兴, 岳敏. 飞机红外辐射的理论计算[J]. 西安电子科技大学学报, 1997(1): 78-81. https://www.cnki.com.cn/Article/CJFDTOTAL-XDKD701.012.htm ZHANG Haixing, YUE Min. Theoretical calculation of the IR radiation of an aeroplane[J]. Journal of XIDIAN University, 1997(1): 78-81. https://www.cnki.com.cn/Article/CJFDTOTAL-XDKD701.012.htm
[7] 未军光, 杨青真, 李岳锋. 飞机发动机排气系统红外辐射强度数值仿真[J]. 计算机仿真, 2011, 28(4): 66-70. DOI: 10.3969/j.issn.1006-9348.2011.04.016 WEI Junguang, YANG Qingzhen, LI Yuefeng. Numerical simulation of infrared radiation for a aroengine exhaust system[J]. Computer Simulation, 2011, 28(4): 66-70. DOI: 10.3969/j.issn.1006-9348.2011.04.016
[8] 冯云松, 李晓霞, 路远, 等. 矩形喷管外尾焰红外辐射特性的数值计算[J]. 兵工学报, 2013, 34(4): 437-442. https://www.cnki.com.cn/Article/CJFDTOTAL-BIGO201304010.htm FENG Yunsong, LI Xiaoxia, LU Yuan, et al. Numerical calculation of infrared radiation characteristics of the exhaust plume outside a rectangular nozzle[J]. Acta Armamentalrii, 2013, 34(4): 437-442. https://www.cnki.com.cn/Article/CJFDTOTAL-BIGO201304010.htm
[9] 冯云松, 杨华, 杨莉, 等. 宽高比对矩形喷管外尾焰红外辐射的影响[J]. 光电工程, 2012, 39(11): 49-54. DOI: 10.3969/j.issn.1003-501X.2012.11.008 FENG Yunsong, YANG Hua, YANG Li, et al. Influence of aspect ratio on infrared radiation of the exhaust plume outside a rectangular nozzle[J]. Opto-Electronic Engineering, 2012, 39(11): 49-54. DOI: 10.3969/j.issn.1003-501X.2012.11.008
[10] 高翔, 杨青真, 施永强, 等. 出口形式对双S弯排气系统红外特性影响研究[J]. 红外与激光工程, 2015, 44(6): 1726-1732. DOI: 10.3969/j.issn.1007-2276.2015.06.009 GAO Xiang, YANG Zhenqing, SHI Yongqiang, et al. Numerical simulation of radiation intensity of double S-shaped exhaust system with different outlet shapes[J]. Infrared and Laser Engineering, 2015, 44(6): 1726-1732. DOI: 10.3969/j.issn.1007-2276.2015.06.009
[11] 刘常春, 吉洪湖, 黄伟, 等. 一种双S弯二元喷管的红外辐射特性数值研究[J]. 航空动力学报, 2013, 28(7): 1482-1488. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201307007.htm LIU Changchun, JI Honghu, HUANG Wei, et al. Numerical simulation on infrared radiation characteristics of serpentine 2-D nozzle[J]. Journal of Aerospace Power, 2013, 28(7): 1482-1488. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201307007.htm
[12] 邓洪伟, 赵春生, 贾东兵, 等. 航空发动机喷管隐身修形设计技术分析[J]. 航空发动机, 2014, 40(2): 10-14. https://www.cnki.com.cn/Article/CJFDTOTAL-HKFJ201402004.htm DENG Hongwei, ZHAO Chunsheng, JIA Dongbin, et al. Analysis of stealthy shape design technology for aeroengine exhaust nozzle[J]. Aeroengine, 2014, 40(2): 10-14. https://www.cnki.com.cn/Article/CJFDTOTAL-HKFJ201402004.htm
[13] 李惠萍, 周起勃, 匡定波. 基于粒子系统和C-G法的尾焰红外模拟[J]. 红外与激光工程, 2008, 37: 601-603. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ2008S2062.htm LI Huiping, ZHOU Qibo, KUANG Dingbo. Simulation of infrared radiance of aircraft plume[J]. Infrared and Laser Engineering, 2008, 37: 601-603. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ2008S2062.htm
[14] 许爱华, 汪中贤, 于坚, 等. 高空高速无人机尾焰红外辐射特性计算研究[J]. 红外与激光工程, 2012, 41(7): 1700-1707. DOI: 10.3969/j.issn.1007-2276.2012.07.005 XU Aihua, WANG Zhongxian, YU Jianxian, et al. Numerical simulation of infrared radiation characteristics for plume of the high-altitude and high-speed UAV[J]. Infrared and Laser Engineering, 2012, 41(7): 1700-1707. DOI: 10.3969/j.issn.1007-2276.2012.07.005
[15] 张建奇. 红外物理[M]. 西安: 西安电子科技大学出版社, 2013. ZHANG Jianqi. Infrared Physics[M]. Xi'an: Xidian University Press, 2013.