带偏航角串列式两风机复杂尾流场数值模拟

在多风机风电场中,通过主动调节上游风机的偏航角度,抑制上游风机尾流对下游风机的影响,减少风力机机组之间的尾流相互干扰,以达到提高整个风电场效率的目的。采用基于开源平台Open FOAM自主开发的FOWT-UALM-SJTU求解器中风电场求解模块ALMWindFarmFoam,将致动线模型与CFD方法相结合,利用大涡模拟(LES)计算研究当上游风机处于不同偏航角度时,两风机之间的复杂尾流干扰效应。对比分析偏航角度改变时,上下游风机气动功率的输出特性,尾流速度变化以及风机的尾涡结构。数值模拟结果表明:在上下游风机沿流向方向距离保持不变的情况下,随着上游风机偏航角度的变化,上下游风机的尾流干扰现象以及下游风机的入流条件会发生明显改变,并会对下游风机的气动功率输出以及两风机风电场的整体流场产生显著影响。     

浮式风机气动—水动—气弹性耦合响应数值模拟

随着海上风电产业的快速发展,大型浮式风机逐渐从概念设计走向工程应用,但仍面临较大的挑战。一方面,在风、浪等环境载荷的作用下,浮式风机的气动载荷和水动力响应之间存在明显的相互干扰作用;另一方面,风力机大型化使得叶片细、长、薄的特点愈发突出,叶片柔性变形十分显著,这会影响到浮式风机的耦合性能。基于两相流CFD求解器naoe-FOAM-SJTU,结合弹性致动线模型和等效梁理论,建立了浮式风机气动—水动—气弹性耦合响应计算模型,并对规则波和剪切风作用下Spar型浮式风机的气动—水动—气弹性耦合响应进行了数值模拟分析。结果表明,风力机气动载荷使得叶片挥舞变形十分显著,而叶片的扭转变形会明显降低风力机的气动载荷。此外,风力机气动载荷会增大浮式平台的纵荡位移和纵摇角,同时,浮式平台运动响应会导致风力机气动载荷产生大幅度周期性变化。进一步地,叶片结构变形响应会使得浮式风机尾流场的速度损失和湍动能有所降低。     

基于非定常面元法的海上浮式风机气动性能研究

对于海上浮式风机而言,由于受到剪切风、塔影效应、浮式基础运动等因素的共同影响,其气动载荷会更加复杂,因此如何准确快速地对海上风力机的气动性能进行预估显得尤为重要。基于速度势的非定常面元法理论,研究海上浮式风机气动载荷特性,编制了相关的计算程序。以NREL 5 MW风机为例,建立了叶片和尾流的三维数值模型,计算得到了不同风速下风机的输出功率以及叶片表面的压力分布,对比数据结果分析了该方法的可靠性。针对非定常流动,模拟了剪切风和塔影效应的作用,并重点分析了浮式基础运动对风机气动载荷的影响。研究表明,浮式基础的纵荡和纵摇会增加输出功率的波动幅值,艏摇运动会导致单个叶片上的气动载荷产生较大的波动,为浮式风机叶片控制提供了参考。     

布设间距对人工鱼礁流场效应影响的数值模拟

人工鱼礁在海中不同的布局方式会产生不同的流场效应.通过水槽与风洞模型实验很难对组合鱼礁流场全局进行观测,为此提出1种基于CFD原理的人工鱼礁流场特性的计算机数值模拟方法,模拟了不同流速条件下布设间距对方形组合礁体流场效应的影响.结果表明:数值模拟较好地反映了人工鱼礁周围上升流和背涡流的分布情况;当两礁体布设间距为1.5倍礁体尺寸时,所产生的上升流高度达到最大值;当布设间距为1倍/时,其上升流的影响面积为最大;当布设间距为1.5倍时所产生的背涡流效果最好;两礁体之间的最佳布设间距应为礁体尺寸的1～1.5倍.数值模拟结果与风洞实验结果基本相符,该结论可为人工鱼礁实际投放提供理论参考.     

五圆柱体结构群绕流特性及互扰效应研究

