AbstractThe conventional thermal network model for satellite surface temperature distribution is directly solved by a new solution method on the basis of the Monte Carlo Ray Tracing (MCRT) method. The solar direct incidence area, the solar radiation transfer coefficient, the infrared radiation transfer coefficient, and the network conduction and radiation coefficients are calculated by Monte Carlo statistical analysis rather than the Gebhart approach. The advantage of the MCRT method is that the surface material characteristics are taken into account in the solving process. The effect of absorptivity and emissivity for temperature distributions is analyzed in detail. Numerical simulation is carried out for the temperature distribution of the satellite surface with different solar incidences of round angle and zenith angle. In the steady stage, the maximum temperature difference reaches up to 200 K and increasing the absorptivity‐emissivity ratio will lead to a rapid temperature rise. In the unsteady stage, the cooling velocity of the main body surface is always smaller than that of the solar array surface and the maximum temperature of the main body is larger than the solar arrays. Under the same absorptivity‐emissivity ratio, the solar incidence angle has a great effect on the uniform character of temperature distribution. © 2010 Wiley Periodicals, Inc. Heat Trans Asian Res, 39(7): 539–553, 2010; Published online 16 July 2010 in Wiley Online Library (wileyonlinelibrary.com). DOI 10.1002/htj.20316