Abstract Submitted for the DPP20 Meeting of The American Physical Society Extending M3D-C1 to stellarator geometry: preliminary results1

Submitted for the DPP20 Meeting of The American Physical Society Extending M3D-C1 to stellarator geometry: preliminary results1 YAO ZHOU, N. M. FERRARO, S. C. JARDIN, Princeton Plasma Physics Laboratory, H. R. STRAUSS, HRS Fusion — Stellarator plasmas have been observed to be nonlinearly stable even when driven beyond linear MHD stability thresholds. Hence, stellarator designs could employ nonlinear stability considerations to relax linear stability constraints, which can often be too restrictive and costly. However, this possibility has not been systematically investigated due to the lack of a state-ofthe-art nonlinear initial-value MHD code for stellarators. We aim to fill this gap by extending the M3D-C1 code from tokamak to stellarator geometry. Our approach introduces a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework. The mapping between the logical and the physical (R,φ, Z) coordinates is then used to calculate derivatives in the latter, in terms of which the existing physics equations are written. This way, no significant changes to the extended-MHD models within M3D-C1 are required. So far, we have successfully implemented this approach in 2D and verified its results against the existing code. Preliminary results from 3D implementation will also be presented and discussed. 1Supported by the US DOE under Contract No. DEAC02-09CH11466. Yao Zhou Princeton Plasma Physics Laboratory Date submitted: 29 Jun 2020 Electronic form version 1.4