Next Article in Journal
Design and Analysis of Broadband LiNbO3 Optical Waveguide Electric Field Sensor with Tapered Antenna
Previous Article in Journal
Novel Methodological Tools for Behavioral Interventions: The Case of HRV-Biofeedback. Sham Control and Quantitative Physiology-Based Assessment of Training Quality and Fidelity
Previous Article in Special Issue
Calculation of Stopping-Power Ratio from Multiple CT Numbers Using Photon-Counting CT System: Two- and Three-Parameter-Fitting Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Frequency Magnetic Induction Tomography System and Algorithm for Imaging Metallic Objects

by
Gavin Dingley
* and
Manuchehr Soleimani
*
Engineering Tomography Laboratory, Department of Electronic and Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY, UK
*
Authors to whom correspondence should be addressed.
Sensors 2021, 21(11), 3671; https://doi.org/10.3390/s21113671
Submission received: 3 February 2021 / Revised: 19 April 2021 / Accepted: 18 May 2021 / Published: 25 May 2021
(This article belongs to the Special Issue Tomography Sensing Technologies)

Abstract

:
Magnetic induction tomography (MIT) is largely focused on applications in biomedical and industrial process engineering. MIT has a great potential for imaging metallic samples; however, there are fewer developments directed toward the testing and monitoring of metal components. Eddy-current non-destructive testing is well established, showing that corrosion, fatigue and mechanical loading are detectable in metals. Applying the same principles to MIT would provide a useful imaging tool for determining the condition of metal components. A compact MIT instrument is described, including the design aspects and system performance characterisation, assessing dynamic range and signal quality. The image rendering ability is assessed using both external and internal object inclusions. A multi-frequency MIT system has similar capabilities as transient based pulsed eddy current instruments. The forward model for frequency swap multi-frequency is solved, using a computationally efficient numerical modelling with the edge-based finite elements method. The image reconstruction for spectral imaging is done by adaptation of a spectrally correlative base algorithm, providing whole spectrum data for the conductivity or permeability.

1. Introduction

Research in the development of magnetic induction tomography (MIT) has largely concentrated on biomedical and industrial process engineering [1,2,3,4]. Studies directed toward biomedical applications often use saline low-conductivity samples to experimentally simulate tissue; these same methods are also applied to imaging pipeline flows [5,6]. Imaging of molten metal flow, and associated trapped gasses, are also investigated [7,8]. The industrial application of MIT in imaging materials with high-conductivity and high-permeability has not received as much attention [9]. Non-destructive testing and eddy-current testing (NDT and ECT, respectively) of metal components has been extensively investigated and successfully adopted as a tool in industry. However, in using ECT to render images for detecting corrosion [10,11,12,13,14], fatigue [15,16], defects [17,18,19] and physical loading [20,21,22], mechanical scanners are often employed [23,24,25,26,27,28]. Furthermore, while ECT is generally regarded as non-contact, close proximity sensor positioning is usually required. In this paper, a recently developed system is described for primarily investigating the application of MIT in metal component imaging. The described unit is designed for simplicity and compactness so that it is easily relocatable to industrial environments. Li et al. have used a similar system for studying defect measurements in hot steel [29]. In this first development, the system operates from 1 to 100 kHz and collects the magnitude of the perturbed excitation field only. A description of the system is given, along with the approach used to characterise a pre-existing sensor array to assess performance. Finally, initial results are presented, demonstrating the operation of the new MIT system. Studies have shown that it is possible, using alternating magnetic fields, to determine the condition of metal components non-destructively without physical contact. Corrosion, mechanical loading and surface discontinuities, such as cracks, are detectable using alternating fields from 1 to 100 kHz. However, ECT is generally the dominant method in this particular application, while MIT-based approaches have not received as much attention. Depending on the electromagnetic properties of the metal under test, the nature of the field perturbation differs. In general, magnetic materials, those with high permeability, enhance the excitation field by focusing the lines of flux, therefore leading to increased signal strength. Materials with high conductivity tend to attenuate the applied field due to the skin effect, resulting in reduced signal strength [30]. In metals, such as many steels, where the object of interest has both high-conductivity and high-permeability, the applied field tends to travel over the surface of the sample, a combination of both the skin effect and magnetic flux focusing [31,32]. Measurement of corrosion in metals is particularly directed toward assessment of steel rebars, which are encased in the concrete of free-standing structures, therefore making non-contact, non-destructive detection of particular interest. Lower frequencies, in the order of 10 to 100 Hz, are often preferred for their greater penetration, where absolute or differential eddy-current probes are often used [10,11,23]. However, there are also studies using higher frequency excitation fields, in the range of 1 to 100 kHz, but these are mostly based on the transient behaviour of coil probes [12,13,14,17,24,33,34]. Although penetration is small at higher frequencies, corrosion is generally a surface effect and therefore detectable [24]. Muttakin and Soleimani demonstrate metallic material characterisation using spatio-spectral analysis over a frequency range of 100 Hz to 100 kHz [35]. Complex impedance data from MIT array sensors are used to interrogate a variety of metallic structures, indicating a change in both phase and magnitude, in response to both sensor position and excitation field frequency. In both testing and monitoring of metal components, it is possible to indirectly measure mechanical stresses in high permeability materials. Ferrous metals exhibit a change in magnetic susceptibility when mechanically loaded (inverse magnetostrictive or Villari effect). While under physical deformation, the alignment of magnetic domains within the material grain structure tends to align with the axis of the applied force [20]. Ferrous metals with helical domains, caused by torsional stress, also show a change in magnetic susceptibility; the Matteucci effect [21]. Variations in stress are therefore reflected as changes in the inductance of coils magnetically coupled to the sample, or the coupling factor of transformer style sensors, as found in MIT systems [36,37,38]. Residual stress in ferrous materials, resulting from quenching and phase transformations, are also indirectly detectable using eddy-current methods [15,16,20,22]. Li et al. describe a method for evaluating the yield stress in cold rolled steel, where the sample is magnetically biased with a low-frequency magnetic field, driving the sample through its hysteresis loop, while recording the complex impedance of a 5 kHz ECT sensor coil [25]. Eddy-current crack detection in metal components, again using high-frequency excitation fields, has also been demonstrated. The abrupt discontinuity in surface conductivity, associated with cracks, interrupts the flow of induced currents [18,19,26,27,28,39,40,41]. Crack depth is reflected in the detected signal magnitude, but is also dependent upon sensor positioning relative to the sample (lift-off effect). Studies have shown that stepped excitation frequency can help indicate crack depth, potentially providing a way of mitigating the effects of sensor location [42,43,44].

