# Time-Varying Incident Velocity

The previously described unsteady airfoil solutions have considered the local free stream velocity to be constant. However, in forward flight a typical blade element on a rotor Figure 8.23 Airloads produced for a series of downstream moving sharp-edged vertical gusts, (a) Lift, (b) Pitching moment about 1 /4-chord.  Figure 8.24 Airloads produced for a series of upstream moving vertical gusts, (a) Lift, (b) Pitching moment about 1 /4-chord.

 Vb

 Yb = 0

 V+ v(t)  X = long

Figure 8.25 A time-varying free-stream velocity results in a nonuniform convection speed of the shed wake.

an arbitrary pitch axis. Essentially, Johnson’s conclusion is that the approximation using the Theodorsen function with the local reduced frequency is appropriate for flow oscillation amplitudes of up to 70% of the mean velocity. For small flow oscillation amplitudes, the Theodorsen function calculated with the reduced frequency argument based on the mean velocity will be accurate enough, which effectively means neglecting the unsteady free stream effects.

The theoretical problem of a time-varying free stream velocity on the unsteady aero­dynamic response has been examined by van der Wall & Leishman (1994), both from the classical frequency domain approach and also from the time domain. This work reempha­sizes that one of the significant effects of an oscillating incident flow velocity is a nonuniform convection velocity of the shed wake vorticity behind the airfoil. Van der Wall & Leishman

As mentioned previously, in the time domain the effect of a time-varying velocity in an incompressible flow can be modeled using Duhamel superposition with the Wagner function. Again, the appropriate apparent mass forces must be included in the solution, even though they will be small for the frequencies and amplitudes of the fluctuating velocity typically found on a helicopter. In this case

jrb /.. d(Va)

C, W = -[h + -4r – aba)

+ f )/а(0)ф(і) + J (j)(s – a)da) . (8.69)

The distance traveled through the flow is now proportional to the area under the velocity versus time curve as given by Eq. 8.4. Van der Wall & Leishman have compared results obtained with this time domain method with Isaac’s frequency domain results, where the agreement is essentially exact.

Results calculated using Eq. 8.69 are shown in Fig. 8.26 for the lift on an airfoil at a constant angle of attack, a, experiencing a harmonic oscillation of velocity of the form V(t) = Voo(l + A sin cot) at a reduced frequency, k, of 0.2 with A = 0, 0.2, 0.4, 0.6, and 0.8. The Duhamel integral was solved by means of a finite-difference approximation using the numerical scheme discussed in Section 8.14.1. The lift in Fig. 8.26 is shown as a fraction of the quasi-steady lift, (C[)qs = 2na. Also shown in this figure are results computed by means of a CFD method based on solutions to the Euler equations. It will be seen from Fig. 8.26 that at a given time the primary effect of unsteadiness is either an increase or attenuation Figure 8.26 The effects of a nonsteady free-stream velocity on the unsteady lift for an airfoil with a constant angle of attack. (CFD calculations courtesy of Arun Jose.)

of the unsteady lift compared to the quasi-steady result, which is a result already noted for oscillations in angle of attack. For 0°«wt <90°, the local velocity (and thus the product Va) is increasing. Here, the shed circulation has the opposite sign to the airfoil circulation and so the unsteady lift is less than the quasi-steady value. The higher local velocity here means that the shed wake is further away from the airfoil in a given time, and this effect contributes to a more rapid buildup in the circulatory lift compared to the case where the wake convection speed would be assumed uniform. For 90° <cot< 270°, the local velocity is decreasing. Here, the shed circulation has the same sign as the airfoil circulation, which tends to increase the lift over the quasi-steady value. Furthermore, for 180° < cot < 360° a velocity lower than the mean incident velocity means that the shed wake is closer to the airfoil; this then enhances the basic unsteady effect and the lift quickly becomes much greater than the quasi-steady values. Finally for 270° < cot < 360°, the local velocity is increasing again. With the change in sign of the wake circulation and the increasing wake convection velocity, the lift drops quickly to become close to the quasi-steady value again. Notice that in Fig. 8.26 there is good agreement between the incompressible theory and the CFD method, despite the slight phase differences. Other theoretical work on modeling unsteady free stream velocity effects is given by Ashley et al. (1952) and Friedmann & Yuan (1977). The extension of this problem to subsonic compressible flow is not without its difficulties, but one approach is discussed later in Section 8.17.