Journal cover Journal topic
Wind Energy Science The interactive open-access journal of the European Academy of Wind Energy
Journal topic
Wind Energ. Sci., 3, 345-352, 2018
https://doi.org/10.5194/wes-3-345-2018
Wind Energ. Sci., 3, 345-352, 2018
https://doi.org/10.5194/wes-3-345-2018

Research articles 07 Jun 2018

Research articles | 07 Jun 2018

# The second curvature correction for the straight segment approximation of periodic vortex wakes

The second curvature correction for periodic vortex wakes
David H. Wood David H. Wood
• Department of Mechanical and Manufacturing Engineering, University of Calgary, Calgary T2N 1N4, AB, Canada
Abstract

The periodic, helical vortex wakes of wind turbines, propellers, and helicopters are often approximated using straight vortex segments which cannot reproduce the binormal velocity associated with the local curvature. This leads to the need for the first curvature correction, which is well known and understood. It is less well known that under some circumstances, the binormal velocity determined from straight segments needs a second correction when the periodicity returns the vortex to the proximity of the point at which the velocity is required. This paper analyzes the second correction by modelling the helical far wake of a wind turbine as an infinite row of equispaced vortex rings of constant radius and circulation. The ring spacing is proportional to the helix pitch. The second correction is required at small vortex pitch, which is typical of the operating conditions of large modern turbines. Then the velocity induced by the periodic wake can greatly exceed the local curvature contribution. The second correction is quadratic in the inverse of the number of segments per ring and linear in the inverse spacing. An approximate expression is developed for the second correction and shown to reduce the errors by an order of magnitude.

1 Introduction

It is common for computational models of the wakes of helicopters, propellers, and wind turbines to use straight vortex segments whose position is iterated until they follow the local flow and the vortex is force free. Solving the Biot–Savart integral gives the induced velocity used in the iteration. Figure 1 shows a representation of a vortex trailing from a two-blade rotor with the straight segment approximation. The labels and symbols on the figure will be defined below. reviewed a range of computational models for wind turbines and describe a recent application of a free-wake vortex model using straight vortex segments.

Figure 1Schematic of two turns of constant radius helical vortices with p=0.1 modelling the far wake of a wind turbine. Nb=2 and the straight segment approximation is shown for Ns=10. The flow is down the page. The “control point” is where the induced velocity is required. The contribution from the “aligned segments” on either side of the control point require the largest second curvature correction.

A well-known difficulty of the straight segment approximation is that it does not reproduce the binormal velocity due to the curvature of the vortex line (e.g. , , and ). This leads to the need for the first curvature correction. To assess the errors of the straight segment approximation and develop a correction, used a vortex ring whose binormal (axial) velocity, U, is given by the well-known Kelvin equation

$\begin{array}{}\text{(1)}& U=\frac{\mathrm{\Gamma }}{\mathrm{4}\mathit{\pi }}\text{log}\left(\frac{\mathrm{8}}{a}\right)-\frac{\mathrm{1}}{\mathrm{4}},\end{array}$

where Γ is the circulation of one vortex, and a is the radius of the vortex core (e.g. ). Figure 2 shows the vortex ring approximation to the helical wake in Fig. 1. Note that by Eq. (1) U increases as a decreases. In this equation, and throughout this paper, all lengths are normalized by the vortex radius (not the core radius a) and all velocities by the wind speed. To reproduce Eq. (1), the numerical evaluation of the Biot–Savart integral for the ring is “cut off” by ignoring the contribution from distances smaller than a from the point at which the velocity is required, the “control point” shown in Figs. 1 and 2. It is emphasized that the cut-off is a heuristic; Kelvin's equation (Eq. 1) is usually derived from impulse considerations or other methods that do not use the Biot–Savart law. documented many factors that alter the vortex velocity from its Biot–Savart value: these include flow along the vortex axis, differing distributions of swirl, etc. Nevertheless, the Biot–Savart prescription is useful and computationally convenient.

Curvature in the wakes of rotors is often associated with vortex periodicity, the “return” of a vortex to the proximity of the control point, which can cause a significant contribution to the binormal velocity. and Wood (2004) used helical line vortices to analyze straight segment errors for this second effect of curvature, but their work has apparently not been considered in subsequent vortex modelling. claimed that the second curvature correction is unnecessary and difficult to implement. The purpose of this paper is to document the importance of the second correction for wind turbine wakes under some operating conditions and to develop an effective and simple correction.

