民用飞机座舱温度场仿真分析与研究
2020-07-20龙正伟汪光文
杨 智 龙正伟 汪光文
民用飞机座舱温度场仿真分析与研究
杨 智1龙正伟2汪光文1
(1.上海飞机设计研究院 上海 201210;2.天津大学环境科学与工程学院 天津 300072)
极端气候条件下民用飞机座舱温度场验证是民用飞机设计的难点,首先基于某型民用飞机座舱几何参数,建立CFD计算模型;其次采用该模型计算常温环境下的座舱温度场,并与试验对比以验证模型;最后再采用验证过的模型计算极端气候条件下的座舱温度场,以验证某型民用飞机极端气候条件下座舱温度场是否满足设计要求。
座舱;温度场;CFD
0 引言
舒适性是民用飞机的关键要求之一,座舱温度是旅客感受民用飞机舒适性的直接因素。同时民用飞机的飞行高度一般在9000~12000m左右,舱外环境温度在可低至-65℃左右。在如此极端的环境下,维持座舱温度场在合适的范围内是民用飞机的基本要求之一。为此,民用飞机的设计要求中规定了座舱温度场的相关设计要求,这些设计要求的验证是民用飞机研制中非常重要的环节。极端气候条件下如地面环境温度为-40℃时,座舱温度设计要求如何验证是民用飞机研制中的难题,一是地面环境温度为-40℃的气候条件较难寻找,且影响民用飞机研制进度;二是若采用试验验证,花费巨大且耗时耗力;三是座舱人员对座舱温度场影响较大,而民用飞机试飞时座舱人员很难满载。因此若可通过计算验证极端气候条件下座舱人员满载时的座舱温度设计要求,将对民用飞机的研制有积极地推动和指导作用。
本文针对极端气候条件下民用飞机座舱温度场设计要求验证问题,基于某型民用飞机座舱真实的几何参数,建立CFD计算模型,采用该模型计算常温环境下的座舱温度场,并与试验对比验证模型后,再用于极端气候条件下的座舱温度场计算,以验证某型民用飞机极端气候条件下座舱温度场是否满足设计要求。
1 民用飞机座舱温度场设计要求
SAE ARP 85[1]中定义了民用飞机座舱温度需要满足的要求,如表1所示。
表1 SAE ARP 85座舱温度要求
SHRAE161-2007[2]也规定了飞机座舱温度设计要求,如表2所示。
表2 ASHRAE 161-2007座舱温度要求
2 座舱温度场计算模型
2.1 物理模型及网格划分
根据某型飞机座舱真实三维几何模型,建立相应的CFD数值模拟模型。座舱天花板附近及行李架附近左右两侧各设置一送风口,即总共有四个送风口。回风口设置于地板附近两侧侧壁面上,几何模型如图1所示。
图1 座舱几何模型示意图
网格生成工具采用Gambit软件,网格划分形式采用非结构网格,生成的网格数量为1000万。为保证计算结果的准确性及计算资源的合理利用,网格划分时分别对座舱内送风口、回风口及壁面附近进行了加密处理,网格示意图如图2所示。
图2 座舱网格示意图
2.2 数学模型
采用RANS湍流模型中稳态RNG-模型对座舱内气流分布进行模拟预测。对于稳态RNG-模型,压力与动量方程的耦合方式采用SIMPLE算法;空间离散方法中,压力离散采用PRESTO算法,动量的差分格式采用有限中心差分准则,其它变量的压力差分格式都采用二阶向前差分准则。当能量的残差低于10-7,其他变量的残差低于10-4时,认为所计算的流场已经达到收敛。
2.3 模型验证
2.3.1 常温天座舱温度场测试
按照ASHRAE161-2007中的要求,分别在某型飞机座舱前部、中部以及后部区域布置三维风速仪,测量了座舱无乘客时每个区域距离地板100mm、610mm、1090mm处的温度和速度;并使用热成像仪测量了座舱天花板、地板、内饰壁面等模型计算中采用的边界参数。试中采用的超声波风速仪如图3所示。
图3 超声波风速仪示意图
2.3.2 常温天座舱温度场计算
采用模型计算了常温天时座舱无乘客时流场与温度场,计算时采用的边界条件如表3所示。
表3 常温天工况座舱内边界条件
座舱流场与温度场计算结果示意如图4~6所示。
图4 座舱内流场分布
图5 座舱内走廊中央沿机身方向流场分布
图6 座舱内走廊中央沿机身方向温度分布
2.3.3 座舱温度场测试与计算结果对比
常温天座舱温度场测试与试验测量结果对比如图7所示。
图7 座舱温度场测试与计算结果对比
根据模拟结果与实测结果的对比情况可以看出,无论是从数值的大小与变化的趋势上看,两种计算工况下的数值计算结果与实验测量数据都能有较好的吻合,证明了现有数值计算方法对座舱温度场进行预测的正确性,即计算模型可以用来验证极端气候条件冷天工况下的座舱温度场。
3 极端气候冷天座舱温度场计算
采用上述计算模型对某型飞机冷天座舱人员满载时座舱流场与温度场进行了计算,以验证极端气候条件冷天(地面环境温度为-40℃)工况下、座舱人员满载时的座舱温度场要求。计算时采用的边界条件设置如表4所示。
表4 冷天地面工况座舱内边界
座舱流场与温度场计算结果示意如图8~13所示。
图8 座舱内流场分布
图9 座舱内走廊中央沿机身方向流场分布
图10 座舱内走廊中央沿机身方向温度场场分布
图11 座舱温度场分布示意图(头靠位置)
图12 座舱温度场分布示意图(腰部位置)
图13 座舱温度场分布示意图(脚部位置)
座舱内空间温度变化计算结果如表5所示。
表5 空间温度变化计算结果
由表5可知,某型飞机冷天座舱空间温度变化满足表2中ASHRAE 161-2007座舱温度要求。
4 结论
基于某型民用飞机座舱真实的几何参数,建立了CFD计算模型。采用计算模型计算了常温天地面环境下的座舱温度,并与试验对比验证了计算模型后;采用验证后的计算模型计算了极端气候条件冷天下的座舱温度,解决了极端气候条件下座舱温度验证问题,并减少试验测试,节约了飞机研制成本。
[1] Society of Automotive Engineers. SAE ARP 85 Air Conditioning Systems For Subsonic Airplanes[S].
[2] American Society of Heating, Refrigerating and Air-Conditioning Engineers. ANSI/ASHRAE Standard 161-2007 Air Quality within Commercial Aircraft[S].
[3] Pang L, Qin Y, Liu D, et al. Thermal comfort assessment in civil aircraft cabins[J]. Chinese Journal of Aeronautics, 2014,24(2):180-186.
[4] Park S, Hellwig R, Grun G, et al. Local and overall thermal comfort in an aircraft cabin and their interrelations[J]. Building and Environment, 2011,46: 1056-1064.
[5] Cui W, Ouyang Q, Zhu Y. Field study of thermal environment spatial distribution and passenger local thermal comfort in aircraft cabin[J]. Building and Environment, 2014,80:213-220.
[6] 林国华,袁修干,杨燕生.人机环境系统中的CFD研究[J].航空学报,1999,20(增):21-23.
[7] 林国华,杨燕生,袁修干.座舱环控系统气流组织的数值研究[J].应用基础与工程科学学报,1998,6(3):302-307.
[8] 林国华.空调座舱流体流动和传热的数值研究[M].北京:北京航空航天大学,1998.
Simulation Research and Analysis of Cabin Temperature Filed for Civil Aircraft
Yang Zhi1Long Zhengwei2Wang Guangwen1
( 1.Shanghai Aircraft Design and Research Institue, Shanghai,201210; 2.School of Environmental Science and Engineering, Tianjin University, Tianjin, 300072 )
The cabin temperature filed verification is difficult at extreme cold and hot day in civil aircraft design. In this paper, the CFD calculation model was built duet to the cabin geometry parameter. Then the cabin temperature filed at stand day was calculated to use the CFD model, and at the same time the calculation result was compared with the test to validated the model. With the validated model, the cabin temperature filed at extreme cold and hot day was calculated to validate the cabin temperature filed requirement is met or not.
Cabin; temperature filed; CFD
V245.3
A
1671-6612(2020)03-312-04
国家重点基础研究发展计划(2012CB720100)
杨 智(1984.10-),男,硕士研究生,高级工程师,E-mail:yangzh@comac.cc
2019-05-25