nJAsA - NASA Technical Reports Server (NTRS)

1 downloads 135 Views 6MB Size Report
East Hartford,. Connecticzlt ..... 5), Widnall and Wolf (ref. 6) and .... point. Also, local radial station in vortex radius of vortex viscous core radius of open .... point of the gust wavefront with the airfoil ..... to the atmosphere through an exhaust tower.
NASA

Contractor

Rotor-Vortex

Robert

CONTRACT OCTOBER

H. Schlinker

Report

3744

Interaction

Noise

and Roy K. Amiet

NASl-16392 1983

LOAN COPY: RETURN TO AFWL TECHNICAL LIBRARY KIRTLAND AFB, N.M. 87117

25th Anniversary 1958-1983

nJAsA

TECH LIBRARY KAFB, NM

NASA

Contractor

Report

Rotor-Vortex

Robert

H. Schlinker

Interaction

and Roy K. Amiet

United Technologies Research East Hartford, Connecticzlt

Prepared for Langley Research Center under Contract NAS l- 16392

runsn

National Aeronautics and Space Administration Scientific and Technical Information Branch 1983

3744

Center

Noise

Rotor-Vortex

Interaction \

Noise

TABLE OF CONTENTS Page SUMMARY............................... INTRODUCTION.

...........................

.......................... Background. .... ; ............... Previous Investigations ..................... Present Investigation LIST OF SYMBOLS. . . . . . . . . . . . . . . . . . . . . . . . . . .

6

ANALYTICAL FORMULATION OF THE BLADE VORTEX INTERACTION MECHANISM . '.

10

Definition of the Problem and Description of Approach General Expression for Far-Field Sound. ............ The Airfoil Lift Response ................... The Velocity Upwash ...................... The Velocity Upwash for an Arbitrary Vortex Orientation Application of the Theory ................... DESCRIPTION OF THE EXPERIMENT. ................... Acoustic Research Tunnel. ................... .......................... ModelRotor Instrumentation ........................ Test Program. ......................... Qualification of Acoustic Test Procedure.

.....

....

21 22 23 26 29

...........

VELOCITY FIELD INCIDENT ON THE ROTOR . . . . . . . . . . . . . . . . Objective ............................ Approach ............................ Vortex Azimuthal and Axial Jlean Velocity Field. ........ Axial Velocity Defect in Two Dimensional Wake ......... Turbulence in Viscous Core. .................. Changes in Velocity Field Due to Contraction of Flow. Unsteady Upwash Due to Blade Vortex Interaction ........

iii

10 11 13 16 18 19 21

.....

31 31 32 32 34 35 35 37

