Home >

news Help

Publication Information


Title
Japanese:GPU-Accelerated General Relativistic Radiation MHD Simulations with M-1 Closure 
English:GPU-Accelerated General Relativistic Radiation MHD Simulations with M-1 Closure 
Author
Japanese: 小林諒平, Hiroyuki R. Takahashi, 額田彰, Yuta Asahina, 朴泰祐, Ken Ohsuga.  
English: Ryohei Kobayashi, Hiroyuki R. Takahashi, Akira Nukada, Yuta Asahina, 朴泰祐, Ken Ohsuga.  
Language English 
Journal/Book name
Japanese:Journal of Information Processing 
English:Journal of Information Processing 
Volume, Number, Page Volume 34        Page 149-160
Published date Feb. 15, 2026 
Publisher
Japanese:情報処理学会 
English:Information Processing Society of Japan 
Conference name
Japanese: 
English: 
Conference site
Japanese: 
English: 
Official URL http://dx.doi.org/10.2197/ipsjjip.34.149
 
DOI https://doi.org/10.2197/ipsjjip.34.149
Abstract Black hole accretion disks serve as the central engines driving high-energy phenomena such as X-ray binaries, active galactic nuclei, and gamma-ray bursts, and give rise to complex physical processes such as turbulence induced by magnetorotational instability, intense heating, and the formation of relativistic jets. Traditional one-dimensional theoretical models and general relativistic magnetohydrodynamic (GR-MHD) simulations have successfully captured many aspects of these phenomena, but they fall short in reproducing the effects of radiative cooling and radiation pressure, which are critical in high accretion rate environments. To overcome these limitations, general relativistic radiation magnetohydrodynamics (GR-RMHD) simulations have been developed to incorporate the interplay between radiation and plasma dynamics. However, the inclusion of radiative transfer significantly increases the computational complexity, requiring advanced computational strategies for high-fidelity multidimensional simulations. In this study, we present a GR-RMHD code implemented using CUDA Fortran and MPI that achieves up to 39-fold speedup over conventional CPU implementations, with nearly linear scalability as the number of GPUs increases. This acceleration enables detailed analysis of the complex dynamics of black hole accretion disks and jet formation, providing new insights into the energy production mechanisms in high-energy astrophysical systems. We discuss the evolution of theoretical models and numerical simulations from GR-MHD to GR-RMHD, describe the algorithms and GPU implementation of our approach, and provide quantitative performance evaluations.

©2007 Institute of Science Tokyo All rights reserved.