航空学报 > 2009, Vol. 30 Issue (12): 2269-2274

流体力学、飞行力学与发动机

基于非结构动态网格的翼型积冰过程数值模拟

王江峰, 夏健, 田书玲   

  1. 南京航空航天大学 航空宇航学院
  • 收稿日期:2008-10-07 修回日期:2008-12-03 出版日期:2009-12-25 发布日期:2009-12-25
  • 通讯作者: 王江峰

Numerical Simulation of Ice Accretion on Airfoil Based on Unstructured Dynamic Grids

Wang Jiangfeng, Xia Jian, Tian Shuling   

  1. College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics
  • Received:2008-10-07 Revised:2008-12-03 Online:2009-12-25 Published:2009-12-25
  • Contact: Wang Jiangfeng

摘要: 讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格库塔方法。采用四阶龙格库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。

关键词: 积冰, 非结构网格, 迎风格式, 动态网格, 数值模拟

Abstract: A numerical method for simulating the ice accretion process on airfoil is presented. The inviscid outer flow around the airfoil is solved using a finite volume upwind splitting VanLeer scheme on unstructured grids with distributed parallel computation, and the time integration is a 5stage RungeKutta method. Then it is coupled with the icing model, with a fourthorder RungeKutta scheme for the description of the trajectories of supercold droplets. A simple method is introduced to calculate the icing quantity on the body surface. And the dynamic mesh technique is implemented to adjust the grids for the iced airfoil. This method can describe the icing process almost without any artificial control parameters. After the validation and verification of the numerical method by the ice accretion problem of the NACA0012 airfoil, it is used to simulate the airfoil ice shape of an airplane under icing weather conditions. The numerical results are obtained and compared well with the limited references.

Key words: ice accretion, unstructured mesh, upwind scheme, dynamic grid, numerical simulation

中图分类号: