An operator-splitting (point implicit) procedure combined with an efficient ordinary differential equation solver has been incorporated into a pressure based code to study finite rate reacting flows. In the predictor step of this procedure, the chemical kinetics terms are integrated implicitly with an ordianry differential equation solver over a characteristic time scale of fluid. In the corrector step, the resulting effective chemical source terms are substituted into fluid dynamics equations and the integration takes place together with convection/diffusion terms in a non-iterative time-marching algorithm. A penalty function method is also presented for the steady-state reacting flow to improve the computational efficiency. The present method has been validated by several benchmark test cases in one and two dimensional reacting flows. Comparisons of the transient calculation between available experimental data and existing numerical results indicate that the combustion wave propagations have been correctly predicted for such highly complicated problems even though some unknown aspects need to be investigated in future studies.
展开▼