尾焰红外辐射特性计算方法研究.docxVIP

尾焰红外辐射特性计算方法研究 1 飞机尾火焰的计算 针对耀斑特征的计算，研究了许多闪闪发光特性的计算。由于近些年来对于飞机红外隐身性能要求的提高,人们也开始关心飞机尾焰的红外辐射特性。飞机尾焰的计算也可以采取类似于计算液体火箭尾焰的方法进行。在液体火箭和飞机尾焰红外辐射特性的计算中,为了简化计算,人们常常忽略粒子的影响,将尾焰作为纯气体计算。 一般说来,尾焰红外辐射特性的计算基本遵循如图1所示的步骤进行。 在每一部分的计算中,都有很多方法可以选择,下文将对每部分使用的主要方法进行概述,重点介绍每种方法的现状及特点,从而给出本文推荐使用的方法。 2 分浓度场的计算方法 尾焰流场的计算主要包括温度、压力场及组分浓度场的计算,归纳起来主要有两类方法,一是采用理想模型、半经验公式进行计算。二是利用计算流体动力学CFD软件计算。 2.1 初始段和主段 该方法假设尾焰为理想的湍流射流,利用理想的湍流射流模型进行计算。如图2所示,该方法的特点是将尾焰分为初始段(x≤x1)和主段(xx1)两个部分,并且假设初始段存在一个组分压强和温度恒定的核心区(斜线区),初始段和主段采用不同公式计算。该计算方法计算公式简单,能很快得到尾焰流场参数,且可以与尾焰辐射算法程序很好的结合,所以很多文献使用了这种方法。 2.2 尾火焰辐射特性计算方法 目前CFD软件有很多,例如CFX、FLUENT、PHO ENICS、STAR-CD等,国内应用较多的是FLUENT。为了简化计算,一般对于简单的圆喷管只进行二维数值模拟,再利用其对称性将尾焰流场数据扩展到三维空间中。 在实际飞行中,尾焰流场受到很多因素的影响,如飞行目标形状、海拔以及喷管形状(如二元喷管)等,所以如果进行更为细致的研究,应采用三维的数值模拟。 利用CFD软件计算的优点是显而易见的,例如可以进行三维数值模拟,可以选择湍流模型,最终可以获得温度场、压力场、组分浓度场等更为详细的流场参数。但是利用CFD软件也有一定缺点,很多尾焰红外辐射特性计算方法不能很好的与CFD软件计算的流场相结合。所以,国内在尾焰红外辐射特性的计算研究中往往利用理想湍流射流模型或者半经验公式进行流场参数的计算。 3 谱带模型法 尾焰辐射参数包括气体辐射参数和粒子辐射参数的计算,粒子辐射参数计算目前广泛应用的是Mie理论。气体辐射参数计算主要有三种方方法:逐线计算法、谱带模型法、κ-分布方法。逐线计算法是逐条计入气体吸收谱线贡献的一种精确计算方法,但是其计算速度慢,计算成本高,所以工程中一般不使用这种方法。谱带模型法是将某一波数间隔内的光谱线排列、重叠性质与单个谱线性质联系起来的计算透过率的方法。该方法计算简单,速度快,如今依然有很多文献中采用该方法。其缺点是:如果不借助其他近似,则不能处理重叠吸收带,而且不能处理具有散射的辐射传输问题。κ-分布方法假定大气吸收系数在不同压力和温度下存在一一对应关系,由此可以将对波数的积分转化为对吸收系数的积分,大大简化了计算。该方法的优点是精度高速度快,可以较好地处理重叠吸收带,以及同时处理吸收与散射问题。所以本文推荐在尾焰辐射相关参数采用κ-分布的方法来计算。 4 尾火焰红外辐射特性模拟方法 通过求解辐射传输方程来研究尾焰红外辐射特性是主要的计算尾焰的红外辐射特性的方法,归纳起来可以分成两大类,一类是不考虑固体粒子吸收与发射的准气体模型,准气体模型使喷焰辐射的计算大为简化;另一类是考虑气体及固体粒子的吸收、发射,还考虑固体粒子散射的模型。 在液体燃料飞行目标尾焰中,粒子成分比较少,人们为了简化起见,往往忽略粒子对于辐射影响。所以很多文献中依然采用准气体模型来解决这类问题,例如目前依然有很多文献中使用NASA推荐的C-G近似谱带模型法计算液体火箭及飞机尾焰红外辐射特性模拟。这在精度要求不高的理想状态下是可行的。 然而,随着对尾焰红外辐射特性测量水平及数值模拟精度要求的提高,人们认识到在液体燃料飞行目标尾焰红外辐射特性的模拟中有必要加入粒子对于辐射的影响,这样得到结果与实际测量结果更相符。如图3和图4为文献中给出的某型Atlas火箭在海拔40 km处飞行中测量与计算结果的比较。可见加入碳粒子后的计算结果与测量的光谱特征更相符合。 因此在后文中我们主要介绍了可以考虑粒子对于辐射影响的三种典型方法,包括热流法、有限体积法及反向蒙特卡洛法。这里的选择是有根据的,美国目前依然在使用的标准红外辐射模型(SIRRM)就是基于双热流法及六热流法发展起来的。近年来发展起来的有限体积法计算工作量较小,计算精度高,而且对非规则计算域的适应性强,所以得到愈来愈多的重视和研究。反向蒙特卡洛法在处理小区域或小立体角接收到的辐射能量问题是很有效,效率远高于正向蒙特卡洛法。它属于概率统计的方法,所以理论上讲,只要模拟的次数足够多,