A New High-Order Magnetohydrodynamic Code with a High-Order Constrained Transport Scheme

Abstract

A New High-Order Magnetohydrodynamic Code with a High-Order Constrained Transport Scheme

Dongsu Ryu^1, Jeongbhin Seo^2

1Depart. of Physics, College of Natural Sciences, UNIST, Rep. of Korea 2Los Alamos National Laboratory, Los Alamos, USA

We introduce a magnetohydrodynamic (MHD) code tailored for astrophysical applications. The code leverages a fifth-order finite-difference weighted essentially non-oscillatory (WENO) scheme for spatial discretization, coupled with a fourth-order strong stability-preserving Runge-Kutta (SSPRK) scheme for temporal integration. In particular, to match the high-order accuracy of the WENO scheme, we implement a novel constrained transport (CT) algorithm of high-order accuracy to enforce div(B) = 0. Here, we outline the MHD code and present its performance. We validate that with the new CT scheme, the code achieves the order of accuracy of the WENO scheme in wave decay tests, including a circularly polarized Alfven wave propagation test. Further, the code's high robustness and fidelity are demonstrated through tests involving complex flows.