北京航空航天大学学报 ›› 2014, Vol. 40 ›› Issue (10): 1442-1450.doi: 10.13700/j.bh.1001-5965.2013.0625

• 论文 • 上一篇    下一篇

基于润滑理论的二维积冰数值模拟

侯硕, 曹义华   

  1. 北京航空航天大学 航空科学与工程学院, 北京 100191
  • 收稿日期:2013-10-31 出版日期:2014-10-20 发布日期:2014-10-29
  • 作者简介:侯硕(1986-),男,山东菏泽人,博士生,houshuo@ase.buaa.edu.cn.
  • 基金资助:

    航空科学基金资助项目(2009ZA51007)

Numerical simulation of two dimensional ice accretion based on lubrication theory

Hou Shuo, Cao Yihua   

  1. School of Aeronautic Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China
  • Received:2013-10-31 Online:2014-10-20 Published:2014-10-29

摘要:

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.

关键词: 润滑理论, 薄膜流, 积冰, 翼型, 传输线, 数值模拟

Abstract:

Water film-ice layer model based on lubrication theory was extended to two-dimensional orthogonal curvilinear coordinate system on arbitrary sections. Two governing partial differential equations of water film flow and ice accretion on curved sections were obtained. The numerical methods for solving these unsteady governing equations were presented. An implicit-explicit discrete scheme was applied to obtain the algebraic form of governing equations. Examples of numerical icing calculation were verified on NACA0012 airfoil and circle cylinder section in typical aero icing and structural icing environment. The results of calculation using current model were compared with ice shapes of simulation using traditional Messinger model and results from ice wind tunnel tests. Ice shapes on airfoil obtained using this numerical method were close to icing curves calculated by Messinger model in low temperature glaze icing conditions. In relatively high temperature icing conditions, more accurate results were obtained than Messinger model compared with results of ice tunnel experiments. Ice accretion on transmission line wires which cannot be predicted properly by traditional Messinger model was also effectively predicted by current model in structural icing conditions.

Key words: lubrication theory, film flow, ice accretion, airfoil, cable, numerical simulation

中图分类号: 


版权所有 © 《北京航空航天大学学报》编辑部
通讯地址:北京市海淀区学院路37号 北京航空航天大学学报编辑部 邮编:100191 E-mail:jbuaa@buaa.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发