DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Computation of resistive instabilities by matched asymptotic expansions

Abstract

Here, we present a method for determining the linear resistive magnetohydrodynamic (MHD) stability of an axisymmetric toroidal plasma, based on the method of matched asymptotic expansions. The plasma is partitioned into a set of ideal MHD outer regions, connected through resistive MHD inner regions about singular layers where q = m/n, with m and n toroidal mode numbers, respectively, and q the safety factor. The outer regions satisfy the ideal MHD equations with zero-frequency, which are identical to the Euler-Lagrange equations for minimizing the potential energy delta W. The solutions to these equations go to infinity at the singular surfaces. The inner regions satisfy the equations of motion of resistive MHD with a finite eigenvalue, resolving the singularity. Both outer and inner regions are solved numerically by newly developed singular Galerkin methods, using specialized basis functions. These solutions are matched asymptotically, providing a complex dispersion relation which is solved for global eigenvalues and eigenfunctions in full toroidal geometry. The dispersion relation may have multiple complex unstable roots, which are found by advanced root-finding methods. These methods are much faster and more robust than the previous numerical methods. The new methods are applicable to more challenging high-pressure and strongly shaped plasmamore » equilibria and generalizable to more realistic inner region dynamics. In the thermonuclear regime, where the outer and inner regions overlap, they are also much faster and more accurate than the straight-through methods, which treat the resistive MHD equations in the whole plasma volume.« less

Authors:
ORCiD logo [1];  [2];  [2]
  1. Fusion Theory and Computation, Inc., Kingston, WA (United States)
  2. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Research Org.:
Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States); Fusion Theory and Computation, Inc., Kingston, WA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES)
OSTI Identifier:
1340151
Alternate Identifier(s):
OSTI ID: 1332585; OSTI ID: 1418992
Report Number(s):
PPPL-5292
Journal ID: ISSN 1070-664X; TRN: US1701243
Grant/Contract Number:  
AC02-09CH11466; FG02-05ER54811; SC0016106
Resource Type:
Accepted Manuscript
Journal Name:
Physics of Plasmas
Additional Journal Information:
Journal Volume: 23; Journal Issue: 11; Journal ID: ISSN 1070-664X
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; 01 COAL, LIGNITE, AND PEAT; resistive MHD instabilities

Citation Formats

