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 »
- Authors:
-
- Fusion Theory and Computation, Inc., Kingston, WA (United States)
- 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}
}
Web of Science
Works referenced in this record:
Linear tearing mode stability equations for a low collisionality toroidal plasma
journal, December 2008
- Connor, J. W.; Hastie, R. J.; Helander, P.
- Plasma Physics and Controlled Fusion, Vol. 51, Issue 1
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
Gravitational instabilities in a compressible collisional plasma layer
journal, October 1964
- Coppi, Bruno
- Annals of Physics, Vol. 30, Issue 1
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
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
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
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
Vacuum calculations in azimuthally symmetric geometry
journal, June 1997
- Chance, M. S.
- Physics of Plasmas, Vol. 4, Issue 6
The direct criterion of Newcomb for the ideal MHD stability of an axisymmetric toroidal plasma
journal, July 2016
- Glasser, A. H.
- Physics of Plasmas, Vol. 23, Issue 7
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
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
Two-dimensional generalizations of the Newcomb equation
journal, April 1990
- Dewar, R. L.; Pletzer, A.
- Journal of Plasma Physics, Vol. 43, Issue 2
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
Convergent solutions of ordinary linear homogeneous differential equations in the neighborhood of an irregular singular point
journal, January 1955
- Turrittin, H. L.
- Acta Mathematica, Vol. 93, Issue 0
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
Resistive instabilities in a tokamak
journal, January 1976
- Glasser, A. H.; Greene, J. M.; Johnson, J. L.
- Physics of Fluids, Vol. 19, Issue 4
Improved poloidal convergence of the MARS code for MHD stability analysis
journal, January 1999
- Liu, D. H.; Bondeson, A.
- Computer Physics Communications, Vol. 116, Issue 1
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
Works referencing / citing this record:
Determination of the non-ideal response of a high temperature tokamak plasma to a static external magnetic perturbation via asymptotic matching
journal, July 2017
- Fitzpatrick, Richard
- Physics of Plasmas, Vol. 24, Issue 7
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
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
Two-fluid nonlinear theory of response of tokamak plasma to resonant magnetic perturbation
journal, November 2018
- Fitzpatrick, Richard
- Physics of Plasmas, Vol. 25, Issue 11
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
Magnetic polarization measurements of the multi-modal plasma response to 3D fields in the EAST tokamak
journal, May 2018
- Logan, N. C.; Cui, L.; Wang, H.
- Nuclear Fusion, Vol. 58, Issue 7
Reduced energetic particle transport models enable comprehensive time-dependent tokamak simulations
journal, August 2019
- Podestà, M.; Bardóczi, L.; Collins, C. S.
- Nuclear Fusion, Vol. 59, Issue 10
Two-Fluid Nonlinear Theory of Response of Tokamak Plasma to Resonant Magnetic Perturbation
text, January 2018
- Fitzpatrick, Richard
- arXiv