Next Article in Journal
Tribological and Thermal Transport of Ag-Vegetable Nanofluids Prepared by Laser Ablation
Next Article in Special Issue
The Influence of Tilt Angle on the Aerodynamic Performance of a Wind Turbine
Previous Article in Journal
Front-Face Fluorescence of Honey of Different Botanic Origin: A Case Study from Tuscany (Italy)
Previous Article in Special Issue
Evaluation of Tip Loss Corrections to AD/NS Simulations of Wind Turbine Aerodynamic Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Characteristics of an Orthopter-Type Vertical Axis Wind Turbine in Shear Flows

by
Rudi Purwo Wijayanto
1,2,*,
Takaaki Kono
3,* and
Takahiro Kiwata
3
1
Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan
2
The Agency of the Assessment and Application of Technology (BPPT), Jakarta 10340, Indonesia
3
Institute of Science and Engineering, Kanazawa University, Kanazawa 920-1192, Japan
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(5), 1778; https://doi.org/10.3390/app10051778
Submission received: 9 January 2020 / Revised: 25 February 2020 / Accepted: 2 March 2020 / Published: 4 March 2020
(This article belongs to the Special Issue Wind Turbine Aerodynamics II)

Abstract

:
To properly conduct a micro-siting of an orthopter-type vertical axis wind turbine (O-VAWT) in the built environment, this study investigated the effects of horizontal shear flow on the power performance characteristics of an O-VAWT by performing wind tunnel experiments and computational fluid dynamics (CFD) simulations. A uniform flow and two types of shear flow (advancing side faster shear flow (ASF-SF) and retreating side faster shear flow (RSF-SF)) were employed as the approaching flow to the O-VAWT. The ASF-SF had a higher velocity on the advancing side of the rotor. The RSF-SF had a higher velocity on the retreating side of the rotor. For each type of shear flow, three shear strengths (Γ = 0.28, 0.40 and 0.51) were set. In the ASF-SF cases, the power coefficients (Cp) were significantly higher than the uniform flow case at all tip speed ratios (λ) and increased with Γ. In the RSF-SF cases, CP increased with Γ. However, when Γ = 0.28, the CP was lower than the uniform flow case at all λ. When Γ = 0.51, the CP was higher than the uniform flow case except at low λ; however, it was lower than the ASF-SF case with Γ = 0.28. The causes of the features of CP were discussed through the analysis of the variation of blade torque coefficient, its rotor-revolution component and its blade-rotation component with azimuthal angle by using the CFD results for flow fields (i.e., horizontal velocity vectors, pressure and vorticity). These results indicate that a location where ASF-SFs with high Γ values dominantly occur is ideal for installing the O-VAWT.

1. Introduction

Since the 2000s, interest in installing small wind turbines (SWTs) in the built environment has been growing [1,2,3,4,5,6,7,8,9]. Wind conditions in the built environment are complex in nature and are characterized by lower wind speeds and higher turbulence because of the presence of obstructions [8,9]. For SWTs to be able to make up their costs within their lifetimes, they should have high efficiency and be placed at sites with high wind speeds, such as coastal sites or high-elevation inland sites. However, in the built environment, keeping the rotational speed of an SWT’s rotor as low as possible is preferable from the viewpoint of aerodynamic noise [10,11]. Therefore, the optimal tip-speed ratio of an SWT in the built environment should be as low as possible, while the maximum power coefficient of the SWT should be as high as possible.
Wind turbines are classified into horizontal-axis wind turbines (HAWTs) and vertical-axis wind turbines (VAWTs), based on the orientation of their rotation axes. Generally, in the built environment, VAWTs are preferable to HAWTs because VAWTs do not suffer, as much as HAWTs, from reduced energy outputs from frequent wind direction changes [12]. Wind turbines are further classified into lift-type wind turbines and drag-type wind turbines, based on the aerodynamic force component that acts on a blade and dominantly contributes to the rotor rotation. With regard to lift-type VAWTs, a lot of research on Darrieus-type VAWTs, including straight-bladed and helical-bladed ones, has been conducted [13,14,15,16]. With regard to drag-type VAWTs, a lot of research on Savonius-type VAWTs has been conducted [17,18,19,20,21]. In general, the optimal tip-speed ratio of a drag-type VAWT is less than 1.0 [11], which is much smaller than that of a lift-type VAWT. Moreover, although the maximum power coefficient of a drag-type VAWT is generally much smaller than that of a lift-type VAWT, the power coefficient of a drag-type VAWT is generally greater than that of a lift-type VAWT at a low tip-speed ratio, of less than 1.0. Therefore, a drag-type VAWT is favorable in the built environment and was researched by our research group.
Our group [22,23] researched a drag-type VAWT called the orthopter-type VAWT (O-VAWT). The O-VAWT is a variable-pitch VAWT; each of the flat-plate blades not only revolves around the main shaft but also rotates around its own blade axis, which is rotationally supported by a pair of connecting arms. We investigated the effects of the number and aspect ratio of the flat-plate blades on the power performance of the O-VAWT in a uniform flow by conducting wind tunnel experiments with an open test section and three-dimensional computational fluid dynamics (CFD) simulations. When the number of the blades was three and the aspect ratio of the blades was 1:1, the maximum power coefficient was 0.25 and the optimal tip-speed ratio was 0.4 [22]. Here, the optimal tip-speed ratios of Savonius-type VAWTs are in the range of 0.45 to 1.0 [21]. Except for several studies that were conducted using a wind tunnel with a very high blocking ratio of the closed test section, the maximum power coefficient of the Savonius-type VAWT was, at most, 0.25 [21]. That is, the O-VAWT has a lower optimal tip-speed ratio than Savonius-type VAWTs, although the maximum power coefficient is relatively high. Therefore, from the viewpoint of aerodynamic noise, the O-VAWT can be more favorable in the built environment as compared to Savonius-type VAWTs. Except for our studies, studies on the power performance of O-VAWTs are very limited. Shimizu et al. [24] investigated the effects of the aspect ratio of the blade on the power performance of an O-VAWT with two blades whose cross-sectional shape was an ellipse by conducting wind tunnel experiments. Bayeul-Line et al. [25] examined the effects of the blade’s cross-sectional shape (elliptical and straight) and the initial blade stagger angle on the performance of an O-VAWT by conducting 2-dimensional CFD simulations. Cooper and Kennedy [26] examined the power performance of an O-VAWT with three blades whose cross-sectional shape was the upstream half of a NACA0010-65 section reflected about the mid-chord by conducting theoretical analysis with a multiple-stream tube model and field measurements. Our group [23] conducted wind tunnel experiments to compare the performance of an O-VAWT with elliptic blades and one with flat-plate blades. By considering the mechanical loss torque, we obtained the maximum power coefficient of 0.246 at a tip-speed ratio of 0.4 for the O-VAWT with elliptic blades and of 0.288 at a tip-speed ratio of 0.4 for the O-VAWT with flat-plate blades. It should be noted that these studies on O-VAWTs were conducted in conditions where the approaching flows had uniform distribution.
To properly conduct a micro-sitting of an O-VAWT in the built environment, it is important to understand the effects of the strong shear approach flow with on the performance of the O-VAWT. Figure 1 illustrates the approaching wind flow to a building. As the wind flow approaches the building, the wind speed decreases, and the pressure increases. Then, the wind flow proceeds along the upwind face of the building and separates at the corners on the roof and the side walls. As the separated wind flow is not obstructed by the building, the pressure decreases, and the wind speed increases. Near the upwind corners, the wind speed increases more than that of the approaching wind. In addition, reverse flow regions are formed between the separated shear layer and the building’s walls. As a result, strong shear flows are formed vertically over the roof surface and horizontally over the side walls. Due to the mixing of momentum, the shear becomes weaker as the flow proceeds downstream. To utilize the increased wind speed over the roof of a building, the effects of building shapes and wind directions on the wind conditions have been investigated (e.g., [27]). Furthermore, the effects of wind conditions, such as wind speed, turbulent intensity, and skew angles, on the potential energy yield and the power performance of a wind turbine have been studied (e.g., [28,29]). In this study, we investigated the effects of horizontal shear flow on the performance of the O-VAWT by conducting wind tunnel experiments and three-dimensional CFD simulations. A uniform flow and two types of shear flow were employed as the approaching flow to the O-VAWT. One type had a higher velocity on the advancing side of the rotor. The other type had a higher velocity on the retreating side of the rotor. For each type of shear flow, we set three different shear strengths.

