Using the concept of impulse in control volume analysis, we derive general expressions for wind turbine thrust in a constant, spatially uniform wind. The absence of pressure in the impulse equations allows for their application in the near wake, where the flow is assumed to be steady in the frame of reference rotating with the blades. The assumption of circumferential uniformity in the near wake – as applies when the number of blades or the tip speed ratio tends to infinity – is needed to reduce these general expressions to the Kutta–Joukowsky (KJ) equation for blade-element thrust. The present derivation improves upon the classical derivation based on the Bernoulli equation by allowing the flow to be rotational in the near wake. The present derivation also yields intermediate expressions for thrust that are valid for a finite number of blades and trailing vortex sheets of finite thickness. For the circumferentially uniform case, our analysis suggests that the magnitudes of the radial velocity and the axial induction factor must be equal somewhere on the plane containing the rotor, and we cite previous studies that show this to occur near the rotor tip across a wide range of thrust coefficients. The derivation reveals one further complication; when deriving the KJ equations using annular control volumes, the existence of vorticity on the lateral control surfaces may cause the local blade loading to differ from the KJ equation, but the magnitude of these deviations is not explored. This complication is not visible to the classical derivation due to its neglect of vorticity.

Introduction

Blade-element theory (BET) for wind turbines uses the fundamental assumption that the forces acting on the elements comprising the rotor blades are given by the Kutta–Joukowsky (KJ) theorem. The thrust and torque are balanced by the change in the axial and angular momentum, respectively, of the flow through a control volume (CV) enclosing the rotor; the combination of BET and momentum theory gives rise to blade-element momentum (BEM) theory. BEM is developed in many texts on wind turbine aerodynamics, such as , , , , and , so it is unnecessary to repeat it here. Herein, the KJ equation for blade-element thrust will be rederived using impulse theory, and the conditions of its validity are examined.

Although the KJ equations are generally introduced as assumptions in BEM theory, they can be derived using the unsteady Bernoulli equation. This was shown as early as and can also be found in and . Assuming the trailing wake to be infinitely thin to permit the use of the Bernoulli equation, and representing the rotor's power extraction by a sudden drop in pressure across the rotor plane, sectional thrust can be expressed as follows: 1ρdTdx=02π12uθ2-λuθxxdθ, where ρ is fluid density, T is the rotor thrust, x is the radial coordinate, λ is the tip speed ratio, and θ is the azimuthal angle. uθ is the circumferential velocity just behind the rotor. For ease of reference, we have included our own version of the classical Bernoulli-based derivation of Eq. () in Appendix A.

In keeping with typical notation in BEM analysis, we define uθ=-2w on this near-wake plane, such that positive w corresponds to positive torque. This substitution gives 1ρdTdx=02π2w2+2λwxxdθ. w is generally interpreted as the circumferential induced velocity at the rotor plane itself, which is assumed in BEM theory to be half that at the plane just behind the rotor. All quantities in Eqs. () and () have been normalized by freestream velocity and rotor radius.

When the flow is circumferentially uniform, Eq. () reduces to the KJ lift theorem, as will be shown later using Stokes' theorem. Here, we merely point out the contradiction in this classical approach. Although the derivation of Eq. () required us to assume irrotational flow everywhere in the near wake – save in the vortical wake of vanishing thickness – the assumption of circumferential uniformity will lead to non-zero vorticity on the near-wake plane except in unrealistic special cases (i.e. uθ0 or uθ=k/x, where k is a constant). Using impulse theory, we offer an improved derivation of the KJ equation for blade-element thrust that avoids this contradiction. The derivation also offers new insights regarding the relative magnitude of radial velocity in the expanding flow upwind and through the rotor.

The remainder of the paper is organized as follows. Since the wind energy research community is likely to be unfamiliar with the concept of “vortical impulse” or “hydrodynamic impulse” – usually referred to simply as “impulse” – we offer a brief introduction in Sect. . Section  then manipulates the general CV equations furnished by impulse theory to yield simplified integral expressions of wind turbine thrust. The reader interested primarily in the consequences of the present impulse analysis for BEM may skip ahead to Sect. . Finally, our conclusions are presented in Sect. .

Background on impulse theory

