Next Article in Journal
Numerical Study on Lateral Response of Offshore Monopile in Sand under Local Scouring Conditions
Previous Article in Journal
Exploring Benthic Scavenger Assemblages, a Multi-Habitat Approach in NW Iberian Peninsula
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Load Characteristics and Extreme Response of Straight-Bladed Floating VAWT Using a Fully Coupled Model

1
Key Laboratory of High Performance Ship Technology, Wuhan University of Technology, Ministry of Education, Wuhan 430063, China
2
Department of Naval Architecture and Ocean Engineering, School of Naval Architecture, Ocean and Energy Power Engineering, Wuhan University of Technology, Wuhan 430063, China
3
China Ship Design and Research Center, Wuhan 430060, China
4
Green & Smart River-Sea-Going Ship, Cruise and Yacht Research Centre, Wuhan University of Technology, Wuhan 430063, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(1), 185; https://doi.org/10.3390/jmse11010185
Submission received: 25 November 2022 / Revised: 29 December 2022 / Accepted: 30 December 2022 / Published: 11 January 2023
(This article belongs to the Topic Advances in Wind Energy Technology)

Abstract

:
Operating Offshore Floating Vertical Axis Wind Turbines (OF-VAWT) have the potential to perform well in the deep-sea area. Some researchers gave performance prediction by developing simplified computing models. However, these models have imperfections in considering load and motion nonlinearity, especially in extreme environments. In this work, a numerical model is developed composed of Computational Fluid Dynamics (CFD) and Dynamic Fluid Body Interaction (DFBI) to acquire the aero-hydrodynamic load and performance of OF-VAWT in general and extreme environments. Unsteady Reynolds-Averaged Navier-Stokes (URANS), SST k-ω and Eulerian Multi-Phase (EMP) models are combined to generate a gas-liquid two-phase flow field; the Volume of Fluid (VOF) model is employed to capture free-surface and make numerical wind-wave. DFBI superposition motion technology is proposed for local motion definition and motion solution, and overset with sliding meshes is introduced to achieve the grid motion. The numerical approach is verified by the tunnel and tank experimental data from the available literature. Simulation results of general cases, such as variable wind speed, wave height and wave length, are compared to discuss the effect of environmental parameters on load and performance. Comparison shows that this straight-bladed OF-VAWT is more susceptible to wind speed. Furthermore, the aerodynamic load generated by the shut-down rotor is still significant in extreme environment, which has implications for the development of OF-VAWT controller.

1. Introduction

Wind energy exploration is a popular research topic in recent years, as it can be harnessed without carbon emission to handle environmental problems. Recently, wind power development has turned its sight from onshore to offshore [1], and coastal fixed-type wind turbines are the beginning. Furthermore, it is well-known that most of the wind energy is far from the coast, so utilizing energy in deep-sea areas interests various scholars and countries today.
Floating Offshore Wind Turbines (FOWT) in which the buoyancy is provided by a floating foundation and the power is generated by a rotor are generally classified as Vertical Axis Wind Turbines (VAWT) and Horizontal Axis Wind Turbines (HAWT) according to the direction of the rotating axis [2]. In addition, the semi-submersible platform with a larger platform is proved to be a better choice for offshore wind turbines in Offshore Code Comparison Collaboration Continuation (OC4) [3]. Further, it is also known the VAWT could capture wind energy in all directions without a steering device, and its nacelle and gearbox are closer to the bottom. Accordingly, the combination of a semi-submersible platform and VAWT is expected to perform well in the deep-sea environment to avoid common damage [4].
Operating Offshore Floating Vertical Axis Wind Turbines (OF-VAWT) can be greatly influenced by the environment. Attentionally, enormous wind turbines are susceptible to being damaged by wind loads. In 2008, five coastal land-based wind turbines in Taiwan were destroyed by super typhoon ‘Jangmi’ [5], and more than 18 wind turbines were damaged during typhoon ‘Weimaxun’ in 2014 [6]. Additionally, numerous offshore floating structures are capsized or damaged by wave load. In 2013, ‘Prestige’ was damaged during an encounter with a large wave [7], and in 2022 a large ship named ‘ALSALMY6’ capsized due to heavy weather conditions in the Persian Gulf. Therefore, it is essential to investigate the load and motion response of OF-VAWT in both general and extreme environments for engineering design implications.
Operating offshore wind turbines generally suffer from aero-hydrodynamic coupling in wind-wave environments. Recently, some scholars made efforts on predicting the performance of several FOWTs considering aero-hydrodynamic coupling, the properties including rotor, platform and simulation module are shown in Table 1. Obviously, floating HAWTs are the primary objects in current research, and OF-VAWTs are rarely investigated especially for the straight-bladed type. In addition, the semi-submersible platform DeepCWind is mostly employed as the floating foundation, and CFD technology is welcomed in solving aero-hydrodynamics simultaneously by more scholars. Some scholars also studied the coupling interaction by other simplified means. For instance, the floating foundation is assumed as a forced motion with regular function, through this way, the effect of pitch and surge motion on the aerodynamic performance of FOWT is investigated [8,9,10]. In a word, the effect of aero-hydrodynamic coupling on performance of straight-bladed OF-VAWT is rarely investigated especially in an extreme environment.
Methodology for investigating the aero-hydrodynamics of OF-VAWTs includes model experiments and numerical simulations [20]. The experimental method is generally used for validation considering its high cost and restricted measurement, thus numerical simulation is more popular in analyzing details and conducting numerous cases. There are several computing codes developed for OF-VAWT in different research institutes. For instance, Tianjin University developed an aero-hydrodynamic code for Darrieus-spar type based on double-multiple-stream-tubes and potential theory [21], as well as the Simo-Riflex-DMS code presented by the Norwegian University of Science and Technology [22] and the Hawc2 code presented by the Technical University of Denmark [23]. However, since too many assumptions are introduced into Blade Element Momentum (BEM), Vortex model and potential theory, there are imperfections in considering nonlinear loads, such as dynamic stall, green water and slamming [24]. Therefore, it is essential to analyze the OF-VAWT using Computational Fluid Dynamics (CFD) that can meet the requirements of accuracy in both aerodynamics and hydrodynamics.
The effect of dreadful environments on FOWT has been investigated by scholars in several ways. However, the prediction of OF-VAWT, especially for the straight-bladed type, is rarely conducted, and most related research has imperfections in consideration of the load and motion nonlinearity due to the adopted simplified computing model. Therefore, the present study aims to investigate the aero-hydrodynamic load and performance characteristics of straight-bladed OF-VAWT in general and extreme environments. To consider the nonlinear aero-hydrodynamics simultaneously, the adopted CFD and Dynamic Fluid Body Interaction (DFBI) superposition motion technology are described in Section 2. We present the validation of this numerical approach, which was carried out by a VAWT wind tunnel experiment and a tank model experiment of free-decay motion, in Section 3. The details of straight-bladed OF-VAWT system, including geometry, properties and numerical model, are present in Section 4. In Section 5, we present several typical cases to analyze the characteristics of aero-hydrodynamic performance and load in time-domain. The differences in load and performance statistics between different environment cases, as well as implications for engineering design, are discussed in Section 6.

