Numerical simulation of two dimensional ice accretion based on lubrication theory
-
摘要: 将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统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
-
[1] Messinger B L. Equilibrium temperature of an unheated icing surface as a function of air speed[J].Journal of the Aeronautical Sciences,1953,20(1):29-42 [2] Ruff G A, Berkowitz B M.User's manual for the NASA lewis ice accretion prediction code (LEWICE)[R].NASA CR 1990-185129,1990 [3] Myers T G, Hammond D W.A mathematical model for atmospheric ice accretion and water flow on a cold surface[J].International Journal of Heat and Mass Transfer,2004,47(25): 54835500 [4] Myers T G, Charpin J P F,Chapman S J.The flow and solidification of a thin fluid film on an arbitrary three-dimensional surface[J].Physics of Fluids,2002,14(8):2788-2803 [5] Gent R W, Dart N P,Cansdale J T.Aircraft icing[J].Philos Trans R Soc London,Ser A,2000,358:2873-2911 [6] Verdin P, Charpin J P F,Thompson C P.Multistep results in ICECREMO2[J].Journal of Aircraft,2009,46(5):1607-1613 [7] Verdin P, Charpin J P F.Multi-stepping ice prediction on cylinders and other relevant geometries[J].Journal of Aircraft,2013,50(3):871-878 [8] Moriarty J A, Schwartz L W,Tuck E O. Unsteady spreading of thin liquid films with small surface tension[J].Physics of Fluids A,1991,3(5):733-742 [9] Brakel T W, Charpin J P F.One-dimensional ice growth due to incoming supercooled droplets impacting on a thin conducting substrate[J].International Journal of Heat and Mass Transfer,2007,50(9/10):1694-1705 [10] Wright W B, Gent R W,Guffond D.DRA/NASA/ONERA collaboration on icing research part IIprediction of airfoil ice accretion[R].NASA CR 1997-202349,1997 [11] Spalart P R. Trends in turbulence treatments[J].AIAA Paper 2000-2306,2000 [12] Shin J, Berkowitz B M.Prediction of ice shapes and their effect on airfoil performance[R].AIAA-91-0264,1991 [13] Beaugendre H, Morency F,Habashi W G.Development of a second generation in-flight icing simulation code[J].Journal of Fluids Engineering,2006,128(2):378-387 [14] 张强,曹义华, 李栋.采用欧拉两相流法对翼型表面霜冰的数值模拟[J].北京航空航天大学学报,2009,35(3): 351355 Zhang Qiang,Cao Yihua,Li Dong.Numerical method to simulate rime ice accretions on an airfoil[J].Journal of Beijing University of Aeronautics and Astronautics,2009,35(3):351-355(in Chinese) [15] Wright W B. Validation results for LEWICE2.0[R].NASA/CR 1999-208690,1999 [16] Fu P, Farzaneh M,Bouchard G.Two-dimensional modelling of the ice accretion process on transmission line wires and conductors[J].Cold Regions Science and Technology,2006,46(2):132-146 [17] 庄礼贤, 尹协远,马晖扬.流体力学[M].合肥:中国科学技术大学出版社,1991:401 Zhuang Lixian,Yin Xieyuan,Ma Huiyang.Fluid mechanics[M].Hefei:Unversity of Science and Technology of China Press,1991:401(in Chinese)
点击查看大图
计量
- 文章访问数: 765
- HTML全文浏览量: 48
- PDF下载量: 636
- 被引次数: 0