# Fish schooling as a basis for vertical axis wind turbine farm design

## Abstract

Most wind farms consist of horizontal axis wind turbines (HAWTs) due to the high power coefficient (mechanical power output divided by the power of the free-stream air through the turbine cross-sectional area) of an isolated turbine. However when in close proximity to neighbouring turbines, HAWTs suffer from a reduced power coefficient. In contrast, previous research on vertical axis wind turbines (VAWTs) suggests that closely-spaced VAWTs may experience only small decreases (or even increases) in an individual turbine’s power coefficient when placed in close proximity to neighbours, thus yielding much higher power outputs for a given area of land. A potential flow model of inter-VAWT interactions is developed to investigate the effect of changes in VAWT spatial arrangement on the array performance coefficient, which compares the expected average power coefficient of turbines in an array to a spatially-isolated turbine. A geometric arrangement based on the configuration of shed vortices in the wake of schooling fish is shown to significantly increase the array performance coefficient based upon an array of wind turbines. Results suggest increases in power output of over one order of magnitude for a given area of land as compared to HAWTs.

## 1 Introduction

Technology for extracting energy from the wind has become a topic of increasing political, economic and social interest in the last fifty years [1, 2]. Nearly all current wind farms use horizontal axis wind turbines (HAWTs) due to their high power coefficient, defined as the mechanical power output of the turbine divided by the power of the free-stream air through the cross-sectional area of the rotor normal to the free-stream direction. However, a major constraint on wind farms is the inter-turbine spacing, since the performance of each HAWT decreases with decreased spacing [3, 4, 5, 6]. When a HAWT extracts energy from the free stream, the speed of the flow behind the turbine is reduced. A well-known aerodynamic result is that the maximum amount of power a rotor can extract is proportional to the cube of the free-stream speed [4]. Horizontal axis wind turbines downstream of other HAWTs have a reduced incoming flow speed and thus are capable of extracting significantly less power than that of a spatially-isolated turbine. Theoretical analysis, computational simulations, and experimental results indicate that decreasing the downstream spacing of a row of multiple HAWTs from ten rotor-diameters to five rotor-diameters can decrease the power generated by roughly 40% [7, 5, 8]. Consequently, most modern HAWT wind farms space turbines five to ten rotor-diameters downstream and three or more rotor-diameters laterally [9, 4].

An alternative paradigm for wind energy extraction is found in vertical axis wind turbine (VAWT) designs. Studies have shown that a single spatially-isolated VAWT has a significantly lower power coefficient and costs approximately a third more than a HAWT with comparable power output [4]. Current literature regarding the performance of VAWT wind farms is scarce when compared to HAWTs, yet some sources indicate that the typical spacing constraints on HAWT wind farms do not apply [10, 11, 12]. Furthermore, a recent U.S. Patent suggests that aerodynamic interactions between a pair of counter-rotating VAWTs contribute to a higher power coefficient per turbine [13]. Previous work by Rajagopalan, Klimas, and Rickerl suggests an increase of up to 4% in the average turbine power coefficient of an array of VAWTs due to steam-tube contraction (induced higher flow speeds around neighbouring VAWTs) [12]. In contrast, Schatzle, Klimas, and Spahr show no change or a mild reduction in power coefficient for a pair of closely-spaced VAWTs [10]. However, neither of these studies considered the possibility of alternating the directions of rotation of neighbouring VAWTs.

For this work, we consider a VAWT array configuration following the arrangement of shed vortices in the wake of schooling fish. These shed vortices form a reversed Kármán vortex street in the wake of the fish [14]. Previous work by Weihs [15] uses a potential flow model to analyze the wake effects of the shed vortices on a school of fish. Motivated by the demonstrated benefits of the reversed Kármán vortex street on the propulsion of the schooling fish, we apply the same configuration and similar modelling tools to analyze VAWT arrays.

We investigate the inter-turbine spacing effects of VAWT arrays through a semi-empirical model. Models are developed, given in section 2, to describe a single-VAWT and arrays of VAWTs using parameters from the field measurements. The field measurements of an installed VAWT are presented in section 3. The results of fitting the model parameters to the field measurements and evaluating effects of array configuration on power output are given in section 4 and discussed in section 5.