2. Numerical Methodology

It is well known that FOWT has complicated aero-hydrodynamic coupling load and multiple motions, thus, simulating the body motion and fluid load with high-fidelity are two keys. Further, the aforementioned CFD technology has advantages in solving nonlinear loads for both aerodynamics and hydrodynamics, and the mooring line could be performed by a quasi-static catenary model. Accordingly, the viscous-flow CFD model, quasi-static catenary model and DFBI superposition motion are combined to predict the aero-hydrodynamic load and performance of OF-VAWT in this paper.

2.1. Computational Fluid Dynamics

Various assumptions are added in traditional theory to simplify problems so that load nonlinearity of wind and wave is not considered. We employed aero-hydrodynamic solvers based on viscous-flow theory, as well as the finite volume method (FVM). The spatial fluid domain is discretized into small-scale control volumes with finite numbers, assuming that each control volume is governed by continuity and momentum equations. However, parameters such as velocity and pressure changed randomly over time, the fluid flow will be more complicated and difficult to solve. On account of accuracy and cost, Unsteady Reynolds-Averaged Navier Stokes (URANS) is mostly used to predict turbulent motion; the idea of RANS method is to average Navier-Stokes equations with time, decompose velocity and pressure into mean and perturbation; the RANS equation is described as follows:
ρ [ u i ¯ t + ( u i u j ) ¯ x j ] = p ¯ x i + μ 2 u i + x j ( ρ u i u j ¯ )
where u is velocity, ρ is density of fluid, p is press, μ is molecular viscosity, u ¯ is average speed of time, p ¯ is average press of time, ( ρ u i u j ¯ ) is an item of Reynolds stress which is a second-order tensor. As new unknowns bring in, a turbulence model should be used to close the equations. Particularly, the approach combining the RANS equation and k-ω turbulence model has been validated in calculating aerodynamics and hydrodynamics [8,25]. Transport equations of turbulent kinetic energy and dissipation rate for SST k-ω model are described as follows:
t ( ρ k ) + x i ( ρ k u i ) = x j ( ( μ + μ t σ k ) k x j ) + G k ˜ Y k + S k
t ( ρ ω ) + x i ( ρ ω u i ) = x j ( ( μ + μ t σ ω ) ω x j ) + G ω Y ω + D ω + S ω
where G k ˜ is a term of turbulent kinetic energy generation for k, G ω is a term of dissipation rate generation, D ω is a term of cross diffusion for ω , Y k is a term of turbulent kinetic energy dissipation for k, Y ω is a term of turbulent kinetic energy dissipation for ω , S k and S ω are user-defined terms. Then, the basic solution of the flow field can be obtained by combining Equations (1)–(3). To capture the gas-liquid two-phase interface precisely, Volume of Fluid (VOF) is adopted in which a phase volume fraction is used to ascertain phase distribution and location of interface, which is defined as follows:
a i = V i V
where V i is volume of phase i in the mesh element, V is volume of mesh element. The sum of the volume fractions of all phases in an element must be one. According to the value of volume fraction, the existence of different phases or fluids in grid cells can be distinguished. a i = 0 means that the corresponding mesh element has no phase at all, a i = 1 means that the corresponding mesh element is completely phase filled, 0 < a i < 1 means that there are interfaces of multiple phases.

2.2. DFBI Superposition Motion

Different from the conventional offshore structures, the OF-VAWT has 7-DOF motion, including global 6-DOF free motion and local 1-DOF rotational motion. Attentionally, this local rotate motion is a superposition motion based on global motion. However, the conventional motion equation can only contain 6-DOF, thus this paper intends to decouple the composite motions by defining the relationship between the global coordinate system and the rotating coordinate system. First, the rotor and the platform are integrated as one body to solve the 6-DOF motion. The overall structure generates three translational motions and three rotational motions based on the global coordinate system, and the equations are shown as follows.
m d v d t = F
I d ω d t + ω × I ω = M
where m is the mass of body, F is the resultant force acting on the body, and v is the velocity. Furthermore, ω is the angular velocity, M is the resultant bending moment, and I is the tensor of the moments of inertia. Gravity and fluid forces are the main external effects except for special conditions. The fluid force F f l u i d , moment M f l u i d and moment of inertia I are calculated by the following formula:
I = [ I x x I x y I x z I x y I y y I y z I x z I y z I z z ]
F f l u i d = i = 1 A p i n i + i = 1 A τ i v i
M f l u i d = i = 1 A l i p i n i + i = 1 A l i τ i v i
where I i i represents principal components along the diagonal and I i j represents off-diagonal components. P i and τ i are, respectively, the pressure and shear stress acting on face i, n i and v i are the area vector of face, l i is the distance vector from the body center of mass to the center of face, A represents the entire body surface area. Figure 1 illustrates the relationship between the local coordinate system and the global coordinate system, so the translation and rotation of the rotating coordinate system can be formulated as follows:
( d θ l x d θ l y d θ l z ) = ( d θ g x d θ g y d θ g z )
( d X l d Y l d Z l ) = ( d X g + L sin d θ y d Y g + L cos d θ y d Z g + L ( 2 cos d θ x cos d θ y ) )
where dX, dY, and dZ represent the tiny translation distance in the three directions, and i represents the tiny rotation angle about the three main axes, subscript L represents local coordinate system (rotation system), and subscript g means global coordinate system. Finally, the Z-axis of the local coordinate system is set as the rotation axis with constant angular speed.

2.3. Mooring Model

The mooring lines in FOWT systems, subjected to their own weight, are hanging between two points (anchors and fairleads). To consider the effect of a mooring system on hydrodynamics and motion, a quasi-static catenary model is employed. Figure 2 illustrates the forces and curve. The shape can be defined as follows:
x = c λ 0 g u + c 2 D L e q λ 0 g sinh ( u ) + α
y = c λ 0 g cosh ( u ) + c 2 2 D L e q λ 0 g sinh 2 ( u ) + β
c = λ 0 L e q g sinh ( u 2 ) sinh ( u 1 )
where g is the gravity acceleration, λ 0 and L e q are the mass per unit length and relaxation length of the catenary, respectively. D is the stiffness of catenary, α and β are integration constants depending on the position of the two end points and the total catenary mass, u is related to inclination angle ϕ as follows:
tan ϕ = sinh ( u )
The parameter values u1 and u2 represent the positions of the catenary end points p1 and p2 in parameter space. Furthermore, the force f1 and f2 acting on the two end points can be expressed by u1 and u2.
f 1 , x = c , f 1 , y = c sinh ( u 1 )
f 2 , x = c , f 1 , y = c sinh ( u 2 )

2.4. Coupling Scheme of Aero-Hydrodynamics

