Explicit integration with GPU acceleration for large kinetic networks

We demonstrate the first implementation of recently-developed fast explicit kinetic integration algorithms on modern graphics processing unit (GPU) accelerators. Taking as a generic test case a Type Ia supernova explosion with an extremely stiff thermonuclear network having 150 isotopic species and 1604 reactions coupled to hydrodynamics using operator splitting, we demonstrate the capability to solve of order 100 realistic kinetic networks in parallel in the same time that standard implicit methods can solve a single such network on a CPU. This orders-of-magnitude decrease in compute time for solving systems of realistic kinetic networks implies that important coupled, multiphysics problems in various scientific and technical fields that were intractible, or could be simulated only with highly schematic kinetic networks, are now computationally feasible.

[1]  Wright-Patterson Afb,et al.  Techniques for solving stiff chemical kinetics on GPUs , 2013 .

[2]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[3]  W. R. Hix,et al.  Computational methods for nucleosynthesis and nuclear energy generation , 1999 .

[4]  Mike Guidry,et al.  Algebraic stabilization of explicit numerical integration for extremely stiff reaction networks , 2011, J. Comput. Phys..

[5]  W. R. Hix,et al.  Thermonuclear kinetics in astrophysics , 2006 .

[6]  Fabian Sewerin,et al.  A methodology for the integration of stiff chemical kinetics on GPUs , 2015 .

[7]  Hans Werner Meuer,et al.  Top500 Supercomputer Sites , 1997 .

[8]  Kyle E. Niemeyer,et al.  Accelerating moderately stiff chemical kinetics in reactive-flow simulations using GPUs , 2013, J. Comput. Phys..

[9]  Jean-Luc Cambier,et al.  GPU-based flow simulation with detailed chemical kinetics , 2013, Comput. Phys. Commun..

[10]  D. Mott,et al.  New Quasi-Steady-State And Partial-Equilibrium Methods For Integrating Chemically Reacting Systems , 1999 .

[11]  W. Green,et al.  Redesigning combustion modeling algorithms for the Graphics Processing Unit (GPU): Chemical kinetic rate evaluation and ordinary differential equation integration , 2011 .

[12]  Jack J. Dongarra,et al.  LU Factorization of Small Matrices: Accelerating Batched DGETRF on the GPU , 2014, 2014 IEEE Intl Conf on High Performance Computing and Communications, 2014 IEEE 6th Intl Symp on Cyberspace Safety and Security, 2014 IEEE 11th Intl Conf on Embedded Software and Syst (HPCC,CSS,ICESS).

[13]  Oluwayemisi O. Oluwole,et al.  Accelerating multi-dimensional combustion simulations using GPU and hybrid explicit/implicit ODE integration , 2012 .