2. Principles of MIT Fundamentals

The following is a description of the general principles used in the design and development of the MIT system, including the method of tomographic image reconstruction. Using magnitude only data, the instrument gives a visual indication of the metal component structure using frequency domain data. However, the coils used to form the sensor array have a frequency-dependent sensitivity. Furthermore, the mutual position and orientation of individual coil pairs also affects relative sensitivity; these two contributing factors affect the sample space dynamic range. To mitigate these effects, the system gain is adapted for the field frequency in use and the particular coil pair in operation during frame acquisition.

2.1. System Design

Assessing the electronics for the MIT system requires a characterisation of the scanning array, determining both the electrical and physical parameters, from which the required sensitivity and dynamic range is established. Current flowing in the excitation coil induces a voltage in the detector coil, where the sample object interacts with the intermediating field; therefore, the electromagnetic properties of the object are reflected in the induced voltage. The time harmonic expression for the voltage v det induced in the detector coil of inductance L det , by a current i ex flowing in the excitation coil, of inductance L ex , in the case of an object-free sample space, is given by:
v det = j ω M i ex
where M is the mutual inductance, which is a function of the magnetic coupling factor k between L det and L ex , dependent upon the distance and relative angle between the two coils, given by M = k L det L ex . However, the excitation coils and detector coils alternate their function during scanning and therefore are of the same inductance ( L det = L ex = L ), such that M = k L . Given that only magnitude information is used in the present system, Equation (1) becomes:
| v det | = ω k L i ex
From Equation (2), the maximum possible signal is determined by setting ω to the maximum system signal frequency ω ( max ) and k = 1, as 0 k 1 by definition. The maximum signal magnitude is usually dictated by the available hardware, for instance the maximum input voltage to the ADC before saturation ( V FS ), combined with initial amplification ( G int ) and tailoring of the excitation coil current i ex , therefore, from Equation (2):
i ex = V FS G int ω ( max ) L
The voltage induced in the detector coil is dependent upon the excitation field frequency ( ω ) and coupling factor (k). So as to maintain a constant output signal magnitude, therefore maintaining the total numerical range, the system gain is configured for different coil pairs and excitation frequencies. The necessary additional system gain ( G sys ) is therefore given by:
G sys ( ω , k ) = ω ( max ) ω k
Setting the system gain to G sys , for a given coil pair coupling factor and frequency, while setting i e x to correspond with the ADC full-scale value (Equation (3)), ensures that the full dynamic range of the ADC is dedicated to capturing data related to the sample space.

2.2. Forward and Inverse Problems

In order to create MIT spectroscopic images, we need to simulate the measurement process in so-called forward modelling. The MIT forward model in each frequency can be described by means of eddy current models, and one can use the magnetic vector potential based formula, Equation (5).
× 1 μ × A + j ω σ A = J s
where A is the magnetic vector potential, J s is the excitation current density, ω is the angular frequency, while σ and μ are the conductivity and permeability of the sample, respectively. The excitation current for each frequency is modelled in transmitting coils as the right-hand side of Equation (5). The measured induced voltage in receiving coils can be writtin as:
V R = j ω V A · J 0 d V
This gives a complex-value-induced voltage, however, as our system is based on amplitude, we use the amplitude of the voltage as the value returned from the forward model. The image reconstruction in MIT spectroscopy is a non-linear inverse problem, but to be able to solve this in a timely manner, a linearisation is used. This is done by calculation of the Jacobian or sensitivity matrix using the results of the forward problem. A sensitivity matrix showing the changes of induced voltage with the electrical conductivity within a given area in the imaging domain is described by, for an excitation coil m and detector coil n:
V m n σ = j ω I m I n Ω E m · E n d Ω
While sensitivity to changes in permeability values can be described by:
V m n μ = j ω I m I n Ω H m · H n d Ω
We linearise the measured changes Δ V as a function of changes in conductivity or permeability, Δ x ; thus, the linear forward model can be written as:
Δ V = J Δ x
where J is the Jacobian, which is essentially a set of sensitivity distributions within the domain, which in this case is not square. The inverse MIT problem is also highly ill-posed, meaning it is sensitive to noise in measured data; in this case, regularisation methods are needed for a stable inverse solution. For the single frequency case, the Tikhonov Regularisation method is used to solve for Δ x , as shown by Equations (10) and (11) [45]:
Δ x = S Δ V
where the sensitivity matrix, S , is given by the Tikhonov Regularisation:
S = J T J + α 2 I 1 J T
where α is the regularisation parameter, usually determined as frequency-dependent, and I is the identity matrix. With a regularisation prior, which accounts for both spatial and spectral correlations between image elements. In spectral MIT, we have measured data in multiple frequencies, and the imaging results in neighbouring frequencies are similar. This allows for additional regulation in the frequency domain, for we can use a spectrally correlated algorithm, providing higher resolution and stability. Specifically, the data frame sequence is concatenated as Δ V f ˜ = Δ V f d ; ; Δ V f ; ; Δ V f + d and the resulting concatenated images Δ σ f ˜ = Δ X f d ; ; Δ X f ; ; Δ X f + d . Therefore, the temporal inverse problem can be rewritten as [46]:
Δ X f ˜ = S ^ Δ V f ˜
For the multi-frequency case, calculating the sensitivity matrix S , using Equation (11) is computationally intensive; therefore, an alternative approach is used, where a sensitivity matrix, S ^ , representing the complete set of frequencies is used, providing a least squares fit approximation. In the following, Δ x is a vector of 2500 elements, reshaped to a 50 by 50 image; however, an interpolation algorithm is used to smooth images. The quality of the rendered image depends on the sensitivity in the area of change, accuracy of measurement (SNR) and regularisation process. In the following sections, we attempt to link these relations for internal and external object imaging. For convenience, the frame data, Δ V , is normalised against the excitation frequency ( Δ V / f ), thus compensating for the array coil transfer function, as described by Equation (2).