Figure 2Vortex ring representation of the helical wake in Fig. 1 and the corresponding vortex segment approximation.

The paper is organized as follows. The next section introduces the vortex ring model of the wake. In the following section, the induced velocity for the periodic component of the wake over a range of vortex spacings is found in terms of its Biot–Savart integral. Section 4 describes the calculation of the induced velocity for the straight segment approximation, determines the second curvature correction, and tests its accuracy. The final section contains the conclusions.

2 The vortex ring model for the wake

For a point with the same radius as a single vortex ring and distance z from it, the Biot–Savart equation for U in the direction of the wind – the binormal direction – is

$\begin{array}{}\text{(2)}& U=\frac{\mathrm{\Gamma }}{\mathrm{4}\mathit{\pi }}\underset{\mathrm{0}}{\overset{\mathrm{2}\mathit{\pi }}{\int }}\frac{\mathrm{1}-\mathrm{cos}\mathit{\theta }}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{z}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}}\mathrm{d}\mathit{\theta },\end{array}$

where θ is the vortex angle in cylindrical polar co-ordinates. If z=0, the integral clearly has a logarithmic singularity as θ→0. The velocity, U1c, requiring the first curvature correction is

$\begin{array}{}\text{(3)}& {U}_{\mathrm{1}\mathrm{c}}=U\left(z=\mathrm{0}\right),\end{array}$

arising from the only ring containing the control point. The integral in Eq. (2) and similar equations will be termed the “influence coefficient” I, which has the same relative error characteristics as U.

The test case used here to investigate the second correction models the far wake of a wind turbine as an infinite row of equispaced vortex rings of constant spacing, s, radius, and Γ, extending to infinity on either side of the control point at z=0. A row of rings is easier to analyze than the helical vortices used by and Wood (2004) but displays the same need at small separation for the second correction. In addition, the discrete nature of the vortex rings helps to localize the correction that is developed in Sect. 4.

The ring vortex wake is consistent with the “Joukowsky” model of the wake, used by ; either the bound vorticity of the blades is constant along their span or all the shed vorticity has rolled up into tip and hub vortices before reaching the far wake. This is clearly a simplification of wind turbine wakes in general, but the linearity of the Biot–Savart law allows more complex wakes to be considered as an assembly of elements such as rings. The velocity associated with the second correction, U2c, is induced by the vortices that do not contain the control point:

$\begin{array}{}\text{(4)}& {U}_{\mathrm{2}\mathrm{c}}=\frac{\mathrm{\Gamma }}{\mathrm{4}\mathit{\pi }}{I}_{\mathrm{2}\mathrm{c}}=\frac{\mathrm{\Gamma }}{\mathrm{4}\mathit{\pi }}\underset{\mathrm{0}}{\overset{\mathrm{2}\mathit{\pi }}{\int }}\mathrm{2}\sum _{j=\mathrm{1}}^{\mathrm{\infty }}\frac{\mathrm{1}-\mathrm{cos}\mathit{\theta }}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+\left(js{\right)}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}}\mathrm{d}\mathit{\theta }.\end{array}$

Equation (4) is not singular as θ→0, which is, possibly, the reason why the need for the second correction has not been appreciated. s can be identified with the pitch p of a helical vortex wake and the number of blades, Nb, according to $s=\mathrm{2}\mathit{\pi }p/{N}_{\mathrm{b}}$. The relationship between p and s can be seen by comparing Figs. 1 and 2.

Figure 3The integrands for the two influence coefficients for a control point on the vortex ring and θ0=0. The stars show the integrand for I1c from Eq. (2) with z=0. The solid line shows I2c. The dashed and dotted lines are the bounds on I2c, Eq. (5) with s=0.2. The sums were evaluated using Eq. (4) for 50 000 rings.

Testing corrections for the straight segment approximation requires an accurate evaluation of the series in Eq. (4) and then an integration in θ. This order is preferred because the integration in θ of the summand results in incomplete elliptic integrals, which are likely to be very difficult to sum. The innocuous looking series in Eq. (4), however, does not appear to have a closed form sum. The standard technique for summing infinite series of algebraic functions is via Laplace transforms (e.g. ). This would be successful if the exponent in the integrand was 1 instead of 3∕2, but for Eq. (4), the method gave a principal value integral that could not be solved in closed form. By the Cauchy integral test for series

$\begin{array}{ll}\frac{\mathrm{1}}{s}& -\frac{\mathrm{1}}{\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}}}\le {I}_{\mathrm{2}\mathrm{c}}\left(\mathit{\theta }\right)\le \frac{\mathrm{1}}{s}-\frac{\mathrm{1}}{\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}}}\\ \text{(5)}& & +\frac{\mathrm{2}\left(\mathrm{1}-\mathrm{cos}\mathit{\theta }\right)}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}},\end{array}$

where

It is easy to show that the average velocity in the direction of the wind at any radius r<1 within the Joukowsky wake is $\mathrm{1}-\mathrm{\Gamma }/s$ (e.g. ), and it is reasonable to assume that the total velocity of the free-wake vortex rings is close to $\mathrm{1}-\mathrm{\Gamma }/\left(\mathrm{2}s\right)$ or the induced velocity, $U\approx \mathrm{\Gamma }/\left(\mathrm{2}s\right)$. Further, the bounds in Eq. (5), which both contain 1∕s, cause U to approach the average of the wake and external velocities, provided the curvature singularity does not contribute significantly to U. There is, unfortunately, only limited experimental information on a and U for wind turbine wakes to guide the assessment of the relative importance of the first and second velocity fields and their corrections. Figure 3 shows the terms in Eq. (5) for s=0.2. This typical value was obtained using the following steps. For modern turbines, Nb=3, and λ≈7 for most of the operating range. For optimal (Betz-Joukowsky) performance, $U=\mathrm{1}/\mathrm{3}$, $p\approx \mathrm{2}/\left(\mathrm{3}\mathit{\lambda }\right)$, where λ is the tip speed ratio, and ${N}_{\mathrm{b}}\mathrm{\Gamma }\mathit{\lambda }/\mathit{\pi }=\mathrm{8}/\mathrm{9}$ (Wood2011). Thus Γ=0.133 and s≈0.2. The sum in Eq. (4) is always zero when θ=0, but, as s decreases, the bounds in Eq. (5) (and hence the sum) tend to 2πs over an increasing range of θ. Integrating over [0,2π] then leads to $U\approx {U}_{\mathrm{2}\mathrm{c}}\approx \mathrm{\Gamma }/\left(\mathrm{2}s\right)$, showing the potential importance of U2c. The integrand of I1c for a small θ is also shown. Its integral and U1c depend on the cut-off, a; to match $U=\mathrm{2}/\mathrm{3}$ for the conditions in Fig. 3, Eq. (1) requires $a\sim {\mathrm{10}}^{-\mathrm{23}}$, which does not seem a reasonable value. Thus it is likely that UU2c at the small p and s typical of the operating conditions of modern wind turbines.

The very limited experimental information on the velocity of the vortices in wind turbine wakes are in general agreement with this argument. measured the wake of a two-bladed turbine in a wind tunnel at λ=4.91 using particle image velocimetry. They determined the vortex velocity in the near wake as 10.8 m s−1 when the wind speed was 12 m s−1. Thus $U=\left(\mathrm{12}-\mathrm{10.8}\right)/\mathrm{12}=\mathrm{0.1}$, which is lower than the value of 1∕3 that follows from assuming optimal power output. Assuming U=U2c and using the general equation $p=\left(\mathrm{1}-U\right)/\mathit{\lambda }$ gives s=0.576 or 360 mm for the rotor of radius 625 mm, which agrees very well with the value read from their Fig. 10. This again implies that U2cU1c. Since most rotor wakes are helices of some form, it is important to note that the equivalent inverse pitch term dominates U for a helical vortex of sufficiently small p . There is a further reason to expect $U\approx {U}_{\mathrm{2}\mathrm{c}}\gg {U}_{\mathrm{1}\mathrm{c}}$ for many turbine wakes: U1c, but not U2c, is associated with the impulse necessary to form a vortex ring. If that impulse and U1c are significant it is unlikely that the wake can be force free.

3 Evaluating the influence coefficient for an infinite array of vortex rings

A closed form sum for U2c in Eq. (4) could not be obtained so the influence coefficients for a range of s values were determined as follows. The Hermite–Hadamard inequality for monotonically decreasing functions that tend to zero at large argument can be used simply to give a tighter bound on I2c(θ). It is

