skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Fully Parallel MHD Stability Code

Authors:
 [1];  [1]
  1. FAR-TECH, Inc., San Diego, CA (United States)
Publication Date:
Research Org.:
FAR-TECH, Inc., San Diego, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24)
OSTI Identifier:
1358076
Report Number(s):
DOE-FAR-TECH-07599
DOE Contract Number:
SC0007599
Type / Phase:
SBIR
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; MHD; stability; Resistive Wall Mode; MARS code; Parallelization; MPI; peeling-ballooning mode; MARS-F; MARS-K

Citation Formats

Svidzinski, Vladimir, and Galkin, Sergei A. Fully Parallel MHD Stability Code. United States: N. p., 2017. Web.
Svidzinski, Vladimir, & Galkin, Sergei A. Fully Parallel MHD Stability Code. United States.
Svidzinski, Vladimir, and Galkin, Sergei A. Sat . "Fully Parallel MHD Stability Code". United States. doi:.
@article{osti_1358076,
title = {Fully Parallel MHD Stability Code},
author = {Svidzinski, Vladimir and Galkin, Sergei A},
abstractNote = {},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Sat May 20 00:00:00 EDT 2017},
month = {Sat May 20 00:00:00 EDT 2017}
}

Technical Report:
This technical report may be protected. To request the document, click here.
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that may hold this item. Keep in mind that many technical reports are not cataloged in WorldCat.

Save / Share:
  • Utilizing a Galerkin procedure to calculate the vacuum contribution to the ideal MHD Lagrangian, the implementation of realistic boundary conditions are described in a linear stability code. The procedure permits calculation of the effect of arbitrary conducting structure on ideal MHD instabilities, as opposed to the prior use of an encircling shell. The passive stabilization of conducting coils on the tokamak vertical instability is calculated within the PEST code and gives excellent agreement with 2-D time dependent simulations of PDX.
  • A nonvariational approach for determining the ideal MHD stability of axisymmetric toroidal confinement systems is presented. The code (NOVA) employs cubic B-spline finite elements and Fourier expansion in a general flux coordinate (psi, theta, zeta) system. Better accuracy and faster convergence were obtained in comparison with the variational PEST and ERATO codes. The nonvariational approach can be extended to problems having non-Hermitian eigenmode equations where variational energy principles cannot be obtained.
  • 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 themore » 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.« less
  • 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 themore » 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.« less
  • Combustion and nitric oxide formation are analyzed for an MHD secondary combustor consisting of multi-parallel jets of air injected into a stream consisting of the primary combustor combustion product containing carbon monoxide. The minimum combustor length required for mixing and complete combustion in the limit of fast chemical kinetics is determined. This length is shorter than that required for uniform mixing. It is a moderate function of the jet spacing and is rather strongly dependent upon the stoichiometric air ratio. In the usual operating temperature range of 1600 to 2000/sup 0/K and for the lengths up to about 5m, themore » NO/sub x/ reactions in the secondary combustors are found to be nearly frozen. Therefore, ignition and combustion, rather than the nitric oxide formation, should be the secondary combustor design criterion.« less