General reconstruction of elastic strain fields from their Longitudinal Ray Transform

Chris Wensrich,Sean Holman,William Lionheart,Matias Courdurier,Anna Polyakova,Ivan Svetov,Ty Doubikin
2024-08-07
Abstract:We develop an algorithm for reconstruction of elastic strain fields from their Longitudinal Ray Transform (LRT) in either two or three dimensions. In general, the LRT only determines the solenoidal part of a symmetric tensor field, but elastic strain fields additionally satisfy mechanical equilibrium, an extra condition that allows for full reconstruction in many cases. Our method provides full reconstruction for general elastic strain fields in connected objects whose boundary only contains one component, while previous results included other requirements such as no residual stress, or zero boundary traction. This allows for full reconstruction in energy resolved neutron transmission imaging for simple objects. Along the way, we prove that the LRT of a potential rank-2 tensor restricted to a bounded set determines the potential on the boundary of the set up to infinitesimal rigid motions on each component of the boundary. The method is demonstrated with numerical examples in two dimensions.
Numerical Analysis,Applied Physics
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the problem of reconstructing the elastic strain field from the Longitudinal Ray Transform (LRT) of the elastic strain field. Specifically, the author has developed an algorithm that can reconstruct the elastic strain field from LRT in two - or three - dimensional space. #### Main problem background 1. **Energy - resolved neutron transmission imaging**: - In energy - resolved neutron transmission imaging, researchers need to reconstruct non - scalar spatial distribution information from projection images. - Bragg - edge strain tomography is a typical example, where the goal is to reconstruct the triaxial elastic strain field (a second - order tensor) within a sample from a series of neutron - transmission - based strain images. 2. **Limitations of existing methods**: - Previous reconstruction methods usually require that the sample has no residual stress or zero boundary traction. - These limitations make the application range of existing methods limited and unable to handle more complex situations. #### Main contributions of the paper 1. **Complete reconstruction without additional conditions**: - The method proposed in this paper can perform complete reconstruction of the elastic strain field satisfying mechanical equilibrium without imposing additional conditions (such as no residual stress or zero boundary traction). - This method is applicable to connected objects whose boundaries contain only one component, thus expanding the scope of application. 2. **Theoretical proof and numerical verification**: - The author has proved that the potential second - order tensor of LRT on a bounded set can determine the potential on the boundary of this set, up to infinitesimal rigid - body motions on each boundary component. - The effectiveness of this method in the two - dimensional case has been verified by numerical examples. 3. **Key formulas and theorems**: - Inversion formulas for reconstructing the elastic strain field have been proposed, such as equations (6) and (8), which utilize LRT and Helmholtz decomposition. - Theorem 1 shows that LRT can uniquely determine the value of the potential part on the boundary, up to infinitesimal rigid - body motions on each boundary component. 4. **Practical applications**: - This method is particularly suitable for energy - resolved neutron transmission imaging of simple objects and can provide complete strain - field reconstruction without requiring additional assumptions on the sample. In summary, this paper has solved the problem of how to reconstruct the elastic strain field from LRT and has broken through the strict limitations of previous methods on sample conditions, expanding the application range of this technology.