前体附面层对轴对称收—扩喷管内外流场数值模拟影响.docVIP

前体附面层对轴对称收—扩喷管内外流场数值模拟影响 　　【摘 要】 在对有较高外流速度的轴对称收-扩喷管内外流场进行数值模拟的过程中，常需考虑实际的模型前体部分所产生的附面层对喷管部分的流场计算结果所可能产生的影响。采用三种常见的模型前体附面层处理方式，通过对包括严重过膨胀在内的三种工况下的喷管进行数值模拟，比较了这三种处理方式对于计算结果的影响。结果表明：在较高外流速度下，模型前体的附面层不可忽略，喷管外壁面附近流场的计算结果受到其很大的影响。对于喷管内部流场，在过膨胀状态下，外流的计算结果差异也会反映到内流激波后的亚音速区域。在模型前体构形较为简单的情况下，采用计算平板湍流附面层的计算公式来模拟前体附面层的速度分布，也可以满足对于后体-喷管数值计算的精度要求。 　　【关键词】 轴对称喷管 附面层 数值模拟 内外流场 前体 　　1 引言 　　在轴对称收-扩喷管的研究中，除模型的风洞试验[1]外，如今CFD技术也得到了大量的应用，其中已经出现了一些对亚音速、跨音速状态下，以及设计工况与非设计工况下喷管内外流场的数值模拟算例[2][3][4]。 　　考虑到风洞试验中，完整的喷管试验模型往往具有较长的前体部分（如图1，喷管收缩段之前的锥形体的大部分区域都可看做是前体），当近壁面的气流到达测量点之前，都需要流经这部分前体。与平板类似，前体的长度会对测量点的近壁面气流即附面层的速度型有明显的影响，因此对于数值模拟而言，前体部分附面层的模拟结果也将会对所要研究的喷管部分的计算结果造成一定的影响。本文对高亚音速下的轴对称收-扩喷管的三种工况进行了数值模拟，对比了三种常用的前体附面层处理方式，并都与试验数据进行了比较。 　　2 数值模拟 　　本文采用的是商用CFD软件Fluent，数值模拟方法为时间推进的有限体积法，控制方程为一般曲线坐标系下强守恒形式的N-S方程。离散格式选用隐式二阶迎风格式，湍流模型为RNG两方程模型。为保证对于非设计工况下喷管内外流场的计算精度，近壁面的处理采用的是增强型壁面函数[3]。 　　研究对象为轴对称喷管[1]（图2），设计落压比???4.25，喷管进口直径为128.63mm，喉道直径为76.2mm，喷管出口直径为85.1916mm。 　　计算采用二维结构化网格，所计算的喷管工作工况为三种情况，分别为：马赫数0.9，NPR=5.09；马赫数0.9，NPR=2.03；马赫数0.8，NPR=1.50。三种工况分别都采用了三种方式进行了计算，即： 　　（1）采用带完整前体的整体模型； 　　（2）外场进口采用UDF进行编制的1/7速度分布律[5]的省略前体的模型，其中1/7速度分布律的公式为： 　　（3）不采用UDF（外场进口速度为单一常数）的省略前体的模型。 　　其中，（1）为单独一套网格（图3），（2）与（3）采用了另一套网格（图4）。另外，针对不同的工况下存在激波的位置进行了网格加密，以确保计算结果的准确性。其余的计算条件为，环境压力101325Pa，环境总温300K。 　　3 结果分析 　　3.1 喷管相关性能参数的定义 　　喷管的流量系数用下式描述[6]： 　　喷管的推力用下式描述[7]： 　　其中为喷管出口处射流在轴向方向上的速度分量，为喷管出口处的压力，为环境压力。 　　推力系数为[6]： 　　压力系数为[1]： 　　其中为壁面静压，为来流速度。 　　压阻系数为[1]： 　　其中为轴向方向上的面积投影，为后体-喷管的最大横截面积。 　　3.2 Ma=0.9，NPR=5.09时的喷管壁面压力 　　图5为外流Ma数为0.9的情况下，后体起始点处的近壁面速度分布。可以看到采用UDF的情况下，由1/7速度分布定律算出的速度分布曲线与通过数值模拟整体模型算出的速度分布曲线相对比较接近。而不使用UDF的情况下，其速度分布曲线则完全不同。 　　图6与图7为Ma数0.9，NPR=5.09的工况下，喷管内外壁面压力分布图。由图来看，外壁面的压力分布，在使用UDF的情况下，计算结果与试验数据在大部分区域吻合得最好，而采用整体模型的计算结果次之，不使用UDF的计算结果最差。三种计算方式的结果差异最大处在压力突变处（激波位置）与其之后的存在逆压梯度的区域，这都是由于收缩段之前的外壁面附面层速度分布由于计算方式的不同而存在着一定差异，这种差异对于大部分区域没有顺压梯度的外壁面收缩段来说，影响是比较大的，这一点在表1中的收缩段压阻系数这一项中也有体现。而内壁面的压力分布计算结果，三者都差不多，这是因为NPR=5.09工况下，落压比已经高于该喷管的设计落压比，喷管内的气流在喉道音速线之后为超音速气流（图8），因而外场流动状况的计算差异无法影响到内场气流，其计算出的流量系数与推力系数也