TABLE OF CONTENTS (Cont'd) Page EXPERIMENTAL ASSESSMENTOF BLADE VORTEX INTERACTION NOISE :. ........................... Objective. ............... Effect of Blade-Vortex Interaction Blade Number Dependence. .................... Vortex Strength Dependence ................... Vortex Intersection Angle Dependence .............. Blade Pitch Angle Dependence .................. Mach Number Dependence ..................... ...................... Directivity Pattern. Assessment of Acoustic Spectrum and Pressure Signature

. . . .

.....

39 39 39 43 44 45 47 48 48 50

CONCLUSIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

53

APPENDIX A - Vortex

. . . . . .

54

Technique for Isolating Blade Vortex Impulsive Signatures . . . . . . . . . . . .

56

Velocity

APPENDIX B - Experimental Interaction

Field

Measurement

Technique.

REFERENCES..............................

57

TABLES...............................

60

FIGURES...............................

61

iv

A theoretical and experimental study was conducted to develop a validated first principles analysis for predicting noise generated by helicopter main-rotor A generalized noise prediction shed vortices interacting with the tail rotor. procedure was formulated in which the incident vortex velocity field, rotor conditions are input variables. The analysis includes geometry, and rotor operating chordwise and spanwise noncompactness, and treats oblique compressibility effects, intersections with the blade planform. Output from the analysis provides acoustic spectra and pressure signatures as a function of observer position. Assessment of the theory involved conducting a model rotor experiment which isolated the blade-vortex interaction noise from other rotor noise mechanisms. The experiment was conducted with a 0.76 m articulated rotor operating in the United Technologies Research Center Acoustic Research Tunnel. An isolated tip vortex, generated by an upstream semispan airfoil, was convected into the model rotor to generate blade-vortex interaction noise. Acoustic spectra, pressure signatures, and directivity were measured for different combinations of rotor tip Since input to the acoustic predicMach number, blade pitch, and blade number. tion required a knowledge of the vortex properties, blade vortex intersection angle and intersection station were determined using smoke visualization. Vortex strength and vortex core radius were documented with hot film probe measurements. Ingestion of the vortex by the rotor was experimentally observed to generate harmonic noise and impulsive waveforms. Assessment of the theory using the measured data verified-the linear dependence on blade number and vortex strength. The absence of noise for vortex trajectories normal to the blade was also predicted. For oblique intersections with a blade, changes in acoustic radiation, due to intersection angle modifications, were calculated to within 2 decibels of the measured A simplified scaling law was also developed to predict these changes. Low data. sensitivity to blade pitch angle modifications for typical rotor operating conditions was confirmed experimentally. The strong dependence on local Mach number was predicted although the sensitivity to Mach number changes was overpredicted at low frequencies. Measured blade vortex interaction directivity patterns confirmed the theoretical dipole acoustic source model including the location of the minimum In the plane of Absolute levels for acoustic spectra and directivity patterns were calthe blade. culated, without the use of empirical or adjustable constants, to within 5 decjlbels At high frequencies, sound pressure levels were underpredicted. at low frequencies. General features of the measured vortex interaction acoustic pressure signature were calculated but the exact waveform shape and amplitude were not predicted.

INTRODUCTION

Background The noise level at a given observer'position relative to a helicopter and within a specified frequency band represents the combined effect of a Some mechanisms occur only number of separate acoustic source mechanisms. others apply to both the main rotor and tail rotor, and on the main rotor, interactions occur linking the wakes of one rotor in more complex situations, The simultaneous occurrence of the with noise generated by the other rotor. acoustic source mechanisms has complicated the design of quiet helicopters. and contributed to empiricism and extrapolation of measured data when predicting full scale helicopter noise signatures. Due to the limitations of empirically based helicopter noise prediction an increasingly important objective of helicopter noise research procedures, has been the development of procedures based on first principles analyses. This effort has been reasonably successful for a number of noise mechanisms such as blade thickness, loading, turbulence ingestion, and boundary layerHowever, in the area of main-rotor shed vortices trailing edge interaction. The interacting with the tail rotor, first principle analyses did not exist. lack of a firm understanding of the tail-rotor vortex interaction mechanism In an created a need for an experimentally validated prediction procedure. the present joint experimental and theoreteffort to achieve this objective, ical research program was conducted to assess this noise mechanism. The blade-vortex interaction mechanism addressed here is considered to As illustrated in Figure la,under forward be potentially important. flight or vertical descent operating conditions vortices shed by the main rotor pass through the tail rotor resulting in the generation of noise. In the strength of this noise mechanism can be reduced by changes in principle, For example, altering the relative position of the the helicopter design. In addition, changing the main rotor and tail rotor is a possible solution. main-rotor blade tip shape may reduce the strength of the tip vortex or diffuse the vortex velocity field in a manner which reduces interaction noise. The decision to pursue such aerodynamic design changes is based on the The major objective ability to predict quantitatively the expected benefits. to develop and validate a first prinof the present study was, therefore, ciples prediction procedure which starts with the incident vortex velocity field as an input.

Previous

Investigations

studies of helicopter tail-rotor Experimental Studies - Few experimental Leverton, et al. (ref. 1) vortex ingestion noise are currently available. observed that the Lynx helicopter produced such noise during approach conditions when the tail rotor intersected tip vortices shed by the main rotor. The experimental results indicated that this noise mechanism can be subjectively annoying and can dominate over the main rotor and engine noise. In addition to full scale measurements, Pegg and Shidler (ref. 2) studied the parameters controlling main-rotor tail-rotor interaction noise using The model permitted varying tail rotor posia small scale helicopter model. tion relative to the main rotor in addition to the tail rotor speed and rotation direction. Experimental results showed that certain features of the tail rotor acoustic spectrum were attributable to interactions with tip vortices shed from the main rotor. A limited study of model scale vortex interaction noise was also conInteraction of an isolated tip vortex ducted by Paterson and Amiet (ref. 3). with the model rotor was observed to produce significant harmonic noise The impulsive acoustic waveform generated which extended to high frequency. by the interaction was found to be sensitive to vortex strength and position. models predicting the blade vortex inTheoretical Studies - Analytical teraction noise mechanism presently treat only the main rotor noise source. To help differentiate these analyses from the present formulation, a brief In this discussion of the main rotor interaction mechanism will be given. case the blade and vortex lie in different (or identical) parallel planes.. Furthermore, the blade leading edge is parallel to the vortex centerline in the most intense interactions (Fig.,lb). This condition is modeled by a two dimensional analysis assuming end effects can be neglected. On the other hand, if the vortex centerline is skewed relative to the airfoil leading edge (but still parallel to rotor plane), the vortex velocity field sweeps across the airfoil. Under this condition, a three dimensional analysis is required.

The above described two dimensional blade vortex interaction was treated The three dimensional interaction by Amiet (ref. 4) assuming linearized flow. was modeled by Widnall (ref. 5), Widnall and Wolf (ref. 6) and Adamczyk (ref. These analyses include the two dimensional problem as a special case 7). but they do not treat the present tail rotor interaction geometry which inIn this case, the rotor disk is approxivolves a nonplanar intersection. One objecstation. mately at 90" to the vortex plane at the intersection to develop the specific analyses tive of the present study was, therefore, needed to treat the tail rotor blade vortex interaction geometry.

3

Present

Investigation

objective of the present study was to develop Objective - The overall The analytan experimentally validated blade-vortex interaction analysis. ical phase of the program required the formulation of a generalized noise prediction procedure in which the incident vortex velocity field, rotor and far field observer position were rotor operating conditions, geometry, Empirical scaling factors were avoided in this first prininput variables. ciples analysis. Assessment of the theory required conducting a tractable acoustic experiment which isolated this noise mechanism from other operative rotor noise Also, detailed measurements of the incident vortex velocity mechanisms. field were obtained to provide the experimental input to the analysis. Confirmation of the theory for the experimental test conditions involved assessment of blade vortex interaction noise sensitivity to changes in the vortex field and the rotor operating conditions. transforms the incident vorTheoretical Approach - The present analysis tex velocity field into a spectrum of sinusoidal gusts incident on an isolated Linearized theory is then used to determine the unsteady blade rotor blade. The response is given by a generalized function lift response function. chordwise and spanwise noncompactwhich includes compressibility effects, ness, and treats skewed or normal gusts in the decomposition of the incident Input to the analysis requires a knowledge of the vortex vortex field. vortex viscous core radius, blade-vortex intersection angle, blade strength, vortex intersection station, blade tip Mach number, blade pitch angle, blade Output from the analyobserver location. chord, blade number, and far field sis provides acoustic spectra and pressure signatures as a function of observer position. arrangeExperimental Approach - Figures 2 and 3 show the experimental Here a tip vortex, generated ment employed to assess the present analysis. by an upstream semispan airfoil model, is convected into a model helicopter Although rotor operating in vertical ascent with a 9 m/set inflow velocity. this test configuration lacks the appearance of a tail rotor, it represents a tractable geometry which permits direct assessment of the blade-vortex inSelection of this geometry was based on an evaluation teraction mechanism. of the parameters controlling the noise mechanism and a detailed discussion is included in the description of the experiment. The experiment was conducted with a 0.76 m diameter articulated rotor. Acoustic spectra, pressure signatures and directivity were measured for different combinations of rotor tip Mach number, blade pitch, and blade number. Blade vortex intersection angle and blade intersection station were documented

with smoke visualization. viscous core turbulence film probe measurements.

Vortex intensity

strength and core radius in addition to and length scale were documented with hot

It was considered important to devise an experiment in which only one This was necessary since the vortex and rotor parameter varied at a time. For example, changing the rotor tip aerodynamic parameters were coupled. Mach number also changed the mean flow contraction ratio for the rotor orienvortex strength, intersection tation shown in Figure 2. Hence, the incident Selection of the test matrix, angle, and intersection station also changed. required documentation of these coupling effects. therefore,

LIST OF SYMBOLSAND NOMENCLATURE

a-

angle

between

b

semichord

B

rotor

blade

BPF

blade

passing

C

chord

cO

sound speed

hot

film

and total

velocity

vector

number frequency

c2

upstream

D

vortex

Ex"

combination

F

point

f

frequency,

f 00

function

g

normalized

G,Gm

single sided and multiple

H

vertical and tip

h

semispan

airfoil

I

distance

from rotor

i

J--i:

Jn'n J ux

sensor

semispar

airfoil

chord

length,

11.4

cm

core diameter of Fresnel

integrals

force

Bessel Eulerian

cycles/set of Mach number; pressure

see Eq. (13)

jump;

see Eq. (2)

spectrum for far field vortex intersections

separation distance of rotor blade

functions time

open jet blade

K

coefficient relating measurements

K1

vortex

tip

penetration tip

of the first scale

between

sound generated

by single

of semispan

airfoil

height

to vortex

intersection

station

and second kind

in autocorrelation velocity

strength

6

and voltage

for

hot

film

probe

normalized

K2

specific

K Y kx,k

Y

value

chordwise

strength,

K = Kl/cU 2

of k ; uS2y/o Y

and spanwise

gust wavenumbers

kxBm2/B2

X=

k

magnitude

i

unit

of wavevector

vector

effective

2 M %

in z direction lift

Mach number

rotor

blade

tip

similarity

P

pressure

AP

pressure

R

observer-vortex

RO

effective potential rotor

R1

function;

local

Mm

see Eq. (6)

Mach number

Mach number for

jump across

airfoil

distance

radius field

radius,

skewed gust

of exponential

decay function

autocorrelation functions for far field single and multiple vortex intersection

RIPS

rotor

r

distance from observer to airfoil-vortex Also, local radial station in vortex radius

rO

S

rB

revolutions

to vortex

sound generated

by a

per second

of vortex

radius of open jet induced by rotor Sears function;

applied

0.38 m

RlgRC0

rAj

0

kx/B2

kX*

k

vortex

viscous

core

flow

before

see Eq. (13)

7

(rB)

intersection

and after

(rA)

point.

contraction

two sided spectrum for far field and multiple vortex intersection SPL

sound pressure

level

in dB relative

dyne/cm2

time

t

time

U

fluid velocity relative to airfoil in present theory; c,M . Also, axial velocity in vortex core or in semispan airfoil two dimensional wake

uT

in Fourier

to 0.0002

by a single

T

'Nl"N2

increment

sound generated

velocity components normal two different orientations total

velocity

vector

V

sweep speed of gust

vp2’v3’v4

vortex

azimuthal

V N

vortex

velocity

73 V

0

vortex

sinusoidal

along

airfoil

film

probe

of the hot relative

sensing

film

element

in

sensor

to airfoil

vectors

component normal

azimuthal

maximum vortex

to hot

in the plane

velocity

sweep speed of gust

'r

summation

along

to blade

airfoil

relative

normal

to plane

planform to fluid

velocity velocity

gust velocity

of airfoil

wg W

Fourier

X,Y,Z

Car'resian

x1‘Yl

airfoil

a

angle

between

rotor

blade

u1 a2 6

amplitude

coordinates;

see Fig.

4

coordinates

semispan 447

of w g

gust

pitch

airfoil

and airfoil. angle

angle

of attack

Also,

constant

used in Eq. (22)

experimentallydeterminedblade vertical plane gF

angle

between hot

film

support

vortex

intersection

prongs

and sensor

angle

in

circulation Y

experimentally determined horizontal plane

6

Dirac

0

1

@2 e

0

8

V

delta

blade

defined

by Eq. (20a)

parameter

defined

by Eq. (20b)

vortex

intersection

function

parameter

observer

vortex

angle; angle;

see Fig.4 see Fig.

4

!J

Cartesian

coordinates;

freestream

see Eq. (37)

density

&P&3

Overbar

velocity

potential

observer

angle;

see Fig.

vortex

angle;

radian

frequency

radian

frequency

quantity

see Fig.

normalized

4 4

increment

in Fourier

by the semichord

9

summation b

angle

in

ANALYTICAL FORMULATION OF THE BLADE VORTEX INTERACTION MECHANISM

Definition

of the Problem

and Description

of Approach

The present analysis was developed for a blade vortex interaction geometry in which the rotor plane of rotation and the vortex plane are skewed. Consequently, the blade cuts through the vortex centerline during its motion as illustrated by the combined airfoil and vortex geometry shown in Figure 4. The orientation of the vortex relative to the airfoil is arbiSmall interaction angles between the airfoil plane and the vortex trary. centerline are, however, not permitted in the analysis. This restriction occurs because the intersection point of the airfoil and the vortex is assumed to be localized on the airfoil allowing end effects to be neglected. treated in references l-4 is excluded from Thus, the blade vortex orientation the present analysis. The present formulation is based on linearized theory blade modeled as a flat plate of infinite span. Because of flow assumption, steady loading is assumed to be small and with the unsteady blade vortex interaction noise mechanism. theory dictates that only the fluid velocity normal to the affects the airfoil loading. Also, because of the linearized the vortex velocity field and its trajectory are not tion, perturbations induced by the airfoil.

with the rotor the linearized is not coupled Standard airfoil airfoil plane flow assumpaltered by velocity

Analytical expressions used for the airfoil response functions are for the case of a flat plate airfoil of infinite span. A sinusoidal gust upwash is assumed incident on the airfoil. Thus, the airfoil response functions are similar to the classical Sears function (ref. 8). The only difference between the airfoil response functions used here and the Sears function is that the present solutions extend the Sears solution to the case of compressible flow. Also, using the similarity rule of Graham (ref. 9>, the gust can be skewed If the Mach number is set equal to zero and the relative to the airfoil. gust is assumed parallel to the airfoil, the present solution reduces exactly to that of Sears. A detailed derivation of the present gust response functions is given Separate solutions are given for the low and high frein references 10-14. The low frequency solution is a series type of solution and quency regime. the high frequency solution is based on a convergent iteration procedure; the solution produced is asymptotically correct as frequency goes to infinity. Comparison with exact numerical solutions of Graham (refs. 9, 15) shows that the present solutions give an error in the lift of at most a few percent. The maximum error occurs near the changeover point from the low to the high frequency solutions.

10

other tance finite

The model for the vortex velocity field is the same one used by several authors (refs. 16, 17). It reduces to a potential vortex at large disand approaches a solid body rotation in the core. This avoids the invelocity existing at r = 0 in a potential vortex.

General

Expression

For Far-Field

Sound

The response of the airfoil to the vortex velocity field will be treated by Fourier decomposing the velocity field and calculating the airfoil response to each resulting gust component. The inverse Fourier transform then gives the time behavior of the gust-vortex interaction. The airfoil is assumed to be a flat plate of infinite span lying in the If the airfoil X,Y plane, with the leading edge at x = -b (see Fig. 1). lift were being calculated, the assumption of infinite span would create problems; since the vortex field has only an r -l decay rate, an infinite lift would be obtained in general. This singularity would become apparent on attempting to perform the Fourier inversion of the expression for lift produced by a gust component; this has a k-l wavenumber behavior for small k. calculates the sound and the sound is However, since the present analysis dependent on the time derivative of the loading, an additional factor of k is introduced which eliminates the singularity at k = 0. There remains, in the calculated value of the average pressure level. however, a difficulty This is considered in the section titled Experimental Assessment of Blade Vortex Interaction Noise. Linearized flow is assumed. This means that the vortex remains fixed relative to the mean flow at infinity. Thus, on Fourier decomposition the z component of the vortex velocity field in the x,y plane can be written in terms of the gusts

This be written

gust will

lead

AP(x,,y,,kx,ky) The far-field x1, y1 is

pressure

P, =

wQ= %(kx,ky)e

-i[k,(x-Ut)+kyy/

to a pressure

jump across

= 2.rrpoUz (kx, ky)g(x,,kx,ky)

of a tranversepoint iFwz 41rcou2

ei[wt+p(Mx-olJ

dipole

(1) the airfoil

which

can

e-i(kyyI-kxUt) or force

F k e

(2)

iwt

at

e-icL[MxI-(xx,+P2yy,)/q] (3)

11

This is obtained from the usual expression for a dipole neglecting terms O(rl/r) where rl is a source (airfoil) the source-observer distance. The far-field ing F in Equation

(ref. 11) by dimension and r is

pressure due to the gust component wg is 3 with AP from Equation 2 noting that

found by replac-

(4)

w= k,U and integrating

over

the airfoil

irp,U2bkxZ p(kx,ky)=

This

surface.

gives i I k,Ut+p(MX-U)](5)

6(ky-~P2y/a)~(kx,ky)~(kx,ky,M)e

COD2

where p( k,, k,, M)= kJbg(x, -b

,k,,k,hl)

e-i@-x’u)xl

(6)

dx,

The exponential phase factor represents an effective lift per unit span. time effects for sound propagating from source in P accounts for retarded if this factor were set equal to 1, 1 would represent the actual to observer; lift per unit span. To find the far-field and ky gust components is

p(t)=-n

bzp,U2 c c2

pressure-time taken giving

response,

the integral

~co[-ikx~(kx,ky)]~(kx,Ky,M)ei

0 -CO

over all

I kxUt+p(MX-u)

k,

1dk,(7)

where

The factor exp[iu(ti-o)] gives only a time shift At = (o-Mx)/~~c to account for the time for the sound to propagate from the airfoil leading'edge to the far field and can be ignored. Thus,

P(t)=-7r

bz poU2 coc2

ja[-ikx~(kx,Ky)]~(kx,Ky,M)eikxutdkx -a-J

12

(8)

The Airfoil

Lift

Response

Graham (ref. 9) has shown that for a flat plate infinite span airfoil in subsonic flow the response function to a skewed gust can be obtained by similarity considerations from either (1) the gust response of the.airfoil to a parallel gust in compressible flow or (2) the gust response of the airfoil to a skewed gust in incompressible flow. The first or second solution should be used depending on whether the intersection point of the gust wavefront with the airfoil leading edge moves supersonically or subsonically relative to the fluid, respectively. For the present case it always moves supersonically a subsonic gust leads to an airfoil. The angle CLof the is

can readily be shown that the intersection point for the gusts of interest. Basically this is because exponentially decaying pressure field away from the gust wavefront relative to the airfoil leading edge

tan a= KY/k,= The intersection

point

sweeps along

the airfoil

My/u with

(9) velocity

V=U Cot Q = Cody Relative

Thus, flow,

to the fluid

the velocity

is

v, = Jm

= coJx2+y2

+pzqy

the solution requires only the airfoil response along with the similarity rules of Graham.

The similarity response functions function g(S, k,, tion is

(10)

2 co to a parallel

(11) gust

in subsonic

9) relates the surface pressure gust rule of Graham (ref. g(E,, k,, ky, M) for a skewed gust in compressible flow to the 0, M) for a parallel gust in a flow of Mach number M. The relag( in Equation 12 gives

(14) Together

with

Equation

6 this

gives

l(kx,Ky,M)=-!-

for P

S(~x~)ei~“X”Y-‘[Jo(~x/o)-iJ,(px/o)] P

(15)

For high frequency the airfoil response can be calculated Case 2: u, 1 0.4 using a method of Landahl (ref. 19) which treats the airfoil as a semiinfinite flat plate with a leading edge but no trailing edge. A correction can be added which treats the airfoil as a semiinfinite flat plate with a trailing edge but no leading An iteration can be carried out between the leading and trailing edges and edge. was shown by Landahl to converge for all frequency. The first two terms of this solution are (refs. 7, 11, 12) g,(x,kX,O,W=

I 7&(l+M)k,(l+X)

14

-I e .[~(I-M)(I+T)+rr/4-Tx]

s

(1W

g,(x,k,,

I

o, M) 2

[( 1+ i) E* (2p(

e’i[~(l-M)(l+~)+n/4-~~]

I-~))-I]

7d27r(l+M)Tix

(Mb)

where

(17) E* is a combination for a skewed gust

of Fresnel

integrals.

Substitution

I

g,(x,kx.Ky,W=

e-

in Equation

12 gives

M~(I+X)-M~~-I I -iv/4 iLx31E 1

&~,(l+MOD)(l+~)

‘J2(x,kx,Ky,f”‘) =

(18a)

I

e

-ii, * [M,(I+E)-M2X-l,-i.rr/4 1

d27&,(l+McD) (18b)

I(I+i)E*(2ix*Mq)(l-jt))-l] From Equation

6 Fis

found

to be

&kx,Ky,W

=

E’(28,)eiQ2

+/?

TTJi;,(l+M&

(Da>

eie2 P2(kx,Ky,W

=

7~,~27r~~(I+M~)

1

i(l-eBi2’I)

(1%)

fie-i20,

+(l-i)[E”(4Ex”M00)-.

E* (2$ F+

+I)]}

JI + (M/M co)(x/u)

where

8,=

ii,+ji(M-X/a)-

15

7~/4

(20b)

The above equations for whereas the computer program o 3 x BPF), tone noise spectra in Figure 35 exhibit A further a directivity pattern which is nearly symmetric about the rotor plane. evaluation of the vortex interaction noise directivity pattern is provided by Figure 36. Here measured tone levels are plotted as a function of B. for increasing

48

Loading noise dominated BPF, 2 x BPF, and 3 x BPF. blade passing frequencies. the directivity associated with the vortex interaction mechanism Consequently, could not be investigated at these frequencies. However, at higher blade passing frequencies experimental evaluation was possible. In this case, the measured radiation patterns were found to be approximately symmetric about the plane of the rotor blade which corresponded to '?, = 80.5'. Predicted directivity patterns are also shown in Figure 36 for selected blade passing frequencies dominated by blade vortex interaction noise (5 x BPF, 8 x BPF, 11 x BPF). Good agreement between theory and experiment exists for the directivity pattern shape. the measured and predicted minima occur In addition, near the plane of the blade. Since the acoustic source dipole axis is aligned normal to the blade planform in the present analysis, the source radiation is zero in the plane of the blade as verified by the experiment. Although generally good agreement between experiment and theory was observed for directivity pattern shape large differences occurred in the absolute tone amplitude. These differences were due to a general decrease in predicted tone amplitude with increasing frequency for all radiation angles, e,. This observation is based on a comparison of measured and predicted acoustic spectra which is discussed later in the subsection titled Assessment of Acoustic Spectrum and Pressure Signature. It should be noted that the value of B. used for the input to the analytical prediction was increased by the value of ~1. This is due to the difference in the reference axes selected for the theoretical and experimental definition of B. (see subsection titled Acoustic Measurements). According to Figure 4, B. in the analysis is measured from the normal to the blade planform. In contrast, the experimental program referenced 8, to the rotor centerline. At the blade vortex intersection station, these angles differed by the blade pitch angle setting. Comparison between theory and experiment, therefore, required increasing B. by This approach is used consistently throughout the al = 9.5" in the predictions. report. Directional Radiation Characteristics ___-- of Blade Vortex Interaction Impulsive Signature - The previous discussion of directivity pattern focused on tone amplitude Additional directivity changes occurring over the range of B. = 60" to 140'. features of blade vortex interaction noise can be identified from acoustic pressure Figure 37 shows ensemble averaged waveforms measured at selected signatures. The data corresponds to the same radiation angles 8, = 60°, 90", 120" and 140'. test condition described in Figure 36. Waveforms in Figure 37 contain the combined loading and blade vortex interaction Dominance of loading noise is evident in the sine wave character of the signatures. In contrast the impulsive signature peaks time history at all radiation angles. at angles outside the plane of rotation but diminishes to zero at 8, = 90" which is close to the plane of the rotor blade. These time history traces provide additional explanations for the directivity features determined from acoustic spectra.

49

Blade vortex interaction acoustic signatures are further isolated in Figure 38. Here measured waveforms are shown for the same microphone angles shown in Figure 37. Loading noise has been removed using the procedure described in Appendix B. Several waveform features are noteworthy. First, waveform shape is observed to be antisymmetric about the plane of rotation. Second, peak amplitude decreases near the plane of rotation. Finally, absolute waveform amplitudes are equal on opposite sides of the rotor. These important characteristics confirm the dipole acoustic source model developed in the present study.

Assessment

of Acoustic

Spectrum

and Pressure

Signature

Objective - Previous discussions focused on the experimental assessment of the sensitivity of blade vortex interaction noise to isolated changes of each parameter in the acoustic theory. Conclusions obtained from this assessment confirmed that the theory models the parametric dependence of each individual parameter. It now remains to determine if the analysis predicts the measured acoustic spectrum shape and absolute sound pressure levels. In addition, the analysis will be tested to determine if the measured acoustic pressure signatures can be predicted. Acoustic Spectrum - Figure 39 shows measured and predicted spectra for Test Condition B. Before comparing these spectra, the presence of steady loading noise must be discussed since this mechanism dominated frequencies given by BPF, 2 x BPF, noise source was removed from the measured spectrum and 3 x BPF. This additional by first using the waveform subtraction technique described in Appendix B.l. A Fast Fourier Transform (FFT) of the resulting isolated blade vortex interaction time history (consisting of 13 consecutive vortex interactions in the continuous time waveform) provided the acoustic spectrum shown in Figure 39. This spectrum coincides with the spectrum in Figure 25b except for the BPF tone and the first two harmonics. The general conclusion from the spectrum comparison in Figure 39 is that tone amplitudes are predicted to within 5 dB at low frequencies. At high frequencies, the present analysis significantly underpredicts the data. acoustic signature for Acoustic Signatures - Figure 40a shows the predicted the experimental test case presented in Figure 38~. Both waveforms are plotted on A comparison of the same scales and correspond to the same measurement angles. the measured and predicted waveforms indicates that the duration of the acoustic In addition, signature is similar as is the sign of the maximum peak pressure. the experimentally observed antisymmetric waveform occurring on opposite sides of Calculated waveforms were found the rotor are predicted by the present analysis. to be identical (not shown) at e. = 50" and 120", except for the sign reversal in the acoustic pressure. (Note that these angles are symmetric around the plane of The qualitative agreement between the rotor blade in the present experiment.) theory and experiment again conforms the dipole acoustic source model.

50

While the general features of the blade vortex interaction signature are predicted, details cannot presently be calculated. For example, the peak amplitude, triangular shape, and combined positive and negative peaks of Figure 38c are not predicted. Assessment of Difference -.-..- - - i . ..-.___. _._ __ _Between Theory and Experiment - The rotor blade interaction with the two dimensional wake of the upstream airfoil was initially considered to be a possible source of the difference between theory and experiment. This mechanism was, however, evaluated experimentally during the blade vortex intersection angle tests. The results indicated that when 81 = 0, the axial velocity defects did not significantly contribute to the tone noise generation. This conclusion confirmed that the vortex azimuthal velocity component dominated the acoustic response during the present study. A second possible source of the difference between theory and experiment is the infinite span airfoil assumption used in the analysis. With this assumption, physical constraints imposed by the finite span rotor blade are not satisfied by the theory. To understand the consequences of the infinite span airfoil model, recall that the far field sound is related to the time derivative of the airfoil loading. Calculating the area under the far field pressure waveform is, therefore, equivalent to a time integral of the pressure signature in Figure 40a. The result of the integral is the difference between the lift when the airfoil is at an infinite distance upstream and downstream of the airfoil

AREA = /

Ldt = Lco-L --m

(41)

For the present infinite span airfoil model and the l/r decay of the azimuthal velocity vector, Equation 42 gives a finite or non-zero difference between L, This occurs because the calculated pressure pulse in Figure 40 is always and L-,. Pressure peaks of opposite sign, as in the experiment, would less than zero. cause the area in Equation 41 to approach zero. The non-zero value obtained from Equation 42 represents a theory limitation since a finite span airfoil situated at an infinite distance from the vortex has To approximate the finite span blade a zero value for the quantity L, - L-,. Equztion 22 was multiplied by an exponential tested in the present experiment, the radius at which the decay function defined as zfp[-(R,/r) 1. Here R, represents field due to exponential amplitude is e . The rapid decay of the vortex velocity the additional exponential decay function ensures that the integral in Equation 42 Introducing the exponential decay function represents more than has a zero value.

51

It is also a mathematical termination of the influence of the vortex velocity field. The acoustic signature is the an approximate simulation of a finite span airfoil. same whether the vortex velocity field is terminated or the airfoil is terminated This equivalence, however, assumes that the blade as in the case of a finite span. end effects can be ignored. Using the exponential decay function to model the finite span blade results of the features of the measured predicted acoustic waveforms which display some This is illustrated in Figure 40b where predicted pressure signals waveforms. to the blade chord. are presented for a range of R,/c values where c corresponds I, in the present experiment, from the Here R,/c = 1 approximates the distance, Terminating the analytical infinite vortex interaction station to the blade tip. pressure signal which span airfoil model at this distance generated an acoustic exhibits combined positive and negative peaks similar to the experiment. The R,/c = 1 curve displays an initial positive pressure pulse, a negative peak, to zero as time increases. and a small positive pressure waveform which returns to return to the one-sided waveform Increasing R,/c causes the acoustic pressure .predicted in Figure 40a.

in

The above discussion illustrated how the waveform shape depends on the infinite span airfoil assumption. This assumption can also result in incorrect predictions of the amplitude of the acoustic signature as illustrated by the following example. Consider a plane normal to the span line of the airfoil. Also let the normal to the airfoil planform and the vorticity vector be in this normal plane (4, = 0" or 180" and 6v arbitrary). Then the unsteady transverse flow sensed by the airfoil will be antisymmetric about this normal plane; i.e., if the flow is upward on one side of this normal plane it is downward on the opposite side with an equal magnitude. An observer stationed anywhere in the normal plane will measure zero sound since a positive pressure perturbation produced on one side of the normal plane is cancelled by the opposing negative pressure. If the infinite span airfoil in the above description is now truncated on one side of the normal plane (as in the case of the finite span rotor blade used in the present experiment) the opposing unsteady pressure responses no longer occur and the above described acoustic cancellation disappears. This finite span effect is not accounted for by the previously described experimental decay factor since the factor retains the symmetry of pressure response on the airfoil surface. Thus, it is possible to also achieve an increase in predicted sound level when using a truncated (finite span) airfoil. The asymmetry associated with the finite span for the discrepancy between predicted and measured in addition to acoustic signatures (Fig. 40). This microphone measurement plane almost coincided with plane. Also, for the small values of Bl tested in vorticity vector was almost in the normal plane.

52

provides a potential explanation acoustic spectra (Fig. 39) follows since the present the above described normal the present experiment, the

CONCLUSIONS

1. Ingestion of a vortex by a rotor generates harmonic noise which extends to high frequencies. Acoustic pressure signatures associated with this noise mechanism display impulsive waveforms superimposed on the blade loading acoustic pressure signals. Linear dependence of measured blade vortex interaction noise on blade number 2. and vortex strength was closely predicted by the present first principles acoustic theory. The absence of noise for vortex trajectories normal to the blade was also predicted. For oblique intersections with a blade, calculated changes in acoustic radiation due to intersection angle modifications were within 2 decibels of the measured data. A simplified scaling law was determined for estimating these changes at small intersection angles. The predicted low sensitivity to blade pitch angle modifications for a typical 3. range of rotor operating conditions was cnnfirmed experimentally. The strong noise dependence on local Mach number was calculated although the sensitivity to Mach number change was overpredicted at low frequencies. Measured blade vortex interaction directivity 4. dipole acoustic source model at all frequencies. observed experimentally in the plane of the blade

patterns confirmed the theoretical The minimum sound pressure level was predicted analytically.

Absolute levels for acoustic spectra and directivity patterns were predicted, 5. without the use of empirical or adjustable constants, to within 5 decibels at low frequencies. sound pressure levels are currently underpredicted. At high frequencies, General features of the measured blade vortex interaction signature were pre6. dicted. Antisymmetry of the measured waveform on opposite sides of the rotor plane confirmed the theoretical dipole directivity model. Duration of the impulsive acoustic signature in addition to the sign of the maximum pressure peak was predicted. Detailed features of the measured acoustic signature, such as peak amplitude, triangular shape, and combined positive and negative peaks, are not presently calculated.

53

APPENDIX A

Vortex

Velocity

Field

Measurement

Technique

Decomposition of Total Velocity Vector - The azimuthal (V,) and axial (U) velocity components of the vortex velocity field are shown in Figure 9. Decomposition of the measured voltage from the single slanted hot film into V0 and U employs the following procedure. it is necessary to recognize that the probe sensing element can be roFirst, tated until it lies in the plane formed by the components Ve and U. (Note that the vortex radial velocity is assumed to be zero in this velocity decomposition Here the Ve and This unique orientation is illustrated in Figure A-l. approach.) U components lie in the plane formed by the hot film sensor and the support prongs. In Figure 9 this corresponds to alignment of the probe sensing plane with a plane In this case, the probe sensing plane is normal to a parallel to the z, y axes. radial line penetrating the vortex center. To decompose the total velocity vector, UT, in Figure A-l into the corresponding values of V, requires two separate measurements. These correspond to two orientations which are 180" apart in Figure A-l. The velocity components are then given by the equations: = UT cos (BF-al)

%

u=

UT sin

(BF-al)

Representing the two separate tion as UN1 and UN2 gives

Since the sensor linearly related

= UT (cos B cos al + sin F = UT (sin

normal

sin

al)

(A-1)

BF cos al - cos BF sin

al)

(A-2)

velocity

components

6

F

at each wire

orienta-

= UT cos al U Nl

(A-3)

U = UT sin N2

(A-4)

al

mean voltage output El and E2 at each of the two orientations to the normal velocity component

is

UNl = KE1 = KE2 U N2 where K represents

the velocity

versus

voltage

54

(A-6)

calibration

coefficient.

Using.Equations A-5 and A-6 in the expressions for Ve and U gives the sum and difference expressions for velocities measured using a standard crosswire %

= 0.707 K(El-E2)

u=

0.707 K(El+E2)

Thus, a measurement at two orientations in a plane trating the vortex center gives Ve and U.

normal

to a radial

line

pene-

Instead of rotating the sensor at each radial velocity station to obtain V, A continuous traverse was conducted on and U the following procedure was used. An analog a radial line through vortex center for each of the two orientations. output from the traverse system was used to drive one axis of a two directional plotter. Two continuous traces of the voltages El and E2 were then obtained as the Graphically summing and differencing the plotter sensor penetrated the vortex. traces decomposed the measured voltages into El-E2 and El+E2 at each radial staA simple application of the calibration factor K then permitted changing tion. the plotter axes notation to Ve and U as shown in Figure 17. radial velocity Identification of the Vortex Center - The above described The followtraverse assumes a knowledge of the line penetrating the vortex center. ing simple procedure was developed to identify the vortex center. there exists an alignment of the As described in the previous subsection, and U are easily determined. This plane is sensor measuring plane for which Ve The plane can be identified normal to the radial line through the vortex center. When the sensor voltage output by rotating the sensor 360" about the probe shaft. components lie in the plane of the sensor. This reaches a maximum, all velocity orientation is designated in Figure A-l as the maximum voltage output orientation. Other orientations can be shown to give smaller voltage values with a minimum occurring at 180" from the maximum orientation. The vortex center then lies on a line perpendicular to the unique measurement Rotating the sensor at numerous stations on the perimeter of the vortex plane. provided a series of radial lines whose intersection point defines the vortex center. Figure A-2 provides a schematic description of this experimental procedure. A digital angular position encoder with a 0.33" resolution was used to sense the probe orientation.

55

APPENDIX B Experimental Technique for Isolating BladexVortex Interaction Impulsive Signatures

Measured acoustic pressure signatures contained both loading noise and the The following procedure was employed blade-vortex interaction impulsive signature. to isolate the impulsive waveform permitting direct comparison of measured and predicted blade-vortex interaction pressure signatures. The experimental technique begins with the ensemble averaged waveform for the A sample time history is shown in combined loading and vortex interaction noise. Figure B-la where the time marker signal from the rotor rig is used to initiate This waveform is obtained with the repetitive sampling of the acoustic waveform. upstream semispan airfoil at a finite angle of attack (~2 = 12" in the present example). A second waveform is obtained with the airfoil at c2 = 0" (Fig. B-lb). Under the vortex strength is forced to zero while the loading noise this condition, the axial remains the same. At the same time, noise generated by intersecting velocity defect in the vortex core and in the two dimensional wake are-considered between Figures B-la and B-lb represents the to remain the same. The difference isolated blade vortex interaction signature which is shown in Figure B-ld. An The above described enlarged version of this waveform is presented in Figure B-le. ensemble averaged pressure signatures and the arithmetic subtraction of these waveforms were obtained on a dual channel Fast Fourier Transform spectrum analyzer with time domain averaging capability. The rationale for using the ct2 for the axial velocity defect which The noise generated by this velocity Using the observed in Figure B-lb. upstream airfoil removed (Fig. B-lc) mechanism.

= 0" waveform is based on the need to account is also present in the waveform of Figure B-la. defect creates the small pressure pulse loading noise pressure signature with the would fail to account for this noise

56

REFERENCES

1.

Leverton, Interaction.

J. Tlr., J. C. Pollard, and C. R. Wills: Main Rotor Vertica, Vol. 1, No. 3, pp. 213-221, 1977.

Wake/Tail

Rotor

2.

Pegg, R. J., and P. A. Shidler: Exploratory Effect of the Main Rotor Wake on Tail Rotor NASA CP-2052, Part 1, pp. 205-219, 1978.

3.

Paterson, Ingestion

4.

Amiet, R. K.: Aerodynamic Sound Production and the Method of Matched Asymptotic Expansions. Ph.D. Thesis, Cornell University, 1969.

5.

Widnall, S.: J. Acoustical

6.

Widnall, S. E. and T. L. Wolf: Effect Noise Due to Blade-Vortex Interaction. 1980.

7.

The Passage of an Infinite Adamczyk, J. J.: Oblique Gust. NASA CR-2395, 1974.

8.

von Karman, T. and W. R. Sears: the Aeronautical Sciences, Vol.

9.

Similarity Graham, J. M. R.: J. Fluid Mech., Vol. Flows.

10.

Amiet, R. K.: AIAA Journal,

11.

Effects of Compressibility in Unsteady Airfoil Lift Theories. Amiet, R. K.: Unsteady Aerodynamics, University of Arizona/AFOSR Symp., Tucson (Ed., R. B. Kinney), pp. 631-653, 1975.

12.

Amiet, R. K.: Journal, Vol.

13.

Low Frequency Approximations in Unsteady Small Amiet, R. K.: J. Fluid Met., Vol. 75, pp. 545-552, 1976. Subsonic Flows.

Wind Tunnel Investigation of the Noise. Helicopter Acoustics,

R. W., R. K. Amiet: Noise of a Model Helicopter of Turbulence. NASA CR-3213, November 1979.

Rotor

Due to

Helicopter Noise Due to Blade-Vortex Interaction. Society of America, Vol. 50, pp. 354-365, 1971.

Airfoil 8, pp.

of Tip Vortex J. Aircraft,

Swept Airfoil

Theory 104-108,

Compressibility Effects in Unsteady Vol. 12, pp. 253-255, 1974.

57

Through

for Nonuniform 1940.

Rules for Thin Airfoils 43, pp. 753-766, 1970.

High Frequency Thin-Airfoil 14, pp. 1076-1082, 1976.

Structure on Helicopter Vol. 17, pp. 705-711,

Theory

an

Motion.

J. of

in Non-Stationary

Thin-Airfoil

for. Subsonic

Theory.

Flows.

Perturbation

AIAA

REFERENCES (Cont'd)

14.

Amiet, R. K.: Airfoil Spanwise Wavenumber.

15.

Graham, J. M. R.: Lifting Surface Theory for the Problem of an Arbitrarily Yawed Sinusoidal Gust Incident on a Thin Aoerfoil in Incompressible-Flow. Aero. Quart., Vol. 21, pp. 182-198, 1970.

16.

Acoustic Ray Paths through a Model Vortex Georges, T. M.: Core. J. Acoust. Sot. Amer., Vol. 51, pp. 206-209, 1972.

17.

Larson, R. S. and K. W. Robbins: Refraction Vortices. AIAA Paper #80-0975, 1980.

18.

Morse, P. M. and K. U. Ingard: Theoretical New York, 1968, Equation 11.2.33.

19.

Landahl,

20.

Gradshteyn, I. S., and I. Academic Press, New York,

21.

Abramowitz, Publications,

22.

Newland, D. E.: Longman, London,

23.

Cooley, and its 1969.

24.

Paterson, R. W., P. G. Vogt, and W. M. Foley: United Aircraft Research Laboratories Acoustic Vol. 10, No. 7, pp. 427-433, 1973.

25.

Hama, F. R.: An Efficient Tripping Vol. 24, pp. 236-237, March 1957.

26.

Schlinker, R. H., and R. K. Amiet: Shear Layer Sound. AIAA Paper No. 80-0973, 1980.

M.:

Response to an Incompressible Skewed Gust of Small AIAA Journal, Vol. 14, pp. 541-542, 1976.

Unstead.y Transonic

Flow.

M. Ryzhik: 1965.

Acoustics.

New York,

of Integrals

to Random Vibrations

Wing Tip

McGraw Hill

M. and I. A. Stegun: Handbook of Mathematical Inc., New York, 1968. An Introduction 1975.

a Viscous

of Sound by Aircraft

Pergamon Press, Tables

with

Series

Book Co.,

1961. and Products.

Functions.

and Spectral

Dover

Analysis,

J. W., P. A. W. Lewis, and P. D. Welch: The Fast Fourier Transofrm Applications. IEEE Transactions on Eductions, Vol. 12, pp. 27-34,

Device.

58

Design and Development of the Research Tunnel. J. of Aircraft,

J. of Aeronautical

Refraction

Sciences,

and Scattering

of

REFERENCES (Cont'd)

L.:

Attaining

a Steady Air

Stream in Wind Tunnels.

NACA TM 726,

27.

Prandtl, 1933.

28.

Batchelor, G. K.: Press, 1976.

29.

Spectrum of Rotor Noise Caused by Atmospheric Turbulence. Hanson, D. B.: J. of the Acoustical Society of America, Vol. 56, No. 1, pp. 110-126, 1974.

An Introduction

to Fluid

59

Dynamics.

Cambridge

University

g

12 12 12

6"

9.5

14"

9.5

9.5

B

C

D

E

3.2"

12 12"

9.5

9.5

I

J

no intersection

5

12

4.4"

9.5

12"

7"

G

0.

4

4

4

2

0.39

0.35

0.39

0.55

0.55

0.55

4

14.6"

2

0.55

2

11.9"

0.39

2

6.8"

0.55

0.55

4 4

-MT

B -

14.6

-B1

H

12"

5"

F

no airfoil

9.5”

A

a2

-5

Test Condition

56

50

56

79

79

79

79

56

79

79

Rps

EXPERIVENTAL TEST CONDITIONS

TABLE I

0.184

0.166

0.173

0.166

0.173

0.186

0.171

0.173

I/R1

-1.1

0.440

0.361

0.388

0.550

0.

0.055

0.277

0.

H/c 2

0.0193

0.0192

0.0195

0.0197

0.00897

0.0185

0.0192

0.0180

Kl/cUo

0.043

0.043

0.043

0.043

0.040

0.044

0.044

0.044

r 0 /c 2

STRONG MAIN ROTOR BLADE VORTEX INTERACTION

FLIGHT DIRECTION

WEAK MAIN ROTOR BLADE VORTEX INTERACTION

ADVANCING

RETREATING SIDE

-

a) ILLUSTRATION OF MAIN ROTOR AND TAIL ROTOR BLADE VORTEX INTERACTION VORTICITY VECTOR NORMAL TO PAGE VORTEX

b) ILLUSTRATION OF MAIN ROTOR BLADE VORTEX INTERACTION WITH VORTEX PARALLEL TO LEADING EDGE

Figure 1 -

Helicopter

Blade Vortex Interactions

61

Figure 2 -

Open Jet Test Section Arrangement

62

in Anechoic Chamber

INLET NOZZL

Figure 3 -

Upstream Model Airfoil and Rotor Blade with Boundary Layer Trip

63

cn .P

Figure 4 - AirfoiLVortex

Geometry

-b CHORD STATION

ZA

ORSERVER

LINED

URBULENCE

ACOUSTIC

CORNERS

ACOUSTIC

WEDG

CONTRACTION COLLECTOR

AIR FLOW +

\I L

HIGH

PRESSURE

AIR

SUPPLY

Es TT1

17m

--EXHAUST TOWER I

SYSTEM

TOP VIEW

CENTRIFUGAL

FAN -..A

SIDE VIEW

Figure 5

-

UTRC Acoustic

65

Research Tunnel

DRIVE

MOTOR

Figure 6 - Acoustic Research Tunnel Control Room with Data Acquisition and Data Reduction Instrumentation

66

I

PLAN VIEW

VORTEX

a) PLAN VIEW WITH UPSTREAM AIRFOIL AND ROTOR BLADES ROTOR BLADE 3 OPEN JET NOZZLE -uo

ri. -iI

r( UPSTREAM AIRFOIL, 11.4

cm CHORD

1

58 cm-

SIDE ELFVATlON

1

I I I

HOLDER FOR VERTICAL ADJUSTMENT

I

II:

FLOOR (COVERED WITH WEDGES DURING ACOUSTIC STUDY)

I //////////////////////////////////////////

b) SIDE VIEW SHOWING OPEN JET TEST FACILITY

Figure 7 -

Detailed Schematic of Vortex Generator and Model Rotor

67

a) PLAN VIEW

60”

f x

Microphone

\

TO MICROPHONE AT/~~=130”.d~=135”

b) VIEW LOOKING UPSTREAM

MICROPHONE

Arrangement

FAR FIELD MICROPHONE STATIONS

Figure 8 -

-4.

I/ n /

I-

OPEN JET INLET NOZZLE

MAXIMUM AZIMUTHAL VELOCITY, v. AZIMUTHAL

AXIAL VELOCITY IN z, x PLANE

UPSTREAM AIRFOIL ANGLE

CENTERLINE

AI

J NOTE 1) VORTEX CENTERLINE, z AXIS, PROBE CENTERLINE AND OPEN JET CENTERLINE ARE ALL CONSIDERED PARALLEL 2) FIGURE ASSUMES DOWNWASH DUE TO AIROFIL NON-EXISTENT

Figure 9 -

Definition of Vortex Velocity Field

69

Figure 10 - Tip Vortex Generating Airfoil and Probe Traverse System in Acoustic Tunnel

Figure

I I

-

Smoke Flow Visualization Showing Tip Vortex for Two Operating Conditions

71

Figure 12 -

Flow Visualization

of Blade Vortex Interaction

72

for Lifting Rotor

TAIL ROTOR CENTERLINE

VORTEX ORIENTATION AFTER SKEWING OF FLOW BY TAIL ROTOR INFLOW

P

t

a) MAIN ROTOR VORTEX TRAJECTORY RELATIVE TO TAIL

x2 TAIL ROTOR

RQTOF( AFTER APPLYING GALILEIAN TRANSFORMATION

VORTEX TRAJECTORY AFTER CONTRACTION

VORTEX TRAJECTORY IN ABSENCE OF CONTRACTION

FIXED AIRFOIL MEAN FLOW NEEDED TO GENERATE VORTEX

b) VORTEX TRAJECTORY IN OPEN JET WIND TUNNEL WITH ROTOR

Figure 13 - Equivalence of Tail Rotor Vortex Interaction Process for Full Scale Helicopter and Open Jet Acoustic Wind Tunnel Model Rotor Experiment

73

20

2 cfY 30

2 4( m -cl z-

0

Figure 14 -

1

3

4

5 FREQUENCY, KHz a) MT = 0.35

6

180”

7

8

Background Noise with Open Jet Operating and Rotor Rig Operating Without Hub or Blades

2

/lo= 120”. $,=

NO HUB OR BLADES

TEST CONDITION: NO AIRFOIL

9

10

9

0

60 --

70

1

2

3

5

Concluded

b) MT = 0.55

FREQUENCY. KHz

Figure 14 -

4

6

7

NO AIRFOIL NO HUB OR BLADES f~o=1200, $,=180”

TEST CONDITION:

8

9

10

60

80

0

Figure 15 -

0.1

0.3

0.5

0.6

0.7

0.8

0.9

Blade Tracking Accuracy and Absence of

BLADES

FREQUENCY, KHz a) B = 2

0.4

Acoustic Spectrum Demonstrating Freestream Turbulence

0.2

n,=90",~,=180"

cY,=14o

MT = 0.39

RPS = 56

B=2

NO AIRFOIL

TEST CONDITION:

1.0

70

80

0

0.1

0.2

FREQUENCY, KHz

0.5

Concluded

b) B = 4 BLADES

