Clustered frequency analysis of shear Alfv́n modes in stellarators

The shear Alfv́n spectrum in three-dimensional configurations, such as stellarators and rippled tokamaks, is more densely populated due to the larger number of mode couplings caused by the variation in the magnetic field in the toroidal dimension. This implies more significant computational requirements that can rapidly become prohibitive as more resolution is requested. Alfv́n eigenfrequencies and mode structures are a primary point of contact between theory and experiment. A new algorithm based on the Jacobi-Davidson method is developed here and applied for a reduced magnetohydrodynamics model to several stellarator configurations. This technique focuses on finding a subset of eigenmodes clustered about a specified input frequency. This approach can be especially useful in modeling experimental observations, where the mode frequency can generally be measured with good accuracy and several different simultaneous frequency lines may be of interest. For cases considered in this paper, it can be a factor of 10² - 10³ times faster than more conventional methods.