[Contents]
Copyright © 2021 jsd

Mach Arrival Times and Angles
John Denker

We consider the times and angles you perceive when a source changes position faster than the speed of sound. The same analysis applies when a source changes position faster than the speed of light, which actually can happen, as discussed below.

*   Contents

1 Setup and Notation

t ‍:= ‍time at which the signal was emitted            
f ‍:= ‍time of flight, from emission to arrival        
c ‍:= ‍speed of sound                                  
v ‍:= ‍speed of emitter (assumed supersonic)            
a ‍:= ‍actual arrival time                             
θ ‍:= ‍angle from which the signal arrives        
x ‍:= ‍observer’s location                             
y ‍:= ‍line of sight to the point of closest approach
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(1)

The time t is measured relative to the point of closest approach, so for the situation shown in the diagram, t is less than zero.

mach-angles-unscaled
Figure ‍1: Mach Angle Construction

2 Key Relationships

At all times:

c2f2 ‍= ‍y2 + v2t2             
f ‍= ‍
y2/c2 + t2v2/c2
             
a ‍= ‍t + f                           
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(2)

sin(θ) ‍= ‍vt / cf              
cos(θ) ‍= ‍y / cf               
tan(θ) ‍= ‍vt / y               
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(3)

See also the scaled versions in section ‍5.

3 Solutions

It is now straightforward to solve for the arrival time a and angle θ (both which you can observe directly) as functions of emission time t (which you can’t). The results are shown in figure ‍2. For this example, we assume y=1, c=1, and v=2.

mach-emission
Figure ‍2: Mach Arrival Time and Angle versus Emission Time

Here’s the graph of what you actually observe: Arrival angles as functions of arrival time. For each time there are two angles.

mach-arrival
Figure ‍3: Mach Arrival Angle versus Arrival Time

When there are two signals arriving at the same time, it is tricky to perceive the angle of either of them using just your two ears. It may be possible to untangle them using frequency information, due to the Doppler shift and/or dispersion. You could easily untangle them using a phased-array passive sonar receiver.

Things are even more interesting when you look at a source that changes its position faster than the speed of light. Your eyes contain millions of photoreceptors, so it is possible to perceive both of the arriving signals. You can’t have a source that physically moves faster than the speed of light, but you can have a long array of sub-sources that light up in turn, by pre-arrangement, creating a pattern that moves faster than the speed of light. As a particularly simple example, you can use to grids to create a moireé pattern that moves faster than the speed of light, even though neither grid is moving very fast, as shown in figure ‍4.

moire
Figure ‍4: Moiré Pattern

4 Special Cases

At the time of onset, when the first signal arrives, [e x b] is a right triangle, and β = θ. We know it’s a right triangle because the line [x b] is tangent to the wavefront, as suggested by the circular arc in the diagram. This can be used as a check on other ways of calculating the minimum time of flight f. Also it tells us that

sin(θ) ‍= ‍cf/vf ‍ ‍ ‍ ‍at onset         
  ‍= ‍c/v
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(4)

which is the famous formula for the Mach angle.

Also, the emission time of the onset signal is:

vt/y ‍= ‍tan(θ) ‍ ‍ ‍ ‍at onset         
t ‍= ‍(y/v) tan(θ) ‍ ‍ ‍ ‍at onset         
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(5)

In the special case where v/c = 2, at onset the triangle [e x b] is a 30-60-90 triangle. This is the case shown in the diagram.

At some later time, we might be tempted to consider the triangle [ex b′], but (1) it doesn’t tell us anything we don’t already know, (2) it bears no relationship to the physics of the Mach angle, and (3) it is not a right triangle.

5 Scaling

It is good practice to think about scaling, to elucidate the structure of a problem, and/or to reduce the number of variables we need to keep track of. The general topic of scaling is discussed in reference ‍1.

For the Mach arrivals problem, all the speeds scale in proportion to c, and all the times scale in proportion to y/c. If we replace v/c by m (the Mach number) and replace y/c by τ (the time of flight from the point of closest approach), then c drops out of the equations entirely. This reduces the number of variables involved. All the legs in figure ‍5 have dimensions of time (not distance).

mach-angles-scaled
Figure ‍5: Scaled Mach Angle Construction

The key equations become:

f2 ‍= ‍τ2 + m2t2             
f ‍= ‍
τ2 + t2m2
             
a ‍= ‍t + f                           
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(6)

sin(θ) ‍= ‍mt / f              
cos(θ) ‍= ‍τ / f               
tan(θ) ‍= ‍mt / τ               
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(7)

6 Inverse Mapping

Equation ‍6 gives us the arrival time a as a function of the emission time t. Sometimes we want the reverse. That requires a bit more work.

at ‍= ‍
τ2 + t2m2
             
a2 −2at + t2 ‍= ‍τ2 + t2m2             
0 ‍= ‍t2 (m2 − 1) + 2at + τ2 − a2          
t ‍= ‍
a ± 
m2a2 − (m2−1)τ2
(m2−1)
       
 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(8)

7 Caustics and Shocks

The black curve in figure ‍6 doesn’t tell us anything new. It is the same as the blue curve in figure ‍2; we have just transposed the axes.

mach-caustic
Figure ‍6: Caustic = Concentration of Pops

Now suppose that the source emits a series of pops, evenly spaced in time, finely spaced along the path of the source. This is shown by the red dots in figure ‍6. The question is, at what time do these pops arrive at the observation point? So let’s look at the spacing of the pops in the arrival-time direction.

You can see that a whole bunch of pops are conentrated at the onset time aon. More generally, for all times we can calculate the density of pops as:

ρ ‍= ‍



de
da



 ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍ ‍(9)

where the sum runs over both branches (upper and lower) of the black curve.

The density is shown by the green curve in figure ‍6. It diverges to infinity for times near onset. It scales like ρ ∝ 1/√(aaon). This sort of accumulation of events is called a caustic.

This means that you observe a very loud boom at onset. You can hear all the other pops, long after onset, but they will be not nearly so loud.

What’s more, the wave equation for sound is slightly nonlinear, and the boom is so loud that it gives rise to a shock, which makes the boom even sharper.

The analysis so far is more-or-less correct if the pops are shaped like delta functions, such as might be produced by a string of firecrackers. It also works for the optical case, if the source is emitting photons incoherently. In contrast, for an actual sonic boom, we would need to account for the fact that the source is emitting a series of delta-prime functions. Accumulating such things is a tricky business.

sr-71-ground-pressure-signature
Figure ‍7: SR 71 Ground Pressure Signature

Not all caustics produce shocks. For example, the bright flickering bands that you see on the bottom of a swimming pool are caustics caused by ripples on the surface. Als the explanation of how a rainbow is formed involves caustics. There is no nonlinearity in these optical examples.

8 References

1.
John Denker,
“Scaling Laws”
www.av8n.com/physics/scaling.htm
[Contents]
Copyright © 2021 jsd