A Newton-Krylov Based Solver for Modeling Finite Rate Chemistry

To date, computational fluid dynamics (CFD) codes aimed at solving practical engineering problems involving chemically reacting flow have incorporated relatively simple descriptions of the chemical mechanisms involved. Techniques are now available to create reduced mechanisms that faithfully represent detailed chemical descriptions over an appropriate range of conditions using many fewer species. However, including reduced mechanisms into a CFD analysis typically leads to numerical difficulties.

As part of National Science Foundation SBIR Phase II contract no. DMI-0216590, a new modeling tool has been created that utilizes a combination of state-of-the-art techniques used by Reaction Engineering International (REI) for modeling finite rate chemistry in chemically reacting flows using reduced mechanisms with emerging Newton-Krylov methods for solving systems of non-linear equations. For tests problems ranging from geometrically simple combustion problems to full-scale utility boiler simulations, the Newton-Krylov solver has reduced the CPU time to achieve a solution by up to 60% compared to our traditional Picard iteration method.