Nested kernel degeneration-based boundary element method solver for rapid computation of eddy current signals

In this article, the nested kernel degeneration (NKD) based integral equation fast solver is proposed to compress the dense system matrix generated by boundary element method (BEM) for calculating eddy current signals in 3D nondestructive evaluation (NDE) problems. The near and diagonal block interactions are computed and stored by full matrices. However, for the far block interactions, not like in the single level kernel degeneration-based method that all the interaction matrices are compressed by the interpolations, the nested approximation framework makes the expression of the higher-level matrices by the ones at leaf levels and the transfer matrices between neighboring levels. Different integral kernels are degenerated and nested to ensure the performance of NKD based method. The robustness and efficiency of the proposed solver are validated and verified by making numerical predictions and comparing them with the ones achieved by analytical method, numerical methods and experiments. Our numerical experiments show that the NKD solver is more efficient than the butterfly solver (multilevel adaptive cross approximation algorithm (MLACA)) for electrically small eddy current NDE problem. With the NKD based method, O(N) computational complexity is achieved for multiscale low frequency eddy current simulations without sacrificing accuracy, where N is the number of unknowns. So far, to our best knowledge, the proposed method with linear complexity is the most efficient integral equation-based BEM forward solver to predict the eddy current signals from cracks and flaws.