Silencer design by using array resonators for low-frequency

band noise reduction

Sang-Hyun Seo
Department of Structure & Materials, KSLV Technology Division, Korea Aerospace Research Institute
(KARI), 45 Eoeun-Dong, Yuseong-Gu, Daejeon, 305-333, Korea

Yang-Hann Kim
Center for Noise and Vibration Control (NOVIC), Department of Mechanical Engineering, Korea Advanced
Institute of Science and Technology (KAIST), Science Town, Daejeon, 305-701, Korea

共Received 1 December 2004; revised 26 July 2005; accepted 26 July 2005兲

Helmholtz resonators are often used to reduce noise. They are particularly useful when noise has a
narrow frequency band. In this study we aim to broaden its narrow band characteristics by
combining many resonators. Serial and parallel arrangements of resonators have been tested to
obtain broader impedance mismatch characteristics in the broader band. Theoretical and
experimental results explain these characteristics in the absence of mean flow. The serial
arrangement mainly increases the peak of TL at the resonance frequency. But the parallel
arrangement logarithmically increases the peak of TL and expands the bandwidth. The change of
acoustic characteristics is explained by introducing an “equivalent impedance analysis.” This shows
that the transmission loss has a maximum value when the distance between resonators is ␭ / 4 of its
wavelength. In this study we propose a novel design method that optimizes the arrangement of
resonators while keeping the volume minimized. Various transmission loss characteristics are
possible when we select different objective functions under constraints. © 2005 Acoustical Society
of America. 关DOI: 10.1121/1.2036222兴
PACS number共s兲: 43.55.Rg, 43.50.Gf 关NX兴 Pages: 2332–2338

I. INTRODUCTION Therefore, we need to study the acoustic characteristics

of array resonators for the design of an effective silencer
The Helmholtz resonator is often used to reduce noise in model. Thus we intend to introduce the equivalent imped-
a narrow frequency band. This type of resonator has a high ance that represents the total acoustic characteristics of the
transmission loss in a narrow band at its resonance fre- silencer and an optimization for determining resonance
quency. It is easy to design this resonator to have a desired frequencies.
resonance frequency, because this is determined by the geo-
metric ratio of the cavity and its neck.
II. THEORY
There have been many studies about the disagreement
between analytic and experimental results of the resonator A. Transmission loss of a branch resonator
due to the uncertainty of damping around the neck or end
Silencer using a resonator reduces noise by an imped-
correction.1–3 With regard to characteristics of the silencer ance mismatch that causes reflection of the incident acoustic
using resonators, Anderson4 studied the effect of flow when a energy and attenuation in the resonator’s neck. When a reso-
single side branch Helmholtz resonator is attached to a cir- nator is attached to a duct by a side branch, as depicted in
cular duct. Chun5 studied the acoustic characteristics of con- Fig. 1, the basic assumption is that plane waves propagate in
centric pipe resonators to obtain a general solution for the a duct and the reflected waves from downstream of a duct do
acoustic field in the resonators by using the Green’s function. not exist by using an anechoic termination in the absence of
Furthermore, many studies about the applications of a reso- mean flow. Considering effects of grazing flow, if the mean
nator are in progress in various fields. flow’s velocity is less than M = 0.1 共M: Mach number兲, its
In this paper, multiple array resonators are used to effect is not serious.4
broaden the narrow band characteristics of a resonator in the The sound pressure 共P兲 and the volume velocity 共U兲 can
low-frequency band. Kim6 and Kuntz7 studied the acoustic be expressed as follows:
characteristics of the panel using multiple array resonators to
P1 = Ae−jkx + Be jkx, P2 = Ce−jkx, 共1a兲
obtain a high absorption coefficient at low frequency. These
studies made use of impedance changes in the resonator 1 1
configuration. U1 = 共Ae−jkx − Be jkx兲, U2 = 共Ce−jkx兲, 共1b兲
Z Z
Also, in the case of a silencer using array resonators, in
order to obtain broadband transmission loss characteristics, where A, B, and C are the magnitude of the incident wave,
we can make use of an additional effect from many different reflected wave, and transmitted wave, respectively, and Z
impedance surfaces. = ␳c / S is the acoustic impedance of the duct. Here k

2332 J. Acoust. Soc. Am. 118 共4兲, October 2005 0001-4966/2005/118共4兲/2332/7/$22.50 © 2005 Acoustical Society of America
 FIG. 1. Helmholtz resonator with a duct.

= 2␲ f / c is the wave number, ␳ is the density of air, and c is

the sound speed.
The transfer matrix between point 1 and point 2 can be
obtained as follows by using the continuity of the sound
pressure and the volume velocity:

冉 冊 冤 冥冉 冊
1 0
P1 P2
= 1
U1 1 U2
− jZc cot kh + Zh

= 冋 T11 T12
T21 T22
册冉 冊 P2
U2
, 共2兲
FIG. 2. The TL change varying the number of resonators. 共a兲 The serial
arrangement 共S: cross-sectional area of duct, L: distance of resonators兲; 共b兲
the parallel arrangement. 共f 0 = 340 Hz, L : 100 mm, l : 25 mm, D : 10 mm,
Sc : ␲共20兲2 mm2, h : 48.5 mm, S : 50*50 mm2兲.
where Zc = ␳c / Sc is the acoustic impedance of a resonator’s
cavity.
The impedance of resonator Zr can be expressed as The TL of a branch resonator has a peak at its resonance
frequency, as is well known. Thus, in order to obtain broad-
Zr = − jZc cot kh + Zh , 共3兲
band characteristics, it is necessary to study the TL change
for the serial and parallel array of resonators.
␳c
Zh = 关0.0072 + jk共l + 0.75兲兴, 共4兲
Sh

where Zh is the hole impedance of a resonator, as suggested B. Serial and parallel array
by Sullivan.8,9 Sullivan obtained the hole impedance of per- Figure 2 illustrates the serial and parallel arrangement of
forated elements in a concentric tube resonator by measure- identical resonators. Considering the TL change varying the
ment. Here, the resistance of Zh was modified in consider- resonator’s number, the serial arrangement mainly increases
ation of the experimental results of this study, and Sh is the the magnitude of TL at the resonance frequency. But the
cross-sectional area of the hole. parallel arrangement logarithmically increases the magnitude
Transmission loss 共TL兲 can be represented as follows by of TL and the bandwidth.
using the transfer matrix: These characteristics give the resonator array high TL in

TL = 20 log10 冏冏 A
C
the broader band. In this sense, we can design the silencer by
using multiple array resonators having different resonance
frequencies, which are included in the objective frequency

= 20 log10 冏 T11 + T12/Z + T21 · Z + T22

2
冏 共dB兲, 共5兲
band.
Figure 3 illustrates a panel-type silencer using each of 8
resonators 共i = 1; 257 Hz, i = 2; 297 Hz, i = 3; 361 Hz, i = 4;
where A is acoustic pressure of the incident wave and C is 413 Hz兲 in the parallel direction. Impedance of each resona-
the acoustic pressure of the transmitted wave. tor is expressed as follows, as given by Eq. 共3兲:
Using Eqs. 共2兲 and 共5兲, the TL of a branch resonator in Zr共i兲 = − jZc共i兲 cot kh共i兲 + Zh共i兲 共i = 1,2,3,4兲. 共7兲
the duct can be obtained as

冉 冊
If the transfer matrix that represents a resonator having

冨 冨
1 resonance frequency 共f i兲 is T f iHz and the direct tube between
2+Z
− jZc cot kh + Zh the resonators is Ttube, the total transfer matrix 共Ttotal兲 of the
TL = 20 log10 . 共6兲
2 panel-type silencer is obtained as follows:

J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array 2333
 FIG. 5. The experiment setup for measurement of transmission loss.

FIG. 3. Panel-type silencer using resonator array 共Li: distance of resonators兲.

In this study, TL was measured by a two-microphone
method from the Seybert and Ross wave decomposition
Ttot = T257 Hz · Ttube · T297 Hz · Ttube · T361 Hz · Ttube · T413 Hz theory.10

冤 冥 冤 冥 冤 冥
1 0 1 0 1 0
= N · Ttube · N · Ttube · N · Ttube B. Results
1 1 1
Zr1 Zr2 Zr3
Figure 6 illustrates the panel-type silencer using each of

