advanced28 min studyLesson 167 in Physics

Longitudinal and Transverse Elastic Waves

Derive the two bulk wave speeds in an isotropic solid and separate polarization, propagation, and interference.

On this page

What you’ll understand

  • Obtain the linear elastodynamic equation from momentum balance.
  • Calculate longitudinal and transverse bulk-wave speeds.
  • Distinguish a bulk solid wave from a slender-rod wave and from a scalar interference display.

An elastic solid can resist both compression and shear, so a disturbance has more than one way to travel. In a longitudinal wave, the displacement points along the direction of propagation. In a transverse wave, it points perpendicular to that direction. These statements describe particle motion near each material point, not the route taken by a marked chunk of material across the specimen. The difference in restoring forces produces two distinct bulk wave speeds.

Momentum balance becomes a field equation

Take a homogeneous isotropic solid initially at rest, with small displacement and small strain. The divergence of the stress tensor is the internal force per volume. Newton’s law equates it to density times displacement acceleration, after any body force has been included. Insert the isotropic elastic constitutive law and the symmetric displacement-gradient strain. Because the elastic constants are spatially uniform, the derivatives collect into a vector Laplacian and the gradient of displacement divergence.

The resulting equation already shows two mechanisms. The shear modulus multiplies spatial bending of the displacement field, while the other term responds to local compression. Density is mass per volume, not a modulus, and supplies inertia. No damping, material anisotropy, or frequency-dependent constitutive response is included here. The continuum approximation also presumes wavelengths large enough compared with the material microstructure. These assumptions matter when comparing the simple prediction with ultrasound in a textured, layered, or strongly dissipative sample.

ρ ∂t2u=μ∇2u+(λ+μ)∇(∇⋅u)\rho\,\partial_t^2\mathbf u=\mu\nabla^2\mathbf u+(\lambda+\mu)\nabla(\nabla\cdot\mathbf u)
This is linear elastodynamics without body forces in a homogeneous isotropic medium. The Lamé constant λ here is not a wavelength.↔ Scroll wide equations horizontally.

Choose a polarization before solving for speed

Consider a plane wave with a wavevector pointing in a fixed direction. If its displacement is parallel to that direction, displacement divergence is nonzero and both restoring terms contribute. Their sum is the longitudinal modulus, the first Lamé constant plus twice the shear modulus. For displacement perpendicular to the wavevector, divergence vanishes. Only the shear contribution remains. There are two independent transverse polarizations in three spatial dimensions, with the same speed in an isotropic solid.

Positive bulk and shear moduli ensure that the squared speeds are positive. The longitudinal modulus can also be written as the bulk modulus plus four thirds of the shear modulus; that form makes the connection with compression and distortion especially clear. Do not substitute Young’s modulus into this bulk longitudinal formula. A slender rod with free lateral surfaces can contract sideways, and its long-wavelength extensional speed instead approaches the square root of Young’s modulus divided by density. The geometry and boundary constraints determine which formula applies.

cL=λ+2μρ=K+4μ/3ρ,cT=μρc_L=\sqrt{\frac{\lambda+2\mu}{\rho}}=\sqrt{\frac{K+4\mu/3}{\rho}},\qquad c_T=\sqrt{\frac{\mu}{\rho}}
These are bulk longitudinal and transverse speeds in a homogeneous three-dimensional isotropic solid. Both ratios under the square roots have units m²/s².↔ Scroll wide equations horizontally.

What a moving wave profile does and does not show

In this ideal bulk medium each polarization has a linear frequency–wavenumber relation, so its phase and group speeds coincide. A localized pulse therefore does not disperse through this particular mechanism. Boundaries can reflect waves and convert one polarization into another; a finite body can also support surface and guided modes that require additional boundary calculations. None of those effects is determined by adding two scalar sine curves. An ordinary fluid has no static shear modulus in this elastic model and consequently has no corresponding propagating bulk shear wave.

In the waves lab, set the amplitude ratio to one and predict the sum for zero phase difference and for a half-turn phase difference. The curves reinforce in the first case and cancel in the second because both displayed components have the same scalar coordinate, frequency, and wavenumber. This is an interference experiment, not a comparison of longitudinal and transverse polarizations. Two orthogonally polarized vector displacements do not cancel simply because their scalar amplitudes have opposite signs. The controls do not change density, elastic moduli, or either physical wave speed derived above.

Two pulse arrivals through an ideal bulk sample

Use ρ = 2,500 kg/m³, μ = 10 GPa, and λ = 20 GPa for an illustrative isotropic solid. Find the bulk speeds and the arrival times over a straight 100 m path, ignoring reflections and attenuation.

  1. For the shear wave, μ/ρ = 10 × 10⁹/2,500 = 4,000,000 m²/s², so cT = 2,000 m/s.

  2. For the longitudinal wave, (λ + 2μ)/ρ = 40 × 10⁹/2,500 = 16,000,000 m²/s², so cL = 4,000 m/s.

  3. The travel times are 100/4,000 = 0.025 s and 100/2,000 = 0.050 s. The longitudinal pulse arrives 0.025 s earlier.

The result

The ideal arrivals are 25 ms for the longitudinal pulse and 50 ms for the transverse pulse. A real measurement must also identify its path and mode.

Explore the model

Interactive lab

Interactive lab · waves

Test interference without changing wave speed

Set the amplitude ratio to one. Predict reinforcement at 0° phase difference and cancellation at 180°, then test both. These are two scalar waves with fixed wavenumber and frequency, not distinct elastic polarizations.

Test interference without changing wave speedTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forPhase controls scalar interference here; the physical bulk speeds instead depend on the appropriate modulus divided by density.

Check your understanding

Why does the ideal ordinary fluid limit μ = 0 lack a propagating bulk transverse elastic wave?

Try it yourself

In the worked material, double μ while leaving λ and density unchanged. By what factors do the two bulk speeds change?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Compare the squared speeds first. Initially λ = 2μ, but hold the original λ fixed during the change.

Compare with the solution

The transverse squared speed doubles, so cT becomes √2 times its old value, about 2,828 m/s. The longitudinal numerator changes from 20 + 2(10) = 40 GPa to 20 + 2(20) = 60 GPa. Its speed therefore gains a factor √(60/40) = √1.5, becoming about 4,899 m/s. The two factors differ because shear stiffness contributes to both modes.

Take it with you

  • Longitudinal and transverse label displacement relative to propagation direction.
  • Bulk isotropic speeds are √((λ + 2μ)/ρ) and √(μ/ρ), within the stated linear continuum assumptions.
  • Changing geometry or boundary conditions can change the applicable wave model; scalar interference does not model polarization conversion.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.