An effective and accurate technique for the numerical solution of 2‐D electromagnetic scattering problems with 3‐D sources is presented. This solution introduces a set of the usual boundary integral equations and uses a scalar Green’s function. In this scalar version, the unknowns of the problem are the boundary values of the longitudinal fields and their normal derivatives in the Fourier domain. A generalization of the usual boundary integral formulation enables us to handle a large class of models composed of piecewise homogeneous domains, including contiguous domains, multiply‐connected domains and unbounded domains. This formulation involves the solution of a system of linear equations, and results in a significant saving in computation time in comparison with other rigorous methods.
The requirements for the numerical implementation of this solution are described in detail. Numerical tests were carried out using the important example of electromagnetic tomography. The specific symmetry properties of the response function in this case are illustrated. Numerical accuracy is verified over a large frequency range, up to 1 MHz.