TY - JOUR
T1 - Magnetohydrodynamic simulation code CANS+
T2 - Assessments and applications
AU - Matsumoto, Yosuke
AU - Asahina, Yuta
AU - Kudoh, Yuki
AU - Kawashima, Tomohisa
AU - Matsumoto, Jin
AU - Takahashi, Hiroyuki R.
AU - Minoshima, Takashi
AU - Zenitani, Seiji
AU - Miyoshi, Takahiro
AU - Matsumoto, Ryoji
N1 - Funding Information:
CANS+ was based on the code originally developed by T. Ogawa at Chiba University. The numerical setup for magnetic reconnection was developed by N. Iwamoto, A. D. Kawamura, J. Sakamoto, and T. Shibayama during the simulation summer school held at Chiba University in 2014. The present simulations used computational resources provided by the Information Technology Center, the University of Tokyo, the Research Institute for Information Technology, Kyushu University, and the RIKEN Center for Computational Science through the HPCI System Research project (Project ID: hp120193, hp120287, hp130027, hp140213, hp140170, hp150263), and Cray XC50 at the Center for Computational Astrophysics, National Astronomical Observatory of Japan. This work was supported in part by MEXT SPIRE, MEXT as “Priority Issue on Post-K computer” (Elucidation of the Fundamental Laws and Evolution of the Universe), JICFuS, the Research Institute of Stellar Explosive Phenomena at Fukuoka University (JM), and JSPS KAK-ENHI Grant Numbers 16H03954 (RM) and 17K14260 (HRT).
Publisher Copyright:
© 2019 The Author(s). Published by Oxford University Press on behalf of the Astronomical Society of Japan.
PY - 2019/8/1
Y1 - 2019/8/1
N2 - We present a new magnetohydrodynamic (MHD) simulation package with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the Harten-Lax-van Leer-discontinuitues (HLLD) approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes has significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems show the advantages of using the high-order scheme by comparing with results from a standard second-order total variation diminishing monotonic upwind scheme for conservation laws (MUSCL) scheme. The present code enables us to explore the long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence causes continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibits a convergent picture of the early nonlinear growth of the MRI in a global model, and indicates that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.
AB - We present a new magnetohydrodynamic (MHD) simulation package with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the Harten-Lax-van Leer-discontinuitues (HLLD) approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes has significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems show the advantages of using the high-order scheme by comparing with results from a standard second-order total variation diminishing monotonic upwind scheme for conservation laws (MUSCL) scheme. The present code enables us to explore the long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence causes continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibits a convergent picture of the early nonlinear growth of the MRI in a global model, and indicates that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.
KW - magnetohydrodynamics (MHD)
KW - methods: numerical
KW - shock waves
KW - turbulence
UR - http://www.scopus.com/inward/record.url?scp=85070825079&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85070825079&partnerID=8YFLogxK
U2 - 10.1093/pasj/psz064
DO - 10.1093/pasj/psz064
M3 - Article
AN - SCOPUS:85070825079
SN - 0004-6264
VL - 71
JO - Publications of the Astronomical Society of Japan
JF - Publications of the Astronomical Society of Japan
IS - 4
M1 - 83
ER -