Nondestructive Testing of Metallic Cables Based on a Homogenized Model and Global Measurements

1 Research Group for Numerical Analysis and Mathematical Modelling, Department of Mathematical Analysis, Galglaan 2, 9000 Gent, Belgium 2 Department of Electrical Energy, Systems and Automation, Ghent University, Sint-Pietersnieuwstraat 41, 9000 Gent, Belgium 3 Department of Electrotechnology, Faculty of Applied Engineering Sciences, University College Ghent, Schoonmeersstraat 52, 9000 Gent, Belgium


Introduction
In 1 , a homogenized model is derived for nondestructive testing NDT of heterogeneous magnetic structures cables . In this contribution, we employ the model for actual NDT.
NDT can be based on eddy currents or on magnetic properties.
In case of eddy current testing, the imposed field has a sufficiently high frequency so that eddy currents are induced in the device under test. Cracks can be detected because they change the path of the eddy currents. The frequency of the excitation determines-via the penetration depth-the thickness of the scanned region: in case of high frequencies, mainly the surface is scanned; in case of lower frequencies, also cracks deeper inside the material will be detected. The paper 2 explains the principle of eddy current sensors and investigates the effect of inhomogeneity of the target material. In 3, 4 , computational techniques are given to solve the magnetic field problem in eddy current testing. Both papers also present eddy current testing experiments: Cardelli et al. consider benchmark problems of an aluminium plate with different cylindrical defects; Li et al. apply the eddy current testing technique to multilayer conductive structures.
In case of NDT based on magnetic properties, there are techniques that exploit the magnetic reluctance variation and techniques that measure magnetic flux leakage. Here, the frequency of excitation can be very low; even excitation with permanent magnets is possible. In 5, 6 , these two techniques for magnetic analysis of NDT are studied for ferromagnetic ropes. The cable is modelled in a 3D numerical model as one solid object with uniform properties. The individual strands are not modelled and-as the magnetic NDT in the cited paper is static-the considered static numerical model has no eddy currents.
The model of the ferromagnetic cable that we derived in 1 is useful for NDT of cables, both for excitation at low frequency where changes in magnetic properties are investigated and at higher frequency eddy current testing . It is valid for a wide range of amplitudes and frequencies.
The motivation of the work in this paper is to improve the performance of nondestructive testing techniques by exploiting information in the magnetic behaviour of the target material. We show that by solving an inverse problem based on the model and on measured hysteresis loops at several frequencies a broader range of defects in the cable can be detected.
In Section 1.2, we recall shortly the model derived in 1 . In Section 2, we introduce the problem we study, which can be classified as a 2D inverse spectral problem. We consider its reduction to 1D. Subsequently in Section 3 we study an inverse problem using the 1D model. It results in a simple two-parameter determination problem. We validate the method in a real data experiment in Section 4.

Detection Setup
Let us have a heterogeneous cable consisting of many thin strands. Usually, the strands are not isolated separately. A coating can be applied to the cable to avoid corrosion. We consider NDT of such a cable based on magnetic hysteresis and eddy currents. The detector has a very simple design. It is made of a single excitation coil and a single measurement coil, which are wound around the cable, see Figure 1. When an alternating current flows through the induction coil, the corresponding alternating magnetic field is induced. Consequently, eddy currents are induced in the cable.
Due to the heterogeneous nature of the problem, very fine space discretizations, for example, very fine meshes for finite element models FEM , are required to obtain a feasible accuracy of a model. A certain model reduction has to be applied to avoid extensive computational costs.

