1. Introduction
The current flows in power electronic devices not only reflect the electromagnetic compatibility, joule heating sources, and impedance performances of devices under test (DUT), but also determine the fault hazard, reliability, as well as stability of the electronic system.
State-of-the-art electronic device detection relies on lumped parameter measurements of the current–voltage (I–V) curve, impedance spectrum, contact resistance and S parameter, etc., by connecting the test probe to the circuit net or node. These lumped parameters are believed to be useful in identifying the running status of DUT, but are uncapable of indicating the abnormal reason and fault position.
Magnetic current imaging (MCI) is a powerful measurement tool to visualize the current flows in DUT by non-invasively sensing the magnetic signature [
1]. This novel technique provides a spatial distribution of the current density (CD) in an image manner to quantitively and straightly localize abnormal positions, as illustrated in
Figure 1.
Magnetic current imaging was firstly proposed for bio-magnetic detection, and applied in an integrated circuit test by Bradley J. Roth and E. F. Fleet, respectively [
2,
3]. Existing MCIs primarily use Fourier transform back-evolution for the vertical component of the magnetic field above the DUT to reconstruct its two-dimensional CD magnitude distribution [
4,
5]. The superiorities of in situ, visibility, and non-invasive conditions promoted MCI applications in the fields of integrated circuit fault diagnosis, superconductor magnet analysis, and lithium-ion battery health prognosis [
1,
6,
7,
8].
Benaiah D. Schrag et al. used magnetic tunnel junction sensors to image the CD magnitude in a self-made application-specific integrated circuit. Current and spatial resolutions of 1 mA and 1 mm, respectively, were achieved over probe-to-sample separation of a 4 mm scanning area of 2 × 2 mm
2. This work discussed the imaging quality dependence on measurement noise and separation distance, however, lacked a technical solution [
9]. J. Gaudestad et al. implemented a
μm-level resolution of CD and successfully localized the die level short-circuit used in giant magneto resistance sensors [
6]. Pauli Kehayias et al. imaged the surface current density magnitude in 555 timer IC by measuring magnetic fields using a quantum diamond microscope, achieving the highest micron-scale spatial resolution and few-micro-Tesla magnetic sensitivity in a 1 × 1 μm
2 pixel at a probe-to-sample separation of 26 μm [
10]. M. Sumi and N. Satoh obtained the CD contribution on a Cockcroft–Walton circuit, a Schottky rectifier and a multi-layered ceramic capacitor. An electric short circuit spot was successfully observed by magneto-impedance sensor within an operation area of 50 × 50 mm
2 [
11]. Mark G. Bason et al. achieved internal CD magnitude imaging in lithium-ion batteries using fluxgate sensor arrays within a scanning area of 150 cm
2 to give a spatial resolution of 5 mm [
7]. Felix Brauchle et al. demonstrated the discovery of cuts and erroneous weld spots inside lithium-ion cells using anisotropic magnetic resistance sensors with a current resolution of 0.15 mA, with a probe-to-sample separation of 7.8 mm. On the other hand, the accuracy of CD in this work is 227 mA/cm
2 at a local resolution of 4 mm
2 [
12]. Their CD imaging results intuitively distinguished CD pattern differences between a normal and faulty status, however, were unable to offer CD directions [
13]. E. Marchiori et al. showed the CD in the superconducting qubit circuit using a superconducting quantum interference device (SQUID) to give a resolution of 3 mA and a probe-to-sample separation of 600 nm [
8].
The emergence and development of MCI enable a new insight to reveal the essential electrical property of the DUT. However, the latest MCI technique is limited in providing CD magnitude pattern, and uncapable of offering vectorial current density (VCD) information, which is of great significance in deriving divergence and rotation signatures to identify current leakage and an equivalent magnetic induction. Moreover, the magnetic field is the fundamental information for a CD reconstruction algorithm; therefore, the magnetic noise will finally couple to the inversed CD results and determine the accuracy and quality of current density images. Filtering, signal windowing, and magnetic shielding were adopted to suppress the equivalent magnetic noise [
14,
15,
16,
17], however, were inapplicable in reducing the common mode noises arising from background electrical equipment and earth magnetic fluctuation [
18,
19].
In this paper, a VCD imaging method is proposed to simultaneously image the VCD distribution based on tensorial magnetic signature, and reduce the common mode background noise using a second-order gradient technique [
20,
21]. A closed-form vectorial current density inversion method is established by combining the difference approximation of differentials and the second-order gradient method. Square, triangular, circle, and snake shapes of current-carrying wires are experimentally adopted as typical DUTs and were vectorially current-imaged using a tri-axial fluxgate sensor. FEA models of the above DUTs are correspondently developed to give coherent VCD images between numerical results and experimental results, with a relative error less than 10%. Although this work is implemented in 16AWG electrical wires, the VCD reconstruction method is fully applicable to other types of current carriers such as superconductor wires, silicon-based chips, and power batteries, etc. It is expected to be migrated to areas of semiconductor, electronic, and electrical energy industries as novel tools for design, characterization and health monitoring. VCD’s promising ability to identify current leakage will greatly facilitate the online inspection of electronic products.
2. The Closed-Form Vectorial Current Density Inversion Method
The magnetic field,
, generated by an electrical current in a conductor is governed by the well-known Ampere Law, of which the differential form can be written as:
where
is the curl operation, and
the current density vector determined by the cross-section and total current flow in the conductor.
By applying the linear constitution relationship of
B–H in air:
, the current density can then be reformulated as:
where
is the magnetic flux density,
the vacuum permeability, and
,
, and
the unit vector for a Cartesian coordinate system, respectively. Equation (2) indicates that the vectorial current density,
, is measurable as long as the spatial variation of the magnetic flux density,
, can be obtained. Applying the physical definition of the magnetic field gradient,
, to Equation (2), the vectorial current density,
, can be linearly expressed in a closed form as a function of the magnetic field gradient tensor elements, described as:
2.1. Measurement of the Magnetic Field Gradient
As shown in
Figure 2, the theoretically defined magnetic field gradient in Equation (3), mathematically also the spatial differential of magnetic flux density
B, can be experimentally obtained by spatial differencing
and
over the baseline (
), given by:
where
and
are common-mode noise-rejected magnetic flux density components at each sampling position obtained by point-by-point background noise (
) subtraction, described as:
2.2. Raw Imaging Data Formatting
Through the data measuring process shown in
Figure 2, we can get a near-field raw dataset, far-field raw dataset, and reference dataset.
Figure 3 illustrates the data flow and processing diagram for a vectorial current density imaging method based on the magnetic gradient tensor. First, a struct buffer is applied to store background noise as a reference raw dataset. Second, far-field and near-field tri-axial magnetic flux density datasets are differenced with respect to the reference raw dataset, respectively, to give a noise-rejected far-field and near-field dataset. Third, the second-order differencing is conducted with respect to baselines (
) to offer a magnetic field gradient dataset. Then, the VCD dataset can be attained based on Equation (3). Finally, by spatially fusing the VCD dataset and location dataset point by point, the vectorial image data of current density are precepted and recorded in an RGB format, where R, G, and B quantitively represent the
x,
y, and
z components of CD. It is practically necessary to mention that each component, as well as the vector strength, can be visualized through a pseudo-color image (e.g., rainbow or traffic cloud picture) with a combination of arrows.
3. Experimental Setup
Figure 4 shows the schematic diagram of the experiment setup for the closed-form vectorial current density inversion method. The experiment system (
Figure 5) consists of seven parts: the non-magnetic DUT substrate, metrology grating module, step-motor control module, host PC, DC current supply, data acquisition device, and a tri-axial fluxgate sensor. Four typical DUTs, namely square, triangular, circle, and snake shapes, made of 16AWG electrical wires (
Figure 6) are excited by the DC constant current of 5A fed by a power supply (RIGOL
® DP831A, Beijing, China). The DUTs are horizontally fixed on a PTFE substrate with a size of 300 mm × 300 mm. The motor control module (CH-Magnetoelectricity
® Technology KZ400C, Beijing, China) communicates with the host PC through an RS-232 protocol with synchronous location data measured via metrology grating. Triaxial magnetic field data at the reference, near, and far positions are precisely measured using a tri-axial fluxgate sensor (Bartington
® Mag-690, Oxon, UK), and recorded by a data acquisition device (CH-Magnetoelectricity Technology
® CH370, China). The imaging data processing program illustrated in
Figure 3 is embedded in a host PC to display the VCD picture. A spatial resolution of 0.01 mm within range of 400 mm, and a magnetic field resolution of 0.1 nT within a range of
100
μT were declared on the manufacture datasheet.
5. Results and Discussion
Figure 7 shows the raw near-field, noise-rejected near-field, and numerical magnetic flux density |
B| images @ z = 2 mm. The magnetic fields are generated by a DC-5A excited snake shape conductor to represent most DUTs without losing its generosity. We see from
Figure 7a that the vector superposition of magnetic fields from DUTs and background magnetic fields are measured to give a relatively high magnetic flux density |
B| up to 69.9
μT.
Figure 7b,c shows that the maximum magnetic flux density |
B| @ z = 2 mm is evaluated to be 41.7
μT, 41.7
μT, and 42.1
μT from Biot–Savart’s solution, numerical simulation and experimentally noise-rejected data, respectively. The results indicate that a highly coherent magnetic field results with a relative error less than 1% can be obtained using a differencing technique with respect to the reference raw dataset based on Equation (5).
Figure 8a,b shows the
x,
y, and
z components of experimentally noise-rejected near-field data and simulated near-field data in
Figure 7b,c correspondingly. Following the positive directions of the imbedded Cartesian coordinate, the feeding DC current from the bottom end of DUT results in obvious
above
y-guided wire,
above
x-guided wire, and
for all wires are governed by famous ‘right-hand rule’. In comparison with magnetic flux density |
B| in
Figure 7, the components of
B intuitively reflect the directions of CD, however, cannot directly offer current flow directions, especially when a non-orthogonal condition occurs between wire directions and
B directions. Benefiting from excellent noise rejection performance endowed using a differencing technique, highly coherent results are obtained between experiment and simulation results, as shown in
Figure 8a,b, respectively. It can be also found that the
image can roughly indicate the wire contours through its cross-zero isogram. This may be the fundamental principle of traditional Fourier Transform back-evolution method for MCI [
2]; however, it is limited in providing the magnitude of CD only.
Figure 9a,b shows the full tensor of magnetic field gradient images of experiments and simulation results, respectively. As the most critical intermediate variable, the full magnetic tensor images are derived from measurement results in
Figure 8 according to Equation (4). An intuitive impression of large zero areas (colored green in the rainbow bar) leaps to the eyes from the tensor images. This can be attributed to the effects of second-order operations, which has exempted magnetic features induced by non-local CD. Moreover, both experiment and simulation results exhibited a high coherence and strictly symmetric pattern (
), which are intrinsic properties of the magnetic tensor. Such property can be practically used as data check criteria. It is also interesting to find that the non-zero value of
or
appears on the corner. Further examination indicates the right-hand rotation CD results in positive
, and vice versa. Such a discovery is not applicable to
and
because the scanning is conducted parallel to the
x–y plane.
Figure 10a,b shows the reconstructed vectorial current density imaging results based on the measured and simulated magnetic gradient tensor, respectively, based on the closed form inversion formula Equation (3). A relative error of less than 3% is achieved between experiment results and simulation results for all shapes of DUTs. The CD magnitude and its directions are successfully imaged with obvious evidence to distinguish different shapes of DUTs. Through the closer examination of the experimental VCD image, the ingoing and outgoing imprints can be recognized as shown in
Figure 10a. The difference between the calculated CD under operating condition result (|
J| = 5 A/1.5 mm
2 3.33 A/mm
2) and reconstructed result |
J| = 3.5 mA/mm
2 can be explained through a probe-to-sample distance induced lift effect. That is, magnetic field attenuated third order inversely with respect to sensing distance, resulting in a smaller CD magnitude and diffusive CD images. The CD magnitude attenuation can be corrected by multiplying a fixing factor (
k = 1050 in this case); however, the diffusion effect can only be alleviated by reducing the probe-to-sample distance. The overall CD resolution is estimated to be 24.15 mA/mm
2 by selecting the maximum of the best resolutions among the four DUT CD experimental images, and by multiplying a fixing factor of 1050. It is worthy to mention that the CD resolution can be improved by simply reducing the probe to sample distance. The spatial resolution is estimated to be 0.69 mm with respect to the wire radius. It is expected to improve the spatial resolution by reducing the size of the sensor to image the CD in a printed circuit board or even a chip.