积冰密度对机翼除冰过程影响的数值研究
|
摘要
针对不同类型的积冰,采用数值方法研究了不同冰密度对于电加热除冰过程中热传导及冰融化过程的影响。选取了典型结冰条件下所生成的两种霜冰及明冰作为分析对象,以标准的三维电加热除冰模型作为数值模拟研究模型,设置温度监测点,取金属蒙皮与冰的交界面作为监测面,在电加热除冰过程中记录监测点及监测面的温度变化情况。通过比较不同密度积冰融化过程中监测点与监测面的温度变化情况,分析积冰密度对冰融化过程的影响。计算结果表明,冰密度对电热除冰过程具有明显的影响,积冰密度降低会提高积冰温度升高的速度,使得积冰初始融化的时刻提前,整个融化过程缩短。
For different types of ice accretion on the airfoil surface under different icing conditions,numerical method is used to study the ice density effect on the heat conduction and the ice melting process during electrothermal de-icing process.A standard three-dimensional electrothermal de-icing model,which is covered with different ice on the airfoil metal surface,is applied in the numerical simulation.The temperature at several points in the simulated model is monitored,and the interface between the airfoil metal surface and the ice is set as the monitoring surface.The temperature variations are recorded at these monitored points and surfaces during the simulation of electrothermal de-icing process.These temperature characteristics can be used to analyze the effect of the ice density on the ice melting.Three kinds of ice with different densities are selected as the research objects.The temperature variations of rime ice during the melting process are compared with glaze ice.The results indicate that the ice density has significant effect on the electrothermal de-icing process.When the ice density decreases,the speed of the ice temperature rising increases,together with early initial time for the ice melting,and the whole melting process is shortened.
For different types of ice accretion on the airfoil surface under different icing conditions,numerical method is used to study the ice density effect on the heat conduction and the ice melting process during electrothermal de-icing process.A standard three-dimensional electrothermal de-icing model,which is covered with different ice on the airfoil metal surface,is applied in the numerical simulation.The temperature at several points in the simulated model is monitored,and the interface between the airfoil metal surface and the ice is set as the monitoring surface.The temperature variations are recorded at these monitored points and surfaces during the simulation of electrothermal de-icing process.These temperature characteristics can be used to analyze the effect of the ice density on the ice melting.Three kinds of ice with different densities are selected as the research objects.The temperature variations of rime ice during the melting process are compared with glaze ice.The results indicate that the ice density has significant effect on the electrothermal de-icing process.When the ice density decreases,the speed of the ice temperature rising increases,together with early initial time for the ice melting,and the whole melting process is shortened.
引文
[1]Fortin G,Perron J.Wind turbine icing and de-icing[C]//47th AIAA Aerospace Sciences Meeting,Orlando,2009.
[2]Varela G G A,Castellano N E,Avila E E.A possible influence of the density of ice accretions on their heat transfer coefficient[J].Quarterly Journal of the Royal Meteorological Society,2005,131:377-379.
[3]Voller V R,Prakash C.A fixed grid numerical modelling methodology for convection-diffusion mushy region phasechange problems[J].International Journal of Heat and Mass Transfer,1987,30(8):1709-1719.
[4]Voller V R,Cross M,Markatos N C.An enthalpy method for convection-diffusion phase change[J].International Journal for Numerical Methods in Engineering,1987,24(1):271-284.
[5]Scanlon T J,Stickland M T,Lacombe M.A PIV and CFDanalysis of natural convection ice melting[C]//9th International Conference on Laser Anemometry,2001.
[6]Tan F L,Hosseinizadeh S F,Khodadadi J M,et al.Experimental and computational study of constrained melting of phase change materials(PCM)inside a spherical capsule[J].International Journal of Heat and Mass Transfer,2009,52(15-16):3464-3472.
[7]Hauk T,Roisman I,Tropea C.Investigation of the melting behaviour of ice particles in an acoustic levitator[C]//11th AIAA/ASME Joint Thermophysics and Heat Transfer Conference,2014.
[8]Xiao C H,Lin G P,Gui Y W,et al.Study on effect of ice shedding on heat transfer characteristics of electrothermal aircraft deicing[J].Acta Aerodynamica Sinica,2012,30(4):551-556.(in Chinese)肖春华,林贵平,桂业伟,等.冰脱落对电热除冰传热特性的影响研究[J].空气动力学学报,2012,30(4):551-556.
[9]Wang K X,Rong J.An investigation of phase-change heat transfer on ice freezing process of a finite dimensional slab of cold energy storage[J].Journal of Engineering Thermophysics,1998,(1):74-76.(in Chinese)王凯旋,荣军.有限厚度冻冰相变过程的传热研究[J].工程热物理学报,1998,(1):74-76.
[10]Guan J C,Zhu H,Adel M S,et al.Numerical analysis of heat transfer with phase-change in cool storage ball[J].Journal of Zhejiang University(Natural Science),1999,(1):85-89.(in Chinese)管建春,朱华,阿迪尔,等.储冰球相变传热数值计算及分析[J].浙江大学学报(自然科学版),1999,(1):85-89.
[11]Zhang M,Xu B,Chai J C,et al.Numerical calculations of phase change about ice and water using moving heat source method(MHSM)[EB/OL].Beijing:Sciencepaper Online[2008-03-25].http://www.paper.edu.cn/releasepaper/content/200803-755.
[12]Yan X.The simulation and experimental study of ice water changing between solid and liquid based on Fluent[D].Qingdao University of Science and Technology,2013.(in Chinese)鄢新.基于Fluent冰水固流转化过程的仿真模拟与实验研究[D].青岛科技大学,2013.
[13]Lei G L,Dong W,Zhu J J,et al.A new melting model in electrothermal de-icing simulation[C]//ASME Turbo Expo2015:Turbine Technical Conference and Exposition,2015.
[14]Lei G L,Dong W,Zhu J J,et al.Numerical investigation of the electrothermal de-icing process of a rotor blade[C]//SAE 2015International Conference,2015.
[15]Huang X M.Study of heat transfer and flow in porous media with phase change for some applications[D].Huazhong University of Science and Technology,2004.(in Chinese)黄晓明.多孔介质相变传热与流动及其若干应用研究[D].华中科技大学,2004.
[16]Zhu Jie.Research on heat and mass transfer in the process of phase change of porous media[D].Dalian University of Technology,2006.(in Chinese)朱杰.多孔介质内的相变传热传质过程研究[D].大连理工大学,2006.
[17]Cui H T,Cao X R,Peng P Y.Numerical simulation on melting phase change heat transfer in heat storage ball filled with metallic foam[J].Journal of Functional Materials,2010,41(8):1439-1442.(in Chinese)崔海亭,曹向茹,彭培英.泡沫金属相变材料熔化传热过程的数值分析[J].功能材料,2010,41(08):1439-1442.
[18]Du Y X,Yuan Y P,Jia D Y,et al.Experimental investigation on melting characteristics of ethanolamine-water binary mixture used as PCM[J].International Communications in Heat and Mass Transfer,2007,34(9-10):1056-1063.
[19]Voller V R,Cross M.Accurate solutions of moving boundary problems using the enthalpy method[J].International Journal of Heat and Mass Transfer,1981,24(3):545-556.
[20]Voller V,Swaminathan C,Thomas B.Fixed grid techniques for phase change problems:a review[J].International Journal for Numerical Methods in Engineering,1990,30(4):875-898.
[21]Voller V R,Swaminathan C R.Real source-based method for solidfication phase change[J].Numerical Heat Transfer,Part B:Fundamentals,1991,19(2):175-189.
[22]Wright W B.Users manual for the improved NASA Lewis ice accretion code LEWICE 1.6[R].NASA CR-198355,1995.
[23]Jones K.The density of natural ice accretions[C]//The Fourth International Conference on Atmospheric Icing of Structure,1988.
[24]Rios M.Icing Simulations Using Jones'Density Formula for Accreted Ice[C]//29th Annual Aerospace Sciences Meeting,1991.
[25]Jones K F.The density of natural ice accretions related to nondimensional icing parameters[J].Quarterly Journal of the Royal Meteorological Society,1990,116(492):477-496.
[26]Levi L,Nasello O B,Prodi F.Morphology and density of ice accreted on cylindrical collectors at low values of impaction parameter.I:Fixed deposits[J].Quarterly Journal of the Royal Meteorological Society,1991,117(500):761-782.
[27]Prodi F,Levi L,Nasello O B,et al.Morphology and density of ice accreted on cylindrical collectors at low values of impaction parameter.Ii:Rotating deposits[J].Quarterly Journal of the Royal Meteorological Society,1991,117(500):783-801.
[28]Liu W,Fan A W,Huang X M.Theory and application of porous medium heat and mass transfer[M].Chinese Science Publishing&Media Ltd.,2006:53-54.(in Chinese)刘伟,范爱武,黄晓明.多孔介质传热传质理论与应用[M].科学出版社,2006:53-54.
[29]Yaslik A D,Dewitt K J,Keith Jr T G,et al.Threedimensional numerical simulation of electrothermal deicing systems[C]//29th AlAA Aerospace Sciences Meeting,1991.
[2]Varela G G A,Castellano N E,Avila E E.A possible influence of the density of ice accretions on their heat transfer coefficient[J].Quarterly Journal of the Royal Meteorological Society,2005,131:377-379.
[3]Voller V R,Prakash C.A fixed grid numerical modelling methodology for convection-diffusion mushy region phasechange problems[J].International Journal of Heat and Mass Transfer,1987,30(8):1709-1719.
[4]Voller V R,Cross M,Markatos N C.An enthalpy method for convection-diffusion phase change[J].International Journal for Numerical Methods in Engineering,1987,24(1):271-284.
[5]Scanlon T J,Stickland M T,Lacombe M.A PIV and CFDanalysis of natural convection ice melting[C]//9th International Conference on Laser Anemometry,2001.
[6]Tan F L,Hosseinizadeh S F,Khodadadi J M,et al.Experimental and computational study of constrained melting of phase change materials(PCM)inside a spherical capsule[J].International Journal of Heat and Mass Transfer,2009,52(15-16):3464-3472.
[7]Hauk T,Roisman I,Tropea C.Investigation of the melting behaviour of ice particles in an acoustic levitator[C]//11th AIAA/ASME Joint Thermophysics and Heat Transfer Conference,2014.
[8]Xiao C H,Lin G P,Gui Y W,et al.Study on effect of ice shedding on heat transfer characteristics of electrothermal aircraft deicing[J].Acta Aerodynamica Sinica,2012,30(4):551-556.(in Chinese)肖春华,林贵平,桂业伟,等.冰脱落对电热除冰传热特性的影响研究[J].空气动力学学报,2012,30(4):551-556.
[9]Wang K X,Rong J.An investigation of phase-change heat transfer on ice freezing process of a finite dimensional slab of cold energy storage[J].Journal of Engineering Thermophysics,1998,(1):74-76.(in Chinese)王凯旋,荣军.有限厚度冻冰相变过程的传热研究[J].工程热物理学报,1998,(1):74-76.
[10]Guan J C,Zhu H,Adel M S,et al.Numerical analysis of heat transfer with phase-change in cool storage ball[J].Journal of Zhejiang University(Natural Science),1999,(1):85-89.(in Chinese)管建春,朱华,阿迪尔,等.储冰球相变传热数值计算及分析[J].浙江大学学报(自然科学版),1999,(1):85-89.
[11]Zhang M,Xu B,Chai J C,et al.Numerical calculations of phase change about ice and water using moving heat source method(MHSM)[EB/OL].Beijing:Sciencepaper Online[2008-03-25].http://www.paper.edu.cn/releasepaper/content/200803-755.
[12]Yan X.The simulation and experimental study of ice water changing between solid and liquid based on Fluent[D].Qingdao University of Science and Technology,2013.(in Chinese)鄢新.基于Fluent冰水固流转化过程的仿真模拟与实验研究[D].青岛科技大学,2013.
[13]Lei G L,Dong W,Zhu J J,et al.A new melting model in electrothermal de-icing simulation[C]//ASME Turbo Expo2015:Turbine Technical Conference and Exposition,2015.
[14]Lei G L,Dong W,Zhu J J,et al.Numerical investigation of the electrothermal de-icing process of a rotor blade[C]//SAE 2015International Conference,2015.
[15]Huang X M.Study of heat transfer and flow in porous media with phase change for some applications[D].Huazhong University of Science and Technology,2004.(in Chinese)黄晓明.多孔介质相变传热与流动及其若干应用研究[D].华中科技大学,2004.
[16]Zhu Jie.Research on heat and mass transfer in the process of phase change of porous media[D].Dalian University of Technology,2006.(in Chinese)朱杰.多孔介质内的相变传热传质过程研究[D].大连理工大学,2006.
[17]Cui H T,Cao X R,Peng P Y.Numerical simulation on melting phase change heat transfer in heat storage ball filled with metallic foam[J].Journal of Functional Materials,2010,41(8):1439-1442.(in Chinese)崔海亭,曹向茹,彭培英.泡沫金属相变材料熔化传热过程的数值分析[J].功能材料,2010,41(08):1439-1442.
[18]Du Y X,Yuan Y P,Jia D Y,et al.Experimental investigation on melting characteristics of ethanolamine-water binary mixture used as PCM[J].International Communications in Heat and Mass Transfer,2007,34(9-10):1056-1063.
[19]Voller V R,Cross M.Accurate solutions of moving boundary problems using the enthalpy method[J].International Journal of Heat and Mass Transfer,1981,24(3):545-556.
[20]Voller V,Swaminathan C,Thomas B.Fixed grid techniques for phase change problems:a review[J].International Journal for Numerical Methods in Engineering,1990,30(4):875-898.
[21]Voller V R,Swaminathan C R.Real source-based method for solidfication phase change[J].Numerical Heat Transfer,Part B:Fundamentals,1991,19(2):175-189.
[22]Wright W B.Users manual for the improved NASA Lewis ice accretion code LEWICE 1.6[R].NASA CR-198355,1995.
[23]Jones K.The density of natural ice accretions[C]//The Fourth International Conference on Atmospheric Icing of Structure,1988.
[24]Rios M.Icing Simulations Using Jones'Density Formula for Accreted Ice[C]//29th Annual Aerospace Sciences Meeting,1991.
[25]Jones K F.The density of natural ice accretions related to nondimensional icing parameters[J].Quarterly Journal of the Royal Meteorological Society,1990,116(492):477-496.
[26]Levi L,Nasello O B,Prodi F.Morphology and density of ice accreted on cylindrical collectors at low values of impaction parameter.I:Fixed deposits[J].Quarterly Journal of the Royal Meteorological Society,1991,117(500):761-782.
[27]Prodi F,Levi L,Nasello O B,et al.Morphology and density of ice accreted on cylindrical collectors at low values of impaction parameter.Ii:Rotating deposits[J].Quarterly Journal of the Royal Meteorological Society,1991,117(500):783-801.
[28]Liu W,Fan A W,Huang X M.Theory and application of porous medium heat and mass transfer[M].Chinese Science Publishing&Media Ltd.,2006:53-54.(in Chinese)刘伟,范爱武,黄晓明.多孔介质传热传质理论与应用[M].科学出版社,2006:53-54.
[29]Yaslik A D,Dewitt K J,Keith Jr T G,et al.Threedimensional numerical simulation of electrothermal deicing systems[C]//29th AlAA Aerospace Sciences Meeting,1991.