Most of the well formulated ice prediction models are capable of predicting rime ice accurately due to its smooth profile. However, glaze ice comprises of horn-ice and feathers which makes it more difficult to predict. Therefore, the aim of the present study is to develop a numerical methodology to predict glaze ice formation on NREL S809 wind turbine blade airfoil. Research shows that ANSYS FENSAP-1CE is the most suitable commercially available software for ice prediction. Hence, the goal of this research is to explore different solver setup conditions on DROP3D and ICE3D modules within FENSAP-ICE and investigate how each of the variables affect the final ice shape. Results from this investigation would enhance the understanding of the solver and help propose the best practice set-up which would be used in the subsequent research on parametric analysis involving glaze ice prediction. Flow conditions from existing wind turbine blade experiment was used for simulation and validation purposes. Over 60 different set-up conditions were tested and validated for this study but only four best cases are presented and examined in this paper based on the ice location, shape and size. Results show that Case 37 with; SLD by body forces splashing and break-up model, Gravity vector of 9.81m/s2 in the Y-axis applied to both DROP3D and ICE3D module, beading enabled in ICE3D module, using six multi-shots proved to be the most accurate set-up for glaze ice.
展开▼