2. Experimental Approach

2.1. Wind Turbine Model

In this paper, we employed a right-handed Cartesian coordinate system (x1, x2, x3) = (x, y, z), in which the z-direction was aligned in the vertical direction. The wind turbine used in this study was an O-VAWT with three flat-plate blades as shown in Figure 2. The blade had a height of h = 4.00 × 10−1 m, chord length of c = 4.00 × 10−1 m and thickness of 4.0 × 10−3 m. Each of the blades not only revolved around the main shaft but also rotated around its own blade axis, which was rotationally supported by a pair of connecting arms. The distance between the main shaft and one of the blade axes was R = 2.55 × 10−1 m. In addition, each of the blade axes was connected with the main shaft by a chain via sprockets. Since the ratio of the number of teeth on the sprocket of the main shaft to that of the blade axis was 1:2, each of the blades rotated around the own blade axis a half time while the rotor revolved around the main shaft one time. When seen from the top, the rotor revolved around the main shaft counterclockwise and each of the blades rotated around the own blade axis clockwise as shown in Figure 2c. Therefore, by using the angular velocity of the rotor revolution (ω), the angular velocity of the blade rotation is expressed as −ω/2. Furthermore, as shown in Figure 2c, according to the azimuthal angle of a blade (φ), we call the range of 90° < φ < 270° the “upwind region” of the rotor, 0° < φ < 90° and 270° < φ < 360° the “downwind region” of the rotor, 180° < φ < 360° the “advancing side of the rotor” and 0° < φ < 180° the “retreating side” of the rotor. The O-VAWT is designed so that the drag force on a blade is large on the advancing side of the rotor while being small on the retreating side of the rotor.
Figure 2d shows a projected swept area of the rotor (A), which is defined as:
A = (2R + 0.5 c) h.
We define the diameter of the rotor as:
D = 2R.
The O-VAWT had a rotor’s diameter of D = 5.1 × 10−1 m and a projected rotor′s swept area of A = 2.84 × 10−1 m2 can be considered as a micro wind turbine. A small scale wind turbine that has a diameter up to 1.25 m and the swept area up to 1.2 m2 is categorized as a micro wind turbine [4].

2.2. Experimental Setup for Uniform Flow Case

Figure 3a shows the experimental setup for the uniform flow case. The experiments were conducted using a closed circuit wind tunnel with an open test section. The size of the cross section of the wind tunnel outlet was 1.25 m × 1.25 m. The blockage ratio which is defined as the ratio of the projected rotor′s swept area to the wind tunnel outlet area was approximately 18%. The O-VAWT was set in the test section so that the rotor center was at the center of the cross-section of the wind tunnel outlet and 0.850 m downwind of the wind tunnel outlet. Here, we defined the rotor center as the point on the rotational axis of the rotor and at the mid-height of the blades. In addition, we set the origin of the coordinate system at the rotor center, as shown in Figure 2c,d. The rotor was driven by a motor (Mitsubishi Electric, GM-S) and its rotational speed (ω) was monitored by using a digital tachometer (Ono Sokki, HT-5500) and controlled by using an inverter (Hitachi, SJ200). The rotor torque was measured by using a torque meter (TEAC, TQ-AR), which was connected to the motor and shaft via couplings. The output signal of the torque meter was converted by a 16-bit analog-to-digital converter with a sampling interval of 0.5°, and 36,000 items (50 revolutions) of data were stored. To measure the reference wind speed U, an ultrasonic anemometer (Kaijo Sonic, DA-650-3TH and TR-90 AH) was set approximately 2 m upwind of the wind tunnel outlet. The value of U was kept at 8 m/s. The value of tip speed ratio λ, which is defined as:
λ = Rω/U,
was varied from 0.1 to 0.8 with an increment of 0.1.

2.3. Experimental Setup for Shear Flow Cases

To generate a horizontal shear flow, a perforated panel and a splitter plate were installed at the outlet of the wind tunnel as shown in Figure 3b,c and Figure 4. The perforated panel had a width of 1.0 m, a height of 1.5 m and a thickness of 2 × 10−3 m, and was set so that it covered half of the wind tunnel outlet in the horizontal wind direction. Due to the existence of the perforated panel, the pressure upwind of the panel increased and the wind flow rate through the wind-tunnel-outlet area covered by the panel decreased while that through the uncovered wind-tunnel-outlet area increased. When the perforated panel covered the wind-tunnel-outlet area upwind of the retreating side of the rotor, the wind speed of the generated shear flow was higher on the advancing side of the rotor. Hereafter, this type of shear flow is referred to as “advancing side faster shear flow” (ASF-SF). On the other hand, when the perforated panel covered the wind-tunnel-outlet area upwind of the advancing side of the rotor, the wind speed of the generated shear flow was higher on the retreating side of the rotor. Hereafter, this type of shear flow is referred to as “retreating side faster shear flow” (RSF-SF). The splitter plate was set vertically, parallel to the wind tunnel wall and at the center of the wind tunnel outlet to avoid the horizontal component of the wind velocity in the generated shear flow becoming significant. The splitter plate had a width of 0.88 m, a height of 1.25 m and a thickness of 4 × 10−3 m.
To investigate the effects of the strength of the shear flow on the performance of the O-VAWT, we generated three kinds of shear flows by using three perforated panels shown in Table 1. With regard to a staggered round-hole perforated panel, the shielding ratio Φ, which is the ratio of the area that shields the airflow to the whole area of the perforated panel, can be computed by:
Φ   =   1 π   d 2 2 3   L 2 ,
where d is a diameter of a hole and L is the distance between the centers of adjacent holes.
Except for the installation of the perforated panels and the splitter plate, the experimental setup for the measurement of the performance of the O-VAWT was the same as the uniform flow case. Prior to the measurement of the performance of the O-VAWT, we measured the horizontal profiles of the generated shear flows at 0.10 m downwind of and at the center height of the wind tunnel outlet by using an x-type hot-wire probe (Kanomax, 0252R-T5).

2.4. Torque and Power Coefficients

Due to the difficulty of evaluating the mechanical losses of the bearings, the sprockets and the chains, this study considers only the aerodynamic torque generated by the blades as the rotor torque of the O-VAWT. The aerodynamic torque generated by the blades was computed by:
TB = TwBTwoB,
where TwB was the measured aerodynamic torque generated by the rotor when the blades were not removed; and TwoB was the measured aerodynamic torque generated by the rotor when the blades were removed. It is worth noting that, generally, when the blades were not removed from the rotor, the rotor generated positive torque while the motor acted as a load to keep the value of ω constant. Conversely, at all tip speed ratios, when the blades were removed from the rotor, the rotor generated negative torque while the motor acted as the driving force of the rotor revolution. Therefore, at all tip speed ratios, TB was higher than TwB.
The power coefficient describes that fraction of the power in the wind that may be converted by the turbine into mechanical work [30] and is defined in this study as:
C P = T B ω 0.5 ρ A U 0 3 ,  
and the torque coefficient is defined as:
C T = T B   0.5 ρ A U 0 2 R .
Here, the U0 is the time-mean stream-wise velocity, u ¯ ( x , y , z ) , averaged over the projected rotor’s swept area at x = 0 and is computed by:
U 0 = ( R + 0.5 c ) R u ¯ ( 0 , y , 0 ) d y   2 R + 0.5 c .

3. Numerical Approach

