Abstract
A first-passage simulation scheme for the computation of effective conductivities of composites with interfacial resistance at the matrix-inclusion interface was developed. The scheme was illustrated with the computation of effective conductivities of two-phase composites containing regularly or randomly distributed spheres. The mean hitting probabilities and mean scaled traveling times of the probing walkers in the close vicinity of the imperfect contact interface were derived by solving boundary value problems. It was observed that the contact between the matrix and inclusions was imperfect such that an interfacial resistance exists at the matrix-inclusion interface, giving rise to a discontinuity in the temperature fields.