Experiments in numerical methods for a problem in combustion modeling (Q760789)

We show that the numerical method of lines can be used effectively to solve time-dependent combustion models in one spatial dimension. By the numerical method of lines (NMOL), we mean the reduction of a system of partial differential equations to a system of ordinary differential equations (ODE's), followed by the solution of this ODE system with an appropriate ODE solver. We use finite differences for the spatial discretization and a variant of the GEAR package for the ODE's. We present various solution methods of interest for the nonlinear algebraic system in this setting; that is, in the corrector iteration section of the GEAR package applied to combustion models. These methods include Newton/block SOR (SOR denotes successive over-relaxation), block SOR/Newton, Newton/block-diagonal Jacobian, Newton/kinetics-only Jacobian, and Newton/block symmetric SOR. These methods have in common their lack of frequent use in ODE software and their easy applicability to partial differential equations in more than one spatial dimension. Finally, we give the results of numerical tests, run of the CDC-7600 and Cray-1 computers. By doing so, we indicate the more promising nonlinear system solvers for the NMOL solution of combustion models.