Parallelization of the P-1 radiation model

The P-1 radiation model is spatially decomposed to solve the radiative transport equation on parallel computers. Mathematical libraries developed by third parties are employed to solve the linear systems that result during the solution procedure. Multigrid preconditioning accelerated the convergence of iterative methods. The parallel performance did not depend strongly on the radiative properties of the medium or the boundary conditions. Predictions from coupling the weighted-sum-of-gray-gases model with the P-1 approximation are compared against benchmarks for model problems. The P-1 approximation resulted in only a moderate loss in accuracy while being significantly faster than the discrete ordinates method.