A non-reflecting wave equation through directional wave-field suppression and its finite difference implementation

Abstract

The acoustic wave equation describes wave propagation directly from basic physical laws, even in heterogeneous acoustic media. When numerically simulating waves with the wave equation, contrasts in the medium parameters automatically generate all scattering effects. For some applications - such as propagation analysis or certain wave-equation based imaging techniques - it is desirable to suppress these reflections, as we are only interested in the transmitted wave-field. To achieve this, a modification to the constitutive relations is proposed, yielding an extra term that suppresses waves with reference to a preferred direction. The scale-factor α of this extra term can either be interpreted as a penetration depth or as a typical decay time. This modified theory is implemented using a staggered-grid, time-domain finite difference scheme, where the acoustic Poynting-vector is used to estimate the local propagation direction of the wave-field. The method was successfully used to suppress reflections in media with bone tissue (medical ultrasound) and geophysical subsurface structures, while introducing only minor perturbations to the transmitted wave-field and a small increase in computation time.