3. MIT Hardware Design

In the following, details of the sensor array and system hardware operation are given. By using the analysis discussed in the previous section, the dynamic range of the instrument, as a function of frequency and array characteristics, is assessed.

3.1. Description of Array

The array is formed from eight coils positioned to form a 105 mm diameter circular sample space, with a coil angular displacement of 45 , as shown in Figure 1. Each coil central axis is positioned vertically 50 mm from the base of the sample space, the combined height of the base mount material and connecting cables.
Each coil has a diameter of 40 mm, wound with 55 turns of 32 swg (0.274 mm diam.) enamelled copper wire, giving a coil length of 15 mm; the axis coil centre to opposite coil centre distance is 140 mm. The measurement of the coils with an LCR instrument gives an average inductance of 126.97 μ H ( ± 0.55 ) and an average equivalent series resistance of 2.61 Ω ( ± 0.49 ). Free space coupling between the coils, where the sample space is empty, was measured for each of the four angular coil pair displacements, as indicated in Figure 1, resulting in eight individual measurements for 45 , 90 , 135 , and four for 180 . Table 1 summarises the average angular coupling factors and standard deviation.
Measuring the transmission between coil pairs is another method of characterising the array. Using a vector network analyser (VNA), S 21 phase and magnitude measurements were taken over 928 frequency points from 10 to 100 kHz. Phase and magnitude data for each coil pair angular displacement are shown in Figure 2. At the lower frequencies, the magnitude follows the trend given by Equation (2), while at the higher frequencies, the impedance of the excitation coil becomes apparent. The coil pair phase data in Figure 2 shows distinct phase inflections from approximately 40 to 50 kHz, depending on the particular coil pair angular displacement, suggesting a resonance at a frequency dependent upon the coupling factor. Impedance measurements of a single coil, derived from VNA S 11 reflection parameters, show a resonance occurring at 1 MHz; however, over the 1 to 100 kHz frequency range of interest, a single coil acts as a pure inductance.
From the array data given in Table 1, the required system gain, as a function of the coil pair angular displacement and excitation field frequency, is calculated with Equation (4), the result of which is shown graphically in Figure 3.

3.2. Hardware Description

Testing of metal samples in an industrial environment requires a physically compact system that is easily deployed; therefore, the hardware was constructed as a single stand-alone unit. System configuration, control and data acquisition were performed by an AVR 8-bit microcontroller (MCU), which communicates with a host PC over a USB bridge, as indicated in the system diagram Figure 4.
During the scanning operation, one coil is selected as the excitation source, while the other coils are multiplexed in sequence as the detector coil. There are two ADG1408 multiplexers, one of which selects the excitation coil, while the other selects the detector; both are controlled by the MCU GPIO port. The common signal input to the excitation multiplexer is fed with the current from a transconductance amplifier, built around an LT1010 power op-amp, set to output a current of 100 mA(rms). An AD9833 direct digital synthesizer (DDS) generates the 1 to 100 kHz tone that is fed to the transconductance amplifier, via a buffer stage to condition and set the sinusoidal signal amplitude. Tone frequency is set by the MCU over an SPI connection to the DDS, while the current amplitude is adjusted manually. The common signal output of the receiver multiplexer is fed to a coil pre-amplifier, which has a manually adjustable pre-set gain of nominally 22 dB and provides the initial gain ( G int ). Signal integrity is maintained by employing separate power and signal grounds connected at a single point. All coils are connected to the power ground plane, therefore providing a clear current return when set as the transmitter. The coil pre-amplifier is built around a Burr-Brown INA103KU instrumentation amplifier; therefore, the differential input is used to re-reference the amplifier output to the signal ground. Op-amp based precision rectifiers, incorporating a low-pass filter network (LPF), are used to convert the AC magnitude of the received signal into a DC voltage. Given the operating frequency range, the LPF of a single detector stage is not sufficient to provide both low-ripple voltage and fast settling time, therefore two selectable stages are used. Given the 40 dB sensitivity range of the detector coil over the total frequency span, a 20 dB amplifier stage precedes the lower frequency absolute detector. The absolute detector DC voltage level outputs are digitised by an ADS1115 ADC, where two of the four multiplexed analogue inputs are used. The ADS1115 also contains a programmable amplifier, with gain settings of effectively 0, 2, 8, 14, 20 and 26 dB, configured by the MCU over an I2C bus. Selecting the appropriate analogue channel of the ADC, for the frequency range of interest, while setting the ADS1115 gain, the effects of coil frequency response and coil pair coupling are compensated, therefore dedicating the dynamic range of the ADC to capturing sample space data. A total system gain of 20 to 66 dB is therefore available; however, from Figure 3, it is apparent that the upper gain is deficient by 21 dB, and this reduction in sensitivity will impact the dynamic range dedicated to the sample space. The ADS1115 has a 16 bit ADC; however, this is for bipolar operation, so only 15 bits are used in the present application, resulting in a total dynamic range of approximately 90 dB (6.02 × N bits ). However, as a result of the deficiency in system gain, 21 dB of the ADC dynamic range is allocated to compensating for array detector coil sensitivity; therefore, 69 dB remains for the sample space.

4. Results and Discussion

At present, the system has a frame acquisition time of 750 ms, partly attributed to the 10 ms dwell time of the absolute detector, therefore contributing 280 ms to the total acquisition period, i.e., 10 ms × N ( N 1 ) / 2 , where N is the number of coils (in the present case, eight). Serial communications and delays resulting from PC activity contribute the remaining 470 ms. Absolute detector dwell times are directly related to the field excitation frequency; in the present system, this is fixed to 10 ms to cater for the lowest expected frequency of 1 kHz. A bank of absolute detectors tailored for sections of the total frequency range, in terms of amplification and dwell time would improve both dynamic range and frame acquisition time. To assess the performance of the system, measurements were taken to determine the quality of frame data and resulting reconstructed images.

