Determination of Magnetisation Conditions in a Double-Core Barkhausen Noise Measurement Set-up

The magnetic Barkhausen effect is useful for assessing 1D and 2D stress states of ferromagnetic steel objects. However, its extension to technically important materials, such as duplex anisotropic steels, remains challenging. The determination of magnetisation inside the studied object and the electromagnet for various geometries, materials and magnetisation angles is a key issue. Three-dimensional, dynamic finite element analysis has been applied to reproduce time-varying fields inside and outside the prototype of a double-core magnetising setup. Useful relationships between characteristics (peak height and location) and magnetic induction vector have been proposed. The qualitative plausibility of simulation has been validated with an experiment and an analytic formula of skin depth. The angular anisotropy of magnetic Barkhausen effect (MBN) in an isotropic sample has been shown in simulation and confirmed experimentally. The numerical model, despite some limitations, seems to be an efficient tool for calibrating stress/MBN relationships at least in isotropic structural steel components.


Introduction
Deducing surface stress state of the ferromagnetic steel structures is one of the most significant goals of the NDT community. The Barkhausen effect has been applied in this context for more than two decades [1][2][3]. Its applicability is easy to demonstrate in simple laboratory set-ups, where controlled uniaxial stress change causes significant monotonic variations of magnetic Barkhausen effect (MBN) [4]. Effect of biaxial stress has been intensively studied as well [5,6], and procedures exist to determine principal stresses from orthogonal MBN measurements, as long as a given steel grade has been properly calibrated in advance.
As in any NDT technique, the inverse, in-situ problem remains a challenge, for several reasons. Firstly, the magnetic field distribution and its time variation within the material cannot be measured directly, so there can be doubts about the repeatability and controllability of the measurements. Because of the trade-off between measurement speed and accuracy, the choice of optimal frequency, amplitude [7] and number of sampled points / probe angles remains open. Secondly, even if the field was well-defined, some uncertainties would remain as to the fundamental sources of MBN. Can its main characteristics (peak location, height and width) be derived from B, dB/dt? Thirdly, the material data, both magnetic and structural, are subject to uncertainty and variation with heat treatment and degradation [3,8]. Finally, there are technically important materials, like duplex steels, which exhibit very significant magnetic anisotropy. Although some work has been started [9], detailed stress assessment procedure is still missing.
The Authors of this paper have mainly dealt with the first and the second issue, using Finite Element modelling strategy which has proven useful and reliable in our previous works [10,11].
Although the application of a nonlinear transient magnetic FEA in the context of MBN-based NDT is our original, unique contribution, the concept of a double-core measurement set-up is not new. It has first been suggested and implemented by Vengrinovic and Tsukerman in 2004 [12], and used to produce co-called directional diagrams (DD), demonstrating significant sensitivity on plane stress state and texture of the tested material.
The magnetisation time-and space characteristics obtained by numerical simulation can provide an input to a feedback control system, thus allowing it to achieve arbitrary shapes and amplitudes of magnetic field waveforms inside the tested structure. The discussion of types and possible enhancements of feedback methods in magnetic circuits can be found in works by Zurek et al., Polik and Kuczman as well as White et al. [13][14][15].
Several abbreviations have been introduced in the article: BSUM resultant magnetic induction-most important output from the finite element analysis; calculated as square root of sum of squares of B x ,B y and B z components at a specified location, usually central point CP Central Point-critical location in the geometric centre of the magnetised area, where a ferrite MBN probe is placed CAx Central Axis-the vertical axis of symmetry of the Set-up, including CP DCMS double-core magnetising set-up-a prototype with four coils wrapped on two perpendicular C-cores (further referred to as X-Core, Y-Core). The applied sinusoidal excitation currents agree in phase, but may have different amplitudes. AX, AY Amplitudes of source current densities on X-Core and Y-Core, expressed in MA/m 2 SCM source current magnitude-useful convention for describing the excitation strength. Calculated as square root of sum of squares of AX and AY AMA actual magnetisation angle-derived from calculated results and/or Hall-probe measurementsrepresents the real induction vector orientation at the central point. NMA nominal magnetisation angle-directly derivable from AX and AY, defined as arctan (AY/AX); in trivial cases, such as AY = 0, NMA = AMA, due to perfect symmetry 2 Model Description

Finite Element Mesh
A series of finite element simulations has been performed in ANSYS software with electromagnetic functionality. The  dimensions of main components, namely X-core, Y-core and the plate have been taken from experimental prototype. The 20-node elements with edge-flux formulation have been applied. Hexahedral, refined discretisation has been defined within the plate and the cores, while the surrounding air space is filled with tetrahedral, relatively coarse elements. There is a gradient of edge length along the plate's depth, with thinnest elements (0.4 mm) at the top surface, and maximum thickness of about 1 mm at the bottom. Dirichlet (flux-parallel) boundary conditions were applied to all the walls of the 6sided hexahedral model space (Figs. 1, 2).

