0   引言

当飞机在云层中飞行时，云层中的过冷水滴会在飞机表面产生结冰现象，破坏其气动外形，影响气动性能，严重威胁飞行安全^[1-2]。鉴于飞机结冰的严重危害，民用航空适航规章要求制造商通过结冰防护等手段，确保飞机在结冰气象条件下的飞行安全。结冰风洞试验是目前重要的飞机结冰适航审定手段，已经广泛应用于飞机适航审定实践^[3]。

为表明结冰风洞具备开展结冰适航审定试验的能力，结冰风洞必须开展全面的流场品质评估，验证其流场品质符合结冰风洞适航应用的要求(需要指出的是：尽管民用航空适航规章并未直接对结冰风洞提出审定要求，但在飞机适航取证前，结冰风洞试验能力必须得到适航审定方的认可，因此结冰风洞仍需满足适航应用的相关要求)，热流场品质是其中一项重要的验证内容。近年来，随着国内大型结冰风洞——3 m×2 m结冰风洞——的建成^[4]、国产C919客机完成首飞并全面进入适航取证阶段，在C919客机适航取证的需求牵引下，结冰风洞流场符合性研究成为目前国内大型结冰风洞适航应用领域内的研究热点，其中热流场符合性是相关研究的重要内容。

针对结冰风洞热流场，世界上主要结冰风洞均开展了相关研究。为改善热流场品质，美国NASA Glenn IRT结冰风洞对制冷系统进行了多次升级改造^[5-11]，建立了系统的热流场品质评估方法，开展了多期热流场品质评估试验，试验结果不仅为制冷系统升级改造提供了数据支撑和改造思路，还为该风洞的适航应用提供了支撑。作为目前国际上尺寸最大、性能最完善的结冰风洞之一，意大利CIRA IWT结冰风洞^[12]也发展了配套的热流场校测设备和方法，开展了全面的热流场校测^[13]，为该风洞的适航应用奠定了基础。此外，美国波音BRAIT结冰风洞和Cox结冰风洞均开展了热流场品质评估试验^[14-15]，试验结果支撑了这些风洞的适航应用。国内目前尚未开展系统的大型结冰风洞热流场符合性验证研究，3 m×2 m结冰风洞缺乏适航应用的基础。

本文建立了结冰风洞热流场符合性验证方法，针对3 m×2 m结冰风洞主试验段构型，开展了热流场符合性验证试验，考察了试验段气流总温、试验段气流速度和喷嘴干空气射流对热流场空间均匀性和时间稳定性的影响，获得了试验段气流总温修正关系，形成了3 m×2 m结冰风洞主试验段热流场控制包线，为其适航应用奠定了基础。

1   3 m×2 m结冰风洞简介

中国空气动力研究与发展中心3 m×2 m结冰风洞是一座闭口回流式高亚声速风洞(图 1)，主要包括结冰喷雾系统、制冷系统、高度模拟系统和风机动力系统。结冰喷雾系统利用喷雾耙和喷嘴产生结冰云雾，其中喷雾耙安装于稳定段蜂窝器下游、收缩段入口位置，共20排，每排设置50个喷嘴安装位置，共计安装1000个喷嘴。制冷系统采用液氨作为制冷媒介，通过热交换器换热，控制气流温度，制冷系统最大功率约为11 MW。该风洞拥有3个可更换试验段，分别为主试验段、次试验段和高速试验段(各试验段尺寸见表 1)。本文仅针对主试验段构型开展热流场符合性验证研究。

图  1  3 m×2 m结冰风洞

Fig.  1  The CARDC icing wind tunnel

下载: 全尺寸图片 幻灯片

表  1  试验段尺寸参数

Table  1  The test section size parameters

Test section	Dimension
	Height/m	Width/m	Length/m
Maintest section	2	3	6.5
Secondarytest section	3.2	4.8	9.0
High speed test section	1.5	2	4.5

下载: 导出CSV 

| 显示表格

结冰风洞的总温探针为四线制PT-100铂电阻总温探针，探针温度测量范围为－50 ℃~50 ℃，精度为±0.2 ℃。8根总温探针安装于热交换器出口处，分别测量热交换器8个模块出口中心处的气流总温，探针位置如图 2所示。

图  2  结冰风洞总温探针安装位置示意图

Fig.  2  The position of total temperature probes mounted in the CARDC icing wind tunnel

下载: 全尺寸图片 幻灯片