The CFD software utilized to simulate the wind flow field was ANSYS Fluent 17.2 [31,32]. The numerical approach was based on our previous paper [22,33], in which the CFD simulations with the delayed detached eddy simulation (DDES) turbulence model of flow around the O-VAWT were conducted and the validities of the grid resolution and the time-step size were confirmed.

3.1. Governing Equations and Discretization Method

The flow field around the wind turbine was assumed to be incompressible and isothermal. The DDES turbulence mode treats near-wall region in a manner like a Reynolds-averaged Navier–Stokes (RANS) turbulence model and treats the rest of the flow field in a manner like a large-eddy simulation (LES) turbulence model [34]. This model has the potential to achieve higher accuracy than RANS models and save a large number of computing resources compared with pure LES models. The governing equations for the CFD simulation with the DDES turbulence model based on the Spalart–Allmaras (SA) model are the continuity equation:
u i x i = 0 ,
the Navier–Stokes equation:
u i t + u i u j x j = 1 ρ p x i + ν x j ( u i x j + u j x i 2 3 δ i j u i x j ) x j [ ν ( u i x j + u j x i ) 2 3 k δ i j ] ,
and the transport equation for the kinematic eddy viscosity ν :
ν t + ν u i x i = C b 1 ν ( S + ν κ 2 d D D E S 2 ( 1 χ 1 + χ f ν 1 ) ) + 1 σ ν [ x j { ( ν + ν ) ν x j } + C b 2 ( ν x j ) 2 ] C w 1 f w ( ν d D D E S ) 2
where ui is the wind-velocity component in the xi direction; p is the pressure; ν is the kinematic viscosity; t is the time; ρ is the air density; k is the turbulence kinetic energy; δij is the Kronecker delta; dDDES is the DDES length scale; χ is ( ν / ν ); S is a scalar measure of the deformation tensor; f ν 1 and f w are damping functions; and C b 1 , C b 2 , C w 1 , σ ν and κ are constants.
The DDES length scale is computed by:
d D D E S = d f d m a x   ( 0 , d c d e s Δ m a x )
where d is the distance to the closest wall; cdes is the empirical constant; and max is the maximum edge length of the local computational cell, i.e., max = max (∆x, ∆y, ∆z). The switching between the RANS and the LES mode depends on the following shielding function;
f d = 1 t a n h ( ( 8 r d ) 3 )
and
r d = ν u i , j u i , j κ 2 d 2
where ui,j is the velocity gradient. The damping functions and closure coefficients are as follows:
f ν 1 = χ 3 c ν 1 3 + χ 3 ,   f w = g [ 1 + c w 3 6 g 6 + c w 3 6 ] 1 / 6 ,   g = r + c w 2 ( r 6 r ) ,   c w 1 = c b 1 κ 2 + ( 1 + c b 2 ) σ ν ,
c b 1 = 0.1355 ,   c b 2 = 0.622 ,   c ν 1 = 7.1 ,   c w 2 = 0.3 ,   c w 3 = 2.0 ,   σ ν = 2 3 ,   c d e s = 0.65 ,   κ = 0.4187 .
The governing equations are discretized by the finite-volume method. The advection terms of the Navier–Stokes equations are discretized by the bounded central-difference scheme. The advection term of the transportation equations for ν is discretized by a second-order upwind scheme. Other spatial derivatives are discretized by the central-difference scheme. The time integration is performed using the second-order implicit method.

3.2. Numerical Setup

The computational domain, the computational meshes and the boundary conditions are shown in Figure 5. As the same as the experimental setup, the origin of the coordinate system is defined at the center of the O-VAWT. The modeled O-VAWT was comprised of three blades, one main shaft and two sets of connecting arms. To reduce the computational cost, other components, such as the chains and sprockets are omitted. The sizes of these components are the same as those used in the experiment. The computational domain consists of three blade domains, one rotor domain and one far-field domain. The blade domain included one of the blades and rotates around each blade axis. The rotor domain included these three blade domains, the main shaft and the connecting arms and rotates around the main shaft. The far-field domain was a stationary domain and its size was 23.5D × 17D × 5D. Except for uniform flow cases, a splitter plate with the same thickness as the experiment was set at y = 0 and x = −11.37D to −1.67D. In all domains, only unstructured meshes were used. Based on our previous mesh-resolution dependency tests [33], the number of computational cells in each domain were set as shown in Table 2. The total number of computational cells was approximately 10 million. All surfaces of the solid components were covered with boundary-layer meshes. The first grid nodes over the surface of the blades were y+ < 1 in all run cases.
At the inlet boundary, the distributions of the stream-wise wind velocity shown in Table 3 were implemented. These distributions were set so that when the O-VAWT is absent, the distributions of the time-mean values of u at x ≈ −0.147 D, which corresponds to 0.1 m downwind of the wind tunnel outlet, matched well with those of the wind tunnel experiment, which were generated by using the three kinds of perforated panels with different shielding ratios Φ, as shown in Figure 6a,b. Here, UH and UL are the maximum and minimum streamwise velocities in the profile of a shear flow measured in the wind tunnel experiment; Γ is the velocity ratio defined as:
Γ = U H U L U H + U L
It is worth noting that amplification factors of Γ = 0.28 and Γ = 0.51 are considered as 1.35 and 1.65, respectively, by computing UH/U. The amplification factor is defined as the ratio of wind speed in the case where there are buildings to wind speed in the case where these buildings are removed. The separated shear flow from a building with an aspect ratio of 1:1:2 reaches an amplification factor of 1.2 [35]. In addition, the separated shear flow from a building with an aspect ratio of 1:1:6 reaches an amplification factor of 1.7 [36]. At x = 0, the values of Γ do no change, as shown in Figure 6c,d. However, due to the momentum diffusion, the horizontal gradients of the time-mean streamwise velocity become weaker. One of the reasons for the relatively large discrepancies between the profiles obtained by the experiment and the CFDs can be that the turbulence intensities of the CFDs are significantly small compared to those of the experiments, as shown in Figure 7. Even though very high values of turbulence intensity were set at the inlet boundary, the turbulence intensities dissipated rapidly and became very small at x = 0 as compared to those of the experiments. As the setting of high turbulence intensity at the inlet boundary leads to computational instability, we set no perturbation condition at the inlet boundary. Table 4 shows the values of U0 computed by Equation (8) and the Reynolds number which is defined as:
R e = U 0 D ν
At the outlet boundary, the pressure outlet condition with p = 0 was imposed. On the surface of the O-VAWT and the splitter plate, the no-slip boundary conditions were set. The sliding mesh technique was used to couple the rotational domains and the stationary domain. The direction of the rotor and the blade rotations are counterclockwise and clockwise, respectively, when viewed from the top (Figure 5a). By changing the rotational speed of the rotor ω, the tip speed ratio λ was set at 0.2, 0.4, 0.5, 0.6, or 0.8. The time step sizes were set as dt = 0.5°/ω.

3.3. Torque and Power Coefficients

As mentioned in sub-Section 2.4, this study considers only the aerodynamic torque generated by the blades as the rotor torque of the O-VAWT. Since each of the blade axes was connected with the main shaft by a chain via sprockets, the aerodynamic torque on each of the blades about each of the blade axes was transmitted through the chain and contributed to the torque about the main shaft. Therefore, the rotor torque generated by a blade at an azimuthal angle φ is expressed as:
T B ( φ ) = T B _ r e v ( φ ) + T B _ r o t ( φ ) ,
where TB_rev(φ) is the conventional blade torque that is calculated by multiplying the rotor radius and the component of the aerodynamic force on the blade at φ in the rotor-revolution direction; and TB_rot(φ) is the torque generated by the component of the aerodynamic force on the blade at φ in the blade-rotation direction about the blade axis. Hereafter, we call TB (φ) the “blade torque,” TB_rev(φ) the “rotor-revolution torque” and TB_rot(φ) the “blade-rotation torque.” The blade torque coefficient (CTB), the rotor-revolution component (CTB_rev) and the blade-rotation component (CTB_rot) are defined as:
C T B ( φ ) = T B   ( φ ) 0.5 ρ A U 0 2 R ,
C T B _ r e v ( φ ) = T B _ r e v   ( φ ) 0.5 ρ A U 0 2 R ,
and
C T B _ r o t ( φ ) = T B _ r o t   ( φ ) 0.5 ρ A U 0 2 R .
It should be noted that CTB, CTB_rev and CTB_rot are coefficients of one blade.
The CFD simulations were conducted for eight revolutions of the rotor. Using the data of the last two rotor revolutions, the torque coefficient CT was computed by the following formula:
C T = n N e N s N = N s + 1 N e 0 2 π C T B   d φ .
Here, n (= 3) is the number of the blades; Ns (= 6) is the number of the rotor revolutions before starting the computation of CT; Ne (= 8) is the number of the rotor revolutions before finishing the computation of CT. The power coefficient of CP was computed by Equation (6).

