Influence of liquid water content onwind turbine blade icing by numericalsimulation
LI Yan1*, SUN Ce1, JIANG Yu1, YI Xian2, GUO Wenfeng1, WANG Shaolong3, FENG Fang1
1. Heilongjiang Provincial Key Laboratory of Technology and Equipment for Utilization of Agricultural Renewable Resources in Cold Region, Northeast Agricultural University, Harbin, Heilongjiang 150030, China; 2. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuang 621000, China; 3. College of Metallurgy & Energy, North China University of Science and Technology, Tangshan, Hebei 063009, China
Influence of liquid water content on wind turbine blade icing by numerical simulation
LI Yan1*, SUN Ce1, JIANG Yu1, YI Xian2, GUO Wenfeng1, WANG Shaolong3, FENG Fang1
1. Heilongjiang Provincial Key Laboratory of Technology and Equipment for Utilization of Agricultural Renewable Resources in Cold Region, Northeast Agricultural University, Harbin, Heilongjiang 150030, China; 2. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang, Sichuang 621000, China; 3. College of Metallurgy & Energy, North China University of Science and Technology, Tangshan, Hebei 063009, China
摘要 In order to research the influence of liquid water content(LWC)on blade icing of wind turbine, a numerical simulation method for blade icing was established. The numerical simulation was based on low speed viscous N-S equation. The trajectory equation of water droplets was established by Lagrangian method. The mass and energy conservation equations of the water droplets impacting on the surface of the blade were solved based on control body theory. Three sections along blade span wise of a 1.5 MW wind turbine were decided to simulate icing. Five kinds of LWC were selected for simulation including 0.2,0.4,0.6,0.8 and 1.0 g/m3 under two ambient temperatures of -10 ℃ and -20 ℃. The medium volume droplet diameter(MVD)was 30 SymbolmA@m. The simulations included icing shape on blade surface, dimensionless icing area and dimensionless maximum stagnation thickness. Furthermore, the flow fields around both the iced blade airfoil and the original one were simulated and analyzed. Accor-ding to the results, the typical icing characteristics of icing shape, icing area and thickness were greatly affected by the difference of LWCs. This study can provide theoretical reference for the research on anti-icing and de-icing of wind turbine blade.
Abstract：In order to research the influence of liquid water content(LWC)on blade icing of wind turbine, a numerical simulation method for blade icing was established. The numerical simulation was based on low speed viscous N-S equation. The trajectory equation of water droplets was established by Lagrangian method. The mass and energy conservation equations of the water droplets impacting on the surface of the blade were solved based on control body theory. Three sections along blade span wise of a 1.5 MW wind turbine were decided to simulate icing. Five kinds of LWC were selected for simulation including 0.2,0.4,0.6,0.8 and 1.0 g/m3 under two ambient temperatures of -10 ℃ and -20 ℃. The medium volume droplet diameter(MVD)was 30 SymbolmA@m. The simulations included icing shape on blade surface, dimensionless icing area and dimensionless maximum stagnation thickness. Furthermore, the flow fields around both the iced blade airfoil and the original one were simulated and analyzed. Accor-ding to the results, the typical icing characteristics of icing shape, icing area and thickness were greatly affected by the difference of LWCs. This study can provide theoretical reference for the research on anti-icing and de-icing of wind turbine blade.
LI Yan*, SUN Ce, JIANG Yu, YI Xian, GUO Wenfeng, WANG Shaolong, FENG Fang. Influence of liquid water content onwind turbine blade icing by numericalsimulation[J]. 排灌机械工程学报, 2019, 37(6): 513-520.
LI Yan*, SUN Ce, JIANG Yu, YI Xian, GUO Wenfeng, WANG Shaolong, FENG Fang. Influence of liquid water content on wind turbine blade icing by numerical simulation. Journal of Drainage and Irrigation Machinery Engin, 2019, 37(6): 513-520.
[1]ZHAO B, MA H, ZHAO Y, et al. Analysis on application of vertical axis wind turbine in Qinghai-Tibet Plateau[J]. Journal of drainage and irrigation machinery engineering, 2018, 36(3): 250-256.(in Chinese)[2]ZHAO Z, QIAN S, ZHENG Y, et al. Enhancement approaches of aerodynamic performance of lift-type vertical axis wind turbine considering small angle of attack[J]. Journal of drainage and irrigation machinery enginee-ring, 2018, 36(2): 146-153.(in Chinese)[3]YAN L, SHOUYANG Z, KOTARO T, et al. Starting performance effect of a truncated-cone-shaped wind gathering device on small-scale straight-bladed vertical axis wind turbine[J]. Energy conversion and management, 2018, 167: 70-80.[4]YAN L, CE S, YU J, et al. Scaling method of the rotating blade of a wind turbine for a rime ice wind tunnel test[J]. Energies, 2019, 12: 626-627.[5]YAN L, KOTARO T, FANG F, et al. A wind tunnel experimental study of icing on wind turbine blade airfoil[J]. Energy conversion and management, 2014, 85: 591-595.[6]YAN L, SHAOLONG W, CE S, et al. Icing distribution of rotating blade of horizontal axis wind turbine based on Quasi-3D numerical simulation[J]. Thermal science, 2018, 22: 1191-1201.[7]RUFF G A. Quantitative comparison of ice accretion shapes on airfoils [J]. Journal of aircraft, 2002, 39(3): 418-426.[8]HOMOLA M C. Effect of atmospheric temperature and droplet size variation on ice accretion of wind turbine blades [J]. Journal of wind engineering and industrial aerodynamics, 2010, 98: 724-729.[9]YI X, WANG K C, MA H L, et al. Computation of icing and its effect of horizontal axis wind turbine [J]. Acta energiae solaris sinica, 2014, 35(6): 1052-1058.(in Chinese)[10]LI Y, WANG S D, LIU Q D, et al, Characteristics of ice accretions on blade of the straight-bladed vertical axis wind turbine rotating at low tip speed ratio [J]. Cold regions science and technology, 2018, 145: 1-13.[11]LI Y, SUN C, JIANG Y, et al. Temperature effect on icing distribution near blade tip of large-scale HAWT by numerical simulation[J]. Advances in mechanical engineering, 2018, 10(11), 1-13.
