低温风洞运行压比相关性研究及应用

黄知龙，王 宁，廖达雄，*

(1. 南京航空航天大学 非定常空气动力学与流动控制工业和信息化部重点实验室，南京 210016；2. 中国空气动力研究与发展中心 设备设计与测试技术研究所，绵阳 621000)

0 引言

风洞试验的雷诺数模拟不足,将会使边界层分离、旋涡流动、激波/边界层干扰、激波/旋涡干扰等黏性起支配作用的流动现象发生变化，使试验数据与真实飞行存在明显差异[1]。研究表明降低试验介质气流温度目前被认为是提高风洞试验雷诺数最有效的技术途径，低温风洞概念的发展极大拓展了亚跨声速风洞的试验能力。美欧建成的大型低温风洞NTF和ETW为波音和空客的创新发展发挥了不可替代的重要作用[2-6]。伴随我国低温高雷诺数风洞设计建设，相关技术研究成为热点[7-12]。

低温风洞具有复杂的风洞控制问题，热力学是一个高度耦合、多变量和非线性过程。主要存在以下困难：一是大范围精确的流场模型难以获得；二是系统存在强烈非线性耦合；三是系统存在大时滞效应，特别是温度。兰利研究中心早期建设的NTF采用了非线性增益调度PI控制[13]，欧洲ETW风洞采用一种具有自学习能力的控制算法不断自我优化[14]。国内祝汝松等探索了鲁棒自适应控制[15-16]，刘为杰等开展了自抗扰控制研究[17]，兰其龙等开展了基于神经网络预测控制[18]，黄丽娟等开展了不协调目标信息系统多变量控制方法仿真[19]。结果表明现代控制理论以及智能控制由于其算法的复杂性，目前结果并不理想，调节时间较长，还不能够完全满足低温风洞运行控制需要。

风洞运行压比定义为风扇段出口截面和入口截面的总压比值，其与风扇转速、液氮喷射流量及排气流量深度耦合，是影响风洞控制的核心参数。本文以中国空气动力研究与发展中心新建成的0.3 m低温跨声速风洞测试数据为研究对象，探索建立了风洞运行压比与不同运行状态参数（压力、温度、马赫数）间统一的动态传递模型，可应用于风洞主控系统优化，通过有效前馈控制，缩短了运行参数调节和稳定时常。

1 研究设备简介

0.3 m低温跨声速风洞是一座连续式单回流增压风洞，利用多级轴流低温风扇驱动。试验段截面尺寸325 mm（宽）× 275 mm（高），运行马赫数范围0.15～1.30，总温范围110～323 K，总压范围100～500 kPa。风扇由1.5 MW的外置电机和变频器驱动，转速范围为500～7200 r/min。洞体材料为不锈钢，表面覆盖绝热材料以减小与环境之间的热交换。通过向洞体内喷入液氮气化吸热平衡风扇运行功率和洞壁传热，液氮喷射流量范围0.6～6.5 kg/s。同时在风洞第三和第四拐角段之间配置排气系统将气化后的氮气排出实现洞体压力控制。风洞回路试验现场照片见图1。

图1 0.3 m低温跨声速风洞Fig. 1 A photo of the 0.3 m cryogenic transonic wind tunnel

2 测试数据

风洞运行压比通过位于风扇段上下游截面周向均布的总压排架测量。入口截面为水平布置的两个5点总压排架，出口截面为周向均布的四个4点总压排架。压力传感器为PSI 9116扫描阀，量程为45 PSI，测量精度0.05%。并按环状面积积分计算得到压力的平均值。压力测量时，试验段内除了轴向静压探测管外为空风洞，同时试验段下游第二喉道完全拉开、支架段再入调节片开度5°。风洞运行状态参数则通过安装于试验段的皮托管获取总静压、总温探针测温度。压力测量采用Mensor CPT6180高精度绝压传感器，量程为70 PSI，测量精度优于0.02%。温度测量采用Lakeshore PT-102-14L传感器，测量精度2.49 mK。