4. Results and Discussion

In this section, the results of the wind tunnel experiments and the CFD simulations for the power performance of the O-VAWT, such as the dependency of the power and torque coefficients on the tip speed ratio, the variations of the torque coefficients with azimuthal angle, are presented for the uniform flow case and the shear flow cases. Subsequently, the causes of the features of the power performance of the O-VAWT are discussed based on the CFD results of the flow fields.

4.1. Performance in Uniform Flow

Figure 8a shows the power and torque coefficients (CP and CT) of the O-VAWT in the uniform flow. The CFD results are in good agreement with the experimental ones. The optimal tip speed ratio at which CP becomes the maximum is less than unity; 0.4 in the experiments and 0.5 in the CFD simulations. As mentioned in the introduction, this low optimal tip speed ratio is a favorable feature for the built environment from the viewpoint of aerodynamic noise. With increasing λ, CT decreases monotonically from a small tip speed ratio (λ = 0.2). These tendencies are commonly found among drag-type wind turbines. The wind-tunnel experimental results by Shimizu et al. [24] for CP and CT of an O-VAWT with two elliptical cross-sectional blades show the same tendencies as our results. The value of the maximum CP is 0.32 in our experiments and CFD simulations, while the value is 0.176 in Shimizu et al.’s experiments. The main factors for the better performance of our O-VAWT as compared to Shimizu et al.’s O-VAWT can be the number of blades and the cross-sectional shape of the blades. In our previous studies, the maximum Cp improved from 0.189 to 0.244 by changing the number of blades from two to three [33] and improved from 0.246 to 0.288 by changing the cross-sectional shape of blades from ellipse to rectangle [23].
Figure 8b shows CT_rev and CT_rot computed by the CFD results. The sum of CT_rev and CT_rot is CT. As well as CT, the value of CT_rev decreases monotonically with an increase in λ. Conversely, the value of CT_rot increases with an increase in λ. The values of CT_rev are positive and larger than those of CT_rot except for λ = 0.8. At λ = 0.8, CT_rev is negative; however, CT_rot is positive and its absolute value is larger than that of CT_rev. As a result, CT is positive at λ = 0.8.
Figure 9 shows the variations of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and blade-rotation component (CTB_rot) with respect to azimuthal angle (φ) at λ = 0.4 and 0.6. The value of CTB is significantly large in the upwind region of the advancing side (φ ≈ 180° to 270°) of the rotor, being the maximum at φ ≈ 210°. In the range of φ where CTB is significantly large, the contribution of CTB_rev is dominant. Except for this range, CTB_rev does not always positively contribute to CTB. At φ where the value of CTB_rev is negative, the value of CTB_rot is generally positive and the rotation of the blade positively contributes to CTB. Due to this positive contribution of CTB_rot to CTB, the value of CTB is positive at almost all φ and the variation of CTB of the O-VAWT with respect to φ is smaller as compared to that of a Savonius-type VAWT. The variation of CTB of a Savonius-type VAWT with respect to φ in Figure 9 is a result of CFD simulation by Tian et al. [37]. The maximum CP and the optimal λ of the Savonius-type VAWT were 0.258 and 1.0, respectively.

4.2. Performance in Shear Flows

Figure 10 compares the experimental results for CP of the O-VAWT in the cases of the advancing side faster shear flow (ASF-SF), the retreating side faster shear flow (RSF-SF) and the uniform flow. In the cases of the ASF-SF, CP is higher than in the case of the uniform flow at all λ and increases with an increase in Γ. In the cases of the RSF-SF, similar to the cases of the ASF-SF, CP increases with an increase in Γ. However, when Γ = 0.28, the values of CP are lower than those of the uniform flow case. Concerning the optimal λ, there is a trend that it shifts to higher λ in both cases of shear flows as compared to the uniform flow case, except for the RSF-SF with Γ = 0.28.
Both in the cases of the ASF-SF (Figure 11a) and the RSF-SF (Figure 11b), the CFD results for Cpλ curves are in good agreement with the experimental ones. In the following discussion, we use the CFD results for the torque of the O-VAWT and the flow field to explain the effects of shear flows on the characteristics of the power performance.
Figure 12 shows the variations of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) with azimuthal angle (φ) in the cases of the ASF-SF. It is confirmed that these profiles are qualitatively the same between the cases of λ = 0.4 and λ = 0.6. Therefore, it is considered that the effects of the shear flow on the characteristics of the blade torque variations with φ do not significantly change around the optimal λ. As compared to the case of the uniform flow, CTB is higher on most of the advancing side of the rotor (φ ≈ 210° to 330°) and lower in the upwind region of the retreating side of the rotor, specifically at φ ≈ 120° to 150°. In particular, with an increase in Γ, CTB increases on most of the advancing side of the rotor (φ ≈ 210° to 330°). The optimal φ at which CTB is the maximum shifts to the downwind direction on the advancing side of the rotor (φ ≈ 240°) as compared to the case of the uniform flow. The effects of the shear flow on the variations of CTB_rev with φ is almost the same as CTB. As compared to the case of the uniform flow, CTB_rev is higher on most of the advancing side of the rotor (φ ≈ 210° to 330°). In contrast, CTB_rot is lower in most of the upwind region of the advancing side (φ ≈ 180° to 240°) and slightly higher in the upwind region of the retreating side of the rotor, specifically at 120° to 150° as compared to the case of the uniform flow. With an increase in Γ, CTB_rot decreases in the upwind region of the advancing side of the rotor (φ ≈ 180° to 240°). Due to its negative values of CTB_rot, the optimal φ at which CTB is the maximum slightly shifts to the downwind direction as compared to CTB_rev.
Figure 13 shows the variation of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) with azimuthal angle (φ) in the cases of the RSF-SF. Since these profiles are qualitatively the same between the cases of λ = 0.4 and λ = 0.6, it is considered that the effects of the shear flow on the characteristics of the blade torque variations with φ do not significantly change around the optimal λ. As compared to the case of the uniform flow, CTB is significantly higher in most of the upwind region of the retreating side (φ ≈ 120° to 180°) and lower on most of the advancing side of the rotor (φ ≈ 210° to 330°). In particular, with an increase in Γ, CTB increases in the upwind region of the retreating side of the rotor. The optimal φ at which CTB is the maximum shifts to the upwind region of the retreating side of the rotor (φ ≈ 150°) as compared to the case of the uniform flow. The effects of the shear flow on the variation of CTB_rev with φ is almost the same as CTB. As compared to the case of the uniform flow, CTB_rev is higher in most of the upwind region of the retreating side (φ ≈ 100° to 160°) and lower on most of the advancing side of the rotor (φ ≈ 210° to 330°). In contrast, CTB_rot is slightly lower in most of the upwind region of the retreating side (φ ≈ 90° to 150°) and slightly higher in most of the upwind region of the advancing side of the rotor (φ ≈ 180° to 240°) as compared to the case of the uniform flow. Furthermore, with an increase in Γ, CTB_rot decreases in most of the upwind region of the retreating side (φ ≈ 100° to 160°) and increases in most of the upwind region of the advancing side of the rotor (φ ≈ 180° to 240°).