冤 冥
1 0 four resonators 共i = 1; 257 Hz, i = 2; 334 Hz, i = 3; 387 Hz, i
· N , 共8兲 = 4; 447 Hz兲 in the parallel direction. Figure 7 shows a com-
1 parison of the measured and predicted transmission loss.
Zr4
where N denotes the number of resonators in the parallel
direction. IV. EQUIVALENT IMPEDANCE ANALYSIS
Using Eqs. 共5兲 and 共8兲, the total TL of the panel-type A. Parallel combination
silencer can be obtained. Figure 4 compares each resonator’s
TL with the total TL of the panel-type silencer. As depicted The panel-type silencer has two important design param-
in Fig. 4, the total TL has each resonance frequency’s peak eters. Each resonator’s resonance frequency and distance be-
and shows the additional effect. tween the resonator’s hole determine the TL of the panel-
Thus we confirmed that the panel-type silencer can type silencer. This means that we have to study the acoustic
obtain high TL characteristics in the broader band, as characteristics of the array resonators to determine an effec-
noted before. This theoretical result should then be verified tive resonator arrangement.
experimentally. Using an equivalent impedance concept that represents
the total acoustic characteristics of the silencer, we can ex-
plain the acoustic characteristics according to the resonator
III. EXPERIMENT
arrangement. Figure 8 illustrates the equivalent impedance of
A. Experimental setup a silencer using array resonators.
The equivalent impedance 共Zeq兲 for a parallel arrange-
Figure 5 shows the experimental setup for the TL mea-
ment is obtained as follows.
surement of the simplified panel-type silencer. The duct’s
As depicted in Fig. 1, the pressure reflection coefficient
cross-sectional area is 50⫻ 50mm2 with a higher cutoff fre-
and the pressure transmission coefficient at x = 0 can be rep-
quency of 3.43 kHz for plane wave propagation. An
resented as follows:
anechoic termination with a cutoff frequency of about
100 Hz is employed in this experiment to avoid the effects of B C
Rcoeff = , Tcoeff = . 共9兲
the reflected waves. A A
Considering the number of resonators in the parallel direc-
tion 共N兲, the transfer matrix between point 1 and point 2 can
be written as

FIG. 4. Transmission loss of panel-type silencer. 共N : 8 , l : 25 mm, D : 10 mm,

S : 50*50 mm2, L : 100 mm兲. FIG. 6. Panel-type silencer using each of four resonators.

2334 J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array
FIG. 7. A comparison between the experimental and theoretical result of the FIG. 9. Impedance change of the parallel arrangement varying the number
panel type using each of four resonators 共L : 100 mm, l : 25 mm, D : 10 mm, of resonators 共Re: resistance, Im: reactance兲.
S : 50*50 mm2兲.
B. Serial combination

冉 冊 冤 冥冉 冊
1 0 As depicted in Fig. 8, the total transfer matrix can be
P1 P2 expressed as follows for a serial array:
共10兲

冉 冊冋 册冉 冊
= N .
U1
Zr
1 U2

That is, Eq. 共10兲 can be rearranged as

冉 冊
P1
U1
= T1 · T2 · T3 · ¯
PN
UN
=
T11 T12
T21 T22
PN
UN
.

共15兲
P 1 = P 2, 共11a兲
Equation 共15兲 can be rewritten as
N · P2 P1 = T11 · PN + T12 · UN , 共16a兲
U1 = + U 2. 共11b兲
Zr
U1 = T21 · PN + T22 · UN . 共16b兲
Using Eqs. 共1兲, 共9兲, and 共11兲 at x = 0, the pressure reflec-
tion coefficient 共Rcoeff兲 and the pressure transmission coeffi- Using Eqs. 共9兲 and 共16兲, the pressure reflection coeffi-
cient 共Tcoeff兲 at the impedance surface can be obtained in cient 共Rcoeff兲 and the pressure transmission coefficient 共Tcoeff兲
terms of the impedance of the resonator 共Zr兲 and the acoustic can be obtained in terms of the element of total transfer
impedance 共Z兲 of the duct: matrix at x = 0 ⬃ L:
−N·Z 2
Rcoeff = , 共12兲 Tcoeff = , 共17a兲
N · Z + 2 · Zr T11 + T12/Z + T21 · Z + T22

2 · Zr T11 + T12/Z − T21 · Z − T22

Tcoeff = . 共13兲 Rcoeff = . 共17b兲
N · Z + 2 · Zr T11 + T12/Z + T21 · Z + T22

