25.07.2013 Views

Second Computational Aeroacoustics (CAA) Workshop on ...

Second Computational Aeroacoustics (CAA) Workshop on ...

Second Computational Aeroacoustics (CAA) Workshop on ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Herec is speedof sound,and 1)themeanflow velocitywhichis assumedto have<strong>on</strong>lytheaxialcomp<strong>on</strong>ent.<br />

Theanalysisis c<strong>on</strong>finedto a steadywavewith theharm<strong>on</strong>ictimedependence -_o,,, and axial angle<br />

dependence e '" _, where m is an integer called the circumferential mode number. Eq. (1)_is then written for<br />

circular cylindrical coordinates (r, _b, x) as<br />

where k = co/c, and M = V/c.<br />

-(<br />

r grt, G-r)+ax -7-p=°<br />

The sound radiati<strong>on</strong> from an unflanged circular duct is schematically displayed in Fig. 1. With reference to this<br />

figure, the entire regi<strong>on</strong> is divided into two: regi<strong>on</strong> 1 for r < a, and regi<strong>on</strong> 2 for r > a, where a is the duct<br />

radius. The subscripts 1 and 2 will be used from now <strong>on</strong> to indicate respectively the regi<strong>on</strong> 1 and 2, unless<br />

specified otherwise. The mean flow velocities are assumed to be uniform in each regi<strong>on</strong>, and denoted by Vj<br />

and Vz. For V_4: V2, there will be the mean flow mismatch at r = a, for x > 0. The sound speed can differ for<br />

the two regi<strong>on</strong>s for reas<strong>on</strong>s such as differences in mean density and temperature. The respective sound<br />

speeds, wave c<strong>on</strong>stants, air densities, and Mach numbers are denoted by q and c2, k_ and k2, Pl and P2,<br />

and M I and M 2, where k t = co / q, k 2 = co / c2, ,,141= Vt / q, and M 2 -- V2 / q.<br />

In a hard-wall circular duct, the general soluti<strong>on</strong> to Eq. (2) is obtained as<br />

p(r,x)= Jr. {A,..e_L'X +B,..e_'_} •<br />

n=l<br />

Here J., is the Bessel functi<strong>on</strong> of order m, /1,.. the n-th zero of J" (x), and A,.. and B,... c<strong>on</strong>stant<br />

coefficients. The wave c<strong>on</strong>stants k_. and k,7,. corresp<strong>on</strong>d to the mode propagati<strong>on</strong>s respectively in the<br />

positive ( to the right) and the negative (to the left) directi<strong>on</strong>s, and are given by<br />

-klM j q-_ kl 2-(1 - M 2)( l'lmn ) 2<br />

kmn ÷ "_l-m,<br />

2 a (4)<br />

The integer n here is called the radial mode number, and the pair (m,n) is used to represent a single duct mode.<br />

C<strong>on</strong>sider the incident wave of a single mode, say (m,g),<br />

Pinc--Jml--I e e .<br />

kay<br />

This wave is incident from x =- _ and propagating towards duct terminati<strong>on</strong> as in Fig. 1. Up<strong>on</strong> arriving at<br />

the duct terminati<strong>on</strong>, it will be partly reflected back into the duct, and partly radiated out of the duct. In<br />

general, the reflected wave c<strong>on</strong>tains many radial modes including propagating and attenuating modes, and is<br />

represented by<br />

28<br />

(2)<br />

(3)<br />

(5)<br />

M<br />

J

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!