CFD technology and catenary model are employed to acquire fluid load and mooring load, respectively, and using the DFBI superposition motion approach to predict the body dynamic responses; thus, a coupling scheme is required to simulate dynamic responses of OF-VAWT in wind-wave environment, and the relationship between each module is shown in Figure 3. In the aerodynamic module, creating the rotate and air domain according to rotor geometry and discretizing them based on FVM, particularly, rotate motion is achieved by using sliding mesh. By importing wind field, the aerodynamic solution is completed by solving URANS equation. In addition, the procedure for creating, discretizing and solution in hydrodynamic module is similar to aerodynamic module, except for importing wave and current rather than wind field, and overset mesh is used to adapt free motion. Among them, Eulerian Multi-Phase (EMP) flow and VOF theory are employed to create gas-fluid two-phase flow and capture the interface precisely, respectively, so that air and water domains are combined together. For the dynamic module, six-degree system motion could be obtained by importing the aero-hydro-mooring dynamic load and moment of inertia into DFBI equation, and single-degree rotate motion is further achieved based on superposition motion method. Finally, the dynamic responses calculated, such as displacement, velocity, angle, and position, are returned to aerodynamic, hydrodynamic and mooring modules to conduct the next-step solution.

3. Validation of Numerical Approach

The accuracy of present numerical approach considering aero-hydrodynamic coupling needs to be validated by model experiments. However, there are rare model tests for straight-bladed OF-VAWT in wind-wave coupled environments, thus verifying the aerodynamic and hydrodynamic modules separately is common for scholars [16].

3.1. Aerodynamic Validation

A suitable and available research containing a three-blade vertical axis wind turbine experiment is supposed to verify the present aerodynamic module. Maeda et al. [26] conducted an aerodynamic model experiment of straight-bladed VAWT to investigate the effect of number of blades on aerodynamic performance. (see Figure 4a) A rotor with a radius of 1 m and a height of 1.2 m is selected as the experimental object, and NACA0021 is chosen as its airfoil section with a chord length of 0.265 m, the number of blades is 3. The power for different tip speed ratios (TSR) was measured in an 8 m/s uniform flow. Furthermore, some scholars have raised numerical research for this experiment based on the improved delayed detached eddy simulation (IDDES) model [9], so the experimental data, the IDDES model data and the present model data are compared in the following parts.
Attentionally, Figure 4b illustrates that the incidence angle of airfoil is determined by wind and angular velocity, so the dimensionless parameter TSR is used as a criterion in the study of aerodynamic loads and performance (see Equation (18)). Furthermore, wind turbines are developed to utilize wind energy so that generated power coefficient CP is a significant index to evaluate the quality, as shown in Equation (19).
T S R = ω R U 0
C P = P 0.5 ρ A U 0 3 = Q ω 0.5 ρ D H U 0 3
where ω represents the rotating angular speed of blades, R is the radius, and U0 is the incoming flow speed. ρ is the density of air. Swept wind area A is equal to the product of height H and diameter D, and power P is equal to the torque Q multiplied by the angular velocity ω.
The power coefficient CP with different TSR of the experiment, IDDES model and present work are compared in Figure 5a, and it is obvious that the trend of these three results is close. However, numerical solution is completed in the ideal environment rather than the experimental environment. This may be a major reason why the numerical results of the present study and IDDES model are slightly larger than the experimental values. To further investigate the period variation of the blade, the power coefficient CP of single blade at the best TSR of 1.78 is shown in Figure 5b. Similarly, the results of the present study and IDDES are identical, and both them match well with the experiment.

3.2. Hydrodynamic Validation

To verify the hydrodynamic module, a tank experiment containing a submersible platform is required. The Offshore Code Comparison Collaboration Continuation (OC4) project, led by the National Renewable Energy Laboratory (NREL), is pretty influential in the wind energy field and carried out a tank test for FOWT floating foundation. Among them, free-decay motion is selected to validate the present hydrodynamic module in this section.
DeepCWind is in a great performance and popular in related research of FOWT. Over the years, many hydrodynamic numerical simulations have been carried out by scholars with different methods based on this platform. Coulling et al. [27] and Tran [28] compared the results of model tests and different versions FAST code. The hydrodynamic analysis of this platform was carried out by Luan et al. [29] based on the traditional potential flow theory. With the development of computer science, CFD technology is more often used for related research, such as commercial software STAR-CCM+ [30] and the open-source platform OpenFOAM [16].
Table 2 shows the calculated periods of free-decay motion including present study, previous related work and experiment. It can be seen that the proposed method in this paper matches well with the experiments and related simulations, and the comparison of pitch motion in free-decay is displayed in Figure 6. It is clear that the periodicity and amplitude of the results of several methods are relatively similar, although with slight differences. Overall, the present method of calculating hydrodynamic forces is effective.

4. Model Description

4.1. OF-VAWT System

The floating vertical axis wind turbine is a hot candidate in deep-sea wind energy, a project is to design an OF-VAWT and place it in the South China Sea. Therefore, the impact of the ocean environment on this OF-VAWT needs to be taken into account in advance in this research. Wind and wave parameter settings are determined according to the environmental statistical report, and this selected full-scale floating wind turbine system is composed of a straight-bladed VAWT and a semi-submersible platform as shown in Figure 7a, and the information of properties is listed in Table 3.

4.2. Wind Turbine

The three blades showed good performance in terms of efficiency and noise in a model experiment [26], so this turbine is selected in this paper to combine a floating system. More detailed information about the full-scale rotor, including number of blades, airfoil section, radius and height, is described in Table 4.

4.3. Floating Foundation

In recent years, spar-type floating platforms are often chosen in related research due to their good stability, but high mooring requirements and small platform area are the major drawbacks. However, the semi-submersible platform does not have these troubles and has been proven to have good performance [27]. A semi-submersible platform named DeepCWind belonging to the OC4 project is selected as the foundation of OF-VAWT. As shown in Figure 7b, one main column and three offset columns contribute to the buoyancy, and several braces are employed to connect each module. Detailed information for geometry is shown in Table 5. In addition, to restrict the displacement, three mooring lines are attached to platform, and properties of the corresponding mooring system refer to Cheng et al. [16].

4.4. Numerical Model of OF-VAWT

4.4.1. Computational Domain

The numerical model based on CFD in this paper is established in the STAR-CCM+ platform. As mentioned above (Section 2.3), the approach combing URANS with k − ω turbulence model is employed to govern fluid flow, creating the two-phase flow and capturing its interface by EMP and VOF, DFBI superposition motion is used to solve the dynamic motion response. Attentionally, establishing and discretizing a computational domain is the precondition of solution. Hence, several geometric regions, including background, overset and rotate domain, are created corresponding to numerical tank, system-body free motion and rotational motion (see Figure 8). Furthermore, volume and surface mesh generators are used to fill these domains with meshes based on FVM.
According to the deep-sea environment, in the numerical tank, the boundary conditions on both two sides and bottom are symmetry, and the top and outlet are pressure outlets. Particularly, the inlet is a velocity inlet where wind and wave are produced. To avoid the potential error induced by drastic motion, the whole system-body free motion domain is filled with overset mesh, and the outline surface of this domain is overset. Additionally, rotating motion is achieved by applying the entire rotate domain to a constant angular velocity, and its outline surface is symmetric.
In addition, to produce a fully-developed wave-wind field, the distance between the velocity inlet and the model is greater than twice the length of the floating structure. While, the distance between the outlet and the model is more than four times the length of the floating structure to consider the effect of wake, especially gas wake. In addition, the distance between the sides and the model is similar to the strategy of velocity inlet. The distance between the top and the model is larger than the rotor height, and the distance between the bottom and the model is greater than twice the platform height.

4.4.2. Strategy and Convergence of Mesh

In the process of meshing, trim and prisms layer mesh generators are used to produce the volume mesh; volume and surface control techniques are combined for mesh refinement. Smaller mesh size is helpful to capture tiny details in the fluid domain, however, it means the computational resource requirement will increase rapidly. Thus, a suitable meshing strategy gives considerations to both accuracy and computational cost.
Figure 9 demonstrates the meshing strategy adopted in this paper by showing a longitudinal section. Modeling the wave shape precisely is essential as wave load is the main contributor to the hydrodynamic load, the entire wave propagation region is contained by the free surface domain. Therefore, using volume control technology refines the mesh in both wave length (x) and wave height (z) directions. In addition, it is also necessary to refine the mesh in wake area as the uniform wind is disturbed significantly after flowing through the rotating rotor. Attentionally, to ensure successful mesh interpolation, the mesh size near the interface between the overset and background domains should be closer, as well as the interface between the overset and sliding domains.
The mesh-independent analysis is required to determine the mesh size under a certain computational model and mesh strategy. The OF-VAWT system in this study will be affected by both uniform wind and regular waves, so there is a complex gas-liquid two-phase domain and energy dissipation during wind-wave propagation. Therefore, the mesh-independent analysis is conducted by checking wave shape. The grids are divided into three types named Mesh1, Mesh2 and Mesh3 according to the difference in the basic size, and the total number of grids are 540, 902 and 1415 million, respectively. Then, wave results of three meshing types are compared with the theoretical values, respectively, and the fully-developed wave results are shown in Figure 10. It is clear that the magnitude of Mesh1 has a large difference from the theoretical value, and there is an obvious phase difference. The results of Mesh2 and Mesh3 are closer to the theoretical value, and the results are acceptable, but the number of meshes of Mesh3 is significantly larger than that of Mesh2. Therefore, Mesh2 is selected as the final grid strategy in this paper, and the total number of grids is 902 million.

4.4.3. Simulation Case

To find the TSR with optimal performance, analyzing the aerodynamic performance of the VAWT separately is essential and then integrating it into the OF-VAWT system. In addition, the VAWT in this paper is the full-scale experimental model [26], and it suffers from tower shadow effect, the same as other HAWTs [31], so the aerodynamic performance reduction is studied briefly. By applying present aerodynamic module, the power and torque coefficients with and without the tower at different TSR are acquired and shown in Figure 11. Obviously, the tower shadow effect makes the value decrease with different degrees, and the reduction is more obvious as the TSR increases. More importantly, the optimal performance of this rotor appears when TSR is approximately 1.75.
To ensure that the rotor is performing well at most times, the rotational speed is determined by optimal TSR of 1.75 in all general cases, including under different wind speed. It is well-known that wave height, wave length, wind speed and propagation direction are the main environmental parameters. However, it can be noticed the inertia of this OF-VAWT model is equal in both the X and Y directions (Table. 3), which means similar motion response in both directions. Thus, the propagation direction of wind-wave is not considered and specified in the X-positive direction.
By referring to the statistics of the South China Sea, it is clear that the average wave height is 6 m and the wind speed is 20 m/s. To investigate the load and performance in different wave height, four cases including GC.1, GC.2, GC.3 and GC.4 are carried out in which wave height is from 2 m to 10 m. We also present five cases to study the effect of different wave length including GC.5, GC.3, GC.6, GC.7 and GC.8. Similarly, the effect of wind speed is considered by conducting four cases including GC.9, GC.10, GC.3 and GC.11. Additionally, it is common that extreme environments may appear suddenly when offshore structures are operating, so two extreme cases are added to compare with general cases. The cases are listed in Table 6.

5. Time-Domain Analysis

Aero-hydrodynamic load and performance in general and extreme cases are obtained by applying the numerical model proposed in this paper. Previous related research mainly concentrates on characteristics in frequency-domain by establishing simplified models [20]. Hence, it is necessary to investigate the time-domain characteristics in typical wind-wave environments before statistical analysis, the time-domain data of GC.3, EC.1 and EC.2 will be presented in this section.
Figure 12 shows the operating state of OF-VAWT including distributions of wind velocity and wave elevation. Obviously, the wave-field near the platform is disturbed slightly and generates reflection, radiation and run-up. Furthermore, there is a significant velocity reduction zone after wind flow through the rotor.
The time-domain forces and bending moments of the axial and lateral directions in a general case are given in Figure 13. Axial thrust attracts more concerns in related works of HAWT [32], however, this VAWT generates the obvious lateral load although axial loads generated by the rotating rotor and platform are significantly greater than lateral loads in head wind-wave conditions, and the period of hydrodynamic load is more than twice of aerodynamics in axial direction while it is similar in lateral direction. More importantly, it is noted that the rotor force is far away from the mass central point. This may be the major reason why aerodynamic bending moment is the dominant component of the total bending moment although the aerodynamic force is relatively small compared to hydrodynamic force in both two directions. In a word, the effect of aerodynamic loads generated by rotor on floating systems presented in this part has implications for engineering design.
Considerable efforts have been made to study the effect of floating foundation motion on turbine aerodynamics in offshore FOWTs by many researchers. It indicates that the thrust and torque values are stable and nearly constant when HAWT is fixed [32,33]. This may be the reason why the aerodynamic frequency is close to motion frequency, such as pitch [34] and surge motion [30,35]. For VAWT, the oscillation of aerodynamics is significant. Cheng et al. [36] noted that the aerodynamic value of two-bladed VAWT varies from zero to twice mean value in a period and revealed the two-peak characteristics. However, the aerodynamics of three-bladed VAWT in this paper is not significantly disturbed by motion on period except for slight fluctuations in amplitude. Figure 14 shows the aerodynamic performance curves in one period, it is evident that the main power of a single blade is generated in the upwind area with a single-peak phenomenon. As a result, the three-peak characteristic of total torque curve is produced due to the sequential passed three blades. Thus, improving the performance in downwind area could be a feasible way to promote application of VAWT.
To date the existing works that can consider the aero- hydrodynamics of FOWT mainly adopted the way modeling and integrating each simplified module, respectively, including FloVAWT [37], HAWC2 [38], SIMO-RIFLEX-AC [36], and SIMO-RIFLEX-DMS [22]. However, the characteristics in frequency-domain are the major concentration due to its limitation in theory. To reveal the performance difference in general and extreme environments, torque, pitch and heave motion are compared in time-domain, as shown in Figure 15. Obviously, the aerodynamic torque peak of operating rotor in EC.1 is larger than twice it in GC.3. However, the expectation of torque reduction does not appear by stopping operation of the rotor in EC.2, and the average, as well as peak and valley, are more serious. The results indicate that the aerodynamic load of the OF-VAWT should be reduced by adding control module adjusting the pitch angle under extreme cases. In addition, motions such as pitch and heave are more violent when wind speed and wave height are larger. Further, two pictures (see Figure 16) are given to demonstrate the floating state in the two extreme environments. In the EC.1 case, the sea level is close to up-deck of OF-VAWT due to obvious pitch motion, and there is a long distinct zone of velocity reduction induced by the operating rotor in air domain. For worse conditions (EC.2), the whole deck is almost submerged under violent pitch motion and rough seas. There is still an obvious decrease in wind speed due to resistance induced by stationary rotor, but with a smaller zone.

