A general scheme is developed to solve many-particle problems of therm
al conduction in two dimensions. Based on the Green's function formali
sm, heat sources induced by discontinuity of temperature fields and te
mperature potentials on inclusion boundaries are explicitly taken into
account, and analytical expressions for temperature distributions in
granular systems. ale established. The unknown coefficients in the ana
lytical expressions for temperature field are determined by a matrix e
quation whose elements are all monomials of the distances between part
icle centers. Decomposition of the matrix equation for the systems con
taining a chain of particles is discussed. Illustrative calculations a
re presented for granular systems with two and three particles.