The development of mathematical and optimization models for reverse supply network design has concerned considerable interest over the past decades. However, the uncertainties that are inherent in the network design and the complex interactions among various uncertain parameters are challenging the capabilities of these developed tools. The aim of this paper is to propose a new mathematical programming model for recycling network design in the iron and steel industry. The considered recycling network is multi-echelon, multi-facility, multi-product, and multi-supplier. Moreover, another objective of this research is to introduce an interval-stochastic robust optimization methodology to deal with various uncertainties in the proposed model. Computational experiments are provided to demonstrate the applicability of the proposed model in recycling network design.