TransROTA: A Code for Solving Collisionality-Extended Braginskii's Closure Formalism for Toroidally-Rotating Plasmas with Nonlinear Successive Overrelaxation

Cheonho Bae,Yifei Jin,Bo Lyu,Jia Fu,Fudi Wang,Hongming Zhang
DOI: https://doi.org/10.1016/j.cpc.2023.108992
IF: 4.717
2023-01-01
Computer Physics Communications
Abstract:TransROTA (Transport analysis based on ROTAtions) is a code that solves the axisymmetric version of StaceySigmar extended plasma rotation model based on a collisionality-extended Braginskii's closure formalism with realistic D-shaped flux surface geometry in toroidally-confined plasmas [Nucl. Fusion 53 (2013) 043,011]. The Stacey-Sigmar model is unique in rotation and transport calculations by retaining all the terms in the angular momentum balance and by considering the first-order poloidal asymmetries in density, velocity, and electrostatic potential, for a higher accuracy in rotation calculations. TransROTA code evolves out of an earlier two-ion code GTROTA [Comp. Phys. Comm. 184 (2013) 2571-2587] and is developed for more practical studies on rotation and transport of up to five ion species (including deuterium and carbon). The main practical application is focused on intrinsic rotation studies. TransROTA requires a well-calibrated radial velocity profile of carbon impurity as one of the 21 required inputs. Once provided, TransROTA stably converges its nonlinear iterations with locally-optimized Successive Overrelaxation (SOR) scheme only, which is a remarkable improvement from GTROTA. This stable convergence is achieved by both modified application of the model equations and improvement in the numerical method. The new numerical scheme successfully removes the numerical instability issues commonly observed in GTROTA to make TransROTA not only to robustly converge the solutions but also to be faster and user-friendly.
What problem does this paper attempt to address?