A non-iterative pressure based scheme for the computation of reacting radiating flows

A non-iterative pressure based algorithm which consists of splitting the solution of momentum energy and species equations into a sequence of predictor–corrector stages was developed for the simulation of transient reacting radiating flows. A semi-discrete approach called the Method of Lines (MOL) which enables implicit time-integration at all splitting stages was used for the solution of conservation equations. The solution of elliptic pressure equation for the determination of the pressure field was performed by a multi-grid (MUDPACK package) solver. Radiation calculations were carried out by coupling an existing radiation code to the algorithm. A first order Arrhenius type rate law expression was utilized to account for the chemistry. The predictions of the algorithm were benchmarked against experimental and numerical data available in the literature. Overall comparisons reveal that numerical results obtained with and without radiation mimic the experimental trends closely. As expected, incorporation of radiation in the simulations leads to better agreement between the predicted and measured velocity and temperature fields when compared to that obtained without radiation. The algorithm developed is an accurate and efficient tool for the simulation of reacting radiating flows and its extension to turbulent flows with the improvement of the existing models is highly promising.