4.3. Flow Characteristics

Figure 14a, Figure 15a and Figure 16a show the temporal sequence of the horizontal distributions of normalized horizontal velocity vectors, normalized pressure and normalized vorticity, respectively, at the mid-height of the O-VAWT at λ = 0.4 in the case of the uniform flow. At φ = 180° to 270°, the approaching flow to the blade has a large velocity component perpendicular to the blade, and the pressure on the upwind side of the blade is high. Due to this high pressure, CTB_rev is significantly high in the range of φ = 180° to 270° in Figure 9. In addition, at φ = 210° to 270°, due to the strong large vortex formed near the outer edge of the downwind side of the blade, the pressure is low near the vortex. This low pressure positively contributes to CTB_rev (see Figure 9 at φ = 210° to 270°). By contrast, this low pressure negatively contributes to CTB_rot (see Figure 9 at φ = 210° to 270°). At φ = 300° and 330°, the approaching flow to the blade has a small velocity component perpendicular to the blade and the pressure on the upwind side of the blade is not high. Therefore, CTB_rev is lower as compared to the upwind region of the advancing side of the rotor (see Figure 9 at φ = 180° to 270°). At φ = 0° to 120°, the attack angle of the blade is positive (here, the counterclockwise direction is defined as positive) and the flow separates over the outer side of the blade. Therefore, on the upwind edge of the blade and on the inner side of the blade near its upwind edge, the pressure is relatively high. In contrast, on the outer side of the blade near its upwind edge, the pressure is relatively lower. This pressure distribution contributes negatively to CTB_rev and positively to CTB_rot.
Figure 14b and Figure 15b and Figure 16b show the temporal sequence of the horizontal distributions of normalized horizontal velocity vectors, normalized pressure and normalized vorticity, respectively, at the mid-height of the O-VAWT at λ = 0.4 in the case of the ASF-SF with Γ = 0.51. At φ = 210° to 330°, the approaching flow to the blade has a larger velocity component perpendicular to the blade, and the pressure on the upwind side of the blade is higher as compared to the uniform flow case. Due to this higher pressure, CTB_rev is higher than the uniform flow case (see Figure 12 at φ = 210° to 330°). In addition, at φ = 210° to 240°, wind speed is significantly increased at the outer edge of the blade and a significantly stronger and larger vortex is formed near the outer edge on the downwind side of the blade. Near the vortex, the pressure is lower as compared to the uniform flow case. This low pressure contributes to higher CTB_rev and lower CTB_rot as compared to the uniform flow case (see Figure 12 at φ = 210° to 240°). At φ = 120° and 150°, the pressure on the outer side of the blade is lower due to the lower speed approaching flow to the blade. Furthermore, the pressure on the inner side of the blade near its upwind edge is higher due to the existence of the blade at φ = 240° or 270°, respectively. This pressure distribution contributes to lower CTB_rev and higher CTB_rot as compared to the uniform flow case (see Figure 12 at φ = 120° to 150°).
Figure 14c, Figure 15c and Figure 16c show the temporal sequence of normalized horizontal distributions of horizontal velocity vectors, pressure and vorticity, respectively, at the mid-height of the O-VAWT at λ = 0.4 in the case of the RSF-SF with Γ = 0.51. At φ = 120° and 150°, the approaching flow to the blade has a larger velocity component perpendicular to the blade and the pressure on the outer side of the blade is higher as compared to the uniform flow case. Due to this higher pressure, CTB_rev is higher than the uniform flow case (see Figure 13 at φ = 120° and 150°). Furthermore, wind speed is increased at the upwind edge of the blade and a stronger vortex is formed near the upwind edge of the blade on its inner side. Near the vortex, the pressure is lower as compared to the uniform flow case. This low pressure contributes to higher CTB_rev and lower CTB_rot as compared to the uniform flow case (see Figure 13 at φ = 120° and 150°). At φ = 210° to 300°, except in the vicinity of the upwind side of the inner edge of the blade at φ = 210°, the approaching flow to the blade has a smaller velocity component perpendicular to the blade and the pressure on the upwind side of the blade is lower as compared to the uniform flow case. Due to this lower pressure, CTB_rev is lower than the uniform flow case (see Figure 13 at φ = 210° and 300°). In addition, at φ = 210° and 240°, the vortex formed near the outer edge of the downwind side of the blade is weaker and the pressure drop becomes smaller as compared to the uniform flow case. This smaller pressure drop contributes to lower CTB_rev and higher CTB_rot as compared to the uniform flow (see Figure 13 at φ = 210° and 240°).

5. Conclusions

We investigated the effects of horizontal shear flow on the performance characteristics of an orthopter-type vertical axis wind turbine (O-VAWT) by conducting wind tunnel experiments and computational fluid dynamics (CFD) simulations. In addition to a uniform flow, two types of shear flow were used as the approaching flow to the O-VAWT. One type was an advancing side faster shear flow (ASF-SF), which had a higher velocity on the advancing side of the rotor. The other type was a retreating side faster shear flow (RSF-SF), which had a higher velocity on the retreating side of the rotor. For each type of shear flow, we set three different velocity ratios (Γ = 0.28, 0.40 and 0.51), which were the ratios of the difference between the highest velocity and the lowest velocity in a shear flow to the sum of the highest and lowest velocities. The main findings are summarized as follows:
  • In the ASF-SF cases, the power coefficients (CP) were significantly higher than the uniform flow case at all tip speed ratios (λ) and increased with Γ. The experimental results for the maximum CP of the ASF-SF case with Γ = 0.51 and the uniform flow case were 0.43 when λ = 0.6 and 0.32 when λ = 0.4, respectively. Around the optimal λ, the blade torque coefficient (CTB) on the advancing side of the rotor was, in general, significantly higher than the uniform flow case and increased with Γ, predominantly contributing to the increase in CP. The CFD results for the maximum discrepancies of CTB on the advancing side of the rotor between the ASF-SF case with Γ = 0.51 and the uniform flow case were 0.38 when λ = 0.4 and 0.44 when λ = 0.6. The high values of CTB of the ASF-SF cases on the advancing side of the rotor were mainly caused by the higher pressure on the upwind side of the blade due to the higher speed of the approaching flow and by the lower pressure near the outer edge of the downwind side of the blade due to the formation of a larger vortex.
  • In the RSF-SF cases, CP increased with Γ. However, when Γ = 0.28, CP was lower than the uniform flow case at all λ. When Γ = 0.51, CP was higher than the uniform flow case except at low λ; however, it was lower than the ASF-SF case with Γ = 0.28. The experimental results for the maximum CP of the RSF-SF case with Γ = 0.28, the RSF-SF case with Γ = 0.51 and the ASF-SF case with Γ = 0.28 were 0.29 when λ = 0.4, 0.33 when λ = 0.5 and 0.35 when λ = 0.5, respectively. Around the optimal λ, the blade torque coefficient (CTB) on the retreating side of the rotor was, in general, higher than the uniform flow case and increased with Γ, predominantly contributing to the increase in CP. The CFD results for the maximum discrepancies of CTB on the retreating side of the rotor between the RSF-SF case with Γ = 0.51 and the uniform flow case were 0.60 when λ = 0.4 and 0.36 when λ = 0.6. The high values of CTB of the RSF-SF cases on the retreating side of the rotor were mainly caused by the higher pressure on the outer side of the blade on the upwind side of the rotor, due to the higher speed of the approaching flow. By contrast, CTB on the advancing side of the rotor was, in general, lower than the uniform flow case, due to the lower pressure on the upwind side of the blade.
  • CTB consists of the rotor-revolution component (CTB_rev) and the blade-rotation component (CTB_rot). In all the shear flow cases, as well as the uniform flow case, the contributions of CTB_rev to CTB were dominant. The dependencies of CTB_rev and CTB_rot on Γ had the opposite tendencies.
These findings are useful for micro-siting of an O-VAWT in the area where shear flows occur. A location where ASF-SFs with high Γ values dominantly occur is ideal for installing the O-VAWT. At a location where not only ASF-SFs but also RSF-SFs occur at high frequencies, higher Γ values are preferable. However, the shear flows utilized in this study are limited in their profiles and the relative positions to the rotor. To properly conduct the micro-siting of an O-VAWT in the area where various kinds of shear flows occur, such as the vicinity of a building, it is essential to understand the performance characteristics of the O-VAWT in the various kinds of shear flows. Therefore, in future research, we plan to investigate the effects of the broadness of the shear layer and the relative position of the shear flow to the rotor on the O-VAWT’s performance characteristics. In addition, we plan to investigate the effects of the turbulence intensity of the approaching flow on the CFD results for the O-VAWT’s performance characteristics by setting obstacles, which emit eddies, upwind of the O-VAWT to avoid the rapid dissipation of high turbulence intensity.

Author Contributions

R.P.W. performed the numerical simulations and prepared this manuscript being supervised by T.K. (Takaaki Kono). All authors contributed to the analyses of the data. T.K. (Takaaki Kono) and T.K. (Takahiro Kiwata) supervised the entire work. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by World Bank Loan No. 8245-ID and a grant from Takahashi Industrial and Economic Research Foundation.

Acknowledgments

This research was supported by the Program for Research and Innovation in Science and Technology (RISET-Pro) Kemenristekdikti. The authors are thankful to technician Kuratani and student Taiki Sugawara for their help with the experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bahaj, A.S.; Myers, L.; James, P.A.B. Urban energy generation: Influence of micro-wind turbine output on electricity consumption in buildings. Energy Build. 2007, 39, 154–165. [Google Scholar] [CrossRef]
  2. Toja-Silva, F.; Colmenar-Santos, A.; Castro-Gil, M. Urban wind energy exploitation systems: Behaviour under multidirectional flow conditions—Opportunities and challenges. Renew. Sustain. Energy Rev. 2013, 24, 364–378. [Google Scholar] [CrossRef]
  3. Ishugah, T.F.; Li, Y.; Wang, R.Z.; Kiplagat, J.K. Advances in wind energy resource exploitation in urban environment: A review. Renew. Sustain. Energy Rev. 2014, 37, 613–626. [Google Scholar] [CrossRef]
  4. Tummala, A.; Velamati, R.K.; Sinha, D.K.; Indraja, V.; Krishna, V.H. A review on small scale wind turbines. Renew. Sustain. Energy Rev. 2016, 56, 1351–1371. [Google Scholar] [CrossRef]
  5. Kumar, R.; Raahemifar, K.; Fung, A.S. A critical review of vertical axis wind turbines for urban applications. Renew. Sustain. Energy Rev. 2018, 89, 281–291. [Google Scholar] [CrossRef]
  6. Stathopoulos, T.; Alrawashdeh, H.; Al-Quraan, A.; Blocken, B.; Dilimulati, A.; Paraschivoiu, M.; Pilay, P. Urban wind energy: Some views on potential and challenges. J. Wind Eng. Ind. Aerodyn. 2018, 179, 146–157. [Google Scholar] [CrossRef]
  7. Toja-Silva, F.; Kono, T.; Peralta, C.; Lopez-Garcia, O.; Chen, J. A review of computational fluid dynamics (CFD) simulations of the wind flow around buildings for urban wind energy exploitation. J. Wind Eng. Ind. Aerodyn. 2018, 180, 66–87. [Google Scholar] [CrossRef]
  8. Anup, K.C.; Whale, J.; Urmee, T. Urban wind conditions and small wind turbines in the built environment: A review. Renew. Energy 2019, 131, 268–283. [Google Scholar]
  9. Carbó Molina, A.; De Troyer, T.; Massai, T.; Vergaerde, A.; Runacres, M.C.; Bartoli, G. Effect of turbulence on the performance of VAWTs: An experimental study in two different wind tunnels. J. Wind Eng. Ind. Aerodyn. 2019, 193, 103969. [Google Scholar] [CrossRef]
  10. Wood, D. Small Wind Turbines; Green Energy and Technology; Springer: London, UK, 2011; ISBN 978-1-84996-174-5. [Google Scholar]
  11. Manwell, J.F.; McGowan, J.G.; Rogers, A.L. Wind Energy Explained: Theory, Design and Application; John Wiley & Sons: Hoboken, NJ, USA, 2010; ISBN 9780470015001. [Google Scholar]
  12. Allen, S.R.; Hammond, G.P.; McManus, M.C. Prospects for and barriers to domestic micro-generation: A United Kingdom perspective. Appl. Energy 2008, 85, 528–544. [Google Scholar] [CrossRef]
  13. Islam, M.; Ting, D.S.K.; Fartaj, A. Aerodynamic models for Darrieus-type straight-bladed vertical axis wind turbines. Renew. Sustain. Energy Rev. 2008, 12, 1087–1109. [Google Scholar] [CrossRef]
  14. 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]
  15. Ghasemian, M.; Ashrafi, Z.N.; Sedaghat, A. A review on computational fluid dynamic simulation techniques for Darrieus vertical axis wind turbines. Energy Convers. Manag. 2017, 149, 87–100. [Google Scholar] [CrossRef]
  16. Kumar, P.M.; Sivalingam, K.; Narasimalu, S.; Lim, T.-C.; Ramakrishna, S.; Wei, H. A review on the evolution of darrieus vertical axis wind turbine: Small wind turbines. J. Power Energy Eng. 2019, 7, 27–44. [Google Scholar] [CrossRef] [Green Version]
  17. Abraham, J.P.; Plourde, B.D.; Mowry, G.S.; Minkowycz, W.J.; Sparrow, E.M. Summary of Savonius wind turbine development and future applications for small-scale power generation. J. Renew. Sustain. Energy 2012, 4, 042703. [Google Scholar] [CrossRef]
  18. Akwa, J.V.; Vielmo, H.A.; Petry, A.P. A review on the performance of Savonius wind turbines. Renew. Sustain. Energy Rev. 2012, 16, 3054–3064. [Google Scholar] [CrossRef]
  19. Roy, S.; Saha, U.K. Review on the numerical investigations into the design and development of Savonius wind rotors. Renew. Sustain. Energy Rev. 2013, 24, 73–83. [Google Scholar] [CrossRef]
  20. Kang, C.; Liu, H.; Yang, X. Review of fluid dynamics aspects of Savonius-rotor-based vertical-axis wind rotors. Renew. Sustain. Energy Rev. 2014, 33, 499–508. [Google Scholar] [CrossRef]
  21. Alom, N.; Saha, U.K. Evolution and progress in the development of savonius wind turbine rotor blade profiles and shapes. J. Sol. Energy Eng. Trans. ASME 2019, 141, 1–15. [Google Scholar] [CrossRef]
  22. Elkhoury, M.; Kiwata, T.; Nagao, K.; Kono, T.; ElHajj, F. Wind tunnel experiments and Delayed Detached Eddy Simulation of a three-bladed micro vertical axis wind turbine. Renew. Energy 2018, 129, 63–74. [Google Scholar] [CrossRef]
  23. Kiwata, T. Vertical axis wind turbine with variable-pitch straight blades. In Proceedings of the International Conferemce on Jets, Wakes and Separated Flows, ICJWSF-2017, Cincinnati, OH, USA, 9–12 October 2017; pp. 1–6. [Google Scholar]
  24. Shimizu, Y.; Maeda, T.; Kamada, Y.; Takada, M.; Katayama, T. Development of micro wind turbine (Orthoptere Wind Turbine). In Proceedings of the International Conference on Fluid Engineering JSME Centennial Grand Congress, Tokyo, Japan, 13–16 July 1997; pp. 1551–1556. [Google Scholar]
  25. Bayeul-Laine, A.-C.; Simonet, S.; Dockter, A.; Bois, G. Numerical study of flow stream in a mini VAWT with relative rotating blades. In Proceedings of the 22nd International Symposium on Transport Phenomena Conference, Delft, The Netherlands, 8–11 November 2011; p. 13. [Google Scholar]
  26. Cooper, P.; Kennedy, O.C.; Cooper, P.; Kennedy, O. Development and analysis of a novel vertical axis wind turbine development and analysis of a novel vertical axis wind turbine. In Proceedings of the Solar 2004: Life, The Universe and Renewables, Perth, Australia, 30 November–3 December 2004; pp. 1–9. [Google Scholar]
  27. Kono, T.; Kogaki, T.; Kiwata, T. Numerical investigation ofwind conditions for roof-mounted wind turbines: Effects of wind direction and horizontal aspect ratio of a high-rise cuboid building. Energies 2016, 9, 907. [Google Scholar] [CrossRef] [Green Version]
  28. Balduzzi, F.; Bianchini, A.; Ferrari, L. Microeolic turbines in the built environment: Influence of the installation site on the potential energy yield. Renew. Energy 2012, 45, 163–174. [Google Scholar] [CrossRef]
  29. Ferreira, C.J.S.; Van Bussel, G.J.W.; Van Kuik, G.A.M. Wind tunnel hotwire measurements, flow visualization and thrust measurement of a VAWT in skew. J. Sol. Energy Eng. Trans. ASME 2006, 128, 487–497. [Google Scholar] [CrossRef]
  30. Burton, T.; Jenkins, N.; Sharpe, D.; Bossanyi, E. Wind Energy Handbook, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011; ISBN 9780470699751. [Google Scholar]
  31. ANSYS, Inc. ANSYS Fluent Theory Guide, Release 17.2; ANSYS, Inc.: Canonsburg, PA, USA, 2016. [Google Scholar]
  32. ANSYS, Inc. ANSYS Fluent User Guide, Release 17.2; ANSYS, Inc.: Canonsburg, PA, USA, 2016. [Google Scholar]
  33. ElCheikh, A.; Elkhoury, M.; Kiwata, T.; Kono, T. Performance analysis of a small-scale orthopter-type vertical axis wind turbine. J. Wind Eng. Ind. Aerodyn. 2018, 180, 19–33. [Google Scholar] [CrossRef]
  34. Spalart, P.R.; Jou, W.H.; Strelets, M.; Allmaras, S.R. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In Proceedings of the First AFOSR International Conference on DNS/LES, Ruston, LA, USA, 4–8 August 1997; pp. 137–147. [Google Scholar]
  35. Tominaga, Y.; Mochida, A.; Yoshie, R.; Kataoka, H.; Nozu, T.; Yoshikawa, M.; Shirasawa, T. AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings. J. Wind Eng. Ind. Aerodyn. 2008, 96, 1749–1761. [Google Scholar] [CrossRef]
  36. Wu, H.; Stathopoulos, T. Wind-tunnel techniques for assessment of pedestrian-level winds. J. Eng. Mech. 1993, 119, 1920–1936. [Google Scholar] [CrossRef]
  37. Tian, W.; Mao, Z.; Zhang, B.; Li, Y. Shape optimization of a Savonius wind rotor with different convex and concave sides. Renew. Energy 2018, 117, 287–299. [Google Scholar] [CrossRef]
