MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

TL;DR

MAS is a new versatile eigenvalue code for plasma stability analysis in general tokamak geometry, incorporating kinetic effects within a Landau-fluid framework, and validated through theoretical benchmarks and application to DIII-D data.

Contribution

The paper introduces MAS, a comprehensive Landau-fluid eigenvalue code that efficiently models plasma stability with kinetic effects in general geometry, bridging fluid and kinetic theories.

Findings

MAS accurately reproduces known plasma wave behaviors.

The code successfully models Alfven eigenmodes in experimental data.

Benchmarking confirms the physical and numerical reliability of MAS.

Abstract

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number n and frequency omega in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field, ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfven wave, ion sound wave, low-n kink, high-n ion temperature gradient mode and kinetic ballooning mode. Moreover, MAS is successfully applied to model the Alfven eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of RSAE, the tunneling damping of TAE, as well as the polarization characteristics of KBAE and BAAE being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.