Integration Strategies for Efficient Multizone Chemical Kinetics Models
Three integration strategies are developed and tested for the stiff, ordinary differential equation (ODE) integrators used to solve the fully coupled multizone chemical kinetics model. Two of the strategies tested are found to provide more than an order of magnitude of improvement over the original, basic level of usage for the stiff ODE solver. One of the faster strategies uses a decoupled, or segregated, multizone model to generate an approximate Jacobian. This approach yields a 35-fold reduction in the computational cost for a 20 zone model. Using the same approximate Jacobian as a preconditioner for an iterative Krylov-type linear system solver, the second improved strategy achieves a 75-fold reduction in the computational cost for a 20 zone model. The faster strategies achieve their cost savings with no significant loss of accuracy. The pressure, temperature and major species mass fractions agree with the solution from the original integration approach to within six significant digits; and the radical mass fractions agree with the original solution to within four significant digits. The faster strategies effectively change the cost scaling of the multizone model from cubic to quadratic, with respect to the number of zones. As a consequence of the improved scaling, the 40 zone model offers more than a 250-fold cost savings over the basic calculation.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 977222
- Report Number(s):
- LLNL-CONF-418126; TRN: US201013%%437
- Resource Relation:
- Journal Volume: 3; Journal Issue: 1; Conference: Presented at: SAE World Congress 2010, Detroit, MI, United States, Apr 13 - Apr 15, 2010
- Country of Publication:
- United States
- Language:
- English
Similar Records
KRYSI, an ODE solver combining a semi-implicit Runge-Kutta method and a preconditioned Krylov method
Reduced-storage techniques in the numerical method of lines