A Virtual Simons' Hour Talk: Models and Methods for Nonlinear Magnetohydrodynamic Simulations of Stellarators

Given by: Nikita Nikulsin
Affiliation: Max Planck Institute for Plasma Physics

Note the time: 8:00am EST / 13:00 UTC.
Zoom link will be sent via email.

Abstract: The JOREK code has recently been extended to allow nonlinear fully 3D stellarator simulations. This is made possible by generalizing the JOREK reduced MHD model to support stellarator geometries, and by allowing the grid to be non-axisymmetric, so that it can be aligned to the flux surfaces in a stellarator. The stellarator reduced model differs mainly in that the magnetic field can be represented as any curl-free field plus a perturbation in the stellarator model, whereas in the tokamak model it is a toroidal field plus a perturbation. It is shown that this model conserves energy, but introduces an error into momentum conservation. An alternate model, which does not guarantee energy conservation, but has a smaller momentum conservation error is also derived. The energy and momentum conservation properties of the main and alternate models are then studied numerically in the tokamak limit. The main model was then tested on a set of l=2 stellarator equilibria based on Wendelstein 7-A. The simulations demonstrate that stable full MHD equilibria are preserved in the reduced model: the flux surfaces do not move throughout the simulation, and closely match the full MHD flux surfaces. Further, both tearing and ballooning modes were simulated, and their growth rates benchmarked against the linear full MHD code CASTOR3D, showing good agreement.