A computationally efficient P1 radiation model for modern combustion systems utilizing pre-conditioned conjugate gradient methods

Utilizing Pre-Conditioners (PC) to the Conjugate Gradients (CG) method, the convergence time of the P1 radiation model reduced by a factor of 30 for modest problem sizes and a factor of 70 for larger sized problems when compared against classical Gauss Seidel sweeps. Further, PC provided 50% computational savings compared to employing CG in a standalone mode. The P1 model calculation times were about 25-30% of the DO model calculation time. The time to solution also scaled linearly with an increase in problem size. The weighted sum of gray gases model employed in this study in conjunction with the P1 model provided good agreement against benchmark data with L2 error norms (defined relative to corresponding DO calculations) improving when isotropic intensities were prevalent.