Direct inversion of the Longitudinal Ray Transform for 2D residual elastic strain fields

TL;DR

This paper introduces a novel method for directly reconstructing 2D elastic strain fields from neutron transmission data using a Helmholtz decomposition and tensor filtered back projection, validated through simulations and experiments.

Contribution

It develops a new inversion approach for the Longitudinal Ray Transform in 2D elastic systems, combining Helmholtz decomposition with tensor filtered back projection.

Findings

The method successfully reconstructs elastic strain fields from simulated and experimental data.

Applying scalar filtered back projection recovers the trace of the solenoidal strain component.

The approach provides physical insights into the reconstructed strain components.

Abstract

We examine the problem of Bragg-edge elastic strain tomography from energy resolved neutron transmission imaging. A new approach is developed for two-dimensional plane-stress and plane-strain systems whereby elastic strain can be reconstructed from its Longitudinal Ray Transform (LRT) as two parts of a Helmholtz decomposition based on the concept of an Airy stress potential. The solenoidal component of this decomposition is reconstructed using an inversion formula based on a tensor filtered back projection algorithm whereas the potential part can be recovered using either Hooke's law or a finite element model of the elastic system. The technique is demonstrated for two-dimensional plane-stress systems in both simulation, and on real experimental data. We also demonstrate that application of the standard scalar filtered back projection algorithm to the LRT in these systems recovers the trace of the solenoidal component of strain and we provide physical meaning for this quantity in the case of 2D plane-stress and plane-strain systems.