$\begin{array}{}\text{(7)}& {I}_{\mathrm{2}\mathrm{c}}\left(\mathit{\theta }\right)=\frac{\mathrm{1}}{s}-\frac{\mathrm{1}}{\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}}}+\frac{\mathrm{1}-\mathrm{cos}\mathit{\theta }}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}}+\mathit{\delta }\left(\mathit{\theta }\right),\end{array}$

where the difference, δ(θ), is always positive but must be determined numerically. The integral of the other terms on the right side of Eq. (7) can be found exactly:

$\begin{array}{ll}& \underset{\mathrm{0}}{\overset{\mathrm{2}\mathit{\pi }}{\int }}\left(\frac{\mathrm{1}}{s}-\frac{\mathrm{1}}{\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}}}+\frac{\mathrm{1}-\mathrm{cos}\mathit{\theta }}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}}\right)\mathrm{d}\mathit{\theta }\\ \text{(8)}& & =\frac{\mathrm{2}\mathit{\pi }}{s}-\frac{\mathrm{2}sE\left(-\mathrm{4}/{s}^{\mathrm{2}}\right)}{{s}^{\mathrm{2}}+\mathrm{4}}-\frac{\mathrm{2}K\left(-\mathrm{4}/{s}^{\mathrm{2}}\right)}{s},\end{array}$

where E(.) and K(.) are the complete elliptic integrals in standard notation. The difference, δ, the integral of δ(θ) over [0,2π], was evaluated using 2000 increments of θ and number of rings, Nr=50 000. This value was chosen using the result obtained from Mathematica, that

$\begin{array}{ll}& \underset{\mathrm{0}}{\overset{\mathrm{2}\mathit{\pi }}{\int }}\mathrm{2}\sum _{j={N}_{\mathrm{r}}}^{\mathrm{\infty }}\frac{\mathrm{1}-\mathrm{cos}\mathit{\theta }}{\left(\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+\left(js{\right)}^{\mathrm{2}}{\right)}^{\mathrm{3}/\mathrm{2}}}\mathrm{d}\mathit{\theta }\\ \text{(9)}& & \to R\left({N}_{\mathrm{r}}\right)=\frac{\mathrm{4}\mathit{\pi }\left(\mathit{\zeta }\left(\mathrm{3}\right)-{H}_{{N}_{\mathrm{r}}}^{\left(\mathrm{3}\right)}\right)}{{s}^{\mathrm{3}}}\end{array}$

for large Nr, where ζ(.) is the zeta function, and H is the harmonic number in standard notation. In later use of this result, R(j) will be called the “remainder”. Using 2πs as an estimate for the integral in Eq. (4), and using Mathematica to evaluate H3(Nr), gave the relative error in truncating the sum at Nr=50 000 as $\mathrm{2}/{s}^{\mathrm{2}}×{\mathrm{10}}^{-\mathrm{10}}=\mathrm{2}×{\mathrm{10}}^{-\mathrm{8}}$ for the smallest value of s considered here, s=0.1. To the number of decimal places used in Table 1, truncation does not alter the integral of the terms in Eq. (7) over [0,2π]. For every calculation up to Sect. 5, Nr=50 000. The θintegral of the sum in Eq. (7) was found using the MATLAB quadrature routine “integral” with an absolute tolerance of 10−8. All integrands are symmetric about θ=π and so were obtained over [0,π]. Table 1 also shows the approximate λ for a Betz–Joukowsky optimal rotor. δ is small: combining Eqs. (6) and (7) leads to

$\begin{array}{}\text{(10)}& \mathrm{0}\le \mathit{\delta }\le \frac{\mathrm{2}K\left(-\mathrm{4}/{s}^{\mathrm{2}}\right)}{s}-\frac{\mathrm{2}sE\left(-\mathrm{4}/{s}^{\mathrm{2}}\right)}{{s}^{\mathrm{2}}+\mathrm{4}},\end{array}$

which is satisfied by all δ values in Table 1. δ is less than 4 % of Ic for the worst case of s=0.8; Eq. (8) is an increasingly good approximation to the sum and the influence coefficient approaches 2πs as s decreases and λ increases. The values of I2c in Table 1 will be compared to the values from the straight segment approximation.