Model
The facts presented in this section are direct consequences of the results obtained in 1 .
Let the circular cross section of a cable be the computational domain Ω. The outer circular boundary of Ω is denoted by Γ E . We introduce into the cable an inclusion Ω I . For a multiplicity of exposition, we assume that Ω I is simply connected, that is, that we have only  Figure 1: a uncoated cable with 7 bundles of strands, b 2D scalar model, c experimental setup: a part of the excitation and the measurement winding around the coated cable. The coating is removed on the left of the figure to clearly show the individual strands, d cable where one bundle of strands is removed over a length of 1 cm. To make the defect visible, the coating is removed and the "defect" bundle of strands is folded in radial direction. Evidently, the method for NDT does not require the defect to be visible. one inclusion. Ω I is filled with nonconductive material air , that is, σ 0. The boundary of Ω I is denoted by Γ I . Let us define Γ : Γ E ∪ Γ I and Ω C : Ω \ Ω I . We suppose that Ω C is filled with a conductive material metal . Further, we assume that the inclusion is buried inside, that is, dist Otherwise it could be detectable by visual inspection, see Figure 2.
If the induction coil is sufficiently uniform, it induces the magnetic field H perpendicular to the cross section of the cable. Let e 1 , e 2 , e 3 be the Cartesian coordinate system and let an infinitely long cable be aligned along e 3 . We have H He 3 . 1.1 In 1 , we assumed that the diffusion of electromagnetic fields in the heterogeneous cable can be described by the Maxwell equations with periodically oscillating coefficients. An accurate and computationally efficient model can be then derived. The idea behind is to replace the heterogeneous material in the cross section by a fictitious homogeneous one, whose behaviour at the macroscopic level is a good approximation of the one of the composite material. We obtained such a homogenized asymptotic model employing the two-scale convergence. We validated the model on real steel cable data. In the isotropic case the model reads in the terms of H: where μ 0 is the permeability of the vacuum, ω is the angular frequency in rad/s, n is the normal outer unit vector with respect to Ω C and H E is the prescribed Dirichlet boundary condition on the exterior boundary Γ E . H I is an a priori unknown boundary condition on the cavity interior boundary Γ I . From physical reasons H I is a complex and H E a real constant. 1.5 is an nonlocal boundary condition on Γ I obtained considering Ampère's law in Ω C and the Faraday's law in Ω I and using the continuity of tangential components of the electric and magnetic field intensities, together with the fact that the conductivity is zero in Ω I . This extra condition is sufficient to determine H I . The homogenized conductivity σ h and homogenized permeability μ h are scalar, generally complex, constants. They are given by certain integral expressions on Y , the socalled Y -cell, from repetition of which Ω C is made. Let us recall here, that, we assume that the heterogeneous cable has a periodic structure. For details, see 1 . Let us remark that σ h is a scalar only in the isotropic case.
The homogenized system has an identical form as the original multiscale system. More importantly, we obtain exactly the system 1.2 -1.5 if we consider a homogeneous cable made from a material with the permeability μ h and the conductivity σ h , that is, μ h and σ h preserve their physical interpretation. Consequently, the results we obtain in this paper are valid for the homogeneous case as well.
If the periodic permeability and the conductivity are bounded and coercive, the homogenized μ h , σ h are bounded and coercive as well. We obtain a unique solution H ∈ H 1 Ω C to 1.2 -1.5 , see 1 .

Mathematical Problems in Engineering 5
We define the permeability in Ω by where μ 0 is the vacuum permeability, and the conductivity by in Ω I . 1.7

Inverse Spectral Problem
In this section, we study the problem of detection of inclusions buried inside the circular cross section of the cable as described in Section 1.2. We want to acquire as much information as possible about the inclusions based on the measurement for different frequencies.
The problem is closely related to the inverse conductivity problem 7 , where one wants to determine an inclusion Ω I , that is, the characteristic function χ Ω I , on the basis of Dirichlet-to-Neuman data on Γ E . Let us prescribe certain Dirichlet boundary condition on Γ E and then let us solve the governing elliptic equation. We can take the Neuman trace of the solution, if this one is regular enough. In this way we define the so-called Dirichletto-Neuman operator. The inverse conductivity problem lies in determining inclusions on the basis of one pair or more pairs of the corresponding Dirichlet-to-Neuman data. This problem is well studied, particularly in two-dimensions, see 7, 8 and references therein. In 8 the authors obtain uniqueness results in 2D for convex polygonal inclusions by two measurements.
Using the Dirichlet-to-Neuman approach for NDT of cables has two disadvantages. First, it is quite expensive to obtain local measurements of the magnetic field and the induction on Γ E : it requires to build quite complicated sensors for local magnetic field and induction measurements 9 , and the gathering of the measurements is a fairly slow process. However, it is the only way to obtain full topological information, that is, the support of an inclusion. Sometimes such a complete information is redundant, for example, it can only be important how big the inclusion is and how deep it is buried inside the cable. The second disadvantage of the Dirichlet-to-Neuman approach is that it does not exploit symmetry.
Let the macroscopic electromagnetic behaviour be described by 1.2 -1.5 , that is, the homogenized model which was validated in 1 . We can define to the solution H to 1.2 -1.5 the corresponding operator H which maps to ω and Ω I this solution, that is, H ω, Ω I : We cannot reach the solution H directly, we can only measure some global quantity dependent on H, for example, the magnetic energy. Consequently, we consider a certain The question is what information can be obtained over Ω I from the range R F of functional F for ω ∈ 0, ∞ . This question is very general and vague. The answer is certainly dependent on functional F, topological restrictions of τ and the governing equation.
We categorize the problem as an inverse spectral problem ISP 7 . The question which ISP asks is if a domain Ω can be determined by the eigenvalues λ resonance frequencies of an elliptic PDE with a homogeneous Dirichlet boundary condition on ∂Ω. The problem, we study does not strictly belong to this category. We do not consider only the eigenvalues but ω ∈ 0, ∞ . We do not try to recover the domain Ω itself, but only its inclusion Ω I . Also the boundary conditions are nonhomogeneous, and the solution has to fulfill 1.5 . The data are not the eigenvalues but certain global measurements F H ω, Ω I . But the main goal is the same: to determine the geometry based on measurements for different frequencies.

Reduction of the Model to 1D
We will not solve the problem in two dimensions, but we will exploit the symmetry of the circular cross section and obtain a one-dimensional problem in radial direction. The reduction will be quite accurate. We will show numerically, that almost all the 2D spectral information can be described by the resulting 1D model. Let us assume that the total magnetic flux through the cross section F ω, Ω I : is the global information available. This is one of the reasonable and physically relevant choices. Another one is the total magnetic energy inside the cross section. First, the volume of the inclusion |Ω I | is uniquely determined by a single measurement for ω 0. Then H H E is the unique solution to 1.2 -1.5 . The permeability inside Ω I is μ 0 and μ μ h / μ 0 in Ω C . Consequently, |Ω I | is indeed determined by F 0, Ω I . It it not possible to obtain measurements for ω 0. By measurements for ω 0, we mean measurements for sufficiently low frequency.
Second, since we hold a constant magnetic field H E on Γ E , we cannot distinguish two identical inclusions with different angular positions buried at the same radial position depth , see Figure 3 a .
Third, the penetration depth of the magnetic field decreases if ω increases. Intuitively, inclusions deep inside can be identified based on low-frequency measurements. Inclusions close to Γ E will alter F ω, Ω I also for high ω. Based on the above heuristic, only the distribution of the "mass" along the radial direction is important. It seems, that there should exist a 1D model, which approximates the 2D case quite accurately.
We employ a simple idea to replace Ω I by an auxiliary radial symmetric inclusion-a ring R I , see Figure 4. Let a be the inner and b be the outer diameter of the ring inclusion. We ask |Ω I | |R I | to be able to copy the case ω 0 and that the mass is buried at the same depth. The physical properties of R I , that is, σ 0, μ μ 0 are the same as for Ω I , that is, we replace Ω I by an isolating ring. In such a way, Ω I is described by only two-parameters, its surface and its depth, that is, the diameter of the ring s a b /2. Another possible approach is to try to derive a model based on averaged physical properties along the radial direction. One could integrate the governing equations in polar coordinates, make some simplifying assumptions on the solution H r, θ as ∂H r, θ /∂θ 0 , and attempt to derive an approximative model using only averages in angular direction. The approach we use is more simple and illustrative. Equation 1.2 in polar coordinates reads Recalling that we have constant Dirichlet boundary condition H H E on Γ E and that the geometry is as depicted in Figure 4 b we immediately argue that H is function of r only. Let 8 Mathematical Problems in Engineering c be the radius of Ω and let R I { r, θ | a < r < b, 0 < θ ≤ 2π}. Then 1.2 -1.5 becomes 1 r dH dr The where J 0 is the Bessel function of the first kind of order 0 and K 0 is the modified Bessel function of the second kind of order 0. Condition 2.6 was used to rule out the fundamental solution K 0 on interval 0, a , since d/dz K 0 z −K 1 z and K 1 z has a singularity in z 0. The constants c 1 , c 2 , and c 3 are determined from 2.4 and 2.5 . In 2.5 , H I is determined from 2.7 , see the proof of Theorem 3.1 in Section 3.1. Now, we compare numerically how well the 1D model 2.3 -2.7 approximates the results of the 2D model for different inclusions Ω I from the point of view of F H ω, Ω I . Real measurements, which are at least two-dimensional, have to be reproducible by the model. Good approximative properties are very important.
The 2D problem is solved by the finite element method. A very fine mesh is used, thus the solution is quasiprecise. To solve the resulting linear problem, a Cholesky solver with accuracy 10 −12 is employed. F ω, Ω I for ω 2π2 i , i 1, . . . , 14 is computed. We compare the resulting F ω, Ω I with the corresponding results of 1D model with the ring inclusion, that is, F ω, R I . First, let Ω I to be a circular inclusion. This case is the most probable in real world applications. It represents a single strand of a cable. Figure 5 shows that visually the results of 1D and 2D model are indistinguishable. The relative error of F ω, Ω I is only 0.724%. Let us note, that the parameters a, b, c, H E , σ h and μ h are here the same as by Figure 4 c , see the corresponding text above.
In Figure 6, there are results for square and box inclusions. For the square inclusion, the relative error of F ω, Ω I is 1.217%. For the box inclusion, the relative error becomes 4.55%. The mass of the box is not concentrated around a point but is more spread along the radial direction. In this case a two-parameter approximation by a ring becomes questionable. However, even in this case the representation error is smaller than the modelling error, see which are computed so that the masses of the object and of the corresponding ring are equal and that a b /2 coincide with the center of mass of the object.

One-Dimensional Inverse Problem
In the previous section, we have shown that quite general inclusions can be reasonably well approximated by a ring inclusion. It means if one has real world data, these can be described by a very simple explicitly solvable model 2.3 -2.7 with only two-parameters a and b. In the case of two disjunct inclusions, we could consider two rings. Let us suppose that a, b are unknown parameters, that is, we do not know if there is an inclusion and if yes where the inclusion is buried. We study the following inverse problem. Re(F(box)) Re(F(ring)) Im(F(box)) Im(F(ring)) b Figure 6: Fingerprints of 1D model for a square inclusion −1.217% relative error of F ω ; b box inclusion −4.55% relative error of F ω .

Uniqueness
Theorem 3.1. The ring inclusion, that is, a ∈ Ê, b ∈ Ê, such that a > 0, c > b > a can be uniquely determined by two measurements of Ê F ω, R I for ω 0 and for ω , where is sufficiently small.
Proof. The proof is quite technical, but the main idea is simple. The area S |R I | of the ring inclusion is uniquely determined by a single measurement for ω 0. Thus it is sufficient to obtain the uniqueness for one of the parameters a, b, say a. Then b S/π a 2 . Consequently we can consider F ω, R I to be a function of a only. We would like to show that if Ê F ω, a a 1 Ê F ω, a a 2 for an ω > 0, then a 1 a 2 . This seems to be, working directly with Bessel functions, quite difficult.
The idea is to exploit the asymptotic behaviour of Bessel functions when ω → 0. We consider the Taylor expansion of F ω, a with respect to ω at ω 0, that is, Mathematical Problems in Engineering

11
The F ω, a is an analytic function in ω, since it is an integral of a linear combination of Bessel functions, which are analytical except for the origin in the case of K 0 10 . The coefficients c i , i 0, . . . , 2 of the linear combination are dependent on ω via 2.7 , but are also analytical. We can indeed use 3.2 . Now, we evaluate F 0, a and ∂F 0, a /∂ω from 3.2 . We obtain immediately that To evaluate ∂F 0, a /∂ω we employ the asymptotic behaviour of Bessel functions. We have

3.4
The solution H to 2. We differentiate the right-hand side of 3.6 with respect to a and take the limit for ω → 0. By a tedious computation we obtain We obtained that for ω small enough we can approximate ∂Ê F ω, a /∂a by a simple expression. This expression is always negative, since 2b 2 ln c/b b 2 − c 2 is always negative for b < c. To see this just we use substitution x b/c. We obtain where −2 ln x x 2 − 1 is negative for x < 1. We obtained that Ê F ω, a is decreasing in a, which gives the uniqueness.