Figure 1. Illustration of the approaching wind flow to the wind turbine near upwind corners of a building. (a) Bird view; (b) enlarged top view.
Figure 1. Illustration of the approaching wind flow to the wind turbine near upwind corners of a building. (a) Bird view; (b) enlarged top view.
Applsci 10 01778 g001
Figure 2. Orthopter-type vertical-axis wind turbine (O-VAWT). (a) A photograph of O-VAWT with three flat blades, main shaft, arm, the chains and sprockets; (b) an isometric view of O-VAWT with three flat blades; (c) motion of rotor and blades viewed from the top; (d) a projected swept area of the rotor viewed from the upwind side.
Figure 2. Orthopter-type vertical-axis wind turbine (O-VAWT). (a) A photograph of O-VAWT with three flat blades, main shaft, arm, the chains and sprockets; (b) an isometric view of O-VAWT with three flat blades; (c) motion of rotor and blades viewed from the top; (d) a projected swept area of the rotor viewed from the upwind side.
Applsci 10 01778 g002
Figure 3. The experimental apparatus and measurement devices; (a) in uniform flow, (b) in shear flows and (c) the porous plate position at the nozzle exit of the wind tunnel in case of shear flows.
Figure 3. The experimental apparatus and measurement devices; (a) in uniform flow, (b) in shear flows and (c) the porous plate position at the nozzle exit of the wind tunnel in case of shear flows.
Applsci 10 01778 g003
Figure 4. The perforated panel and splitter plate set at the outlet of the wind tunnel.
Figure 4. The perforated panel and splitter plate set at the outlet of the wind tunnel.
Applsci 10 01778 g004
Figure 5. Computational domain, computational mesh and boundary conditions. (a) Top view of the computational domain; (b) bird’s eye view of the computational domain with computational mesh and boundary conditions; (c) modeled rotor; (d) blade domains and (e) rotor domain.
Figure 5. Computational domain, computational mesh and boundary conditions. (a) Top view of the computational domain; (b) bird’s eye view of the computational domain with computational mesh and boundary conditions; (c) modeled rotor; (d) blade domains and (e) rotor domain.
Applsci 10 01778 g005
Figure 6. Horizontal distribution of the time-mean values of streamwise velocity at the mid-height of the wind tunnel outlet at x ≈ −0.147D, which corresponds to 0.1 m downwind of the wind tunnel outlet, and x = 0, which corresponds to the position of the rotational axis of the O-VAWT, when the O-VAWT is absent. (a) ASF-SF cases at x ≈ −0.147D, (b) RSF-SF cases at x ≈ −0.147D, (c) ASF-SF cases at x = 0 and (d) RSF-SF cases x = 0.
Figure 6. Horizontal distribution of the time-mean values of streamwise velocity at the mid-height of the wind tunnel outlet at x ≈ −0.147D, which corresponds to 0.1 m downwind of the wind tunnel outlet, and x = 0, which corresponds to the position of the rotational axis of the O-VAWT, when the O-VAWT is absent. (a) ASF-SF cases at x ≈ −0.147D, (b) RSF-SF cases at x ≈ −0.147D, (c) ASF-SF cases at x = 0 and (d) RSF-SF cases x = 0.
Applsci 10 01778 g006
Figure 7. Horizontal distribution of the turbulence intensity of ASF-SF cases at the mid-height of the wind tunnel outlet at x = 0, which corresponds to the position of the rotational axis of the O-VAWT, when the O-VAWT was absent. (a) Linear scale on the horizontal axis and (b) logarithmic scale in the horizontal axis.
Figure 7. Horizontal distribution of the turbulence intensity of ASF-SF cases at the mid-height of the wind tunnel outlet at x = 0, which corresponds to the position of the rotational axis of the O-VAWT, when the O-VAWT was absent. (a) Linear scale on the horizontal axis and (b) logarithmic scale in the horizontal axis.
Applsci 10 01778 g007
Figure 8. Performance of the O-VAWT in the uniform flow; (a) the variation of power coefficient (CP) and torque coefficient (CT) with tip speed ratios (λ) by the experiments and CFD simulations, (b) rotor-revolution (CT_rev) and blade-rotation (CT_rot) components of torque coefficients (CT) computed by the CFD simulations. The wind-tunnel experimental results by Shimizu et al. [24] for CP and CT of an O-VAWT with two elliptical cross-sectional blades are added for reference.
Figure 8. Performance of the O-VAWT in the uniform flow; (a) the variation of power coefficient (CP) and torque coefficient (CT) with tip speed ratios (λ) by the experiments and CFD simulations, (b) rotor-revolution (CT_rev) and blade-rotation (CT_rot) components of torque coefficients (CT) computed by the CFD simulations. The wind-tunnel experimental results by Shimizu et al. [24] for CP and CT of an O-VAWT with two elliptical cross-sectional blades are added for reference.
Applsci 10 01778 g008
Figure 9. The variation of blade torque coefficients (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) with azimuthal angle (φ) at; (a) λ = 0.4 and (b) λ = 0.6. Note that CTB, CTB_rev and CTB_rot are coefficients of one blade. The CFD simulation result by Tian et al. [37] for the variation of CTB of a Savonius-type VAWT with φ is added for comparison.
Figure 9. The variation of blade torque coefficients (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) with azimuthal angle (φ) at; (a) λ = 0.4 and (b) λ = 0.6. Note that CTB, CTB_rev and CTB_rot are coefficients of one blade. The CFD simulation result by Tian et al. [37] for the variation of CTB of a Savonius-type VAWT with φ is added for comparison.
Applsci 10 01778 g009
Figure 10. Power coefficient (CP) of O-VAWT in shear flows by the experiments.
Figure 10. Power coefficient (CP) of O-VAWT in shear flows by the experiments.
Applsci 10 01778 g010
Figure 11. Power coefficient (CP) of O-VAWT in shear flows by the experiments and the CFD simulations in the cases of; (a) the ASF-SF and (b) the RSF-SF.
Figure 11. Power coefficient (CP) of O-VAWT in shear flows by the experiments and the CFD simulations in the cases of; (a) the ASF-SF and (b) the RSF-SF.
Applsci 10 01778 g011
Figure 12. The variation of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) in the cases of the ASF-SF computed by the CFD simulations; (a) for λ = 0.4 and (b) for λ = 0.6. Note that CTB, CTB_rev and CTB_rot are coefficients of one blade.
Figure 12. The variation of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) in the cases of the ASF-SF computed by the CFD simulations; (a) for λ = 0.4 and (b) for λ = 0.6. Note that CTB, CTB_rev and CTB_rot are coefficients of one blade.
Applsci 10 01778 g012
Figure 13. The variation of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) in the cases of the RSF-SF computed by the CFD simulations; (a) for λ = 0.4 and (b) for λ = 0.6. Noted that CTB, CTB_rev and CTB_rot are coefficients of one blade.
Figure 13. The variation of blade torque coefficient (CTB), its rotor-revolution component (CTB_rev) and its blade-rotation component (CTB_rot) in the cases of the RSF-SF computed by the CFD simulations; (a) for λ = 0.4 and (b) for λ = 0.6. Noted that CTB, CTB_rev and CTB_rot are coefficients of one blade.
Applsci 10 01778 g013
Figure 14. The temporal sequence of the horizontal distributions of normalized horizontal velocity vectors of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Figure 14. The temporal sequence of the horizontal distributions of normalized horizontal velocity vectors of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Applsci 10 01778 g014
Figure 15. The temporal sequence of normalized pressure at the mid-height of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Figure 15. The temporal sequence of normalized pressure at the mid-height of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Applsci 10 01778 g015
Figure 16. The temporal sequence of normalized vorticity at the mid-height of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Figure 16. The temporal sequence of normalized vorticity at the mid-height of the O-VAWT for λ = 0.4 in the case of; (a) the uniform flow, (b) the ASF-SF with Γ = 0.51 and (c) the RSF-SF with Γ = 0.51.
Applsci 10 01778 g016
Table 1. Perforated panels. Here, d is the diameter of a hole, L is the distance between the centers of adjacent holes and Φ is the shielding ratio.
Table 1. Perforated panels. Here, d is the diameter of a hole, L is the distance between the centers of adjacent holes and Φ is the shielding ratio.
Named [m]L [m]Φ [-]Enlarged View
Perforated panel A3 × 10−34 × 10−30.49 Applsci 10 01778 i001
Perforated panel B3 × 10−34.5 × 10−30.60 Applsci 10 01778 i002
Perforated panel C3 × 10−35 × 10−30.67 Applsci 10 01778 i003
Table 2. Numbers of grid sizes.
Table 2. Numbers of grid sizes.
DomainMesh on This ResearchRegular Mesh Used by ElCheikh [33]
Blade 12,300,5301,610,660
Blade 22,303,7251,614,903
Blade 32,304,6071,588,363
Blade 4-1,614,661
Rotor1,756,6352,127,497
Far end1,770,583441,196
TOTAL10,436,0808,997,280
Table 3. Distribution of u at the inlet boundary. Here, Γ is the velocity ratio; Φ is the shielding ratio of a perforated panel; UH and UL are the maximum and minimum velocities in a shear flow.
Table 3. Distribution of u at the inlet boundary. Here, Γ is the velocity ratio; Φ is the shielding ratio of a perforated panel; UH and UL are the maximum and minimum velocities in a shear flow.
Flow TypeΓ Φ UH [m/s]Range of UH Region [m]Velocity Distribution in Transition RegionRange of Transition Region [m]UL [m/s]Range of UL Region [m]
Uniform0-8 8 8
ASF-SF 0.510.6713.2y ≤ −0.24UH ∙(−y/0.24)0.3−0.24 ≤ y ≤ −0.014.3y ≥ −0.01
ASF-SF0.400.6012.1y ≤ −0.22UH ∙(−y/0.22)0.35−0.22 ≤ y ≤ 0.025.2y ≥ −0.02
ASF-SF0.280.4910.8y ≤ −0.18UH ∙(−y/0.18)0.2−0.18 ≤ y ≤ −0.016.1y ≥ −0.01
RSF-SF 0.510.6713.2y ≥ 0.24UH ∙(y/0.24)0.30.01 ≤ y ≤ 0.244.3y ≤ 0.01
RSF-SF0.400.6012.1y ≥ 0.22UH ∙(y/0.24)0.350.02 ≤ y ≤ 0.225.2y ≤ 0.02
RSF-SF0.280.4910.8y ≥ 0.18UH ∙(y/0.24)0.20.01 ≤ y ≤ 0.186.1y ≤ 0.01
Table 4. Experimental and computational fluid dynamics (CFD) results for the values of U0 and Re for the uniform flow, ASF-SF and RSF-SF cases.
Table 4. Experimental and computational fluid dynamics (CFD) results for the values of U0 and Re for the uniform flow, ASF-SF and RSF-SF cases.
Flow TypeΓU0 [m/s]Re
ExperimentCFDExperimentCFD
Uniform0882.79 × 1052.79 × 105
ASF-SF0.5110.089.783.52 × 1053.41 × 105
ASF-SF0.409.389.223.28 × 1053.22 × 105
ASF-SF0.288.738.723.05 × 1053.05 × 105
RSF-SF0.517.767.482.71 × 1052.61 × 105
RSF-SF0.407.727.362.70 × 1052.57 × 105
RSF-SF0.287.527.332.62 × 1052.56 × 105

Share and Cite

MDPI and ACS Style

Wijayanto, R.P.; Kono, T.; Kiwata, T. Performance Characteristics of an Orthopter-Type Vertical Axis Wind Turbine in Shear Flows. Appl. Sci. 2020, 10, 1778. https://doi.org/10.3390/app10051778

AMA Style

Wijayanto RP, Kono T, Kiwata T. Performance Characteristics of an Orthopter-Type Vertical Axis Wind Turbine in Shear Flows. Applied Sciences. 2020; 10(5):1778. https://doi.org/10.3390/app10051778

Chicago/Turabian Style

Wijayanto, Rudi Purwo, Takaaki Kono, and Takahiro Kiwata. 2020. "Performance Characteristics of an Orthopter-Type Vertical Axis Wind Turbine in Shear Flows" Applied Sciences 10, no. 5: 1778. https://doi.org/10.3390/app10051778

APA Style

Wijayanto, R. P., Kono, T., & Kiwata, T. (2020). Performance Characteristics of an Orthopter-Type Vertical Axis Wind Turbine in Shear Flows. Applied Sciences, 10(5), 1778. https://doi.org/10.3390/app10051778

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