Modeling thermal radiation in combustion environments can be extremely challenging for two main reasons. First, the radiative transfer equation (RTE), which is the cornerstone of modeling radiation in such environments, is a five-dimensional integro-differential equation. Second, the absorption and scattering coefficients of molecular gases and particulates prevalent in combustion environments oscillate strongly with the wavenumber (or wavelength), i.e., the medium is strongly nongray, requiring the solution of the RTE for a large number of wavenumbers. This article reviews the progress that has been made in this area to date with an emphasis on the work performed over the past three decades. Progress in both deterministic and stochastic (Monte Carlo) solutions of the RTE is reviewed, in addition to the review of the treatment of the spectral properties of gases, soot, and fuel droplets that dominate combustion environments, i.e., spectral or nongray models. The application of the various state-of-the-art nongray models and RTE solution methods to flames (particularly turbulent), fires, combustors, and other combustion systems are summarized along with a critical discussion of the pros and cons of the models and methods. Finally, the challenges that remain in modeling thermal radiation in combustion systems are highlighted and future outlooks are shared.