垂直轴风轮是风能转换的主要设备之一，CFD方法是对其复杂绕流场进行数值模拟的有效手段。实度是影响垂直轴风轮性能的主要参数之一，由叶片数、翼型弦长及转子半径三个参数共同决定。为研究不同参数组合工况下垂直轴风轮的气动力学行为，采用SST-DDES方法对垂直轴风轮二维模型进行了数值模拟，并与已有研究结果对比分析，评估了SST-DDES预测垂直轴风轮气动性能的准确性。其次，对模拟结果有关键影响的计算域、网格精度以及时间步长进行了无关性研究，减少计算资源消耗并增进数值模拟的可靠性。最后，对多组不同实度的工况进行数值模拟，得到每种工况下功率系数与尖速比的关系曲线。结果表明，风轮实度在0.2 ~ 0.5之间时其功率系数峰值最大；两叶片风轮比多叶片风轮具有更高的功率系数峰值和更宽的有效尖速比范围；增大风轮的半径可显著提高风轮的功率峰值系数；风轮实度及半径一定时，不同叶片数和弦长组合具有相似的功率系数曲线
Vertical axis wind turbine (VAWT) is one of the main equipment for wind energy conversion, and computational fluid dynamics (CFD) is the most effective method to study its aerodynamic behavior. Rotor solidity, which is defined by three physical parameters, blade number, chord length and rotor radius, has been identified as a key parameter on the performance of a VAWT. However, there is a lack of systematic studies on these three parameters in the literature. To examine the effect of the solidity, two-dimensional numerical simulations were performed based on the shear stress transport-delayed detached eddy simulation (SST-DDES) method. The SST-DDES model was validated against published experiment data by comparing the power coefficient (CP) and drag coefficient (CD). Besides, domain size, mesh density and time step were analyzed and optimized in order to reduce the computational time cost and to improve the reliability of the simulations. Cases of different solidities were computed and curves of CP were obtained. There exists a maximum CP when the solidity value varies between 0.2 and 0.5. VAWTs with two blades achieve a higher peak of CP and a wider range of effective tip speed ratio. Increasing the radius of the rotor can significantly increase the peak CP of VAWTs. By keeping the solidity and rotor radius constant, VAWTs with different blade numbers have similar CP curves.
[1] PARASCHIVOIU I. Wind turbine design: With emphasis on Darrieus concept[M]. Montreal, Canada: Polytechnic International Press, 2002.
[2] FERREIRA C S, VAN KUIK G, VAN BUSSEL G, et al. Visualization by PIV of dynamic stall on a vertical axis wind turbine[J]. Experiments in fluids, 2009, 46(1): 97-108. DOI: 10.1007/s00348-008-0543-z.
[3] LIANG X T, FU S, QU B X, et al. A computational study of the effects of the radius ratio and attachment angle on the performance of a Darrieus-Savonius combined wind turbine[J]. Renewable energy, 2017, 113: 329-334. DOI: 10.1016/j.renene.2017.04.071.
[4] HOWELL R, QIN N, EDWARDS J, et al. Wind tunnel and numerical study of a small vertical axis wind turbine[J]. Renewable energy, 2010, 35(2): 412-422. DOI: 10.1016/j.renene.2009.07.025.
[5] LI Q A, MAEDA T, KAMADA Y, et al. Effect of number of blades on aerodynamic forces on a straight-bladed Vertical Axis Wind Turbine[J]. Energy, 2015, 90: 784-795. DOI: 10.1016/j.energy.2015.07.115.
[6] LI S M, LI Y. Numerical study on the performance effect of solidity on the straight-bladed vertical axis wind turbine[C]//Proceedings of 2010 Asia-Pacific Power and Energy Engineering Conference. Chengdu, China: IEEE, 2010: 1-4. DOI: 10.1109/APPEEC.2010.5449269.
[7] STRICKLAND J H. Darrieus turbine: A performance prediction model using multiple streamtubes[R]. Albuquerque, United States: Sandia Laboratory, 1975.
[8] HOLME O. A contribution to the aerodynamic theory of the vertical-axis wind turbine[C]//International Symposium on Wind Energy Systems. Cambridge, England: British Hydromechanics Research Association, 1977.
[9] WANG S Y, INGHAM D B, MA L, et al. Numerical investigations on dynamic stall of low Reynolds number flow around oscillating airfoils[J]. Computers & fluids, 2010, 39(9): 1529-1541. DOI: 10.1016/j.compfluid.2010. 05.004.
[10] SIMÃO FERREIRA C J, BIJL H, VAN BUSSEL G, et al. Simulating dynamic stall in a 2D VAWT: Modeling strategy, verification and validation with Particle Image Velocimetry data[J]. Journal of physics: conference series, 2007, 75(1): 012023. DOI: 10.1088/1742-6596/75/1/012023.
[11] SPALART P R, DECK S, SHUR M L, et al. A new version of detached-eddy simulation, resistant to ambiguous grid densities[J]. Theoretical and computational fluid dynamics, 2006, 20(3): 181-195. DOI: 10.1007/s00162- 006-0015-0.
[12] MENTER F R, KUNTZ M, LANGTRY R. Ten years of industrial experience with the SST turbulence model[M]//HARJALI? K, NAGANO Y, TUMMERS M. Turbulence, Heat and Mass Transfer 4. Antalya, Turkey: Begell House, Inc., 2003.
[13] STRELETS M. Detached eddy simulation of massively separated flows[C]//39th Aerospace Sciences Meeting and Exhibit, Aerospace Sciences Meetings. Reno, NV, USA: AIAA, 2001: 1-18. DOI: 10.2514/6.2001-879.
[14] MCLAREN K W. A Numerical and experimental study of unsteady loading of high solidity vertical axis wind turbines[D]. Hamilton: McMaster University, 2011.
[15] AMET E, MAÎTRE T, PELLONE C, et al. 2D Numerical simulations of blade-vortex interaction in a Darrieus turbine[J]. Journal of fluids engineering, 2009, 131(11): 111103. DOI: 10.1115/1.4000258.
[16] LANEVILLE A, VITTECOQ P. Dynamic stall: The case of the vertical axis wind turbine[J]. Journal of solar energy engineering, 1986, 108(2): 140-145. DOI: 10.1115/ 1.3268081.
