A particular application of the Volume Integral Method (VIM) dedicated to the efficient simulation of eddy current testing configurations involving several objects with very different dimensions is presented in this paper. Other problems addressed by this development are the modeling of flaws with complex shapes and the efficient computation of interactions between remote objects. Simulation results obtained with this direct model have been processed by the statistical module of the CIVA platform in order to carry out Probability of Detection (POD) calculations in the case of the ECT of a rivet neighbored by a flaw. Two parametric models used for the POD estimation are presented and simulation results are discussed.