Bilanin 1973
Bilanin 1973
The physical processes that sustain discrete-frequency oscillations of cavities adjacent to compressible flow
are modeled analytically, yielding a formula which predicts excitation frequencies as a function of Mach number
and cavity geometry. These physical processes are similar to those used by Powell in describing the mechanism
underlying the production of edge tones. The empirically determined constants appearing in Rossiter's formula for
excitation frequencies are computed from the model. It appears that instability of the shear layer as well as
Downloaded by YORK UNIVERSITY LIBRARIES on August 12, 2014 | http://arc.aiaa.org | DOI: 10.2514/3.6747
interaction between the shear layer and the cavities' trailing edge is required to sustain discrete-frequency
oscillation. It is suggested that the simultaneous excitation of two or more discrete frequencies (which are not
harmonic), as have been observed in practice, correspond to the simultaneous participation of two or more
vortex sheet displacement modes. The model analyzed yields qualitatively correct acoustic mode shapes in the
cavity. Theoretical results include the calculation of possible excitation frequencies over the range 0.8 ^ M ^ 3
for rectangular cavities and show their dependence on cavity depth. Analytic results are in general agreement
with experimental data.
Nomenclature = S/a + K
= correction factor for edge effects
A = magnitude of harmonic force acting on vortex sheet at = fluid density
\x\ ^ e, y = 0 = — ^ + S2, phase of H 0 (1) for large argument
a = sound speed = lag of vortex sheet displacement after forcing at x = 0
B = constant in mass-flux equation [Eq. (27)] = velocity potential
C(u) = Fresnal integral = circular excitation frequency
D = cavity depth = Laplace transform variable
/ = frequency = Fourier-Laplace transformed dependent variable
H(u) = Hankel function
i = (-l)" 2 Subscripts
K = complex wave number in x direction and Fourier transform
variable ( ) + = outside the cavity, y > 0
L = cavity length ( ) _ = inside the cavity, y < 0
M = u/a, Mach number ( ) r = real part of ( )
m = mass flux at cavity trailing edge (),- = imaginary part of ( )
n, m = summation indices in x, y direction, respectively
P = perturbation pressure
R = image-correction sum Introduction
r = thermal recovery factor
rj, 7
S
= source and observation point position vectors, respectively
= -i(tt-KU}
D ATA from many experiments conducted over the years have
shown that cut-outs and cavities exposed to a flowing
fluid sometimes produce intense, periodic, but not necessarily
51 = fL/U+, exterior Strouhal number
simple harmonic pressure fluctuations. Data from deep and
52 = coL/a_, interior Strouhal number
S(u) = Fresnal integral shallow cavities, from rectangular, triangular, and circular
t = time cavities, suggest that the unsteady motion may be treated as if it
U = convection velocity (note U = 0 inside the cavity) were two-dimensional (i.e., the oscillations are independent of the
x, y = coordinate system (see Fig. 1) transverse or spanwise direction) even though the steady motion
y = ratio of specific heats is three-dimensional.1"3 Clearly two-dimensional motion is
A = dispersion relation much simpler to analyze. These analyses do not provide a means
e = small physical length of estimating the amplitude of the pressure fluctuations and,
Y\ = shear-layer displacement except for Ref. 1, do not attempt to determine when the pressure
K = fractional vortex convection speed
in the cavity will oscillate. By studying possible modes one can
Aa = wavelength of acoustic wave
AB = wavelength of shed vortices estimate the frequency of the oscillation.
Any investigation of this problem must start from the idea
Received February 2, 1972; revision received September 15, 1972. of two flowfields. (See Fig. 1 for the geometry of the situation.)
This material is an abridgement of A. Bilanin's thesis (S.M.) submitted One is an outside flowfield, characterized by the forced con-
to the Department of Aeronautics and Astronautics at M.I.T. vection. The second is an inside flowfield, characterized by an
A. Bilanin would like to thank NASA and NSF for fellowship energy-storing mechanism. The energy-storing mechanism may
support during this investigation. The authors would like to thank be either acoustic (volume process) or elastic (surface process).
C. Tarn of the Department of Aeronautics and Astronautics, M.I.T. These two flowfields are separated by a mixing layer that is
for his suggestions.
Index categories: Aerodynamic and Powerplant Noise; Nonsteady usually thin. This mixing layer may be approximated by a
Aerodynamics. vortex sheet. Consequently, flow is assumed inviscid everywhere.
* NSF Fellow, Department of Aeronautics and Astronautics. In Ref. 1 the problem is considered in terms of hydrodynamic
t Professor, Department of Aeronautics (on leave); presently Chief stability. The vortex sheet is subjected to an infinitesimal dis-
Scientist, USAF. Associate Fellow AIAA. turbance and the temporal growth or decay of the perturbation
348 A. J. BILANIN AND E. E. COVERT AIAA JOURNAL
to the combined flowfield determined. Results of this analysis The results presented here are derived in the same pragmatic
suggest the classification of deep or shallow cavities depends spirit that characterizes Rossiter's model.! In fact Eq. (3) will be
upon the effect of the vertical or horizontal boundaries on the rederived and the constants will be determined from the calcula-
vortex induced velocity field. For deep cavities the vertical tion. The calculation is not based upon "first principles" but
walls seem to cause vertical motion to dominate. Asymptotic based upon an analytical model. Following Ref. 1, it is assumed
results for deep cavities in Ref. 1 support this. However, no that the motion of the vortex sheet is more strongly dependent
results are given there for shallow cavities. upon the bottom boundary than the ends. Hence this part of
Downloaded by YORK UNIVERSITY LIBRARIES on August 12, 2014 | http://arc.aiaa.org | DOI: 10.2514/3.6747
Shallow cavities have been studied at supersonic speeds by the solution will be based on the study of motion of an infinite
Karamcheti2 and at subsonic and transonic speeds by Rossiter.3 vortex sheet separating the inviscid convecting flow with its
In Ref. 3, Rossiter presents results of his study of shadow- velocity, U + , from a quiescent inviscid fluid field of depth D.
graphs. During oscillation, periodic density variations were This vortex sheet is excited at the origin x = 0, corresponding
observed across the mouth of the cavity and were assumed to to the leading edge of the cavity, by a periodic pressure pulse
be vortices shed from the leading edge of the cavity. The of length 2s in the x direction (Fig. 2a). The dominant pressure
shadowgraphs also show acoustic waves whose primary source fluctuations at the trailing edge of the cavity (x = L) are modeled
is near the trailing edge of the cavity. Rossiter used these by a single periodic acoustic monopole (Fig. 2b). This is a
observations to construct a model for the unsteady process, simplified representation of a much more complex phenomenon
which is based upon two conditions. The first condition is that that is illustrated in Fig. 3. The separation of leading and trailing
the frequency / of the vortex shedding process is equal to the edge effects together with the assumption that the pressure field
frequency of the acoustic process—that is, from the trailing edge has little effect on the vortex sheet
except near the leading edge, allows the leading and trailing
edge processes to be computed separately. The computed results
Where /„ is the wavelength of the vortex shedding process, U+ are combined following Rossiter.3 This leads to a phase shift of
is the forced convection velocity outside the cavity, and K is 2nn, a condition that selects the allowable frequencies. The rms
the fraction of this speed at which the vortex is convected pressure distribution along the cavity may also be computed,
across the cavity mouth. Taking the frequency of the acoustic and this result is compared with Heller's measurements4 (Fig. 4).
wave to be equal to a_/A fl implies a negligibly small recirculating The qualitative agreement lends support to the idea that this
velocity in the cavity. Rossiter3 finds this to be the case for simple model contains many of the features that occur in practice.
shallow cavities (L/D > 1) with transonic outside flowfields. Data Finally comparisons of predicted Strouhal numbers with data are
taken by Heller et al.4 suggest this approximation is valid at offered.
higher Mach numbers if L/D > 4. The second condition for
excitation is that the phase shift of the convected vortex traveling
from the leading edge to the trailing edge of the cavity, plus Analysis
the phase shift of the acoustic wave traveling the reverse direction For simplification, the analysis will be divided into three
in the cavity, is 2nn plus a correction factor for edge effects; parts. The first will be a study of the stability of a vortex sheet
thus subjected to a periodic pressure impulse. The second will be an
approximate study of the effects of an acoustic source or
Combining Eqs. (1) and (2) gives monopole located on the downstream face of the cavity on the
pressure field in the cavity. The third will be coupling the first
Si =• (3)
and second such that conditions for excitation are determined.
TjU.t)
MASS FLUX
FROM CAVITY
u_-o
00 e
HELLER 1
HANKEL
= lim-H P_ e~iKxdx +
Fig. 4 Comparison FCN. APPROX. -»0 (JO U -00
of normalized rms
pressure. P,e-iK
The kinematic and dynamic boundary conditions at the shear chosen such that Re(^L + ) ^ 0 in order to satisfy the radiation
layer q(x, t) can be combined to give condition as y —> oo.
Vortex sheet displacement function, q(x, t) is obtained from the
inverse transform of Eq. (15). A brief description of the inversion
y-»0_ p_ d2r]/dt2 = -dP_/dy (7) procedure is given in the Appendix. It is simply noted here that
P + (x,o,t) = P_(x,o,t) (8) the inversion in general requires mappings, which must be
where P_(x,o, t) contains the harmonic forcing (A/2e) e~io)t. The obtained numerically, of the zeros of the transcendental
boundary conditions as y -> +00 require outgoing radiation and dispersion relation
the rigid wall requires that A(Q,K)-( J u + 2 -K 2 )/> + +[(/^_ 2 -K 2 )/ A i _]ctnh/i_D (16)
in the Q and K planes. Unfortunately only the limiting case of
—— Eq. (16) as D -> oo allows analytic determination of the zeros as
dy y=-D
shown by Miles.9 Tarn10 has identified the amplifying and
Assuming the system can have and has achieved a steady state evanescent roots for this limit by using a stability criterion
implies that the initial conditions are not important. Thus they described by Briggs.12
can be set equal to zero. The general form of the steady state amplifying solution is,
The solution for the shear-layer displacement will be obtained for x > 0,
by the method of Fourier-Laplace Transforms. Transforming ] (17)
Eqs. (4-6, 8, and 9)§ yields
d2P+/dy2 = (S + 2/a + 2 + K2)P+ (10) The phase of the constant of proportionality is 0 and
represents the lag of vortex sheet displacement behind forcing
where at x = 0. 0 is the phase of \_-i(d/dK)^(o})K)~]~1\^(^K)=0. Only
amplifying solutions are considered here since they will dominate
d2P_/dy2 = (S_ the response downstream. The vortex sheet displacement for all
(11) waves can be shown to be exponentially small for x < 0.
where
1
This integral is evaluated exactly yielding
1
(2e-^4/S2) cos (S 2 )i
)
N f .
v^;
1
(24)
— 2L ——————°
Y Note that C(u) and S(u) are the Fresnal integrals. These functions
t
1 have been tabulated by Abramowitz and Stegun.14
1 1
Ultimately the phase of the pressure [Eq. (20)], evaluated at
the leading edge of the cavity, will be needed. The phase of
— iPei(0t has been evaluated for N = 2 and L/D ratios of 4 and oo.
The results are presented in Fig. 6 along with the results obtained
by neglecting R [Eq. (22)]. U N
It is apparent that the error incurred for the phase of the
O ACOUSTIC
MONOPOLE
pressure at the leading edge of the cavity by neglecting all
sources except the one located at the trailing edge is small for
cavities with L/D ratios in excess of 4. Therefore the phase shift
of an acoustic disturbance traveling from the trailing edge to the
o-2L*o o o
leading edge of a shallow cavity is taken to be that given by a
single source located at the trailing edge.
Downloaded by YORK UNIVERSITY LIBRARIES on August 12, 2014 | http://arc.aiaa.org | DOI: 10.2514/3.6747
•- - f -^ (31)
10.
This is just Rossiter's formula [Eq. (3)] with £ = -((/>/27i+f)
and K = a)/Kr U+. Note that 0 and f come from the leading
8.0 and trailing edge processes, respectively.
SOLID 8 DASHED
| 6.0 LINES-PHASE
OF EON. 21
O NEGLECTING R, Results
2
- 4.0
b
In the limit, D->oo the dispersion relation, Eq. (16), takes
a greatly simplified form. Tarn10 has shown that one amplifying
2.0
t The phase of the pressure at the leading edge of the cavity has
the same value for L/D ratios 0 and oo.
2.0 4.0 6.0 8.0 10. ** It is known that near a source at a distance of O(T/) the flow-
field to lowest order may be treated as incompressible if J///1 < 1.
Experimental results show that A is 0(L). Therefore the inequality
Fig. 6 Phase of Eq. (21). holds for small vortex sheet deflections.
MARCH 1973 EXCITATION FREQUENCIES FOR RECTANGULAR CAVITIES 351
1.0 I- PLANE
0.8
3.0
0.4
K-PLANE
0.2
1. Farzin Ghanadi, Maziar Arjomandi, Ben Cazzolato, Anthony Zander. 2014. Interaction of a flow-excited Helmholtz
resonator with a grazing turbulent boundary layer. Experimental Thermal and Fluid Science 58, 80-92. [CrossRef]
2. Ryan F. Schmit, James E. Grove, Frank Semmelmayer, Mitchell Haverkamp. Nonlinear Feedback Mechanisms Inside a
Rectangular Cavity. AIAA Journal, ahead of print1-16. [Abstract] [Full Text] [PDF] [PDF Plus]
3. Taro Handa, Hiroaki Miyachi, Hatsuki Kakuno, Takaya Ozaki, Shinji Maruyama. Modeling of a Feedback Mechanism
in Supersonic Deep-Cavity Flows. AIAA Journal, ahead of print1-6. [Abstract] [Full Text] [PDF] [PDF Plus]
4. Rajarshi Das, Heuy-Dong Kim, Job Kurian. Experimental Study of Supersonic Entrainment Using a Cavity. Journal of
Propulsion and Power, ahead of print1-9. [Abstract] [Full Text] [PDF] [PDF Plus]
5. Daming Sun, Ya Xu, Haijun Chen, Qie Shen, Xuejun Zhang, Limin Qiu. 2013. Acoustic characteristics of a mean flow
acoustic engine capable of wind energy harvesting: Effect of resonator tube length. Energy 55, 361-368. [CrossRef]
6. Sami Yamouni, Denis Sipp, Laurent Jacquin. 2013. Interaction between feedback aeroacoustic and acoustic resonance
mechanisms in a cavity flow: a global stability analysis. Journal of Fluid Mechanics 717, 134-165. [CrossRef]
7. Weipeng Li, Taku Nonomura, Kozo Fujii. 2013. On the feedback mechanism in supersonic cavity flows. Physics of Fluids
Downloaded by YORK UNIVERSITY LIBRARIES on August 12, 2014 | http://arc.aiaa.org | DOI: 10.2514/3.6747
34. Clarence W. Rowley, David R. Williams. 2006. DYNAMICS AND CONTROL OF HIGH-REYNOLDS-NUMBER
FLOW OVER OPEN CAVITIES. Annual Review of Fluid Mechanics 38:1, 251-276. [CrossRef]
35. W. Koch. 2005. Acoustic Resonances in Rectangular Open Cavities. AIAA Journal 43:11, 2342-2349. [Citation] [PDF]
[PDF Plus]
36. Lawrence Ukeiley, Nathan Murray. 2005. Velocity and surface pressure measurements in an open cavity. Experiments in
Fluids 38:5, 656-671. [CrossRef]
37. B. Umesh Chandra, S. R. Chakravarthy. 2005. Experimental Investigation of Cavity-Induced Acoustic Oscillations in
Confined Supersonic Flow. Journal of Fluids Engineering 127:4, 761. [CrossRef]
38. K. Takeda, C.M. Shieh. 2004. Cavity Tones by Computational Aeroacoustics. International Journal of Computational
Fluid Dynamics 18:6, 439-454. [CrossRef]
39. Nathan Murray, Lawrence UkeileyLow-Dimensional Estimation of Cavity Flow Dynamics . [Citation] [PDF] [PDF Plus]
40. Donald P. Rizzetta, Miguel R. Visbal. 2003. Large-Eddy Simulation of Supersonic Cavity Flowfields Including Flow
Control. AIAA Journal 41:8, 1452-1462. [Citation] [PDF] [PDF Plus]
41. S. Ziada, H. Ng, C.E. Blake. 2003. Flow excited resonance of a confined shallow cavity in low Mach number flow and
its control. Journal of Fluids and Structures 18:1, 79-92. [CrossRef]
42. K. Chung. 2003. Characteristics of Compressible Rectangular Cavity Flows. Journal of Aircraft 40:1, 137-142. [Citation]
[PDF] [PDF Plus]
43. Ganesh Raman, Edmane Envia, Timothy J. Bencic. 2002. Jet-Cavity Interaction Tones. AIAA Journal 40:8, 1503-1511.
[Citation] [PDF] [PDF Plus]
44. Nathan Murray, Lawrence UkeileyEstimating the Shear Layer Velocity Field Above an Open Cavity from Surface Pressure
Measurements . [Citation] [PDF] [PDF Plus]
45. H. KOOK, L. MONGEAU. 2002. ANALYSIS OF THE PERIODIC PRESSURE FLUCTUATIONS INDUCED BY
FLOW OVER A CAVITY. Journal of Sound and Vibration 251:5, 823-846. [CrossRef]
46. Jianbo Zhang, Etsuo Morishita, Takeo Okunuki, Hiroshi Itoh. 2002. Experimental Investigation on the Mechanism
of Flow-Type Changes in Supersonic Cavity Flows. TRANSACTIONS OF THE JAPAN SOCIETY FOR
AERONAUTICAL AND SPACE SCIENCES 45:149, 170-179. [CrossRef]
47. Kung-Ming Chung. 2000. Characteristics of Transonic Rectangular Cavity Flows. Journal of Aircraft 37:3, 463-468.
[Citation] [PDF] [PDF Plus]
48. C. NOGER, J.C. PATRAT, J. PEUBE, J.L. PEUBE. 2000. AEROACOUSTICAL STUDY OF THE TGV
PANTOGRAPH RECESS. Journal of Sound and Vibration 231:3, 563-575. [CrossRef]
49. Kung-Ming Chung. 1999. Pressure Fluctuations in Rectangular Cavity Flows. Journal of Mechanics 15:03, 97-102.
[CrossRef]
50. Drazen Fabris, David WilliamsExperimental measurements of cavity and shear layer response to unsteady bleed forcing .
[Citation] [PDF] [PDF Plus]
51. Y. Kang, A. R. Karagozian, O. I. Smith. 1998. Transport Enhancement in Acoustically Excited Cavity Flows, Part 1:
Nonreactive Flow Diagnostics. AIAA Journal 36:9, 1562-1567. [Citation] [PDF] [PDF Plus]
52. A.V. Savin. 1996. On mechanism of self-oscillatory interaction of plane subsonic jet with acoustic resonator. Experimental
Thermal and Fluid Science 12:4, 388-398. [CrossRef]
53. R. A. SMITH, E. GUTMARK, K. C. SCHADOW. 1992. Mitigation of pressure oscillations induced by supersonic flow
over slender cavities. Journal of Aircraft 29:6, 999-1004. [Citation] [PDF] [PDF Plus]
54. R. SMITH, E. GUTMARK, K. SCHADOWMitigation of pressure oscillations induced by supersonic flow over slender
cavities . [Citation] [PDF] [PDF Plus]
55. DONALD P. RIZZETTA. 1988. Numerical simulation of supersonic flow over a three-dimensional cavity. AIAA Journal
26:7, 799-807. [Citation] [PDF] [PDF Plus]
56. DONALD RIZZETTANumerical simulation of supersonic flow over a three-dimensional cavity . [Citation] [PDF]
[PDF Plus]
57. Shear Layer Instabilities, Flow Tones, and Jet Noise 1, 130-218. [CrossRef]
58. D. ROCKWELL. 1983. Oscillations of impinging shear layers. AIAA Journal 21:5, 645-664. [Citation] [PDF] [PDF Plus]
59. V. Sarohia, P. F. Massier. 1977. Control of Cavity Noise. Journal of Aircraft 14:9, 833-837. [Citation] [PDF] [PDF Plus]
60. Yung H. Yu. 1977. Measurements of Sound Radiation from Cavities at Subsonic Speeds. Journal of Aircraft 14:9, 838-843.
Downloaded by YORK UNIVERSITY LIBRARIES on August 12, 2014 | http://arc.aiaa.org | DOI: 10.2514/3.6747