A novel method for implementing the statistical narrowband model into the radiative transfer equation was devised to perform nongrey radiative transfer calculationsin three-dimensional absorbing,emitting, and scattering media. In this new method the radiation intensity is split into two parts: the nonscattered part and the scattered part. The nonscattered part is solved accurately using the statistical narrowband model and a ray-tracing technique and does not require iteration. The scattered part is solved with approximation by using the gray-band model to estimate the band correlation between the scattered intensity and the gas absorption coef cient. The accuracy of this method for narrowband radiation intensity prediction was evaluated in a homogeneous H 2 O-N 2 -Al 2 O 3 mixture under both isothermal and nonisothermal conditions by comparing its results with those of the statistical narrowband correlated-K method. This new implementation method alleviates to some extent the dif culty of the band model-scattering incompatibility and offers good results provided the particle loading is not too high. Nomenclature f = species molar fraction I nm = nonscattered part of spectral radiation intensity, W m ¡2 cm sr ¡1 I sm = scattered part of spectral radiation intensity, W m ¡2 cm sr ¡1 I m = spectral radiation intensity, W m ¡2 cm sr ¡1 N I m = mean radiation intensity over a narrowband, W m ¡2 cm sr ¡1 k jm = absorption coef cient associated with the jth quadrature point, m ¡1 N k m = mean line intensity to spacing ratio, cm ¡1 atm ¡1 L = path length, m l m = local mean path length, m n = particle number density, m ¡3 n = unit normal vector of a wall surface pointing toward the gas side p = pressure, atm S = surface area of a computational grid, m 2 s, s 0 = position variables along a line of sight, m T = temperature, K V = volume of a computational grid, m 3 w jm = weight parameter associated with the j th quadrature point x, y, z = Cartesian coordinates, m N b m = mean line width to spacing ratio N c m = mean half width of an absorption line, cm ¡1 D m = wave-number interval of a narrowband, cm ¡1 N d m = equivalent line spacing, cm ¡1 ² w = wall surface emissivity j gm = gas spectral absorption coef cient, m ¡1 N j gm = narrowband averaged gas absorption coef cient, m ¡1 pm = particulate spectral absorption coef cient, m ¡1 j sm = particulate spectral scattering coef cient, m ¡1 m = wave number, cm ¡1 n , g , l = direction cosines t gm = gas spectral transmissivity t nm = medium spectral transmissivity associated with I nm U = scattering phase function \ = solid angle, sr X = direction of radiation propagation Subscripts b = blackbody g = gas i = spatial discretization (along a line of sight) index n = angular discretization index p = particulate w = wall m = spectral