A partir de l'établissement d'une approximation in stationnaires en 3D des potentiels retardés pour des courants électromagnétiques sur des polyèdres non-nécessairement convexes, une méthode de résolution pour la simulation de la diffraction par un ensemble non-connexe d'objets est formulée. Une partition de ce dernier est effectuée suivant les inhomogénéités présentes. Le problème est alors traduit en un système d'équations de Maxwell couplées, chacune étant homogène hors d'un élément correspondant de la partition, qui induit la construction d'une solution du problème initial. Par approximation des termes de couplage, il s'en suit une méthode naturellement hybride et parallèle sur un système stable et bien-posé. La restriction de chaque sous-système à un voisinage du support de ses inhomogénéités est obtenue par introduction de conditions aux limites absorbantes de type "PML" dont un formalisme généralisé est étudié. Des exemples numériques illustrent l'ensemble de ces développements.