摘要
数值模拟是研究飞机结冰的重要手段,为获取水滴撞击特性,通常采用欧拉法求解水滴运动方程。液态水含量在空间分布上存在间断,会导致水滴欧拉方程的求解极易发散;若采用一阶格式,虽然稳定性得到提高,但计算精度难以保证满足工程设计需求;有研究者通过添加扩散项以提升收敛性,但会人为增加数值耗散。针对这一问题,本文采用一阶迎风格式对水滴欧拉方程的对流通量进行离散,并通过线性重构与限制器函数保证计算精度达到二阶,建立一种基于迎风离散格式的水滴欧拉方程二阶求解方法。通过对比NACA0012、NACA23012和NLF0414翼型水滴撞击特性计算结果与FENSAP-ICE软件的结果,验证了该方法的有效性与计算精度。
Numerical simulation is an important means to study aircraft icing.In order to obtain water droplet impingement characteristics,Euler method is usually used to solve the control equation of droplet movement.Due to the discontinuity in the spatial distribution of liquid water content,the solution of the Euler equation of water droplet is easily divergent.If the first-order format is adopted,although the stability is improved,the calculation accuracy is difficult to meet the engineering design requirements.Some researchers increase convergence by adding diffusion terms,but artificially increase numerical dissipation.To solve this problem,the convective flux of the water droplet Euler equation is discretized by a first-order upwind scheme,and the calculation accuracy is guaranteed to reach the second-order by linear reconstruction and limiter functions.The effectiveness and accuracy of the method are verified by comparing the calculation results of NACA0012 NACA23012 and NLF0414 airfoils droplet impingement characteristics with FENSAP-ICE software.
作者
宁义君
顾洪宇
朱东宇
乔龙
李艳亮
Ning Yijun;Gu Hongyu;Zhu Dongyu;Qiao Long;Li Yanliang(Liaoning Provincial Key Laboratory of Aircraft Ice Protection,AVIC Aerodynamics Research Institute,Shenyang 110034,China;Shenyang Key Laboratory of Computational Fluid Dynamics,AVIC Aerodynamics Research Institute,Shenyang 110034,China)
出处
《航空科学技术》
2024年第1期105-114,共10页
Aeronautical Science & Technology
基金
中国-意大利绿色航空技术“一带一路”联合实验室建设与飞机低能耗防冰技术联合研究(2022YFE0203700)。
关键词
水滴撞击特性
水滴运动轨迹
欧拉法
迎风格式
结冰软件
water droplet impingement characteristics
water droplet trajectory
Euler method
upwind scheme
icing software