4.1. SNR Measurements of Frame Data

Figure 5 shows the signal to noise ratio (SNR) as a function of field excitation frequency. A total of 100 sets of frame data were taken to determine the mean and standard deviation for each coil pair, from which the average SNR over a frame was taken for each frequency point. A total of 600 frequency points were taken from 1 to 100 kHz—the nominal operating range of the system.
Two artefacts in Figure 5 are apparent; first, a notch at 42 kHz with a depth of 6 dB, which corresponds with the phase inflection in Figure 2b, identified as a resonance. Although small, the resonance creates a local instability in the frequency sweep, resulting in slightly decreased SNR. A second larger artefact in Figure 5 occurs at 84 kHz, with a decrease in SNR of approximately 14 dB. On further investigation, it was found that at this frequency, the excitation field increases by 6 dB, inducing an instability in the transconductance amplifier, therefore decreasing SNR. The harmonic relation between the smaller artefact at 42 kHz and the larger decrease in SNR at 84 kHz suggests that the larger artefact results from the interaction of the transconductance amplifier with the coil pair resonance. Figure 6 shows the measured SNR over a single frame, i.e., SNR as a function of the coil pair, for 1, 10 and 100 kHz excitation fields.
From Figure 6, it is apparent that the SNR is dependent upon the coil pair, which directly corresponds with the coupling factor. However, the frequency of the excitation field has a greater influence, as also indicated by Figure 5.

4.2. Measurements of Reconstructed Images

System performance measurements based on reconstructed images from frame data were made using the algorithms described by Adler et al. [47]. Of particular interest are the relative resolution and relative image deformation with frequency, particularly when other conducting objects are introduced into the sample space with the sample of interest. Quantitative assessment of image reconstruction begins with determining the optimum amplitude threshold of the column vector x ^ representing the reconstructed image. As the objects are clearly defined and known, the threshold value is determined empirically and defined by a factor k t h , where 1 k t h 0 , resulting in the amplitude-limited vector x ^ q , from Equation (13):
[ x ^ q ] i = 1 if [ x ^ ] i k t h · max ( x ^ ) 0 otherwise
Image resolution (RES) is determined by the area ratio of the amplitude-limited image, represented by vector x ^ q and the total image area ( A 0 ), measured in pixels; therefore, it has a value 1 RES 0 . A square root is used to represent this ratio in terms of radius; see Equation (14):
RES = k [ x ^ q ] k A 0
Shape deformation (SD) is defined as the fraction of the threshold-limited, reconstructed image ( [ x ^ q ] ) outside a circular boundary of equal area; therefore, it is a ratio with a value 1 SD 0 , given by Equation (15):
SD = k C [ x ^ q ] k k [ x ^ q ] k
Measurements were taken of an aluminium rod sample, 100 mm in length and 10 mm in diameter. In the first set of measurements, the sample rod was positioned at a radius of 40 mm from the centre of the sample space. An identical rod was then positioned symmetrically opposite the sample, located at an equal radial distance from the sample space centre, see Figure 7. Data were collected for 100 frequency points, over a range of 1 to 100 kHz. Background data, where the sample space is empty, were taken by averaging 100 frames ( V bg ). With the rods positioned in the sample space, each frequency point was averaged over five frames of data ( V both ), from which the background data were subtracted ( Δ V = V both V bg ). Images were reconstructed using Tikhonov regularisation, as described previously in Equation (11), then the threshold condition was applied, as described by Equation (13). It is noted that the regularisation parameter α in Equation (11) should be frequency dependent; however, in this case, the same value was used across the frequency range. Figure 8 shows the reconstructed images for 5, 50 and 95 kHz. From Equations (14) and (15), the resolution and shape deformation of the reconstructed image are calculated for each frequency point, as shown in Figure 9.
Figure 10 shows the variation in image maximum ( max ( Δ x ) ) as a function of excitation field frequency for both the single rod in the sample space ( Δ V = V sgl V bg ) and with the addition of the identical external rod ( Δ V = V both V sgl ).
Figure 8 suggests that there is very little variation in image quality with field excitation frequency, although Figure 8a indicates some degradation at the lower frequencies, caused by lower SNR; see also Figure 5. Figure 9 gives a better indication of the image quality as a function of frequency, where below 20 kHz, both the resolution and deformation are relatively poor; however, they improve with increased field frequency, again following the SNR trend in Figure 5. There is a slight degradation in resolution as a result of the previously identified 84 kHz resonance, but this is a direct consequence of the associated decrease in SNR. Figure 10 indicates that there is very little difference in the image maximum with or without the inclusion of an identical rod sample in the sample space, which corresponds with the image quality observations drawn from Figure 8. Resonance effects are apparent, and a system frequency response below approximately 10 kHz is clearly shown, resulting in the lower SNR and, hence, image quality at these frequencies.
Reconstructed concealed object image measurements were also taken of a partially enclosed aluminium rod, centrally positioned within a C-shaped cylindrical shroud, 42 mm in diameter and with a height of 110 mm; Figure 11 shows the arrangement. The conducting material forming the shroud is 40 μ m thick copper, so minimal attenuation was provided to the applied excitation field. Both the rod and shroud were placed in the centre of the sample space, where again 100 frequency points were recorded, each averaged over five frames. Reconstructed images with the rod and shroud present are shown in Figure 12, where the background data are subtracted from the frame data ( Δ V = V both V bg ). Figure 13 is reconstructed from the frame data with the rod and shroud in position, but with the frame data of the shroud subtracted ( Δ V = V both V shrd ), thus leaving the rod. Figure 14 shows plots of both the resolution (RES) and shape deformation (SD) as a function of frequency, derived from the images in Figure 13 using Equations (13)–(15), i.e., the data from both rod and shroud, with the data from the shroud subtracted ( Δ V = V both V shrd ).
Figure 15 shows the variation in image maximum ( max ( Δ x ) ) as a function of frequency for both the rod ( Δ V = V both V shrd ) and shroud ( Δ V = V shrd V bg ). With increasing excitation field frequency, the internally-included shroud has increased visibility as a result of the skin effect, therefore decreasing the visibility of the rod by its shielding action.
Figure 12 indicates that there is little variation in the reconstructed image with frequency; however, the shroud is dominant and hides the image of the rod. Subtracting the frame data of the shroud alone, without the rod, reveals the internal object, as shown in Figure 13. It is apparent from Figure 13 that at higher field frequencies the skin effect of the shroud shields the rod, thus attenuating the image. Figure 15 further indicates how the presence of the shroud in rendered images is more dominant at higher frequencies, while that of the sample rod is reduced. Below approximately 10 kHz, the frequency response is low, as a result of the limitations of the system; however, beyond this point, the rod has a strong presence in the rendered image, but rapidly decreases with frequency, as it is shielded by the shroud. Figure 14 shows the resolution and shape deformation of the reconstructed images in Figure 13 as a function of frequency, which indicates that the rod image resolution improves slightly with increased frequency, while the shape deformation degrades. Increased resolution results from the improved SNR with frequency, as indicated by Figure 5, while the shape deformation is dependent upon the increasing influence of the conducting shroud with frequency, due to the skin effect.

