Mapping corrosion depths along pipeline sections using guided-wave-based tomographic methods is a challenging task. Accurate defect sizing depends heavily on the precision of the forward model in guided wave tomography. This model is fitted to measured data using inversion techniques. This study evaluates the effectiveness of a recursive extrapolation scheme for tomography applications and full waveform inversion. It employs a table-driven approach, with precomputed extrapolation operators stored across a spectrum of wavenumbers. This enables fast modelling for extensive pipe sections, approaching the speed of ray tracing while accurately handling complex velocity models within the full frequency band. This ensures an accurate representation of diffraction phenomena. The study examines the assumptions underlying the extrapolation approach, namely, the negligible reflection and conversion of modes at defects. In our tomography approach, we intend to use multiple wave modes— (Formula presented.), (Formula presented.), and (Formula presented.) —and helical paths. The acoustic extrapolation method is validated through numerical studies for different wave modes, solving the 3D elastodynamic wave equation. Comparison with an experimentally measured single-mode wavefield from an aluminium plate with an artificial defect reveals good agreement.

Over the last 15 years, literature on nondestructive testing has shown that the generation of higher harmonics and nonlinear mixing of waves could be used to obtain the nonlinearity parameters of an elastic medium and thereby gather information about its state, e.g., aging and fatigue. To design ultrasound measurement setups based on these phenomena, efficient numerical modeling tools are needed. In this paper, the iterative nonlinear contrast source method for numerical modeling of nonlinear acoustic waves is extended to the one-dimensional elastic case. In particular, nonlinear mixing of two collinear bulk waves (one compressional, one shear) in a homogeneous, isotropic medium is considered, taking into account its third-order elastic constants ( A,  B,  and  C). The obtained results for nonlinear propagation are in good agreement with a benchmark solution based on the modified Burgers equation. The results for the resonant waves that are caused by the one-way and two-way mixing of primary waves are in quantitative agreement with the results in the literature [Chen, Tang, Zhao, Jacobs, and Qu, J. Acoust. Soc. Am. 136(5), 2389-2404 (2014)]. The contrast source approach allows the identification of the propagating and evanescent components of the scattered wavefield in the wavenumber-frequency domain, which provides physical insight into the mixing process and explains the propagation direction of the resonant wave.