图2给出了风洞三个运行总压（115 kPa、200 kPa、450 kPa）、五个总温(110 K、150 K、200 K、280 K、323 K)和12个马赫数的组合状态校测点数据，共计119个组合，马赫数范围0.20～1.30，阶梯ΔMa= 0.10。其中，各孤立标示点的纵坐标为压比测量值，以变量ε表示；横坐标为基于风洞运行状态参数和试验段水力直径为特征尺寸的雷诺数，以Rec表示。每条曲线则为同一试验段马赫数下压比随特征雷诺数变化的多项式拟合曲线，所有测试状态雷诺数最小约为1.23×106，最大约55.3×106，变化范围达到45倍。

从图2中可以看出，对于所有给定的试验段马赫数，风洞运行压比总体上呈现出随着试验雷诺数增加而逐渐降低的趋势，且在雷诺数较低区域压比下降较快，而随着雷诺数的增加其变化率逐渐减小。笔者认为风洞运行压比的降低主要是由于雷诺数增大时回路内气体的黏性损失降低所致。且当雷诺数增大到一定程度后洞壁内附面层的气流流动全部发展为湍流附面层，管道黏性损失进一步降低则变得困难。各马赫数下测试数据随着雷诺数增加缓慢渐变还表明，导致这种结果最主要原因来自于气体黏性损失的降低，而不是由于风洞内部段局部气流流动状态的改善所致。

图2 风洞运行压比与雷诺数关系Fig. 2 The variation of the running pressure ratio with Reynolds number at different Mach numbers

另外需要指出虽然图中马赫数Ma= 1.1、1.2和1.3的测试结果也表现出与亚声速相似的变化趋势，但此时风洞试验段、支架段和第二喉道等部段气流已达到超声速，运行出现阻塞，进一步增加风扇转速和驱动功率并不会使试验段马赫数再增加。在这种情况下，风洞运行压比与雷诺数的关系则不再符合上述分析。

3 动态模型建立

3.1 数据拟合分析

为了更全面的理解和预测低温风洞运行特性，为低温风洞运行参数的耦合自动控制方案优化提供技术支撑，需要以尽可能简洁的形式表达出风洞运行压比与试验参数的关系。为此对测试数据尝试进行拟合分析。

采用多项式拟合将风洞运行压比与风洞其它变量联系起来是一种最常用数据处理方法[20-21]。然而，由于低温风洞运行状态（温度、压力和马赫数）变化非常宽广，要将各参数之间的关系以方程的形式表示且达到理想精度，方程将变得非常复杂，实现起来也较为困难，因此需要探索新的解决途径。而现有的测试数据结果表明了风洞运行压比是试验段马赫数和特征雷诺数这两个变量的函数，能否构造出一个简单的函数表达式描述这三个变量之间的关系，就成为进一步分析的方向。

因为气流流动的压力损失与速压成正比，也就是与速度的平方成正比，因此可将测试数据进一步整理为运行压比与试验段马赫数平方的关系，见图3所示。其中横坐标为马赫数平方，纵坐标为运行压比。从图中可以看出风扇的压比与马赫数的平方成近似线性关系。假设压比与雷诺数为简单的幂函数，则可通过改变幂函数指数大小获得相对最佳数据拟合结果[22]。

图3 风洞运行压比与马赫数平方关系Fig. 3 The variation of the running pressure ratio with the square of the Mach number

图4 风洞运行压比与马赫数-雷诺数组合关系Fig. 4 The variation of the running pressure ration as a function of Ma2Rec−0.095

我们再针对上图中线性度较好的区域数据（Ma≤1.0），采用最小二乘法进行了线性拟合[23]，见图5所示。直线斜率为0.9218，与纵坐标截距为1.0033。拟合优度可决系数R2为0.9987，表明拟合程度良好。因此，可给出风洞运行压比与组合函数的构造关系表达式如下：

图5 测试数据线性拟合Fig. 5 A linear fitting of test data

式中，Ma为试验段马赫数，Rec为特征雷诺数。

依据风洞试验数据的不确定度评估理论，对上述测试马赫数和压比的不确定度进行分析。以风洞典型运行测试状态（总压pt= 115000 Pa，总温T0= 110 K，马赫数Ma= 0.90）为例进行分析。主要运行参数包括：驻 室 参 考 静 压pc= 67994 Pa，风 扇 段 入 口 总 压pFI= 99958 Pa，风扇段出口总压pFO= 116121 Pa。

自由流马赫数Ma是测量值总压pt和驻室参考静压pc的函数，通过状态方程式对各参数求导得到偏导数关系式，带入参数可得到偏导数值如下：

