Techniques, Tricks and Algorithms for Efficient GPU-Based Processing of Higher Order Hyperbolic PDEs

DOI10.1007/S42967-022-00235-9MaRDI QIDQ6418316zbMATH OpenFDO

Authors S. Subramanian, Dinshaw S. Balsara, Deepak Bhoriya, Harish Kumar

Publication date 23 November 2022

PDEs in connection with fluid mechanics (35Q35) Maxwell equations (35Q61) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08) Magnetohydrodynamics and electrohydrodynamics (76W05) Models and methods for concurrent and distributed computing (process algebras, bisimulation, transition nets, etc.) (68Q85) Finite volume methods applied to problems in fluid mechanics (76M12) PDEs in connection with computer science (35Q68) PDEs in connection with astronomy and astrophysics (35Q85) Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs (65M22)

Abstract: GPU computing is expected to play an integral part in all modern Exascale supercomputers. It is also expected that higher order Godunov schemes will make up about a significant fraction of the application mix on such supercomputers. It is, therefore, very important to prepare the community of users of higher order schemes for hyperbolic PDEs for this emerging opportunity. We focus on three broad and high-impact areas where higher order Godunov schemes are used. The first area is computational fluid dynamics (CFD). The second is computational magnetohydrodynamics (MHD) which has an involution constraint that has to be mimetically preserved. The third is computational electrodynamics (CED) which has involution constraints and also extremely stiff source terms. Together, these three diverse uses of higher order Godunov methodology, cover many of the most important applications areas. In all three cases, we show that the optimal use of algorithms, techniques and tricks, along with the use of OpenACC, yields superlative speedups on GPUs! As a bonus, we find a most remarkable and desirable result: some higher order schemes, with their larger operations count per zone, show better speedup than lower order schemes on GPUs. In other words, the GPU is an optimal stratagem for overcoming the higher computational complexities of higher order schemes! Several avenues for future improvement have also been identified. A scalability study is presented for a real-world application using GPUs and comparable numbers of high-end multicore CPUs. It is found that GPUs offer a substantial performance benefit over comparable number of CPUs, especially when all the methods designed in this paper are used.

This page was built for publication: Techniques, Tricks and Algorithms for Efficient GPU-Based Processing of Higher Order Hyperbolic PDEs

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6418316)