新型升阻混合型垂直轴风力机气动特性研究

考虑S型与H型垂直轴风力机的特点,设计了一种新型升阻混合型垂直轴风力机,采用CFD法计算其启动与气动性能。结果表明,原始H型垂直轴风力机数值结果与试验值在各工况下吻合良好;新型升阻混合型垂直轴风力机不同方位角下的启动力矩均大于原始H型风力机,最小及最大值分别提升232%和83.3%;S型风轮输出功率随叶片重叠比增加而减小,完全重叠时输出功率基本为0;新型升阻混合型垂直轴风力机最大风能利用率为0.298,具有更复杂的流场特性。     

风力机气动特性的浸入边界法模拟1)

风力发电机的空气动力学性能是决定风力机安全与效率的最重要因素之一。但由于影响风力机气动性能的参数众多,更加高效精确地模拟风力机气动特性一直是风力机的重要发展方向。本文提出了基于浸入边界法的风力机建模,网格离散,以及数值模拟的统一性框架。利用同伦变形来生成光滑的叶片模型,并且使用仿射变换来处理叶片的渐缩与扭转问题。首先,针对二维翼型的升阻力,检验了算法的数值精度。表明此方法对于阻力的模拟具有非常严格的一阶精度,进而提出采用理查森外推法来精确高效修正升阻力模拟结果。同时,模拟研究了拱曲度以及厚度对二维翼型升阻力的影响。随后,模拟研究了单风力机(包含塔架)在不同尖速比下的功率系数,并对塔架与叶片间的相互气动作用进行了初步分析。最后,模拟研究了双风力机在风场中不同前后间隔距离下的气动干涉问题。本文主要意义在于验证建模,离散,与数值模拟的一体化框架的有效可行性,进而为后续研究(给定约束下风力机自动优化选型)提供坚实基础。     

基于涡尾迹方法的风力机气动性能分析

基于(势流)涡尾迹方法开发了水平轴风力机叶片气动性能分析程序,采用固定尾迹涡模型和自由尾迹模型分别对气动设计性能进行计算分析,得到风力机设计工况下的涡位置、诱导因子、功率系数及扭矩系数等气动性能参数,并与设计结果对比。结果表明,涡尾迹方法能够快速准确地计算风力机叶片气动性能参数,对风力机叶片气动分析,固定尾迹涡模型较自由尾迹模型计算时间短,具有较好的实用性。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε（RNG）湍流模型,数值模拟美国可再生能源实验室（NREL）5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε(RNG)湍流模型,数值模拟美国可再生能源实验室(NREL)5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

可变形飞行器新概念升阻特性分析

飞行器的变形方式可分为三大类:局部小变形、中等变形和大尺度变形.针对大尺度变形,提出了两种可变形概念:变前掠翼布局和可伸缩翼布局.通过数值方法计算分析了不同变形外形在不同马赫数下的升力、阻力和升阻特性,从中总结出适合于各个速度段的能够获得最大升阻比的最佳变形方式.结果表明:大尺度变形能够显著改变飞行器的升阻特性,使可变形飞行器能适应较广范围飞行条件的变化,因而在全飞行周期中比传统的固定外形飞行器具有更优的性能.     

某型电动飞机气动特性数值模拟与风洞实验

为了验证某型大展弦比电动飞机气动设计的合理性以及为飞行性能及品质计算提供数据,采用有限体积法离散求解三维可压雷诺平均Navier-Stokes方程,并选用Spalart-Allmaras湍流模型对该电动飞机流场进行CFD数值模拟。结果表明,该型电动飞机气动设计合理,巡航速度升阻比最高可达23,具有较高气动效率;通过CFD数值模拟得到了全机升力系数、阻力系数和升阻比。为了验证CFD计算结果,对该型电动飞机进行了缩比模型风洞实验,结果显示,CFD数值模拟法计算结果与风洞实验结果高度吻合,说明CFD计算结果准确。该方法可为大展弦比电动飞机气动设计提供指导。     

风剪切对风力机叶片气动性能及尾迹形状的影响

风力机通常运行在非定常工况中,其气动性能及尾迹会随着工况的变化而变化.风剪切是风力机长期所处的环境,它会影响到叶片气动载荷、尾迹形状、总体性能等,分析风剪切作用下的叶片气动性能对风力机的设计有重要意义.本文采用一种时间步进自由涡尾迹(free vortex wake,FVW)方法,耦合FVW方法与风剪切模型,计算不同风...     

3MW叶片翼型设计与气动特性分析

对于大功率叶片的设计研究，翼型的选择和气动特性分析是重要内容。考虑到叶片的整体结构和气动特性，叶片在主要输功区域布置了风力机专用的DU系列翼型，叶尖部分布置NACA64系列翼型，能满足雷诺数较大，阻力小的要求。在确定截面翼型的基础上，结合Wilson理论设计模型初步算出翼型的轴向、周向诱导因子，叶尖、轮毂损失因子等具体...     

气动载荷下防热材料剥离颗粒输运特性的直接数值模拟研究

高超声速飞行器防热材料在气动载荷下发生机械剥蚀,进而影响绕流流态、气动性能、热载荷等,相关颗粒剥离动力学是高超声速热防护系统设计及防热材料体系评价中的共性基础性科学问题.研究通过近壁流动量纲分析,将烧蚀颗粒剥离过程模化为单个圆球惯性烧蚀颗粒在Couette流动中的动力学问题,并采用颗粒解析的直接数值模拟方法开展数值研究,获得了烧蚀颗粒关键特征参量对颗粒输运动力学的影响规律.研究发现,随着颗粒/流体密度比ρr越大,颗粒惯性St越大,则颗粒水平和法向输运速度均减小;随着颗粒粒径d p越大,颗粒惯性St越大,则颗粒水平输运速度减小,但是,法向输运速度和位移均因大颗粒受到更大的Saffman升力而增大.此外,烧蚀颗粒法向位移远小于水平位移,颗粒以水平输运为主.本研究最终建立了颗粒启动速度归一化表达式,发现归一化颗粒启动速度是颗粒和流体惯性的函数,即颗粒水平输运速度等于流体微团或中性浮力颗粒的速度减去惯性修正项.研究结果为烧蚀颗粒调制边界层作用机理研究提供支撑.     

