-
Flexible Stellarator Physics Facility
Authors:
F. I. Parra,
S. -G. Baek,
M. Churchill,
D. R. Demers,
B. Dudson,
N. M. Ferraro,
B. Geiger,
S. Gerhardt,
K. C. Hammond,
S. Hudson,
R. Jorge,
E. Kolemen,
D. M. Kriete,
S. T. A. Kumar,
M. Landreman,
C. Lowe,
D. A. Maurer,
F. Nespoli,
N. Pablant,
M. J. Pueschel,
A. Punjabi,
J. A. Schwartz,
C. P. S. Swanson,
A. M. Wright
Abstract:
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Robustness of high-$β$ W7-X plasmas against ideal ballooning instability
Authors:
Yao Zhou,
K. Aleynikova,
Chang Liu,
N. M. Ferraro
Abstract:
To examine the robustness of the designed 5% $β$-limit for high-performance operation in the W7-X stellarator, we undertake nonlinear magnetohydrodynamic (MHD) simulations of pressure-driven instabilities using the M3D-$C^1$ code. Consistent with linear analyses, ideal ballooning instabilities occur as $β$ exceeds 5% in the standard configuration. Nonetheless, the modes saturate nonlinearly at rel…
▽ More
To examine the robustness of the designed 5% $β$-limit for high-performance operation in the W7-X stellarator, we undertake nonlinear magnetohydrodynamic (MHD) simulations of pressure-driven instabilities using the M3D-$C^1$ code. Consistent with linear analyses, ideal ballooning instabilities occur as $β$ exceeds 5% in the standard configuration. Nonetheless, the modes saturate nonlinearly at relatively low levels without triggering large-scale crashes, even though confinement degradation worsens as $β$ increases. In contrast, in an alternative configuration with nearly zero magnetic shear, ideal interchange modes induce a pressure crash at $β= 1\%$. These results suggest that the standard W7-X configuration might have a soft $β$-limit that is fairly immune to major MHD events, and that interchange instabilities tend to be more detrimental than ideal ballooning in stellarators.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
Nonlinear magnetohydrodynamic modeling of current-drive-induced sawtooth-like crashes in the W7-X stellarator
Authors:
Yao Zhou,
K. Aleynikova,
N. M. Ferraro
Abstract:
Sawtooth-like core electron temperature crashes have been observed in W7-X experiments with electron cyclotron current drive. We present nonlinear single-fluid magnetohydrodynamic simulations of this phenomenon using the newly developed stellarator modeling capability of the M3D-$C^1$ code. The near-axis current drive gives rise to two $ι=1$ resonances in the equilibrium rotational transform profi…
▽ More
Sawtooth-like core electron temperature crashes have been observed in W7-X experiments with electron cyclotron current drive. We present nonlinear single-fluid magnetohydrodynamic simulations of this phenomenon using the newly developed stellarator modeling capability of the M3D-$C^1$ code. The near-axis current drive gives rise to two $ι=1$ resonances in the equilibrium rotational transform profile so that two consecutive $(1,1)$ internal kink modes are seen in the simulations. A small-amplitude crash at the inner resonance occurs first, which may correspond to the sawtooth precursors observed in the experiments. A bigger crash at the outer resonance then flattens the core temperature profile, which shows semi-quantitative agreements with experimental measurements on certain metrics such as the crash amplitude and the inversion radius of the temperature change. These results illustrate a likely mechanism of the current-drive-induced sawtooth-like crashes in W7-X and, to some extent, validate the stellarator modeling capability of M3D-$C^1$.
△ Less
Submitted 16 March, 2023; v1 submitted 18 July, 2022;
originally announced July 2022.
-
Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration
Authors:
Chang Liu,
Stephen C. Jardin,
Hong Qin,
Jianyuan Xiao,
Nathaniel M. Ferraro,
Joshua Breslau
Abstract:
The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing…
▽ More
The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the $δf$ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes has been achieved.
△ Less
Submitted 28 July, 2021;
originally announced July 2021.
-
Modeling of carbon pellets disruption mitigation in an NSTX-U plasma
Authors:
C. F. Clauser,
S. C. Jardin,
R. Raman,
B. C. Lyons,
N. M. Ferraro
Abstract:
Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted in an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a $1\,$mm rad…
▽ More
Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted in an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a $1\,$mm radius vitreous carbon pellet (~ 3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first' simulations for disruption mitigation in NSTX-U and other devices, such as ITER.
△ Less
Submitted 28 June, 2021;
originally announced June 2021.
-
Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry
Authors:
Yao Zhou,
N. M. Ferraro,
S. C. Jardin,
H. R. Strauss
Abstract:
The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-$C^1$ code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-$C^1$. A $C^1$ coordinate mapping connects the logical domain to…
▽ More
The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-$C^1$ code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-$C^1$. A $C^1$ coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-$C^1$ extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of stellarator plasmas with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.
△ Less
Submitted 10 July, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Self-consistent simulation of resistive kink instabilities with runaway electrons
Authors:
Chang Liu,
Chen Zhao,
Stephen C. Jardin,
Nathaniel Ferraro,
Carlos Paz-Soldan,
Yueqiang Liu,
Brendan C. Lyons
Abstract:
A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris…
▽ More
A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.
△ Less
Submitted 21 April, 2021;
originally announced April 2021.
-
A Community Plan for Fusion Energy and Discovery Plasma Sciences
Authors:
Scott Baalrud,
Nathaniel Ferraro,
Lauren Garrison,
Nathan Howard,
Carolyn Kuranz,
John Sarff,
Wayne Solomon
Abstract:
This document is the final report of the Community Planning Process (CPP) that describes a comprehensive plan to deliver fusion energy and to advance plasma science. The CPP was initiated by the executive committee of the American Physical Society Division of Plasma Physics (APS DPP) to help the Fusion Energy Sciences Advisory Committee (FESAC) fulfill a charge from the U.S. Department of Energy (…
▽ More
This document is the final report of the Community Planning Process (CPP) that describes a comprehensive plan to deliver fusion energy and to advance plasma science. The CPP was initiated by the executive committee of the American Physical Society Division of Plasma Physics (APS DPP) to help the Fusion Energy Sciences Advisory Committee (FESAC) fulfill a charge from the U.S. Department of Energy (DOE) to develop a strategic plan for the DOE Office of Fusion Energy Sciences (FES). In this charge, dated Nov 30, 2018, DOE Deputy Director for Science Dr. Stephen Binkley requested that FESAC "undertake a new long range strategic planning activity for the Fusion Energy Sciences (FES) program. The strategic planning activity to encompass the entire FES research portfolio (namely, burning plasma science and discovery plasma science) should identify and prioritize the research required to advance both the scientific foundation needed to develop a fusion energy source, as well as the broader FES mission to steward plasma science." The CPP represents the first phase in developing a long range strategic plan for FES, and will serve as the basis for the second phase activity conducted by FESAC. It is worth noting that enacting the full scope of the recommendations in the strategic plan in this document will require suitable partnerships with other offices and governmental agencies, as well as with private industry and international partners.
△ Less
Submitted 9 November, 2020;
originally announced November 2020.
-
Structure and overstability of resistive modes with runaway electrons
Authors:
Chang Liu,
Chen Zhao,
Stephen C. Jardin,
Amitava Bhattacharjee,
Dylan P. Brennan,
Nathanial M. Ferraro
Abstract:
We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge wh…
▽ More
We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfvén velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C1 linear simulation, and good agreement is found in both cases.
△ Less
Submitted 9 September, 2020; v1 submitted 15 June, 2020;
originally announced June 2020.
-
Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry
Authors:
Robert Hager,
C. S. Chang,
N. M. Ferraro,
R. Nazikian
Abstract:
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking…
▽ More
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.
△ Less
Submitted 16 March, 2020;
originally announced March 2020.
-
Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C$^1$, NIMROD and JOREK
Authors:
I. Krebs,
F. J. Artola,
C. R. Sovinec,
S. C. Jardin,
K. J. Bunkers,
M. Hoelzl,
N. M. Ferraro
Abstract:
A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magneto-hydrodynamic codes M3D-C$^1$, JOREK and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wa…
▽ More
A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magneto-hydrodynamic codes M3D-C$^1$, JOREK and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. The benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C$^1$ simulations as well as from the linear phase of axisymmetric nonlinear JOREK, NIMROD and M3D-C$^1$ simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared and excellent agreement is found regarding plasma location and plasma currents as well as eddy and halo currents in the wall.
△ Less
Submitted 19 August, 2019; v1 submitted 6 August, 2019;
originally announced August 2019.
-
Observations of Heteroclinic Bifurcations in Resistive MHD Simulations of the Plasma Response to Resonant Magnetic Perturbations
Authors:
T. E. Evans,
W. Wu,
G. P. Canal,
N. M. Ferraro
Abstract:
A new class of static magnetohydrodynamic (MHD) magnetic island bifurcations is identified in rotating spherical tokamak plasmas during single- and two-fluid resistive MHD simulations. As the magnitude of an externally applied non-axisymmetric magnetic field perturbation is increased in these simulations, the internal flux surfaces that make up a sub-set of the resonant helical magnetic islands in…
▽ More
A new class of static magnetohydrodynamic (MHD) magnetic island bifurcations is identified in rotating spherical tokamak plasmas during single- and two-fluid resistive MHD simulations. As the magnitude of an externally applied non-axisymmetric magnetic field perturbation is increased in these simulations, the internal flux surfaces that make up a sub-set of the resonant helical magnetic islands in the plasma gradually elongate and undergo heteroclinic bifurcations. The bifurcation results in the creation of a new set of hyperbolic-elliptic fixed points as predicted by the Poincaré-Birkoff fixed point theorem. Field line calculations without including the resistive MHD plasma response to the applied perturbation field do not undergo this class of bifurcations indicating the importance of plasma self-organization in the bifurcation process.
△ Less
Submitted 12 June, 2018; v1 submitted 25 May, 2018;
originally announced May 2018.
-
Magnetic flux pumping in 3D nonlinear magnetohydrodynamic simulations
Authors:
I. Krebs,
S. C. Jardin,
S. Günter,
K. Lackner,
M. Hoelzl,
E. Strumberger,
N. Ferraro
Abstract:
A self-regulating magnetic flux pumping mechanism in tokamaks that maintains the core safety factor at $q\approx 1$, thus preventing sawteeth, is analyzed in nonlinear 3D magnetohydrodynamic simulations using the M3D-C$^1$ code. In these simulations, the most important mechanism responsible for the flux pumping is that a saturated $(m=1,n=1)$ quasi-interchange instability generates an effective ne…
▽ More
A self-regulating magnetic flux pumping mechanism in tokamaks that maintains the core safety factor at $q\approx 1$, thus preventing sawteeth, is analyzed in nonlinear 3D magnetohydrodynamic simulations using the M3D-C$^1$ code. In these simulations, the most important mechanism responsible for the flux pumping is that a saturated $(m=1,n=1)$ quasi-interchange instability generates an effective negative loop voltage in the plasma center via a dynamo effect. It is shown that sawtoothing is prevented in the simulations if $β$ is sufficiently high to provide the necessary drive for the $(m=1,n=1)$ instability that generates the dynamo loop voltage. The necessary amount of dynamo loop voltage is determined by the tendency of the current density profile to centrally peak which, in our simulations, is controlled by the peakedness of the applied heat source profile.
△ Less
Submitted 20 June, 2017;
originally announced June 2017.