6. Discussion

The load and performance characteristics in time-domain were analyzed in the previous section, particularly, the aerohydrodynamic load attract attention from structural safety and survival. Some researchers have conducted a series of wind-wave cases investigating the OF-VAWT performance [39,40], but the effect of each environmental variable has not yet been concluded. In this section, we compared the differences in motions and power under different cases to discuss the effect of environmental parameters, including wave height, wave length and wind speed, on OF-VAWT performance.
Wave height as an important part of the wave parameters represents the elevation of the wave peak from the wave valley, the wave shape is steeper as the wave height increases. Figure 17 shows the statistics for variable wave height cases, mean and standard deviation (STD) values are combined to analyze the trend and destabilization of data. Obviously, the pitch motion is always significant compared to heave and roll, and roll motion increased rapidly in extreme environments although it is negligible in general cases. The motion mean value increases with wave height, especially for pitch motion, and the amplitude range of all motion grows with wave height by observing the STD. In addition, there is a slight reduction of the power coefficient when wave height increases, particularly, the power coefficient in extreme cases is unacceptable.
Similarly, wave length represents the distance between two wave peaks, the wave shape is smoother as the wave length increases. By analyzing Figure 18, it is obvious that the motions increase first and then decrease in general cases, however, the STD is of a slight difference. The trend of power coefficient is opposite to motion in variable wave length cases. Attentionally, the mean value of both motion and power have obvious changes induced by wave, but the STD value is stable. However, both of them are susceptible to an extreme environment.
The wind speed indicating wind strength is important to influence aerodynamics. Figure 19 shows the motion and power at different wind speeds. Different from wave impact, it is obvious that both mean and STD values of motion increase rapidly with wind speed, and the amplitude of motion is more susceptible to wind than waves, for example, the max mean value of pitch in general wind speed cases is approximately 7 degrees, while it is 5 degrees in general wave cases. Furthermore, the power coefficient is in decline with increased wind speed, and the maximum reduction of power coefficient induced by wind speed is significant compared to it induced by waves. However the STD is relatively independent of wind speed in general cases.
It is obvious that there is a large difference in motion and in the power of extreme cases and general cases. With wave height or length increases, the maximum mean value of motion in general cases is merely two-thirds of it in extreme cases. However, this value in general wind speed cases is close to extreme cases. A similar pattern appears in the power coefficient variation. This comparison indicates that wind speed is more crucial than wave parameters to the OF-VAWT system. In addition, extreme environments change both the mean and STD value apparently, especially for power coefficient. For instance, the pitch motion in extreme cases is approximately 7.5 degrees while it is approximately 4 in general cases. Particularly, the power coefficient in extreme cases is only 0.075, which means the wind generator is not effective any longer.
Researchers found that VAWTs perform better in wind farms considering wake effect [41]. Further, it is noted that several offshore VAWTs were developed or produced recently, such as 1MW-S2x designed by SeaTwirl and 30kW-S1 operated in Lysekil. However, these products or concepts are mainly developed for tests and validation, and almost no floating VAWTs for deep sea exist. Therefore, the findings in this paper reveal the physics in advance and have implications for designing OF-VAWTs in the future.
It is worth noting that the prediction of aero-hydrodynamics using CFD is accurate and reliable, however, the time cost is a trouble for promotion and application. Therefore, there are still challenges for VAWT development. For instance, the detailed analysis in aerodynamics is important for VAWT improvement, such as dynamics stall effect [42,43,44,45]. Some scholars have also considered small VAWT development [46]. To balance the time cost and accuracy, researchers proposed and improved aerodynamic models for VAWT, such as double-multiple streamtubes [47,48]. These fields are still the key point in VAWT in the future.

7. Conclusions

In this study, a numerical model of OF-VAWT considering aero-hydro-mooring dynamics was developed by combining CFD technology, DFBI superposition motion method and the catenary model. First, the validation of this method was conducted by comparing the aero-hydrodynamic data with two model experiments, respectively. Then, the time-domain load and performance of OF-VAWT in both general and extreme environments were presented and analyzed. Finally, we discussed the performance distributions with different environmental parameters by the statistics and compared them with extreme responses. Several conclusions were obtained as follows.
  • The aerodynamic bending moment is always dominant in a total bending moment, and the aerodynamic torque induced by rotor resistance in extreme environments is significant. These results indicate that OF-VAWT has the requirement for a control module adjusting the pitch angle under extreme cases to reduce aerodynamic loads although VAWT has the ability to harness the energy in each direction without steering.
  • In a wind-wave environment, the trend of power coefficient is opposite to motion. The performance variation is different in each case, quantitative analysis indicates the performance of OF-VAWT is more susceptible to wind compared to waves.
  • In extreme environments, the power coefficient generated by rotors is unacceptable under a violent motion, nearly a quarter of the normal. In addition, there is severe green water induced by violent motion in enormous waves, and even the deck is almost submerged in the worst condition.
In this paper, we concluded the loads and performance of straight-blade OF-VAWT in both general and extreme environments; a comparison shows that OF-VAWT is more susceptible to wind speed than waves. The findings of this study have implications for engineering design. However, the present work concentrates on aero-hydrodynamic performance, and the structural responses are not considered, thus related fluid–structure interaction research will be conducted in the future.

Author Contributions

