Accelerating moderately stiff chemical kinetics in reactive-flow simulations using GPUs

Abstract

The chemical kinetics ODEs arising from operator-split reactive-flow simulations were solved on GPUs using explicit integration algorithms. Nonstiff chemical kinetics of a hydrogen oxidation mechanism (9 species and 38 irreversible reactions) were computed using the explicit fifth-order Runge-Kutta-Cash-Karp method, and the GPU-accelerated version performed faster than single- and six-core CPU versions by factors of 126 and 25, respectively, for 524,288 ODEs. Moderately stiff kinetics, represented with mechanisms for hydrogen/carbon-monoxide (13 species and 54 irreversible reactions) and methane (53 species and 634 irreversible reactions) oxidation, were computed using the stabilized explicit second-order Runge-Kutta-Chebyshev (RKC) algorithm. The GPU-based RKC implementation demonstrated an increase in performance of nearly 59 and 10 times, for problem sizes consisting of 262,144 ODEs and larger, than the single- and six-core CPU-based RKC algorithms using the hydrogen/carbon-monoxide mechanism. With the methane mechanism, RKC-GPU performed more than 65 and 11 times faster, for problem sizes consisting of 131,072 ODEs and larger, than the single- and six-core RKC-CPU versions, and up to 57 times faster than the six-core CPU-based implicit VODE algorithm on 65,536 ODEs. In the presence of more severe stiffness, such as ethylene oxidation (111 species and 1566 irreversible reactions), RKC-GPU performed more than 17 times faster than RKC-CPU on six cores for 32,768 ODEs and larger, and at best 4.5 times faster than VODE on six CPU cores for 65,536 ODEs. With a larger time step size, RKC-GPU performed at best 2.5 times slower than six-core VODE for 8192 ODEs and larger. Therefore, the need for developing new strategies for integrating stiff chemistry on GPUs was discussed.

27 pages, LaTeX; corrected typos in Appendix equations A.10 and A.11

Cited by in corpus (15)

• pyJac: analytical Jacobian generator for chemical kinetics
• PeriPy -- A High Performance OpenCL Peridynamics Package
• A Neural Network Inspired Formulation of Chemical Kinetics
• 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
• Concentration-Flux-Steered Mechanism Exploration with an Organocatalysis Application
• Effects of Langmuir Turbulence on Upper Ocean Carbonate Chemistry
• Using SIMD and SIMT vectorization to evaluate sparse chemical kinetic Jacobian matrices and thermochemical source terms
• GPU-Based Parallel Integration of Large Numbers of Independent ODE Systems
• Accelerating reactive-flow simulations using vectorized chemistry integration
• SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems
• Optimized thread-block arrangement in a GPU implementation of a linear solver for atmospheric chemistry mechanisms
• Solving large number of non-stiff, low-dimensional ordinary differential equation systems on GPUs and CPUs: performance comparisons of MPGOS, ODEINT and DifferentialEquations.jl
• Performance assessment of CUDA and OpenACC in large scale combustion simulations
• Adapting Atmospheric Chemistry Components for Efficient GPU Accelerators