Material Model
The steel has constant electrical conductivity of 5 MS/m. The cores are perfectly laminated, i.e. they do not exhibit any eddy

Solution Scheme and Stability Check
The simulation is transient, and the magnetic induction distribution results from time-varying current density loading on four excitation coils. The amplitudes of current densities (max. of 0.8 MA/m 2 ) are based on the preliminary setting of the prototype. Two coils wound on a common core work in series, i.e. they have got identical electromotive forces (EMFs) which add up. The X-core and Y-core loadings are controlled independently, but there is no phase shift between them, i.e. the zero-current and maximum-current occurs at the same instant. A convention of phase shifts "ϕ" was adopted, according to which the zero-current corresponds to ϕ = 0, minimumcurrent to ϕ = −90 and maximum current to ϕ = +90. The variation of magnetic induction may be thus related to this reference angle.
During test runs, the convergence of the Newton-Raphson scheme was invariably poor at the first time-step. Success was achieved only when an extra static step was introduced, the allowable number of substeps was increased and the tolerance of solution error was temporarily loosened. Nevertheless, high magnetic permeabilities of order of 50,000 (corresponding e.g. to permalloy cores) still produced convergence failure.
The solution stability was another issue on which results accuracy depended (see Fig. 4). The second calculated period differed significantly from the first one, so the results had to be treated as numerical, unrealistic transients. The 3rd, 4th and 5th period turned out to be similar to the second one and could be treated as a stabilised enough, reliable solution.
The applied software ANSYS has no intrinsic B-H hysteresis feature, but it can represent eddy-current, 'apparent' hysteresis. This notion can be defined as a hysteretic dependency between the source excitation current and the magnetic induction within the magnetised object. The carbon steels typically exhibit intrinsic hysteresis of order of 100 A/m. The eddy-current, apparent coercivity H ce obtained from the simulation at 10 Hz was of order of 500 A/m. This value has been derived from the phase shift between the zero-crossings of source current and the magnetic induction, related to 90 • and multiplied by maximum value of magnetic field intensity, H max .  Figures 7 and 8 show the BSUM induction distribution in the topmost layer of the plate for NMA = 0 • ,atϕ = 0 (zero excitation) and 90 • (maximum excitation) respectively. The Y-core disturbs the field distribution even though its coils are unloaded all the time. It is argued, that the laminated core material leads to local decrease of flux density, providing an extra low-reluctance pathway for flux lines. Again it is interesting to note, that at zero-current phase, the induction values are far from being negligible, e.g. there are regions of B>1.0T.

Analysis of Time-Variation of Magnetic Induction Along the Central Axis (CAx)
The magnetic field time-dependency does not vary with angle when the single core is being mechanically located at different angles. Experience with DCMS prototype shows however some Simulation shows that it is harder to magnetise the sample at 0 • than at 45 • or 60 • . The effect may be accounted for by the saturation/nonlinearity of the core, which is the most significant when a single core is loaded with the entire RSC, i.e.at0 • .
Actual Magnetisation Angles (AMA) are calculated as 0 • , 46 • , and 48 • , for Nominal Magnetisation Angles of 0 • ,45 • and 60 • , respectively. The surprising discrepancy between AMA and NMA at 60 • may again be associated with the saturation effects. The core which is more intensely loaded saturates more significantly, so larger proportion of its flux is lost in form of stray fields, and the flux within the circuit is reduced due to lower effective magnetic permeability. Consequently the flux values penetrating the sample from both cores become closer. This results in BX/BY ratio at the Central Point (CP) tending to 1.0, and Actual Magnetisation Angle (AMA) approaching 45 • .
Although the first top layer generates most of the MBN, it is useful to study the B time-through-depth dependency. Activity of eddy currents within the deeper layers entails retardation of magnetic induction with regard to the source current. This phenomenon is termed as "eddy-current hysteresis". The conservation of Gauss theorem requires that the sum of fluxes in all the layers be approximately proportional to the source currents. In particular, the total flux must zeroout at ϕ = 0. Consequently, the surface layers exhibit field of opposite sign, appearing as "advanced in phase" (Fig. 10).
Analysis of the numerical data behind Fig. 10 indicates a non-trivial phenomenon of hysteresis loop reversal between surface (1) and deeper (especially 7+) layers. The intermediate (4th) layer exhibits unusual hysteresis loop with two self-crossings.

Through-Depth Field Decay
The distribution of magnetic flux inside a ferromagnetic bulk sample depends significantly on the source current frequency. The well-known skin-depth formula (1) has been found to be relevant in similar set-up [16], so it can be reapplied in the double-core set-up.
where ω − angular frequency µ r − relative magnetic permeability (dependent on local H) µ 0 − absolute magnetic permeability of the vacuum σ − electrical conductivity of the sample The skin-depth "d" corresponds to the distance from the surface at which the induction amplitude decreases by the factor of 2.71 ("e").   It is expected, that the magnetic induction B reaches highest values at the surface, and its amplitude (maximum value in time) decreases with depth. The simulation run at nominal 60 • with relatively "high" source currents (source current amplitude of 6 MA/m 2 ) has been selected to examine the through-depth induction variation and compare it against analytical formula.
The assumed B(H) curve leads to varying estimates of d as a function of local H. At H = 100, 1000 and 5000 A/m, d∼2, 6, 15 mm respectively.
The calculated induction amplitudes at varying depths are shown in Fig. 11. The decrease by factor ∼3 at 6 mm correlates well with the analytical prediction.

Comparison of Low-Current Versus High-Current Magnetisation
Two series including the total of 6 simulations have been performed. Three of these had the elevated source current magnitudes (SCM) of 6 MA/m 2 , three others had the SCM of 0.3 MA/m 2 . In both series the nominal magnetisation angle (NMA) was set to 0 • ,45 • and 60 • . The induction magnitude (BSUM) was recorded at the Central Point as a function of the phase of the magnetising current. The comparison at 60 • , for high SCM versus low SCM is shown in Fig. 12.
The low magnetising currents allow operating within the linear portion of the B(H) curve, producing induction magnitudes not exceeding 0.6 T, whereas the selected high SCM causes a very prominent non-linearity due to saturation of the cores and the central area of the plate, with BSUM approaching 2 T. In such a case, the time-variation of any component of the induction is no longer monotonic when the magnetising current increases from its negative peak value up to the positive peak value. This strong non-linearity and lack of monotonic behaviour make the setting unsuitable for interpretable Barkhausen Noise measurement.
The comparison of nominal magnetisation angle (NMA) and actual magnetisation angle (AMA) is shown in Table  1. The difference is negligible at NMA = 0 • , regardless the SCM. At 45 • the discrepancy, due to the inequality of lengths of the C-Cores, does not exceed 3 • , and still can be deemed negligible in practice. However, a significant difference between NMA and AMA (11 • ) is observed for NMA = 60 • , and high source current Magnitude. The induction vector at the central point deviates from expected 60 • towards the 45 • direction. It is argued, that one of the C-Cores penetrates deeper into the magnetic saturation region. As its differential magnetic reluctance increases, it generates more stray field, which does not enter the plate. Consequently, uneven saturation of the C-Cores leads to a decrease of difference between BX and BY at the Central Point of the plate. This, in turn, results in AMA approaching 45 • . Unequal air gaps at the poles in the actual system (e.g. operating over a curved surface) may play some role too, calling for an application of some flux control system, in order to balance the flux contributions from both cores.

Summary and Perspectives
A practically useful, 3D nonlinear dynamic Finite Element model of the double-core Barkhausen noise (DCBN) measuring setup has been constructed. The convergence of calculations has been achieved for the default operational frequency of 10 Hz.
Original definitions (NMA-nominal magnetisation angle, AMA-actual magnetisation angle, SCM-source current magnitude) have been proposed for the sake of unambiguous characterisation and repeatability of measuring conditions. Insight has been obtained into the time/space distribution of the magnetic induction within the studied plate as well as the C-Cores, which is essential for planning and interpreting any magnetic-field based NDT procedure. The following observations have been made: -The eddy current-related hysteresis loop is strongly dependent on the depth of the plate at which it is recorded. It exhibits a characteristic advancement in phase at the top surface, which turns into a phase lag at depths of approx. 2 mm and greater. The positive phase shift at the surface was demonstrated experimentally in a previous study [17]. -The spatial decay of the time-maximum of the magnetic induction is consistent with the standard skin-depth formula. This suggests, that this simple and practical formula can be reliably used in the initial prediction of magnetic field inside a DCBN set-up. -The AMA is acceptably close to the expected NMA provided that the set-up operates within a linear portion of the B(H) material characteristics. Otherwise there is a non-negligible dependency of all the studied magnetisation parameters on NMA. This effect is accounted for by the uneven saturation of both C-Cores. As expected, the maximum induction generated in the plate for a constant, high current magnitude, is higher at 45 • than at 0 • /90 • . -When comparing the low-current and high-current source magnitudes, the resulting magnetisation patterns have been shown to differ significantly, both quantitatively and qualitatively.
From a practical point of view, the presented simulations allow determination of an optimum workpoint of DCBN set-up, at which the magnetic induction in the plate is as high as possible without introducing uncontrollable distortions due to B(H) saturation.
Transforming the calculated B(x,y,z,t) function into numerical local MBN patterns, and comparing them against experimental results, is a pending task.