7.3: The Speed of Sound
( \newcommand{\kernel}{\mathrm{null}\,}\)
The physics of sound waves is obviously a three-dimensional problem. However, we can learn a lot about sound by considering motion of air in only one-dimension. Consider, for example, standing waves in the air in a long narrow tube like an organ pipe, shown in cartoon form in Figure 7.3. Here, we will ignore the motion of the air perpendicular to the length of the pipe, and consider only the one-dimensional motion along the pipe. As we will see later, when we can deal with three-dimensional problems, this is a sensible thing to do for low frequencies, at which the transverse modes of oscillation cannot be excited. If we consider only one-dimensional motion, we can draw an analogy between the oscillations of the air in the pipe and the longitudinal waves in a massive spring.
Figure 7.3: An organ pipe.
It is clear what the analog of ρL is. The linear mass density of the air in the tube is ρL=ρA
where A is the cross-sectional area of the tube. The question then is what is Kℓ for a tube of air?
Consider putting a piston at the top of the tube, as shown in Figure 7.4. With the piston at the top of the tube, there is no force on the piston, because the pressure of the air in the tube is the same as the pressure of the air in the room outside. However, if the piston is moved in a distance dz, as shown Figure 7.5, the volume of the air in the tube is decreased by −dV=Adz.
Figure 7.4: The organ pipe with a piston at the top. The air in the tube acts like a spring.
Figure 7.5: Pushing in the piston changes the volume of the air in the tube.
If the piston were moved in slowly enough for the temperature of the gas to stay constant, then the pressure would simply be inversely proportional to the volume. However, in a sound wave, the motion of the air is so rapid that almost no heat has a chance to flow in or out of the system. Such a change in the volume is called “adiabatic.” When the volume is decreased adiabatically, the temperature goes up (because the force on the piston is doing work) and the pressure increases faster than 1/V, like p∝V−γ
where γ is a positive constant that depends on the thermodynamic properties of the gas. More precisely, γ is the ratio of the specific heat at constant pressure to the specific heat at constant volume:3 CP/CV
In air, at standard temperature and pressure γair ≈1.40
Now we can write from (7.34), dpp=−γdVV
or dp=−γpdVV≈γAp0Vdz=γp0ℓdz
where p0 is the equilibrium (room) pressure. Then the force on the piston is dF=Adp=γA2p0Vdz=γAp0ℓdz
so that K=dFdz=γAp0ℓ
and Kℓ is Kℓ=γAp0.
Thus we expect the dispersion relation to be ω2=v2sound k2=KℓρLk2=γp0ρk2
where we have defined the “speed of sound”, vsound , as v2sound =γp0ρ
For air at standard temperature and pressure, vsound ≈332ms.
As we will see in the next chapter, this is actually the speed at which sound waves travel. For now, it is just a parameter in our calculation of the normal modes.
In the pipe shown in (7.3), the displacement of the air, which we will call ψ(z,t), must vanish at z=0, because the bottom of the tube is closed and there is nowhere for the gas to go.
The z derivative of ψ must vanish at z=ℓ, because the excess pressure is proportional to −∂∂zψ. The pressure is proportional to the force in our analogy with longitudinal waves in ∂z the massive spring. Using (7.41) and (7.23), we expect the longitudinal force to be ±γAp0∂∂zψ
or the excess pressure to be p−p0=−γp0∂∂zψ.
We want the negative sign because for ∂∂zψ>0, the air is spreading out and has lower pressure.
Thus for a standing wave in the pipe, (7.3), we expect the boundary conditions ψ(0,t)=0,∂∂zψ(z,t)|z=ℓ=0,
for which the solution is ψ(z,t)=sinkzcosωt
k=(n+1/2)πℓ,ω=vk,
where v=vsound , for nonnegative integer n. In particular, the lowest frequency mode of the tube corresponds to n=0, ω=vπ2ℓ,ν=ω2π=v4ℓ.
Helmholtz Approximation
Let’s consider a slightly different problem. What is the lowest frequency mode of a one-liter soda bottle, shown in Figure 7.6? A typical set of parameters is given below: A≈2.85 cm2: area of neck ℓ≈5.7 cm: length of neck L≈25 cm: length of bottle V0≈1000 cm: volume of body
Figure 7.6: A one liter soda bottle.
Putting the length, L, of the bottle into (7.50) gives ν≈332 hertz . In American standard pitch (see Table 7.1), this is an E above middle C.
This is obviously wrong. If you have ever blown into your soda bottle, you know that the frequency of the lowest mode is much lower than that. The problem, of course, is that the soda bottle is not shaped anything like the tube. To determine the modes is a complicated three-dimensional problem. It turns out, however, that we can find the lowest mode to a decent approximation rather easily.
The idea is that in the lowest mode, the air in the neck of the bottle is moving rapidly, but in the body of the bottle, the air quickly spreads out so that it is not moving much at all. The idea of the Helmholtz approximation to try is to treat the air in the neck as a single chunk with mass ρAℓ,
and to treat the body as a spring, that contributes restoring force but no inertia (because the air is not moving much). Then all we must do is to compute the K of the “spring.” That is easy, using (7.38). In this case, dV=Adz,
so dp=−γpAdzV≈−γp0AdzV0
Table 7.1: American standard pitch (A440) — frequencies are in Hertz.
Note | ν | Note | ν | Note | ν |
---|---|---|---|---|---|
A | 880 | A | 440 | A | 220 |
G# | 831 | \boldsymbol{G #} | 415 | \boldsymbol{G #} | 208 |
G | 784 | G | 392 | G | 196 |
F# | 740 | \boldsymbol{F #} | 370 | \boldsymbol{F #} | 185 |
F | 698 | F | 349 | F | 175 |
E | 659 | E | 330 | E | 165 |
Eb | 622 | Eb | 311 | Eb | 156 |
D | 587 | D | 294 | D | 147 |
C# | 554 | C# | 277 | C# | 139 |
C | 523 | C | 262 | C | 131 |
B | 494 | B | 247 | B | 123 |
Bb |
466 |
Bb | 233 | Bb | 117 |
and F≈−γp0A2dzV0
or "K"=γp0A2V0.
Then using ω2=K/m, we expect ω=√γA2p0/V0ρAℓ=v√AℓV0.
For the soda bottle, (7.6), this gives ν≈118 hertz
or roughly a Bb below low C. This is just about right (see problem 7.5).
Corrections to Helmholtz
There are many possible corrections to (7.57) that might be considered. One is to include the so-called “end effect.” The point is that the velocity of the air in the lowest mode does not drop to zero immediately when you go past the ends of the neck. Thus the actual mass is somewhat larger than ρAℓ. The lore is that you can do better by replacing ℓ→ℓ+0.6r
where r is the radius of the neck.
Here we will discuss another correction that can be dealt with systematically using the methods of space translation invariance and local interactions. If the bottle has a long neck, it is probably not a good idea to treat the air in the neck as a solid mass. Furthermore, there is a simple alternative. A better analogy for the neck is a massive spring with Kℓ=γAp0. Because the neck is a space translation invariant, essentially one-dimensional system, we expect a displacement of the form ycosωzv
in the neck, where z=0 is the open end and y is the displacement of the air at z=0. Thus, where the neck attaches to the body, the displacement is ycosωℓv.
The force at this point from the compression of the air in the neck is (from (7.45)) Fneck =−γAp0∂ψ∂z=γAp0ωvysinωℓv.
This must be the negative of the force from the air in the body, from (7.39), −Fbody =γA2p0V0ycosωℓ/v,
or ωV0Avtanωℓv=1.
You will explore the consequences of this in problem 7.5.
This analysis does not distinguish between the area of the top and bottom of the neck. Perhaps the area at the bottom is more appropriate. What matters is the area at the bottom that determines the force per unit area where the wave in the neck matches onto the body.
___________________
3See, for example, Halliday and Resnick.