4.3. Discussion of Results

Resolution and shape deformation are dependent upon frame data SNR, where any decrease in signal or increase in noise results in diminished performance. Reduction in coil pair transmission, from either low excitation frequency, coil coupling, or the presence of shielding, decreases signal strength. Although a lower excitation frequency corresponds with a smaller signal induced in the detector coil, it also provides greater shield penetration, thus presenting a trade-off. Noise primarily originates from solid-state devices in detector amplifiers in the form of Gaussian noise. However, the excitation coil driver, in this case, the transconductance amplifier, can also introduce noise in the form of instabilities. Self-resonance of array coils can cause instability noise when either acts as a detector or field generator; in the present case, the latter is most apparent. Although characterisation of an array provides insight into the dynamic-range dedicated to the sample space, frequency domain measurements of coil impedance and coil pair transmission are also important, indicating the presence of excitation frequencies offering reduced performance.

5. Conclusions

Although the application of MIT for metal component testing is growing, the methods described for eddy-current testing suggest the potential for generating cross-sectional images indicating corrosion and fatigue. A compact magnitude-only MIT system was developed for investigating metal component imaging in industrial test facilities. Given the levels of excitation-field perturbation over a range of metals, having either high conductivity, high permeability or a combination of both, a significant dynamic range is needed, which was analytically quantified for a given sensor array. The ability to alter the excitation field frequency opens the possibility of spectral based image reconstruction, as well as greater image rendering stability, for which an algorithm was presented. Initial test results were given indicating the performance of the new system, particularly in discriminating metallic samples within the sample space. While image quality improves with frequency for external objects, there is a degradation in image quality for concealed objects. Characterisation of a chosen array not only gives an indication of required system gain, but also identifies the presence of parasitic resonances that can potentially degrade rendered image quality.

Author Contributions

