swMATH ID: 1030
Software Authors: H. Lütjens; J. Luciani
Description: XTOR-2F solves a set of extended magnetohydrodynamic (MHD) equations in toroidal tokamak geometry. In the original XTOR code, the time stepping is handled by a semi-implicit method. Moderate changes were necessary to transform it into a fully implicit one using the NITSOL library with Newton-Krylov methods of solution for nonlinear system of equations. After addressing the sensitive issue of preconditioning and time step tuning, the performances of the semi-implicit and the implicit methods are compared for the nonlinear simulation of an internal kink mode test case within the framework of resistive MHD including anisotropic thermal transport. A convergence study comparing the semi-implicit and the implicit schemes is presented. Our main conclusion is that on one hand the Newton-Krylov implicit method, when applied to basic one fluid MHD is more computationally costly than the semi-implicit one by a factor 3 for a given numerical accuracy. But on the other hand, the implicit method allows to address challenging issues beyond MHD. By testing the Newton-Krylov method with diamagnetic modifications on the dynamics of the internal kink, some numerical issues, to be addressed further, are emphasized.
Homepage: http://dl.acm.org/citation.cfm?id=1850839.1850885
Keywords: plasma; tokamak; magnetohydrodynamics; nonlinear simulations; fully implicit; Newton-Krylov; internal kink; thermal transport; bifluid; diamagnetic rotations
Related Software: XTOR; BOUT++; CHEASE; NITSOL; Boost C++ Libraries; GYRO; cppmhd; phipm; NIMROD; FLASH; Boost; EPIC; Algorithm 919; CVODE; SUNDIALS; Expokit; MKL; ScaLAPACK; JOREK; FINESSE
Cited in: 8 Publications

Citations by Year