8. Jet Sound 8 pts

Waves · Doppler effect, Sound propagation, Supersonic flow

Decode a supersonic flyby spectrogram to extract the Mach number and closest-approach distance from the Doppler asymptotes and the sonic-boom delay.

Solution by Teo Kai Wen and Jaan Kalda.

Part i) (1.5 points)

The spectrogram shows four qualitative features.

(a) Silence before t0t_0. Since sound travels at a finite speed, it takes some time to hear the sound (this feature would also hold for a non-supersonic plane).

(b) Broadband spike at t0t_0. An entire interval of past emissions reaches the observer simultaneously, forming a shockwave front — the “boom” — that registers as a short, strong pressure peak at the sensor and represents the superposition of a wide range of frequencies.

(c) Two simultaneous tones for t>t0t > t_0. For each reception time t>t0t > t_0 (let us set t=0t = 0), two distinct past emission events 1 and 2 contribute. Suppose the sound from event 1, emitted at time t1-t_1 from distance L1L_1, reaches us now. Since the jet is supersonic, there is a second emission — event 2 at a later time t2-t_2, distance L2L_2 from us, separated from event 1 by distance L0L_0 along the trajectory — such that

L0v=L1L2c,\frac{L_0}{v} = \frac{L_1 - L_2}{c},

i.e., the jet’s traversal time from 1 to 2 equals the sound delay (L1L2)/c(L_1 - L_2)/c. Both arrivals are Doppler-shifted, but to different degrees: the sound coming from behind (event 1, jet far upstream and approaching the observer) is blue-shifted (upper ridge), while the sound coming from ahead (event 2, jet receding past closest approach) is red-shifted (lower ridge).

(d) Ridges decrease in time. Both ridges diverge at the boom (the broadband spike at t0t_0) and decrease monotonically afterward, tending toward finite asymptotic values as tt \to \infty.

Graph: Geometry sketch: the supersonic jet flies along a horizontal line; observer "obs" is below. The earlier event 1 (red, upper) is upstream and emits sound that travels distance L_1 to the observer; the later event 2 (lower) is downstream and emits sound that travels distance L_2. Distance L_0 = vt along the trajectory separates the two emission points.

Grading

  • Silence before t0t_0, explained as the sound not reached the receiver yet: 0.2 pts
  • Sonic boom at t0t_0, explained as the shock-wave front producing a short broadband pressure peak: 0.5 pts
  • Associating the changes in frequency with the Doppler effect: 0.2 pts
  • Explanation of two simultaneous tones for t>t0t > t_0, with two emission events reaching the observer at the same time: 0.6 pts
  • Explanation of ridges diverging at the boom and tending asymptotically to finite values as tt \to \infty: 0.5 pts

Part ii) (3 points)

As tt \to \infty, the jet’s velocity becomes nearly aligned with the line of sight, so the radial velocity (component along the line from source to observer) approaches ±v\pm v. The classical Doppler formula for a source moving radially at speed vrv_r relative to a stationary observer is

f=f01±vr/c,f = \frac{f_0}{|1 \pm v_r/c|},

with the - sign for an approaching source (vr>0v_r > 0, source moving toward observer) and ++ for a receding source. Applied to the two asymptotic limits with vrv|v_r| \to v:

f2=f0M1(upper ridge, approaching),f_2 = \frac{f_0}{M-1} \quad \text{(upper ridge, approaching)}, f1=f0M+1(lower ridge, receding).f_1 = \frac{f_0}{M+1} \quad \text{(lower ridge, receding)}.

Combining: 1f11f2=2f0\frac{1}{f_1} - \frac{1}{f_2} = \frac{2}{f_0} and 1f1+1f2=2Mf0\frac{1}{f_1} + \frac{1}{f_2} = \frac{2M}{f_0}, so

f0=2f1f2f2f1,M=f2+f1f2f1.f_0 = \frac{2f_1f_2}{f_2 - f_1}, \qquad M = \frac{f_2 + f_1}{f_2 - f_1}.

Reading from the spectrogram: f1320Hzf_1 \approx 320\,\text{Hz} and f21600Hzf_2 \approx 1600\,\text{Hz}. Therefore

f0=232016001280=800Hz,f_0 = \frac{2\cdot 320\cdot 1600}{1280} = 800\,\text{Hz}, M=19201280=1.5.M = \frac{1920}{1280} = 1.5.

Grading

  • Recognising that the radial velocity tends to ±v\pm v in the asymptotic limits: 0.5 pts
  • Correct Doppler formula f=f0/1±vr/cf = f_0/|1 \pm v_r/c|: 0.5 pts
  • Asymptotic upper-ridge value f2=f0/(M1)f_2 = f_0/(M-1): 0.2 pts
  • Asymptotic lower-ridge value f1=f0/(M+1)f_1 = f_0/(M+1): 0.2 pts
  • Reading f1[310,350]Hzf_1 \in [310, 350]\,\text{Hz} and f2[1500,1700]Hzf_2 \in [1500, 1700]\,\text{Hz} from the spectrogram (0.2p each): 0.4 pts
  • Correct symbolic form M=(f2+f1)/(f2f1)M = (f_2 + f_1)/(f_2 - f_1): 0.8 pts
  • Numerical answer M[1.35,1.65]M \in [1.35, 1.65]: 0.4 pts

