Asymptotics of the resonances for a continuously stratified layer

Ultrasound wave propagation in a nonhomogeneous linearly elastic layer of constant thickness is considered. The resonances for the corresponding acoustic propagator are studied. It is shown that the distribution of the resonances depends on the smoothness of the coefficients. Namely, if the coefficients have jump discontinuities at the boundaries, then the resonances are asymptotically distributed along a straight line parallel to the real axis on the unphysical sheet of the complex frequency plane. In the contrary, if the coefficients are continuous, then it is shown that the resonances are asymptotically distributed along a logarithmic curve. The spacing between two successive resonances turns out to be sensitive to articular cartilage degeneration. The application of the obtained results to ultrasound testing of articular cartilage is discussed.

💡 Research Summary

The paper investigates the high‑frequency asymptotic distribution of resonances (also called scattering poles) for a one‑dimensional acoustic propagator that models an ultrasound wave traveling through a linearly elastic layer of constant thickness whose material properties vary with depth. The authors formulate the problem as a Sturm–Liouville eigenvalue problem for the complex frequency ω, introduce the transfer matrix of the layer, and define resonances as the complex ω‑values at which the (1,1) entry of the transfer matrix vanishes on the non‑physical sheet of the complex ω‑plane.

Two distinct regimes are examined, depending on the smoothness of the density ρ(x) and elastic modulus E(x) at the two interfaces (x = 0 and x = h).

1. Jump discontinuities at the boundaries – When ρ and E have finite jumps at the interfaces, a WKB‑type high‑frequency approximation yields a resonance condition of the form
   Im k(ω)·h ≈ π n, Re k(ω) ≈ constant,
   where k(ω)≈ω/c_eff and c_eff is an effective wave speed determined by the impedance contrast. Consequently the resonances lie asymptotically on a straight line parallel to the real axis on the unphysical sheet, with a nearly constant spacing ΔIm ω ≈ π/(c_eff h). This reflects the fact that a sharp impedance mismatch forces the wave to be almost completely reflected, producing high‑Q quasi‑standing‑wave modes whose decay rates are almost identical.

2. Continuous coefficients – If ρ(x) and E(x) are continuous (and sufficiently smooth) across the boundaries, the transfer matrix acquires a logarithmic factor in its high‑frequency expansion:
   M(ω) ≈ exp(i ω τ)·(log ω)^α·