A high performance spectral code for nonlinear MHD stability
A new spectral code, NSTAB, has been developed to do nonlinear stability and equilibrium calculations for the magnetohydrodynamic (MHD) equations in three dimensional toroidal geometries. The code has the resolution to test nonlinear stability by calculating bifurcated equilibria directly. These equilibria consist of weak solutions with current sheets near rational surfaces and other less localized modes. Bifurcated equilibria with a pronounced current sheet where the rotational transform crosses unity are calculated for the International Thermonuclear Experimental Reactor (ITER). Bifurcated solutions with broader resonances are found for the LHD stellarator currently being built in Japan and an optimized configuration like the Wendelstein VII-X proposed for construction in Germany. The code is able to handle the many harmonics required to capture the high mode number of these instabilities. NSTAB builds on the highly successful BETAS code, which applies the spectral method to a flux coordinate formulation of the variational principle associated with the MHD equilibrium equations. However, a new residue condition for the location of the magnetic axis has been developed and implemented. This condition is based on the weak formulation of the equations and imposes no constraints on the inner flux surfaces.
- Research Organization:
- New York Univ., NY (United States). Courant Inst. of Mathematical Sciences
- Sponsoring Organization:
- USDOE, Washington, DC (United States); Department of Defense, Washington, DC (United States); National Science Foundation, Washington, DC (United States)
- DOE Contract Number:
- FG02-86ER53223
- OSTI ID:
- 10113306
- Report Number(s):
- DOE/ER/53223-193; ON: DE93004541; CNN: Grant: AFOSR-91-0042; Grant: DMS-8922805; Grant: DMS-900002P
- Resource Relation:
- Other Information: PBD: Sep 1992
- Country of Publication:
- United States
- Language:
- English
Similar Records
A high-performance spectral code for nonlinear MHD stability
MHD stability of the MHH2 stellarator