SIESTA is a new ideal MHD 3D-equilibrium solver, capable of dealing with magnetic islands and stochastic regions in a fast, accurate and scalable manner. [1] The SIESTA code will complement other existent 3D MHD island solvers and is particularly suited for applications that require not only accuracy but speed of evaluation as well, such as experimental 3D equilibrium reconstruction or stellarator design. SIESTA will also be useful to calculate MHD equilibria at very high spatial resolutions, such as those that might be required for the investigation of NTMs at ITER-relevant temperature and resistivity conditions. SIESTA is based on a preconditioned, iterative algorithm that takes advantage of a pre-existent VMEC solution to provide both a background coordinate system that guarantees a compact representation as well as a good initial guess for the iterative procedure.

References