Radiative transfer of acoustic waves in continuous complex media: Beyond the Helmholtz equation

TL;DR

This paper derives a radiative transfer equation for acoustic waves in complex media, highlighting the significant impact of operator terms on transport parameters, especially the mean-free path, across frequency regimes.

Contribution

It introduces a comprehensive radiative transfer model that includes operator contributions, extending beyond the traditional Helmholtz equation approach for acoustic wave scattering.

Findings

Operator terms significantly affect phase function and anisotropy factor at low frequencies.

Disregarding operator terms causes large errors in the transport mean-free path, up to 300%.

Results are validated through numerical wave simulations and Monte Carlo methods.

Abstract

Heterogeneity can be accounted for by a random potential in the wave equation. For acoustic waves in a fluid with fluctuations of both density and compressibility (as well as for electromagnetic waves in a medium with fluctuation of both permittivity and permeability) the random potential entails a scalar and an operator contribution. For simplicity, the latter is usually overlooked in multiple scattering theory: whatever the type of waves, this simplification amounts to considering the Helmholtz equation with a sound speed $c$ depending on position $\mathbf{r}$. In this work, a radiative transfer equation is derived from the wave equation, in order to study energy transport through a multiple scattering medium. In particular, the influence of the operator term on various transport parameters is studied, based on the diagrammatic approach of multiple scattering. Analytical results are obtained for fundamental quantities of transport theory such as the transport mean-free path $\ell^*$, scattering phase function $f$ and anisotropy factor $g$. Discarding the operator term in the wave equation is shown to have a significant impact on $f$ and $g$, yet limited to the low-frequency regime i.e., when the correlation length of the disorder $\ell_c$ is smaller than or comparable to the wavelength $\lambda$. More surprisingly, discarding the operator part has a significant impact on the transport mean-free path $\ell^*$ whatever the frequency regime. When the scalar and operator terms have identical amplitudes, the discrepancy on the transport mean-free path is around $300\,\%$ in the low-frequency regime, and still above $30\,\%$ for $\ell_c/\lambda=10^3$ no matter how weak fluctuations of the disorder are. Analytical results are supported by numerical simulations of the wave equation and Monte Carlo simulations.