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

Title: Axisymmetric benchmarks of impurity dynamics in extended-magnetohydrodynamic simulations

Abstract

A verification benchmark has been carried out between the M3D-C1 and NIMROD extended-magnetohydrodynamics (xMHD) codes for simulations of impurity-induced disruption mitigation. Disruptions are a significant concern for future tokamaks, and high-fidelity simulations are required in order to ensure the success of disruption-mitigation techniques (e.g., shattered-pellet injection) in large-scale fusion reactors. Both MHD codes have been coupled to the KPRAD code for impurity dynamics. The codes show excellent agreement in four axisymmetric, nonlinear simulations, particularly during the thermal quench. This agreement is seen in the time histories of global plasma quantities such as thermal energy, radiated power, and total number of electrons, as well as 2D contours of temperature and current density. The simulations predict that, given the same number of atoms injected, argon quenches the plasma two-to-three times as fast as neon. Furthermore, the inclusion of temperature-dependent Spitzer resistivity causes the current to diffuse and to decay, inducing axisymmetric MHD instabilities that result in a current quench. As a result, this work represents an important verification of the coupled impurity and MHD models implemented in M3D-C1 and NIMROD, giving greater confidence in the ability of both codes to perform more sophisticated disruption-mitigation simulations.

Authors:
ORCiD logo [1];  [2];  [1];  [3];  [3];  [1];  [1];  [1]
  1. General Atomics, San Diego, CA (United States)
  2. SLS2 Consulting, San Diego, CA (United States)
  3. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Research Org.:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); General Atomics, San Diego, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21). Scientific Discovery through Advanced Computing (SciDAC)
OSTI Identifier:
1505779
Alternate Identifier(s):
OSTI ID: 1569257
Grant/Contract Number:  
AC02-05CH11231; FC02-04ER54698; FG02-95ER54309; SC0018109
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Plasma Physics and Controlled Fusion
Additional Journal Information:
Journal Volume: 61; Journal Issue: 6; Journal ID: ISSN 0741-3335
Publisher:
IOP Science
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; magnetohydrodynamics; impurities; disruptions; disruption mitigation; integrated modeling; benchmark

Citation Formats

Lyons, B. C., Kim, C. C., Liu, Y. Q., Ferraro, N. M., Jardin, S. C., McClenaghan, J., Parks, P. B., and Lao, L. L. Axisymmetric benchmarks of impurity dynamics in extended-magnetohydrodynamic simulations. United States: N. p., 2019. Web. doi:10.1088/1361-6587/ab0e42.
Lyons, B. C., Kim, C. C., Liu, Y. Q., Ferraro, N. M., Jardin, S. C., McClenaghan, J., Parks, P. B., & Lao, L. L. Axisymmetric benchmarks of impurity dynamics in extended-magnetohydrodynamic simulations. United States. doi:10.1088/1361-6587/ab0e42.
Lyons, B. C., Kim, C. C., Liu, Y. Q., Ferraro, N. M., Jardin, S. C., McClenaghan, J., Parks, P. B., and Lao, L. L. Tue . "Axisymmetric benchmarks of impurity dynamics in extended-magnetohydrodynamic simulations". United States. doi:10.1088/1361-6587/ab0e42. https://www.osti.gov/servlets/purl/1505779.
@article{osti_1505779,
title = {Axisymmetric benchmarks of impurity dynamics in extended-magnetohydrodynamic simulations},
author = {Lyons, B. C. and Kim, C. C. and Liu, Y. Q. and Ferraro, N. M. and Jardin, S. C. and McClenaghan, J. and Parks, P. B. and Lao, L. L.},
abstractNote = {A verification benchmark has been carried out between the M3D-C1 and NIMROD extended-magnetohydrodynamics (xMHD) codes for simulations of impurity-induced disruption mitigation. Disruptions are a significant concern for future tokamaks, and high-fidelity simulations are required in order to ensure the success of disruption-mitigation techniques (e.g., shattered-pellet injection) in large-scale fusion reactors. Both MHD codes have been coupled to the KPRAD code for impurity dynamics. The codes show excellent agreement in four axisymmetric, nonlinear simulations, particularly during the thermal quench. This agreement is seen in the time histories of global plasma quantities such as thermal energy, radiated power, and total number of electrons, as well as 2D contours of temperature and current density. The simulations predict that, given the same number of atoms injected, argon quenches the plasma two-to-three times as fast as neon. Furthermore, the inclusion of temperature-dependent Spitzer resistivity causes the current to diffuse and to decay, inducing axisymmetric MHD instabilities that result in a current quench. As a result, this work represents an important verification of the coupled impurity and MHD models implemented in M3D-C1 and NIMROD, giving greater confidence in the ability of both codes to perform more sophisticated disruption-mitigation simulations.},
doi = {10.1088/1361-6587/ab0e42},
journal = {Plasma Physics and Controlled Fusion},
issn = {0741-3335},
number = 6,
volume = 61,
place = {United States},
year = {2019},
month = {4}
}

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

