Relativistic hydrodynamics on graphics processing units

14th International Conference on Strangeness in Quark Matter

Sikorski, Jan^1 ; Cygert, Sebastian^2 ; Porter-Sobieraj, Joanna^2 ; Slodkowski, Marcin^1 ; Krzyanowski, Piotr^1 ; Ksiek, Natalia^1 ; Duda, Przemyslaw^1

Faculty of Physics, Warsaw University of Technology, Koszykowa 75, 00-662 Warsaw, Poland^1

Faculty of Mathematics and Information Science, Warsaw University of Technology, Koszykowa 75, 00-662 Warsaw, Poland^2

关键词: Graphics Processing Unit (GPU);    Hydrodynamics calculations;    Hydrodynamics simulations;    Order of accuracy;    Orders of magnitude;    Quark-gluon plasma;    Relativistic hydrodynamics;    Runge-kutta schemes;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/509/1/012059/pdf
DOI  :  10.1088/1742-6596/509/1/012059

来源: IOP

【 摘 要 】

Hydrodynamics calculations have been successfully used in studies of the bulk properties of the Quark-Gluon Plasma, particularly of elliptic flow and shear viscosity. However, there are areas (for instance event-by-event simulations for flow fluctuations and higher-order flow harmonics studies) where further advancement is hampered by lack of efficient and precise 3+1D program. This problem can be solved by using Graphics Processing Unit (GPU) computing, which offers unprecedented increase of the computing power compared to standard CPU simulations. In this work, we present an implementation of 3+1D ideal hydrodynamics simulations on the Graphics Processing Unit using Nvidia CUDA framework. MUSTA-FORCE (MUlti STAge, First ORder CEntral, with a slope limiter and MUSCL reconstruction) and WENO (Weighted Essentially Non-Oscillating) schemes are employed in the simulations, delivering second (MUSTA-FORCE), fifth and seventh (WENO) order of accuracy. Third order Runge-Kutta scheme was used for integration in the time domain. Our implementation improves the performance by about 2 orders of magnitude compared to a single threaded program. The algorithm tests of 1+1D shock tube and 3+1D simulations with ellipsoidal and Hubble-like expansion are presented.