This article reports results of the four popular and widely used numerical methods, viz., the Monte Carlo method (MCM), the discrete transfer method (DTM), the discrete ordinates method (DOM) and the finite volume method (FVM) used to calculate radiative information in any thermal problem. Different classes of problems dealing with radiation and/or conduction heat transfer problems in a 2-D rectangular absorbing, emitting and scattering participating medium have been considered. In radiative equilibrium and non-radiative equilibrium cases, the MCM results have been used as the benchmark data for comparing the performances of the DIM, the DOM and the FVM. In the combined radiation and conduction mode problem, the energy equation has been formulated using the lattice Boltzmann method (LBM). To compare the performance of the DTM, the DOM and the FVM, the required radiative field data computed using these methods have been provided to the LBM formulation. Temperature distributions obtained using the four methods and those obtained from the LBM in conjunction with the DIM, the DOM and the RIM have been compared for different parameters such as the extinction coefficient, the scattering albedo, the conduction-radiation parameter, the wall emissivity, the aspect ratio and heat generation rate. In all the cases, results of these methods have been found in good agreements. Computationally, the DTM was found the most time consuming, and the DOM was computationally the most efficient. (C) 2012 Elsevier Ltd. All rights reserved.