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=L1−L2c,\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 (L1−L2)/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 t→∞t \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 t→∞t \to \infty: 0.5 pts

Part ii) (3 points)

As t→∞t \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=f0∣1±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 ∣vr∣→v|v_r| \to v:

f2=f0M−1(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: 1f1−1f2=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=2f1f2f2−f1,M=f2+f1f2−f1.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: f1≈320 Hzf_1 \approx 320\,\text{Hz} and f2≈1600 Hzf_2 \approx 1600\,\text{Hz}. Therefore

f0=2⋅320⋅16001280=800 Hz,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/c∣f = f_0/|1 \pm v_r/c|: 0.5 pts
  • Asymptotic upper-ridge value f2=f0/(M−1)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)/(f2−f1)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=2f1f2f2−f1=2⋅320⋅16001280=800 Hz.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=800 Hzf_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∣=∣OA∣2∣JA∣=ctM,|BA| = \frac{|OA|^2}{|JA|} = \frac{ct}{M}, d=∣OA∣⋅∣JO∣∣JA∣=ctM2−1M.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

τ=∣BA∣v+dc−t,\tau = \frac{|BA|}{v} + \frac{d}{c} - t, τ=tM2+tM2−1M−t.\tau = \frac{t}{M^2} + \frac{t\sqrt{M^2-1}}{M} - t.

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

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

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

1−1−1/2.25=1−0.5556≈0.255,1 - \sqrt{1 - 1/2.25} = 1 - \sqrt{0.5556} \approx 0.255, d≈340⋅1.10.255≈1.5 km.d \approx \frac{340\cdot 1.1}{0.255} \approx 1.5\,\text{km}.

Grading

  • Computing f0=2f1f2/(f2−f1)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