La présence de petites inclusions dans un domaine de référence modifie la solution de l'équation de Laplace dans ce domaine. Les cas d'une inclusion isolée ou de plusieurs bien séparées ont été largement étudiés. Dans cet exposé, nous considérons le cas où la distance entre deux inclusions tend vers zéro mais moins vite que la taille caractéristique des inclusions. Nous donnons un développement asymptotique multi-échelle complet de la solution de l'équation de Laplace dans la situation de deux inclusions parfaitement isolantes. Nous présentons également quelques simulations numériques basée sur une méthode de superposition multi-échelle provenant du développement au premier ordre. Nous appliquerons ensuite les mêmes techniques pour étudier l'équation de l'élasticité linéaire afin de prédire le comportement à rupture de certains matériaux présentant des micro-défauts. Ce travail a été réalisé en collaboration avec D. Brancherie, M. Dambrine, F. Hérau, S. Tordeux et G. Vial.