垂直轴自由变偏角水轮机的多体耦合数值求解方法

基于达朗伯原理建立了自由变偏角水轮机结构动力学模型,并结合粘性CFD方法,形成了一种用于分析自由变偏角潮流能水轮机水动力特性的多体耦合数值计算方法.将该方法与试验及非耦合算法进行了对比分析,并研究了自启动、流速变化和负载变化对水轮机性能的影响规律.研究表明:本文建立的数值计算方法能够有效模拟自由变偏角水轮机在各种非定常过程中的运动特性和动力特性.     

Progresses in application of computational fluid dynamic methods to large scale wind turbine aerodynamics

The computational fluid dynamics (CFD) methods are applied to aerodynamic problems for large scale wind turbines. The progresses including the aerodynamic analyses of wind turbine profiles, numerical flow simulation of wind turbine blades, evaluation of aerodynamic performance, and multi-objective blade optimization are discussed. Based on the CFD methods, significant improvements are obtained to predict two/three dimensional aerodynamic characteristics of wind turbine airfoils and blades, and the vortical structure in their wake flows is accurately captured. Combining with a multi-objective genetic algorithm, a 1.5 MW NH-1500 optimized blade is designed with high efficiency in wind energy conversion.     

Status of large scale wind turbine technology development abroad

To facilitate the large scale (multi-megawatt) wind turbine development in China, the foreign efforts and achievements in the area are reviewed and summarized. Not only the popular horizontal axis wind turbines on-land but also the offshore wind turbines, vertical axis wind turbines, airborne wind turbines, and shroud wind turbines are discussed. The purpose of this review is to provide a comprehensive comment and assessment about the basic work principle, economic aspects, and environmental impacts of turbines.     

Analogy between a flapping wing and a wind turbine with a vertical axis of revolution

Based on an analysis of available experimental data, the hypothesis about an analogy between a flapping wing and a wind turbine of the Darrieus rotor type is justified. It is demonstrated that the torque on the shaft of the Darrieus rotor is generated by thrust forces acting on the blades in a pulsed flow. A conclusion is drawn that it is necessary to perform aerodynamic calculations of blades on the basis of the nonlinear theory of the wing in an unsteady flow with allowance for the airfoil thickness. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 152–155, March–April, 2009.     

漂浮运动对风力机气动特性的影响分析

海上浮式风机为蕴藏丰富的深海风能开发提供了有效的解决方案,浮式基础存在的大幅度纵荡、纵摇和艏摇运动可以改变风轮与流场间的相互作用,从而影响风力机的气动特性。基于叶素-动量理论及其修正方法,以NREL 5MW风力机为研究对象,考虑浮式基础运动对叶片不同径向位置处相对入流风速的影响,提出了风载荷的计算模型,通过编程计算获得了叶轮转矩和风力机功率,并比较了不同运动形式对风力机功率波动的影响。结果表明,纵摇对其功率特性影响最大,这为海上浮式风机的优化设计提供理论依据与数据基础。     

Implementation of all‐Mach Roe‐type schemes in fully implicit CFD solvers – demonstration for wind turbine flows

This paper presents the implementation of all‐Mach Roe‐type schemes in a fully implicit CFD solver. Simple 2D cases, such as the flow around inviscid and viscous aerofoils, were used for an initial validation of these methods, along with more challenging computations consisting of the 3D flow around the Model Experiments in Controlled Conditions wind turbine, in parked and rotating conditions. This work is motivated by the increased interest of the wind turbine industry in larger diameter wind turbines where compressibility effects near the blade tips may be important. Instead of using an incompressible flow solver, this paper explores the option of modifying an existing, efficient, compressible flow solver for use at lower Mach numbers. The good performance of the Roe solver and its popularity influenced the selection of schemes for this work. The results suggest that effective all‐Mach solutions are possible with implicit solvers, and the paper defines the implementation of the new fluxes and Jacobian, including an investigation of some numerical parameters, using as platform the Helicopter Multi‐Block solver of Liverpool University. Copyright © 2013 John Wiley & Sons, Ltd.     

Investigation of horizontal-axis wind turbine (HAWT) blade threedimensional rotational effect based on field experiments

Field experiments are performed on a two-bladed 33 kW horizontal-axis wind turbine (HAWT). The pressures are measured with 191 pressure sensors positioned around the surfaces of seven spanwise section airfoils on one of the two blades. Three-dimensional (3D) and two-dimensional (2D) numerical simulations are performed, respectively, on the rotor and the seven airfoils of the blade. The results are compared with the experimental results of the pressure distribution on the seven airfoils and the lift coefficients. The 3D rotational effect on the blade aerodynamic characteristics is then studied with a numerical approach. Finally, some conclusions are drawn as follows. From the tip to the root of the blade, the experimental differential pressure of the blade section airfoil increases at first and then decreases gradually. The calculated 3D result of the pressure distribution on the blade surface is closer to that of the experiment than the 2D result. The 3D rotational effect has a significant impact on the blade surface flow and the aerodynamic load, leading to an increase of the differential pressure on the airfoils and their lift coefficient than that with the 2D one because of the stall delay. The influence of the 3D rotational effect on the wind turbine blade especially takes place on the sections with flow separation.