因为测试pt和pc的传感器通过校准系统得到的校准不确定度是固定的，可作为测量值的偏差极限进行传递，pt和pc的相关偏差极限为：

马赫数的不确定度UMa由 总压偏差极限Bpt、驻室参考静压偏差极限Bpc， 及相应的精度极限Ppt和Ppc确定，计算方程式如下：

DM的校准不确定度也作为偏差极限传递，没有伴随产生的精度极限，根据校准结果BDM取值0.00177。将各参量带入式（4）中，可得到马赫数不确定度为：UMa=±0.0019。

测试压比的不确定度Uε可表示为：

其中，

将各参量值代入式（5），则可得到压比不确定度为：Uε=0.002376。

3.2 理论分析

虽然上文中风洞运行压比与马赫数和雷诺数组合函数的关系是通过风洞的实测数据归纳得出，但从传统的流体力学理论也同样可以推导出类似的方程式，下文进行简单阐述。

风洞管道流动压力的无因次损失系数k定义式如下：

式中，v1为入口截面速度，p01为入口截面总压,p02为出口截面总压。

由于低温下的气体为不完全气体，需考虑真实气体效应，相关参数计算公式如下[2，24]：

气体密度：

速度：

动压：

雷诺数：

式中，R为气体常数，γ为比热比，Z为压缩性因子，µ为动力黏性系数。

利用连续性方程，并将式（10）带入式（7），则风扇段出口至入口的压力损失系数k的表达式为：

式中，γ为比热比，p1为风扇入口截面静压，p2为风扇出口截面静压，Ma为试验段马赫数。

假设风扇段前后静压比值用变量 ε表示，则有：

上式可变形为：

来流马赫数Ma< 1的无激波平板边界层流动阻力研究表明摩擦损失系数可表示为雷诺数的幂函数倒数和比例常数的乘积[25]，如下式：

式中，ko为比例常数（与风洞回路几何外形有关，而与雷诺数无关），指数n是雷诺数的函数。

将式（15）带入式（14），则得到风洞运行压比的理论表达式如下：

可以看出理论推导得到的关系式（16）和测试数据的拟合方程式（1）具有相同的表达型式。测试数据拟合得到的雷诺数指数（−1/n）值约为−0.095，介于−1/10和−1/11之间，基于平板附面层表面摩擦阻力得到的值介于−1/3到−1/9[25]，两者存在差异可能是因为风洞回路流道的几何外形更为复杂，且包含了拐角导流片、开槽壁板和支架等流道内部件，需进行修正。另外，常数项偏差0.0033可能是由于液氮喷射所致或数据系统误差引起的，需要开展进一步精确测量研究解决，但该动态模型对于工程应用已有足够精度。

4 动态模型应用

低温风洞运行控制多变量耦合关系示意见图6所示。最困难的是气流温度精确控制，它与压缩机转速、液氮喷射流量都是强耦合，动力学和热力学关系复杂。

图6 运行控制多变量耦合关系Fig. 6 A schematic of the Multi-variable coupled control

利用前文获得的压比与运行状态之间的传递函数（式（1））使得风洞运行的液氮需求和压缩机功率以相对简单的方式与试验运行状态相关联，表达式如下：

式中，Pf为风扇轴功率，Tt,i为风扇入口温度，β为液氮吸热能力，m˙LN2为液氮流量。

将该动态模型应用于0.3 m低温风洞控制方案，将使控制系统能够根据设定的风洞运行目标状态，提前获得风洞运行所需的精确液氮流量和压缩机转速，通过前馈控制则可提升控制系统效率，缩短运行状态参数稳定时长。

5 结论

利用0.3 m低温风洞的流场校测数据，成功构造马赫数和雷诺数的组合幂函数，并建立了马赫数小于1.0时风洞各运行工况条件下运行压比与马赫数-雷诺数的线性关联方程式，并利用传统的流体力学方程式推导获得了相同的数学表达型式。结果表明：风洞运行压比与马赫数-雷诺数组合函数存在简单线性关系，该研究成果作为动态模型应用于低温风洞多变量控制，可简化参数传递过程和优化方案，缩短状态参数的稳定时间。同时，由于马赫数大于1.0时风洞流道内存在激波、阻塞等流动状态，本文所得到压比关联式在超声速区域并不适用，对此还将进一步开展针对性研究。