Elastic Waves, and Where the Tensors Lead

One master equation, one substitution

Recall the toolkit we have built: the inertia tensor for rotation, and the stress and strain tensors for deformation — all symmetric rank-2 machines. To find how a solid moves, we close the loop: take Cauchy's equation of motion \rho\,\ddot{u}_i = \partial_j\sigma_{ij} and substitute the isotropic Hooke's law \sigma_{ij}=\lambda\,\varepsilon_{kk}\delta_{ij}+2\mu\,\varepsilon_{ij} for the stress. Everything reduces to the displacement field alone.

\rho\,\ddot{\mathbf{u}} = (\lambda + \mu)\,\nabla(\nabla\!\cdot\mathbf{u}) + \mu\,\nabla^2\mathbf{u}

The Navier equation of elastodynamics — Cauchy's law with Hooke's law folded in, a closed wave equation for the displacement field.

Two waves hiding in one equation

The Navier equation looks like one equation, but it secretly contains two waves. Split the displacement (a Helmholtz decomposition) into a curl-free part, which compresses the material, and a divergence-free part, which shears it. Each part obeys an ordinary wave equation — with its own speed:

v_P = \sqrt{\dfrac{\lambda + 2\mu}{\rho}}, \qquad v_S = \sqrt{\dfrac{\mu}{\rho}}

The longitudinal (P) and transverse (S) elastic wave speeds: P compresses the material along its direction of travel, S shears it sideways.

Because \lambda and \mu are both positive, v_P > v_S always — the compression wave outruns the shear wave in every solid. Their ratio depends on nothing but Poisson's ratio:

\frac{v_P}{v_S} = \sqrt{\frac{\lambda + 2\mu}{\mu}} = \sqrt{\frac{2(1-\nu)}{1-2\nu}} > 1

The P-wave always leads the S-wave; the size of the lead is set entirely by Poisson's ratio \nu.

X-raying the Earth with two speeds

This is the foundation of seismology. An earthquake radiates both waves at once. The faster P-wave (Primary) arrives first, the S-wave (Secondary) later; the delay between them tells a seismograph how far away the quake struck. Timing these arrivals across a global network reconstructs the layered interior of a planet no drill could ever reach.

In honesty, the real Earth is layered, anisotropic and slightly dissipative, so the simple two-speed isotropic model is only the first approximation — but it is the one that opened the whole subject, and everything else is a refinement of it.

The whole track, assembled

Step back and see the arc. From the point mass we grew two theories of extended matter. Rigid rotation: the inertia tensor, its principal axes, Euler's equations, the precessing gyroscope. Deformation: the strain and stress tensors, Hooke's law, Cauchy's equation, and now elastic waves. The common thread is unmistakable — symmetric rank-2 tensors and the linear machines they define.

The reach of these ideas is enormous: spacecraft attitude control and reaction wheels (gyroscopes); the stability of a spinning rifle bullet and the tumble of a thrown wrench; why beams and bridges bend the way they do; the speed of sound in a steel rail; ultrasound imaging; and the seismology that maps our planet's core.

Where the tensors lead next

Let the shear modulus fall to zero and allow the material to flow, and continuum mechanics becomes fluid mechanics. Make the stress depend on the rate of strain rather than the strain itself — that is viscosity — and you arrive at the Navier–Stokes equations and the unsolved riddles of turbulence.

The stress tensor itself keeps generalizing. In electromagnetism the Maxwell stress tensor carries the momentum of the field; in relativity the stress–energy tensor is the very source of gravity in Einstein's equations. The rank-2 tensor you built for a steel bar is, mathematically, the same object that curves spacetime.