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 2n   +   1. In other words, they are in odd multiples, and the even ones are missing. The fundamental resonance frequency is given by f   =   c/(4l)   =   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/(2l)   =   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
(2w2+1w·w+k2)˜p(w)=0.
image
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
˜p(w)=˜pJ0(kw),
image
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
w˜p(w)|w=0=0.
image
At the cylinder wall, where w   =   a, we have
1jkρ0c·w˜p(w)|w=a=˜u0,
image
which leads to
˜p=jρ0c˜u0J1(ka).
image
Hence,
˜p(w)=jρ0c˜u0J0(kw)J1(ka).
image (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)
(2w2+1w·w+k2)˜p(w)=δ(wa)w˜p(w)|x=a=jkρ0cδ(wa)˜u0,
image (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),
image (A3.3)
where β n is the solution to
J1(βn)=0.
image
When the wall is stationary, this satisfies the boundary conditions
w˜p(w)|w=0=0,w˜p(w)|w=a=0.
image
Inserting Eq. (A3.3) into Eq. (A3.2) and multiplying both sides by J0(βmw/a) image while integrating over the radius of the cylinder gives
n=0˜An(k2β2na2)a0J0(βmw/a)J0(βnw/a)wdw=jkρ0c˜u0a0J0(βmw/a)δ(wa)wdw.
image
The two integrals have the following identities:
a0J0(βmw/a)J0(βnw/a)wdw={0,mna2J20(βn)/2,m=n
image
a0J0(βmw/a)δ(wa)wdw=aJ0(βm).
image
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) image 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)
image
so that the solution of Eq. (2), given by Eq. (3), becomes
˜p(w)=2jkaρ0c˜u0n=0J0(βnw/a)(k2a2β2n)J0(βn).
image (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.
image
Hence,
Q(x)=2J1(x)xJ0(x)=(n=0x2x2β2n)1.
image
For large n, we can substitute βn(n+1/4)π,n. image
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=22×10(SPLmax205)=2.83Pa
image
and the maximum volume displacement require to produce that pressure
Vmax=rpmax2πf2Bρ0=1×2.832×3.14×44.22×1.18=195cm3.
image
We let the port volume be V P   =   10V max   =   1.95   L. The approximate tube length t is given by Eq. (7.97)
tcωBVPVAB=3452×3.14×44.21.9530.3=31.5cm
image
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   cm2 and the diameter is dP=2SP/π= image 8.9   cm. If we use a 9   cm diameter tube with a cross-sectional area S P   =   ¼π   ×   92   =   63.6   cm2, the exact length is given by Eq. (7.98)
t=SPc2VABω2B0.84SP=63.6×104×345230.3×103×(2×3.14×44.2)20.8463.6×104=25.7cm.
image
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   m3 and, from the table, the total absorptive area is S tot   =   100   +   125   +   32   =   257   m2. 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.
image
The reverberation time T is given by Eq. (10.50)
T=0.161VStotαtot=0.161Vαs,iSi=0.161×1000257×0.626=1s,
image
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×(10.626)257×0.626)=11dB.
image
The distance r rev from a source to equal direct and reverberant fields is given by Eq. (10.68)
rev=14QStotαtotπ(1αtot)=141×257×0.6263.14×(10.626)=2.93m.
image
The average Eyring absorption coefficient α ey is given by Eq. (10.70)
αey=1e0.161V/(StotT)=12.720.161×1000/(257×1)=0.466.
image
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×(10.466)257×0.466)1=0.938W.
image
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),
image
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).
image
Applying the Wronskian yields,
Zs=ρ0c(2πka(J21(ka)+Y21(ka))+jJ0(ka)J1(ka)+Y0(ka)Y1(ka)J21(ka)+Y21(ka)).
image
Problem 12.2.
The spherical cap is set in a rigid sphere of radius R and moves with a radial velocity ˜u0 image such that the velocity distribution is described by
˜u(R,θ)={˜u0,0θα0,α<θπ,
image (A3.1)
where α is the half angle of the arc formed by the cap. The total volume velocity is given by
˜U0=˜u0R22π0α0sinθdθdφ=S˜u0,
image (A3.2)
where S is the effective surface area of the cap given by
S=2πR2(1cosα)=4πR2sin2(α/2).
image (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˜u0n=0Anh(2)n(kr)Pn(cosθ).
image (A3.4)
Applying the velocity boundary condition gives
˜u(R,θ)=1jkρ0crp(r,θ)|r=R=˜u0jkn=0Anh(2)n(kR)Pn(cosθ)={˜u0,0θα0,α<θπ,
image (A3.5)
where the derivative of the spherical Hankel function h(2)n(kR) image 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
1jkn=0Anh(2)n(kR)π0Pn(cosθ)Pm(cosθ)sinθdθ=α0Pm(cosθ)sinθdθ,
image (A3.6)
from which we obtain the coefficients as follows:
An=jk(2n+1)sinαP1n(cosα)2h(2)n(kR),
image (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˜u0n=0(n+12)sinαP1n(cosα)Pn(cosθ)h(2)n(kr)h(2)n(kR).
image (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˜u0ejkr4πrD(θ),
image (A3.9)
where S is the dome effective area given by Eq. (A3.3) and
D(θ)=sinα2k2R2sin2(α/2)n=0jn+1(2n+1)2P1n(cosα)Pn(cosθ)nh(2)n1(kR)(n1)h(2)n+1(kR).
image (A3.10)
The far field on axis response is given by
D(0)=sinα2k2R2sin2(α/2)n=0jn+1(2n+1)2P1n(cosα)nh(2)n1(kR)(n1)h(2)n+1(kR).
image (A3.11)
Radiation Impedance. The total radiation force ˜F image is given by
˜F=R22π0α0˜p(R,θ)sinθdθdφ=4πR2jρ0c˜u0n=0(n+12)2(sinαP1n(cosα))2h(2)n(kR)nh(2)n1(kR)(n1)h(2)n+1(kR)
image (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(P1n(cosα))2h(2)n(kR)nh(2)n1(kR)(n1)h(2)n+1(kR),
image (A3.13)
where we have used the expression for ˜U0 image from Eq. (A3.2).
Problem 13.1.
The far-field pressure is given by Eq. (13.6)
˜p(x,y,z)=2jkρ0c˜u0ly/2ly/2lx/2lx/2g(x,y,x|x0,y0,z0)|z0=0dx0dy0,
image
where the Green's function is given by Eqs. (13.33) and (13.34)
g(x,y,x|x0,y0,z0)=ejk(x(xx0)+y(yy0)+z(zz0))/R4πR,
image
R=x2+y2+z2.
image
By combining the above equations, we get
˜p(x,y,z)=2jkρ0c˜u0ejkR4πRlx/2lx/2ejkxx0/Rdx0ly/2ly/2ejkyy0/Rdy0=jklxlyρ0c˜u0ejkR2πR·sin(klxx/2R)klxx/2R·sin(klyy/2R)klyy/2R.
image
We note that sin θ x   =   x/R and sin θ y   =   y/R which yields
˜p(x,y,z)=jklxlyρ0c˜u0ejkR2πRD(θx,θy),
image
where
D(θx,θy)=sin(12klxsinθx)12klxsinθx·sin(12klysinθy)12klysinθy.
image
Problem 13.2.
The near-field pressure is given by Eq. (13.6)
˜p(x,y,z)=2jkρ0c˜u0ly/2ly/2lx/2lx/2g(x,y,x|x0,y0,z0)|z0=0dx0dy0,
image
where the Green's function is given by Eqs. (13.33) and (13.34)
g(x,y,x|x0,y0,z0)=j8π2ej(kx(xx0)+ky(yy0)+kz|zz0|)kzdkxdky
image
kz={k2k2xk2y,k2x+k2yk2jk2x+k2yk2,k2x+k2y>k2.
image
By combining the above equations, we get
˜p(x,y,z)=kρ0c4π2˜u0ej(kxx+kyy+kz|z|)kzlx/2lx/2ejkxx0dx0ly/2ly/2ejkyy0dy0dkxdky.
image
The total radiation force is then given by
˜f=ly/2ly/2lx/2lx/2˜p(x,y,z)|z=0dxdy.
image
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˜u0lx/2lx/2cos(kxx)dxly/2ly/2cos(kyy)dy×lx/2lx/2cos(kxx0)dx0ly/2ly/2cos(kyy0)dy01kzdkxdky.
image
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π2sinc2(kxlx/2)sinc2(kyly/2)1kzdkxdky,
image
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
Zs=RsjXs,
image
where
Rs=klxlyρ0c4π2π/2010sinc2(klx2tcosϕ)sinc2(kly2tsinϕ)tdtdϕ1t2,
image
Xs=klxlyρ0c4π2π/201sinc2(klx2tcosϕ)sinc2(kly2tsinϕ)tdtdϕt21.
image
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ϕ,
image
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ϕ)
image
Problem 14.1.
Let the Green's function be of the form
G(x,y|x0,y0)=m=1n=1Amnsin(mπx/lx)sin(nπy/ly).
image
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=1n=1Amn(k2Dm2π2l2xn2π2l2y)lx0sin(mπxlx)sin(pπxlx)dxly0sin(nπyly)sin(qπyly)dy=lx0δ(xx0)sin(pπxlx)dxly0δ(yy0)sin(qπyly)dy.
image
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π2l2yk2D)sin(mπx0lx)sin(nπy0ly).
image
Hence,
G(x,y|x0,y0)=4m=1n=1sin(mπx/lx)sin(nπy/ly)sin(mπx0/lx)sin(nπy0/ly)(k2mnk2D)lxly,
image
where
kmn=πm2l2x+n2l2y.
image
Problem 14.2.
Let the Green's function be of the form
G(w,ϕ|w0,ϕ0)=m=0n=1Amncos(mϕ)Jm(αmnw/a).
image
Inserting this into Eq. (14.44), multiplying both sides by cos()J p (α pq w/a), and integrating over the surface of the membrane gives
m=0n=1Amn(k2Dα2mna2)2π0cos(mϕ)cos(pϕ)dφa0Jm(αmnw/a)Jp(αpqw/a)wdw=2π0δ(ϕϕ0)cos(pϕ)dϕa0δ(ww0)Jp(αpqw/a)dw.
image
Applying the integrals of Eqs. (14.38) and (14.39) together with the property of the Dirac delta function of Eq. (A2.154) from Appendix II yields
Amn=2cos(mϕ0)Jm(αmnw0/a)π(1+δm0)(α2mnk2Da2)J2m+1(αmn).
image
Hence,
G(w,ϕ|w0,ϕ0)=2πm=0n=1cos(mϕ)Jm(αmnw/a)cos(mϕ0)Jm(αmnw0/a)(1+δm0)(α2mnk2Da2)J2m+1(αmn).
image
Problem 14.3.
f01=α201h4πa2Y3(1ν2)ρD=2.22152×0.2×1034π×(12.5×103)269×1093×(10.32)×2700=1.54kHz,f02=(α02α01)2f01=(5.45162.2215)2668=9.26kHz,f03=(α03α01)2f01=(8.61142.2215)2668=23.1kHz,
image
Problem 14.4.
Let the Green's function be of the form
G(w,φ|w0,φ0)=m=0n=1Amnηmn(w,φ),
image
where η mn (w, ϕ)   =   cos() (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,ϕ) image , and integrating over the surface of the shell gives
m=0n=1Amn(α4mna4k4D)2π0a0ηmn(w,φ)ηpq(w,φ)wdwdϕ=2π0a0δ(ww0)δ(φφ0)ηpq(w,φ)dwdϕ.
image
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).
image
Hence,
G(w,φ|w0,φ0)=a2πm=0n=1ηmn(w,φ)ηmn(w0,φ0)Δmn(k4Da4α4mn).
image
Problem 14.5.
f01=h4πa2Y3ρD(β411ν2+48H2h2)=40×1064π×(12.5×103)269×1093×2700(109.2×10310.32+48×(2×10340×106)2)=0Hz,f02=40×1064π×(12.5×103)269×1093×2700(5.9034410.32+48×(2×10340×106)2)=20.7kHz,f03=40×1064π×(12.5×103)269×1093×2700(9.1898410.32+48×(2×10340×106)2)=21.3kHz,f04=0.2×1034π×(12.5×103)269×1093×2700(12.386410.32+48×(2×10340×106)2)=22.7kHz.
image
Problem 15.1.
prms(r,0)=2πfε0a2(EPd)2122rc=2×3.14×103×8.85×1012×0.052×(2000/103)22×2×1×345=0.570Pa
image
Hence, the sound pressure level is 20 log10(0.57/(20   ×   10 6))   =   89.1   dB SPL.
Problem 15.2.
From Eq. (15.119), the polarization voltage is
EP=(EPd)d=2000103×0.224×103=448V.
image
From Eq. (15.86), the minimum tension is
T=ε0a2d(EPd)2=8.85×1012×0.0520.224×103(2000103)2=395N/m.
image
From Eq. (15.114), the resonance frequency is
f0=13.26aTaρ0=13.26×0.053950.05×1.18=502Hz.
image
From Eq. (15.117), the resistance is
RS=2Qρ0Ta=211.18×3950.05=193rayls.
image
From Eq. (15.113), the maximum tension is
T=0.6σmaxh=0.6×55×106×12×106=396N/m.
image
..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset