高超声速风洞轴对称喷管收缩段设计

Φ240 mm高超声速风洞扩压器设计研究范孝华;杨波;朱涛;钟俊【摘要】采用一维管流方法对口径为Φ240 mm、马赫数6.0的自由射流高超声速风洞扩压器进行设计,运用数值模拟方法对扩压器结构形式进行了筛选,并通过试验验证扩压器设计方法的有效性.结果表明:综合考虑风洞流场品质、扩压器抗反压能力以及扩压器总压恢复等指标,确定单级15°-0.8D e扩压器方案为最优设计方案;该方案扩压器的试验与数值模拟结果吻合较好,能够使风洞正常起动,并长时间维持风洞正常运行;筛选得到的扩压器方案可用于风洞建设,其设计方法可为同类高超声速风洞扩压器的研制提供参考.%1-D tube flow method was adopted for the design of Φ 240 mm hypersonic wind tunnel diffuser at Mach Number6.0,and numerical simulation was carried out in order to select the configuration of the diffuser.Then,experimental test was used to verify the design method of the diffuser. The results show that, in comprehensive consideration of the flow-field performance of the wind tunnel,the resistance pressure ability,and total pressure recover,the diffuser has single stage compression half-cone angle being 15°and second throat diameter being 0.8D e (the exit diameter of the nozzle).The experimental results of the optimal diffuser agree well with corresponding numerical ones.The hypersonic wind tunnel can start smoothly and maintain long-time operation with this diffuser.The optimal diffuser design can be implemented in the construction of new wind tunnel,and the design method of the diffuser provides a reference for the development of similar diffusers in hypersonic wind tunnel.【期刊名称】《空气动力学学报》【年(卷),期】2017(035)005【总页数】8页(P633-639,644)【关键词】自由射流;高超声速风洞;扩压器;抗反压能力【作者】范孝华;杨波;朱涛;钟俊【作者单位】中国空气动力研究与发展中心超高速空气动力研究所,四川绵阳621000;中国空气动力研究与发展中心超高速空气动力研究所,四川绵阳 621000;中国空气动力研究与发展中心超高速空气动力研究所,四川绵阳 621000;中国空气动力研究与发展中心超高速空气动力研究所,四川绵阳 621000【正文语种】中文【中图分类】V235.21高超声速风洞是进行高超声速飞行气动试验和研究的重要地面设备，对于暂冲吹吸式高超声速风洞，扩压器的主要作用是在总压损失尽可能小的条件下将超声速气流减速增压以提高扩压器出口压力，减小由于下游背压不断升高对风洞流场带来的影响，维持高超声速风洞在尽可能长的时间内正常运行[1]。

前体附面层对轴对称收 【摘 要】 在对有较高外流速度的轴对称收-扩喷管内外流场进行数值模拟的过程中，常需考虑实际的模型前体部分所产生的附面层对喷管部分的流场计算结果所可能产生的影响。采用三种常见的模型前体附面层处理方式，通过对包括严重过膨胀在内的三种工况下的喷管进行数值模拟，比较了这三种处理方式对于计算结果的影响。结果表明：在较高外流速度下，模型前体的附面层不可忽略，喷管外壁面附近流场的计算结果受到其很大的影响。对于喷管内部流场，在过膨胀状态下，外流的计算结果差异也会反映到内流激波后的亚音速区域。在模型前体构形较为简单的情况下，采用计算平板湍流附面层的计算公式来模拟前体附面层的速度分布，也可以满足对于后体-喷管数值计算的精度要求。 【关键词】 轴对称喷管 附面层 数值模拟 内外流场 前体 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），因而外场流动状况的计算差异无法影响到内场气流，其计算出的流量系数与推力系数也都基本相同（表1）。 3.3 ma=0.9，npr=2.03时的喷管壁面压力与流场结构 此状态下，后体起始点处的近壁面速度分布依旧如图5所示。而图9与图10则为ma数0.9，npr=2.03的工况下，喷管内外壁面压力分布图。与之前的工况类似，外壁面压力分布的计算结果与试验数据的对比由好到差，依次为采用udf、采用整体模型、无udf，其原因也与之前基本相同；对于内壁面的压力分布，同样出现与外壁面类似的情况，即在采用不同计算方法，压力突变位置以及之后的区域的压力分布也是不同的，具体来说，无udf的情况下，计算结果与试验数据差距较大，而采用udf或整体模型的情况下，计算结果与试验数据较为接近，其原因则是在npr=2.03的情况下，喷管处于较为严重的过膨胀状态，喷管的扩张段会出现一道较强的正 激波（图11），使得波后的扩张段气流降至亚音速，而外流场的计算结果则可以影响到这部分亚音速气流。由内外流场的压力分布图可以看出，在靠近喷口的外壁面处，计算结果的压力较高时，内壁面的激波后部压力也相对较高，从而使得激波的位置也更加靠近喉道。由于激波存在于扩张段，因而喷管喉道处的音速线都存在，而从表2与表1的对比可以看出，流量系数未受到影响，但是推力系数则降低不少，这是因为正激波的存在使得气流的总压明显降低，这将直接影响到喷管的推力。而表2内推力系数的差异则是由于各个前体附面层处理方式计算出的激波位置和波后压力的不同而产生的。 3.4 ma=0.8，npr=1.50时的喷管壁面压力与流场结构 图6为外流ma数为0.9的情况下，后体起始点处的近壁面速度分布。可以看到，在此速度下，由1/7速度分布定律算出的速度分布曲线同样可以和通过数值模拟整体模型算出的速度分布曲线较好地匹配。 图13与图14为ma数0.8，npr=1.50的工况下，喷管内外壁面压力分布图。在此工况下，外壁面已经没有了很激烈的压力突变点，从与试验数据的对比来看，依然是在无udf的情况下计算结果最差。而采用整体模型和采用udf的情况下，计算结果则都和试验数据相差不远，区别主要是在靠近喷口的壁面附近，整体模型计算出的压力值要高于采用udf计算出的压力值。同样这种压力差也影响到了内壁面的压力分布计算结果，即：喷管扩张段的壁面在激波后的区 域里，采用整体模型的压力计算结果同样高于采用udf的计算结果，而且激波位置也相对更加靠近喉道。由于激波依然存在于扩张段（图15），与npr=2.03的情况类似，从表3可以看出流量系数三者依旧基本相同，而推力系数由于正激波造成的总压损失而与表1相比下降不少。其在表3中的三个算例之间的推力系数差异则因为计算出的激波位置不同而差别明显。 4 结语 本文通过对喷管在不同工况下的内外流场进行数值模拟，对比了三种不同的模型前体附面层处理方法，进一步得出了以下结论： （1）在对喷管内外流场进行数值模拟时，对于外壁面的压力分布，由于只处于零压梯度、逆压梯度或者很小的顺压梯度下，因而附面层对其计算的结果具有较为明显的影响； （2）对于内壁面的压力分布，当喷管处于正常工况或是欠膨胀工况时，喷管的扩张段的流场为超音速，此时外流场的计算结果几乎不会影响到喷管内流场；而当喷管处于过膨胀特别是严重过膨胀工况时，喷管扩张段出现激波以及激波之后的亚音速流场区域，这两者的计算结果都会较为明显地受到外壁面附近流场特别是靠近喷口处的流场的影响，从而在激波位置和壁面压力方面都存在比较明显的差异。 （3）在外流速度较大的情况下，忽略前体部分附面层的影响将无法得到令人满意的计算结果。若是考虑前体的附面层，当前体的构型较为简单时，使用1/7速度分布律来定义外场的进口速度，其 计算结果已经可以达到比较满意的精度。 （4）当通过采用完整的前体模型来引入附面层的影响时，本文的计算结果反而不如使用1/7速度分布律时的结果，其中的原因可能涉及到湍流模型的计算精度问题。在本文的算例中，前体占了整个计算模型的很大一部分，而计算所关注的后体收缩段所占比例很小，因而单纯用于工程计算用的湍流模型在整个模型的计算上或许满足精度要求，但是只抽取一小部分来看的话可能还是存在精度不足的现象。这一猜测有待进一步更深入的研究。 参考文献： [1]george t.carson，and edwin e.lee.experimental and analytical investigation of axisymmetric supersonic cruise nozzle geometry at mach numbers from 0.60 to 1.30 [r]. nasa tp-1953，1981. [2]黄宏艳，王强.过膨胀状态下轴对称收-扩喷管内外流场计算及分析[j].航空动力学报，2007，22（7）：1069—1073. [3]余铭，刘友宏.轴对称收-扩喷管内外流场一体化数值模拟[j].科学技术与工程，2011，11（32）：7979—7984. [4]reza ghias.simulation of flow through supersonic cruise nozzle： a validation study [r].tfaws paper session，2011. [5]陈懋章.粘性流体动力学基础[m].北京：高等教育出版社，2002.

基于轴线马赫数分布的喷管扩张段无粘型面设计胡振震＊，李震乾，石义雷，陈爱国【摘要】摘要：针对高超声速风洞轴对称喷管设计问题，开展了喷管扩张段无粘型面设计研究。

介绍了基于预设轴线马赫数分布的直接设计方法，改进了基于面积比的轴线马赫数分布预设方法，提出了一种方便多点控制的轴线特征点分布方法。

对设计喷管流场进行特征线网三角化，与数值模拟结果进行比较，并分析了影响喷管无粘型面的设计因素。

表明：改进的面积比方法可以保证轴线马赫数分布预设的合理性；轴线马赫数分布、轴线上特征点分布和边界特征点数明显影响喷管无粘型面。

【期刊名称】实验流体力学【年(卷),期】2016(030)004【总页数】8【关键词】高超声速风洞；扩张段无粘型面；直接设计；轴线马赫数分布0 引言高超声速风洞喷管设计目前普遍采用无粘型面加边界层修正的直接设计方法，相比于优化设计［1－3］，该方法原理简单计算便捷，在广泛应用中取得了成功。

然而在实践中有时也发现一个问题，即风洞运行的实际流场与喷管设计流场存在差异。

该问题的产生因素有很多，包括来流非均匀非稳态、非平衡膨胀、边界层转捩、热化学非平衡、气体组分冷凝、加工和组装精度、喷管内表面侵蚀和型面的力热变型等，但喷管的无粘型面设计则是其中最重要的一个方面，因而分析其差异应首先确认喷管无粘型面设计方面的问题。

无粘型面直接设计方法根据原理可分为3类：基于预设初始膨胀段型面的方法［4］；基于泉流（或源流）假设的方法（如Foelsch［5］方法、Cresci［6］方法和Sivells［7－8］方法）和基于超声速区预设轴线马赫数分布的方法（如ASN［9］方法）。

研究表明，第2类基于泉流假设的方法（特别是Sivells方法）在选择合适的设计参数情况下可以得到良好的无粘喷管设计，但该类方法设计的喷管长度较长，调整不易。

第1类和第3类方法均需计算完整的超声速区特征线网，区别在于第1类方法依据预设的初始膨胀段型面确定轴线上核心区的马赫数分布，而第3类方法则根据预设的轴线马赫数分布反求喷管型面。

轴对称收-扩喷管内外流场一体化数值模拟余铭;刘友宏【摘要】基于求解Navier-Stokes方程组对某航空发动机的轴对称固定几何收-扩喷管进行了六种工况的内外流场一体化数值模拟,计算出了喷管的流量系数和推力系数.结果表明,收-扩喷管超音速射流,在地面工况下表现为过膨胀特征,在高空不同工况下表现为欠膨胀特征,形成了复杂的激波系和膨胀波系.在欠膨胀状态,随着落压比的增大,喷管的流量系数无明显变化,但推力系数显著减小.另外,从理论解出发,计算了前面各种工况的喷管出口气动性能参数,并与数值模拟结果进行了对比分析,验证了数值模拟结果的可靠性.%Entity with a tunnel project, an integrated approach to the working principle and technical characteristics of TST geological prediction and LTD2100 ground penetrating radar system is described. The rationality of design and construction guidance and recommendations in the construction are given. Integrated tunnel geological prediction and efficient implementation of construction are provided a scientific basis.【期刊名称】《科学技术与工程》【年(卷),期】2011(011)032【总页数】6页(P7979-7984)【关键词】超音速射流;轴对称收-扩喷管;数值模拟;推力系数;流量系数【作者】余铭;刘友宏【作者单位】北京航空航天大学能源动力工程学院,航空发动机气动热力国家级重点实验室,北京 100191;北京航空航天大学能源动力工程学院,航空发动机气动热力国家级重点实验室,北京 100191【正文语种】中文【中图分类】V231.3轴对称收-扩喷管的超音速射流一直是喷管技术研究中的一个热点。