Table 1Values of the influence coefficient for varying s values with θ0=0.

4 Straight segment approximation of vortex rings and its accuracy

Each of the rings not containing the control point was approximated by an even number, Ns, of straight segments. The ith segment of each ring started at $\mathit{\theta }=\mathrm{2}\mathit{\pi }i/{N}_{\mathrm{s}}+{\mathit{\theta }}_{\mathrm{0}}$ when measured from the control point and finished at $\mathrm{2}\mathit{\pi }\left(i+\mathrm{1}\right)/{N}_{\mathrm{s}}+{\mathit{\theta }}_{\mathrm{0}}$. θ0 is the angular displacement between the control point and the start of the first (i=0) segment. A straightforward application of Eq. (10.115) of gives I2c(i,j), the contribution of the ith segment on the two jth vortex rings to U2c, as

and

$\begin{array}{ll}\text{(12)}& {A}_{i}& =\mathrm{8}\mathrm{sin}\left(\frac{\mathit{\pi }i}{{N}_{\mathrm{s}}}+\frac{{\mathit{\theta }}_{\mathrm{0}}}{\mathrm{2}}\right)\mathrm{sin}\left(\frac{\mathit{\pi }\left(i+\mathrm{1}\right)}{{N}_{\mathrm{s}}}+\frac{{\mathit{\theta }}_{\mathrm{0}}}{\mathrm{2}}\right),{B}_{i,j}& =\mathrm{sin}\frac{\mathit{\pi }}{{N}_{\mathrm{s}}}\left(\frac{\mathrm{1}}{{b}_{i+\mathrm{1},j}}+\frac{\mathrm{1}}{{b}_{i,j}}\right)\\ \text{(13)}& & +\mathrm{sin}\left(\frac{\left(\mathrm{2}i+\mathrm{1}\right)\mathit{\pi }}{{N}_{\mathrm{s}}}+{\mathit{\theta }}_{\mathrm{0}}\right)\left(\frac{\mathrm{1}}{{b}_{i+\mathrm{1},j}}-\frac{\mathrm{1}}{{b}_{i,j}}\right),{C}_{i,j}& =\mathrm{cos}\frac{\mathrm{2}\mathit{\pi }}{{N}_{\mathrm{s}}}+{b}_{i,j}^{\mathrm{2}}+{b}_{i+\mathrm{1},j}^{\mathrm{2}}\\ \text{(14)}& & +\mathrm{cos}\left(\frac{\mathrm{2}\mathit{\pi }\left(\mathrm{2}i+\mathrm{1}\right)}{{N}_{\mathrm{s}}}+\mathrm{2}{\mathit{\theta }}_{\mathrm{0}}\right)-\mathrm{2},\end{array}$

and

$\begin{array}{}\text{(15)}& {b}_{i,j}=\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\left(\frac{\mathrm{2}\mathit{\pi }i}{{N}_{\mathrm{s}}}+{\mathit{\theta }}_{\mathrm{0}}\right)+\left(js{\right)}^{\mathrm{2}}}.\end{array}$

For the first calculations, the junction of the first and Nsth segment was aligned with the control point so θ0=0. The influence coefficients calculated from Eqs. (11)–(14) are compared to the results from Table 1 in Fig. 4 in terms of the relative error using the integral over [0,π] as the denominator since the integrand must be symmetric about θ=π. The error is defined as the exact integral minus the straight segment approximation so the latter is always an underestimate. Two separate ranges of θ are considered: for $\mathrm{0}\le |\mathit{\theta }|\le \mathit{\pi }/\mathrm{2}$ the error is 1 to 2 orders of magnitude higher than for $\mathit{\pi }/\mathrm{2}\le |\mathit{\theta }|\le \mathit{\pi }$ and the errors in the first range increase with decreasing s whereas in the second range they decrease. Figure 5 shows the reason. For this test case, no aligned segment contributes to U2c as the Biot–Savart velocity must lie in the plane containing the segment and the control point. Otherwise, both errors scale as $\mathrm{1}/{N}_{\mathrm{s}}^{\mathrm{2}}$, as was found in the helix simulations of Wood (2004). As s decreases, however, the aligned segment error increases proportionally to 1∕s, but the error for the remaining range of θ decreases at constant Ns.