Conceptualization, G.D. and M.S.; methodology, G.D. and M.S.; software, G.D. and M.S.; validation, M.S.; formal analysis, G.D. and M.S.; investigation, G.D.; resources, G.D. and M.S.; data curation, G.D.; writing–original draft preparation, G.D. and M.S.; writing–review and editing, G.D. and M.S.; visualization, G.D. and M.S.; supervision, M.S.; project administration, M.S.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ma, L.L.; Soleimani, M. Magnetic induction tomography methods and applications: A review. Meas. Sci. Technol. 2017, 28, 072001. [Google Scholar] [CrossRef] [Green Version]
  2. Wei, H.Y.; Soleimani, M. Electromagnetic Tomography for Medical and Industrial Applications: Challenges and Opportunities [Point of View]. Proc. IEEE 2013, 101, 559–565. [Google Scholar] [CrossRef]
  3. Griffiths, H. Magnetic induction tomography. Meas. Sci. Technol. 2001, 12, 1126–1131. [Google Scholar] [CrossRef]
  4. Igney, C.H.; Pinter, R.; Such, O. Magnetic Induction Tomography System and Method. U.S. 20080258717A1, 23 October 2008. [Google Scholar]
  5. Wei, H.Y.; Soleimani, M. Hardware and software design for a National Instrument-based magnetic induction tomography system for prospective biomedical applications. Physiol. Meas. 2012, 33, 863–879. [Google Scholar] [CrossRef]
  6. Watson, S.; Williams, R.J.; Gough, W.; Griffiths, H. A magnetic induction tomography system for samples with conductivities below 10 S/m. Meas. Sci. Technol. 2008, 19, 045501. [Google Scholar] [CrossRef]
  7. Terzija, N.; Yin, W.; Gerbeth, G.; Stefani, F.; Timmel, K.; Wondrak, T.; Peyton, A.J. Use of electromagnetic induction tomography for monitoring liquid metal/gas flow regimes on a model of an industrial steel caster. Meas. Sci. Technol. 2011, 22. [Google Scholar] [CrossRef]
  8. Ma, X.; Peyton, A.J.; Higson, S.R.; Drake, P. Development of multiple frequency electromagnetic induction systems for steel flow visualization. Meas. Sci. Technol. 2008, 19, 094008. [Google Scholar] [CrossRef]
  9. Wei, H.Y.H.Y.; Soleimani, M. A Magnetic Induction Tomography System for Prospective Industrial Processing Applications. Chin. J. Chem. Eng. 2012, 20, 406–410. [Google Scholar] [CrossRef]
  10. Bailey, J.; Long, N.; Hunze, A. Eddy Current Testing with Giant Magnetoresistance (GMR) Sensors and a Pipe-Encircling Excitation for Evaluation of Corrosion under Insulation. Sensors 2017, 17, 2229. [Google Scholar] [CrossRef] [Green Version]
  11. Guilizzoni, R.; Finch, G.; Harmon, S. Subsurface Corrosion Detection in Industrial Steel Structures. IEEE Magn. Lett. 2019, 10, 1–5. [Google Scholar] [CrossRef]
  12. De Alcantara, N.; da Silva, F.; Guimarães, M.; Pereira, M. Corrosion Assessment of Steel Bars Used in Reinforced Concrete Structures by Means of Eddy Current Testing. Sensors 2015, 16, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Minesawa, G.V.; Sasaki, E. Eddy current inspection of corrosion defects for concrete embedded steel members. AIP Conf. Proc. 2014, 781, 781–786. [Google Scholar] [CrossRef]
  14. Frankowski, P.P.K. Corrosion detection and measurement using eddy current method. In Proceedings of the 2018 International Interdisciplinary PhD Workshop (IIPhDW), Swinoujscie, Poland, 9–12 May 2018; pp. 398–400. [Google Scholar] [CrossRef]
  15. Oka, M.; Yakushiji, T.; Tsuchida, Y.; Enokizono, M. Evaluation of fatigue damage in an austenitic stainless steel (SUS304) using the eddy current probe. In Proceedings of the INTERMAG Asia 2005, Digests of the IEEE International Magnetics Conference, Nagoya, Japan, 4–8 April 2005; pp. 427–428. [Google Scholar] [CrossRef]
  16. Oka, M.; Tsuchida, Y.; Yakushiji, T.; Enokizono, M. Fatigue Evaluation for a Ferritic Stainless Steel (SUS430) by the Eddy Current Method Using the Pancake-Type Coil. IEEE Trans. Magn. 2010, 46, 540–543. [Google Scholar] [CrossRef]
  17. Zheng, D.; Dong, Y. A Dual-parameter Oscillation Method for Eddy Current Testing with The Aid of Impedance Nonlinearity. IEEE Trans. Instrum. Meas. 2019, 69, 4476–4486. [Google Scholar] [CrossRef]
  18. Tsukada, K.; Tomioka, T.; Wakabayashi, S.; Sakai, K.; Kiwa, T. Magnetic detection of steel corrosion at a buried position near the ground level using a magnetic resistance sensor. IEEE Trans. Magn. 2018, 54, 1. [Google Scholar] [CrossRef]
  19. Gao, P.; Wang, X.; Han, D.; Zhang, Q. Eddy current testing for weld defects with different directions of excitation field of rectangular coil. In Proceedings of the 2018 4th International Conference on Control, Automation and Robotics (ICCAR), Auckland, New Zealand, 20–23 April 2018; pp. 486–491. [Google Scholar] [CrossRef]
  20. Uchanin, V.; Minakov, S.; Nardoni, G.; Ostash, O.; Fomichov, S. Nondestructive Determination of Stresses in Steel Components by Eddy Current Method. Strojniški Vestnik-J. Mech. Eng. 2018, 64, 690–697. [Google Scholar] [CrossRef]
  21. Charubin, T.; Nowicki, M.; Szewczyk, R. Spectral analysis of Matteucci effect based magnetic field sensor. AIP Conf. Proc. 1996, 2018, 020019. [Google Scholar] [CrossRef]
  22. Chady, T.; Sikora, R.; Psuj, G.; Enokizono, M.; Todaka, T. Fusion of electromagnetic inspection methods for evaluation of stress-loaded steel samples. IEEE Trans. Magn. 2005, 41, 3721–3723. [Google Scholar] [CrossRef]
  23. Lo, C.C.H.; Nakagawa, N. Evaluation of eddy current and magnetic techniques for inspecting rebars in bridge barrier rails. AIP Conf. Proc. 2013, 1511, 1371–1377. [Google Scholar] [CrossRef] [Green Version]
  24. Miller, G.; Gaydecki, P.; Quek, S.; Fernandes, B.T.; Zaid, M.A. Detection and imaging of surface corrosion on steel reinforcing bars using a phase-sensitive inductive sensor intended for use with concrete. NDT E Int. 2003, 36, 19–26. [Google Scholar] [CrossRef]
  25. Li, K.; Li, L.; Wang, P.; Liu, J.; Shi, Y.; Zhen, Y.; Dong, S. A fast and non-destructive method to evaluate yield strength of cold-rolled steel via incremental permeability. J. Magn. Magn. Mater. 2020, 498. [Google Scholar] [CrossRef]
  26. Postolache, O.; Pereira, M.D.; Ramos, H.G.; Ribeiro, A.L. NDT on Aluminum Aircraft Plates based on Eddy Current Sensing and Image Processing. In Proceedings of the 2008 IEEE Instrumentation and Measurement Technology Conference, Victoria, BC, Canada, 12–15 May 2008; pp. 1803–1808. [Google Scholar] [CrossRef]
  27. Bo, L.; Feilu, L.; Zhongqing, J.; Jiali, L. Eddy Current Array Instrument and Probe for Crack Detection of Aircraft Tubes. In Proceedings of the 2010 International Conference on Intelligent Computation Technology and Automation, Changsha, China, 11–12 May 2010; Volume 2, pp. 177–180. [Google Scholar] [CrossRef]
  28. Arismendi, N.; Pacheco, E.; Lopez, O.; Espina-Hernandez, J.; Benitez, J. Classification of artificial near-side cracks in aluminium plates using a GMR-based eddy current probe. In Proceedings of the 2018 28th International Conference on Electronics, Communications and Computers, CONIELECOMP 2018, Cholula, Mexico, 21–23 February 2018; pp. 31–36. [Google Scholar] [CrossRef]
  29. Li, F.; Spagnul, S.; Odedo, V.; Soleimani, M. Monitoring Surface Defects Deformations and Displacements in Hot Steel Using Magnetic Induction Tomography. Sensors 2019, 19, 3005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Peyton, A.J.; Beck, M.S.; Borges, A.R.; De Oliveira, J.E.; Lyon, G.M.; Yu, Z.Z.; Brown, M.W.; Ferrerra, J. Development of Electromagnetic Tomography (EMT) for Industrial Applications. Part 1: Sensor Design and Instrumentation. In Proceedings of the 1st World Congress on Industrial Process Tomography, Buxton, Greater Manchester, UK, 14–17 April 1999; pp. 306–312. [Google Scholar]
  31. Rubinacci, G.; Tamburrino, A.; Ventre, S. Concrete rebars inspection by eddy current testing. Int. J. Appl. Electromagn. Mech. 2007, 25, 333–339. [Google Scholar] [CrossRef]
  32. Borges, A.R.; de Oliveira, E.; Velez, J.; Peyton, A.J.; Tavares, C. Development of electromagnetic tomography (EMT) for industrial application: Part 2: Image reconstruction and software framework. World Congr. Ind. Process. Tomogr. 1999, 2, 219–225. [Google Scholar]
  33. Miller, G.; Gaydecki, P.; Quek, S.; Fernandes, B.; Zaid, M. A combined Q and heterodyne sensor incorporating real-time DSP for reinforcement imaging, corrosion detection and material characterisation. Sens. Actuators Phys. 2005, 121, 339–346. [Google Scholar] [CrossRef]
  34. Zheng, D.; Dong, Y. A Novel Eddy Current Testing Scheme by Transient Oscillation and Nonlinear Impedance Evaluation. IEEE Sens. J. 2018, 18, 4911–4919. [Google Scholar] [CrossRef]
  35. Muttakin, I.; Soleimani, M. Magnetic Induction Tomography Spectroscopy for Structural and Functional Characterization in Metallic Materials. Materials 2020, 13, 2639. [Google Scholar] [CrossRef] [PubMed]
  36. Schonekess, H.; Ricken, W.; Becker, W.J. Improved multi-sensor for force measurement on pre-stressed steel cables by means of eddy current technique. In Proceedings of the IEEE Sensors, 2004, Vienna, Austria, 24–27 October 2004; pp. 260–263. [Google Scholar] [CrossRef]
  37. Cao, B.; Iwamoto, T.; Bhattacharjee, P.P. An experimental study on strain-induced martensitic transformation behavior in SUS304 austenitic stainless steel during higher strain rate deformation by continuous evaluation of relative magnetic permeability. Mater. Sci. Eng. A 2020, 774, 138927. [Google Scholar] [CrossRef]
  38. Xu, H.; Lu, M.; Avila, J.R.S.; Zhao, Q.; Zhou, F.; Meng, X.; Yin, W. Imaging a weld cross-section using a novel frequency feature in multi-frequency eddy current testing. Insight-Non-Destr. Test. Cond. Monit. 2019, 61, 738–743. [Google Scholar] [CrossRef]
  39. Obeid, S.; Dogaru, T.; Tranjan, F.M. Rotational GMR magnetic sensor based eddy current probes for detecting buried corner cracks at the edge of holes in metallic structures. In Proceedings of the IEEE SoutheastCon 2008, Huntsville, AL, USA, 3–6 April 2008; pp. 314–317. [Google Scholar] [CrossRef]
  40. Dalal Radia, T.; Ahmed, D.; Bachir, H.; Khaldoun, L.I. Detection of Defects Using GMR and Inductive Probes. In Smart Energy Empowerment in Smart and Resilient Cities; Springer: Berlin/Heidelberg, Germany, 2020; Volume 102, pp. 617–622. [Google Scholar] [CrossRef]
  41. Sen, T.; Anoop, C.S.; Sen, S. Study and analysis of two GMR-based eddy-current probes for defect-detection. In Proceedings of the 2017 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Turin, Italy, 22–25 May 2017; pp. 1–6. [Google Scholar] [CrossRef]
  42. Lopes Ribeiro, A.; Pasadas, D.; Ramos, H.; Rocha, T. Determination of crack depth in aluminum using eddy currents and GMR sensors. AIP Conf. Proc. 2015, 1650, 361–367. [Google Scholar] [CrossRef]
  43. Bernieri, A.; Betta, G.; Ferrigno, L.; Laracca, M. Multi-frequency ECT method for defect depth estimation. In Proceedings of the 2012 IEEE Sensors Applications Symposium Proceedings, Brescia, Italy, 7–9 February 2012; pp. 1–6. [Google Scholar] [CrossRef]
  44. Bernieri, A.; Betta, G.; Ferrigno, L.; Laracca, M. Crack Depth Estimation by Using a Multi-Frequency ECT Method. IEEE Trans. Instrum. Meas. 2013, 62, 544–552. [Google Scholar] [CrossRef]
  45. Neubauer, A. Tikhonov regularisation for non-linear ill-posed problems: Optimal convergence rates and finite-dimensional approximation. Inverse Probl. 1989, 5, 541–557. [Google Scholar] [CrossRef]
  46. Soleimani, M.; Mitchell, C.N.; Banasiak, R.; Wajman, R.; Adler, A. Four-Dimensional Electrical Capacitance Tomography Imaging Using Experimental Data. Prog. Electromagn. Res. 2009, 90, 171–186. [Google Scholar] [CrossRef] [Green Version]
  47. Adler, A.; Arnold, J.H.; Bayford, R.; Borsic, A.; Brown, B.; Dixon, P.; Faes, T.J.C.; Frerichs, I.; Gagnon, H.; Gärber, Y.; et al. GREIT: A unified approach to 2D linear EIT reconstruction of lung images. Physiol. Meas. 2009, 30, S35–S55. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Diagram of array, indicating coil position and dimensions.
