Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling
Abstract
Neoclassical tearing modes are macroscopic ($$L$$ ~ 1 m) instabilities in magnetic fusion experiments; if unchecked, these modes degrade plasma performance and may catastrophically destroy plasma confinement by inducing a disruption. Fortunately, the use of properly tuned and directed radiofrequency waves ($$λ$$ ~ 1 mm) can eliminate these modes. Numerical modeling of this difficult multiscale problem requires the integration of separate mathematical models for each length and time scale [T. G. Jenkins, S. E. Kruger, Phys. Plasmas 19 (2012) 122508]; the extended MHD model captures macroscopic plasma evolution while the RF model tracks the flow and deposition of injected RF power through the evolving plasma profiles. The scale separation enables use of the eikonal (ray-tracing) approximation to model the RF wave propagation. In this work we demonstrate a technique, based on methods of computational geometry, for mapping the ensuing RF data (associated with discrete ray trajectories) onto the finite-element/pseudospectral grid that is used to model the extended MHD physics. In the new representation, the RF data can then be used to construct source terms in the equations of the extended MHD model, enabling quantitative modeling of RF-induced tearing mode stabilization. Though our specific implementation uses the NIMROD extended MHD [C. R. Sovinec et al., J. Comp. Phys. 195 (2004) 355] and GENRAY RF [A. P. Smirnov, R. W. Harvey, K. Kupfer, Bull. Amer. Phys. Soc. 39 (1994) 1626] codes, the approach presented can be applied more generally to any code coupling requiring the mapping of ray tracing data onto Eulerian grids.
- Authors:
-
- Tech-X Corporation, Boulder, CO (United States)
- Utah State Univ., Logan, UT (United States)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Univ. of California, Oakland, CA (United States); UT-Battelle LLC/ORNL, Oak Ridge, TN (United States); Univ. of Wisconsin, Madison, WI (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC)
- OSTI Identifier:
- 1565301
- Alternate Identifier(s):
- OSTI ID: 1249981
- Grant/Contract Number:
- AC02-05CH11231; AC05-00OR22725; FC02-06ER54899
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 297; Journal Issue: C; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; computer science; physics; computational geometry; ray tracing; finite element methods; pseudospectral methods; magnetohydrodynamics; radiofrequency waves
Citation Formats
Jenkins, Thomas G., and Held, Eric D. Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling. United States: N. p., 2015.
Web. doi:10.1016/j.jcp.2015.05.035.
Jenkins, Thomas G., & Held, Eric D. Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling. United States. https://doi.org/10.1016/j.jcp.2015.05.035
Jenkins, Thomas G., and Held, Eric D. Tue .
"Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling". United States. https://doi.org/10.1016/j.jcp.2015.05.035. https://www.osti.gov/servlets/purl/1565301.
@article{osti_1565301,
title = {Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling},
author = {Jenkins, Thomas G. and Held, Eric D.},
abstractNote = {Neoclassical tearing modes are macroscopic ($L$ ~ 1 m) instabilities in magnetic fusion experiments; if unchecked, these modes degrade plasma performance and may catastrophically destroy plasma confinement by inducing a disruption. Fortunately, the use of properly tuned and directed radiofrequency waves ($λ$ ~ 1 mm) can eliminate these modes. Numerical modeling of this difficult multiscale problem requires the integration of separate mathematical models for each length and time scale [T. G. Jenkins, S. E. Kruger, Phys. Plasmas 19 (2012) 122508]; the extended MHD model captures macroscopic plasma evolution while the RF model tracks the flow and deposition of injected RF power through the evolving plasma profiles. The scale separation enables use of the eikonal (ray-tracing) approximation to model the RF wave propagation. In this work we demonstrate a technique, based on methods of computational geometry, for mapping the ensuing RF data (associated with discrete ray trajectories) onto the finite-element/pseudospectral grid that is used to model the extended MHD physics. In the new representation, the RF data can then be used to construct source terms in the equations of the extended MHD model, enabling quantitative modeling of RF-induced tearing mode stabilization. Though our specific implementation uses the NIMROD extended MHD [C. R. Sovinec et al., J. Comp. Phys. 195 (2004) 355] and GENRAY RF [A. P. Smirnov, R. W. Harvey, K. Kupfer, Bull. Amer. Phys. Soc. 39 (1994) 1626] codes, the approach presented can be applied more generally to any code coupling requiring the mapping of ray tracing data onto Eulerian grids.},
doi = {10.1016/j.jcp.2015.05.035},
journal = {Journal of Computational Physics},
number = C,
volume = 297,
place = {United States},
year = {Tue Sep 15 00:00:00 EDT 2015},
month = {Tue Sep 15 00:00:00 EDT 2015}
}
Web of Science
Works referenced in this record:
Neoclassical tearing modes and their control
journal, May 2006
- La Haye, R. J.
- Physics of Plasmas, Vol. 13, Issue 5
Observation of Nonlinear Neoclassical Pressure-Gradient–Driven Tearing Modes in TFTR
journal, June 1995
- Chang, Z.; Callen, J. D.; Fredrickson, E. D.
- Physical Review Letters, Vol. 74, Issue 23
Transport effects of low (m,n) MHD modes on TFTR supershots
journal, October 1994
- Chang, Z.; Fredrickson, E. D.; Callen, J. D.
- Nuclear Fusion, Vol. 34, Issue 10
Influence of neoclassical tearing modes on energy confinement
journal, January 1999
- Günter, S.; Gude, A.; Maraschek, M.
- Plasma Physics and Controlled Fusion, Vol. 41, Issue 6
MHD instabilities near the β limit in the doublet III-D tokamak
journal, March 1989
- Strait, E. J.; Lao, L.; Kellman, A. G.
- Physical Review Letters, Vol. 62, Issue 11
Beta limits in long-pulse tokamak discharges
journal, May 1997
- Sauter, O.; La Haye, R. J.; Chang, Z.
- Physics of Plasmas, Vol. 4, Issue 5
High beta tokamak operation in DIII-D limited at low density/collisionality by resistive tearing modes
journal, March 1997
- Haye, R. J. La; Lao, L. L.; Strait, E. J.
- Nuclear Fusion, Vol. 37, Issue 3
Chapter 1: Overview and summary
journal, December 1999
- Editors, ITER Physics Basis; Co-Chairs, ITER Physics Expert Group Chairs an; Unit, ITER Joint Central Team and Physics
- Nuclear Fusion, Vol. 39, Issue 12
Experiments on neoclassical tearing mode stabilization by ECCD in ASDEX Upgrade
journal, May 1999
- Zohm, H.; Gantenbein, G.; Giruzzi, G.
- Nuclear Fusion, Vol. 39, Issue 5
Complete Suppression of Neoclassical Tearing Modes with Current Drive at the Electron-Cyclotron-Resonance Frequency in ASDEX Upgrade Tokamak
journal, August 2000
- Gantenbein, G.; Zohm, H.; Giruzzi, G.
- Physical Review Letters, Vol. 85, Issue 6
Complete stabilization of a tearing mode in steady state high-β p H-mode discharges by the first harmonic electron cyclotron heating/current drive on JT-60U
journal, November 2000
- Isayama, A.; Kamada, Y.; Ide, S.
- Plasma Physics and Controlled Fusion, Vol. 42, Issue 12
Control of neoclassical tearing modes in DIII–D
journal, May 2002
- La Haye, R. J.; Günter, S.; Humphreys, D. A.
- Physics of Plasmas, Vol. 9, Issue 5
Active control for stabilization of neoclassical tearing modes
journal, May 2006
- Humphreys, D. A.; Ferron, J. R.; La Haye, R. J.
- Physics of Plasmas, Vol. 13, Issue 5
Stabilization and prevention of the 2/1 neoclassical tearing mode for improved performance in DIII-D
journal, April 2007
- Prater, R.; Haye, R. J. La; Luce, T. C.
- Nuclear Fusion, Vol. 47, Issue 5
Achievement of high fusion triple product, steady-state sustainment and real-time NTM stabilization in high-Â p ELMy H-mode discharges in JT-60U
journal, October 2003
- Isayama, A.; Kamada, Y.; Hayashi, N.
- Nuclear Fusion, Vol. 43, Issue 10
Enhancement of the Stabilization Efficiency of a Neoclassical Magnetic Island by Modulated Electron Cyclotron Current Drive in the ASDEX Upgrade Tokamak
journal, January 2007
- Maraschek, M.; Gantenbein, G.; Yu, Q.
- Physical Review Letters, Vol. 98, Issue 2
A closure scheme for modeling rf modifications to the fluid equations
journal, November 2009
- Hegna, C. C.; Callen, J. D.
- Physics of Plasmas, Vol. 16, Issue 11
Fluid and drift-kinetic description of a magnetized plasma with low collisionality and slow dynamics orderings. I. Electron theory
journal, August 2010
- Ramos, J. J.
- Physics of Plasmas, Vol. 17, Issue 8
Fluid and drift-kinetic description of a magnetized plasma with low collisionality and slow dynamics orderings. II. Ion theory
journal, October 2011
- Ramos, J. J.
- Physics of Plasmas, Vol. 18, Issue 10
Fluid equations in the presence of electron cyclotron current drive
journal, December 2012
- Jenkins, Thomas G.; Kruger, Scott E.
- Physics of Plasmas, Vol. 19, Issue 12
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
Dynamics of the major disruption of a DIII-D plasma
journal, May 2005
- Kruger, S. E.; Schnack, D. D.; Sovinec, C. R.
- Physics of Plasmas, Vol. 12, Issue 5
Kinetic Effects of Energetic Particles on Resistive MHD Stability
journal, March 2009
- Takahashi, R.; Brennan, D. P.; Kim, C. C.
- Physical Review Letters, Vol. 102, Issue 13
Edge localized linear ideal magnetohydrodynamic instability studies in an extended-magnetohydrodynamic code
journal, March 2010
- Burke, B. J.; Kruger, S. E.; Hegna, C. C.
- Physics of Plasmas, Vol. 17, Issue 3
Runaway electron confinement modelling for rapid shutdown scenarios in DIII-D, Alcator C-Mod and ITER
journal, May 2011
- Izzo, V. A.; Hollmann, E. M.; James, A. N.
- Nuclear Fusion, Vol. 51, Issue 6
First-order finite-Larmor-radius effects on magnetic tearing in pinch configurations
journal, April 2011
- King, J. R.; Sovinec, C. R.; Mirnov, V. V.
- Physics of Plasmas, Vol. 18, Issue 4
Nonlocal closures for plasma fluid simulations
journal, May 2004
- Held, E. D.; Callen, J. D.; Hegna, C. C.
- Physics of Plasmas, Vol. 11, Issue 5
Lower hybrid current drive experiments on Alcator C-Mod: Comparison with theory and simulation
journal, May 2008
- Bonoli, P. T.; Ko, J.; Parker, R.
- Physics of Plasmas, Vol. 15, Issue 5
Investigation of lower hybrid physics through power modulation experiments on Alcator C-Mod
journal, May 2011
- Schmidt, A.; Bonoli, P. T.; Meneghini, O.
- Physics of Plasmas, Vol. 18, Issue 5
Off-midplane launch of electron Bernstein waves for current drive in overdense plasmas
journal, May 2000
- Forest, C. B.; Chattopadhyay, P. K.; Harvey, R. W.
- Physics of Plasmas, Vol. 7, Issue 5
Power deposition by mode converted electron Bernstein waves in the DIII-D `heat pinch' experiments
journal, May 2001
- Forest, C. B.; Harvey, R. W.; Smirnov, A. P.
- Nuclear Fusion, Vol. 41, Issue 5
Finite-Larmor-radius kinetic theory of a magnetized plasma in the macroscopic flow reference frame
journal, August 2008
- Ramos, J. J.
- Physics of Plasmas, Vol. 15, Issue 8
Computational modeling of fully ionized magnetized plasmas using the fluid approximation
journal, May 2006
- Schnack, D. D.; Barnes, D. C.; Brennan, D. P.
- Physics of Plasmas, Vol. 13, Issue 5
Analysis of a mixed semi-implicit/implicit algorithm for low-frequency two-fluid plasma modeling
journal, August 2010
- Sovinec, C. R.; King, J. R.
- Journal of Computational Physics, Vol. 229, Issue 16
Calculating electron cyclotron current drive stabilization of resistive tearing modes in a nonlinear magnetohydrodynamic model
journal, January 2010
- Jenkins, Thomas G.; Kruger, Scott E.; Hegna, C. C.
- Physics of Plasmas, Vol. 17, Issue 1
Geometric optics in plasmas characterized by non-Hermitian dielectric tensors
journal, October 1980
- Friedland, L.; Bernstein, I. B.
- Physical Review A, Vol. 22, Issue 4
The quickhull algorithm for convex hulls
journal, December 1996
- Barber, C. Bradford; Dobkin, David P.; Huhdanpaa, Hannu
- ACM Transactions on Mathematical Software, Vol. 22, Issue 4
Development of ITER-relevant plasma control solutions at DIII-D
journal, July 2007
- Humphreys, D. A.; Ferron, J. R.; Bakhtiari, M.
- Nuclear Fusion, Vol. 47, Issue 8
Algorithm 660: QSHEP2D: Quadratic Shepard Method for Bivariate Interpolation of Scattered Data
journal, June 1988
- Renka, Robert J.
- ACM Transactions on Mathematical Software, Vol. 14, Issue 2
Algorithm 792: accuracy test of ACM algorithms for interpolation of scattered data in the plane
journal, March 1999
- Renka, Robert J.; Brown, Ron
- ACM Transactions on Mathematical Software, Vol. 25, Issue 1
Creating an Asymmetric Plasma Resistivity with Waves
journal, September 1980
- Fisch, N. J.; Boozer, A. H.
- Physical Review Letters, Vol. 45, Issue 9
Verification of continuum drift kinetic equation solvers in NIMROD
journal, March 2015
- Held, E. D.; Kruger, S. E.; Ji, J. -Y.
- Physics of Plasmas, Vol. 22, Issue 3
Exact linearized Coulomb collision operator in the moment expansion
journal, October 2006
- Ji, Jeong-Young; Held, Eric D.
- Physics of Plasmas, Vol. 13, Issue 10