US Dept. of Energy 08.18.17
New supercomputing capabilities help understand how to cope with large-scale instabilities in tokamaks.
A cutaway of the lower divertor region of the National Spherical Torus Experiment in a simulation of a disruption. The red and blue colors indicate electrical “halo” currents into the walls of the tokamak. The rope-like structures show the paths of three magnetic field lines that intersect the walls.
The Science
Inside a fusion experiment, where scientists study the reactions at the heart of our sun, disruptions — large-scale instabilities of the plasma — cause rapid and complete loss of magnetic confinement. Models of fusion plasmas now combine advanced numerical methods with high-performance computing capabilities. The result? Scientists can explore the causes and dynamics of disruptions in unprecedented detail.
The Impact
Disruptions pose one of the most significant challenges to designing a fusion reactor. During these events, electrical currents arising in the walls create significant forces that can damage the walls of the tokamak vessel. Now, scientists can model these currents in a fully three-dimensional geometry, with realistic plasma parameters. The results can lead to strategies that avoid and mitigate disruptions in future reactor-sized devices.
Summary
The tokamak is an efficient design for confining superheated plasmas with magnetic fields because much of the magnetic field is produced by electrical currents in the plasma. This advantage can become a liability, because perturbations to the plasma current can reduce the magnetic field in a self-reinforcing cycle, causing rapid loss of confinement. Moreover, these disruptions impose strong electromagnetic forces and heat loads, posing a major challenge to successful operation of a tokamak reactor.
Researchers are now carrying out fully three-dimensional simulations of large-scale instabilities in the NSTX and DIII-D tokamaks. These simulations use the M3D-C1 code, which models the plasma as an electrically conducting fluid. New high-fidelity capabilities in the code show the electrical “halo” currents that can lead to disruptions flowing into and through the walls of the tokamak. And further simulations of vertical displacement events, which often cause or accompany disruptions, show that violent secondary instabilities may develop as the plasma is pushed against the vessel wall.
These secondary instabilities generally lead to a three-dimensional distribution of halo current, which consists of symmetric and asymmetric components. Asymmetric currents may produce forces that are particularly damaging to the tokamak vessel. Fortunately, in these simulations the asymmetric component remains localized and strongly subdominant to the symmetric component, even in cases that exhibit a strongly growing secondary instability. The simulations also show that cooling the plasma before or during the vertical displacement event may further suppress the instabilities that lead to asymmetric current. Future work will model disruptions initiated by other instabilities in which the asymmetric component of the halo currents is expected to be larger.
Contact
Nate Ferraro
Princeton Plasma Physics Laboratory
nferraro@pppl.gov
Funding
This work was supported by the U.S. Department of Energy (DOE) under contract DE-AC02-09-CH11466. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science user facility supported by the DOE Office of Science under contract DE-AC02-05CH11231. The research was performed as part of the activities of the Scientific Discovery through Advanced Computing Center for Extended-Magnetohydrodynamic Modeling.
Publications
N.M. Ferraro, S.C. Jardin, L.L. Lao, M.S. Shephard, and F. Zhang, “Multi-region approach to free-boundary three-dimensional tokamak equilibria and resistive wall instabilities” Physics of Plasmas 23, 056114 (2016). [DOI: 10.1063/1.4948722]
D. Pfefferlé, N. Ferraro, S. Jardin, and A. Bhattacharjee, “Fully 3D modeling of tokamak vertical displacement events with realistic parameters.” APS Division of Plasma Physics Meeting, abstract GP10.086. San Jose, California (2016).