Figure 1. Diagram of array, indicating coil position and dimensions.
Sensors 21 03671 g001
Figure 2. Coil pair transmission measurements ( S 21 ), for the four different angular displacements.
Figure 2. Coil pair transmission measurements ( S 21 ), for the four different angular displacements.
Sensors 21 03671 g002
Figure 3. Graph showing required system gain as a function of frequency for different coil pair angles.
Figure 3. Graph showing required system gain as a function of frequency for different coil pair angles.
Sensors 21 03671 g003
Figure 4. System diagram of the MIT instrument.
Figure 4. System diagram of the MIT instrument.
Sensors 21 03671 g004
Figure 5. SNR as a function of excitation field frequency.
Figure 5. SNR as a function of excitation field frequency.
Sensors 21 03671 g005
Figure 6. Frame data SNR for excitation field frequencies 1, 10 and 100 kHz.
Figure 6. Frame data SNR for excitation field frequencies 1, 10 and 100 kHz.
Sensors 21 03671 g006
Figure 7. Photograph of inclusion measurements using aluminium rod samples 10 mm in diameter and 100 mm in length, with a 105 mm diameter circular sample space.
Figure 7. Photograph of inclusion measurements using aluminium rod samples 10 mm in diameter and 100 mm in length, with a 105 mm diameter circular sample space.
Sensors 21 03671 g007
Figure 8. Reconstructed images from aluminium external object inclusion measurements ( Δ V = V both V bg ), acquired using 5, 50 and 95 kHz excitation field frequencies. Images represent a circular sample space of 105 mm diameter.
Figure 8. Reconstructed images from aluminium external object inclusion measurements ( Δ V = V both V bg ), acquired using 5, 50 and 95 kHz excitation field frequencies. Images represent a circular sample space of 105 mm diameter.
Sensors 21 03671 g008
Figure 9. Reconstructed image resolution (RES) and shape deformation (SD) of dual aluminium external object inclusion data, applying a threshold value of k t h = 0.3 in Equations (13)–(15), as a function of excitation field frequency 2 to 100 kHz.
Figure 9. Reconstructed image resolution (RES) and shape deformation (SD) of dual aluminium external object inclusion data, applying a threshold value of k t h = 0.3 in Equations (13)–(15), as a function of excitation field frequency 2 to 100 kHz.
Sensors 21 03671 g009
Figure 10. Image maximum ( max ( Δ x ) ) as a function of frequency, for the rendered images of the aluminium rod sample individually ( Δ V = V sgl V bg ) and with the including object ( Δ V = V both V sgl ).
Figure 10. Image maximum ( max ( Δ x ) ) as a function of frequency, for the rendered images of the aluminium rod sample individually ( Δ V = V sgl V bg ) and with the including object ( Δ V = V both V sgl ).
Sensors 21 03671 g010
Figure 11. Photographs of internal object measurements using an aluminium rod and copper shroud.
Figure 11. Photographs of internal object measurements using an aluminium rod and copper shroud.
Sensors 21 03671 g011
Figure 12. Reconstructed image measurements of aluminium rod partially concealed in a copper shroud, acquired using 5, 50 and 95 kHz excitation field frequencies. Image constructed from frame data of both rod and shroud present in the sample space ( Δ V = V both V bg ).
Figure 12. Reconstructed image measurements of aluminium rod partially concealed in a copper shroud, acquired using 5, 50 and 95 kHz excitation field frequencies. Image constructed from frame data of both rod and shroud present in the sample space ( Δ V = V both V bg ).
Sensors 21 03671 g012
Figure 13. Reconstructed image of aluminium rod partially encapsulated in a copper shroud, with shroud frame data subtracted ( Δ V = V both V shrd ), leaving the rod; acquired using 5, 50 and 95 kHz excitation field frequencies. Image (a) clearly shows the aluminium rod, while the influence of the shroud at higher frequencies is apparent in (b,c), where the rod is weaker and distorted.
Figure 13. Reconstructed image of aluminium rod partially encapsulated in a copper shroud, with shroud frame data subtracted ( Δ V = V both V shrd ), leaving the rod; acquired using 5, 50 and 95 kHz excitation field frequencies. Image (a) clearly shows the aluminium rod, while the influence of the shroud at higher frequencies is apparent in (b,c), where the rod is weaker and distorted.
Sensors 21 03671 g013
Figure 14. Reconstructed image resolution (RES) and shape deformation (SD) of internal object measurement using a conducting shroud, where shroud frame data are subtracted ( Δ V = V both V shrd ). Using a threshold value of k t h = 0.3 in Equations (13)–(15), frame data were collected over a frequency range of 2 to 100 kHz.
Figure 14. Reconstructed image resolution (RES) and shape deformation (SD) of internal object measurement using a conducting shroud, where shroud frame data are subtracted ( Δ V = V both V shrd ). Using a threshold value of k t h = 0.3 in Equations (13)–(15), frame data were collected over a frequency range of 2 to 100 kHz.
Sensors 21 03671 g014
Figure 15. Image maximum as a function of frequency, for the rendered images for both the aluminium rod sample and the internally-included copper shroud independently.
Figure 15. Image maximum as a function of frequency, for the rendered images for both the aluminium rod sample and the internally-included copper shroud independently.
Sensors 21 03671 g015
Table 1. Average coupling factors of coil pairs with the four angular displacements.
Table 1. Average coupling factors of coil pairs with the four angular displacements.
AngleAverage kSTD k
45°0.0440.0037
90°0.00750.0029
135°0.00460.0023
180°0.00460.00218
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dingley, G.; Soleimani, M. Multi-Frequency Magnetic Induction Tomography System and Algorithm for Imaging Metallic Objects. Sensors 2021, 21, 3671. https://doi.org/10.3390/s21113671

AMA Style

Dingley G, Soleimani M. Multi-Frequency Magnetic Induction Tomography System and Algorithm for Imaging Metallic Objects. Sensors. 2021; 21(11):3671. https://doi.org/10.3390/s21113671

Chicago/Turabian Style

Dingley, Gavin, and Manuchehr Soleimani. 2021. "Multi-Frequency Magnetic Induction Tomography System and Algorithm for Imaging Metallic Objects" Sensors 21, no. 11: 3671. https://doi.org/10.3390/s21113671

APA Style

Dingley, G., & Soleimani, M. (2021). Multi-Frequency Magnetic Induction Tomography System and Algorithm for Imaging Metallic Objects. Sensors, 21(11), 3671. https://doi.org/10.3390/s21113671

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop