迎风面三维积冰过程中水膜流动的计算方法

曹广州; 吉洪湖; 斯仁

doi:10.13224/j.cnki.jasp.2015.03.019

迎风面三维积冰过程中水膜流动的计算方法

doi: 10.13224/j.cnki.jasp.2015.03.019

1. 南京航空航天大学无人机研究院, 南京 210016;
2. 南京航空航天大学能源与动力学院, 南京 210016

基金项目:

中央高校基本科研业务费专项资金(NS2015096)

详细信息

作者简介:
曹广州(1981-),男,江苏赣榆人,助理研究员,博士,主要从事航空发动机应用和传热传质研究.

中图分类号: V231.1
计量
- 文章访问数: 1473
- HTML浏览量: 86
- PDF量: 799
- 被引次数: 0
出版历程
- 收稿日期: 2014-07-01
- 刊出日期: 2015-03-28

Computational methodology of water film flow in three-dimensional ice accretion on upwind surface

1. Research Institute of Unmanned Aircraft, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China;
2. College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

摘要

摘要: 对影响迎风面三维积冰过程的冰层表面水膜流动进行分析,建立描述其流动的数学模型,引入在流动前锋下游预设薄水膜的方法处理其在表面干区的流动推进过程,提出超过临界厚度的水膜被吹除的方法处理角状冰下游局部水膜被吹离表面的现象.通过对翼型-平板结构中翼型上的冰形对比和水膜流动规律分析,验证了该水膜流动数学模型的合理性和计算方法的可行性.研究表明:与Messinger模型相比,积冰模型中引入水膜流动后,可以较好地模拟三维明冰积冰,同时明冰极限的模拟结果也更接近试验值;冰层表面的水膜厚度约为10^-5m,流动速度约为10^-2m/s.
- 三维积冰 /
- 水膜流动 /
- 数学模型 /
- 流动前锋 /
- 水膜吹离
Abstract: The thin water film flow on the ice layer affecting the three-dimensional (3-D) ice accretion on the upwind surface was analyzed. A mathematic model for simulating the flow was developed. A very thin precursor water film beyond the contact line was introduced to deal with the film flow on the dry surface. To treat with the phenomenon that some water film was blew off the surface behind the ice horn, a computational methodology was introduced that the top part of the water film which exceeding the critical height was blew away. The ice shape contrast and water film flow analysis on the airfoil-plane indicates the rationality of the mathematic model and the feasibility of the computational methodology. The research indicates: icing model involving the water film flow can simulate 3-D glaze ice well and provide the ice layer boundary closer to the experimental result than Messinger model. The water film height on the ice layer is about 10^-5m and its velocity is about 10^-2m/s.
- three-dimensional ice accretion /
- water film flow /
- mathematic model /
- flow contact line /
- film blow off

HTML

参考文献(22)

[1]	裘燮刚,韩凤华.飞机防冰系统[M].北京:航空专业教材编审组,1996.
[2]	Bourgault Y,Habashi W G,Beaugendre H.Development of a shallow water icing model in FENSAP-ICE[R].AIAA 99-0246,1999.
[3]	Myers T G,Charpin J P,Thompson C P.Slowly accreting ice due to supercooled water impacting on a cold surface[J].Physics of Fluids,2002,14(1):240-256.
[4]	Messinger B L.Equilibrium temperature of an unheated icing surface as a function of airspeed[J].Journal of Aeronautical Science,1953,20(1):29-42.
[5]	Wright W B,Gent R W,Guffond D.DRA/NASA/ONERA collaboration on icing research:Part Ⅱ prediction of airfoil ice accretion[R].NASA-CR-202349,1997.
[6]	Hedde T,Guffond D.Development of a three-dimensional icing code,comparison with experimental shapes[R].AIAA 92-0041,1992.
[7]	Myers T G,Charpin J P,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.
[8]	陈维建.飞机机翼结冰的数值模拟研究[D].南京:南京航空航天大学,2007. CHEN Weijian.The investigation on the numerical simulation for airfoil ice accretion[D].Nanjiing:Nanjing University of Aeronautics and Astronautics,2007.(in Chinese)
[9]	曹广州.迎风表面三维积冰的数学模型与计算方法研究[D].南京:南京航空航天大学,2011. CAO Guangzhou.An icing model and the corresponding calculational method for simulating three dimensional ice accretion on the upwind surface[D].Nanjiing:Nanjing University of Aeronautics and Astronautics,2011.(in Chinese)
[10]	Wang G,Rothmayer A P.Thin water films driven by air shear stress through roughness[J].Computer and Fluid,2009,38(2):235-246.
[11]	Rothmayer A P,Hu H.Linearized solutions of three-dimensional condensed layer films[R].San Diego:5th AIAA Atmospheric and Space Environments Conference,2013.
[12]	Zhang K,Johnson B,Rothmayer A P,et al.An experimental investigation on wind-driven rivulet/film flow over a NACA0012 airfoil by using digital image projection technique[R].AIAA-2014-0741,2014.
[13]	Silva G A,Silvares O M,Zerbini E J.Water film breakdown and rivulets formation effects on thermal anti-ice operation simulation[R].AIAA-2006-3758,2006.
[14]	Silva G A,Silvares O M,Zerbini E J.Differential boundary-layer analysis and runback water flow model applied to flow around airfoils with thermal anti-ice[R].AIAA-2009-3967,2009.
[15]	于意奇.大尺度平板水膜流动行为的数值模拟和试验研究[D].上海:上海交通大学,2012. YU Yiqi.Experimental and numerical study on film flow behavior on large scale flat plate[D].Shanghai:Shanghai Jiao Tong University,2012.(in Chinese)
[16]	周孝清,陈沛霖.竖直壁面上的水膜流动[J].中南工学院学报,1996,13(1):1-7. ZHOU Xiaoqing,CHEN Peilin.Characteristics of thin water film flow along a vertical plate[J].Journal of Central-South Institute of Technology,1996,13(1):1-7.(in Chinese)
[17]	张群,童一俊,任海刚,等.竖壁降膜流动特性与稳定性的数值研究[J].西北工业大学学报,2013,31(5):733-736. ZHANG Qun,TONG Yijun,REN Haigang,et al.Numerical investigating flow characteristics and stability of film falling down a vertical plate[J].Journal of Northwestern Polytechnical University,2013,31(5):733-736.(in Chinese)
[18]	Hansman R H,Turnock S T.Investigation of surface water behavior during glaze ice accretion[J].Jounal of Aircraft,1989,26(2):140-147.
[19]	曹广州,吉洪湖,胡娅萍,等.模拟飞机迎风面三维积冰的数学模型研究[J].航空动力学报,2011,26(9):1953-1963. CAO Guangzhou,JI Honghu,HU Yapping,et al.An icing model for simulating three dimensional ice accretion on the upwind surface of a plane[J].Journal of Aerospace Power,2011,26(9):1953-1963.(in Chinese)
[20]	Moriarty J A,Schwartz L W.Unsteady spreading of thin liquid films with small surface tension[J].Physics of Fluids A,1991,3(5):733-742.
[21]	Ruff G A,Berkowitz B M.User's manual for the NASA Lewis ice accretion prediction code (LEWICE)[R].NASA CR-185129,1990.
[22]	Shin J,Bond T H.Results of an icing test on a NACA0012 airfoil in the NASA Lewis icing research tunnel[R].AIAA 92-0647,1992.