基于Fluent流体计算平台,运用大涡模拟方法对亚临界雷诺数Re=3900下“X”形排列五圆柱体结构群三维绕流特性进行研究,主要分析来流攻角α与间距比L/D两个关键参数对五圆柱体结构群的尾流区三维涡结构演化与流体力系数的影响,并揭示其内在流动互扰机理。研究表明:来流攻角和间距比的变化对五圆柱体结构群流动控制及互扰效应的影响显著。在小间距比工况下,观察到柱体群间隙区域内流体高速流动的现象,导致五圆柱体之间的互扰作用十分强烈。间隙流对中间圆柱体和下游圆柱体有较强的冲击作用,对其表面的流体力分布特性有显著的影响。另外,大间距比工况下,当α=0°与L/D≥5.0工况时,柱体群尾流效应强于其间隙流效应。当α=22.5°与L/D=7.0时,位于下游与中间处的圆柱体流体绕流特性存在较大差异。而当α=45°与L/D≥6.0时,位于上游与中间处的圆柱体尾流区均会产生正负交替的漩涡结构。     

近海海上风电场水下噪声传播模型适用性研究

通过现场采集近海海上风电场工程区运营期风机水下噪声和背景噪声数据,计算了噪声信号的倍频带声压级,功率谱级和峰值声压级,确定了海上风电场水下噪声总声源级为148.3 d B,以此开展近海海上风电工程风机水下噪声频域特性、功率密度谱特性等研究。在此基础上使用Kraken简正波模型和Bellhop射线模型对风电场运营期风机水下噪声在水平与垂直方向上的传播进行模拟,模拟了噪声在不同频带内的衰减程度,结果显示模型模拟结果在不同频率下的衰减趋势有着很大差异,产生了明显的多途干涉现象,通过实测数据对建立的噪声传播模型进行验证,发现Kraken简正波模型在500 Hz以下,Bellhop射线模型在500 Hz以上适合模拟实际水下噪声传播情形,同时海区本身背景噪声的存在会对预测的准确性产生影响。这些结论可用于进一步对近海海上风电场水下噪声传播的研究。     

偏航工况下潮流能水轮机性能及尾流场特性分析

潮流能发电机组在工作中,水轮机与潮流流向之间的偏航角对于机组的性能有较大影响。基于ANSYS Fluent软件采用SST k-ω湍流模型和滑移网格旋转模型对不同偏航角下水轮机进行数值模拟,分析偏航工况下水轮机性能及尾流场特性,并进行试验验证。结果表明:随着偏航角的增大,轴向力系数与功率系数降低,尾流场流速变化不对称现象趋于明显,同时造成周围流场湍流强度增大。     

潮流能水轮机尾流场及涡特性DDES模拟

为研究潮流能水轮机尾流场流动特性及涡结构组成,基于DDES(Delayed Detached Eddy Simulation)模型对不同流速和转速的4种工况下水轮机尾流场进行数值模拟,并进一步探究其尾流场空间涡结构的变化特性。结果表明,水轮机的数值模拟结果与试验结果能够较好吻合。对比不同工况下的尾流场模拟结果可知:水轮机尾流区域流动复杂,延迟分离涡模拟方法能有效模拟水轮机旋转过程中产生的叶尖涡、叶尖脱落涡、轮毂涡等不同涡结构,并能完整观察到叶尖涡的产生、脱落、失稳、破碎过程。转速一定时,流速越大,叶尖脱落涡、轮毂涡的发展距离越远;流速一定时,转速越大,涡的发展距离越短。本文数值模拟计算结果可为实际海况中潮流能阵列水轮机的布局提供可靠依据。     

湍流对潮流能发电水轮机性能影响数值仿真研究

潮流能发电水轮机的实际工作海域往往存在不同程度的湍流,而湍流会对潮流能发电水轮机的获能系数、轴向力系数和尾流场性能等产生影响。研究湍流对潮流能发电水轮机性能的影响规律,对于实海况下潮流能发电水轮机的性能预测、可靠性和安全性的提高以及潮流能发电场多机组排布优化等具有一定的参考价值。通过对潮流能发电水轮机试验模型进行建模,并用CFD(Computational Fluid Dynamics)分析软件Fluent对处于不同湍流强度下的潮流能发电水轮机性能进行数值模拟,得到其获能系数、轴向力系数及尾流场特性。通过分析数值模拟结果,并与相关参考文献的试验结果进行对比。研究结果表明:湍流强度越大,水轮机获能系数和轴向力系数越小,尾流场速度恢复越快;水轮机后方尾流场纵向和横向影响区域更大。     

叶片变桨失效过程中空气动力失衡对海上风机影响

叶片桨距角之间的角度差异产生的空气动力失衡是海上风机的主要动力问题之一。基于海上风机分析程序FAST和水动力计算程序WADAM开发的一种时域数值模拟程序,可计算海上风机系统在风浪载荷作用下的耦合动力响应。应用此数值工具,模拟一个叶片上变桨控制系统失效的情况,研究空气动力失衡对浮式海上风机系统运动响应的影响。分析表明,空气动力载荷失衡引起的激振不仅激发了浮式基础的横向运动,而且增大了基础的纵荡运动和首摇运动。同时,空气动力失衡还大幅增加了风机塔柱底部受到的横向剪切力,对风机系统的安全性造成了威胁。     

基于尾流模型的风场偏航控制优化研究

在风机大尺度化与风场大型化的趋势下,如何通过合适的控制策略以降低尾流损失成为关键问题。以包括30台NREL-5MW风机并采用5行6列平行四边形布置方式的小型风场为研究对象,基于显示尾流模型,以各风机偏航角度为优化参数,风场总功率为目标函数,使用粒子群优化算法对比分析了偏航控制对不同风速、风向、湍流强度下的风场性能提升效果。结果表明,偏航控制优化可在风向与风机行或列方向平行时发挥明显效果,当风机行列间距为4倍风轮直径且湍流强度为5%时,在不同风速下偏航控制可分别将风场总体发电量提升15%~20%,但对于布置间距大于7倍风轮直径或湍流强度高于15%时的风场,其作用十分有限,总体发电量提升在5%以内。     

Modeling of twin-turbine systems with vertical axis tidal current turbines: Part I—Power output

Recent interest in the tidal current industry has driven development of the prototype from the stand-alone turbine to the twin-turbine system. In this paper, we develop a numerical model to systematically analyze the relationship between the power output and the configuration of a twin-turbine system. First, we present the design principle of the twin-turbine system. We then develop the numerical model for simulating the operation of the system, and validate the model by conducting towing tank experimental tests. We then use the model to predict the power output of the system. The results of this study show that the total power output of a twin-turbine system with optimal layout can be about 25% higher than two times that of a stand-alone turbine. We also discuss the hydrodynamic interaction between the two turbines under different configurations of the system. We conclude that the optimally configured counter-rotating system should be a side-by-side system, and that the optimally configured co-rotating system should have the downstream turbine partially in the wake of the upstream turbine, depending on the detailed configuration of the turbine.     

Modeling of twin-turbine systems with vertical axis tidal current turbine: Part II—torque fluctuation

We recently showed the advantage of using a numerical system to extract energy from tidal currents by developing a new twin-turbine model (Li and Calisal, 2010a). Encouraged by this result, we decided to use this model to study another important characteristic of the turbine system, torque fluctuation. This effort is summarized in this paper. The torque fluctuation is expected to reduce the fatigue life of tidal current turbines, though potentially it also may deteriorate the power quality of tidal current turbines. In this paper, after reviewing the twin-turbine model, we use it to predict the torque fluctuation of the system with the same configurations as we used to study the power output in Li and Calisal (2010a). Specifically, we investigate the torque fluctuation of twin-turbine systems with various turbine parameters (e.g., relative distance between two turbines and incoming flow angle) and operational condition (e.g., tip speed ratio). The results suggest that the torque of an optimally configured twin-turbine system fluctuates much less than that of the corresponding stand-alone turbine, under the same operating conditions. We then extensively compare the hydrodynamic interaction’s impact on the torque fluctuation and the power output of the system. We conclude that the hydrodynamic interactions pose more constructive impacts on the torque fluctuation than on the power output. The findings indicate that the optimally configured counter-rotating system should be a side-by-side system, and that the optimally configured co-rotating system should have the downstream turbine partially in the wake of the upstream turbine depending on the detailed configuration of the turbines. Furthermore, one must balance the optimal torque fluctuation against the optimal power output.     

Tidal Turbine Array Optimization Based on the Discrete Particle Swarm Algorithm