2   结冰风洞热流场符合性验证方法

图 3为结冰风洞热流场符合性验证方法流程图，主要包括试验标准、试验内容、试验仪器、试验方法和数据处理等5部分内容，本节分别进行介绍。

图  3  结冰风洞热流场符合性验证方法流程图

Fig.  3  The flow diagram of the compliance verification method of thermodynamic flowfield of the icing wind tunnel

下载: 全尺寸图片 幻灯片

2.1   试验标准及内容

目前，国际上普遍依据国际自动机工程师学会(SAE International)发布的标准文件 Calibration and acceptance of icing wind tunnels(SAE ARP 5905-2003)验证结冰风洞流场符合性^[16]。该文件明确给出了结冰风洞在适航应用中需要满足的流场品质指标，其中的热流场品质指标如表 2所示。表中针对气流静温，分别给出了测试设备最大不确定度、空间均匀性和时间稳定性等3个符合性指标。本文依据SAE ARP 5905-2003标准，以气流总温为试验对象，开展热流场符合性研究。需要指出的是：由于气流总温可通过总温探针直接测量得到，进而可以根据气流速度计算得到气流静温，因此本文选择气流总温作为试验对象，更具有代表性。

表  2  结冰风洞热流场品质指标

Table  2  The quality index of thermodynamic flowfield of icing wind tunnel

Items	Static air temperature
	Ts < －30 ℃	－30 ℃ < Ts < 5 ℃
Measurement instrumentation maximum uncertainty	±2.0 ℃	±0.5 ℃
Spatial uniformity	±2.0 ℃	±1.0 ℃
Temporal stability	±2.0 ℃	±0.5 ℃

下载: 导出CSV 

| 显示表格

2.2   试验仪器

采用总温探针和温度格栅装置测量试验段内的气流总温，如图 4所示。

图  4  总温探针和温度格栅

Fig.  4  Total temperature probe and temperature grid

下载: 全尺寸图片 幻灯片

总温探针为四线制PT-1000铂电阻总温探针，如图 4(a)所示，探针长100 mm、外径6 mm，温度测量范围为－50 ℃~50 ℃，精度为±0.2 ℃。试验前对总温探针进行标定^[8]，获得了总温恢复系数，在马赫数0~0.6范围内，探针恢复率范围为0.994~1。

温度格栅装置由7个垂直栅条和4个水平栅条组成，如图 4(b)所示。垂直栅条厚20 mm、宽100 mm，栅条间距为375 mm。

垂直栅条前缘设置7个探针安装孔，如图 5所示。图中x轴从试验段左壁(沿流向左侧为左壁)指向右壁，y轴从试验段下壁面指向上壁面。试验段内共设置49个测点，其中红点表示试验段中心线处测点，测点的横向(x方向)、纵向(y方向)间距分别为375和250 mm，覆盖75%横纵向试验段区域。

图  5  温度格栅测点位置矩阵

Fig.  5  Measurement point position matrix of the temperature grid device

下载: 全尺寸图片 幻灯片

2.3   试验方法及数据处理

气流总温试验工况如表 3所示，表中给出了试验段名义总温(T_t)、试验段名义气流速度(v_TS)、喷嘴喷气压力(p_a)和喷雾耙入口干空气温度(T_a)的试验参数。为防止喷雾系统在低温环境下损坏，应保持热气在喷雾耙内持续流动以维持耙内温度，因此最小喷气气压选择0.03 MPa。

表  3  温度场校测工况

Table  3  Test conditions of airflow total temperature measurement

Parameter	Value
Tt/ ℃	－30, －20, －15, －10, －5, 0, 5
vTS/(m·s-1)	40, 80, 120, 140
pa/MPa	0.03, 0.40, 0.90
Ta/ ℃	80

下载: 导出CSV 

| 显示表格

试验时，当气流参数稳定后，同时采集试验段和风洞气流总温数据，各测点参数采样时间为60 s，采样频率1 Hz。

试验后，对试验段气流总温数据进行总温恢复系数修正，获得当地气流总温T_{t, local}：

式中，T_{t, probe}为总温探针测量总温(单位为℃)，R_r为探针总温恢复率。

气流总温空间分布和时间分布以试验段中心线处气流总温T_{tc, TS}为基准，采用气流总温空间偏差ΔT_{t, S}和时间偏差ΔT_{t, T}表征：

