Appendix III
Answers to problems
Problem 2.1.
From
Eq. (2.112), we see that the eigenfrequencies of the pan flute are in multiples of 2
n
+
1. In other words, they are in
odd multiples, and the even ones are missing. The fundamental resonance frequency is given by
f
=
c/(4
l)
=
345/(4
×
0.294)
=
293
Hz.
From
Eq. (2.94), we see that the eigenfrequencies of the pan flute are in multiples of
n. In other words, they are in
odd and
even multiples. The fundamental resonance frequency is given by
f
=
c/(2
l)
=
345/(2
×
0.294)
=
587
Hz.
Problem 2.2.
In the steady state, the
homogenous wave equation in cylindrical coordinates from
Eq. (2.23) becomes
(∂2∂w2+1w·∂∂w+k2)˜p(w)=0.
Although a solution is given by Eqs. (2.125) and (2.126), we omit the Y
0 function because of continuity at the center so that the solution becomes
where the J
0 function is analogous to the cosine function in the case of a closed tube and satisfies the boundary condition at the center
At the cylinder wall, where w
=
a, we have
1−jkρ0c·∂∂w˜p(w)|w=a=−˜u0,
which leads to
Hence,
˜p(w)=jρ0c˜u0J0(kw)J1(ka).
(A3.1)
Next, we shall rewrite the
inhomogeneous Eqs. (2.80) and (2.81) in cylindrical coordinates using the Laplace operator of
Eq. (2.23)
(∂2∂w2+1w·∂∂w+k2)˜p(w)=−δ(w−a)∂∂w˜p(w)|x=a=jkρ0cδ(w−a)˜u0,
(A3.2)
which includes the sound source at w
=
a on the right-hand side, where δ is the Dirac delta function. Let the solution be in the form of an eigenfunction expansion
˜p(w)=∑∞n=0˜AnJ0(βnw/a),
(A3.3)
where β
n
is the solution to
When the wall is stationary, this satisfies the boundary conditions
∂∂w˜p(w)|w=0=0,∂∂w˜p(w)|w=a=0.
Inserting Eq. (A3.3) into Eq. (A3.2) and multiplying both sides by
J0(βmw/a)
while integrating over the radius of the cylinder gives
∑∞n=0˜An(k2−β2na2)∫a0J0(βmw/a)J0(βnw/a)wdw=jkρ0c˜u0∫a0J0(βmw/a)δ(w−a)wdw.
The two integrals have the following identities:
∫a0J0(βmw/a)J0(βnw/a)wdw={0,m≠na2J20(βn)/2,m=n
∫a0J0(βmw/a)δ(w−a)wdw=aJ0(βm).
The first is the property of
orthogonality from Eq. (A2.101a) of Appendix II, and the second is a property of the Dirac delta function from Eq. (A2.154) of Appendix II. Hence,
J0(βmw/a)
in this case is the
orthogonal (or normalizing) function and we can eliminate the summation to yield
˜An=2jkaρ0c˜u0(k2a2−β2n)J0(βn)
so that the solution of Eq. (2), given by Eq. (3), becomes
˜p(w)=2jkaρ0c˜u0∑∞n=0J0(βnw/a)(k2a2−β2n)J0(βn).
(A3.4)
By comparing the pressures given by Eqs. (A3.1) and (A3.4) at w
=
a, we obtain
J0(ka)J1(ka)=∑∞n=02kak2a2−β2n.
Hence,
Q(x)=2J1(x)xJ0(x)=(∑∞n=0x2x2−β2n)−1.
For large
n, we can substitute
βn≈(n+1/4)π,n→∞.
Problem 7.1.
From
Eq. (6.10), we obtain
Q
TS
=
0.46
×
7.58/(0.46
+
7.58)
=
0.43. This is within 7% of the
Q
TS
value in
Table 7.4 required for a 0.01
dB Chebyshev alignment. Multiplying the
V
AB
/
V
AS
figure from the table by
V
AS
of the drive unit gives the box volume
V
AB
=
1.5511
×
19.5
=
30.3
L. Also, from the table, we see that the required port resonance frequency is
f
B
=
f
B
/
f
S
×
f
S
=
0.8838
×
50
=
44.2
Hz. From the “Summary of bass-reflex design” on p. 334, we obtain the maximum peak pressure
pmax=2√2×10(SPLmax20−5)=2.83Pa
and the maximum volume displacement require to produce that pressure
Vmax=rpmax2πf2Bρ0=1×2.832×3.14×44.22×1.18=195cm3.
We let the port volume be
V
P
=
10
V
max
=
1.95
L. The approximate tube length
t is given by
Eq. (7.97)
t≈cωB√VPVAB=3452×3.14×44.2√1.9530.3=31.5cm
where the cross-sectional area of the tube from
Eq. (7.96) is
S
P
=
V
P
/
t
=
1.96
×
10
−3/0.315
=
62.2
cm
2 and the diameter is
dP=2√SP/π=
8.9
cm. If we use a 9
cm diameter tube with a cross-sectional area
S
P
=
¼π
×
9
2
=
63.6
cm
2, the exact length is given by
Eq. (7.98)
t=SPc2VABω2B−0.84√SP=63.6×10−4×345230.3×10−3×(2×3.14×44.2)2−0.84√63.6×10−4=25.7cm.
Also, from
Table 7.4, we see that the cutoff frequency is
f
3dB
=
f
3dB/
f
S
×
f
S
= 0.8143
×
50
=
40.7
Hz.
Problem 10.1.
The total volume is simply
V
=
12.5
×
10
×
8
=
1000
m
3 and, from the table, the total absorptive area is
S
tot
=
100
+
125
+
32
=
257
m
2. The average Sabine absorption coefficient
α
tot
is given by
Eq. (10.52)
αtot=∑αs,iSiStot=0.7×100+0.6×125+0.5×32257=161257=0.626.
T=0.161VStotαtot=0.161V∑αs,iSi=0.161×1000257×0.626=1s,
which according to
Fig. 10.13 is fairly optimum for the size of auditorium. The sound strength
G is given by
Eq. (10.65)
G=10log10(100r2+1600π(1−αtot)Stotαtot)=10log10(100102+1600×3.14×(1−0.626)257×0.626)=11dB.
The distance
r
rev
from a source to equal direct and reverberant fields is given by
Eq. (10.68)
rev=14√QStotαtotπ(1−αtot)=14√1×257×0.6263.14×(1−0.626)=2.93m.
The average Eyring absorption coefficient
α
ey
is given by
Eq. (10.70)
αey=1−e−0.161V/(StotT)=1−2.72−0.161×1000/(257×1)=0.466.
The required acoustic power is given by
Eq. (10.71) using the maximum peak SPL of a large orchestra from column three of
Table 10.2
W=4×10(SPL/10)−10ρ0c(Q4πr2ref+4(1−αey)Stotαey)−1=4×10(102.5/10)−101.18×345(14π×102+4×(1−0.466)257×0.466)−1=0.938W.
If the loudspeakers have an efficiency of 0.5%, the required amplifier power is 0.938
×
100/0.5
=
188
W.
Problem 12.1.
The specific radiation impedance is given by
Zs=˜p(a)˜u0=jρ0cH(2)0(ka)H(2)1(ka),
which is separated into real and imaginary parts as follows:
Zs=jρ0cH(2)0(ka)H(1)1(ka)H(2)1(ka)H(1)1(ka)=ρ0cY0(ka)J1(ka)−J0(ka)Y1(ka)+j(J0(ka)J1(ka)+Y0(ka)Y1(ka))J21(ka)+Y21(ka).
Applying the Wronskian yields,
Zs=ρ0c(2πka(J21(ka)+Y21(ka))+jJ0(ka)J1(ka)+Y0(ka)Y1(ka)J21(ka)+Y21(ka)).
Problem 12.2.
The spherical cap is set in a rigid sphere of radius
R and moves with a radial velocity
˜u0
such that the velocity distribution is described by
˜u(R,θ)={˜u0,0≤θ≤α0,α<θ≤π,
(A3.1)
where α is the half angle of the arc formed by the cap. The total volume velocity is given by
˜U0=˜u0R2∫2π0∫α0sinθdθdφ=S˜u0,
(A3.2)
where S is the effective surface area of the cap given by
S=2πR2(1−cosα)=4πR2sin2(α/2).
(A3.3)
Near-field Pressure. Again, assuming that the pressure field generated is a general axisymmetric solution to
Eq. (2.180), the Helmholtz wave equation in spherical coordinates is
˜p(r,θ)=ρ0c˜u0∑∞n=0Anh(2)n(kr)Pn(cosθ).
(A3.4)
Applying the velocity boundary condition gives
˜u(R,θ)=1−jkρ0c∂∂rp(r,θ)|r=R=˜u0−jk∑∞n=0Anh′(2)n(kR)Pn(cosθ)={˜u0,0≤θ≤α0,α<θ≤π,
(A3.5)
where the derivative of the spherical Hankel function
h′(2)n(kR)
is given by
Eq. (12.32). We now multiply both sides of
Eq. (A3.5) with the orthogonal function
P
m
(cos
θ) and integrate over the surface of the sphere, where the area of each surface element is given by
δS
=
2
πR
2
sin
θ δθ, so that
1−jk∑∞n=0Anh′(2)n(kR)∫π0Pn(cosθ)Pm(cosθ)sinθdθ=∫α0Pm(cosθ)sinθdθ,
(A3.6)
from which we obtain the coefficients as follows:
An=−jk(2n+1)sinαP−1n(cosα)2h′(2)n(kR),
(A3.7)
where we have used the integral identities from Eqs. (A2.66) and (A2.69) in
Appendix II. Finally, by inserting Eq. (A3.7) in Eq. (A3.4), we can write the near-field pressure as
˜p(r,θ)=−jkρ0c˜u0∑∞n=0(n+12)sinαP−1n(cosα)Pn(cosθ)h(2)n(kr)h′(2)n(kR).
(A3.8)
Far-field Pressure. In the far field, we can use the asymptotic expression for the spherical Hankel function from
Eq. (12.18), which when inserted into Eq. (A3.8) gives
˜p(r,θ)|r→∞=˜p(r,θ)=−jkρ0cS˜u0e−jkr4πrD(θ),
(A3.9)
where S is the dome effective area given by Eq. (A3.3) and
D(θ)=sinα2k2R2sin2(α/2)∑∞n=0jn+1(2n+1)2P−1n(cosα)Pn(cosθ)nh(2)n−1(kR)−(n−1)h(2)n+1(kR).
(A3.10)
The far field on axis response is given by
D(0)=sinα2k2R2sin2(α/2)∑∞n=0jn+1(2n+1)2P−1n(cosα)nh(2)n−1(kR)−(n−1)h(2)n+1(kR).
(A3.11)
Radiation Impedance. The total radiation force
˜F
is given by
˜F=R2∫2π0∫α0˜p(R,θ)sinθdθdφ=−4πR2jρ0c˜u0∑∞n=0(n+12)2(sinαP−1n(cosα))2h(2)n(kR)nh(2)n−1(kR)−(n−1)h(2)n+1(kR)
(A3.12)
using the identity of Eq. (A2.69). The specific impedance Z
s
is then given by
Zs=˜F˜U0=−jρ0csin2αsin2(α/2)∑∞n=0(n+12)2(P−1n(cosα))2h(2)n(kR)nh(2)n−1(kR)−(n−1)h(2)n+1(kR),
(A3.13)
where we have used the expression for
˜U0
from Eq. (A3.2).
Problem 13.1.
˜p(x,y,z)=2jkρ0c˜u0∫ly/2−ly/2∫lx/2−lx/2g(x,y,x|x0,y0,z0)|z0=0dx0dy0,
g(x,y,x|x0,y0,z0)=e−jk(x(x−x0)+y(y−y0)+z(z−z0))/R4πR,
By combining the above equations, we get
˜p(x,y,z)=2jkρ0c˜u0e−jkR4πR∫lx/2−lx/2ejkxx0/Rdx0∫ly/2−ly/2ejkyy0/Rdy0=jklxlyρ0c˜u0e−jkR2πR·sin(klxx/2R)klxx/2R·sin(klyy/2R)klyy/2R.
We note that sin θ
x
=
x/R and sin θ
y
=
y/R which yields
˜p(x,y,z)=jklxlyρ0c˜u0e−jkR2πRD(θx,θy),
D(θx,θy)=sin(12klxsinθx)12klxsinθx·sin(12klysinθy)12klysinθy.
Problem 13.2.
˜p(x,y,z)=2jkρ0c˜u0∫ly/2−ly/2∫lx/2−lx/2g(x,y,x|x0,y0,z0)|z0=0dx0dy0,
g(x,y,x|x0,y0,z0)=−j8π2∫∞−∞∫∞−∞e−j(kx(x−x0)+ky(y−y0)+kz|z−z0|)kzdkxdky
kz={√k2−k2x−k2y,k2x+k2y≤k2−j√k2x+k2y−k2,k2x+k2y>k2.
By combining the above equations, we get
˜p(x,y,z)=kρ0c4π2˜u0∫∞−∞∫∞−∞e−j(kxx+kyy+kz|z|)kz∫lx/2−lx/2ejkxx0dx0∫ly/2−ly/2ejkyy0dy0dkxdky.
The total radiation force is then given by
˜f=∫ly/2−ly/2∫lx/2−lx/2˜p(x,y,z)|z=0dxdy.
The imaginary parts of the complex exponents cancel over the negative and positive values of k
x
and k
y
so that this simplifies to
˜f=kρ0c4π2˜u0∫∞−∞∫∞−∞∫lx/2−lx/2cos(kxx)dx∫ly/2−ly/2cos(kyy)dy×∫lx/2−lx/2cos(kxx0)dx0∫ly/2−ly/2cos(kyy0)dy01kzdkxdky.
The specific radiation impedance is the ratio of the total radiation force to the total volume velocity. Hence, after integrating over the surface of the piston, we have
Zs=˜f˜U0=˜flxly˜u0=klxlyρ0c4π2∫∞−∞∫∞−∞sinc2(kxlx/2)sinc2(kyly/2)1kzdkxdky,
where sinc(z)
=
sin(z)/z. By using polar coordinates, where k
x
=
kt cosϕ, k
y
=
kt sinϕ, and dk
x
dk
y
=
k
2
tdtdϕ, we reduce the double infinite integral to a single infinite one over t and a single finite one over ϕ. Also, the infinite integral is split to produce the real resistance
R
s
and the imaginary reactance X
s
so that
where
Rs=klxlyρ0c4π2∫π/20∫10sinc2(klx2tcosϕ)sinc2(kly2tsinϕ)tdtdϕ√1−t2,
Xs=klxlyρ0c4π2∫π/20∫∞1sinc2(klx2tcosϕ)sinc2(kly2tsinϕ)tdtdϕ√t2−1.
Substituting t
=
sin θ in the above equation yields
Rs=klxlyρ0c4π2∫π20∫π20sinc2(klx2sinθcosϕ)sinc2(kly2sinθsinϕ)sinθdθdϕ,Xs=klxlyρ0c4π2∫π20∫π2+j∞π2+j0sinc2(klx2sinθcosϕ)sinc2(kly2sinθsinϕ)sinθdθdϕ,
where we note that sin(π/2
+
j∞)
=
cos j∞
=
cosh ∞
=
∞.
Problem 13.3.
D(θ,ϕ)=2J1(kasinθ)kasinθ·sin(2kdsinθsinϕ)4sin(12kdsinθsinϕ)=2J1(kasinθ)kasinθ·cos(12kdsinθsinϕ)cos(kdsinθsinϕ)
Problem 14.1.
Let the Green's function be of the form
G(x,y|x0,y0)=∑∞m=1∑∞n=1Amnsin(mπx/lx)sin(nπy/ly).
Inserting this into
Eq. (14.24), multiplying both sides by sin(
pπx/
l
x
) sin(
qπy/
l
y
), and integrating over the surface of the membrane gives
∑∞m=1∑∞n=1Amn(k2D−m2π2l2x−n2π2l2y)∫lx0sin(mπxlx)sin(pπxlx)dx∫ly0sin(nπyly)sin(qπyly)dy=−∫lx0δ(x−x0)sin(pπxlx)dx∫ly0δ(y−y0)sin(qπyly)dy.
Applying the integral of
Eq. (14.18) together with the property of the Dirac delta function of Eq. (A2.154) from
Appendix II yields
Amn=4lxly(m2π2l2x+n2π2l2y−k2D)sin(mπx0lx)sin(nπy0ly).
Hence,
G(x,y|x0,y0)=4∑∞m=1∑∞n=1sin(mπx/lx)sin(nπy/ly)sin(mπx0/lx)sin(nπy0/ly)(k2mn−k2D)lxly,
where
Problem 14.2.
Let the Green's function be of the form
G(w,ϕ|w0,ϕ0)=∑∞m=0∑∞n=1Amncos(mϕ)Jm(αmnw/a).
Inserting this into
Eq. (14.44), multiplying both sides by cos(
pϕ)
J
p
(
α
pq
w/a), and integrating over the surface of the membrane gives
∑∞m=0∑∞n=1Amn(k2D−α2mna2)∫2π0cos(mϕ)cos(pϕ)dφ∫a0Jm(αmnw/a)Jp(αpqw/a)wdw=−∫2π0δ(ϕ−ϕ0)cos(pϕ)dϕ∫a0δ(w−w0)Jp(αpqw/a)dw.
Amn=2cos(mϕ0)Jm(αmnw0/a)π(1+δm0)(α2mn−k2Da2)J2m+1(αmn).
Hence,
G(w,ϕ|w0,ϕ0)=2π∑∞m=0∑∞n=1cos(mϕ)Jm(αmnw/a)cos(mϕ0)Jm(αmnw0/a)(1+δm0)(α2mn−k2Da2)J2m+1(αmn).
Problem 14.3.
f01=α201h4πa2√Y3(1−ν2)ρD=2.22152×0.2×10−34π×(12.5×10−3)2√69×1093×(1−0.32)×2700=1.54kHz,f02=(α02α01)2f01=(5.45162.2215)2668=9.26kHz,f03=(α03α01)2f01=(8.61142.2215)2668=23.1kHz,
Problem 14.4.
Let the Green's function be of the form
G(w,φ|w0,φ0)=∑∞m=0∑∞n=1Amnηmn(w,φ),
where
η
mn
(
w,
ϕ)
=
cos(
mϕ) (
J
m
(
α
mn
w/
a)
−
B
mn
I
m
(
α
mn
w/
a)
+
C
mn
). Inserting this into
Eq. (14.149), multiplying both sides by
η∗pq(w,ϕ)
, and integrating over the surface of the shell gives
∑∞m=0∑∞n=1Amn(α4mna4−k4D)∫2π0∫a0ηmn(w,φ)η∗pq(w,φ)wdwdϕ=−∫2π0∫a0δ(w−w0)δ(φ−φ0)η∗pq(w,φ)dwdϕ.
Applying the integrals of
Eq. (14.142) together with the property of the Dirac delta function of Eq. (A2.154) from
Appendix II yields
Amn=a2η∗mn(w0,φ0)πΔmn(k4Da4−α4mn).
Hence,
G(w,φ|w0,φ0)=a2π∑∞m=0∑∞n=1ηmn(w,φ)η∗mn(w0,φ0)Δmn(k4Da4−α4mn).
Problem 14.5.
f01=h4πa2√Y3ρD(β411−ν2+48H2h2)=40×10−64π×(12.5×10−3)2√69×1093×2700(−109.2×1031−0.32+48×(2×10−340×10−6)2)=0Hz,f02=40×10−64π×(12.5×10−3)2√69×1093×2700(5.903441−0.32+48×(2×10−340×10−6)2)=20.7kHz,f03=40×10−64π×(12.5×10−3)2√69×1093×2700(9.189841−0.32+48×(2×10−340×10−6)2)=21.3kHz,f04=0.2×10−34π×(12.5×10−3)2√69×1093×2700(12.38641−0.32+48×(2×10−340×10−6)2)=22.7kHz.
Problem 15.1.
prms(r,0)=2πfε0a2(EPd)212√2rc=2×3.14×103×8.85×10−12×0.052×(2000/10−3)22×√2×1×345=0.570Pa
Hence, the sound pressure level is 20 log
10(0.57/(20
×
10
−6))
=
89.1
dB SPL.
Problem 15.2.
EP=(EPd)d=200010−3×0.224×10−3=448V.
T=ε0a2d(EPd)2=8.85×10−12×0.0520.224×10−3(200010−3)2=395N/m.
f0=13.26a√Taρ0=13.26×0.05√3950.05×1.18=502Hz.
RS=2Q√ρ0Ta=21√1.18×3950.05=193rayls.
T=0.6σmaxh=0.6×55×106×12×10−6=396N/m.