农业工程学报
農業工程學報
농업공정학보
2015年
16期
78-85
,共8页
石亚丽%左红梅%杨华%周捍珑%沈文忠
石亞麗%左紅梅%楊華%週捍瓏%瀋文忠
석아려%좌홍매%양화%주한롱%침문충
风力机%数值分析%模型%偏航角%计算流体力学%动态失速%速度分布
風力機%數值分析%模型%偏航角%計算流體力學%動態失速%速度分佈
풍력궤%수치분석%모형%편항각%계산류체역학%동태실속%속도분포
wind turbines%numerical analysis%models%the yaw angle%computational fluid dynamics%dynamic stall%the velocity distribution
偏航工况水平轴风力机存在典型的动态特性,为了提高动态载荷特性的预测精度,该文采用计算流体力学方法(computational fluid dynamics,CFD)研究了MEXICO(model experiments in controlled conditions)风轮在偏航角0、15°、30°、45°工况下的整机气动性能。数值模拟得到的叶片截面压力系数分布、载荷系数随方位角变化规律以及轴向入流时速度分布与试验测量值均吻合较好。当偏航角在30°以内时,采用CFD方法计算的轴向载荷系数的相对误差在±5%以内,切向载荷系数的相对误差在±15%以内;当偏航角达到45°时,轴向载荷系数的相对误差超过±15%,切向载荷系数的相对误差接近±30%,同时偏航运行时速度分布与试验测量相差较大。偏航运行时叶根处的翼型升阻力迟滞特性较叶尖处显著,但叶根处攻角变化范围小于叶尖处。采用动量叶素法进行风力机性能预测时必需充分考虑该特性。该研究为工程预测模型的建立和偏航工况风力机设计运行提供了参考。
偏航工況水平軸風力機存在典型的動態特性,為瞭提高動態載荷特性的預測精度,該文採用計算流體力學方法(computational fluid dynamics,CFD)研究瞭MEXICO(model experiments in controlled conditions)風輪在偏航角0、15°、30°、45°工況下的整機氣動性能。數值模擬得到的葉片截麵壓力繫數分佈、載荷繫數隨方位角變化規律以及軸嚮入流時速度分佈與試驗測量值均吻閤較好。噹偏航角在30°以內時,採用CFD方法計算的軸嚮載荷繫數的相對誤差在±5%以內,切嚮載荷繫數的相對誤差在±15%以內;噹偏航角達到45°時,軸嚮載荷繫數的相對誤差超過±15%,切嚮載荷繫數的相對誤差接近±30%,同時偏航運行時速度分佈與試驗測量相差較大。偏航運行時葉根處的翼型升阻力遲滯特性較葉尖處顯著,但葉根處攻角變化範圍小于葉尖處。採用動量葉素法進行風力機性能預測時必需充分攷慮該特性。該研究為工程預測模型的建立和偏航工況風力機設計運行提供瞭參攷。
편항공황수평축풍력궤존재전형적동태특성,위료제고동태재하특성적예측정도,해문채용계산류체역학방법(computational fluid dynamics,CFD)연구료MEXICO(model experiments in controlled conditions)풍륜재편항각0、15°、30°、45°공황하적정궤기동성능。수치모의득도적협편절면압력계수분포、재하계수수방위각변화규률이급축향입류시속도분포여시험측량치균문합교호。당편항각재30°이내시,채용CFD방법계산적축향재하계수적상대오차재±5%이내,절향재하계수적상대오차재±15%이내;당편항각체도45°시,축향재하계수적상대오차초과±15%,절향재하계수적상대오차접근±30%,동시편항운행시속도분포여시험측량상차교대。편항운행시협근처적익형승조력지체특성교협첨처현저,단협근처공각변화범위소우협첨처。채용동량협소법진행풍력궤성능예측시필수충분고필해특성。해연구위공정예측모형적건립화편항공황풍력궤설계운행제공료삼고。
A typical dynamic characteristic of horizontal axis wind turbine shows up under yaw condition. Prediction accuracy is low for momentum-blade element theory and related engineering prediction model. In order to improve the prediction accuracy of dynamic load characteristics, the whole wind turbine models, based on the experiment about MEXICO (model experiments in controlled conditions) rotor in 2006, are established by three-dimensional software called Pro/E. under different yaw conditions, i.e. yaw angle of 0, 15, 30 and 45 degree. ICEM CFD (integrated computer engineering and manufacturing code for computational fluid dynamics) is applied to grid division. The rotating domain containing rotor part is meshed into hexahedral grids, and the static domain containing part of wheel hub, tower and outflow field is meshed into tetrahedral grids. When the grid size of the first layer of blade surface is set as 5×10-6 m to ensure the first dimensionless size near the wall Y+<0.5 on the wall, the 2 numbers of grids are determined by the error of axial load on the airfoil in the 60%section of blades, which respectively are 6 572 451 and 2 961 385. The aerodynamic performance of models under rated condition is simulated by ANSYS CFX with the turbulence model of SST (shear stress transport), high resolution is chosen as advection scheme, and transient rotor stator as the domain interface method. The results are converted into data, processed and analyzed by MATLAB. Finally the following conclusions are drawn. The distributions of pressure coefficients along the airfoil chord in different blade sections calculated by CFD method are in good agreement with the experimental measurements, and the error on the suction surface of airfoil is mainly caused by stall separation occurring on the pressure surface of airfoil. With the increasing of yaw angle, the pressure coefficients of the suction side are increasing and the location of minimum pressure coefficient moves to airfoil trailing edge slightly. For the pressure side, the pressure coefficients increase at first and then decrease, and the location of maximum pressure coefficient moves to airfoil leading edge slightly. The axial load coefficients and tangential load coefficients of blades first decrease and then increase and then decrease again with the increase of the azimuthal angle. With the increase of the yaw angle, the axial and tangential load coefficients are both reduced. When the yaw angle is within 30°, the relative error of axial load coefficients is in the range of ±5% and the relative error of tangential load coefficients is in the range of ±15%. CFD method is higher than BEM (blade element momentum) method in forecasting accuracy of dynamic load calculation. Under yaw condition, the hysteresis characteristic of airfoil lift and drag in blade root is more remarkable than blade tip, while the variation range of the angle of attack in blade root is much less than that in blade tip. This characteristic must be considered when BEM method is used to predict wind turbine performance. For axial inflow condition, CFD method can well predict the average speed, but restricted by turbulence model and the wake model, CFD calculation did not show the velocity characteristics of rotating vortex shedding from wind turbine impeller under yaw condition. The study provides a data support to build up the forecast model on the engineering and provides the basis for wind turbine design under yaw condition.