Save / Share:

Works referenced in this record:

Experimental study of disruption mitigation using massive injection of noble gases on Tore Supra
journal, July 2010


LSODE and LSODI, two new initial value ordinary differnetial equation solvers
journal, December 1980


Nonlinear magnetohydrodynamics simulation using high-order finite elements
journal, March 2004

  • Sovinec, C. R.; Glasser, A. H.; Gianakon, T. A.
  • Journal of Computational Physics, Vol. 195, Issue 1
  • DOI: 10.1016/j.jcp.2003.10.004

MHD stability in X-point geometry: simulation of ELMs
journal, June 2007


Impurity mixing and radiation asymmetry in massive gas injection simulations of DIII-D
journal, May 2013


Magnetohydrodynamic simulations of massive gas injection into Alcator C-Mod and DIII-D plasmas
journal, May 2008

  • Izzo, V. A.; Whyte, D. G.; Granetz, R. S.
  • Physics of Plasmas, Vol. 15, Issue 5
  • DOI: 10.1063/1.2841526

Nonlinear reconnecting edge localized modes in current-carrying plasmas
journal, May 2017


Disruption mitigation studies in DIII-D
journal, May 1999

  • Taylor, P. L.; Kellman, A. G.; Evans, T. E.
  • Physics of Plasmas, Vol. 6, Issue 5
  • DOI: 10.1063/1.873445

Chapter 3: MHD stability, operational limits and disruptions
journal, December 1999

  • Mhd, ITER Physics Expert Group on Disrup; Editors, ITER Physics Basis
  • Nuclear Fusion, Vol. 39, Issue 12
  • DOI: 10.1088/0029-5515/39/12/303

Gas jet disruption mitigation studies on Alcator C-Mod
journal, October 2006


Demonstration of Tokamak Discharge Shutdown with Shell Pellet Payload Impurity Dispersal
journal, February 2019


The role of MHD in 3D aspects of massive gas injection
journal, June 2015


Runaway electron confinement modelling for rapid shutdown scenarios in DIII-D, Alcator C-Mod and ITER
journal, May 2011


3D two-temperature magnetohydrodynamic modeling of fast thermal quenches due to injected impurities in tokamaks
journal, November 2018


Low-Z Shell Pellet Experiments on DIII-D
conference, January 2009

  • Hollmann, E. M.; James, A. N.; Parks, P. B.
  • ATOMIC PROCESSES IN PLASMAS: Proceedings of the 16th International Conference on Atomic Processes in Plasmas, AIP Conference Proceedings
  • DOI: 10.1063/1.3241211

Dynamo-driven plasmoid formation from a current-sheet instability
journal, December 2016


Disruptions in ITER and strategies for their control and mitigation
journal, August 2015


Modeling of rapid shutdown in the DIII-D tokamak by core deposition of high-Z material
journal, June 2017

  • Izzo, V. A.; Parks, P. B.
  • Physics of Plasmas, Vol. 24, Issue 6
  • DOI: 10.1063/1.4990447

First demonstration of rapid shutdown using neon shattered pellet injection for thermal quench mitigation on DIII-D
journal, March 2016


Multiple timescale calculations of sawteeth and other global macroscopic dynamics of tokamak plasmas
journal, January 2012


Numerical Studies of Impurities in Fusion Plasmas
journal, March 1983


Study of plasma termination using high- Z noble gas puffing in the JT-60U tokamak
journal, April 2005


Analysis of a mixed semi-implicit/implicit algorithm for low-frequency two-fluid plasma modeling
journal, August 2010


Chapter 3: MHD stability, operational limits and disruptions
journal, June 2007


Bézier surfaces and finite elements for MHD simulations
journal, August 2008