Numerical simulation for the MHD system in 2D using OpenCL

Abstract

Abstract. In this work we compute the MHD equations with divergence cleaning on GPU. The method is based on the finite volume approach and Strang dimensional splitting. The simplicity of the approach makes it a good candidate for a GPU implementation with OpenCL. With adequate memory optimization access, we achieve very high speedups, compared to a classical sequential implementation.Résumé. Dans ce travail, nous résolvons les équations de la MHD avec correction de divergence sur carte graphique. La méthode est basée sur les volumes finis et le splitting directionnel de Strang. La simplicité de l'algorithme en fait un bon candidat pour la programmation sur carte graphique sous OpenCL. Avec de bonnes optimisations des accès mémoire, nous obtenons de très bonnes accélérations, comparé à une programmation séquentielle classique.

Asynchronous OpenCL/MPI Numerical Simulations of Conservation Laws

“…The numerical simulations thus require heavy computations and are well adapted to GPU hardware. For more details and bibliography on the MHD equations, we refer to [11].…”

“…It is not the objective of this work to detail the MHD equations. For this we refer for instance to [11]. In Section 2, we present numerical results obtained with the MHD equations.…”