Part iii) (3 points)

The closest-approach signal is the sound emitted when the jet’s radial velocity is zero (jet directly overhead BB). At that emission, vr=0v_r = 0, so the received frequency equals f0f_0. From part (ii)‘s asymptotic relations we have

f0=2f1f2f2f1=232016001280=800Hz.f_0 = \frac{2f_1f_2}{f_2 - f_1} = \frac{2\cdot 320\cdot 1600}{1280} = 800\,\text{Hz}.

The observer hears this frequency on the lower ridge (the receding branch crosses f0f_0 as the jet passes overhead). Read t1t_1 from the spectrogram as the time when the lower ridge crosses f0=800Hzf_0 = 800\,\text{Hz}.

Set t=0t = 0 at the moment the jet emits the sound that the observer eventually hears as the boom; call this position AA. Let JJ denote the jet’s position when the observer hears the boom, and OO the observer. Then

AJ=vt,AO=ct,|AJ| = vt, \qquad |AO| = ct,

and the boom-tangency condition fixes AOJ=π/2\angle AOJ = \pi/2 (the wavefront from AA just grazes the observer). In particular, sin(AJO)=AO/AJ=1/M\sin(\angle AJO) = |AO|/|AJ| = 1/M, identifying AJO\angle AJO as the Mach angle α\alpha.

Graph: The supersonic jet at speed v moves along trajectory from A (emission of the boom-causing wave) past B (closest-approach foot) to J (jet's position when boom is heard). Observer O is below, at perpendicular distance d = |OB|. Right triangles △JOA (right angle at O) and △OBA (right angle at B) share angle α at A; the Mach cone has half-angle α with sin α = 1/M.

Let BB be the foot of perpendicular from OO onto trajectory AJAJ, so OB=d|OB| = d (closest-approach distance). The right triangles JOA\triangle JOA and OBA\triangle OBA are similar (both right-angled, sharing the angle at AA). From this similarity:

BA=OA2JA=ctM,|BA| = \frac{|OA|^2}{|JA|} = \frac{ct}{M}, d=OAJOJA=ctM21M.d = \frac{|OA|\cdot|JO|}{|JA|} = \frac{ct\sqrt{M^2-1}}{M}.

The observer hears the boom at time tt and the closest-approach signal at time BA/v+d/c|BA|/v + d/c. The lag is

τ=BAv+dct,\tau = \frac{|BA|}{v} + \frac{d}{c} - t, τ=tM2+tM21Mt.\tau = \frac{t}{M^2} + \frac{t\sqrt{M^2-1}}{M} - t.

Eliminating tt via t=dM/(cM21)t = dM/(c\sqrt{M^2-1}) and simplifying:

τc=d(11M2),\tau c = d\left(1 - \sqrt{1 - M^{-2}}\right), d=τc11M2.\boxed{\,d = \frac{\tau c}{1 - \sqrt{1 - M^{-2}}}.\,}

Reading from the spectrogram: t03.3st_0 \approx 3.3\,\text{s} (boom onset) and t14.4st_1 \approx 4.4\,\text{s} (lower ridge crosses 800Hz800\,\text{Hz}), so τ1.1s\tau \approx 1.1\,\text{s}. With M=1.5M = 1.5:

111/2.25=10.55560.255,1 - \sqrt{1 - 1/2.25} = 1 - \sqrt{0.5556} \approx 0.255, d3401.10.2551.5km.d \approx \frac{340\cdot 1.1}{0.255} \approx 1.5\,\text{km}.

Grading

  • Computing f0=2f1f2/(f2f1)f_0 = 2f_1f_2/(f_2 - f_1) with f0[720,910]Hzf_0 \in [720, 910]\,\text{Hz} using the asymptotes from part (ii): 0.4 pts
  • Using Mach cone and the corresponding condition sinα=1/M\sin\alpha = 1/M: 0.3 pts
  • Recognising that at the point of closest approach, f=f0f = f_0 (no Doppler shift): 0.3 pts
  • Expressing the time t0t_0 between the emission (event A) and reception of the sonic boom in terms of dd: 0.5 pts
  • Expressing the time t1t_1 between event AA and reception of the wave emitted at BB in terms of dd: 1.0 pts
  • Reading t0t_0, t1t_1 from the spectrogram and computing τ[0.85,1.15]s\tau \in [0.85, 1.15]\,\text{s}: 0.2 pts
  • Numerical answer if approach correct d[1.0,1.7]kmd \in [1.0, 1.7]\,\text{km}: 0.3 pts