In consideration of the resource wasted by unreasonable layout scheme of tidal current turbines, which would influence the ratio of cost and power output, particle swarm optimization algorithm is introduced and improved in the paper. In order to solve the problem of optimal array of tidal turbines, the discrete particle swarm optimization (DPSO) algorithm has been performed by re-defining the updating strategies of particles’ velocity and position. This paper analyzes the optimization problem of micrositing of tidal current turbines by adjusting each turbine’s position, where the maximum value of total electric power is obtained at the maximum speed in the flood tide and ebb tide. Firstly, the best installed turbine number is generated by maximizing the output energy in the given tidal farm by the Farm/Flux and empirical method. Secondly, considering the wake effect, the reasonable distance between turbines, and the tidal velocities influencing factors in the tidal farm, Jensen wake model and elliptic distribution model are selected for the turbines’ total generating capacity calculation at the maximum speed in the flood tide and ebb tide. Finally, the total generating capacity, regarded as objective function, is calculated in the final simulation, thus the DPSO could guide the individuals to the feasible area and optimal position. The results have been concluded that the optimization algorithm, which increased 6.19% more recourse output than experience method, can be thought as a good tool for engineering design of tidal energy demonstration.     

0903号台风对我国海上风力发电的影响

基于CCMP(cross-calibrated multi-platform)再分析风场数据,研究了0903号台风经过前后21天对我国东部5个海上风电场的影响。发现:风机的输出功率受其所处台风场的相对位置影响,台风北进过程中前沿(后沿)在2009年6月18—19(27—28)日使得风机输出功率显著增加,而6月22—23日,台风中心进入东海时, 5个站位风机输出功率最低。唐山市海港区S1站单机运营能力最低(64.3%)且额定容量最低(0);嘉兴市平湖县S4站单机运营能力最高(88.1%);南通市如东县S3站单机额定容量最大(13.1%)。0903号台风过境整个过程中,风速未超过风机的切出风速,对我国东部海洋风力发电场,提高了发电量。     

A three-dimensional wake impingement model and applications on tandem oscillating foils

Potential flow based vortex numerical methods have been widely used in aerodynamics and hydrodynamics. In these methods, vortices shed from lifting bodies are traced by using vortex filaments or dipole panels. When the wake elements encounter a downstream body, such as a rudder behind a propeller or a stator behind a rotor, a treatment is necessary to divert the wake elements to pass by the body. This treatment is vital to make wake simulations realistic and to satisfy the non-penetration condition during wake body interaction. It also helps to avoid pure numerical disturbances such as when a vortex filament or an edge of a dipole panel passes through the collection point of a body element; this is a singularity for induced velocity and it will introduce a large numerical disturbance. This necessary treatment for three-dimensional problems with geometrical complexity has not been found to date. In this study, a wake impingement model was developed to divert wake elements to slip over the body surface, model the vortex/body interaction, and predict forces on fluctuating components. The model was also tested on configurations of oscillating foils in tandem with an existing panel method code. Simulation results with the wake impingement model are shown to be in closer agreement with limited published experimental data than those without the model. With the established wake impingement model, force fluctuations on the after body due to the wake vortex impingement were investigated based on a series of simulations. The series varied several parameters including distance between two foils, oscillating frequency, span, rear foil pitch angle, swap angle and vertical position.     

Experimental analysis of the flow field around horizontal axis tidal turbines by use of scale mesh disk rotor simulators

Understanding the flow field around horizontal axis marine current turbines is important if this new energy generation technology is to advance. The aim of this work is to identify and provide an understanding of the principal parameters that govern the downstream wake structure and its recovery to the free-stream velocity profile. This will allow large farms or arrays of devices to be installed whilst maximising device and array efficiency. Wake characteristics of small-scale mesh disk rotor simulators have been measured in a 21 m tilting flume at the University of Southampton. The results indicate that wake velocities are reduced in the near wake region (close behind the rotor disk) for increasing levels of disk thrust. Further downstream all normalised wake velocity values converge, enforcing that, as for wind turbines, far wake recovery is a function of the ambient flow turbulence. Varying the disk proximity to the water surface/bed introduces differential mass flow rates above and below the rotor disk that can cause the wake to persist much further downstream. Finally, the introduction of increased sea bed roughness whilst increasing the depth-averaged ambient turbulence actually decreases downstream wake velocities. Results presented demonstrate that there are a number of interdependent variables that affect the rate of wake recovery and will have a significant impact on the spacing of marine current turbines within an array.