The reflection coefficient of a fractional reflector

Laurent Demanet,Olivier Lafitte

2023-05-06

Abstract:This paper considers the question of characterizing the behavior of waves reflected by a fractional singularity of the wave speed profile, i.e., of the form \[ c(x_1, x_2, x_3) = c_0 \left(1 + \left( \frac{x_1}{\ell}\right)_{+}^\alpha \right)^{-1/2}, \] for $\alpha > 0$ not necessarily integer. We first focus on the case of one spatial dimension and a harmonic time dependence. We define the reflection coefficient $R$ from a limiting absorption principle. We provide an exact formula for $R$ in terms of the solution to a Volterra equation. We obtain the asymptotic limit of this coefficient in the large $\ell \omega / c_0$ regime as \[ R = \frac{\Gamma(\alpha + 1)}{(2 i)^{\alpha + 2}} \left( \frac{c_0}{\ell \omega} \right)^{\alpha} + \mbox{lower order terms.} \] The amplitude is proportional to $\omega^{-\alpha}$, and the phase rotation behavior is obtained from the $i^{-(\alpha+2)}$ factor. The proof method does not rely on representing the solution by special functions, since $\alpha > 0$ is general. In the multi-dimensional layered case, we obtain a similar result where the nondimensional variable $\ell \omega / c_0$ is modified to account for the angle of incidence. The asymptotic analysis now requires the waves to be non-glancing. The resulting reflection coefficient can now be interpreted as a Fourier multiplier of order $- \alpha$. In practice, the knowledge of the dependency of both the amplitude and the phase of $R$ on $\omega$ and $\alpha$ might be able to inform the kind of signal processing needed to characterize the fractional nature of reflectors, for instance in geophysics.

Analysis of PDEs,Classical Analysis and ODEs