Massively parallel simulations of relativistic fluid dynamics on graphics processing units with CUDA

Relativistic fluid dynamics is a major component in dynamical simulations of the quark-gluon plasma created in relativistic heavy-ion collisions. Simulations of the full three-dimensional dissipative dynamics of the quark-gluon plasma with fluctuating initial conditions are computationally expensive and typically require some degree of parallelization. In this paper, we present a GPU implementation of the Kurganov-Tadmor algorithm which solves the 3 + 1d relativistic viscous hydrodynamics equations including the effects of both bulk and shear viscosities. We demonstrate that the resulting CUDA-based GPU code is approximately two orders of magnitude faster than the corresponding serial implementation of the Kurganov-Tadmor algorithm. We validate the code using (semi-)analytic tests such as the relativistic shock-tube and Gubser flow. Program summary Program Title: GPU-VH Program Files doi: http://dx.doi.org/10.17632/dhnmtfpz9k.1 Licensing Provisions: GPLv3 Programming language: CUDA C External routines/libraries: Google Test, GNU Scientific Library (GSL) Nature of problem: Dynamical evolution of the fluid dynamic stage of the quark-gluon plasma produced in nuclear collisions Solution method: Kurganov-Tadmor algorithm (C) 2017 Elsevier B.V. All rights reserved.