IRUS Total

A methodology for the integration of stiff chemical kinetics on GPUs

File Description SizeFormat 
Combustion and Flame_2014.pdfAccepted version452.67 kBAdobe PDFView/Open
Title: A methodology for the integration of stiff chemical kinetics on GPUs
Authors: Sewerin, F
Rigopoulos, S
Item Type: Journal Article
Abstract: Numerical schemes for reacting flows typically invoke the method of fractional steps in order to isolate the chemical kinetics model from diffusion/convection phenomena. Here, the reaction fractional step requires the solution of a collection of independent ODE systems which may be severely stiff. Recently, researchers have begun to explore the highly parallel structure of graphics processing units (GPUs) in order to accelerate integration schemes for these ODE systems. However, much of the existing work concentrates on explicit integration algorithms which may fall short in the presence of stiffness. In this light, we have carefully reimplemented in OpenCL C the Fortran 77 program of the 3-stage/5th order implicit Runge–Kutta method Radau5 by Hairer and Wanner (1991) and tested it extensively in the context of a transient equilibrium scheme for the flamelet model. Our implementation can easily be integrated with any existing reactive flow software in order to solve the reaction fractional step on an OpenCL-enabled GPU. Moreover, it is suited for any Chemkin-format reaction mechanism with ≲200≲200 species without incurring a loss in occupancy and it reaches its limit speedup (which is largely independent of the mechanism size) at a small problem size (≈500 ODE systems). In view of memory constraints, we include an optimized scheme for splitting the ODE systems across several kernel invocations and overlapping the kernel execution with data transfers. An in-depth evaluation is based upon runtime measurements of the CPU and the GPU implementation on a user level and a high-end CPU/GPU for an increasing number of ODE systems, reduced and detailed reaction mechanisms and a range of time step sizes.
Issue Date: 1-Apr-2015
Date of Acceptance: 3-Nov-2014
URI: http://hdl.handle.net/10044/1/18757
DOI: 10.1016/j.combustflame.2014.11.003
ISSN: 0010-2180
Publisher: Elsevier
Start Page: 1375
End Page: 1394
Journal / Book Title: Combustion and Flame
Volume: 162
Issue: 4
Copyright Statement: © 2014 The Combustion Institute.Copyright. NOTICE: this is the author’s version of a work that was accepted for publication in Combustion and Flame. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Combustion and Flame, 2 December 2014 DOI 10.1016/j.combustflame.2014.11.003
Keywords: Science & Technology
Physical Sciences
Energy & Fuels
Engineering, Multidisciplinary
Engineering, Chemical
Engineering, Mechanical
Chemical kinetics
Implicit Runge-Kutta methods
0902 Automotive Engineering
0904 Chemical Engineering
0913 Mechanical Engineering
Publication Status: Published
Online Publication Date: 2014-12-02
Appears in Collections:Mechanical Engineering
Faculty of Engineering