飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

三维翼面结冰过程及其影响数值研究

为得到三维机翼及翼身尾构型在结冰条件大气中的结冰过程和影响,采用结构化网格对M6机翼和翼身尾构结冰过程进行数值模拟,对飞机结冰过程的数值模拟采用欧拉两相流法求解水滴撞击特性,基于Messinger热力学模型模拟冰形增长过程,对其水滴撞击特性、结冰形状及结冰后对气动特性影响进行研究.分析结果可得:机翼展向,随着弦长变小,结冰量逐渐增多,尤其对大展弦比机翼更为明显.同时机身前端迎风面也是比较容易积冰的区域.     

大型商用运输机机翼增升构型水滴撞击特性计算

基于欧拉-拉格朗日数值模拟框架,发展了一套适合求解三维多体模型水滴撞击特性的计算方法。采用粒子统计与面积比率计算相结合的方式,实现了水滴收集率的高效精确预测。在常规粒径条件下,该方法计算结果与国外标模试验数据吻合良好。应用上述算法,开展了某大型商用运输机增升构型的水滴撞击特性数值模拟研究,计算结果显示:水滴在缝翼背面缝道处出现聚集和轨迹交汇现象,聚集的水滴撞击在缝翼背面后缘,致使当地的水滴收集率出现了严重的激增,这一现象对飞机防冰系统设计提出了新的要求;另外,三维效应对增升构型的撞击特性影响显著,其中主翼下翼面水滴收集率沿展向分布存在较大差异。     

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

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

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

CJ818客机机翼防冰范围的确定及防冰热载荷计算

通过数值计算的方法,对CJ818客机机翼周围的流场和水滴运动轨迹进行了计算;介绍了水滴撞击特性的计算方法,编写程序计算了机翼的水滴撞击特性,并据此确定出机翼的防冰范围;分析了防冰表面热流的组成,并给出了每种热流的计算方法,在此基础上编写程序计算得到了CJ818客机机翼的防冰热载荷。     

三维笛卡尔网格的结冰数值模拟

<正>飞机结冰现象会对飞机安全飞行带来严重威胁,因此对飞机结冰现象研究的开展具有必要性。目前研究结冰的方法主要有两种:试验法和数值模拟法。数值模拟法,则是利用计算机建立机翼几何模型、划分计算网格,计算水滴撞击特性,进行机翼表面的传热、传质计算得到冰形。数值模拟不需要在现实中去模拟真实的飞机飞行情况,仅仅需要在数学上构建相关模型,这在成本上相对低廉,同时也无需承担试验风险,目前在结冰研究中成为了不可或缺的重要方法。国外三维结冰数值模拟研究开     

风洞侧壁干扰对机翼表面霜冰的影响

采用数值模拟的方法研究了不同阻塞度下侧壁干扰对NACA0012翼型表面霜冰生长的影响。通过对流场和水滴撞击特性的研究发现:侧壁干扰会压缩来流空气,致使更多水滴跟随气流向机翼前缘汇集,增大了翼面水滴收集系数,且翼面水滴总收集量与阻塞度成正比、与液滴雷诺数成反比;同时会改变水滴的撞击方向,使水滴撞击极限和最大水收集系数的位置发生偏移。采用多步法对结冰特性的研究发现:在基准计算工况下,当阻塞度为20%时,侧壁干扰效应带来结冰量的相对附加增量约为16%;在同一阻塞度下,各时间步长内的结冰相对附加增量基本相等。提出了一种表征侧壁干扰对水滴收集影响强弱的因子,该影响因子与水滴收集相对附加增量具有较高的线性关系。      

飞机热气防冰系统与冰脊预测的数值模拟

基于流固耦合传热的思想建立了一套飞机热气防冰系统的的数值模拟方法,并将其与积冰热力学模型结合起来,实现了热气防冰系统开启时的机翼积冰预测.采用格心格式有限体积法求解N-S方程获得防冰腔与外流场;通过欧拉法在外流场的基础上获得过冷水滴撞击特性;求解三维热传导偏微分方程获得蒙皮的传热特性;采用交接面插值的方法实现防冰腔到外流场的热量传递;建立了考虑三维溢流效应的积冰热力学模型并在此基础上开展了机翼冰脊的数值预测.数值模拟结果表明:热气防冰系统开启时加热机翼表面温度最高可达308K,加热区后的上下机翼表面均有冰脊形成,通过对结果的分析表明该方法是合理可行的.      

