Comparative turbulent three-dimensional Navier-Stokes hydrodynamic analysis and performance assessment of oscillating wings for renewable energy applications

Oscillating wings can extract energy from an oncoming water or air stream, and first large-scale marine demonstrators are being tested. Oscillating wing hydrodynamics is highly unsteady, may feature dynamic stall and leading edge vortex shedding, and is significantly three-dimensional due to finite-wing effects. Understanding the interaction of these phenomena is essential for maximizing power generation efficiency. Much of the knowledge on oscillating wing hydrodynamics stemmed from two-dimensional low-Reynolds number computational fluid dynamics studies and laboratory testing; real installations, however, will feature Reynolds numbers higher than 1 million and unavoidable finite-wing-induced losses. This study investigates the impact of flow three-dimensionality on the hydrodynamics and the efficiency of a realistic aspect ratio 10 device in a stream with Reynolds number of 1.5 million. The improvements achievable by using endplates to reduce finite-wing-induced losses are also analyzed. Three-dimensional time-dependent Navier-Stokes simulations using the shear stress transport turbulence model and a 30 million-cell grid are performed. Detailed comparative hydrodynamic analyses of the finite and the infinite wings reveal that flow three-dimensionality reduces the power generation efficiency of the finite wing with sharp tips and that with endplates by about 17% and 12% respectively. Presented analyses suggest approaches to further reducing these power losses

Massively parallel time- and frequency-domain Navier-Stokes Computational Fluid Dynamics analysis of wind turbine and oscillating wing unsteady flows

Increasing interest in renewable energy sources for electricity production complying with stricter environmental policies has greatly contributed to further optimisation of existing devices and the development of novel renewable energy generation systems. The research and development of these advanced systems is tightly bound to the use of reliable design methods, which enable accurate and efficient design. Reynolds-averaged Navier-Stokes Computational Fluid Dynamics is one of the design methods that may be used to accurately analyse complex flows past current and forthcoming renewable energy fluid machinery such as wind turbines and oscillating wings for marine power generation. The use of this simulation technology offers a deeper insight into the complex flow physics of renewable energy machines than the lower-fidelity methods widely used in industry. The complex flows past these devices, which are characterised by highly unsteady and, often, predominantly periodic behaviour, can significantly affect power production and structural loads. Therefore, such flows need to be accurately predicted. The research work presented in this thesis deals with the development of a novel, accurate, scalable, massively parallel CFD research code COSA for general fluid-based renewable energy applications. The research work also demonstrates the capabilities of newly developed solvers of COSA by investigating complex three-dimensional unsteady periodic flows past oscillating wings and horizontal-axis wind turbines. Oscillating wings for the extraction of energy from an oncoming water or air stream, feature highly unsteady hydrodynamics. The flow past oscillating wings may feature dynamic stall and leading edge vortex shedding, and is significantly three-dimensional due to finite-wing effects. Detailed understanding of these phenomena is essential for maximising the power generation efficiency. Most of the knowledge on oscillating wing hydrodynamics is based on two-dimensional low-Reynolds number computational fluid dynamics studies and experimental testing. However, real installations are expected to feature Reynolds numbers of the order of 1 million and strong finite-wing-induced losses. This research investigates the impact of finite wing effects on the hydrodynamics of a realistic aspect ratio 10 oscillating wing device in a stream with Reynolds number of 1.5 million, for two high-energy extraction operating regimes. The benefits of using endplates in order to reduce finite-wing-induced losses are also analyzed. Three-dimensional time-accurate Reynolds-averaged Navier-Stokes simulations using Menter's shear stress transport turbulence model and a 30-million-cell grid are performed. Detailed comparative hydrodynamic analyses of the finite and infinite wings highlight that the power generation efficiency of the finite wing with sharp tips for the considered high energy-extraction regimes decreases by up to 20 %, whereas the maximum power drop is 15 % at most when using the endplates. Horizontal-axis wind turbines may experience strong unsteady periodic flow regimes, such as those associated with the yawed wind condition. Reynolds-averaged Navier-Stokes CFD has been demonstrated to predict horizontal-axis wind turbine unsteady flows with accuracy suitable for reliable turbine design. The major drawback of conventional Reynolds-averaged Navier-Stokes CFD is its high computational cost. A time-step-independent time-domain simulation of horizontal-axis wind turbine periodic flows requires long runtimes, as several rotor revolutions have to be simulated before the periodic state is achieved. Runtimes can be significantly reduced by using the frequency-domain harmonic balance method for solving the unsteady Reynolds-averaged Navier-Stokes equations. This research has demonstrated that this promising technology can be efficiently used for the analyses of complex three-dimensional horizontal-axis wind turbine periodic flows, and has a vast potential for rapid wind turbine design. The three-dimensional simulations of the periodic flow past the blade of the NREL 5-MW baseline horizontal-axis wind turbine in yawed wind have been selected for the demonstration of the effectiveness of the developed technology. The comparative assessment is based on thorough parametric time-domain and harmonic balance analyses. Presented results highlight that horizontal-axis wind turbine periodic flows can be computed by the harmonic balance solver about fifty times more rapidly than by the conventional time-domain analysis, with accuracy comparable to that of the time-domain solver

Three-dimensional Navier-Stokes turbulent hydrodynamic analysis and performance assessment of oscillating wings for power generation

A wing simultaneously heaving and pitching can extract energy from an oncoming water or air stream. First large-scale commercial demonstrators are being installed and tested. The operating conditions of this device is likely to feature Reynolds numbers in excess of 500,000. Strong finite wing effects spoiling the power generation efficiency are also expected. This paper thoroughly investigates the hydrodynamics of oscillating wings at a Reynolds number of 1,500,000 considering finite wing effects for an aspect ratio 10 wing with either sharp tips or endplates to reduce tip vortex losses. The study of these periodic flows uses three-dimensional time-dependent Navier-Stokes simulations with grids featuring more than 30 million cells The shear stress turbulence model of Menter is used for the turbulence closure. Main contributions include: a) the quantification of the efficiency improvement achievable by using wings with endplates rather than bare tips, and b) detailed comparisons of the wing hydrodynamics with and without endplates, and the infinite wing

Comparative assessment of the harmonic balance Navier Stokes technology for horizontal and vertical axis wind turbine aerodynamics

Several important wind turbine unsteady flow regimes, such as those associated with the yawed wind condition of horizontal axis machines, and most operating conditions of all vertical axis machines, are predominantly periodic. The harmonic balance Reynolds-averaged Navier-Stokes technology for the rapid calculation of nonlinear periodic flow fields has been successfully used to greatly reduce runtimes of turbomachinery periodic flow analyses in the past fifteen years. This paper presents an objective comparative study of the performance and solution accuracy of this technology for aerodynamic analysis and design applications of horizontal and vertical axis wind turbines. The considered use cases are the periodic flow past the blade section of a utility-scale horizontal axis wind turbine rotor in yawed wind, and the periodic flow of a H-Darrieus rotor section working at a tip-speed ratio close to that of maximum power. The aforementioned comparative assessment is based on thorough parametric time-domain and harmonic balance analyses of both use cases. The paper also reports the main mathematical and numerical features of a new turbulent harmonic balance Navier-Stokes solver using Menter’s shear stress transport model for the turbulence closure. Presented results indicate that a) typical multi-megawatt horizontal axis wind turbine periodic flows can be computed by the harmonic balance solver about ten times more rapidly than by the conventional time-domain analysis, achieving the same temporal accuracy of the latter method, and b) the harmonic balance acceleration for Darrieus rotor unsteady flow analysis is lower than for horizontal axis machines, and the harmonic balance solutions feature undesired oscillations caused by the wide harmonic content and the high-level of stall predisposition of this flow field type

Load balance and Parallel I/O: Optimising COSA for large simulations

This paper presents the optimisation of the parallel functionalities of the Navier-Stokes Computational Fluid Dynamics research code COSA, a finite volume structured multi-block code featuring a steady solver, a general purpose time-domain solver, and a frequency-domain harmonic balance solver for the rapid solution of unsteady periodic flows. The optimisation focuses on improving the scalability of the parallel input/output functionalities of the code and developing an effective and user-friendly load balancing approach. Both features are paramount for using COSA efficiently for large-scale production simulations using tens of thousands of computational cores. The efficiency enhancements resulting from optimising the parallel I/O functionality and addressing load balance issues has provided up to a 4x performance improvement for unbalanced simulations, and 2x performance improvements for balanced simulations

Cross-comparative analysis of loads and power of pitching floating offshore wind turbine rotors using frequency-domain Navier-Stokes CFD and blade element momentum theory

Reliable predictions of the aero- and hydrodynamic loads acting on floating offshore wind turbines are paramount for assessing fatigue life, designing load and power control systems, and ensuring the overall system stability at all operating conditions. However, significant uncertainty affecting both predictions still exists. This study presents a cross-comparative analysis of the predictions of the aerodynamic loads and power of floating wind turbine rotors using a validated frequency-domain Navier-Stokes Computational Fluid Dynamics solver, and a state-of-the-art Blade Element Momentum theory code. The considered test case is the National Renewable Energy Laboratory 5 MW turbine, assumed to be mounted on a semi-submersible platform. The rotor load and power response at different pitching regimes is assessed and compared using both the high- and low-fidelity methods. The overall qualitative agreement of the two prediction sets is found to be excellent in all cases. At a quantitative level, the high- and low-fidelity predictions of both the mean rotor thrust and the blade out-of-plane bending moments differ by about 1 percent, whereas those of the mean rotor power differ by about 6 percent. Part of these differences at high pitching amplitude appear to depend on differences in dynamic stall predictions of the approaches

Harmonic balance Navier–Stokes aerodynamic analysis of horizontal axis wind turbines in yawed wind

Multi–megawatt horizontal axis wind turbines often operate in yawed wind transients in which the resulting periodic loads acting on blades, drive–train, tower and foundation adversely impact on fatigue life. Accurately predicting yawed wind turbine aerodynamics and resulting structural loads can be challenging, and would require the use of computationally expensive high–fidelity unsteady Navier–Stokes Computational Fluid Dynamics. The high computational cost of this approach can be significantly reduced by using a frequency–domain framework. The paper summarizes the main features of the COSA Harmonic Balance Navier–Stokes solver for the analysis of open rotor periodic flows, presents initial validation results based on the analysis of the NREL Phase VI experiment, and it also provides a sample application to the analysis of a multi–megawatt turbine in yawed wind. The reported analyses indicate that the harmonic balance solver determines the considered periodic flows from 30 to 50 times faster than the conventional time–domain approach with negligible accuracy penalty to the latter. Copyright c2017 John Wiley & Sons, Ltd

Darrieus wind turbine blade unsteady aerodynamics:a three-dimensional Navier-Stokes CFD assessment

Energized by the recent rapid progress in high-performance computing and the growing availability of large computational resources, computational fluid dynamics (CFD) is offering a cost-effective, versatile and accurate means to improve the understanding of the unsteady aerodynamics of Darrieus wind turbines, increase their efficiency and delivering more cost-effective and structurally sound designs. In this study, a Navier-Stokes CFD research code featuring a very high parallel efficiency was used to thoroughly investigate the three-dimensional unsteady aerodynamics of a Darrieus rotor blade. Highly spatially and temporally resolved unsteady simulations were carried out using more than 16,000 processor cores on an IBM BG/Q cluster. The study aims at providing a detailed description and quantification of the main three-dimensional effects associated with the periodic motion of this turbine type, including tip losses, dynamic stall, vortex propagation and blade/wake interaction. Presented results reveal that the three-dimensional flow effects affecting Darrieus rotor blades are significantly more complex than assumed by the lower-fidelity models often used for design applications, and strongly vary during the rotor revolution. A comparison of the CFD integral estimates and the results of a blade-element momentum code is also presented to highlight strengths and weaknesses of low-fidelity codes for Darrieus turbine design. The reported CFD results may provide a valuable and reliable benchmark for the calibration of lower-fidelity models, which are still key to industrial design due to their very high execution speed

Harmonic balance Navier-Stokes analysis of tidal stream turbine wave loads

ARCTIC, a novel incompressible Reynolds–averaged Navier–Stokes finite volume code for the hydrodynamic analysis of open rotor unsteady loads is presented. One of its unique features is a harmonic balance solver enabling high–fidelity analyses of turbine periodic hydrodynamic loads with runtimes reduced by more than one order of magnitude over conventional time–domain CFD, and with negligible accuracy penalty. The strength of the new technology is demonstrated by analyzing with both harmonic balance and time–domain solvers the load fluctuations of a realistic tidal stream turbine. Such fluctuations are caused by a harmonic perturbation of the freestream velocity similar to that due to surface gravity waves