## 2 Modelling

### 2.1 Potential Flow Model

In the interest of examining the effects of spacing and orientation of VAWTs in an array, a potential flow model is first constructed for a single turbine. The model consists of uniform flow with a dipole and a point vortex superimposed to represent one VAWT. The point vortex accounts for the flow rotation caused by the turbine and the dipole represents flow blockage due to the presence of the turbine in the flow. The complex velocity potential, , and the velocity vector field, , are given by:

(1) |

(2) |

where refers to the free-stream speed, is the strength of the point vortex, is the strength of the dipole, and are the unit vectors in the and directions, respectively, and ( and are real numbers). To create the velocity potential corresponding to an array of VAWTs, the model in equation 1 is summed over the VAWT positions in the array,

(3) |

where corresponds to the location of the th VAWT. The velocity field is obtained using equation 2. To account for a velocity-deficit in each VAWT wake, the velocity-deficit curves given by Hau [4] are applied:

(4) |

where is the resulting velocity vector after accounting for wake-effects and is a spatial function representing the velocity deficit. The function is a normal probability density function (PDF) in the lateral direction and a beta PDF in the downstream direction. The term adjusts the magnitude of the downstream velocity to reflect the velocity deficit in the wake. The wake extends approximately six rotor-diameters downstream and 20 bilaterally from the free-stream wind direction. This velocity deficit is included for each VAWT in the array. Accounting for the wake in this manner implies that is no longer a solution to the potential flow equations.

In order to evaluate the expected power produced by a turbine in an array, we define a parameter representative of the power output. By considering a simple drag-based VAWT, the tangential force on each turbine blade is , where is the velocity tangent to the blade motion. Furthermore, we assume that angular velocity of the rotor is . Putting these together, the model for the power output of a single, spatially-isolated VAWT can be represented by:

(5) |

where is the power output of the turbine, is the density of air, is a non-dimensional parameter specific to a VAWT design, is the radius of the wind turbine, is the unit vector tangent to the contour of integration, and is the contour of integration that is a circle with radius centered on the VAWT in the direction of rotation. The non-dimensional parameter is included to allow for the variety of wind turbine designs and accounts for effects such as blade shape, height of turbine, and other factors that affect wind turbine performance. The same expression for (with a different ) is obtained for a lift-based VAWT model as the lift is also and it is still reasonable to assume the angular velocity of the rotor is .

Having defined we introduce the array performance coefficient, , to evaluate the effects of closely spaced turbines in an array:

(6) |

where is the average performance parameter for an array. It is implicitly assumed that all turbines in an array are homogeneous (i.e. same ). For , the average VAWT in the array performs worse than an equivalent, spatially-isolated VAWT. However, if , the model suggests that the array configuration provides beneficial, inter-VAWT, aerodynamic interactions and increases the average power output of the array.

Although comparing is a means for understanding the flow physics, as a means to compare the performance of an array configuration in terms of power output per unit land area occupied by the VAWT farm, a power density parameter, , is defined:

(7) |

where is the land area of the entire wind farm, is the land area of a single turbine, taken to be arbitrarily , and is the number of turbines in the array. This parameter is most useful in the practical design of a wind farm and determining the appropriate spacing for maximizing the power output per a given area.

### 2.2 BioInspiration

A mathematical model similar to that described in the previous section is used by Weihs [15] to study the wake effects of fish schooling arrangements on fish propulsion. Weihs’ potential flow model uses uniform flow and superimposes point-vortices at the location of shed vortices from fish as they swim, depicted in figure a. The location of these shed vortices can be described in the complex plane by

(8) | ||||

(9) |

(10) | ||||

(11) |

where and correspond to positive (anticlockwise rotating) vortices, and correspond to negative (clockwise rotating) vortices, is the downstream distance between two vortices of the same line, is the lateral distance between the two lines of shed vortices from a given fish, is the distance between two adjacent fish, and and are integers. Fish in a school typically swim with , depending on species, and 2 [15]. Incorporating equation 9 and equation 11 into equation 3, the resulting velocity potential for an array of VAWTs following this particular arrangement is

(12) | |||||

### 2.3 Numerical Simulation

equation 12 is evaluated numerically for a VAWT array configuration where , corresponding to an array of turbines (256 turbines total). The values of , , and are obtained from the field measurements, as described in section 3.1. The model is interrogated for the range of parameters: , -, and , where corresponds to the turbine diameter. The range of is chosen to extend from turbines immediately downstream of each other to where the turbines have sufficiently cleared the downstream wake of each other. The range of is chosen to avoid overlapping turbines.

## 3 Field Measurements

Field measurements of a Savonius-type VAWT located in Glendora, CA were collected in spring 2009 in order to generate input data for the potential flow model. The turbine was a prototype model with rotor radius of m and rotor height of m. In order to characterize the flow around the VAWT, two cup-type flow anemometers were used. A reference anemometer consisted of an adjustable-length pole with an anemometer (InSpeed Vortex, with a range of 1.4 - 55.6 m s and m s accuracy) mounted on the end and connected to a dedicated computer that recorded the wind speed every second. The pole was fixed to the ground using guy wires for untended operation. The data from the reference anemometer was averaged over the entire visit to produce the free-stream speed during the field visit.

The second anemometer was mounted to a second pole that was manually positioned at 18 locations around the VAWT, as indicated in figure 2. In addition to measuring wind speed, the second anemometer contained a wind vane (InSpeed E-Vane) that measured the direction of the incoming wind ( accuracy). The wind vane was wired to a computer data acquisition system (Measurement Computing USB-1208FS) which recorded the voltage output of the wind vane at 100 Hz. This was then down-sampled through averaging the 100 measurements in one-second to obtain 1 Hz recordings that synchronized with the wind speed data. The anemometer recorded wind data at each station for approximately three minutes. For our field measurements, both instrument setups were extended to the mid-height of the VAWT and were adjusted according to the underlying terrain ( m). At the end of the measurements, the mobile instrument setup was moved near the reference setup in order to obtain the free-stream direction, which was assumed constant during the one-hour measuring period. This is a reasonable assumption as local weather stations report only westerly wind during our measurement time period. Data from both anemometers was synchronized for post-processed based on the recorded time-stamps.

### 3.1 Data Assimilation

The model requires the calculation of three parameters to describe a single, spatially-isolated VAWT: , , and . To obtain from the experimental measurements, the average free-stream speed from the reference anemometer is used. Using this measured free-stream velocity, the values of and are determined from a least-squares fit of a single-turbine model (one point vortex and one dipole) to the measured VAWT flow field.

## 4 Results

The measured flow field from the field visit to the VAWT along with the flow field calculated from equation 1, using the model parameters described in section 3.1, are given in figure 2. The measured average free-stream speed, , is 3.00 m s and appears to be well-described by a normal distribution having a standard deviation of 0.86 m s. The fitting of the measured flow field to the single-VAWT model, as given in equation 1, by least-squared difference yielded a unimodal minimum within the parameter space of and . The measured results are well-represented by a range of values given by =7.413.82 m s and = 0.000.77 m s. Despite the range of , and , in this letter we only present results using m s, = 7.41 m s, and m s.

The calculated for the fish-schooling VAWT-array models, given by equation 12 and using the range of values of and described in section 2.3, is shown in figure a. There is a clear region for which is greater than unity, corresponding to close-spaced pairs of counter-rotating VAWTs. The maximum in this parameter space suggests the potential for inter-turbine interactions that yield more power output than spatially-isolated turbines. figure b shows the normalized power density for the fish-schooling VAWT-array models. Comparing figures a and b shows that the optimal array configuration in terms of power extracted per unit land area does not coincide with the largest values, which implies improving the average performance of individual turbines might not be the desired goal of wind farm design. In figure 3, the dashed lines show the permissible parameter values for schooling fish. As seen in figure a, this region is below unity in the results. Array configurations of turbines with and turbines were also explored and show results similar to the presented result.

