-
Growing Neoclassical Tearing Modes Seeded via Transient-Induced-Multimode Interactions
Authors:
E. C. Howell,
J. R. King,
J. D. Callen,
R. J. La Haye,
R. S. Wilcox,
S. E. Kruger
Abstract:
Nonlinear extended MHD simulations demonstrating seeding of neoclassical tearing modes (NTMs) via MHD-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inher…
▽ More
Nonlinear extended MHD simulations demonstrating seeding of neoclassical tearing modes (NTMs) via MHD-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge [R.J. La Haye, et al., Proceedings IAEA FEC 2020]. The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase like that observed experimentally. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth which is well described by the nonlinear modified Rutherford equation. This work highlights the role of nonlinear mode coupling in seeding NTMs.
△ Less
Submitted 22 October, 2021;
originally announced October 2021.
-
The deteriorating effect of plasma density fluctuations on microwave beam quality
Authors:
Alf Köhn,
Max E. Austin,
Michael W. Brookman,
Kenneth W. Gentle,
Lorenzo Guidi,
Eberhard Holzhauer,
Rob J. La Haye,
Jarrod B. Leddy,
Omar Maj,
Craig C. Petty,
Emanuele Poli,
Terry L. Rhodes,
Antti Snicker,
Matthew B. Thomas,
Roddy G. L. Vann,
Hannes Weber
Abstract:
Turbulent plasma edge density fluctuations can broaden a traversing microwave beam degrading its quality. This can be a problem for scenarios relying on a high spatial localization of the deposition of injected microwave power, like controlling MHD instabilities. Here we present numerical estimations of the scattering of a microwave beam by density fluctuations over a large parameter range, includ…
▽ More
Turbulent plasma edge density fluctuations can broaden a traversing microwave beam degrading its quality. This can be a problem for scenarios relying on a high spatial localization of the deposition of injected microwave power, like controlling MHD instabilities. Here we present numerical estimations of the scattering of a microwave beam by density fluctuations over a large parameter range, including extrapolations to ITER. Two codes are used, the full-wave code IPF-FDMC and the wave kinetic equation solver WKBeam. A successful comparison between beam broadening obtained from DIII-D experiments and corresponding full-wave simulations is shown.
△ Less
Submitted 24 May, 2018;
originally announced May 2018.
-
Feedforward and feedback control of locked mode phase and rotation in DIII-D with application to modulated ECCD experiments
Authors:
Wilkie Choi,
R. J. La Haye,
M. J. Lanctot,
K. E. J. Olofsson,
E. J. Strait,
R. Sweeney,
F. A. Volpe
Abstract:
The toroidal phase and rotation of otherwise locked magnetic islands of toroidal mode number n=1 are controlled in the DIII-D tokamak by means of applied magnetic perturbations of n=1. Pre-emptive perturbations were applied in feedforward to "catch" the mode as it slowed down and entrain it to the rotating field before complete locking, thus avoiding the associated major confinement degradation. A…
▽ More
The toroidal phase and rotation of otherwise locked magnetic islands of toroidal mode number n=1 are controlled in the DIII-D tokamak by means of applied magnetic perturbations of n=1. Pre-emptive perturbations were applied in feedforward to "catch" the mode as it slowed down and entrain it to the rotating field before complete locking, thus avoiding the associated major confinement degradation. Additionally, for the first time, the phase of the perturbation was optimized in real-time, in feedback with magnetic measurements, in order for the mode's phase to closely match a prescribed phase, as a function of time. Experimental results confirm the capability to hold the mode in a given fixed-phase or to rotate it at up to 20 Hz with good uniformity. The control coil currents utilized in the experiments agree with the requirements estimated by an electromechanical model. Moreover, controlled rotation at 20 Hz was combined with Electron Cyclotron Current Drive (ECCD) modulated at the same frequency. This is simpler than regulating the ECCD modulation in feedback with spontaneous mode rotation, and enables repetitive, reproducible ECCD deposition at or near the island O-point, X-point and locations in between, for careful studies of how this affects the island stability. Current drive was found to be radially misaligned relative to the island, and resulting growth and shrinkage of islands matched expectations of the Modified Rutherford Equation for some discharges presented here. Finally, simulations predict the as designed ITER 3D coils can entrain a small island at sub-10 Hz frequencies.
△ Less
Submitted 15 January, 2018;
originally announced January 2018.
-
Resolving ECRH deposition broadening due to edge turbulence in DIII-D by heat deposition measurement
Authors:
Michael W Brookman,
Matthew B Thomas,
Jarrod Leddy,
C Craig Petty,
Robert J La Haye,
K Barada,
Terry L Rhodes,
Zheng Yan,
Max E Austin,
Roddy G L Vann
Abstract:
Interaction between microwave power, used for local heating and mode control, and density fluctuations can produce a broadening of the injected beam, as confirmed in experiment and simulation. Increased power deposition width could impact suppression of tearing mode structures on ITER. This work discusses the experimental portion of an effort to understand scattering of injected microwaves by turb…
▽ More
Interaction between microwave power, used for local heating and mode control, and density fluctuations can produce a broadening of the injected beam, as confirmed in experiment and simulation. Increased power deposition width could impact suppression of tearing mode structures on ITER. This work discusses the experimental portion of an effort to understand scattering of injected microwaves by turbulence on the DIII-D tokamak. The corresponding theoretical modeling work can be found in M.B. Thomas et. al.: Submitted to Nuclear Fusion (2017)[Author Note - this paper to be published in same journal]. In a set of perturbative heat transport experiments, tokamak edge millimeter-scale fluctuation levels and microwave heat deposition are measured simultaneously. Beam broadening is separated from heat transport through fitting of modulated fluxes. Electron temperature measurements from a 500 kHz, 48-channel ECE radiometer are Fourier analyzed and used to calculate a deposition-dependent flux. Consistency of this flux with a transport model is evaluated. A diffusive and convective transport solution is linearized and compared with energy conservation-derived fluxes. Comparison between these two forms of heat flux is used to evaluate the quality of ECRF deposition profiles, and fitting finds a significant broadening of 1D equilibrium ray tracing calculations from the benchmarked TORAY-GA ray tracing code is needed. The physical basis, cross-validation, and application of the heat flux method is presented. The method is applied to a range of DIII-D discharges and finds a broadening factor of the deposition profile width which scales linearly with edge density fluctuation level. These experimental results are found to be consistent with the full-wave beam broadening measured by the 3D full wave simulations in the same discharges.
△ Less
Submitted 10 October, 2017;
originally announced October 2017.
-
Resolving ECRH deposition broadening due to edge turbulence in DIII-D by 3D full-wave simulations
Authors:
Matthew B Thomas,
Michael W Brookman,
Max E Austin,
Alf Köhn,
Robert J La Haye,
Jarrod B Leddy,
Roddy G L Vann,
Zheng Yan
Abstract:
Edge plasma density fluctuations are shown to have a significant effect on the electron cyclotron resonance heating (ECRH) beam in the DIII-D tokamak. Experimental measurements of the ECRH deposition profile have been taken in three operating scenarios: L-mode, H-mode and negative triangularity. Each scenario corresponds to distinct turbulence characteristics in the edge region through which the b…
▽ More
Edge plasma density fluctuations are shown to have a significant effect on the electron cyclotron resonance heating (ECRH) beam in the DIII-D tokamak. Experimental measurements of the ECRH deposition profile have been taken in three operating scenarios: L-mode, H-mode and negative triangularity. Each scenario corresponds to distinct turbulence characteristics in the edge region through which the beam must propagate. The measured ECRH deposition profile is significantly broadened by comparison to the profile predicted by the ray tracing code TORAY-GA and has been shown to scale with the severity of edge turbulence. Conventional ray tracing does not include the effects of turbulence and therefore a 3D full-wave cold plasma finite difference time domain code EMIT-3D is presented and used for the simulations. The turbulence is generated through the Hermes model in the BOUT++ framework which takes as input the measured time averaged electron density, temperature and magnetic field profiles for the specific shot in question. The simulated turbulence is constrained to match the experimentally measured (by use of the BES and DBS systems) correlation length and normalised fluctuation levels. The predictions of the beam broadening from the simulations are found to agree very well with the experimentally-observed broadening in all cases: L-mode, H-mode and negative triangularity. Due to the large gradients within the H-mode edge, the resolution uncertainty and error in the measurement from Thomson scattering and BES diagnostics result in a spread in the simulated turbulence amplitude. In light of this a parameter scan through the range in experimental diagnostic measurement uncertainty has been conducted to explore the impact on beam broadening predictions.
△ Less
Submitted 9 October, 2017;
originally announced October 2017.
-
Statistical analysis of $m/n$ = 2/1 locked and quasi-stationary modes with rotating precursors at DIII-D
Authors:
R. Sweeney,
W. Choi,
R. J. La Haye,
S. Mao,
K. E. J. Olofsson,
F. A. Volpe,
the DIII-D Team
Abstract:
A database has been developed to study the evolution, the nonlinear effects on equilibria, and the disruptivity of locked and quasi-stationary modes with poloidal and toroidal mode numbers $m=2$ and $n=1$ at DIII-D. The analysis of 22,500 discharges shows that more than 18% of disruptions are due to locked or quasi-stationary modes with rotating precursors (not including born locked modes). A para…
▽ More
A database has been developed to study the evolution, the nonlinear effects on equilibria, and the disruptivity of locked and quasi-stationary modes with poloidal and toroidal mode numbers $m=2$ and $n=1$ at DIII-D. The analysis of 22,500 discharges shows that more than 18% of disruptions are due to locked or quasi-stationary modes with rotating precursors (not including born locked modes). A parameter formulated by the plasma internal inductance $l_i$ divided by the safety factor at 95% of the poloidal flux, $q_{95}$, is found to exhibit predictive capability over whether a locked mode will cause a disruption or not, and does so up to hundreds of milliseconds before the disruption. Within 20 ms of the disruption, the shortest distance between the island separatrix and the unperturbed last closed flux surface, referred to as $d_{edge}$, performs comparably to $l_i/q_{95}$ in its ability to discriminate disruptive locked modes. Out of all parameters considered, $d_{edge}$ also correlates best with the duration of the locked mode. Within 50 ms of a locked mode disruption, average behavior includes exponential growth of the $n=1$ perturbed field, which might be due to the 2/1 locked mode. Even assuming the aforementioned 2/1 growth, disruptivity following a locked mode shows little dependence on island width up to 20 ms before the disruption. Separately, greater deceleration of the rotating precursor is observed when the wall torque is large. At locking, modes are often observed to align with a residual error field. Timescales associated with the mode evolution are also studied and dictate the response times necessary for disruption avoidance and mitigation.
△ Less
Submitted 28 November, 2016; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Avoiding Tokamak disruptions by applying static magnetic fields that align locked modes with stabilizing wave-driven currents
Authors:
F. A. Volpe,
A. Hyatt,
R. J. La Haye,
M. J. Lanctot,
J. Lohr,
R. Prater,
E. J. Strait,
A. Welander
Abstract:
Non-rotating (`locked') magnetic islands often lead to complete losses of confinement in tokamak plasmas, called major disruptions. Here locked islands were suppressed for the first time, by a combination of applied three-dimensional magnetic fields and injected millimetre waves. The applied fields were used to control the phase of locking and so align the island O-point with the region where the…
▽ More
Non-rotating (`locked') magnetic islands often lead to complete losses of confinement in tokamak plasmas, called major disruptions. Here locked islands were suppressed for the first time, by a combination of applied three-dimensional magnetic fields and injected millimetre waves. The applied fields were used to control the phase of locking and so align the island O-point with the region where the injected waves generated non-inductive currents. This resulted in stabilization of the locked island, disruption avoidance, recovery of high confinement and high pressure, in accordance with the expected dependencies upon wave power and relative phase between O-point and driven current.
△ Less
Submitted 29 October, 2015;
originally announced October 2015.
-
Power requirements for electron cyclotron current drive and ion cyclotron resonance heating for sawtooth control in ITER
Authors:
I. T. Chapman,
J. P. Graves,
O. Sauter,
C. Zucca,
O. Asunta,
R. J. Buttery,
S. Coda,
T. Goodman,
V. Igochine,
T. Johnson,
M. Jucker,
R. J. La Haye,
M. Lennholm,
JET-EFDA Contributors
Abstract:
13MW of electron cyclotron current drive (ECCD) power deposited inside the q = 1 surface is likely to reduce the sawtooth period in ITER baseline scenario below the level empirically predicted to trigger neo-classical tearing modes (NTMs). However, since the ECCD control scheme is solely predicated upon changing the local magnetic shear, it is prudent to plan to use a complementary scheme which di…
▽ More
13MW of electron cyclotron current drive (ECCD) power deposited inside the q = 1 surface is likely to reduce the sawtooth period in ITER baseline scenario below the level empirically predicted to trigger neo-classical tearing modes (NTMs). However, since the ECCD control scheme is solely predicated upon changing the local magnetic shear, it is prudent to plan to use a complementary scheme which directly decreases the potential energy of the kink mode in order to reduce the sawtooth period. In the event that the natural sawtooth period is longer than expected, due to enhanced alpha particle stabilisation for instance, this ancillary sawtooth control can be provided from > 10MW of ion cyclotron resonance heating (ICRH) power with a resonance just inside the q = 1 surface. Both ECCD and ICRH control schemes would benefit greatly from active feedback of the deposition with respect to the rational surface. If the q = 1 surface can be maintained closer to the magnetic axis, the efficacy of ECCD and ICRH schemes significantly increases, the negative effect on the fusion gain is reduced, and off-axis negative-ion neutral beam injection (NNBI) can also be considered for sawtooth control. Consequently, schemes to reduce the q = 1 radius are highly desirable, such as early heating to delay the current penetration and, of course, active sawtooth destabilisation to mediate small frequent sawteeth and retain a small q = 1 radius.
△ Less
Submitted 19 June, 2013;
originally announced June 2013.