A 3D PISO-MHD was implemented on Kepler-class graphics processing units (GPUs) using CUDA technology.
A consistent and conservative scheme is used in the code which was validated by three basic benchmarks in a rectangular and round ducts.
Parallelized of CPU and GPU acceleration were compared relating to single core CPU in MHD problems and non-MHD problems.
Different preconditions for solving MHD solver were compared and the results showed that AMG method is better for calculations.