aersp 305w aerospace technology...

14
AERSP 305W Aerospace Technology Laboratory Laboratory Section 17 Laboratory Experiment Number 3 2D Airfoil Testing of S805 Airfoil February 10, 2011 Performed in Room 8 Hammond Building Brett Davis Lab Partner’s Names : Titos Gosalvez Kevin Harrigan James Trexler Zachary Watson Lab TA: Brian Wallace Course Instructor: Richard Auhl

Upload: doanthien

Post on 26-Jun-2018

224 views

Category:

Documents


0 download

TRANSCRIPT

AERSP 305W Aerospace Technology Laboratory

Laboratory Section 17

Laboratory Experiment Number 3 2D Airfoil Testing of S805 Airfoil

February 10, 2011 Performed in Room 8 Hammond Building

Brett Davis

Lab Partner’s Names: Titos Gosalvez Kevin Harrigan James Trexler

Zachary Watson

Lab TA: Brian Wallace

Course Instructor: Richard Auhl

Brett Davis Lab 3 Page 2 of 14

Abstract

The objective of this experiment was to determine the 2D lift and drag coefficients of an S805

airfoil. Pressure taps in the airfoil were used to measure the local pressures on the airfoil’s

surface. These were used to calculate the lift coefficients. A hot-wire anemometer was used

to take a velocity survey of the airfoil’s wake in order to calculate the drag coefficients. The

coefficients were calculated from -18o to 18o at every two degrees. The angle of attach at

which stall and that at which the flow reattaches were also determined. The data from this

experiment coincided with the data presented in the NREL report on the S805 airfoil. This

experiment successfully determined the 2D lift and drag coefficients for the S805 airfoil and

was consistent with the NREL report.

Brett Davis Lab 3 Page 3 of 14

Introduction

The objective of this experiment was to determine the 2D lift and drag coefficients of an S805

airfoil. The data obtained in this experiment was compared to prior experiments performed

by the National Renewable Energy Laboratory (NREL). With pressure readings from the

pressure taps in the airfoil, the local pressure coefficients were determined using Equations

(1a, 1b) for the upper and lower surface of the airfoil.

The normal force coefficient was calculated using the equations below.

(2)

(3)

Equation (3) was obtained by assuming a thin airfoil. From the normal force coefficient, the

2D lift coefficient was obtained with Equation (4) by assuming a small angle of attack.

(4)

To determine the 2D drag coefficient, Equation (5) was evaluated using Simpson’s Rule,

Equation (6), and the velocities measured by a hot-wire anemometer

.

(5)

(6)

During Lab 1, the wind tunnel’s venturi was calibrated using Equation (7). This was used to

hold the wind tunnel at a constant Reynolds number.

(7)

Brett Davis Lab 3 Page 4 of 14

Experimental Procedure

To determine the 2D lift and drag coefficients for the S805 airfoil a series of experiments were

run in the 2x3 foot wind at a Reynolds number of 1,000,000 and at angles of attach ranging

from negative 18o to positive 18o. Time constraints only allowed each lab group to run tests

at two angles of attack. The data collected was shared between all lab groups. This report

concerns the experiments performed at 1o and -18o angle of attack. The details of the

procedures are outlined below.

A solid aluminum airfoil with a constant chord of 19.685” and a span of 24” was used to

perform these experiments. The airfoil had 48 static pressure ports on the upper surface and

49 ports on the lower surface. 17 of the upper and 17 of the lower ports were connected to a

manometer bank as seen in Table 1.

Table 1. S805 Pressure Tap Locations

X in. Y in. X in. Y in.

1 0.00 0.00 19

2 0.16 0.27 20 0.04 -0.08

3 0.39 0.46 21 0.24 -0.22

4 0.63 0.59 22 0.47 -0.30

5 1.26 0.86 23 0.79 -0.37

6 3.94 1.48 24 1.97 -0.55

7 6.89 1.76 25 4.92 -0.79

8 8.86 1.76 26 7.87 -0.88

9 9.45 1.73 27 9.84 -0.87

10 10.04 1.68 28 12.80 -0.77

11 10.63 1.61 29 14.17 -0.65

12 1.22 1.52 30 14.76 -0.58

13 11.81 1.42 31 15.35 -0.50

14 13.78 1.09 32 15.94 -0.41

15 16.73 0.59 33 16.34 -0.34

16 18.11 0.33 34 16.93 -0.25

17 19.29 0.07 35 17.72 -0.13

18 36 19.29 -0.01Atmosphere

Manometer

Port Number

Manometer

Port Number

Tap Location

Upper Surface Lower Surface

Tap Location

T.S. Static

Pressure transducer A was connected to the pitot static probe in the wind tunnel test section,

and transducer B measured the static pressure drop over the venturi. Transducers A and B

were set at a span of 4.3 and previously calibrated to 2.945 PSF/volt and 2.985 PSF/volt

Brett Davis Lab 3 Page 5 of 14

respectively. Downstream from the airfoil a hot-wire anemometer was mounted on a sting

and motorized traverse system to measure the wake velocities as seen in Figure 1.

Figure 1. Experimental set up, location of hot-wire probe and airfoil

The hot-wire anemometer was previously calibrated to yield Equation (8).

(8)

To run the wind tunnel at a constant Reynolds number of 1 million, the venturi pressure drop

was held constant at 3.47 volts on transducer B. This pressure was determined from the

venturi calibration performed in Lab 1. The calibration constants in Lab 1 were determined

for the front, middle and rear of the test section as seen in the Appendix A. This experiment

was performed between the rear and middle locations, so these values were averaged to get

K = 0.899.

The airfoil was set at -17.95o (as close to -18o as possible) angle of attack and LabVIEW was

used to traverse the hot-wire probe vertically across the test section to measure the velocities

in the airfoil’s wake at every 0.125” over the 24” wind tunnel height. The probe made 2000

measurements at each location with a sample rate of 2 kHz. With the manometer bank tilted

at 19.8o the height of the fluid in each bank was recorded. When the hot-wire wake survey

was completed, the airfoil was rotated to 0.96o (approximately 1o). The hot-wire survey was

run again, and the fluid levels in each bank were recorded while the manometer bank was

tilted at 19.7o.

Hot-wire probe

S805 Airfoil

Traverse

direction

Brett Davis Lab 3 Page 6 of 14

General observations about the stall point of the airfoil were made with the wind tunnel

running at 50%. The airfoil was rotated from -18o to stall while observing the manometer

bank. After the airfoil stalled at 20o, the airfoil was rotated back through lower angles of

attack until the flow reattached at 17.6o.

Results and Discussion

The data collected during this experiment and during the experiments performed by other lab

groups was combined to determine the 2D coefficients of lift and drag for the S805 airfoil.

The results of these experiments were compared to the data in the NREL Design and

Experimental Results for the S805 Airfoil report to assess the findings of this experiment.

During the lab, the room temperature, pressure and air density were recorded in Table 2.

Table 2. Room conditions during experiments

During Experiment Temperature (F) Pressure (in Hg) Density (slug/ft^3)

α = 1 deg 66.4 29.42 0.002305

α = -18 deg 67.1 29.42 0.002302

A time trace of the velocity at the start of surveys at 1 and negative 18 degrees was taken at

the top of the test section. As seen in Figure 2, the turbulence intensity was much larger

when the airfoil was at negative 18 degrees angle of attack than at positive 1 degree. These

results would be expected, as the larger angle of attack should produce more turbulence

downstream from the airfoil.

Figure 2. Velocity time trace downstream from S805 airfoil

70.00

80.00

90.00

100.00

110.00

120.00

130.00

0.00 0.20 0.40 0.60 0.80 1.00

Ve

loci

ty (

ft/s

)

Time (sec)

1 Deg

Re = 1x10^6 Ti (1o) = 0.003464

Ti (18o) = 0.016203

Brett Davis Lab 3 Page 7 of 14

From the hot-wire anemometer’s measurements, the average velocity at each position was

calculated. These velocities were plotted as a function of vertical position to form the wake

profile at 1o in Figure 3 and at -18o in Figure 4. These profiles show that the -18o angle of

attach created a much larger velocity deficit than the 1o orientation. This large wake results in

a larger coefficient of drag, 0.2412, than the smaller wake at 1o, which had a cd of 0.0083.

These differences are consistent with what is expected.

Figure 3. Wake velocity profile at 1o angle of attack

Figure 4. Wake velocity profile at -18o angle of attack

0

5

10

15

20

88 90 92 94 96 98 100

Ve

rtic

al P

osi

tio

n (

inch

es)

Velocity (ft/s)

Re = 1x10^6 α = 0.96o

cd = 0.0083

0

5

10

15

20

0 20 40 60 80 100 120

Ve

rtic

al P

osi

tio

n (

inch

es)

Velocity (ft/s)

Re = 1x10^6 α = -17.95o

cd = 0.2412

Brett Davis Lab 3 Page 8 of 14

Using the pressure readings taken from the manometer bank, the pressure distribution over

the airfoil was plotted. In Figure 5, it can be seen that at α = 1o the upper surface of the airfoil

experienced negative pressures and the lower surface had positive pressure, indicating

positive lift resulting in a cl of 0.3344. The data from this experiment closely follows the data

from the NREL report as can be seen from Figure 5. There was a favorable pressure

gradient to about 50% of the chord on the upper surface and 65% on the lower surface,

followed by adverse pressure gradients, which is close to the design specifications in the

NREL report.

Figure 5. Pressure distribution over S805 airfoil at 1o angle of attack

The pressure distribution over the airfoil at -18o can be seen in Figure 6. From this graph, it

can be seen that the airfoil generated negative lift, as the upper surface experienced positive

pressures, and the lower surface had negative pressures. There was not a large pressure

gradient on the lower surface in this orientation, while the upper surface had a large adverse

pressure gradient up to 0.5% of the chord followed by a mostly favorable pressure gradient.

At this angle, the cl was -0.6714. This data did not follow the NREL data as closely, although

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

0.000 0.200 0.400 0.600 0.800 1.000

Cp

x/c

Experimental Upper Surface

Experimental Lower Surface

NREL upper

NREL lower

Re = 1x10^6 cl = 0.3344

αExperiment = 0.96o

αNREL = 1o

LE TE

Brett Davis Lab 3 Page 9 of 14

the general trend was consistent. The angles of attack did not exactly match the NREL data

in this case which could have contributed to the differences in the data.

Figure 6. Pressure distribution over S805 airfoil at -18o angle of attack

The compiled data from each lab group was used to generate a plot of the coefficient of lift as

a function of the angle of attack in Figure 7. The linear portion of the cl versus α curve was a

near perfect match to the NREL results. In the non-linear regions of the curve, the

coefficients of lift from this experiment were slightly greater than those from the NREL report

for positive angles of attack. At -16o and -14o the coefficients of lift are slightly less than

expected possibly due to measurement errors. Overall, the cl versus α curve is consistent

with expectations and the NREL findings.

-1.500

-1.000

-0.500

0.000

0.500

1.000

1.500

0.000 0.200 0.400 0.600 0.800 1.000

Cp

x/c

Experimental Upper Surface

Experimental Lower Surface

NREL lower

NREL upper

Re = 1x10^6 cl = -0.6714

αExperiment = -17.95o

αNREL = -18.1o

TE LE

Brett Davis Lab 3 Page 10 of 14

Figure 7. Coefficient of lift vs. angle of attack of S805 airfoil

The compiled data was also used to create a plot of the coefficient of lift versus the coefficient

of drag in Figure 8. The plot shows the expected relationship between cl and cd, and the

NREL data appears to be consistent with the findings of this experiment. However, the NREL

data was over a much smaller range of drag coefficients. Figure 9 shows the data over the

range of cd presented in the NREL report. At small values of cd this experiment’s data is

relatively close to the NREL data. The differences between the two sets of data are likely

due to the difficulty in determining these very small drag coefficients. In this experiment, the

cd was evaluated by integrating the wake profile using Simpson’s Rule. The bounds of

integration for each angle of attack had to be manually selected because the velocity varied

slightly over the height of the wind tunnel. The approximations in this approach could be

responsible for the differences seen in Figure 9.

-1

-0.5

0

0.5

1

1.5

-20.00 -15.00 -10.00 -5.00 0.00 5.00 10.00 15.00 20.00

Cl

Angle of Attack (deg)

Experimental

NREL

Re = 1x10^6

Brett Davis Lab 3 Page 11 of 14

Figure 8. Coefficient of lift vs. coefficient of drag of S805 airfoil

Figure 9. cl vs. cd*1000 of S805 airfoil at small cd values

-1

-0.5

0

0.5

1

1.5

0 0.05 0.1 0.15 0.2 0.25

cl

cd

Experimental

NREL

Re = 1x10^6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

1.2

0 5 10 15 20

cl

cd *1000

Experimental

NREL

Re = 1x10^6

Brett Davis Lab 3 Page 12 of 14

With the tunnel running at 50%, the airfoil stalled at approximately 20o, however, the

protractor maxed out at 20o so this measurement could be slightly off. Figure 10 shows the

manometer bank when the airfoil was at 6o to 20o angle of attack.

Figure 10. Manometer bank showing airfoil pressures from 6o until stall

The flow reattached when the airfoil was rotated back to 17.6o. The flow separated at a

higher angle of attack than it reattached at because the flow needs more energy to reattach.

Once the flow is attached, its momentum resists separation until the adverse pressure

gradient over comes it. Once separated, the flow needs to gain enough energy to initially

reattach which is why the two angles were different.

6.0o 11.2

o 13.5

o 18.9

o ~20

o

(stall)

Brett Davis Lab 3 Page 13 of 14

Conclusions

The 2D lift and drag coefficients for the S805 airfoil found from this experiment were

consistent with those presented in the NREL report. The coefficient of lift as a function of the

angle of attack yielded the same relationship as that in the NREL report with only a slight

error at large angles of attack. The lift and drag coefficients’ relationship was also consistent

with the NREL data. Some errors occurred at small coefficients of drag likely due to the

difficulty in computing such small amounts of drag. These differences were insignificant

when considered over the full range of data taken. General observations indicated that the

angle of stall is greater than the angle required to reattach the flow, which is consistent with

what was expected. Overall, the experiment successfully determined the 2D lift and drag

coefficients for the S805 airfoil.

Brett Davis Lab 3 Page 14 of 14

The Appendix

Appendix A. 2x3 ft. wind tunnel venturi calibration

∆P = 0.9976q - 0.0937

∆P = 0.8976q + 0.1094

∆P = 0.8766q + 0.0078

∆P = 0.9897q

0

5

10

15

20

25

30

0 5 10 15 20 25 30

Co

ntr

acti

on

Pre

ssu

re D

rop

(P

SF)

Test Section Dynamic Pressure (PSF)

Front

Middle

Rear

Theoretical