Glasser, A. H., Wang, Z. R., and Park, J. -K. Computation of resistive instabilities by matched asymptotic expansions. United States: N. p., 2016. Web. doi:10.1063/1.4967862.
Glasser, A. H., Wang, Z. R., & Park, J. -K. Computation of resistive instabilities by matched asymptotic expansions. United States. https://doi.org/10.1063/1.4967862
Glasser, A. H., Wang, Z. R., and Park, J. -K. Thu . "Computation of resistive instabilities by matched asymptotic expansions". United States. https://doi.org/10.1063/1.4967862. https://www.osti.gov/servlets/purl/1340151.
@article{osti_1340151,
title = {Computation of resistive instabilities by matched asymptotic expansions},
author = {Glasser, A. H. and Wang, Z. R. and Park, J. -K.},
abstractNote = {Here, we present a method for determining the linear resistive magnetohydrodynamic (MHD) stability of an axisymmetric toroidal plasma, based on the method of matched asymptotic expansions. The plasma is partitioned into a set of ideal MHD outer regions, connected through resistive MHD inner regions about singular layers where q = m/n, with m and n toroidal mode numbers, respectively, and q the safety factor. The outer regions satisfy the ideal MHD equations with zero-frequency, which are identical to the Euler-Lagrange equations for minimizing the potential energy delta W. The solutions to these equations go to infinity at the singular surfaces. The inner regions satisfy the equations of motion of resistive MHD with a finite eigenvalue, resolving the singularity. Both outer and inner regions are solved numerically by newly developed singular Galerkin methods, using specialized basis functions. These solutions are matched asymptotically, providing a complex dispersion relation which is solved for global eigenvalues and eigenfunctions in full toroidal geometry. The dispersion relation may have multiple complex unstable roots, which are found by advanced root-finding methods. These methods are much faster and more robust than the previous numerical methods. The new methods are applicable to more challenging high-pressure and strongly shaped plasma equilibria and generalizable to more realistic inner region dynamics. In the thermonuclear regime, where the outer and inner regions overlap, they are also much faster and more accurate than the straight-through methods, which treat the resistive MHD equations in the whole plasma volume.},
doi = {10.1063/1.4967862},
journal = {Physics of Plasmas},
number = 11,
volume = 23,
place = {United States},
year = {Thu Nov 17 00:00:00 EST 2016},
month = {Thu Nov 17 00:00:00 EST 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 24 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Linear tearing mode stability equations for a low collisionality toroidal plasma
journal, December 2008


Numerical solution of the resistive magnetohydrodynamic boundary layer equations
journal, January 1984

  • Glasser, A. H.; Jardin, S. C.; Tesauro, G.
  • Physics of Fluids, Vol. 27, Issue 5
  • DOI: 10.1063/1.864730

Gravitational instabilities in a compressible collisional plasma layer
journal, October 1964


LAPACK Users' Guide
software, January 1999


Resistive instabilities in general toroidal plasma configurations
journal, January 1975

  • Glasser, A. H.; Greene, J. M.; Johnson, J. L.
  • Physics of Fluids, Vol. 18, Issue 7
  • DOI: 10.1063/1.861224

Linear Stability of Resistive MHD Modes: Axisymmetric Toroidal Computation of the Outer Region Matching Data
journal, December 1994

  • Pletzer, A.; Bondeson, A.; Dewar, R. L.
  • Journal of Computational Physics, Vol. 115, Issue 2
  • DOI: 10.1006/jcph.1994.1215

Full toroidal plasma response to externally applied nonaxisymmetric magnetic fields
journal, December 2010

  • Liu, Yueqiang; Kirk, A.; Nardon, E.
  • Physics of Plasmas, Vol. 17, Issue 12
  • DOI: 10.1063/1.3526677

A new general method for solving the resistive inner layer problem
journal, September 2002

  • Galkin, S. A.; Turnbull, A. D.; Greene, J. M.
  • Physics of Plasmas, Vol. 9, Issue 9
  • DOI: 10.1063/1.1501092

Vacuum calculations in azimuthally symmetric geometry
journal, June 1997


The direct criterion of Newcomb for the ideal MHD stability of an axisymmetric toroidal plasma
journal, July 2016


Finite-Resistivity Instabilities of a Sheet Pinch
journal, January 1963

  • Furth, Harold P.; Killeen, John; Rosenbluth, Marshall N.
  • Physics of Fluids, Vol. 6, Issue 4
  • DOI: 10.1063/1.1706761

A new formulation of the resistive magnetohydrodynamics stability problem for finite β toroidal plasmas
journal, January 2000

  • Galkin, S. A.; Turnbull, A. D.; Greene, J. M.
  • Physics of Plasmas, Vol. 7, Issue 10
  • DOI: 10.1063/1.1289688

Two-dimensional generalizations of the Newcomb equation
journal, April 1990


Modification of Δ′ by magnetic feedback and kinetic effects
journal, September 2012

  • Liu, Yueqiang; Hastie, R. J.; Hender, T. C.
  • Physics of Plasmas, Vol. 19, Issue 9
  • DOI: 10.1063/1.4754281

Feedback stabilization of nonaxisymmetric resistive wall modes in tokamaks. I. Electromagnetic model
journal, September 2000

  • Liu, Y. Q.; Bondeson, A.; Fransson, C. M.
  • Physics of Plasmas, Vol. 7, Issue 9
  • DOI: 10.1063/1.1287744

Resistive instabilities in a tokamak
journal, January 1976

  • Glasser, A. H.; Greene, J. M.; Johnson, J. L.
  • Physics of Fluids, Vol. 19, Issue 4
  • DOI: 10.1063/1.861490

Improved poloidal convergence of the MARS code for MHD stability analysis
journal, January 1999


Calculation of the vacuum Green’s function valid even for high toroidal mode numbers in tokamaks
journal, January 2007

  • Chance, M. S.; Turnbull, A. D.; Snyder, P. B.
  • Journal of Computational Physics, Vol. 221, Issue 1
  • DOI: 10.1016/j.jcp.2006.06.025

Works referencing / citing this record:

A Riccati solution for the ideal MHD plasma response with applications to real-time stability control
journal, March 2018

  • Glasser, Alexander S.; Kolemen, Egemen; Glasser, A. H.
  • Physics of Plasmas, Vol. 25, Issue 3
  • DOI: 10.1063/1.5007042

A robust solution for the resistive MHD toroidal Δ′ matrix in near real-time
journal, August 2018

  • Glasser, Alexander S.; Kolemen, Egemen
  • Physics of Plasmas, Vol. 25, Issue 8
  • DOI: 10.1063/1.5029477

Two-fluid nonlinear theory of response of tokamak plasma to resonant magnetic perturbation
journal, November 2018


Asymptotic solutions and convergence studies of the resistive inner region equations
journal, January 2020

  • Glasser, A. H.; Wang, Z. R.
  • Physics of Plasmas, Vol. 27, Issue 1
  • DOI: 10.1063/1.5134999

Magnetic polarization measurements of the multi-modal plasma response to 3D fields in the EAST tokamak
journal, May 2018


Reduced energetic particle transport models enable comprehensive time-dependent tokamak simulations
journal, August 2019