Using Eq. 共12兲, the equivalent impedance 共Zeq兲 for a parallel Figure 10 illustrates an electrical analogy11 of a simple
arrangement can be evaluated as model using two resonators.
First, considering the second resonator and a direct tube
P1 Z共A + B兲 Z共1 + Rcoeff兲 Z · Zr between the first and the second resonator, the impedance at
Zeq = = = = . 共14兲
U1 A−B 1 − Rcoeff N · Z + Zr x = L can be expressed as
Figure 9 shows the change of equivalent impedance Z · Z2
ZL = . 共18兲
共Zeq / Z兲, which is normalized by the acoustic impedance of Z + Z2
the duct while varying the number of resonators in the par-
allel direction. As the number of resonators increases, the As depicted in Fig. 10, considering the wave propagation
resistance 共Re兲 of the impedance broadens. The reactance having magnitudes A and B in the direct tube, the sound
共Im兲 of the impedance has a value of zero at resonance fre- pressure can be represented as
quency 共200 Hz兲.

FIG. 8. Equivalent impedance of a silencer using array resonators 共P: sound

pressure, U: volume velocity, A: incident wave, B: reflected wave, C: trans- FIG. 10. Electrical analogy of the simple model using two resonators 共P:
mitted wave兲. sound pressure, U: volume velocity兲.

J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array 2335
 P共x兲 = Ae−jk共x−L兲 + Be jk共x−L兲 . 共19兲
Using Eq. 共19兲, the impedance at x = 0, x = L can be obtained
as
P共L兲 A+B
ZL = =Z· , 共20兲
U共L兲 A−B

P共0兲 Ae jkL + Be−jkL

Z0 = = Z · jkL . 共21兲
U共0兲 Ae − Be−jkL
Using Eqs. 共20兲 and 共21兲, the impedance 共Z0兲 at x = 0 can be
evaluated as
共ZL/Z兲 + j tan共kL兲
Z0 = Z · . 共22兲
1 + j共ZL/Z兲tan共kL兲
Finally, considering the first resonator, the equivalent
impedance 共Zeq兲 of a serial array can be obtained as follows:

Zeq = 冉 1
+
1
Z1 Z0
冊 −1
. 共23兲

The change of acoustic characteristics is explained by using

the equivalent impedance analysis. Thus, we can determine
the effective resonator arrangement for obtaining the TL
characteristics for low-frequency band noise reduction.

C. Characteristics change by the distance between

resonators
As depicted in Fig. 10, in the case of using two resona-
tors with 202 Hz 共front兲 and 258 Hz 共rear兲, the power trans-
mission coefficient 共T⌸ = 兩Tcoeff兩2兲 and the power reflection
coefficient 共R⌸ = 兩Rcoeff兩2兲, the absorption coefficient 共D⌸ = 1
− R⌸ − T⌸兲 for varying the distance between resonators can
be obtained as the contour graph in Fig. 11.
It has a periodic characteristic of ␭ / 2 共␭: wavelength兲 as
a whole and the TL between the first and second resonance
frequency has a maximum value when the distance between
resonators is ␭ / 4 of its wavelength.
Figure 11 shows that the front resonator reflects sound
power uniformly and the rear resonator attenuates sound
power dominantly. When the arrangement is changed to FIG. 11. Acoustic characteristics varying the distance between resonators
with 202 Hz 共front兲 and 258 Hz 共rear兲. 共a兲 Power transmission coefficient;
258 Hz 共front兲 and 202 Hz 共rear兲, the power reflection coef-
共b兲 power reflection coefficient; 共c兲 absorption coefficient.
ficient 共R⌸兲 and the absorption coefficient 共D⌸兲 have a sym-
metric shape of 202 Hz 共front兲 and 258 Hz 共rear兲. However,
the TL is not changed in this case. D. Effective arrangement of resonators
Thus, the change of distance between resonators is more
According to the equivalent impedance analysis, we can
important than the arrangement order of the resonators from
determine the optimal distance between resonators to obtain
the viewpoint of TL. These characteristics are identified by
a high TL in the objective frequency band.
the equivalent impedance analysis, as noted earlier.
In the model shown in Fig. 3 共i = 1; 257 Hz, i = 2;
Figure 12 shows the resistance 共Re兲 and reactance 共Im兲
297 Hz, i = 3; 361 Hz, i = 4; 413 Hz兲, the distance between
of equivalent impedance that is normalized by the acoustic
resonator’s hole can be determined as
impedance of the duct while varying the distance between

冉 冊
resonators for the case of two resonators, 202 Hz 共front兲 and
␭i ␭i+1 1 1
258 Hz 共rear兲. The resistance of impedance is related to the + c +
energy dissipation and the reactance of impedance is related 4 4 f i f i+1
Li = = . 共24兲
to the energy reflection. Figure 11 and Fig. 12 show that the 2 8
absorption coefficient is high at the high resistance and the
reflection coefficient is high at the zero-crossing point of the Figure 13 shows a comparison of TL curves when the
reactance. distances of the resonators are 100 mm and ␭ / 4.

2336 J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array
 FIG. 14. The silencer model for the minimization of the design parameters
共W ⫻ H ⫻ Li: volume of cavity, l: neck length of resonator, D: hole diameter
of resonator, ⌬si: distance between resonators, Ltot: total length of silencer兲.

nators: ⌬si兲. If all design parameters are considered, the

silencer model becomes very complex. Therefore, it is nec-
essary to minimize the design parameters. Figure 14 illus-
trates the silencer model that minimizes the design param-
eters. Here l, W, H, and D can be fixed by geometric shape.
The resonance frequency 共f i兲 can be expressed as

fi =
c
2␲
冑 ␲D2/4
共l + 0.75D兲WHLi
共i = 1,2, . . . ,n兲, 共25兲

where n is the number of resonators in the serial direction.

Thus the function of TL can be defined as
TL = fun共Li,⌬si兲. 共26兲
FIG. 12. Equivalent impedance varying the distance between resonators
with 202 Hz 共front兲 and 258 Hz 共rear兲. 共a兲 Resistance; 共b兲 reactance. From Eq. 共25兲, Li is expressed in terms of f i as follows:
c2 · D2 A
When the distances between resonators are ␭ / 4, the TL Li = =
has a higher value in the objective frequency band. But the 16␲共l + 0.75D兲W · H · f 2i f 2i
silencer must have a long length for in order to effectively
control noise at low frequency because the wavelength is
relatively long. Therefore we intend to propose a new design
冉 A=
c2 · D2
16␲共l + 0.75D兲W · H
, 冊 共27兲

method that optimizes the arrangement of the resonators. where A is constant.

The distance between holes 共⌬si兲 also can be expressed
V. DESIGN PROBLEM in terms of f i on the assumption that the resonator’s hole is
located at the center of the cavity:
A. Problem definition
Li + Li+1 A/f 2i + A/f i+1
2
A共f 2i + f i+1
2
兲
The silencer using array resonators has many design pa- ⌬si = = = 2 2 , 共28兲
2 2 2f i · f i+1
rameters 共volume of cavity: W ⫻ H ⫻ Li, neck length of reso-
nator: l, hole diameter of resonator:D, distance between reso- From Eqs. 共26兲–共28兲, the function of TL can be repre-
sented as
TL = fun共f 1, f 2, f 3, . . . , f n兲. 共29兲
If the total volume of all resonators is given, the constraint of
the silencer design can be expressed as
n

兺
i=1
Li = Ltot . 共30兲

Using Eq. 共27兲, this constraint can be rewritten as

n
A
兺
i=1 f
2 = Ltot . 共31兲
i

Thus, we can obtain the multivariable optimization problem

FIG. 13. Transmission loss of panel-type silencer when the distances be-
tween resonators are 100 mm and ␭ / 4 共N : 8, l : 25 mm, D : 10 mm, in terms of resonance frequency with a total volume
S : 50*50 mm2兲. constraint.

J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array 2337
 VI. CONCLUSION

Helmholtz resonators are used to reduce noise in the

