Accelerating Euler Equations Numerical Solver on Graphics Processing Units

Pierre Kestener, Frederic Chateau, Romain Teyssier
CEA, Centre de Saclay, DSM/IRFU/SEDI, F-91191 Gif-Sur-Yvette, France
Algorithms and Architectures for Parallel Processing, Lecture Notes in Computer Science, 2010, Volume 6082/2010, 281-288


   title={Accelerating Euler Equations Numerical Solver on Graphics Processing Units},

   author={Kestener, P. and Ch{^a}teau, F. and Teyssier, R.},

   journal={Algorithms and Architectures for Parallel Processing},





Download Download (PDF)   View View   Source Source   



Finite volume numerical methods have been widely studied, implemented and parallelized on multiprocessor systems or on clusters. Modern graphics processing units (GPU) provide architectures and new programing models that enable to harness their large processing power and to design computational fluid dynamics simulations at both high performance and low cost. We report on solving the 2D compressible Euler equations on modern Graphics Processing Units (GPU) with high-resolution methods, i.e. able to handle complex situations involving shocks and discontinuities. We implement two different second order numerical schemes, a Godunov-based scheme with quasi-exact Riemann solver and a fully discrete second-order central scheme as originally proposed by Kurganov and Tadmor. Performance measurements show that these two numerical schemes can achieves x30 to x70 speed-up on recent GPU hardware compared to a mono-thread CPU reference implementation. These first results provide very promising perpectives for designing a GPU-based software framework for applications in computational astrophysics by further integrating MHD codes and N-body simulations.
No votes yet.
Please wait...

* * *

* * *

HGPU group © 2010-2021 hgpu.org

All rights belong to the respective authors

Contact us: