Due to their distinct textures, porous cellular materials have been of interest to many engineering applications. Energy conversion through such materials, especially at high temperatures, is governed by naturally coupled conduction-radiation physics. Therefore, a thorough understanding of the conduction-radiation behavior within these materials is important to properly design and optimize these materials. Several numerical methods developed over the last few decades allow to solve and study the influence of different textural parameters on coupled conductive-radiative heat transfers. These numerical methods require 3D digitized images of the sample of interest obtained either using 3D imagery technique or generated using numerical algorithms. To better represent the complex geometry, the 3D image of the sample requires high spatial resolution, and for a sample which is representative of involved physics might contain hundreds of millions of voxels which is complex to solve and computationally expensive.To cope with this issue, we developed a parallelized discrete scale numerical approach using cell centered Finite Volume Method (FVM) and deterministic ray tracing to solve coupled heat transfer within highly porous large scale complex cellular materials. At the heart of this method rests a decomposition approach based on modified zonal method which significantly reduces the coupling efforts, memory requirements, and computation time. The results of two test cases presented in this study are cross-verified with those presented in literature and computed using Star-CCM+ software. Finally, the method is applied to virtual fibrous sample. These results present the ability of the solver to consider different boundary conditions, large temperature gradient across boundaries, and dealing with large samples consisting hundreds of million of voxels while providing accurate results. We also present and discuss the sensitivity of other associated parameters on the final results.