W.L. (Wenping Luo): simulation, validation, writing—original draft, writing—review and editing, visualization; W.L. (Weiqin Liu): funding acquisition, supervision, methodology; M.Y.: supervision, methodology; S.C.: investigation, data curation, software, information retrieval; X.S.: supervision, methodology; W.W.: supervision, methodology. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No.: 52071243, 52101371) and the National Defense Basic Research Program of China (No.: JCKY2020206B037).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data is in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaldellis, J.K.; Apostolou, D.; Kapsali, M.; Kondili, E. Environmental and social footprint of offshore wind energy. comparison with onshore counterpart. Renew. Energy 2016, 92, 543. [Google Scholar] [CrossRef]
  2. Ismail, M.F.; Vijayaraghavan, K. The effects of aerofoil profile modification on a vertical axis wind turbine performance. Energy 2015, 80, 20–31. [Google Scholar] [CrossRef]
  3. Borg, M.; Collu, M.; Kolios, A. Offshore floating vertical axis wind turbines, dynamics modelling state of the art. part ii: Mooring line and structural dynamics. Renew. Sustain. Energy Rev. 2014, 39, 1226–1234. [Google Scholar] [CrossRef] [Green Version]
  4. Rosenthal, W.; Lehner, S. Rogue waves: Results of the maxwave project. J. Offshore Mech. Arct. Eng. 2008, 130, 21006. [Google Scholar] [CrossRef]
  5. Chou, J.S.; Chiu, C.K.; Huang, I.K.; Chi, K.N. Failure analysis of wind turbine blade under critical wind loads. Eng. Fail. Anal. 2013, 27, 99–118. [Google Scholar] [CrossRef]
  6. Zhang, X.; Wang, M.; Wang, Y.; Liu, Y. Research on Failure of Offshore Wind Turbine Tower in Typhoon. In Proceedings of the 2018 OCEANS-MTS/IEEE Kobe Techno-Oceans (OTO), Kobe, Japan, 28–31 May 2018. [Google Scholar]
  7. Lechuga, A. Were freak waves involved in the sinking of the tanker”prestige”? Nat. Hazards Earth Syst. Sci. 2006, 6, 973–978. [Google Scholar] [CrossRef]
  8. Su, J.; Li, Y.; Chen, Y.; Han, Z.; Zhou, D.; Zhao, Y.; Bao, Y. Aerodynamic performance assessment of φ-type vertical axis wind turbine under pitch motion. Energy 2021, 225, 1. [Google Scholar] [CrossRef]
  9. Lei, H.; Zhou, D.; Bao, Y.; Chen, C.; Ma, N.; Han, Z. Numerical simulations of the unsteady aerodynamics of a floating vertical axis wind turbine in surge motion. Energy 2017, 127, 1. [Google Scholar] [CrossRef]
  10. Lei, H.; Zhou, D.; Lu, J.; Chen, C.; Han, Z.; Bao, Y. The impact of pitch motion of a platform on the aerodynamic performance of a floating vertical axis wind turbine. Energy 2017, 119, 369. [Google Scholar] [CrossRef]
  11. Wu, J.; Meng, L.; Zhao, Y.; He, Y. Coupled aerodynamic and hydrodynamic analysis of floating offshore wind turbine using cfd method. Trans. Nanjing Univ. Aeronaut. Astronaut. 2016, 33, 80–87. [Google Scholar]
  12. Tran, T.T.; Kim, D.H. A cfd study of coupled aerodynamic-hydrodynamic loads on a semisubmersible floating offshore wind turbine. Wind Energy 2018, 21, 70–85. [Google Scholar] [CrossRef]
  13. Wang, H.; Hu, Z.Q.; Meng, X.Y. Dynamic performance investigation of a spar-type floating wind turbine under different sea conditions. China Ocean Eng. 2018, 32, 256–265. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Kim, B. A fully coupled computational fluid dynamics method for analysis of semi-submersible floating offshore wind turbines under wind-wave excitation conditions based on oc5 data. Appl. Sci. 2018, 8, 2314. [Google Scholar] [CrossRef] [Green Version]
  15. Guo, Y.; Liu, L.Q.; Li, Y.; Xiao, C.S.; Tang, Y.G. The surge-heavepitch coupling motions of the ϕ-type vertical axis wind turbine supported by the truss spar floating foundation. J. Hydrodyn. 2019, 31, 669–681. [Google Scholar] [CrossRef]
  16. Cheng, P.; Huang, Y.; Wan, D. A numerical model for fully coupled aero-hydrodynamic analysis of floating offshore wind turbine. Ocean Eng. 2019, 173, 183–196. [Google Scholar] [CrossRef]
  17. Li, X.; Zhu, C.; Fan, Z.; Chen, X.; Tan, J. Effects of the yaw error and the wind-wave misalignment on the dynamic characteristics of the floating offshore wind turbine. Ocean Eng. 2020, 199, 106960. [Google Scholar] [CrossRef]
  18. Gao, J.; Griffith, D.T.; Sakib, M.S.; Boo, S.Y. A semi-coupled aeroservo-hydro numerical model for floating vertical axis wind turbines operating on tlps. Renew. Energy 2022, 181, 692–713. [Google Scholar] [CrossRef]
  19. Lee, H.; Poguluri, S.K.; Bae, Y.H. Development and verification of a dynamic analysis model for floating offshore contra-rotating vertical axis wind turbine. Energy 2022, 240, 122492. [Google Scholar] [CrossRef]
  20. Borg, M.; Shires, A.; Collu, M. Offshore floating vertical axis wind turbines, dynamics modelling state of the art. part i: Aerodynamics. Renew. Sustain. Energy Rev. 2014, 39, 1214–1225. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, L.; Guo, Y.; Zhao, H.; Tang, Y. Motions of a 5 mw floating vawt evaluated by numerical simulations and model tests. Ocean Eng. 2017, 144, 21–34. [Google Scholar] [CrossRef]
  22. Wang, K.; Hansen, M.O.L.; Moan, T. Dynamic analysis of a floating vertical axis wind turbine under emergency shutdown using hydrodynamic brake. Energy Procedia 2014, 53, 56–69. [Google Scholar] [CrossRef]
  23. Vita, L. Offshore Floating Vertical Axis Wind Turbines with Rotating Platform. Ph.D. Thesis, Technical University of Denmark (DTU), Denmark, 2011. [Google Scholar]
  24. Jin, X.; Zhao, G.; Gao, K.; Ju, W. Darrieus vertical axis wind turbine: Basic research methods. Renew. Sustain. Energy Rev. 2015, 42, 212–225. [Google Scholar] [CrossRef]
  25. Bianchini, A.; Balduzzi, F.; Ferrara, G.; Ferrari, L.; Persico, G.; Dossena, V.; Battisti, L. Detailed analysis of the wake structure of a straight blade h-darrieus wind turbine by means of wind tunnel experiments and computational fluid dynamics simulations. J. Eng. Gas Turbines Power 2018, 140, 32604. [Google Scholar] [CrossRef]
  26. Maeda, T.; Kamada, Y.; Murata, J.; Furukawa, K.; Yamamoto, M. Effect of number of blades on aerodynamic forces on a straight-bladed vertical axis wind turbine. Energy 2015, 90, 784. [Google Scholar]
  27. Coulling, A.J.; Goupee, A.J.; Robertson, A.N.; Jonkman, J.M.; Dagher, H.J. Validation of a fast semi-submersible floating wind turbine numerical model with deepcwind test data. J. Renew. Sustain. Energy 2013, 5, 23116. [Google Scholar] [CrossRef]
  28. Tran, T.T.; Kim, D.H. The coupled dynamic response computation for a semi-submersible platform of floating offshore wind turbine. J. Wind Eng. Ind. Aerodyn. 2015, 147, 104–119. [Google Scholar] [CrossRef]
  29. Luan, C.; Gao, Z.; Moan, T. Modelling and Analysis of a Semi-Submersible Wind Turbine with a Central Tower with Emphasis on the Brace System. In Proceedings of the Asme International Conference on Ocean, Nantes, France, 9–14 June 2013. [Google Scholar]
  30. Tran, T.T.; Kim, D.H. A cfd study into the influence of unsteady aerodynamic interference on wind turbine surge motion. Renew. Energy 2016, 90, 204–228. [Google Scholar] [CrossRef]
  31. Mohammadi, E.; Fadaeinedjad, R.; Naji, H.R. Using a new wind turbine emulator to analyze tower shadow and yaw error effects. Energy Convers. Manag. 2018, 174, 378–387. [Google Scholar] [CrossRef]
  32. Huang, Y.; Cheng, P.; Wan, D. Numerical analysis of a floating offshore wind turbine by coupled aero-hydrodynamic simulation. J. Mar. Sci. Appl. 2019, 18, 82. [Google Scholar] [CrossRef]
  33. Borg, M.; Collu, M. A comparison between the dynamics of horizontal and vertical axis offshore floating wind turbines. Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci. 2015, 373, 20140. [Google Scholar] [CrossRef] [Green Version]
  34. Fang, Y.; Duan, L.; Han, Z.; Zhao, Y.; Yang, H. Numerical analysis of aerodynamic performance of a floating offshore wind turbine under pitch motion. Energy 2020, 192, 116. [Google Scholar] [CrossRef]
  35. Farrugia, R.; Sant, T.; Micallef, D. A study on the aerodynamics of a floating wind turbine rotor. Renew. Energy 2016, 86, 770–784. [Google Scholar] [CrossRef]
  36. Cheng, Z.; Madsen, H.A.; Gao, Z.; Moan, T. A fully coupled method for numerical modeling and dynamic analysis of floating vertical axis wind turbines. Renew. Energy 2017, 107, 604. [Google Scholar] [CrossRef]
  37. Collu, M.; Borg, M.; Shires, A.; Brennan, F.P. FloVAWT: Progress on the development of a coupled model of dynamics for floating offshore vertical axis wind turbines. In Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering, American Society of Mechanical Engineers, Nantes, France, 9–14 June 2013; p. V008T09A045. [Google Scholar]
  38. Larsen, T.J.; Madsen, H.A. On the way to reliable aeroelastic load simulation on vawts. Proc. EWEA 2013, 9, 2013. [Google Scholar]
  39. Cheng, Z.; Wang, K.; Gao, Z.; Moan, T. A comparative study on dynamic responses of spar-type floating horizontal and vertical axis wind turbines. Wind Energy 2017, 20, 305–323. [Google Scholar] [CrossRef] [Green Version]
  40. Cheng, Z.; Wang, K.; Gao, Z.; Moan, T. Dynamic response analysis of three floating wind turbine concepts with a two-bladed Darrieus rotor. J. Ocean Wind Energy 2015, 2, 213–222. [Google Scholar] [CrossRef] [Green Version]
  41. Hansen, J.T.; Mahak, M.; Tzanakis, I. Numerical modelling and optimization of vertical axis wind turbine pairs: A scale up approach. Renew. Energy 2021, 171, 1371–1381. [Google Scholar] [CrossRef]
  42. Bangga, G.; Hutomo, G.; Wiranegara, R.; Sasongko, H. Numerical study on a single bladed vertical axis wind turbine under dynamic stall. J. Mech. Sci. Technol. 2017, 31, 261–267. [Google Scholar] [CrossRef]
  43. Bangga, G.; Hutani, S.; Heramarwan, H. The Effects of Airfoil Thickness on Dynamic Stall Characteristics of High-Solidity Vertical Axis Wind Turbines. Adv. Theory Simul. 2021, 4, 2000204. [Google Scholar] [CrossRef]
  44. Ferreira, C.S.; Bijl, H.; Van Bussel, G.; Van Kuik, G. Simulating dynamic stall in a 2D VAWT: Modeling strategy, verification and validation with particle image velocimetry data. J. Phys. Conf. Ser. 2007, 75, 012023. [Google Scholar] [CrossRef] [Green Version]
  45. Buchner, A.J.; Soria, J.; Honnery, D.; Smits, A.J. Dynamic stall in vertical axis wind turbines: Scaling and topological considerations. J. Fluid Mech. 2018, 841, 746–766. [Google Scholar] [CrossRef]
  46. Bianchini, A.; Bangga, G.; Baring-Gould, I.; Croce, A.; Cruz, J.I.; Damiani, R.; Erfort, G.; Simao Ferreira, C.; Infield, D.; Nayeri, C.N.; et al. Current status and grand challenges for small wind turbine technology. Wind Energy Sci. 2022, 7, 2003–2037. [Google Scholar] [CrossRef]
  47. Bangga, G.; Dessoky, A.; Lutz, T.; Krämer, E. Improved double-multiple-streamtube approach for H-Darrieus vertical axis wind turbine computations. Energy 2019, 182, 673–688. [Google Scholar] [CrossRef]
  48. Paraschivoiu, I. Double-multiple streamtube model for studying vertical-axis wind turbines. J. Propuls. Power 1988, 4, 370–377. [Google Scholar] [CrossRef]
