Three-Dimensional Fem-Bem Computation of Electromagnetic Responses of Flaws

Electromagnetic (EM) nondestructive testing of a lossy dielectric material involves exciting the material with an external EM source, measuring responses on a measuring surface, and inferring the location, shapes and properties of possible flaws. But the interaction of electromagnetic fields with the material is a highly complex process. Three-dimensional computation is therefore required to accurately predict the overall interaction of EM fields with the material and flaws. The problem is unbounded, and strong inhomogeneity exists inside the material, due to the presence of flaws. The coupling of the finite element method (FEM) and the boundary element method (BEM) provides an efficient tool for modeling this type of problems [1]. Two questions must be answered before a successful use of this method. i) The interior EM fields are governed by a weak integral form based on the curlcurl equation of a field variable, leaving the divergence free equation unspecified, and the interface and boundary conditions to be forced separately. Spurious solutions may occur, when the standard nodal based finite element method is applied. Therefore, special care should be taken to prevent the so-called vector parasites [2]. ii) The exterior EM fields are governed by a surface integral representation, which takes care of the infinite extent of the exterior domain. This integral equation invloves tangential surface currents. The standard nodal based boundary element method does not provide a rigorous description of the tangential surface currents especially near corners and edges, because of the ambiguity in normal directions. Semi-discontinuous superparametric elements can be used [3]. Also, discretization of the integral equation is central to couple the source fields into the materials and flaws