Figure 4Relative error of the straight segment approximation for the conditions in Table 1 (θ0=0). +: s=0.1; : s=0.2; ×: s=0.4; : s=0.8. The top group of symbols show the errors for the first half of the rings that contain the aligned segments, $\mathrm{0}\le |\mathit{\theta }|\le \mathit{\pi }/\mathrm{2}$, and the bottom group covers $\mathit{\pi }/\mathrm{2}\le |\mathit{\theta }|\le \mathit{\pi }$. ${N}_{\mathrm{s}}=\mathrm{20},\mathrm{40},\mathrm{80}$, and 160.

Figure 5Angular contribution of straight segments to the influence coefficient for z=0.2, Ns=40, and θ0=0. The value of θ is the midpoint of the segment. The solid line shows the numerical solution of the exact integral. shows the straight segment contribution, × the contribution from Eq. (17), and + the contribution from Eq. (18).

Figure 5 shows the angular contribution to the influence coefficients for s=0.2 and Ns=40. The value of θ used for plotting is the midpoint of each segment. The solid line shows the exact integral from θi to θi+1. As was found by and Wood (2004), the errors are localized near θ=0. A correction for the error for the aligned segments can be developed from the small-θ expansion of the series in Eq. (4):

and ζ(3)=1.2026. Equation (16) has two important implications. First, the best possible error for periodic straight segments scales as (Nss)−3 but it is likely that an unrealistically high value of Ns would be required to achieve this. Second, 1∕ζ(3) or over 80 % of the correction to U2c is due to the two rings (j=1) on either side of the control point. This is the justification for pointing out the aligned segments in Figs. 1 and 2. A general form of the correction, therefore, can be based on the returned vortex on either side of the control point. Since the distance from the control point to the vortex segments must be calculated in a free-wake simulation, it should not be difficult to determine the proximity in terms of θz and apply a correction. A more general correction is Δ(θs), where ${\mathit{\theta }}_{\mathrm{s}}=\mathrm{2}\mathit{\pi }/{N}_{\mathrm{s}}$ is obtained by integrating in θ only for j=1 and then using ζ(3) to correct approximately for the remaining rings. The result is

$\begin{array}{ll}\mathrm{\Delta }\left({\mathit{\theta }}_{\mathrm{s}}\right)& \approx \mathrm{2}\mathit{\zeta }\left(\mathrm{3}\right)\left[\frac{F\left({\mathit{\theta }}_{\mathrm{s}}/\mathrm{2},-\mathrm{4}/{s}^{\mathrm{2}}\right)}{s}-\frac{sE\left({\mathit{\theta }}_{\mathrm{s}}/\mathrm{2},-\mathrm{4}/{s}^{\mathrm{2}}\right)}{{s}^{\mathrm{2}}+\mathrm{4}}\\ \text{(17)}& & -\frac{\mathrm{2}\mathrm{sin}{\mathit{\theta }}_{\mathrm{s}}}{\left({s}^{\mathrm{2}}+\mathrm{4}\right)\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}{\mathit{\theta }}_{\mathrm{s}}+{s}^{\mathrm{2}}}}\right],\end{array}$

where E(.) and F(.) are the incomplete elliptic integrals. Equation (17) is shown in Fig. 3 to give a better estimate for the aligned segments. An alternative, simpler correction than Eq. (17) can be found by using $\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }\sim {\mathit{\theta }}^{\mathrm{2}}$ for small θ values to give

$\begin{array}{}\text{(18)}& \mathrm{\Delta }\left({\mathit{\theta }}_{\mathrm{s}}\right)\approx \mathrm{2}\mathit{\zeta }\left(\mathrm{3}\right)\left[\mathrm{log}\left(\frac{{\mathit{\theta }}_{\mathrm{s}}+\sqrt{{\mathit{\theta }}_{\mathrm{s}}^{\mathrm{2}}+{s}^{\mathrm{2}}}}{s}\right)-\frac{{\mathit{\theta }}_{\mathrm{s}}}{\sqrt{{\mathit{\theta }}_{\mathrm{s}}^{\mathrm{2}}+{s}^{\mathrm{2}}}}\right],\end{array}$

which gives almost the same correction for the aligned segments; Fig. 5. These results for the application of Eqs. (17) and (18) to the aligned segments are similar at the other values of s as well, but are not shown in the interests of brevity. It is noted that the correction developed here is simple in the sense that the vortex curvature is known a priori. As pointed out by , however, and shown by the analysis of , the modelling of three-dimensional wakes of varying geometry can be considerably more complex.