Impulse theory expresses fluid-dynamic forces in terms of the first moment of vorticity (e.g. ): 1N-1Vx×ΩdV, where V is some fluid volume, N is the dimension of the space, x is the position vector, and Ω is the vorticity vector. The use of the word “impulse” may confuse the unfamiliar reader; it is an imprecise yet well-established nomenclature. showed that the first moment of vorticity is equal to the impulse (i.e. the integral of force applied over time) necessary to establish an unbounded vortical flow from rest (the domain is unbounded, but the region of vortical fluid remains bounded). Since then, the term “impulse” has been co-opted to refer specifically to the first moment of vorticity, as recounted in Sect. 2.5 of , and we will continue to use this terminology here.

The introduction of the concept of impulse removes the pressure and introduces vorticity to the equations of linear momentum conservation. This approach can be traced at least as far back as , who sought to determine the speed of a vortex ring. In their review, recount that exact impulse-based expressions for aerodynamic force were derived independently by , , and . Discussions of impulse formulations can also be found in , , and . They have recently gained popularity for use with planar or volume measurements of fluid velocities from particle image velocimetry and related techniques (e.g. ).

Derivation of an impulse-based force formulation begins with the conventional application of the Reynolds transport theorem to the momentum equation and proceeds by manipulating the equations using vector calculus identities. For a full presentation, we recommend the doctoral thesis, as only a brief overview is given here.

For an inertial CV in an incompressible fluid, the equation for force, F, is given in many fluid mechanics texts and by , Eq. (3.36): Fρ=-ddtVUdV+Sn-PρI-UU+TdS, where t is time, and U is the velocity vector. S is the control surface bounding the control volume, V; n is the outward-facing normal on S; I is the unit tensor; P is the pressure; and T is the viscous stress tensor. The nUU term gives the conventional momentum deficit when the equation is used to determine thrust.

