Wind Energy Handbook. Michael Barton Graham
Чтение книги онлайн.

Читать онлайн книгу Wind Energy Handbook - Michael Barton Graham страница 49

Название: Wind Energy Handbook

Автор: Michael Barton Graham

Издательство: John Wiley & Sons Limited

Жанр: Физика

Серия:

isbn: 9781119451167

isbn:

СКАЧАТЬ Schematic illustration of the radial and axial variation of axial velocity in the vicinity of an actuator disc, a=13.

      Note that the bound vorticity (being the circulation on the rotor blades in response to the incident and induced flow) induces zero rotation at the disc and decays axially up and downstream. The discontinuity in tangential velocity at the disc is because the idealised changes are assumed to take place through a disc of zero thickness. In reality the azimuthal velocity rises rapidly but continuously as the flow passes through the rotor blades, which sweep through a disc and influence region of finite thickness as shown in Figure 3.5.

      At the disc itself, because the bound vorticity induces no rotation and the wake cylinder induces no rotation within the wake cylinder either, it is only the root vortex that does induce rotation, and that value is half the total induced generally in the wake. Hence the root vortex induced rotation that is only half the rotational velocity is used to determine the flow angle at the disc. At a radial distance equal to half the disc radius, as an example, the axial variation of the three contributions is shown in Figure 3.10.

      The rotational flow is confined to the wake, that is, inside the cylinder, and tends asymptotically to 2a′Ω well downstream of the rotor. There is no rotational flow anywhere outside the wake, neither upstream of the disc nor at radial distances outside the wake cylinder. Because of this there is no first order transverse effect of the proximity of a ground plane on the downstream convection of the vortex wake of a wind turbine as there is on the trailing vortices of a fixed wing aircraft. The rotational flow within the wake cylinder decreases radially from the axis to the wake boundary but is not zero at the outer edge of the wake, therefore there is an abrupt fall of rotational velocity across this cylindrical wake surface vortex sheet.

Graph depicts the axial variation of tangential velocity in the vicinity of an actuator disc at 50% radius, a=13, λ = 6. Graph depicts the axial variation of tangential velocity in the vicinity of an actuator disc at 101% radius, a=13, λ = 6.

      3.4.8 Axial thrust

      The axial thrust T on the disc can be determined using the Kutta–Joukowski theorem:

StartFraction italic d upper T Over italic d r EndFraction equals rho normal upper Gamma upper V

      where V is the tangential velocity component at the disc. If V = rΩ(1 + a′), then, using Eq. (3.32) ignoring any additional inflow at the disc caused by the centrifugal pressure reduction due to wake swirl discussed at the end of Section 3.3.2:

      Integration of Eq. (3.41) over the entire disc gives the thrust coefficient as

      (3.42)upper C Subscript upper T Baseline equals 4 a left-parenthesis 1 minus a right-parenthesis

      That is, the same as for the simple momentum theory and so in balance with the rate of change of axial momentum. Note that if the induced tangential velocity arΩ is included in V as it is in blade‐element/momentum (BEM) theory and the blade circulation is constant from the axis to the tip, there is a singularity in the axial force on the blade section at the axis as there is also at the outer tip. This points to the failure of a simple constant strength bound vortex model at the blade ends as discussed in the section on tip‐loss corrections.

      3.4.9 Radial flow and the general flow field

      The radial velocity field that is predicted is largest on any given streamline at the actuator disc rising from zero at the axis to a weak logarithmic infinite value at the edge of the disc, which is the path of the blade tips. The infinite radial velocity at the edge is associated with non‐zero disc loading right up to the edge. This is not realistic, being a consequence of assuming the rotor to consist of an infinite number of blades whose effect is ‘smeared’ uniformly over the disc, but being a weak singularity does not significantly affect the rest of the flow field. In applying the more detailed BEM theory the tip region is corrected by a tip correction factor to recognise that in reality the blade loading must fall to zero at the blade tips.

      An alternative method of deriving the velocity field of the actuator disc has been given more recently by Conway (1998). This method takes the approach of building up the flow field from a sum of Bessel functions that are fundamental СКАЧАТЬ