基于CFD的机翼结冰过程分析

采用一种适用于分析机翼结冰过程的数值模拟方法,该方法在飞机结冰气象条件下采用计算流体力学（Computational Fluid Dynamics,CFD）对机翼流场进行计算,获得流场中二维机翼翼型的水滴撞击特性,并遵循质量和能量守恒方程,计算出机翼表面的结冰量,预测出机翼的积冰形状。所得冰形与文献中的试验数据吻合良好,表明数值模拟方法有效。另外还分析了结冰对机翼升力、阻力和压力系数的影响,其结果可以指导飞机机翼防／除冰系统的设计研究     

翼型表面明冰的数值模拟

提出一种翼型表面明冰的数值模拟方法.基于欧拉两相流理论,建立空气、水滴两相流控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击;采用考虑了粗糙度影响的边界层积分法求解结冰表面的对流换热特性,基于经典Messinger模型提出明冰数值模拟的热力学模型;利用冰层时间推进法模拟积冰过程,采用积冰法向生长假设生成积冰外形.对NACA 0012翼型在不同环境温度下的明冰冰形进行了预测,并与实验结果进行了对比.      

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

发动机旋转整流帽罩的水滴撞击特性分析

基于计算流体力学软件FLUENT,采用拉格朗日法,提出了一种求取旋转部件表面水滴收集系数的方法,并采用用户自定义函（UDF）予以实现。将该方法运用于发动机旋转整流帽罩的水滴撞击特性研究,结果表明：随着帽罩锥体锥角及转速的增大,帽罩表面的局部水收集系数沿其表面外形线下降较快,且撞击区域变小。     

横向气流中液体圆柱射流的破碎特性和表面波现象

采用高速摄像仪对横向气流场中的圆形液体射流的圆柱破碎过程(低Weber数)进行了实验研究.实验采用直射式喷嘴,喷孔直径为0.3mm和0.5mm,喷嘴长径比均为40.实验工质采用水.Weber数为1.7~7,液/气动量通量比为3.4~83.实验观察了表面波现象及射流破碎形成的液滴的尺寸及其速度.观测发现,表面波波长随着气流Weber数的对数的增加而线性减小.破碎后形成的液滴直径 d_p/d₀与Weber数的对数成正比.液滴的初始x方向速度与y方向速度大小与液滴直径无关.在实验范围内,液滴y方向速度均不随Weber数和液/气动量通量比的变化,约为初始射流速度的0.94倍;而液滴x方向速度约为0.085倍的气流速度.对表面波和破碎后形成的液滴尺寸及速度的研究有助于构建更精确的初始雾化模型.      

Effects of flow parameters on thermal performance of an inner-liner anti-icing system with jets impingement heat transfer

The multiple jets impingement heat transfer is widely applied in the wing anti-icing system. It is challenging to apply the similarity criterion to carry out the anti-icing experiments due to the complex flow and heat transfer behavior. In the present study, the full-scale slat model is used to carry out anti-icing experimental researches in a 2 m × 3 m icing wind tunnel of China Aerodynamics Research and Development Center. The effects of icing parameters Liquid Water Content (LWC) and Median Volume Diameter (MVD) and hot air parameters (mass flow rate and temperature) on the thermal performance of an inner-liner anti-icing system with jets impingement heat transfer are studied. The effects of the experimental parameters are analyzed in detail by combining impingement and evaporation heat transfer mechanisms. The impingement hot air mass flow rate dramatically affects the heat transfer performance of the impingement stagnation region within the range of the experimental parameters. The temperature of impingement hot air and that of wing skin are approximately linear correlated. The experimental results show the effects of LWC and MVD on water film formation and runback ice accretion. The formation of water film is analyzed by an analytical method based on the wing skin temperature difference of dry and wet air conditions.     

发动机进气帽罩防冰热载荷的数值模拟研究

对某发动机进气帽罩防冰热载荷的计算方法进行研究。利用FLUENT软件对帽罩周围的流场进行模拟,同时采用FLUENT中的离散相模型计算帽罩周围的水滴运动轨迹,得到帽罩表面的水滴撞击效率;利用数值计算结果及流动换热经验关系式分析帽罩表面的各项热流,获得该发动机帽罩的防冰热载荷。计算中还比较了不同湍流模型对防冰热载荷的影响。