Figure 6Variation in relative error with θ0 for s=0.1. Ns=20, +; Ns=40, ; Ns=80, ×; Ns=160, .

One of these complexities is that the control point may not align with the junction of segments on (in this case) adjacent rings. The effect of this can be investigated by using non-zero θ0 in Eqs. (12)–(15). The results are shown in Fig. 6 for $\mathrm{20}\le {N}_{\mathrm{s}}\le \mathrm{160}$ and $\mathrm{0}\le {\mathit{\theta }}_{\mathrm{0}}\le \mathit{\pi }/{N}_{\mathrm{s}}$ and s=0.1. As was found for other values of s, there is remarkably little variation in the error with θ0 except for the lowest Ns, suggesting that the correction derived above for the aligned case (θ0=0) is also applicable to other values. This is not an immediately obvious result from Eqs. (12)–(15). For ${\mathit{\theta }}_{\mathrm{0}}=\mathit{\pi }/{N}_{\mathrm{s}}$ and $\mathit{\theta }=\mathrm{2}\mathit{\pi }/{N}_{\mathrm{s}}$:

$\begin{array}{ll}{I}_{\mathrm{2}\mathrm{c}}\left(\mathit{\theta },\mathrm{1}\right)& =\frac{\mathrm{8}{\mathrm{sin}}^{\mathrm{2}}\left(\mathit{\theta }/\mathrm{2}\right)\mathrm{sin}\mathit{\theta }}{\phantom{\rule{0.33em}{0ex}}\left(\mathrm{3}-\mathrm{4}\mathrm{cos}\mathit{\theta }+\mathrm{cos}\mathrm{2}\mathit{\theta }+\mathrm{2}{s}^{\mathrm{2}}\right)\sqrt{\mathrm{2}-\mathrm{2}\mathrm{cos}\mathit{\theta }+{s}^{\mathrm{2}}}}\\ \text{(19)}& & \to \left(\frac{\mathit{\theta }}{s}{\right)}^{\mathrm{3}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\text{as}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathit{\theta }/s\to \mathrm{0},\end{array}$

which suggests a difference from the case when θ0=0.

Figure 7Variation in the straight segment approximation to I2c(θ) for θ0=0 and Nr= 5, ×; 10, ; 20, ; +, 50; ×, 50 000. Ns=40, s=0.2.

5 Using a finite array of rings to determine the influence coefficient

The second curvature error was shown in the last section to be caused largely by aligned segments on the rings either side of the control point. For increasing θ, I2c(θ) becomes dominated by rings at a larger distance from the control point. This is shown in Fig. 7, which implies that Nr either must be large to ensure an accurate determination of I2c or a suitable remainder term be used. This allows an approximate determination of the influence coefficient for the case in which θ0=0 according to

$\begin{array}{}\text{(20)}& {I}_{\mathrm{2}\mathrm{c}}\approx \mathrm{2}\mathrm{\Delta }\left({\mathit{\theta }}_{\mathrm{s}}\right)+\sum _{j=\mathrm{1}}^{{N}_{\mathrm{r}}}\sum _{i=\mathrm{1}}^{{N}_{\mathrm{s}}}{I}_{\mathrm{2}\mathrm{c}}\left(i,j\right)+R\left({N}_{\mathrm{r}}\right),\end{array}$

where one possibility for the remainder R(Nr) is given by Eq. (9).

Table 2Terms in Eq. (20) for s=0.2 and θ0=0. Exact value of I2c=26.889628. The error is the relative error.

The terms in Eq. (20) are listed in Table 2 for s=0.2. A significant number of vortex rings, Nr, or equivalently a large stream-wise distance is needed to make the remainder, R(Nr), accurate. Typically, Nr≥20 for this s, and then R(Nr) is comparable to Δ(θs). For Ns=20, for example, after applying an accurate remainder, the second curvature correction changes the relative error from 3.4 % to less than 0.2 %, which is a reduction by 2 orders of magnitude.

6 Conclusions

The widely used straight segment approximation for approximating the curved and periodic vortex wakes of wind turbines, propellers, and helicopters can have two errors associated with the wake curvature. The first is the well-known error in reproducing the locally induced binormal velocity. This is usually accommodated by a cut-off in the Biot–Savart determination of the vortex velocity using Eqs. (2) and (3) at a distance comparable with the radius of the vortex core. The second, less well-known error is the subject of this paper. It arises from the alignment of the segments of the periodic vortex returning to the proximity of the point at which the velocity is being determined.

By modelling the far wake of a wind turbine as an infinite row of equispaced vortex rings, two important results were obtained. First, it was shown that the velocity associated with the second error dominates at the small spacings typical of modern wind turbine operation. The available experimental evidence on wake structure is consistent with this finding. Then it is shown that the second error is quadratic in the number of segments per revolution and inversely proportional to the spacing of the rings, which is proportional to the pitch of a more realistic, but more difficult, helical wake. The model to investigate the second correction is artificial in that a single, infinite row of vortex rings of constant spacing, radius, and circulation is not applicable to the near wake. Nevertheless the model demonstrated the general importance of the rings adjacent to the control point at which the velocity is being calculated. These adjacent rings contribute over 80 % of the correction that is needed because the straight segment approximation does not correctly determine the contribution to the induced velocity from the closest parts of the adjacent rings, called the aligned segments.

It was also shown that the best behaviour possible for the second error is cubic in the product of the number of segments per revolution and the vortex spacing. It is likely, however, that larger numbers of vortex segments would be needed to achieve this error than are used in practice. This result was generalized to develop a second correction that improves the computed induced velocity by nearly 1 order of magnitude.

Code availability
Code availability.

The MATLAB codes used in this study are available from the author.

Competing interests
Competing interests.

The author declares that he has no conflict of interest.

Acknowledgements
Acknowledgements.

This work is part of a research project on wind turbine aerodynamics funded by the NSERC Discovery Grants Program.

Edited by: Alessandro Bianchini
Reviewed by: Joseph Saverin and Wang Xiaodong

References

Bhagwat, M. J. and Leishman, J. G.: Self-Induced Velocity of a Vortex Ring Using Straight-Line Segmentation, AIAA Journal, 59, 1–7, 2014. a, b

Govindarajan, B. M. and Leishman, J. G.: Curvature Corrections to Improve the Accuracy of Free-Vortex Methods, J. Aircraft, 53, 378–386, 2016. a, b, c

Katz, J. and Plotkin, A.: Low-Speed Aerodynamics, 2nd edition, C.U.P., Cambridge, 2001. a

Kim, C. J., Park, S. H., Sung, S. K., and Jung, S. N.: Dynamic modeling and analysis of vortex filament motion using a novel curve-fitting method, Chinese J. Aeronaut., 29, 53–65, 2016. a, b

Kuibin, P. A. and Okulov, V. L.: Self-induced motion and asymptotic expansion of the velocity field in the vicinity of a helical vortex filament, Phys. Fluids, 10, 607–614, 1998. a

O'Brien, J. M., Young, T. M. O'Mahoney, D. C., and Griffin, P. C.: Horizontal axis wind turbine research: A review of commercial CFD, FE codes and experimental practices, Prog. Aerosp. Sci., 97, 1–24, 2017. a

Saffman, P. G.: Vortex Dynamics, C.U.P., Cambridge, 1992. a, b

Sarmast, S., Segalini, A., Mikkelsen, R.  F., and Ivanell, S.: Comparison of the near-wake between actuator-line simulations and a simplified vortex model of a horizontal-axis wind turbine, Wind Energy, 19, 471–481, 2016. a, b

Wheelon, A. D.: On the Summation of Infinite Series in Closed Form, J. Appl. Phys., 113, 113–118, 1954. a

Wood, D. H.: Method to Improve the Accuracy of Straight-Segment Representation of Helical Vortices, AIAA Journal, 41, 256–262, 2004.  a, b, c, d

Wood, D. H.: Small Wind Turbines: Analysis, Design, and Application, Springer, London, 2011. a, b

Wood, D. H. and Li, D.: Assessment of the accuracy of representing a helical vortex by straight segments, AIAA Journal, 40, 647–651, 2002. a, b, c

Xiao, J. P., Wu, J., Chen, L., and Shi, Z. Y.: Particle image velocimetry (PIV) measurements of tip vortex wake structure of wind turbine, Appl. Math. Mech.-Engl., 32, 729–738, 2011. a