The momentum integral – the first integral on the right-hand side of Eq. () – can be related to the impulse as given in Eq. () by the so-called “impulse-momentum identity” (Noca's Eq. 3.1) at the expense of additional terms. By then removing P using the Navier-Stokes equations and vector identities, Noca obtained his Eqs. (3.55) and (3.56), which we combine as Fρ=-1N-1ddtVx×ΩdV+Sn(12U2I-UU-1N-1U(x×Ω)5+1N-1Ω(x×U))dS+Fv, where N=3 for our analysis, 12U2=12UU is the kinetic energy per unit mass, and Fv is shorthand for viscous terms evaluated on S. In common with other CV analyses of wind turbines, Fv will be ignored. Equation (), or its variations in , , and , is the complete, general form of the impulse formulation for fluid-dynamic force for an impermeable body in steady motion. It is as exact as the more familiar Eq. (), and it serves as the starting point for our analysis.

The impulse equation for rotor thrust

We consider a wind turbine rotating steadily at a tip speed ratio, λ, defined as the ratio of the circumferential velocity of the blade tips to the wind speed. The latter velocity will be used to normalize all velocities. All lengths are normalized by rotor radius, R. We start by using two variants of a cylindrical CV, both of whose upwind face is well ahead of the blades, and whose radius RCVR. This ensures that the streamwise velocity along the surface at RCV is equal to the wind speed. The downstream control surface is either just downstream or just upstream of the rotor plane; these variants are called CV1 and CV2, respectively, as shown in Fig. . In both cases, the origin of the cylindrical coordinate system (x,θ,z) is centred on the downstream face of the CV, as shown on the left of Fig. . x is the radial coordinate, z is the axial coordinate, and θ is the azimuthal angular coordinate (and the subscript θ refers to the circumferential direction).

Control volumes (CVs) for the present analysis. In both variants, the upstream face is well upstream of the rotor, where the velocity is equal to the freestream, and RCVR. The downstream control surface is just downstream and just upstream of the rotor plane in CV1 and CV2, respectively, and the corresponding downstream control surfaces (CSs) are labelled SD and SU.

We now consider the contributions of each term in Eq. () to the thrust, defined as T=Fez, where ez is the unit vector in the z direction. The time derivative of the impulse integral – the first integral in Eq. () – does not necessarily vanish identically, but its thrust contribution does. Since the vortical wake is assumed to rotate with the blades, only the direction (but not the magnitude) of the impulse integral can change, and its time derivative reduces to ddtCV1x×ΩdV=Λ×CV1x×ΩdV, where Λ is the rotation vector. For steady rotation, Λ=λez, and the thrust contribution of Eq. () vanishes since Λ and ez are parallel.

The first two terms in the surface integral on the right-hand side of Eq. () can be referred to jointly as the velocity terms. The axial component of the nUU term gives the conventional momentum-deficit term, derived in any textbook on wind turbine aerodynamics: -ezCS1nUUdS=SDu(1-u)dS=SDa(1-a)dS where CS1 is the external control surface bounding CV1, and SD denotes the downwind face of CV1. u=1-a is the axial velocity magnitude, and a is the well-known axial induction factor. The n(12U2I) term, which we will call the kinetic energy term, contains contributions from all three components of velocity; the thrust contribution from axial velocity is -12SD(1-u2)dS=-12SD(2a-a2)dS. The right-hand sides of Eqs. () and () can be summed to give -12SDa2dS. The kinetic energy term also has a contribution from the radial velocity, v: 12SDv2dS, and the contribution from the circumferential velocity is 2SDw2dS, where we have continued the convention to express the circumferential velocity on SD as uθ=-2w. The velocities on the lateral cylindrical control surface do not contribute to the kinetic energy term because its normal vector is perpendicular everywhere to ez. Vorticity can be non-zero only on the downstream face of the CV, and the last two terms on the right-hand side of Eq. (), which we call the vortex terms, simplify to 12SD(1-a)xΩθdS-SDΩzwxdS. Because the radial components of velocity and vorticity are parallel to x=xer on SD, where er is the radial unit vector, they lend no thrust contribution through the vortex terms.

So, we have seen that the impulse term's thrust contribution vanishes, the two velocity terms in Eq. () are replaced by the sum of Eqs. () through (), and the vortex terms lend the contribution given in Eq. (). Summing these contributions, the thrust becomes Tρ=12[SDv2dS-SDa2dS+4SDw2dS14+SD(1-a)xΩθdS-SD2ΩzwxdS], where all velocities and vorticities are expressed relative to an inertial frame.

Velocity (in the rotating frame of reference) and vorticity vectors in the near wake. If diffusion is neglected and the wake is assumed to be rigid, velocity and vorticity triangles on each orthogonal plane are geometrically similar.

Equation () has incorporated the assumptions of steady rotation, a rigid wake, and vanishingly small viscosity. To proceed towards the KJ equation for blade-element thrust, let us assume that the effect of diffusion is negligible on SD, such that Kelvin's circulation theorem applies to individual fluid elements. When combined with the assumption of a rigid wake, it is clear that vortex lines must coincide with streamlines in the rotating frame of reference, or else the structure of the vortical wake would appear to deform in that frame. This alignment of streamlines and vortex lines in three dimensions is illustrated in Fig. . On any two-dimensional plane cutting through our three-dimensional volume, the projected velocity and vorticity components generate similar triangles. We are specifically interested in the resulting relationship between axial and circumferential vector components: 2w+λx1-a=ΩθΩz. This becomes useful if we rewrite the vortex terms as 12SD(1-a)Ωθ-(2w+λx)ΩzxdS+12SDΩzλx2dS, because the terms in the square brackets of the first integrand cancel by Eq. (). The remaining term in Eq. () can be interpreted as a thrust due to wake rotation. We now manipulate this term to express it in terms of velocity. Substituting dS=xdxdθ into the previous expression, the term becomes 12λ02π0Ωzx3dxdθ. The streamwise vorticity in polar coordinates is Ωz=-2dwdx+wx-1xvθ, where the negative sign on the w terms is due to our chosen sign convection, whereby positive w corresponds to positive torque. When substituted into Eq. (), the contribution of v/θ vanishes: -12λ0x202πvθdθdx=0. The vortex terms thus become -λ02π01dwdxx3dx+01wx2dxdθ. After integrating the first term by parts, and assuming wx30 as x, the whole expression becomes 2λ02π0wx2dxdθ=2λSDwxdS. Substituting the expression on the right-hand side into Eq. (), we arrive at Tρ=12SDv2-a2dS+2SDw2+λwxdS. It is now clear that our derivation is approaching Eq. (), the precursor to the KJ thrust equation. The next step is to investigate the relationship between v and a, which is done in the next section.

Since a future goal of this research program is unsteady turbine modelling, an alternative derivation of Eq. () in the rotating frame of reference is offered in Appendix B.

If the downwind face of the CV is moved to be just ahead of the rotor – replacing integrals on SD with integrals on SU in our first use of CV2 shown in Fig.  – then T in Eq. () becomes zero, since no body is enclosed by the CV. The flow on SU is irrotational, and so the vortex terms in Eq. () also vanish, and our CV impulse analysis reduces to 0=SUv2-a2+uθ2dS. If there are azimuthal variations of v on SU, then there must be corresponding radial variations in uθ to ensure that the axial vorticity remains zero. In the special case of azimuthal uniformity, uθ must vanish on SU so that the circulation around any circle of radius x remains zero . Thus, 0=SUv2-a2dS for the circumferentially uniform case. directly corroborates Eq. (), showing it to hold for his model of rotor flow. Since v and a must be continuous across any actuator disk, the equation also holds immediately behind the rotor, but that is the limit to its validity; it cannot hold, for example, in the far wake, where v is everywhere zero but a0.

Equation () implies that v2=a2 somewhere on SU if both v and a are C0-continuous. In a real flow with finite viscosity, we indeed expect both variables to vary continuously, even if we allow for arbitrarily large (but finite) gradients. At the rotor axis, v=0 by symmetry, but in general a>0 at the axis for a loaded turbine rotor. If a2>v2 over the inner rotor, Eq. () requires that v2>a2 over some other part of SU, and the restriction that v and a are continuous then requires v2=a2 for at least one radial location. Previous actuator-disk simulations show that v>a in the flow external to the rotor (x>1), and that v=a occurs near the tip of the actuator disk for a wide range of thrust coefficients, e.g. Figs. 15 and 16 of , Figs. 3.2 and 3.3 of , and Fig. 4 of . These results indicate a significant radial deflection of the streamlines in the tip region as they pass through a loaded rotor. We are unaware of any study of the effect of this “crossflow” on blade-element forces, but simple expressions for crossflow alterations to airfoil lift and drag are developed by .

For a turbine with a finite number of blades, the situation is more complex. The early near-wake measurements of and the recent ones of , when converted to co-ordinates rotating with the blades, show large positive and negative values of v(θ) particularly in the tip region. The circumferential average of v or the value at the blades can, therefore, be small while the average of v2 can be significant. The flow in the tip region of a rotor with a finite number of blades can include motion of the tip vortex towards, rather than away from, the axis of rotation .

Assuming that a and v are continuous through the rotor, Eq. () applies on SD as well as SU, and Eq. () simplifies to Tρ=2SDw2+λwxdS. All that remains to recover Eq. () is to determine when the blade-element version of Eq. () is also valid.

To determine the blade-element version of Eq. (), we consider two variants of a cylindrical annular CV: one with the downstream face just upwind of the rotor (V1) and one with it just downwind (V2), as shown in Fig. . In the former case, the flow everywhere in the control volume and on the bounding control surface (S1) is irrotational, and since there is no body enclosed by this CV, the thrust is zero and we obtain 0=S1n12U2I-UUdS.

Annular CVs used for the blade-element force analysis. The downstream face of width dx is just upstream and just downstream of the rotor for V1 and V2, respectively.

As noted above, the circumferential velocity ahead of the rotor vanishes for a circumferentially uniform flow, in which case only radial and axial velocities contribute to the integral in Eq. (). Assuming these velocities to be continuous across the rotor, their contribution to the same integrand over S2 also vanishes and we obtain S2n12U2I-UUdS=2SDw2dS, where SD is the downstream face of the annular CV, V2. Since the trailing vortices pierce SD, the vortex terms contribute 2λwx to the integrand of the thrust integral, as in Eq. (). Vorticity also exists on the lateral boundaries of V2, which we denote SL- and SL+, and the complete thrust contribution of the vortex terms is 12ezS2-nU(x×Ω)+nΩ(x×U)dS=12SL+vΩθ-uθΩrxdS-12SL-vΩθ-uθΩrxdS28=+2SDλwxdS. The integrals over SL- and SL+ do not vanish if the vortex lines and streamlines are coincident (as they did on SD) because diffusion is not negligible in the immediate vicinity of the blades where vorticity is generated. We will nonetheless neglect these contributions without claiming their insignificance, as this simplification facilitates the recovery of the KJ equation. Thereafter combining Eqs. () and (), the expression for blade-element thrust becomes 1ρdTdx=202πw2+λwxxdθ, which is, at last, identical to Eq. ().

Recall that the classical derivation of Eq. () requires the following assumptions: the trailing wake is thin, i.e. it occupies zero volume; it rotates rigidly with the rotor; and the axial and radial velocities are continuous across the disk. We have arrived at Eq. () using the latter two of these three assumptions, but we have avoided the thin-wake assumption; the new derivation permits the flow to be rotational in the wake. We have, however, invoked the circumferential uniformity assumption early; this was not required in the classical derivation of Eq. () but is necessary to recover the KJ blade-element thrust equation.

Kutta–Joukowsky equation

The sum of blade circulations at radial station x can be calculated as an area integral of radial vorticity over Sa, or a closed line integral over Ca, where only the downstream circular section of Ca lends a non-zero contribution.

Circumferential uniformity is approached as N, but Eq. () is approximately true for finite N when λxw. This occurs for the power-producing region of the blades at sufficiently high λ.

In blade-element analysis, the Kutta–Joukowsky equation is usually introduced as an assumption, but this can only be valid as N and/or λ. This result is directly deducible from Eq. () regardless of the method of its derivation.

Recall that we neglected the contribution of the vortex terms on the lateral boundaries of the annular CV to arrive at Eq. (). Although this assumption led to a result in agreement with the classical approach, this agreement is not sufficient justification. Rather, the classical derivation is simply blind to this effect due to its neglect of vorticity. Future work is needed to investigate the possible significance of the neglected terms.

Summary and conclusions

This paper describes the application of vortical impulse theory to determine the thrust on a steadily rotating wind turbine in a steady, spatially uniform wind. The principal attraction of the impulse approach is the absence of pressure in the force equations, which allows them to be applied anywhere in the flow, including immediately behind the rotor. We assume that the vorticity field is steady when viewed by an observer rotating with the blades, so the vortex lines and streamlines in the near wake are aligned in this frame. Subsequently, we rederive the Kutta–Joukowsky (KJ) equation for blade-element thrust under the special condition of circumferential uniformity in the wake (i.e. as the number of blades or the tip speed ratio tends to infinity). Our derivation improves upon the classical derivation – which is based on the Bernoulli equation – by allowing the wake to be rotational and three-dimensional and is exact for any amount of flow expansion through the rotor. In addition, intermediate steps en route to the KJ equation provide thrust expressions which allow for circumferential non-uniformity and trailing vortex sheets of finite thickness.

The derivation also yields insight into the radial velocity based on Eq. (). When the wake is circumferentially uniform, the magnitudes of the radial velocity and the axial induction factor must be equal somewhere on the plane containing the rotor, and we cite previous studies that show this to occur near the rotor tip across a wide range of thrust coefficients. Equation () was verified in the simulations of .

An important avenue for future work is an investigation of the effect of vorticity on the lateral boundaries of the annular control volumes in the immediate proximity of the blades. The blade-element KJ equation was recovered only when these terms were neglected, but the validity of this neglect was not established. This complication is not visible to the classical derivation of the KJ equation due to its neglect of vorticity.

Assume that, at all times, the fluid elements along a streamline passing through the rotor remain irrotational. The trailing vorticity in the wake is assumed to be infinitely thin. The unsteady Bernoulli equation along this streamline is ϕt+12(u2+v2+uθ2)+p=C, where C is a constant, and ϕ is the scalar potential. If the wake rotates rigidly with the blades, the unsteady term can be expressed in terms of the circumferential velocity and the rotation rate, λ: ϕt=θtϕθ=-λxuθ, where the negative sign arises because θ is evaluated in the moving frame, i.e. θ/t=-λ, and the relationship uθ=(1/x)ϕ/θ comes directly from the definition of the scalar potential, U=ϕ. We now apply Eq. () at locations just up- and downstream of the rotor to solve for the pressure difference across the rotor, Δp. Assuming radial and axial velocities to be continuous across the disk, only the circumferential velocity contributes to Δp: Δp=12uθ2-λxuθ. Substituting uθ=-2w, and integrating Δp over the rotor, we recover Eq. (): Tρ=SDΔpdS=2SDw2+λwxdS. Since the Bernoulli equation along any one streamline is independent of neighbouring streamlines, the blade-element version of this equation should also hold, i.e. 1ρdTdx=202πw2+λwxxdθ, which is Eqs. () and () in the main text.

The impulse equation for thrust in steady rotating coordinates

When the polar co-ordinates are attached to the blades rotating at λ, the velocities corresponding to those used in the main text are (0,λx,1) well upwind, (v,λx,1-a) immediately upwind of the rotor, and (v,uθ,1-a)=(v,2w+λx,1-a) downwind. Equation () requires modification to be expressed in terms of velocities and vorticities evaluated in the rotating frame of reference, U and Ω, respectively. The momentum integral in Eq. () is replaced by ddtVUdV=ddtVUdV+V2Λ×UdV+VΛ×(Λ×x)dV.

In Noca's derivation of Eq. () from Eq. (), the pressure, P, was removed by using the Navier–Stokes equations in the form of his Eq. (3.44): Ut=-P+12U2+U×Ω+T.

For a steadily rotating CV, Eq. () becomes Ut=-P+12U2+U×Ω-2Λ×U-Λ×(Λ×x), as given, for example, by Eq. (3.2.10) of . Clearly, the Coriolis term is the penultimate one and the last is the centrifugal. These have now to be included in Eq. (). The right side of Eq. () appears explicitly in the second term of Noca's Eq. (3.49) as Sx×n×[(P+12U2)+U×Ω+T]dS, which means that, in addition to swapping U for U and Ω and Ω in Eq. (), use of a noninertial CV requires the extra terms -Sx×n×[2Λ×U+Λ×(Λ×x)]dS. The general form of the impulse formulation, Eq. (), in a steadily rotating frame of reference thus becomes Fρ=-1N-1ddtVx×ΩdV+Sn12U2I-UU-1N-1U(x×Ω)+1N-1Ω(x×U)dS-V2Λ×UdV-VΛ×(Λ×x)dV-1N-1Sx×n×2Λ×UB6+Λ×(Λ×x)dS+Fv. Fv will be neglected as before, as we expect viscous effects to be no more important for rotating than stationary wakes. It is assumed that the vortical wake appears stationary relative to the rotating blades, such that the time derivative of impulse – the first integral on the right-hand side of Eq. () – vanishes. The Coriolis force can only have radial and circumferential components and the centrifugal force is radial. The radial components cannot contribute to an axial (or indeed to a circumferential) force, and thus the volume integrals of the pseudo-forces at the end of the first line and start of the second line of Eq. () do not contribute to T=ezF. The axial component of the centrifugal contribution to the pseudo-force surface integral in the second line of Eq. () also vanishes by symmetry, and the remaining Coriolis term can be evaluated using the vector identity a×(b×c)=(ac)b-(ab)c. With SD denoting the downwind face of the CV, the net contribution, which lies in the axial direction, is -2SDλwxdS, because the term involving λ2x2 will be equal and opposite at the upwind and downwind faces of the CV.

Thus the reduced form of Eq. () is Tρ=ez[12SnU2dS-SnUUdS-12SnU(x×Ω)dSB8+12SnΩ(x×U)dS-2SDezλwxdS]. Equation () is applied to the CVs shown in Fig. , but now rotating at the same λ as the rotor. Rotation does not alter the MD term in Eq. () or the contributions to the kinetic energy term from the axial and radial velocities, Eqs. () and () respectively. The circumferential velocity contribution is altered to 12Suθ2dS=2SD(w2+λwx)dS. It will again be assumed that w=0 outside the wake.

The vortex terms are better expressed in the inertial frame, Ω=Ω+2Λ, because non-zero values of Ω indicate fluid elements that have experienced viscous shear, whereas the apparent vorticity Ω has no such physical interpretation. The vortex terms are equivalently expressed as [-12SnU(x×Ω)dS+12SnΩ(x×U)dS]B10-SnΛ(x×U)dS. For the vortical wake to appear stationary in the rotating frame, Ω must be parallel to U everywhere that vorticity is non-zero – that is, vortex lines and streamlines are aligned in the rotating frame – by which nΩ(x×U)=nU(x×Ω) and the two integrals in square brackets identically cancel, leaving only -SnΛ(x×U)dS=SD2λwxdS for T. Combining these results leads to B12Tρ=12SDv2dS-12SDa2dS+2SDw2+λwxdS+2SDλwxdS-2SDλwxdS.

Curiously, three integrals with the integrand 2λwx have arisen from three different origins: the kinetic energy term and the vortex terms each yield the same integral in the positive sense, and the Coriolis term lends an equal and opposite contribution. The choice of which two to cancel is arbitrary, but either choice leads to the final equation: Tρ=12SDv2-a2xdx+2SDw2+λwxdS, which is identical to Eq. ().

Data availability

Author contributions

The work was conceived, executed, and written up equally by both authors.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

David H. Wood's contribution to this work is part of a research project on wind turbine aerodynamics funded by the NSERC Discovery Program. Eric J. Limacher acknowledges receipt of an NSERC Post-Doctoral Scholarship. The authors thank the anonymous referees and Gijs van Kuik for their valuable comments.

Financial support

This research has been supported by the NSERC (Discovery Program and PDS).

Review statement

This paper was edited by Jens Nørkær Sørensen and reviewed by three anonymous referees.