An investigation of GPU-based stiff chemical kinetics integration methods

Abstract

A fifth-order implicit Runge-Kutta method and two fourth-order exponential integration methods equipped with Krylov subspace approximations were implemented for the GPU and paired with the analytical chemical kinetic Jacobian software pyJac. The performance of each algorithm was evaluated by integrating thermochemical state data sampled from stochastic partially stirred reactor simulations and compared with the commonly used CPU-based implicit integrator CVODE. We estimated that the implicit Runge-Kutta method running on a single GPU is equivalent to CVODE running on 12-38 CPU cores for integration of a single global integration time step of 1e-6 s with hydrogen and methane models. In the stiffest case studied---the methane model with a global integration time step of 1e-4 s---thread divergence and higher memory traffic significantly decreased GPU performance to the equivalent of CVODE running on approximately three CPU cores. The exponential integration algorithms performed more slowly than the implicit integrators on both the CPU and GPU. Thread divergence and memory traffic were identified as the main limiters of GPU integrator performance, and techniques to mitigate these issues were discussed. Use of a finite-difference Jacobian on the GPU---in place of the analytical Jacobian provided by pyJac---greatly decreased integrator performance due to thread divergence, resulting in maximum slowdowns of 7.11-240.96 times; in comparison, the corresponding slowdowns on the CPU were just 1.39-2.61 times, underscoring the importance of use of an analytical Jacobian for efficient GPU integration. Finally, future research directions for working towards enabling realistic chemistry in reactive-flow simulations via GPU\slash SIMD accelerated stiff chemical kinetic integration were identified.

arXiv comment34 pages, 6 figures; pdfLaTeX

References in corpus (6)

• Recent progress and challenges in exploiting graphics processors in computational fluid dynamics
• pyJac: analytical Jacobian generator for chemical kinetics
• Accelerating moderately stiff chemical kinetics in reactive-flow simulations using GPUs
• An investigation of GPU-based stiff chemical kinetics integration methods
• Accelerating finite-rate chemical kinetics with coprocessors: comparing vectorization methods on GPUs, MICs, and CPUs
• How to obtain efficient GPU kernels: an illustration using FMM & FGT algorithms

Cited by in corpus (7)

• Reducing the complexity of chemical networks via interpretable autoencoders
• An investigation of GPU-based stiff chemical kinetics integration methods
• Accelerating finite-rate chemical kinetics with coprocessors: comparing vectorization methods on GPUs, MICs, and CPUs
• Using SIMD and SIMT vectorization to evaluate sparse chemical kinetic Jacobian matrices and thermochemical source terms
• Accelerating reactive-flow simulations using vectorized chemistry integration
• The Stabilized Explicit Variable-Load Solver with Machine Learning Acceleration for the Rapid Solution of Stiff Chemical Kinetics
• Modular, general purpose ODE integration package to solve large number of independent ODE systems on GPUs