式中：下标mt表示时间平均，T_{t, TS}为试验段各测量点处气流总温。进而采用总温空间偏差和时间偏差的标准差(σ(ΔT_{t, S})、σ(ΔT_{t, T}))以及最大绝对值(|ΔT_{t, S}|_max和|ΔT_{t, T}|_max)评估试验段内气流总温的空间均匀性和试验段中心处的时间稳定性。

最后，利用试验段中心线处的气流总温T_{tc, TS}与风洞测量气流总温T_{t, WT}，获得试验段总温修正关系：

式中，K_Tt为总温修正函数。

3   试验结果与分析讨论

3.1   气流总温

3.1.1   气流总温的空间均匀性

为评估试验段气流总温的空间均匀性，图 6给出了试验段气流总温空间分布云图，包括40、80、120和140 m/s的气流速度条件下，气流总温0 ℃、－5 ℃和－20 ℃对应的典型试验结果(喷嘴气压均为0.03 MPa)。图中，红色虚线框表示模型区，范围为－750 mm≤x≤750 mm、－500 mm≤y≤500 mm。可以看出：随着试验段气流总温和气流速度的降低，试验段上下壁面附近逐渐出现总温偏差超过1 ℃的非均匀峰值区域，且非均匀峰值区域不断向模型区内部发展，特别是在140 m/s工况下，非均匀峰值区域已经发展到模型区的内部核心区域。由此可见，基于目前的制冷系统性能，降低试验段气流总温和提高试验段气流速度均会减弱试验段内热流场的空间均匀性。

图  6  试验段气流总温空间分布云图

Fig.  6  Spatial distribution of the total temperature in the test section

下载: 全尺寸图片 幻灯片

为进一步定量评估试验段热流场的空间均匀性，图 7给出了试验段模型区内的气流总温空间偏差标准差和最大绝对值。图中黑色虚线表示气流总温空间均匀性指标(见表 2)。可以看出：在各工况下，气流总温空间偏差标准差均满足空间均匀性指标要求，且在40~120 m/s、－20 ℃~0 ℃工况范围内标准差一致性较好，均小于0.6 ℃；随着试验段气流速度增大及气流总温降低，气流总温空间偏差最大绝对值整体增大，且逐渐超出标准范围，这与试验段上下壁面附近的非均匀峰值区域发展扩大相对应。在40、80和120 m/s条件下，非均匀峰值区主要集中于模型区边界处(如图 6所示)，当气流总温高于－20 ℃后，温度偏差最大绝对值均小于1.5 ℃；但在140 m/s条件下，当温度低于－5 ℃后，非均匀峰值区发展扩大至模型区核心区域，温度偏差峰值整体增大(普遍超过1.5 ℃)；此外，喷嘴干空气射流对气流总温空间均匀性并无显著影响。由此可见，模型区内气流总温空间均匀性在大部分范围内均满足指标要求，其中40 ~120 m/s、－20 ℃~0 ℃范围内的热流场空间均匀性最优，但是随着试验段气流速度增大以及气流总温降低，模型区内逐渐出现总温超标的峰值区域，尤其是在140 m/s条件下，模型区内普遍存在超标区域。

图  7  试验段模型区气流总温空间偏差标准差和最大绝对值

Fig.  7  The standard deviation and maximum absolute value of the spatial deviation of airflow total temperature in the model area of the test section

下载: 全尺寸图片 幻灯片

3.1.2   气流总温的时间稳定性

为评估试验段气流总温的时间稳定性，图 8给出了试验段中心线处的气流总温时间偏差变化曲线，包括喷嘴气压0.03、0.40和0.90 MPa状态下，气流速度40、80、120和140 m/s、气流总温0 ℃和－20 ℃对应的试验结果。可以看出：各工况下，气流总温时间偏差均在±0.5 ℃范围内；－20 ℃条件下的总温偏差脉动程度比0 ℃条件下更加剧烈，而试验段气流速度和喷嘴干空气射流对气流总温时间偏差并无显著影响。

图  8  试验段中心线处气流总温时间偏差变化曲线

Fig.  8  Profiles of the time deviation of the airflow total temperature in the centerline of the test section

下载: 全尺寸图片 幻灯片

进一步地，图 9给出了试验段中心线处的气流总温时间偏差标准差和最大绝对值。可以看出：试验段气流速度和气流总温以及喷嘴干空气射流对气流总温的时间稳定性并无显著影响，各工况下的气流总温时间偏差标准差和最大绝对值分别小于0.2 ℃和0.4 ℃。

图  9  试验段中心线处气流总温时间偏差标准差和最大绝对值

Fig.  9  The standard deviation and maximum absolute value of the time deviation airflow total temperature in the centerline of test section

下载: 全尺寸图片 幻灯片

3.1.3   气流总温修正关系和不确定度

图 10为试验段中心线处的气流总温修正关系和不确定度，包括喷嘴气压0.03、0.40和0.90 MPa对应的试验结果、拟合曲线和不确定度。从图 10(a)可以看出，风洞采集的气流总温与试验段中心线处的气流总温具有显著的线性关系，其拟合公式为：

图  10  试验段中心线处气流总温修正关系和不确定度

Fig.  10  The correction relationship and uncertainty of the airflow total temperature in the centerline of the test section

下载: 全尺寸图片 幻灯片

式中，K_T和B_T分别为斜率和截距函数，均为喷嘴气压的函数。

在图 10(b)中，纵轴为采用拟合公式计算的总温(T_{ta, TS})与试验段中心线处的气流总温之差，黑色虚线为不确定度指标。从图中可以看出，在0.03和0.40 MPa条件下，拟合温度偏差均满足标准要求，但在0.90 MPa条件下，存在个别超标的温度点，且这些点主要集中在40 m/s速度条件下，如图中的阴影区所示。考虑到主要试验工况下喷嘴气压均小于0.40 MPa，因此整体而言温度拟合公式的不确定度满足标准要求。

图 11给出了试验段中心线处气流总温与风洞测量总温之差。从图中可以看出：喷嘴干空气射流会提高气流温差，随着喷嘴气压增大，温差不断增大，但增大幅度会随着试验段气流速度的增大而减小，其中，在40 m/s、－30 ℃工况下最大温差超过2 ℃，而超过80 m/s后，温差普遍小于1 ℃(在140 m/s、－20 ℃工况下，温度场均匀性较差，则此时温差超过1 ℃)。这主要是因为：喷嘴向试验气流中注入的高温空气将自身携带的能量传递给试验气流，增大了试验段气流总温，且随着喷嘴气压增大，喷嘴热气流量增大，进而提高了试验段气流总温的增大幅度；但随着试验段气流速度的增大，试验气流流量不断增大，减弱了喷嘴干空气射流对试验气流的影响，导致气流总温的增大幅度减小。

图  11  试验段中心线处气流总温与风洞测量总温之差

Fig.  11  The difference between the total temperature in the centerline of the test section and the wind tunnel total temperature

下载: 全尺寸图片 幻灯片

3.2   结冰风洞热流场控制包线

根据结冰风洞主试验段热流场品质评估结果，给出了结冰风洞主试验段热流场控制包线，如图 12所示，其中横轴为试验段名义气流速度，纵轴为试验段名义气流总温，黑点为试验工况点。图中综合考虑试验段模型区内热流场的空间均匀性和时间稳定性，给出了控制包线(红线区域)，在该包线范围内，模型区内热流场品质基本满足SAE ARP 5905-2003的指标要求，但在模型区边界处仍存在温度偏差超标的峰值区。需要指出的是，在高风速(超过140 m/s)和低总温(小于－20 ℃)条件下，试验段模型区内存在明显的不均匀峰值区，该区域严重影响气流总温的空间均匀性，制约着结冰风洞的试验模拟范围，因此，需要对3 m×2 m结冰风洞制冷系统进行升级改造，提高高风速和低总温条件下的温度空间均匀性，扩展试验模拟范围。

图  12  3 m×2 m结冰风洞主试验段热流场控制包线

Fig.  12  The thermodynamic flowfield operating envelop of CARDC icing wind tunnel in the main test section

下载: 全尺寸图片 幻灯片

4   结论

本文建立了结冰风洞热流场符合性验证方法，开展了3 m×2 m结冰风洞主试验段热流场符合性验证试验，得到以下结论：

(1) 基于目前制冷系统性能，降低试验段气流总温和提高试验段气流速度会减弱试验段内热流场的空间均匀性，但对时间稳定性并无显著影响。

(2) 喷嘴干空气射流对试验段模型区内气流总温的空间均匀性和时间稳定性并无显著影响，但会提高试验段气流总温，尤其对40 m/s条件下的热流场的影响最为显著。

(3) 在主要试验条件下，结冰风洞热流场品质基本满足SAE ARP 5905-2003指标要求。