Abstract: The central-weighted essentially non-oscillatory (CWENO) method is commonly employed for magnetohydrodynamics (MHD) simulations due to its robust numerical stability and high resolution. However, MHD simulations typically rely on Cartesian coordinate systems, whereas laboratory plasmas are often based on curved coordinate systems. To address this, we propose an enhanced version of the CWENO method, which demonstrates third-order accuracy and is specifically adapted to 3D curved coordinate systems. We utilize this method to conduct MHD simulations in both 3D Cartesian and spherical coordinate systems. While the benchmark results of our approach align with those obtained using Zhou’s extension of the CWENO method and the Lax–Friedrichs splitting scheme, our method exhibits more stability than Zhou’s method. Thus, this study offers a more efficient approach for experimental plasma-physics simulations.