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

Abstract

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 summaryProgram Title: GPU-VHProgram Files doi:http://dx.doi.org/10.17632/dhnmtfpz9k.1Licensing Provisions: GPLv3Programming language: CUDA CExternal 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 collisionsSolution method: Kurganov-Tadmor algorithm