Cloud calibration method of 3m×2m icing wind tunnel
-
摘要: 为综合检验结冰风洞的性能,并为飞机适航取证中结冰合格审定工作提供测试依据,依照SAE ARP5905的相关要求,根据冰积聚原理、相位多普勒干涉原理和热线原理,采用格栅、PDI-FPDR、冰刀、LWC-200在3m×2m结冰风洞主试验段完成了云雾参数校测工作,总结了测试方法与流程以提高校测效率。代表性云雾参数校测结果表明,云雾场均匀区能够覆盖试验段横截面积的60%以上,云雾容积平均直径(MVD)稳定性优于±10%,液态水含量(LWC)稳定性优于±20%,云雾参数满足ARP 5905相关指标要求,能够为3m×2m结冰风洞的标检提供数据支撑。Abstract: The cloud calibration has been carried out under SAE ARP5905 standard for 3m×2m icing wind tunnel, results of which can be used for certificating icing wind tunnel performance and supplying data proof for obtaining icing flight certification. Based on the ice accretion mechanism, the phase doppler interference principle, and the hot-wire principle, the calibration work was completed in the main test section by using the icing grid, PDI-FPDR, the icing blade, and LWC-200. Methods and procedures are summarized to increase the efficiency of calibration. Typical cloud calibration results indicate that the icing cloud uniformity is good, 60% of the main test section could remain target LWC within ±20%. The MVD maximum deviation is smaller than ±10%. The LWC stability is greater than ±20%. All the cloud parameters could meet the requirements of ARP5905, which could be used for supplying data proof of 3m×2m icing wind tunnel certification.
-
Keywords:
- icing wind tunnel /
- icing test /
- airworthiness /
- cloud parameter /
- calibration
-
0 引言
结冰会破坏飞机气动外形,影响飞机的气动性能与操稳特性,是影响飞行安全的重要隐患之一。因此,世界各国的民用航空适航规章中都明确规定,飞机适航取证过程中必须进行结冰的适航合格审定工作。
结冰风洞作为地面模拟设备,其试验条件连续可控、不受外界气候条件限制,较易实现FAR/CCAR 25部附录C的结冰云雾包线模拟,且试验成本低、周期短、可靠性高,常被选择作为结冰适航审定的主要符合性验证手段之一。
意大利CIRA-IWT和美国NASA-IRT是国际上最具代表性的2座结冰风洞,先后参与完成了波音、空客等多个型号飞机的结冰合格审定工作,在国际结冰研究机构中占有重要地位。这些结冰风洞在正式承接型号试验前均要经过相关标准的验证与标检,目前国际上大多数结冰风洞依照SAE颁布的ARP5905[1]进行校测,以综合检验结冰风洞的相关性能指标是否满足结冰和防除冰试验要求。
3m×2m结冰风洞建成后,参照FAR/CCAR 25部附录C的结冰云雾包线范围CanKaoWenXian_3 ,按照SAE ARP5905相关标准要求完成了云雾参数校测工作。本文详细介绍了3m×2m结冰风洞在主试验段完成的云雾参数校测流程与方法,并给出代表性结果。
我国大飞机、宽体客机、通用飞机等民用飞机未来走向国际市场,必须经过适航规章的合格审定。因此,本文旨在通过介绍3m×2m结冰风洞校测与验证的流程、方法和结果,为我国航空飞行器研制过程中结冰合格审定的取信工作提供技术支撑。
1 试验设备
3m×2m结冰风洞是一座拥有3个可更换试验段的闭口、回流式风洞,能够实现低温、低压和高湿的模拟结冰条件。
结冰风洞的辅助系统包括:制冷系统、喷雾系统、高度模拟系统、加湿系统、发动机进气模拟系统和防/除冰系统等,分别可实现温度模拟、云雾模拟、压力模拟、湿度控制、发动机功率状态模拟和模型热气供气等功能。风洞的基本性能参数如表 1所示。
表 1 结冰风洞基本性能参数Table 1 Main performance parameters of icing wind tunnelContent Main test
sectionSecond test
sectionHigh speed
test sectionSize 3m×2m×6.5m 4.8m×3.2m×9m 2m×1.5m×4.5m Speed 21~210m/s 8~78m/s 26~256m/s Temperature Normal~-40℃ Humidity 70% ~100% Altitude 0~20000m Cloud MVD: 10~300μm
LWC:0.2~3g/m3
Uniformity: 60% of section喷雾系统由1000个水气内混合雾化喷嘴(Spray 98818)组成,可通过改变喷嘴的供水、供气压力生成具有特定液滴尺寸和液态水含量的结冰云雾,喷嘴可独立控制开关。
本文的校测工作在主试验段开展。
2 云雾参数校测
云雾参数是衡量结冰云雾场性能的重要参数,是影响结冰试验结果精准度的关键,其主要包括均匀度、容积平均直径(MVD,Median Volume Diameter)和液态水含量(LWC,Liquid Water Content)等。3m×2m结冰风洞建成以来,依据SAE(美国机动车工程师学会)、FAA(美国联邦航空管理局)、中国民用航空局等颁布的若干标准[1-6],并主要参考ARP5905的《结冰风洞校准与验收》,开展了主试验段的云雾场校测。此外,校测流程与方法还参照了美国NASA-IRT结冰风洞和意大利CIRA-IWT结冰风洞的校测报告[7-10]。
云雾参数中的3个重要参数可通过调节水压、气压、喷嘴开度实现不同条件的模拟。其中水压通过水处理系统调节,气压通过供气系统调节,喷嘴开度通过喷雾控制系统开关喷嘴水路的电磁阀进行调节,通常采用全开、1/2开和1/4开等3个开度。
2.1 云雾均匀度校测
在进行MVD和LWC校测前,首先需要进行云雾均匀度校测,用以确定喷雾系统喷嘴的布局与开度。目前国内外大多数结冰风洞均以格栅校测为主。
3m×2m结冰风洞主试验段云雾均匀度校测采用1955mm×2705mm的格栅,格栅与上下洞壁相连,其面积覆盖了试验段中心88%的横截面积。格栅网格为150mm×150mm,格栅单元深60mm,迎风面宽度为5mm。
校测中,结冰试验按霜冰条件进行,试验静温低于-18℃,MVD为20μm(喷嘴地面测试时初步获得初始的水压、气压匹配),格栅前缘目标厚度6.4mm左右。图 3为典型条件下均匀度格栅前缘结冰情况。
获得格栅结冰后,采用经过预冷的数显千分尺测量格栅迎风面的结冰厚度,依据编号,记录所有格栅横、竖栅条中心位置的结冰厚度。
结冰风洞的云雾均匀度表征的是试验段液态水含量的均匀程度,通常用相对量表示。根据ARP5905,液态水含量均匀度可通过格栅上结冰的相对厚度表示,其关系式如下:
(1) 式中:LWC(x, y)为格栅上某位置的液态水含量,g/m3;LWCC为格栅上中心位置(交叉点)的液态水含量,g/m3;τ(x, y)为格栅上某位置的冰厚度,mm;τgrid为格栅厚度,mm;τC为格栅上中心位置的冰厚度,mm;τiceR为相对厚度。
格栅上的测量位置一般取水平和竖直栅条的中心位置,如图 4所示。格栅中心位置的冰厚度计算公式如下:
(2) 记录不同条件下均匀度格栅上冰的相对厚度,绘制二维厚度云图(等高线图),用以评估试验段的云雾场均匀度。
图 5和6给出3m×2m结冰风洞主试验段典型条件下的云雾场均匀度校测结果。结果表明,试验段截面内80%以上区域液态水含量偏差在±20%内,云雾场的均匀性良好,满足ARP5905中关于LWC的空间分布要求。靠近上洞壁中部位置LWC整体偏小,而两侧一致性较好,根据格栅支撑结构位置情况,分析认为格栅上部支撑的3块矩形连接装置是造成该现象的主要原因,平面阻挡带来了前缘气流扰动,改变了水滴的撞击特性,使得格栅前缘的水收集率降低,表征为结冰厚度减小。该连接装置需要进一步优化,减少结构迎风面积。
云雾场均匀度校测是一项十分复杂且工作量巨大的测试工作,采用良好的校测流程会大大提高校测效率。根据已完成的校测工作建议结冰风洞校测流程如下:
(1) 根据喷嘴的地面测试结果,确定喷嘴的良好工作区间(即水气压组合条件),在良好的工作条件下,喷嘴之间的个体差异通常会随着水气压的增大而降低;
(2) 建立不同喷嘴开度与LWC的对应关系图,每个开度的LWC下限通常为有效的喷嘴水气压组合状态,即最小的水气压组合。图 7为喷嘴开度与LWC的对应关系的示意图;
(3) 减小喷嘴开度时,以喷嘴之间近似等间距、距离最短为原则进行配置;
(4) 影响均匀度的因素主要包括:喷嘴个体差异、水气压组合、喷嘴开度、风速等。在不更换喷嘴或改变喷嘴开闭的情况下,均匀度与影响因素之间的关系如图 8所示。在低LWC范围内,建议优先选择更高的水气压组合条件、更低的开度;
(5) 在保证均匀度良好的前提下,喷嘴开度与试验段LWC近似正比关系,风速与LWC近似反比关系,校测前可根据此原则进行估算。
2.2 MVD校测
云雾中悬浮的液态小水滴大小不一,其变化范围可由几微米至几百微米,通常采用单一参数描述云雾全部水滴的总体特性,常采用容积平均直径MVD代表[11-12]。目前各国结冰风洞对于MVD校测采用的仪器有所不同,主要设备包括前向散射分光仪(FSSP)、光学阵列测量仪(OAP)、马尔文测量仪(Malvern Spraytec)、相位多普勒粒子分析仪(PDPA)和机载式相位多普勒干涉仪(PDI-FPDR)等。
3m×2m结冰风洞MVD校测采用Malvern Spraytec和PDI-FPDR进行多仪器匹配测量。Malvern Spraytec主要进行喷嘴的地面测试,测量范围为0.5~2000μm,测量精度为±1%(Dv50标准粒子)。PDI-FPDR主要用于洞内校测,测量范围为0.5~1000μm,测量精度为0.5μm。
PDI-FPDR是在相位多普勒干涉仪(PDI)基础上发展的新一代相位多普勒颗粒分析系统,其工作原理是采用激光器作为光源,将颗粒看作一个微小的透镜,测量颗粒对平行入射激光的散射光变化,以其不同空间接收位置散射光的相位变化来反映该微小透镜的焦距大小,即颗粒的粒径大小;通过频率变化反映其运动速度的大小[11]。
MVD洞内校测是将PDI-FPDR测量仪置于风洞中心,给定固定的水压(pw)、气压(pa)匹配工况,连续采集云雾水滴直径30s以上,待测试结果稳定后记录MVD、水压和气压。改变不同水压、气压匹配状态,重复测量过程,建立MVD、pw和pa数据库。
(1) 云雾MVD模拟方法
基于校测过程中所建立的MVD、pw和pa数据库,通过数值解析法分析MVD与pw、pa的相关性,采用最小二乘误差分析法建立MVD与pw、pa的关系函数MVD=f(pw, pa)。该关系函数作为结冰云雾的控制依据,可以根据试验前所需的MVD和LWC快速查询系统控制所需的pw、pa。
(2) 误差分析
绘制MVD=f(pw, pa)函数曲线,并将实际校测值以离散点形式绘制于同一曲线图中,分析校测值与理论分析结果的偏差是否满足±10%的偏差要求(见ARP5905 P16)。此外,校测标准中备注指出,当MVD小于30μm时,以±3μm的不确定度替代10%偏差要求。
(3) 重复性校测
开展相同模拟条件下的重复性校测,检验是否满足±10%的偏差要求。
(4) 针对性校测
针对不同试验的模拟条件需求,试验前可根据试验条件进行查表与计算获得喷雾系统所需的pw、pa,再抽检典型状态进行风洞内实际校测,以确保云雾参数控制的准确度。
图 10给出了3m×2m结冰风洞主试验段典型MVD校测结果。其中(a)为试验测量结果与拟合曲线,(b)为试验测量结果与数值拟合结果的差异性分析。结果表明,MVD随水压增大而增加,随气压增大而减小;在气压大于某一临界值时,将无法形成稳定的喷雾条件;对于同一水气压组合条件下的MVD重复性偏差小于10%,在30μm以下偏差值小于3μm;水滴直径模拟范围能够覆盖FAR 25部附录C规定的15~50μm的范围,且相对于数值拟合分析结果,实际测试结果在±10%的偏差范围内。
2.3 LWC校测
LWC代表单位体积内的总含水量。目前各国结冰风洞采用的测试仪器主要包括冰刀、旋转圆柱、LWC热线测量仪等[11-12]。
3m×2m结冰风洞LWC校测主要有LWC-200热线测量仪(见图 11)和冰刀装置(见图 12)。LWC-200热线测量仪因能够快速获得LWC测量值,主要用于提高LWC的校测效率。测量范围为0~3.0g/m3,响应频率25Hz。校测过程中,测量仪安装在试验段中心,即可以根据仪器测量值建立起LWC与其他影响参数之间的关系。
冰刀法是国内外常用LWC校测方法,但总体效率较低,主要用于校测结果间的相互对比验证。校测过程中,将冰刀支撑于风洞试验段中心,在指定试验条件下运行风洞并开启喷雾系统,喷雾稳定后,将冰刀防护罩打开,待冰生长一定时间后,将防护罩关闭。使用预冷的数显游标卡尺测量冰刀工作面上的冰厚度(不少于3个测点位置),记录平均值,并按下式计算液态水含量[14]。
(3) 式中:ρice为冰的密度,kg/m3;τ为冰刀上平均结冰厚度,m;Eb为收集系数,由CFD计算获得;V为模型区中心风速,m/s;t为结冰时间,s。
改变不同水压、气压、水滴直径、风速(V)、模拟高度(H)、喷嘴开度(N_D)等控制条件,重复测量过程,建立LWC与影响参数的数据库。
(1) 云雾LWC模拟方法
通过多步数值分析法建立LWC与影响参数的关系函数如下:
(4) 该关系函数同样作为结冰云雾的控制依据,可以根据试验前所需的MVD和LWC查询或求解系统控制所需的pw、pa。
(2) 误差分析
设置其他影响参数为定值,绘制LWC=f(pw, MVD)函数曲线,并将实际校测值以离散点形式绘制于同一曲线图中,分析校测值与理论分析结果的偏差是否满足±20%的偏差要求(见ARP5905 P16)。
(3) 重复性校测
开展相同模拟条件下的重复性校测,检验是否满足±20%的偏差要求。
(4) 针对性校测
针对不同试验的模拟条件需求,首先根据试验条件进行查表与计算获得喷雾系统所需的pw、pa,再进行试验条件抽样校测,确保云雾参数控制的准确度。
图 13描述了当前结冰风洞部分已校测的数据点与FAR 25部附录C规定的结冰云雾包线间的关系,这表明风洞已经完成了一定范围的包线内典型状态校测,部分状态点的模拟能力已超出结冰云雾包线范围。理论上,在风速不受限定的情况下,3m×2m结冰风洞能够模拟大部分的结冰云雾包线,但仍需关注高风速、低液态水含量下的均匀性。图 14给出了3000m高度下的校测结果,其中(a)为LWC校测试验结果与数值拟合曲线,(b)为试验测量结果与数值拟合结果的不确定度分析。结果表明,在保持MVD不变的情况下,LWC随着水压的增大而增加;对于同一水气压组合条件,LWC的重复性偏差小于20%,且部分工况下优于±10%;LWC校测结果与数值拟合结果能够满足±20%的指标要求。
2.4 小结
3m×2m结冰风洞云雾参数校测主要依据ARP5905校测标准要求,以实际校测结果数据为输入条件,分析建立喷雾系统控制策略与方法,并进行误差分析与可靠性分析。具体流程如图 15所示。
结冰风洞实际能力能够覆盖大部分FAR/CCAR 25部附录C的结冰云雾包线范围,但校测前应统筹规划好需要测试的离散结冰状态点,以提高校测工作的效率,确保云雾参数校测的数据库完整性,并需在风洞运行过程中定期开展校测与检验。
3 结论与建议
3m×2m结冰风洞已经完成了主试验段的云雾场校测,云雾场均匀度能够覆盖试验段横截面积的60%以上,云雾容积平均直径(MVD)稳定性优于±10%,液态水含量(LWC)稳定性优于±20%,云雾参数满足ARP 5905相关指标要求。然而,结冰风洞的云雾参数校测是一个复杂且长期的工作,ARP5905中指出结冰风洞每5年需要进行一次全面校测,期间还要进行中间校测和检验性校测,其中检验性校测需每6个月进行一次,主要检测云雾参数是否有大的变动。因此,结冰风洞需要通过持续不断地开展定期与不定期校测工作,以确保试验参数与试验结果的重复性和稳定性。
由于国内大多数MVD和LWC测试设备均为引进国外成熟产品,在设备本身的标定与检测方面还存在一定困难和局限性,建议我国各类结冰风洞在建立正规的测试仪器标定流程前,先进行单仪器的重复性对比与多仪器间的互对比研究,以确保测试仪器的精准度。
-
表 1 结冰风洞基本性能参数
Table 1 Main performance parameters of icing wind tunnel
Content Main test
sectionSecond test
sectionHigh speed
test sectionSize 3m×2m×6.5m 4.8m×3.2m×9m 2m×1.5m×4.5m Speed 21~210m/s 8~78m/s 26~256m/s Temperature Normal~-40℃ Humidity 70% ~100% Altitude 0~20000m Cloud MVD: 10~300μm
LWC:0.2~3g/m3
Uniformity: 60% of section -
[1] SAE Aerospace. SAE ARP5905 Calibration and acceptance of icing wind tunnels[S]. SAE, 2003. http://www.sae.org/standards/content/arp5905/
[2] Federal Aviation Administration. FAR Part25 Airworthiness standards: Transport category airplanes[S]. Federal Aviation Agency, 2008.
[3] 中国民用航空局. CCAR-25-R4中国民用航空规章第25部: 运输类飞机适航标准[S]. 北京: 中国民用航空局, 2011. [4] SAE Aerospace. SAE AIR4906 Droplet sizing instrumentation used in icing facilities[S]. SAE, 1995. https://www.sae.org/standards/content/air4906/
[5] SAE. SAE AIR5320 Summary of icing simulation test facilities[S]. SAE, 1999. https://www.sae.org/standards/content/air5320/
[6] SAE Aerospace. SAE AIR5624 Aircraft inflight icing terminology[S]. SAE, 2008. http://articles.sae.org/2917/
[7] Sten L E, Ide R F. NASA glenn icing research tunnel: 2014 and 2015 cloud calibration procedures and results[R]. NASA/TM-2015-218758, 2015. https://www.sciencedirect.com/science/article/pii/S0376042118300058
[8] Van Zante J F, Ide R F. NASA glenn icing research tunnel: 2014 cloud calibration procedure and results[R]. NASA/TM-2014-218392, 2014. DOI: 10.2514/6.2012-2933
[9] Bellucci M, Esposito B M. Calibration of the CIRA IWT in the low speed configuration[R]. AIAA-2007-1092, 2007. http://edu.wanfangdata.com.cn/Periodical/Detail/kjdb201325007
[10] Ragni A, Esposito B. Calibration of the CIRA IWT in the high speed configuration[R]. AIAA-2005-471, 2005. http://www.academia.edu/28736885/Calibration_of_the_CIRA_IWT_in_the_High_Speed_Configuration
[11] 林贵平, 卜雪琴, 申晓斌, 等.飞机结冰与防冰技术[M].北京:北京航空航天大学出版社, 2016. [12] 朱春玲, 朱程香.飞机结冰及其防护[M].北京:科学出版社, 2016. [13] 程尧, 张平涛, 郭向东, 等.机载式相位多普勒干涉仪在结冰风洞的应用[J].兵工自动化, 2017, 36(9):55-57. http://www.cqvip.com/QK/90272A/201502/664576270.html Cheng Y, Zhang P T, Guo X D, et al. Application of PDI-FPDR in icing wind tunnel[J]. Ordnance Industry Automation, 2017, 36(9):55-57. http://www.cqvip.com/QK/90272A/201502/664576270.html
[14] 赖庆仁, 郭龙, 李明, 等.结冰风洞液态水含量测量装置设计与实现[J].空气动力学学报. 2016, 24(6):750-755. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_kjdb200921018 Lai Q R, Guo L, Li M, et al. Design and implementation of the device for liquid water content measurement in icing wind tunnel[J]. Acta Aerodynamica Sinica, 2016, 24(6):750-755. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_kjdb200921018
-
期刊类型引用(19)
1. 高郭池,张波,全敬泽,尹崇,丁丽,姜裕标. 正常类飞机自然结冰试飞适航审定技术. 航空学报. 2024(01): 195-216 . 百度学术
2. 陈方备,王畅,戴铮,崔燚,吴健. 空气雾化喷嘴喷雾特性实验研究. 工程热物理学报. 2024(03): 873-878 . 百度学术
3. 赵照,熊建军,冉林,易贤. 大型结冰风洞热气供气防除冰试验技术. 航空动力学报. 2024(05): 16-22 . 百度学术
4. 赵照,王梓旭,熊建军,冉林,易贤. 大型结冰风洞双供水供气喷雾控制系统研制. 测控技术. 2024(07): 48-54 . 百度学术
5. 魏龙涛,刘森云,王桥,郭奇灵. 0.75m×0.50m结冰风洞气动——热流场品质评估. 航空工程进展. 2024(04): 162-170+2 . 百度学术
6. 魏龙涛,左承林,郭向东,刘森云,郭奇灵. 三维冰形激光测量技术在结冰风洞中的应用. 航空动力学报. 2024(11): 246-253 . 百度学术
7. 赵照,易贤,冉林,熊建军. 多路热气流量控制及其在防除冰试验中的应用. 测控技术. 2023(10): 18-23+37 . 百度学术
8. 赵照,熊建军,冉林,郭向东,李自雨. 基于ZigBee技术的结冰风洞温度无线测量系统设计. 电子设计工程. 2022(05): 51-55+60 . 百度学术
9. 熊建军,刘锡,冉林,赵照. 基于控制律的电加热防除冰系统设计与验证. 测控技术. 2021(02): 130-134+139 . 百度学术
10. 赵照,熊建军,张平涛,冉林,李自雨. 基于分布式光纤测温的结冰风洞喷雾耙温度场测量. 计算机测量与控制. 2021(03): 14-17+22 . 百度学术
11. 郭向东,柳庆林,赖庆仁,杨升科,赵照. 大型结冰风洞气流场适航符合性验证. 空气动力学学报. 2021(02): 184-195 . 百度学术
12. 赵照,熊建军,冉林,何苗. 基于动态神经网络的喷雾耙供气温度数据重构. 仪表技术与传感器. 2021(04): 116-121 . 百度学术
13. 陈舒越,郭向东,王梓旭,刘森云,吴迎春. 结冰风洞过冷大水滴粒径测量初步研究. 实验流体力学. 2021(03): 22-29 . 本站查看
14. Zhiqiang GUO,Mei ZHENG,Qian YANG,Xiaofeng GUO,Wei DONG. Effects of flow parameters on thermal performance of an inner-liner anti-icing system with jets impingement heat transfer. Chinese Journal of Aeronautics. 2021(09): 119-132 . 必应学术
15. 王先炜,樊晓锋,林森什,陈彪. 结冰对直升机进气系统压力损失影响的试验研究及分析. 直升机技术. 2021(04): 62-66 . 百度学术
16. 战培国. 结冰云小水滴粒径测量设备综述. 测控技术. 2020(06): 1-7 . 百度学术
17. 郭向东,张平涛,赵献礼,杨升科,林伟. 大型结冰风洞热流场符合性验证. 实验流体力学. 2020(05): 79-88 . 本站查看
18. 郭向东,张平涛,赵照,赖庆仁,郭龙. 大型结冰风洞云雾场适航应用符合性验证. 航空学报. 2020(10): 205-219 . 百度学术
19. 倪章松,刘森云,王桥,王梓旭,郭龙. 3m×2m结冰风洞试验技术研究进展. 实验流体力学. 2019(06): 46-53 . 本站查看
其他类型引用(7)