The parameter sensitivities affecting the flutter speed of the NREL (National Renewable Energy Laboratory) 5-MW baseline HAWT (horizontal axis wind turbine) blades are analyzed. An aeroelastic model, which compris...The parameter sensitivities affecting the flutter speed of the NREL (National Renewable Energy Laboratory) 5-MW baseline HAWT (horizontal axis wind turbine) blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases (i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area) will decrease the flutter speed. Therefore, the judicious selection of the four parameters (the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density) can increase the relative inflow speed at the blade section associated with the onset of flutter.展开更多
It is important to study the pressure distribution on the blade and in the adjacent area while searching the power augmentation theory with adding a tip vane to the wind turbine. This paper shows the CFD simulation re...It is important to study the pressure distribution on the blade and in the adjacent area while searching the power augmentation theory with adding a tip vane to the wind turbine. This paper shows the CFD simulation relationship of the pressure distribution on the rotor blade and in the adjacent area, after calculating the pressure of the different chordwise and spanwise point on the blade with the tip vane-V(8.8×8) and without the tip vane under tip speed ratio λ 4.5. Combining the isobaric section figure in certain location, it can be seen that the tip vane improve the pressure difference between pressure and suction surface. The most influenced zone is found and these can further display the power augmentation theory of the wind turbine using the tip vane. The simulation calculation was based on N-S equations. 3-D, steady, implicit solver was chosen. Turbulence model was k-ω SST. Discretization scheme is SECOND ORDER UPWIND. Pressure-velocity coupling was a typical SIMPLE scheme. In the whole grid system, two-divided grid formation was adopted, that is, inner region and outer region. Inner region including rectangular solid blade and neighboring, outer region is semi-cylinder. There were together 720,000 nodes with tetra-prism unstructured mesh.展开更多
基金Project(2015B37714)supported by the Fundamental Research Funds for the Central Universities of ChinaProject(51605005)supported by the National Natural Science Foundation of China+1 种基金Project(ZK16-03-03)supported by the Open Foundation of Jiangsu Wind Technology Center,ChinaProject([2013]56)supported by the First Group of 2011 Plan of Jiangsu Province,China
文摘The parameter sensitivities affecting the flutter speed of the NREL (National Renewable Energy Laboratory) 5-MW baseline HAWT (horizontal axis wind turbine) blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases (i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area) will decrease the flutter speed. Therefore, the judicious selection of the four parameters (the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density) can increase the relative inflow speed at the blade section associated with the onset of flutter.
基金Project 50566001 supported by NSFCProject 200308020207 supported by Inner Mongolia Autono- mous Region Natural Science Foundation of China.
文摘It is important to study the pressure distribution on the blade and in the adjacent area while searching the power augmentation theory with adding a tip vane to the wind turbine. This paper shows the CFD simulation relationship of the pressure distribution on the rotor blade and in the adjacent area, after calculating the pressure of the different chordwise and spanwise point on the blade with the tip vane-V(8.8×8) and without the tip vane under tip speed ratio λ 4.5. Combining the isobaric section figure in certain location, it can be seen that the tip vane improve the pressure difference between pressure and suction surface. The most influenced zone is found and these can further display the power augmentation theory of the wind turbine using the tip vane. The simulation calculation was based on N-S equations. 3-D, steady, implicit solver was chosen. Turbulence model was k-ω SST. Discretization scheme is SECOND ORDER UPWIND. Pressure-velocity coupling was a typical SIMPLE scheme. In the whole grid system, two-divided grid formation was adopted, that is, inner region and outer region. Inner region including rectangular solid blade and neighboring, outer region is semi-cylinder. There were together 720,000 nodes with tetra-prism unstructured mesh.