Transient tomography for defect detection

Lead Research Organisation: University of Leeds
Department Name: Applied Mathematics

Abstract

The ability to locate, estimate the size and shape, along with the properties of a hidden defect concealed in an exterior enclosure is of importance to the sectors of Security (landmine detection or detecting explosive material or illicit substances in vehicles stationary or in motion), Energy (detecting underground pipe blockages or, monitoring the structural integrity of a reactor while reducing human exposure in harmful and hostile environments), Health (detecting anomalies/ tumours/viruses), etc. Unlike similar tomographic techniques, e.g., electrical impedance/ resistive/capacitance tomography, which are stationary methods, time-dependent tomography uses transient data information to detect hidden defects which may be static or moving in time. In comparison with its stationary boundary potential + current flux formulation, the transient boundary temperature + heat flux formulation provides more temporal information to retrieve the unknown physical properties of the defect that is imaged. As such, the proposed research enhances non-destructively monitoring the integrity of structures, and practically it can be used to detect foreign obstacles concealed inside other objects.

Any tomography technique has at its heart a difficult inverse problem that needs to be solved, hence the mathematical analysis on the well- or ill-posedness of the model is necessary to be undertaken for scientific justification, as well as to be able to improve it. Difficulties arise due to the non-existence, non-uniqueness or the instability of solution, and establishing the degree of ill-posedness of the operator that needs inverting, is non-trivial. Moreover, concerning the actual inversion, to detect buried/hidden objects based on the transient approach one has to solve a difficult nonlinear and ill-posed moving boundary problem, which leads to a non-convex multi-dimensional optimization that needs to be further regularized to achieve stable results.

The thermographic principle of taking surface temperature measurements as we dynamically heat or cool an object offers an interesting transformative idea, but the approach is yet to be tested in an uncontrolled environment, and current understanding of the applicability of the technique to industrial scenarios is yet to be apprehended. Therefore, this mathematical modelling project using thermal-waves will provide a solid platform on which improved instrumentation for imaging can be built. The adventure of the research is to verify and validate the appropriateness of the new model by inverting both numerically simulated and experimental data in order to ultimately become available for real life application.
 
Description 1. Determination of a space-dependent source in the thermal-wave model of bio-heat transfer
This piece of research has considered the inverse problems of identifying the unknown space-dependent source in the thermal-wave model of bio-heat transfer from final-time or time-average temperature measurements. First, the uniqueness of the solutions of these inverse source problems has been proved using the separation-of-variables method.
Then, the inverse problems have been solved numerically by minimizing the least-squares objective functionals using the CGM combined with the discrepancy principle. To show the accuracy and stability of the numerical results, examples concerning the reconstruction of sources of different forms in one and two dimensions have been presented and thoroughly discussed.

2. Determination of the space-dependent blood perfusion coefficient in the thermal-wave model of bio-heat transfer
The reconstruction of the space-dependent perfusion coefficient and the temperature in the thermal-wave model of bio-heat transfer from final time temperature measurement has been investigated. For the numerical discretization, an unconditionally stable finite-difference method based on the Crank-Nicolson scheme has been used as a direct solver. This has been combined with a constrained regularized minimization problem. The resulting objective functional, penalized by a Tikhonov regularization term to restore the stability of the solution, has been minimized iteratively using the MATLAB optimization toolbox routine lsqnonlin. Accurate and stable numerical solutions for the unknown perfusion coefficient and the temperature, from both exact and noisy data, have been successfully achieved using the proposed computational method which has been verified for three benchmark numerical examples. Moreover, a dimensional blood perfusion rate of a biological tissue subjected to an external source of laser irradiation has been successfully identified.

3. Inverse problems of heterogeneous blood perfusion determination from non-intrusive boundary measurements
The inverse coefficient problem of recovering the unknown space-dependent blood perfusion coefficient in the hyperbolic thermal-wave model of bio-heat transfer from boundary temperature measurements has been investigated. Uniqueness and conditional Lipschitz stability have been established using the technique of Carleman estimates. These have been found valid over a time interval that is shorter than the usual previously reported in the literature. Moreover, the micro-local analysis performed leads to an improved stability estimate. Further, the problem has been reformulated as a nonlinear least-squares minimization problem and has been numerically solved using the conjugate gradient method (CGM) combined with the discrepancy principle for achieving stability. Numerical results associated with a physical example have been presented and discussed. Accurate and stable solutions have been obtained for both exact and noisy data using the proposed iterative CGM.
Exploitation Route There is a lot to follow from this investigation. In particular:
1. Further work associated with the thermal-wave model of bio-heat transfer is possible, e.g. determining the shape, size and location of tumours within biological tissues using non-intrusive boundary observations.
2. Extensions to more complex and potentially more accurate governing equations of bio-heat transfer, e.g. models of fractional-order or in integro-differential form.
3. Simultaneous recovery of the space-dependent source and perfusion coefficient in the thermal-wave model of bio-heat transfer can also be attempted
Sectors Digital/Communication/Information Technologies (including Software),Energy,Healthcare,Manufacturing, including Industrial Biotechology,Pharmaceuticals and Medical Biotechnology,Security and Diplomacy,Other

 
Description Immediate target markers/end users for this innovative interdisciplinary project, which will enhance the UK's position in tomographic reconstruction research, include nuclear, aerospace, security/defence and medical sectors. It is hoped that this emerging technology could lead to a new, safe, reliable, cost-effective and real-time procedure for monitoring and detecting concealed defects hidden in large containers or buried (at various depths and in various soils). It may also result in a more accurate prediction and assessment of the number, location, size and shape of anomalies (e.g. tumours, cracks), and a more accurate calculation of their evolvement in time.
First Year Of Impact 2023
Sector Aerospace, Defence and Marine,Healthcare
Impact Types Economic