## 5 Discussion

A qualitative comparison between the flow fields from the VAWT model and the collected field data, as shown in figure 2, demonstrates that the current single-VAWT model captures the dominant flow characteristics for a single turbine. This verifies the model’s ability to reduce a complicated flow field to a three-parameter model. However from the discrepancy between velocity vectors of the single-VAWT model and collected field data in figure 2, it is evident that the wake of a spatially-isolated turbine is not accurately captured. These discrepancies motivated the introduction of an empirically-based wake from HAWTs equation 4, which is considered to be an appropriate means of accounting for an expected velocity deficit in the wake when calculating .

Measurements were limited to the area immediately surrounding the VAWT. Since the dominant flow effects of the VAWT are most pronounced in the closer locations, the vortex and dipole contributions to the velocity field drop off as and , respectively, where is the distance from the vortex or dipole center. This explains our model’s ability to capture the measured flow field despite the aforementioned limitation. Furthermore, measurements farther from the VAWT would likely reflect ambient fluctuations to a greater degree more than VAWT aerodynamics.

In order to understand the practical implications of this work for wind farm design, a comparison of the array power density (total array power output over array area) between VAWTs and HAWTs is made, as given by . For this analysis, we assume that the parameters, and , are the same for all VAWT designs. Additionally, we restrict ourselves to the parameter space representative of schooling fish, and chose the array configuration that maximizes following this constraint based on figure b. figure 4 demonstrates that significantly higher power densities are obtained with VAWT farms compared to operational HAWT farms, with increases of over one order of magnitude. Although close spacing of VAWTs potentially leads to higher power outputs of the individual turbines (regions of in figure a), the major advantage of VAWT wind farms is their significantly smaller inter-turbine spacing, compared to HAWT wind farms, without a substantial decrease in performance as demonstrated by figures a and b.

The proposed array configuration adapted from the geometry of fish schooling provides an interesting arrangement for study. figure a indicates decreases in the average turbine power output if turbines are geometrically arranged following the same configuration of shed vortices in the wake of schooling fish. However the increase in array power density of close-spaced arrays overcomes these mild reductions in expected power output, as seen in figure 4. If one allows for a broader range of parameter space than seen in schooling fish, there is the opportunity for significant improvements () in the power output of these turbines. This is not entirely surprising due to differences in the criteria of beneficial arrangements. The fish aim to align themselves to optimize their forward propulsion whereas spatial configurations for turbine arrays aim to maximize energy extraction.

While only one case is presented here, that of a VAWT array, two other finite arrays were explored with dimensions of and turbines. These did not show any significant deviation from the case presented, and thus explains their omission. For all turbine array sizes explored, the wind farm configurations corresponding to the highest values of , shown in figure a, consist of pairs of closely-spaced, counter-rotating VAWTs with their adjoining axis perpendicular to the free stream. A streamline analysis shows that stream-tube contraction between pairs of counter-rotating VAWTs, and the concomitant flow acceleration, is primarily responsible for improvements in expected power output. This is consistent with Rajagopalan et al. [12], who studied Darrieus-type VAWT arrays and showed a 4% increase in power coefficient due to the same mechanism. From our results it appears as though vortex interaction is the dominating mechanism for the stream-tube contraction, rather than blockage effects from neighbouring wind turbines. However, with further study, it is possible that this conclusion is dependent on VAWT design or array configuration, thus allowing VAWT designers to consider blockage effects in future design iterations.

## 6 Conclusion

