Researcher profile

William R. B. Lionheart

William R. B. Lionheart contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Recovering the second moment of the strain distribution from neutron Bragg edge data

Point by point strain scanning is often used to map the residual stress (strain) in engineering materials and components. However, the gauge volume and hence spatial resolution is limited by the beam defining apertures and can be anisotropic for very low and high diffraction (scattering) angles. Alternatively, wavelength resolved neutron transmission imaging has a potential to retrieve information tomographically about residual strain induced within materials through measurement in transmission of Bragg edges - crystallographic fingerprints whose locations and shapes depend on microstructure and strain distribution. In such a case the spatial resolution is determined by the geometrical blurring of the measurement setup and the detector point spread function. Mathematically, reconstruction of strain tensor field is described by the longitudinal ray transform; this transform has a non-trivial null-space, making direct inversion impossible. A combination of the longitudinal ray transform with physical constraints was used to reconstruct strain tensor fields in convex objects. To relax physical constraints and generalise reconstruction, a recently introduced concept of histogram tomography can be employed. Histogram tomography relies on our ability to resolve the distribution of strain in the beam direction, as we discuss in the paper. More specifically, Bragg edge strain tomography requires extraction of the second moment (variance about zero) of the strain distribution which has not yet been demonstrated in practice. In this paper we verify experimentally that the second moment can be reliably measured for a previously well characterised aluminium ring and plug sample. We compare experimental measurements against numerical calculation and further support our conclusions by rigorous uncertainty quantification of the estimated mean and variance of the strain distribution.

preprint2021arXiv

Crystalline phase discriminating neutron tomography using advanced reconstruction methods

Time-of-flight neutron imaging offers complementary attenuation contrast to X-ray computed tomography (CT), coupled with the ability to extract additional information from the variation in attenuation as a function of neutron energy (time of flight) at every point (voxel) in the image. In particular Bragg edge positions provide crystallographic information and therefore enable the identification of crystalline phases directly. Here we demonstrate Bragg edge tomography with high spatial and spectral resolution. We propose a new iterative tomographic reconstruction method with a tailored regularisation term to achieve high quality reconstruction from low-count data, where conventional filtered back-projection (FBP) fails. The regularisation acts in a separated mode for spatial and spectral dimensions and favours characteristic piece-wise constant and piece-wise smooth behaviour in the respective dimensions. The proposed method is compared against FBP and a state-of-the-art regulariser for multi-channel tomography on a multi-material phantom. The proposed new regulariser which accommodates specific image properties outperforms both conventional and state-of-the-art methods and therefore facilitates Bragg edge fitting at the voxel level. The proposed method requires significantly shorter exposure to retrieve features of interest. This in turn facilitates more efficient usage of expensive neutron beamline time and enables the full utilisation of state-of-the-art high resolution detectors.

preprint2020arXiv

Lipschitz stability at the boundary for time-harmonic diffuse optical tomography

We study the inverse problem in Optical Tomography of determining the optical properties of a medium $Ω\subset\mathbb{R}^n$, with $n\geq 3$, under the so-called diffusion approximation. We consider the time-harmonic case where $Ω$ is probed with an input field that is modulated with a fixed harmonic frequency $ω=\frac{k}{c}$, where $c$ is the speed of light and $k$ is the wave number. We prove a result of Lipschitz stability of the absorption coefficient $μ_a$ at the boundary $\partialΩ$ in terms of the measurements in the case when the scattering coefficient $μ_s$ is assumed to be known and $k$ belongs to certain intervals depending on some a-priori bounds on $μ_a$, $μ_s$.

preprint2020arXiv

Scanning electron diffraction tomography of strain

Strain engineering is used to obtain desirable materials properties in a range of modern technologies. Direct nanoscale measurement of the three-dimensional strain tensor field within these materials has however been limited by a lack of suitable experimental techniques and data analysis tools. Scanning electron diffraction has emerged as a powerful tool for obtaining two-dimensional maps of strain components perpendicular to the incident electron beam direction. Extension of this method to recover the full three-dimensional strain tensor field has been restricted though by the absence of a formal framework for tensor tomography using such data. Here, we show that it is possible to reconstruct the full non-symmetric strain tensor field as the solution to an ill-posed tensor tomography inverse problem. We then demonstrate the properties of this tomography problem both analytically and computationally, highlighting why incorporating precession to perform scanning precession electron diffraction may be important. We establish a general framework for non-symmetric tensor tomography and demonstrate computationally its applicability for achieving strain tomography with scanning precession electron diffraction data.

preprint2019arXiv

Histogram Tomography

In many tomographic imaging problems the data consist of integrals along lines or curves. Increasingly we encounter "rich tomography" problems where the quantity imaged is higher dimensional than a scalar per voxel, including vectors tensors and functions. The data can also be higher dimensional and in many cases consists of a one or two dimensional spectrum for each ray. In many such cases the data contain not just integrals along rays but the distribution of values along the ray. If this is discretized into bins we can think of this as a histogram. In this paper we introduce the concept of "histogram tomography". For scalar problems with histogram data this holds the possibility of reconstruction with fewer rays. In vector and tensor problems it holds the promise of reconstruction of images that are in the null space of related integral transforms. For scalar histogram tomography problems we show how bins in the histogram correspond to reconstructing level sets of function, while moments of the distribution are the x-ray transform of powers of the unknown function. In the vector case we give a reconstruction procedure for potential components of the field. We demonstrate how the histogram longitudinal ray transform data can be extracted from Bragg edge neutron spectral data and hence, using moments, a non-linear system of partial differential equations derived for the strain tensor. In x-ray diffraction tomography of strain the transverse ray transform can be deduced from the diffraction pattern the full histogram transverse ray transform cannot. We give an explicit example of distributions of strain along a line that produce the same diffraction pattern, and characterize the null space of the relevant transform.