结垢对飞机空调换热器影响的试验研究

为研究换热器结垢对于飞机空调性能的影响,并得出规律及结论,搭建了换热器性能试验台,包括高压引气单元、冲压空气单元、换热器连接单元、测控系统等。该试验台试验了飞机制冷系统初级散热过程,通过干净及结垢换热器的试验对比结果,从换热及阻力特性两方面得出换热器结垢影响的相应规律:换热器冷边流量越大,结垢对于换热器换热特性影响越小;换热器冷边流量越大,结垢对于换热器阻力特性影响越大。与理论分析相符,表明了试验的有效性。     

电脉冲除冰系统电磁脉冲力仿真分析

为了获得更精确的电脉冲除冰(EIDI)系统中电磁力在飞机蒙皮上的分布及其随时间变化的情况,针对系统的工作原理与工作过程进行了仿真研究。计算分析了EIDI系统电路中电流随时间变化的关系;通过麦克斯韦方程组的求解完成了EIDI系统电磁场的仿真;计算分析了电磁力随时间变化的规律,获得了蒙皮上的磁感应强度分布、涡电流密度分布和瞬态电磁力密度分布。考虑线圈工作时蒙皮对其电感的影响,建立了更接近真实电路的激励加载方法,得到了准确的电磁力变化趋势;考虑趋肤效应对蒙皮厚度方向上涡电流密度、电磁力密度分布的影响,用有限元分析方法使蒙皮上电磁力分布的计算更为全面合理。      

三维发动机进气道水滴撞击特性分析

用带数值扩散项的欧拉法对某型三维发动机进气道水滴撞击特性进行了研究.实现了可压流场中三维复杂形体的局部水滴收集系数的计算.数值扩散项的引入解决了欧拉法计算产生奇异结果发散的问题.通过对不同进气道流量的算例结果的比较,发现进气道流量对进气道唇口附近的水滴撞击区和撞击量都有影响.进气量大,最大局部水滴收集系数也大.通过对有无整流罩的计算结果进行比较,发现此型发动机整流罩对水滴撞击特性的影响很小.      

二维翼型大尺寸过冷水滴撞击特性及冰形分析

基于现有大尺寸过冷水滴(supercooled large droplet, SLD)动力学特性,分析水滴变形对阻力的影响.并根据几种典型的反弹/飞溅模型,分析了SLD的阻力变化、反弹、飞溅等对水滴撞击特性的影响,采用软件FENSAP-ICE的飞溅模型和LEWICE 2.0的反弹模型研究了反弹及飞溅现象对冰形的影响.计算结果显示:水滴撞击前的破碎现象对水滴尺寸分布有较大的影响,进行撞击特性以及冰形计算的时候需进行考虑;SLD破碎、飞溅、反弹降低了局部水收集系数、减小了水滴撞击范围;飞溅现象主要发生在机翼前缘附近区域,反弹主要在撞击边缘区域;水滴直径增加,飞溅现象逐渐减弱,但边缘位置的反弹现象一直很明显.SLD变形带来的阻力影响对冰形及结冰区域影响很小;与未考虑飞溅及反弹现象得到的冰形比较,考虑飞溅及反弹得到冰形前缘区域形状变化不大,但是整体结冰区域减小.      

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.      

Roe格式在水滴欧拉方程数值求解中的应用

目前普遍采用欧拉法求解水滴运动控制方程以模拟结冰表面水滴撞击特性。水滴欧拉方程近似没有压力项的空气欧拉方程,但已应用于空气欧拉方程的空间离散格式多数不能直接应用于水滴欧拉方程,有研究者将经典Roe格式应用于水滴欧拉方程时出现了奇点问题。通过对水滴欧拉方程进行形式改造,基于中心格式验证改造方法的有效性,并应用Roe格式结合线性重构与限制器函数实现对水滴欧拉方程二阶精度的求解。将NACA23012翼型计算结果与FENSAP-ICE软件进行对比,两者水滴收集系数分布基本一致。     

基于特征正交分解法的翼型结冰冰形快速预测

为加快翼型结冰冰形的计算速度,提出了一套基于特征正交分解(POD)法的结冰冰形快速预测算法.通过计算流体动力学(CFD)数值模拟计算得到的冰形结果作为样本,以结冰温度单参数变化时为例,详细介绍了POD法预测结冰冰形的实现步骤.考虑结冰温度、结冰时间以及液态水含量的影响,完成了单参数、两参数与三参数变化时的POD结冰冰形快速预测.通过算例结果发现:完成POD预测计算只需几秒钟且POD法与CFD法得到的冰形吻合较好,仅在冰角附近有较小差别,表明POD方法能够快速、准确地得到翼型的结冰冰形.      

飞机旋转表面结冰数值模拟

在Shallow-Water结冰热力学模型的基础上, 搭建了一套适用于三维旋转表面的非稳态结冰模型, 采用Jacobi迭代法和Gauss-Seidel迭代法对表面结冰进行非稳态迭代求解。用所提模型对简化后的旋转桨叶模型进行计算, 并与FENSAP软件结果进行对比, 验证了所提模型的准确性, 分析了转速、水滴直径和液态水含量等因素对旋转表面结冰冰形和水膜流动的影响。结果表明:随着转速的增加, 结冰范围和水膜覆盖范围偏移愈加明显;结冰范围和水膜覆盖范围随水滴直径增长逐渐增加, 水膜厚度也逐渐增大;结冰厚度和水膜厚度随液态水含量增长而相应增加, 水膜覆盖范围也明显变大。      

发动机进气道短舱前缘结冰三维模拟研究

为计算发动机进气道短舱前缘的结冰冰形,在Messinger结冰热力学模型的基础上发展了一套纯三维的表面溢流水流动结冰模型。考虑空气摩擦力为表面溢流水流动的主要驱动力,用空气对进气道表面的剪切力来确定溢流水的流动方向及流量分配。为求解溢流水结冰模型,发展了一套重复查找表面控制体状态的方法,能够快速完成整个三维表面的计算。用该方法对某三维发动机进气道进行计算,得到了三维结冰冰形,并将计算结果与FENSAP-ICE计算结果进行了对比,结果显示两者的冰形轮廓基本一致,仅在冰角处存在差异,表明本文三维发动机结冰计算模型与计算方法是有效的,其计算精度与FENSAP-ICE结果相当。     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。