Remark 3.2.
One of the parameters in the above proof is identified on the basis of the measurement at zero frequency. Let us note, that the second one could be roughly identified using the skin depth as long as a defect much deeper than ρ cannot be detected. The higher the frequency, the smaller the skin depth. We have the following estimates based on the asymptotic expansions of Bessel functions of a very large argument see, e.g., 10 :

3.12
where ν ∈ Ê is the order of the Bessel function, α is any positive real number and |z| is sufficiently large. Thus we have that H, given by 2.8 , relatively quickly drops in amplitude. The parameter k in the argument of the Bessel function in 2.8 is inversely proportional to the skin depth: |k| √ 2ρ −1 . Consequently, when using a sufficiently high angular frequency ω, skin depth ρ is small, the field penetrate not deep enough and we are not able to see a difference in the flux F for a cable with and without an inclusion as long as this inclusion is buried much deeper than ρ. This can be observed in Figure 5.
To confirm numerically Theorem 3.1 we computed F ω, a for small f 0.1 Hz. The inclusion is given by a 0.0010 m and b 0.0015 m. F ω, a , depicted in Figure 7 a , is indeed a decreasing function of a. Let us define a least square problem by the following functional: 3.13 where F i , i 1, . . . , n, are measurements of the real part of the magnetic induction through the cross section of the cable for different frequencies f i , i 1, . . . , n. In Figure 7 b we depicted l a corresponding to F ω, a of Figure 7 a . We have two exact measurements for f 0 0 Hz and f 1 0.1 Hz, that is, n 2. We see that l a really achieves its minimum at a 0.0010 m.

Stability Analysis
Relation 3.9 from the proof of Theorem 3.1 could be used to analyze the sensitivity of Problem 1 for very low frequencies. But since the problem is ill-posed, one usually tries to sample for a very broad range of frequencies to ensure stability. We will analyze the stability of Problem 1 numerically.
To calculate a and b on the basis of measurements we use a hybrid minimization algorithm. First we run a global search on a rough grid. Then we employ the Levenberg-Marquardt 11, 12 algorithm to find the local minimum. We implemented the problem in C using 13 . The global search on a rough grid is computationally feasible since the evaluation of the explicit function 3.1 takes only approximately 0.3 ms on a regular PC Intel Core2 Duo CPU T7250 @ 2.00 GHz .
First, we will consider only two measurements, one for f 0 0 Hz and one for a variable frequency f 1 . We want to see how the choice of f 1 influences the accuracy. We consider uniform distributed noise of 1%. In the global search we take a uniformly distributed grid of 20 × 20 points. The maximum number of iterations for the Levenberg-Marquardt algorithm is set to 100. We compute the relative standard deviation s for samples of N 1000 computations, that is, for example for a The first observation is that only the measurements for a quite narrow interval of frequencies yield the best possible accuracy. The second observation is that for two measurements with 1% level of noise s never falls under 5% level, that is, two measurements are certainly not enough. And the third observation is that the inclusions buried deeper inside, that is, a 0.0006 m, are recovered less accurately than the inclusions close to the surface. Now, let us take a fixed frequency f 500 Hz, which lies in the optimal interval according to the above experiment. We consider the sensitivity of the inverse problem with two measurements for different levels of noise. The global search is performed on a uniformly distributed grid of 50 × 50 points. The relative standard deviation s for recovered parameters, computed for samples of N 100, is depicted in Figure 8 b . The results again confirm that the inclusions buried deep inside are much more sensitive to the noise in the measurements than those close to the surface. Even for very low levels of noise is the accuracy of the recovery of a deep buried inclusion quite poor.
Let us analyze how the number of measurements influences the accuracy. We considered N 2, 4, 6, 8 measurements where f 0 0 Hz in all the cases and f i 500 100 i − N/2 Hz, i 1, . . . , N − 1, that is, the measurements are regularly distributed around "optimal" value f 500 Hz. We again sample 100 times and we also take a uniformly distributed grid 50 × 50 for the global search as before. Figure 9 depicts standard deviations s of a and of b − a for different numbers of measurements as a function of noise. Considering more measurements improves the accuracy significantly. For N 8 and the noise level of 5% which is our expected modeling error, see Section 4.3 of 1 , s a is around 10% and s b − a even drops under 10%. Slightly better accuracy for b − a with respect to a is due to the fact that b − a is computed based on all the measurements, but the position a is indifferent to the measurement f 0 0 Hz. Figure 10 depicts the averages a and b − a. Let us recall we sample 100 times. The error of recovery lies for N 8 within 10% error band for both a and b − a. Notice, that the corresponding s a and s b − a are around 0.3. We see that the averages a and b − a are biased estimators of a and b − a, respectively. We assume that the reason is of algorithmic nature.

Real Data Experiment
We present a real case study to confirm the usefulness of the proposed inclusion detection method.

Experiments
The experimental setup is identical to the setup from 1 . The considered cable has 7 bundles of strands, each bundle consists of 19 wires with a diameter of 0.3 mm: Figure 1 a . The inspection of the cross section of the cable shows that the center s of all bundles except the central one is situated at a radial distance of approximately 1.6 mm.
To obtain measurements for validating the model, three ring core sensors for hysteresis loops were made by winding an excitation and a measurement winding around the material to test.
i For the first ring core, the material to test is a set of 18 strands that were electrically insulated by giving them an acryl coating.
ii For the second ring core, the material to test is a piece of undamaged coated cable: Figure 1 c .
iii For the third ring core, the material to test is the same cable where one bundle of strands is removed over 10 mm length: Figure 1 d .
To make the ring cores, firstly the strands for the first ring core or the cable for the second and the third ring core were folded to a ring toroid with 300 mm circumference. A ferrite core was added around the connection of the cable ends to close the magnetic lines. This simulates an infinitely long cable. Then, an excitation and a measurement winding were added, uniformly distributed along the circumference. Finally, an isolating tape was added: see Figures 1 c and 1 d . Let us note that it is also possible to make a straight sensor that can be moved over a cable to continuously scan new sections of a long cable.
We measured on each ring core on the one hand 40 quasistatic BH-loops 0.5 Hz between 0 and 5150 A/m, and on the other hand dynamic BH-loops and corresponding losses for four levels of magnetic induction 50 mT, 100 mT, 200 mT, 500 mT and for frequencies f of 50, 100, 200, 300, 500, 1000, 2000, 3000, and 5000 Hz. Generating the current for the excitation winding is done by a National Instruments PCI-6110 data acquisition card in combination with a Kepco 50 V/8 A linear amplifier. The voltage waveform induced in the measurement coil-proportional to dφ/dt-is integrated via an analog integrator. For each induction level and for each frequency, this integrated waveform obtained from the measurement coil is sampled simultaneously with the enforced current waveform in the excitation coil. The sampling is done by the same data acquisition card at a sampling rate of 200 samples per period for the loops at 3 and 5 kHz, and at a rate of 1000 samples per period for the other frequencies 0.05 to 2 Megasamples per second . In the LabVIEW software, the integrated waveform of the measurement coil is converted to an induction waveform B t and the current waveform is converted to the externally applied field waveform H t via Ampère's law. This results in the dynamic BH hysteresis loops, used by the inverse problem. The LabVIEW program has a feedback control algorithm on the waveform shape of the induction: it adjusts the excitation current waveform in order to obtain a sinusoidal induction waveform.

Identification of the Parameters for the Homogenized Models
The homogenized models used in the inverse problem, see Section 4.3, require the homogenized permeability μ h , and the homogenized conductivity σ h .
The homogenized permeability is a function of the magnetic induction. It is found via quasistatic 0.5 Hz hysteresis loop measurements on ring core 1, that consists of individually isolated strands. The modelling of the hysteresis loops is done using a complex magnetic permeability. In a time harmonic model, this complex permeability produces ellipses in the BH-plane. The determination of the complex permeability for a given frequency and amplitude is done by a least squares fitting of the ellipses on the measured quasistatic hysteresis loops. Evidently, the ellipses are only an approximation of the experimental hysteresis loops: the deviation between the ellipses and the measured loops is small if the material is not saturated Figure 11 a ; the shape of the experimental hysteresis loops deviates more from the fitted ellipses if the material is saturated Figure 11 b . Once the complex permeability is found for a broad range of induction values, the homogenized permeability is found using the homogenization result-Theorem 3.1 in 1 . We obtained μ h 50 mT 8.43799 · 10 −5 , −9.72202 · 10 −6 , μ h 100 mT 0.00010495, −2.79189 · 10 −5 , μ h 200 mT 0.000160883, −9.77414 · 10 −5 , and μ h 500 mT 0.000318907, −0.00028671 . The homogenized conductivity σ h is found for each of the models via an inverse problem with as input the measured dynamic loops on ring core 2 the cable without an inclusion , as explained in 1 . See Section 4.3 for the values of σ h for a particular model from that section.
The physical electrical conductivity not the homogenized conductivity σ h is found to be 4.1 MS/m by enforcing a DC current in a strand and by measuring the voltage over a well-known length of the strand.

Identification of Defects by Solving the Inverse Problem
In 1 we derived three successive versions of the homogenized model. All the models incorporate magnetic hysteresis via the complex homogenized permeability μ h .
In the case of Model 1, σ h is a real function of the induction level. Model 2 considers σ h to be a complex function of the induction level. Consequently, it can describe electric hysteresis effects as well. The most accurate model, Model 3, assumes that σ h is a complex function  are only functions of the induction level. The models achieve accuracies around 15%, 10%, and 5%, respectively. In Tables 1, 2, and 3 we present the homogenized σ h for all the models as obtained in 1 . The achieved modeling errors are therein presented as well.
We test how the models can comply with a real data inclusion detection problem described above. The results are summarized in Table 4. Two numbers are presented: i N-the number of wires detected defined by N : where 0.00015 m is the radius of the small wire and χ is the fill factor, that is, how much from the cable area is occupied by the wires. We have We expect that the optimal values are approximately N 19 and s 1.6 mm. Each row of Table 4 represents a different set of measurements which is taken into account.   We have altogether 40 measurements for 4 different induction levels. We consider all the measurements N 40 , low frequency measurements for f ∈ 0.5, 300 Hz N 20 , middle frequencies for f ∈ 100, 1000 Hz N 20 and only the high frequency measurements for f ∈ 500, 5000 Hz N 20 . From the results of the first two models, that is, Model 1 and Model 2, we conclude that the measurements for the complete spectrum of frequencies have to be included to obtain accurate results. Taking all the measurements yields very good agreement with the expected values. Model 1 found that 23.03 small wires were removed from the radial position r 1.5979 mm. Model 2 found that 17.90 small wires were removed from the radial position r 1.5897 mm. There is very good agreement of the recovery between the models, especially of the position.
Model 2 recovers the number of removed wires more accurately. It suggests that this model is better in predicting low frequency behaviour of the problem. The low frequencies magnetic effects are responsible for accurate recovery of the removed mass.
The results of the most accurate Model 3 are not consistent even for the whole set of the measurements, see Table 4 a . Based on the good results of the other two models, we can conclude that there are two possible reasons for the behaviour. First, probably 4 parameters per induction level are too many to fit. Let us note, the problem of finding σ h 1 and σ h 2 is illposed. A sufficient number of measurements and/or some kind of regularization are needed to get feasible σ h 1 and σ h 2 . Second, the results of the model from Figure 10 in 1 show that the losses are not monotonically decreasing, which we think is not physical. One probably has to enforce some constraints on σ h 1 and σ h 2 to correct the model.

Conclusions
We considered a problem of detection of inclusions in a metallic cable using a homogenized model derived in 1 and global magnetic hysteresis loops measurements for both low and high frequency up to 5000 Hz . We reduced the original 2D problem to a one-dimensional problem by symmetrization of a possible inclusion to a ring. We showed numerically that this reduction is rather accurate with respect to the flux through the cross section. At the same time, the inclusion is described by only two-parameters. The resulting 1D model has a piecewise analytical solution, which allowed us to evaluate explicitly the flux through the cross section of the cable. Thus the model is very fast.
We studied an inverse problem of detection of the ring inclusion. We obtained a uniqueness result, studied sensibility and considered a real data experiment.
The results show that-despite of its simplicity-the model is applicable in real world conditions. The most accurate adaptation of the model achieves an accuracy around 5% when estimating the number of defect strands and the radial position of the defect.