0.4

Figure 15 -

0.3

$,=90”.

$+,=180”

0.6

n, = 14”

MT = 0.39

RPS = 56

NO AIRFOIL B=4

0.7

TEST CONDITION: IDENTICAL

0.8

0.9

1 .o

60

80

0

Figure 16 -

0.2 0.6

FREQUENCY, KHz

1 .o

1.2

1.4

1.6

Blade Tracking Accuracy and Absence of

a) FREQUENCY RANGE o-2 KHz

0.8

Acoustic Spectrum Demonstrating Freestream Turbulence

0.4

TEST CONDITION NO AIRFOIL B=4 ‘Y, = 9.5” n,=120”, c#1~=180’ RPS = 79 MT = 0.55

1 .a

2.0

40

90

0 1 2 3

FREQUENCY, KHz

5

IDENTICAL

6

Figure 16 -

Concluded

b) FREQUENCY RANGE O-10 KHz

4

TEST CONDITION -

7

8

9

10

-0.5

-0.4

-0.3

-0.2

-0.1

0

0.1

0.2

0.3

0.4

0.5

-0.75

Figure 17 -

-1 .oo

f-

0

r

0 .25

VORTEX CENTER

_

0

a) H/c2 = 0.88, -1 I r/c2 I 1, AR = 3.1, CY~= 12’

0.5

ANALYTICAL MODEL

v,J% MEASURED

AIRFOIL PRESSURE SIDE

NORMALIZED RADIAL DISTANCE, r/c2

-0.25

0

--t

0.75

Azimuthal Velocity Dependence on Radial Distance from Vortex Center

-0.50

--I Pcm

AIRFOIL SUCTION SIDE

1 .oo

Z

8

2

_

-

-0.5 -0.25

-0.4

-0.3

-0.2

-0.1

0

5 cl :

i

0.1

0.2I

3” -k ? c G 2

0.3,

0.4I-

0.5,

0.Ei-

-0.20

0

O

}-

-

1 cm

-0.05

0 0 0

0

0

0

-

0

0

0

0

0

0 0

0.05

O&O

0 ’

c

Figure 17 -Continued

AR=3.1,

0.10

0

~~~=12’

NORMALIZED RADIAL DISTANCE, r/c2

-0.10

-1

o-o

D/c, = 0.088 --I

b) H/c2 =0.88, -0.251 rlc2s0.25,

-0.15

0

0

0

0

l---

0

0.15

0

0

0.20

0

0

0.25

-1 .oo

-0.5 1

-0.4

-0.3

-0.2

-0.1

0

0.1

0.2

0.3

0.4

0.5

0

I -0.75

0

0

o

0

0

0

I

0

0 0

0

0

0.25

I

0

VORTEX CENTER

Figure 17 -

Continued

c) H/c2 = -0.55, -1 I r/c2 5 1, AR = 1.6, cyy2= 12”

NORMALIZED RADIAL DISTANCE, r/c2

I

-0.25

I

-0.50

--I1 cmt-

O

0

0

0

0

000

0.50

I

0

I

0.75

O

0

0

1

1.00

0.1

0

-0.1

-0.2

c 0 0

2 n I4 2

z g

L

-0.25

-0.5

-0. ’

-0.3

0.2

?I .>

0.3

0.4

0.5

-0.20

I

0

O

I

1 cm

0

-0.05

I

0 0 0 0 0

0

0

I

0

0 0 0

I

0.05

Oo&

--I

O

Figure 17 -

Concluded

ty2

0

0.10

I

= 12’

NORMALIZED RADIAL DISTANCE, r/c2

-0.10

I

-4

0

OCWD

D/c2 = 0.091

d) H/c2 = -0.55, -0.251 r/c2 I 0.25, AR = 1.6,

-0.15

I-

0

0

I-

0

0.15

I

0

0

0.20

I

0

0.25

0

I

III

n

VORTEX CENTER

”/

NORMALIZED

Figure 16 -

NOTE - VORl ‘EX CENTEfi SHIFTED BETWEEN TRACES TO ACCOMODATE ALL TRACES ON SAME PAGE

RADIAL DISTANCE, r/c2

Radial Traverses Showing Vortex Symmetry, H/c2 =0, AR =2.22, ~~‘12”

84

0.14 0” z

ff2, deg

0.12

0

12

0

6

5 kl 2 n cl ;

A v KO