In this work we present the experimental data used to create a potential flow model to describe the flow field around a spatially-isolated VAWT. This single-VAWT model consisted of three parameters: free-stream speed, vortex strength, and dipole strength. A model for an array of VAWTs, based on the single-VAWT model, is developed and an empirical wake is introduced to calculate the expected power output of the array. To evaluate the potential benefits of a particular array configuration, an array performance coefficient, , is defined to compare the average power output of an array configuration to that of a spatially-isolated VAWT. Configurations inspired by the arrangement of shed vortices in the wake of schooling fish are used for interrogating values with varied geometric parameters. Results show the possibility of values up to 1.4, for the VAWT array, when not restricted to parameters corresponding to actual schools of fish. The leading mechanism for this increase is stream-tube contraction and concomitant flow acceleration between counter-rotating turbines in close proximity. However, values of are approximately or less than unity for configurations following the schooling of fish regardless of array dimensions considered. These configurations significantly reduced the land use for VAWT wind farms, resulting in array power density increases of over one order of magnitude compared to operational HAWT wind farms.

Although the current model does not explicitly allow for complicated flow behavior, such as vortex shedding, turbulence, and three-dimensional effects, the current vortex and dipole model yields a flow field sufficient for assessing the performance of VAWTs. Furthermore, the field data is only from a small segment of time (one hour for the particular visit reported) and does not fully describe the variety of conditions in which VAWTs can produce power. It is of interest to conduct further measurements on a broad range of full-scale VAWT designs and over different operating conditions in order to refine the single-VAWT model. Similar measurements and considerations are required to refine the model for arrays of full-scale VAWTs. Such refinements include investigation of the dependence of and on VAWT design, wind speed, and surrounding environment (e.g. turbulence and the atmospheric boundary layer). Including these effects may potentially impact the results in terms of expected power output and array power density.

## Acknowledgments

We thank J Meier and P Abad-Manterola for assisting with the field work. We also thank California Energy & Power for allowing us the use of their wind turbine. This work has been funded by the National Science Foundation: Energy for Sustainability Program. The technology described in this paper is protected under both US and international pending patents filed by the California Institute of Technology.

### Footnotes

- The technology described in this paper is protected under both US and international pending patents filed by the California Institute of Technology.

### References

- T. Ackermann. Wind power in power systems. John Wiley and Sons, 2005.
- 20% wind energy by 2030: Increasing wind energy’s contribution to u.s. electricity supply, July 2008.
- PBS Lissaman. Energy effectiveness of arbitrary arrays of wind turbines. AIAA Paper, 79:323–328, 1979.
- Erich Hau. Wind Turbines. Springer, 2nd edition, 2006.
- M.K. Liu, MA Yocke, and TC Myers. Mathematical model for the analysis of wind-turbine wakes. Journal of Energy, 7:73–78, 1983.
- RJ Barthelmie, O. Rathmann, ST Frandsen, KS Hansen, E. Politis, J. Prospathopoulos, K. Rados, D. Cabezón, W. Schlez, J. Phillips, et al. Modelling and measurements of wakes in large wind farms. In Journal of Physics: Conference Series, volume 75, page 012049 (9pp). Institute of Physics Publishing, 2007.
- NO Jensen. A note on wind generator interaction. Risø-M-2411, Risø National Laboratory, Roskilde, Denmark, 1984.
- S. Ivanell, J.N. Sørensen, and D. Henningson. Numerical computations of wind turbine wakes. pages 259–263, 2007.
- B. Sørensen. Renewable energy: its physics, engineering, use, environmental impacts, economy, and planning aspects. Academic Press, 2004.
- PR Schatzle, PC Klimas, and HR Spahr. Aerodynamic interference between two Darrieus wind turbines. Wind Energy, pages 8–13, 1980.
- T. Nguyen. A vortex model of the Darrieus turbine. Master’s thesis, Texas Tech. University, Jan 1978.
- R. Rajagopalan, P. Klimas, and T. Rickerl. Aerodynamic interference of vertical axis wind turbines. Journal of Propulsion and Power, 6:645–653, 1990.
- R. N. Thomas. Patent #6,784,566 coupled vortex vertical axis wind turbine. United States Patent, 2004.
- M. S. Triantafyllou, G. S. Triantafyllou, and D. K. P. Yue. Hydrodynamics of fishlike swimming. Annual Review of Fluid Mechanics, 32:33–53, 2000.
- Daniel Weihs. Some Hydrodynamical Aspects of Fish Schooling, volume II, pages 703–18. Plenum Press, 1975.