Figure 1. Definition of local and global coordinates in superposition motion.
Figure 1. Definition of local and global coordinates in superposition motion.
Jmse 11 00185 g001
Figure 2. Force analysis of a mooring line in quasi-static catenary model.
Figure 2. Force analysis of a mooring line in quasi-static catenary model.
Jmse 11 00185 g002
Figure 3. Coupling scheme and modules in this fully coupled model.
Figure 3. Coupling scheme and modules in this fully coupled model.
Jmse 11 00185 g003
Figure 4. Aerodynamic experiment of a VAWT: (a) Testing arrangement (reproduced from [26], with permission from Elsevier, 2022) (b) Parameters in rotation.
Figure 4. Aerodynamic experiment of a VAWT: (a) Testing arrangement (reproduced from [26], with permission from Elsevier, 2022) (b) Parameters in rotation.
Jmse 11 00185 g004
Figure 5. Comparison aerodynamic data of simulation against experiment: (a) Power coefficient against TSR (b) Power coefficient against azimuth angle [26,9]
Figure 5. Comparison aerodynamic data of simulation against experiment: (a) Power coefficient against TSR (b) Power coefficient against azimuth angle [26,9]
Jmse 11 00185 g005
Figure 6. Comparison hydrodynamic data of simulation against experiment: (a) Testing arrangement (reproduced from [27], with permission from AIP Publishing, 2022) (b) Time-domain pitch curve of free−decay motion.
Figure 6. Comparison hydrodynamic data of simulation against experiment: (a) Testing arrangement (reproduced from [27], with permission from AIP Publishing, 2022) (b) Time-domain pitch curve of free−decay motion.
Jmse 11 00185 g006
Figure 7. Geometric model of this OF-VAWT: (a) OF-VAWT system (b) OC4 foundation.
Figure 7. Geometric model of this OF-VAWT: (a) OF-VAWT system (b) OC4 foundation.
Jmse 11 00185 g007
Figure 8. Computational domain and boundary.
Figure 8. Computational domain and boundary.
Jmse 11 00185 g008
Figure 9. Grid partitioning strategy in longitudinal section.
Figure 9. Grid partitioning strategy in longitudinal section.
Jmse 11 00185 g009
Figure 10. Time−domain wave−making amplitude curve of different mesh type.
Figure 10. Time−domain wave−making amplitude curve of different mesh type.
Jmse 11 00185 g010
Figure 11. Power and torque coefficient of rotor with tower CP1, CQ1 and without tower CP2, CQ2.
Figure 11. Power and torque coefficient of rotor with tower CP1, CQ1 and without tower CP2, CQ2.
Jmse 11 00185 g011
Figure 12. The operating OF−VAWT in the gas-liquid flow field of case 3.
Figure 12. The operating OF−VAWT in the gas-liquid flow field of case 3.
Jmse 11 00185 g012
Figure 13. Time−domain curve of aero−hydrodynamic force and bending moment in case 3: (a) Axial force and bending moment (X−axes); (b) Lateral force and bending moment (Y−axes).
Figure 13. Time−domain curve of aero−hydrodynamic force and bending moment in case 3: (a) Axial force and bending moment (X−axes); (b) Lateral force and bending moment (Y−axes).
Jmse 11 00185 g013
Figure 14. Time−domain curve of aerodynamic performance with azimuth angle in case 3.
Figure 14. Time−domain curve of aerodynamic performance with azimuth angle in case 3.
Jmse 11 00185 g014
Figure 15. Time−domain curve of aerodynamic performance and hydrodynamic motion in case 3.
Figure 15. Time−domain curve of aerodynamic performance and hydrodynamic motion in case 3.
Jmse 11 00185 g015
Figure 16. Floating state of OF-VAWT in extreme cases: (a) EC.1; (b) EC.2.
Figure 16. Floating state of OF-VAWT in extreme cases: (a) EC.1; (b) EC.2.
Jmse 11 00185 g016
Figure 17. Aero-hydrodynamic performance in variable wave height case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Figure 17. Aero-hydrodynamic performance in variable wave height case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Jmse 11 00185 g017
Figure 18. Aero-hydrodynamic performance in variable wave length case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Figure 18. Aero-hydrodynamic performance in variable wave length case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Jmse 11 00185 g018
Figure 19. Aero-hydrodynamic performance in variable wind speed case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Figure 19. Aero-hydrodynamic performance in variable wind speed case: (a) Mean and STD value of motion response; (b) Mean and STD value of power coefficient.
Jmse 11 00185 g019
Table 1. Related efforts of FOWT considering aero-hydrodynamic coupling.
Table 1. Related efforts of FOWT considering aero-hydrodynamic coupling.
Author (Year)RotorPlatformAerodynamic ModuleHydrodynamic Module
Wu et al. (2016) [11]NERL-5MWSpar-typeCFDCFD
Tran and Kim (2018) [12]NERL-5MWDeepCWindCFDCFD
Wang et al. (2018) [13]NERL-5MWSpar-typeBEMPotential flow
Zhang and Kim (2018) [14]NERL-5MWDeepCWindCFDCFD
Guo et al. (2019) [15]Φ-typeSpar-typeStreamtubePotential flow
Cheng et al. (2019) [16]NERL-5MWDeepCWindCFDCFD
Li et al. (2020) [17]NERL-5MWDeepCWindBEMPotential flow
Gao et al. (2022) [18]Darrieus-typeTLP-typeVortexPotential flow
Lee et al. (2022) [19]Darrieus-typeTLP-typeStreamtubePotential flow
Table 2. Period comparison of free decaying motion.
Table 2. Period comparison of free decaying motion.
DOFExp [27]Simo/Riflex [29]Fast V8.1 [28]Unsteady CFD [30]naoe-FOAM-SJTU [16]Present
Pitch26.825.825.025.225.825.6
Roll26.926.025.025.3-25.7
Heave17.517.117.717.817.5817.7
Table 3. Properties of this system.
Table 3. Properties of this system.
PropertyParameter
Mass1.37 × 107 kg
CM location below SWL12.53 m
Total structure roll inertia bout CM8.6 ×109 kg × m2
Total structure pitch inertia bout CM8.6 ×109 kg × m2
Total structure yaw inertia bout CM1.226 × 1010 kg × m2
Table 4. Wind turbine geometry.
Table 4. Wind turbine geometry.
PropertyParameter
Aerofoil sectionNACA0021
Number of blades(n)3
Length of chord(c)13.25 m
Rotor diameter(D)100 m
Height of straight-bladed(H)60 m
Pitching angle(β)8 deg
Height of center above main column48 m
Table 5. Floating foundation geometry.
Table 5. Floating foundation geometry.
PropertyParameter
Depth of platform base below SWL (total draft)20 m
Elevation to platform top (tower base) above SWL10 m
Distance between offset columns50 m
Height of upper columns26 m
Height of base columns6 m
Diameter of main column6.5 m
Diameter of offset (upper) column12 m
Diameter of offset (bottom) column24 m
Diameter of braces1.6 m
Table 6. General and extreme cases of OF-VAWT in wind-wave coupling environment.
Table 6. General and extreme cases of OF-VAWT in wind-wave coupling environment.
CaseH (m)λ/LV (m/s)
GC.12120
GC.24120
GC.36120
GC.410120
GC.560.7520
GC.661.2520
GC.761.520
GC.86220
GC.96110
GC.106115
GC.116130
EC.17.5-40
EC.2 (parked)14-60
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, W.; Liu, W.; Yang, M.; Chen, S.; Song, X.; Wu, W. Load Characteristics and Extreme Response of Straight-Bladed Floating VAWT Using a Fully Coupled Model. J. Mar. Sci. Eng. 2023, 11, 185. https://doi.org/10.3390/jmse11010185

AMA Style

Luo W, Liu W, Yang M, Chen S, Song X, Wu W. Load Characteristics and Extreme Response of Straight-Bladed Floating VAWT Using a Fully Coupled Model. Journal of Marine Science and Engineering. 2023; 11(1):185. https://doi.org/10.3390/jmse11010185

Chicago/Turabian Style

Luo, Wenping, Weiqin Liu, Meng Yang, Shuo Chen, Xuemin Song, and Weiguo Wu. 2023. "Load Characteristics and Extreme Response of Straight-Bladed Floating VAWT Using a Fully Coupled Model" Journal of Marine Science and Engineering 11, no. 1: 185. https://doi.org/10.3390/jmse11010185

APA Style

Luo, W., Liu, W., Yang, M., Chen, S., Song, X., & Wu, W. (2023). Load Characteristics and Extreme Response of Straight-Bladed Floating VAWT Using a Fully Coupled Model. Journal of Marine Science and Engineering, 11(1), 185. https://doi.org/10.3390/jmse11010185

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop