超声速总压管气动特性数值研究.docVIP

超声速总压管气动特性数值研究 　　摘要：建立了一种总压管的计算模型，0度攻角时，在马赫为1.1～2.4的超声速均匀来流条件下，计算值与理论值吻合得很好。在此基础上，模拟了两支总压管的角度特性，研究了来流马赫数及头部形状的影响，计算结果也与已有实验值吻合得较好。超声速来流条件下，该计算模型能够很好地捕捉总压管前的激波，刻画激波的强间断特性。 　　关键词：超声速流场；总压管；角度特性；激波。 　　前言 　　总压管是测量流场总压的有效工具之一。NACA在90年代中期就对总压的角度特性展开了大量研究[1][2][3]，得到了不同头部几何形状及来流马赫数对总压管角度特性影响的一些规律。国内对总压管在超声速流场中的应用也进行了一些研究，文献[4]提出了一种超声速条件下利用不加修正总压探针测量值峰值连线近似总压分布测量值的上包络线法；文献[5]研究了总压管在超声速流场测量中的影响。但是，在国内关于总压管在超声速流场中的角度特性的相关文献还不多见。 　　本文建立了一种总压管的计算模型，首先模拟了0度攻角时，某一总压管在不同马赫数的超声速均匀来流条件下的“测量”情况，计算结果与理论值吻合得很好。在此基础上模拟了两支总压管的角度特性，研究了不同来流马赫数及头部形状的影响，计算结果也与已有实验值[1]吻合得较好。超声速来流条件下，该计算模型能够很好地捕捉总压管前的激波，刻画激波的强间断特性。 　　1、计算模型的建立 　　1.1总压管的模化 　　图1给出了研究对象总压管A和B头部的几何参数。A、B总压管均具有圆柱形外形和圆锥形内腔，外径D都为2mm，内径d都为0.375mm（d/D=0.1875），两者唯一的区别是锥形内腔的半锥角，A、B总压管锥形内腔的半锥角分别为15°和25°。 　　图1 总压管A和B的头部几何形状 　　计算之前对总压管进行了模化，计算中去除了总压管支杆，保持总压管的头部几何不变，总压管长度取为管外径的5倍。总压管的测量原理就是将气流等熵滞止，并将滞止压力连接到传感器上而得到气流的总压，此处的做法是保留了一段引压通道并将通道的末端封死。 　　1.2计算域的选取 　　计算域的进口与探针前端距离为总压管外径的20倍，这是考虑到所计算的来流为超声速，扰动对上游的影响较小。计算域出口与总压管尾部的距离是在试算后认为不再有边条反射而确定的，最终计算域出口距总压管尾部的距离为40倍的管外径。考虑超音速风洞的堵塞比，同时考虑实际应用尺寸，计算域的外径定为总压管外径的75倍，图2为最终选取的计算域。 　　图 2 计算域 　　2、计算网格 　　采用激波捕捉和激波装配型两种不同类型的计算网格，如图3所示。激波装配型网格是在捕捉型网格计算结果的基础上建立起来的，具体为：根据激波捕捉模型网格的计算结果，初步判断激波的位置，并在捕捉模型网格的基础上进行局部加密得到。实践证明激波捕捉模型的网格能够更准确地反映激波形态，从而获得正确的计算流场[3]。为满足计算湍流模型SST的需要，网格在壁面处都进行了加密，第一层网格为0.005mm，Y-plus小于20。 　　图 4 计算网格 　　3、边界条件设置 　　进口：给定总温、总压和速度的大小与方向。 　　出口：1）亚声速时，给定出口静压；2）超声速时，直接设为超音（不给边界条件）。 　　计算域周向：给滑移壁面边条。 　　总压管的壁面：给无滑移壁面边条。 　　4、结果分析 　　4.1零度攻角计算验证 　　选择了总压管A为计算对象，来流马赫数分别取1.1、1.2、1.3、1.4、1.62、2.0和2.4，来流攻角均取为0度。从后面图7中部分马赫数条件下数据提取线AB上的总压分布情况也可以看出，引压通道入口不远气流的总压就保持不变。为了方便处理，直接读取总压管引压通道后半段AB线上的总压作为总压管“测量值”。总压的理论值定义为来流经过一道正激波后的总压。 　　图5为总压管“测量值”在总压与马赫数理论关系线上的分布情况。从图中可以看出“测量值”与理论值几乎是重合的。总体上看，计算结果与理论值吻合较好。 　　图5 总压管“测量值”在总压与马赫数理论关系线上的分布 　　为了说明总压管在超声速流场中实现总压测量的过程，提取了总压管轴线从计算域进口到总压管引压通道末端之间（设为AB线上）的马赫数和总压，并作了分析。 　　图6是不同来流马赫数下AB线上的马赫数分布情况。从图中可以看出，超声速来流的马赫数在总压管上游的某个位置发生了突变，从超声速变为亚声速，说明此处产生了一道正激波。马赫数变化的程度反映出了激波强度的大小，对比不同来流马赫数下AB线上的马赫数变化情况可以看出，随着来流马赫数的增大激波逐渐增强。气流经过正激波后没有直接滞止下来，而是经过一段距离直到引压通道入口才实现完全滞止。对比不同来流马赫数下AB线上的马赫数分布也可