Simulations study of fluid-structure interaction of a wind turbine

This report presents the study of fluid-structure interaction (FSI) of a 5kW horizontal axis wind turbines (HAWT). The motivation of the current investigation is to address the FSI problem of HAWT in the whole flow field analysis with hybrid turbulence model method which is used to find out the flow and structure characteristics of wind turbine in operation under different tip speed ratio (TSR). The target of current work is to obtain detailed information of flow field and structure response simultaneously using the finite element analysis (FEA) method as well as providing new reliable and accurate approach which could be quickly and economically employed in wind turbine FSI problem analysis. In this report, a detailed and comprehensive literature review related to FSI problem of wind turbine has been presented and includes aerodynamic, wind turbine theory, computational fluid dynamic, fluid structure interaction and aeroelastic stability. A 5kW wind turbine is designed according to blade element momentum (BEM) method and a 3D wind turbine modeling is generated in SolidWorks. The 3D wind turbine modeling is imported into FEA software by coupling with the structure domain and the fluid domain to achieve FSI analysis and different TSR are applied to the computational fluid dynamics (CFD) analysis to gain the wind turbine operation information under different scenarios. The simulation results are investigated to reveal the influence of the FSI phenomenon on wind turbine performance. Analysis of wind turbine aeroelastic stability is conducted to prevent blade from experiencing flutter problem. A multiple objective optimization method for wind turbine design including aerodynamic and structure has been developed to improve wind turbine performance. The flow field analysis results showed that the strength of vortex at the tip and root of the blade increases as wind speed increases when wind turbine is under operating conditions. With wind speed of 10m⁄s, most of the flow over the lower surface remains attached and the wind turbine blade was running under design condition. As the wind speed increases, flow separation emerges. When the wind speed increasing to 16m⁄s, the wind turbine blade works under deep stall condition. The wind turbine power coefficient drops continuously from 0.46 to 0.22, as wind speed increased from 10m⁄s to 20m⁄s. The comparison of wind turbine power and torque coefficient under different CFD configuration shows the advantage of whole flow field analysis with hybrid turbulence model. The structural analysis results demonstrated that most of the pressure force acts at a region around the 1⁄4 chord from the leading edge of the blade under operating condition and the displacement of blade tip in flap direction is bigger than those in the other directions. Modal analysis results deduced that flap is the main type of vibration in first and second modes. The difference of natural frequency between single blade and blade rotor becomes larger as vibration order increases. The dynamic performances of wind turbine blade and rotor are totally different and show the whole flow field method is better than the one third domain method. The wind turbine aeroelastic analysis results illustrate that the bending frequency and torsion frequency are around 23.94Hz and 38.46Hz respectively under TSR=4 and both move towards 30Hz under TSR=1.7. When TSR is smaller than 1.7 the wind turbine will experience flutter due to the overlapping of bending and torsion frequencies. The multiple objectives optimization method takes the aerodynamic performance and blade loading into account to gain better design parameter. The comparison between initial and optimal design shows that optimized method could provide higher electrical energy around 1073kW•h per annum under rated wind speed. The optimized wind turbine critical flutter wind speed will increase from initial design speed of 35m/s to 38m/s under optimization design.Doctor of Philosoph