Abstract
This paper presents numerical simulation results of eddy current signals to determine the appropriated dimensions of imitative stress-corrosion cracks (SCCs) in viewpoint of similar eddy current testing signals. The imitative SCCs were modelled as multiple tiny slits in three dimensions. A Finite Element and Boundary Element hybrid method (FEM-BEM) was utilized to simulate the eddy current signals in an impedance form of an absolute pancake coil. Equivalent conductivity (
Introduction
Stress corrosion crack (SCC) is an essential service failure in engineering materials of the components and structures. The SCC occurs by environmentally induced crack propagation [1, 2, 3]. It can form as only a single crack or a group of colonies depended on many factors such as material, shape, type of stress, pH environment and type of corrosion.
Evaluation of SCC by using the eddy current testing (ECT) is an important issue in nuclear power plants, petrochemical and aerospace industries [4, 5, 6]. Although the ECT has high performance to treat the surface breaking SCC in length, it is still difficult to quantify the actual profiles with the measured signals. At the presents, SCCs have been researched and simulated to be conductive cracks in the viewpoint of ECT. There were some reports to explore the equivalent conductivity inside the SCC regions by numerical simulation of the uniform conductive cracks [7, 8, 9]. The results indicated that there was no correlation between conductivity and SCC width, but the SCC size related to its equivalent electrical resistance. The equivalent conductivity of SCC was applied to reconstruct the crack profile using inversion methods by authors [10, 11]. There is also a report of the local conductivity determination in the SCC region by using tested signals of the direct current potential drop (DCPD) [12]. The SCC conductivity showed as a linear distribution along the crack depth. In addition, imitative SCCs in perspective of ECT using multiple slits fabricated by different techniques were presented [13, 1, 4, 15]. They could be corresponding to the conductive distribution of the actual SCCs by introducing the conductive bridge between tiny slits.
As in literatures, the imitative SCCs are not exactly identical form to the actual SCCs, but they provide the same response to ECT signal. In order to determine the appropriated dimensions of imitative SCCs for simulating practical SCCs, this paper studied on the characteristics of the ECT signals due to different sizes of the imitative SCCs through numerical analyses. The specimens of imitative SCCs were modelled as a plate with tiny rectangular slits by varying the dimensions of gap between the tiny slits, length of the tiny slits and width of the tiny slits respectively. An absolute pancake coil was utilized in the numerical simulation as an ECT probe with different testing frequencies. Based on the numerical results, the equivalent conductivity (
Numerical simulations
Basic scheme of FEM-BEM
To predict impedance of ECT probe, the solution of the magnetic vector potential (
The air region can be governed as in Eq. (3), where
Based on the Galerkin’s discretization procedure, the following system of linear equations can be deduced to solve the solution of vector potential
where,
As shown in Fig. 1a, the imitative SCCs fabricated using the diffusion bounding technique were modelled using the multiple tiny slits separated from each other to allow the eddy current pass through the crack region. An austenitic stainless steel 316 L plate was chosen as the specimen for 3D simulation with size of 20
Models for simulations. (a) Imitative SCC model for simulation and definition of variable, (b) the model of uniform conductive notch with its dimensions.
Details of the imitative SCC dimensions
Referring to works of Yusa et al. [8, 13], it is appropriate to simulate the ECT signal of the SCC using a notch with uniform conductive region. Thus, the equivalent modeling utilized in this paper was chosen as a uniform conductive notch of width, length and depth being 0.2 mm, 20 mm and 3 mm respectively. This model was used to estimate the equivalent conductivity (
For FEM mesh of the models, both of the imitative SCCs and the uniform conductive notch were subdivided into a unity mesh. FEM and BEM nodes and elements were set covering the volume and the surface of the model. Dense elements were specified in the flaw zone with 0.2 mm side length. The numbers of elements of the FEM and BEM were 5,512 and 3,708 respectively. Figure 2 shows details of meshes for numerical simulation.
Meshing details of the imitative SCCs and the uniform conductive notch models.
Relations between the imitative SCC dimensions and the equivalent conductivity were investigated by simulation of the absolute pancake coil with a diameter of 3.2 mm and height of 0.8 mm as the ECT probe. The turn number of the coil was set to 140 turns. Total excitation current was set to 1 A, and the probe was discretized to 10
Validation of the imitative SCC model
An ECT experiment was conducted to validate the simulation results by using the tiny slit dimension of 0.4 mm, 0.8 mm and 0.2 mm in width length and gap between tiny slit respectively. They were fabricated to a thickness of 10 mm of SS316L plate by using the diffusion bonding technique as illustrated in Fig. 3a. The tolerances of tiny slit sizes were
(a) The imitative SCC dimensions of tiny slit width 
The average impedance signals obtained from the experimental data are compared to the simulation results as shown in Fig. 4. The impedances were plotted between resistance (R) and inductive reactance (XL) of the utilized ECT probe with error bars of standard deviation. The correlation between the experiment and the numerical results were in good agreement for all tested frequencies. The phases of signals could be similarly observed from both of the experimental and simulation results. The errors of the peak signal amplitudes were less than 3% of all cases.
Comparison of ECT signals between numerical results and experiment results of the imitative SCC dimension of the tiny slit width 
The simulation results of ECT signals corresponding to imitative SCCs are categorized to the impedance signals of varying gap between tiny slits, varying length of tiny slits and varying width of tiny slits. Then, the equivalent conductivities are estimated for different imitative SCCs. In summary, the prediction results of appropriated imitative SCC dimensions are established based on the specified equivalent conductivity for the ECT probe utilized.
Numerical results of varying gap between tiny slits
The ECT impedance amplitudes and phases affected by the varying gap between tiny slits at the constant tiny slit width
ECT impedance amplitudes and phases of varying gap between tiny slits at constant width of tiny slit 0.4 mm and length of tiny slit 0.8 mm for testing frequencies at 50 kHz, 100 kHz and 300 kHz.
The results of ECT impedance signals by varying length of tiny slits are shown in Fig. 6 by keeping constant width of tiny slit and gap between tiny slits at 0.4 mm and 0.2 mm respectively. The impedance amplitudes rose up owing to the increasing slit length from 0.4 mm to 1.2 mm for all of test frequencies in logarithm scale. Afterward, only small amplitude changes could be indicated for varying the tiny slit length from 1.6 mm to 2.0 mm. These results were caused from the spatial resolution of the ECT probe utilized. For the impedance phase results, there were small phase angle changes of all of test frequencies. However, shifting in phase of simulated frequency at 300 kHz was higher than those other frequencies. The results of further tiny slit widths and gap between tiny slits showed the similar tendency.
ECT impedance amplitudes and phases of varying length of tiny slit at constant width of tiny slit 
As shown in Fig. 7, the amplitudes of ECT impedance signals linearly increased in logarithm scale when the tiny slit widths were enlarged from 0.2 mm to 1.0 mm for all of test frequencies. These results were obtained from fixing the length of tiny slit at 0.8 mm and the gap between tiny slits at 0.2 mm. On the other hand, the impedance phases tended to decrease for all of the test frequencies. There was a significant phase decline in testing frequency of 300 kHz. The simulation results of other tiny slit lengths and the gap between tiny slits showed in the same trend.
ECT impedance amplitudes and phases of varying width of tiny slit at constant length of tiny slit 
Correlations between the absolute values of impedance (
Approximated coefficients of correlation between impedance (Z) and the equivalent conductivity
(a) Curve fitting between the impedance (Z) and the equivalent conductivity from the uniform notch at width of 0.2 mm, (b) Residual from fitting.
Figure 9 shows the calculated
The equivalent conductivity obtained from (a) varying gap between tiny slits, (b) varying length of tiny slit, (c) varying width of tiny slit.
The results of varying length of tiny slit are shown in Fig. 9b by setting the constant width of tiny slits and the gap between tiny slits as 0.4 mm and 0.2 mm respectively. The value of
From the estimation results, it could be indicated that there were small variations in the equivalent conductivity for frequencies of 50 kHz, 100 kHz and 300 kHz due to varying gap between tiny slits
In order to design the imitative SCC dimensions, statistical methodology was utilized to treat the significant variables affecting the equivalent conductivity including frequency and dimensions of slits. All variables were primarily analyzed by using the analysis of variance (ANOVA). The results could be indicated that there was no significant effect to the response of
Estimated parameters of response surface at specified different width
Estimated parameters of response surface at specified different width
Contour plots for prediction the length of tiny slit and gap between tiny slit for the imitative SCC at fixed width of tiny slit to (a) 0.2 mm, (b) 0.4 mm, (c) 0.6 mm.
As RSM scheme, the equivalent conductivity was fitted by a full quadratic model in terms of the gap between tiny slits (
In practice, the width of tiny slit was specified to design the imitative SCCs dimensions as given equivalent conductivity. Four typical observations of the values, for example cases, 20%, 15%, 10% and 5% of the base material, were chosen as target values. By setting the specified widths of tiny slit, the optimized RSM results of the gap between tiny slits and length of tiny slits could be determined and presented as in Table 4. These dimensions also gave the accurate prediction of the
Prediction results of the imitative SCC dimensions at specified the equivalent conductivity
This work studied the characteristics of ECT signals due to imitative SCCs of different dimensions through numerical simulation by using an FEM-BEM code. The SS316L specimens were modelled as a plate with multiple tiny slits. An absolute pancake ECT probe was utilized in the numerical simulation. The ECT signals of different excitation frequencies were simulated by scanning the probe along the length of the crack region with constant lift off, and scanning steps. The simulation results from the conductive notch model width of 0.2 mm were used to estimate the equivalent conductivity. The dependences of the
Increasing gap between tiny slits decreased the impedance amplitude of the ECT signals. On the other hand, increasing tiny slit length and width led to increase the impedance amplitude of the ECT signals. Phase of the ECT signals changes significantly at a higher frequency. Therefore, the test frequency of the pancake probe should be set to produce the eddy current penetration deeper than the crack depth to avoid the effect of phase shifting. The significant dimensions of the imitative SCC should be smaller than the utilized ECT probe to prevent the effects from the phase shifting and spatial resolution of the probe. Following the simulation results obtained by an ECT probe with a diameter of 3.2 mm and height of 0.8 mm, the gap between tiny slits and width of tiny slit should be set The maximum amplitude of ECT impedance signals could be used effectively to estimate the equivalent conductivity. In addition, there was no significant effect from changing the frequencies on the value. The statistical methodology is efficient to predict the appropriated dimensions of imitative SCCs, which is important to design an imitative SCC specimens, and helpful to SCC profile reconstruction through inverse analysis of ECT signals.
Footnotes
Acknowledgments
The authors would like to thank the National Magnetic Confinement Fusion Program of China (2013GB113005), National Science Foundation of China (No. 51577139 and 51407132) for funding this study.
