We estimate the heat conducted by a cluster of many small cavities. We show that the dominating heat is a sum, over the number of cavities, of the heat generated by each cavity after interacting with each other. This interaction is described through densities computable as solutions of a closed, and invertible, system of time domain integral equations of second kind. As an application of these expansions, we derive the effective heat conductivity which generates approximately the same heat as the cluster of cavities, distributed in a three-dimensional bounded domain, with explicit error estimates in terms of that cluster. At the analysis level, we use time domain integral equations. Doing that, we have two choices. First, we can favor the space variable by reducing the heat potentials to the ones related to the Laplace operator (avoiding Laplace transform). Second, we can favor the time variable by reducing the representation to the Abel integral operator. As the model under investigation has time-independent parameters, we follow the first approach here.