narrow frequency band. The transmission loss characteristics
of a silencer through the use of many resonators are obtained
by theoretical and experimental results. The results show that
it is possible to design a silencer to have broadband charac-
teristics through serial and parallel arrangements of resona-
tors at low frequency.
FIG. 15. The objective function for optimization 共f lb: low bound frequency, The change of acoustic characteristics of the silencer
f ub: upper bound frequency兲. model can be explained by using an equivalent impedance
analysis. Considering the equivalent impedance analysis, the
B. Optimization change of the distance between resonators is more significant
than the arrangement order of resonators from the viewpoint
For the desired TL, it is necessary to select a suitable
of TL. As a result, it has a periodic characteristic of ␭ / 2 as a
objective function. The determination of the objective func-
whole and the transmission loss has a maximum value when
tion for the design of TL shape is very important.
the distance between resonators is ␭ / 4 of its wavelength.
The optimization problem for TL shape having broad-
In this study we propose a new design method that op-
band characteristics at low frequency can be expressed as
timizes the arrangement of resonators for a TL shape that has
minimize max共TL兲 − min共TL兲:objective function broadband characteristics in low frequency. Thus we can ob-
tain various transmission loss shape by selecting a suitable
n
A objective function under given constraints.
subject to 兺
i=1 f
2 = Ltot:equality constraint
i

f lb 艋 f i 艋 f ub:bound constraint. 共32兲 ACKNOWLEDGMENTS

Figure 15 illustrates the objective function in Eq. 共32兲. This study was partly supported by the NRL 共National
When the bound constraint is 300 Hz艋 f i 艋 550 Hz and Research Laboratory兲 project of KISTEP 共Korea Institute Of
each parameter is SD = 50⫻ 50 mm2, W = 70 mm, H = 30 mm, Science and Technology Evaluation and Planning兲 and the
D = 10 mm, l = 10 mm, Ltot = 30 mm, the results of optimiza- BK21 共Brain Korea 21兲 project initiated by Ministry of Edu-
tion are as presented in Fig. 16. cation and Human Resources Development of Korea.
Figure 16 shows that each resonance frequency is 336, 1
342, 358, 385, 414, 448, and 481 Hz, and the difference K. Uno Ingard, “On the theory and design of acoustic resonators,” J.
Acoust. Soc. Am. 25, 1037–1062 共1953兲.
between the maximum and minimum of total TL is less than 2
J.-G. Ih, “On the inertial end correction of resonators,” Acustica 78, 1–15
4.7 dB. 共1993兲.
3
From these results, we can obtain various TL character- S.-T. Haa and Y.-H. Kim, “The limitation and applicability of Helmholtz
istics by selecting different objective functions under con- resonator, regarding as equivalent single-degree of vibration system,”
Journal of the Korean Society for Noise and Vibration Engineering 4,
straints. Finally, this optimization method shows that it is 209–219 共1994兲.
4
possible to design a transmission loss shape of a silencer J. S. Anderson, “The effect of an air flow on a single side branch Helm-
using array resonators. holtz resonator in a circular duct,” J. Sound Vib. 52, 423–431 共1977兲.
5
Y.-D. Chun and Y.-H. Kim, “A study on the characteristics of concentric
pipe resonators,” Journal of the Korean Society for Noise and Vibration
Engineering 8, 49–56 共1998兲.
6
S.-R. Kim, M.-S. Kim, Y.-H. Kim, and Y.-W. Kim, “Absorptive charac-
teristics of resonator panel for low frequency noise control, and its appli-
cations,” Inter-noise 2003, Seogwipo, Korea, N932, 2003.
7
R. A. Prydz, L. S. Wirt, and H. L. Kuntz, “Transmission loss of a
multilayer panel with internal tuned Helmholtz resonators,” J. Acoust.
Soc. Am. 87, 1597–1602 共1990兲.
8
J. W. Sullivan and M. J. Crocker, “Analysis of concentric-tube resonators
having unpartitioned cavities,” J. Acoust. Soc. Am. 64, 207–215 共1978兲.
9
J. W. Sullivan, “A method for modeling perforated tube muffler compo-
nents,” J. Acoust. Soc. Am. 66, 772–788 共1979兲.
10
A. F. Seybert and D. F. Ross, “Experimental determination of acoustic
properties using a two-microphone random-excitation technique,” J.
Acoust. Soc. Am. 61, 1362–1370 共1977兲.
11
F. A. Firestone, “A new analogy between mechanical and electrical sys-
FIG. 16. The result of optimization for the broadband TL shape. tems,” J. Acoust. Soc. Am. 4, 249–267 共1932兲.

2338 J. Acoust. Soc. Am., Vol. 118, No. 4, October 2005 S. Seo and Y. Kim: Silencer design by using array