-
Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach
Authors:
Erol Balkovic,
Joaquim Loizu,
Jonathan P. Graves,
Yi-Min Huang,
Christopher B. Smiet
Abstract:
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson e…
▽ More
We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al. 2020; Loizu & Bonfiglio 2023), we make use of SPEC (Hudson et al. 2012), an equilibrium solver based on the variational principle of the Multi-Region relaxed MHD, featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model $J_\textrm{bs} = C \nabla p$ to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter $Δ'$ and the bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHD code (Huang & Bhattacharjee 2016) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical Modified Rutherford Equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and neoclassical tearing modes, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Equilibrium $β$-limits dependence on bootstrap current in classical stellarators
Authors:
A. Baillod,
J. Loizu,
Z. Qu,
H. P. Arbez,
J. P. Graves
Abstract:
While it is important to design stellarators with high magneto-hydrodynamic (MHD) stability $β$-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of $β$ values. As $β$ increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another ki…
▽ More
While it is important to design stellarators with high magneto-hydrodynamic (MHD) stability $β$-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of $β$ values. As $β$ increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another kind of $β$-limit, the so-called equilibrium $β$-limit. In this paper, we explore numerically the dependence of the equilibrium $β$-limit on the bootstrap current strength using the Stepped Pressure Equilibrium Code (SPEC). We develop a diagnostic to determine whether or not magnetic islands are expected to participate significantly to radial transport, and we build an analytical model to predict the expected equilibrium $β$-limit, which recovers the main features of the numerical results. This research opens the possibility to include additional targets in stellarator optimization functions, provides additional understanding on the existence of magnetic surfaces at large $β$, and is a step forward in the understanding of the equilibrium $β$-limit.
△ Less
Submitted 16 May, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
Stellarator optimization for nested magnetic surfaces at finite $β$ and toroidal current
Authors:
A. Baillod,
J. Loizu,
J. P. Graves,
M. Landreman
Abstract:
Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this pap…
▽ More
Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this paper, we build upon their work to show the first finite-$β$, fixed- and free-boundary optimization of SPEC equilibria for good magnetic surfaces. The objective function is constructed with the Greene's residue of selected rational surfaces and the optimization is driven by the SIMSOPT framework [M. Landreman et al., J. Open Source Software 6, 3525 (2021)]. We show that the size of magnetic islands and the consequent regions occupied by chaotic field lines can be minimized in a classical stellarator geometry by optimizing either the injected toroidal current profile, the shape of a perfectly conducting wall surrounding the plasma (fixed-boundary case), or the coils (free-boundary case), in a reasonable amount of computational time. This work shows that SPEC can be used as an equilibrium code both in a two-step or single-step stellarator optimization loop.
△ Less
Submitted 30 November, 2021;
originally announced November 2021.
-
Verification and validation of gyrokinetic and kinetic-MHD simulations for internal kink instability in DIII-D tokamak
Authors:
G. Brochard,
J. Bao,
C. Liu,
N. Gorelenkov,
G. Choi,
G. Dong,
P. Liu,
J. Mc. Clenaghan,
J. H. Nicolau,
F. Wang,
W. H. Wang,
X. Wei,
W. L. Zhang,
W. Heidbrink,
J. P. Graves,
Z. Lin,
H. Lütjens
Abstract:
Verification and validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kin…
▽ More
Verification and validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kink mode when all kinetic effects are suppresed. The simulated radial mode structures agree quantitatively with the electron cyclotron emission measurement after adjusting, within the experimental uncertainty, the safety factor q=1 flux-surface location in the equilibrium reconstruction. Compressible magnetic perturbations strongly destabilize the kink, while poloidal variations of the equilibrium current density reduce the growth rate of the kink. Furthermore, kinetic effects of thermal ions are found to decrease the kink growth rate in kinetic-MHD simulations, but increase the kink growth rate in gyrokinetic simulations, due to the additional drive of the ion temperature gradient and parallel electric field. Kinetic thermal electrons are found to have negligible effects on the internal kink instability.
△ Less
Submitted 20 September, 2021;
originally announced September 2021.
-
Hybrid guiding-centre/full-orbit simulations in non-axisymmetric magnetic geometry exploiting general criterion for guiding-centre accuracy
Authors:
David Pfefferlé,
Jonathan P. Graves,
Wilfred A. Cooper
Abstract:
To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details a…
▽ More
To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details are provided on how to implement it in cylindrical coordinates, and in flux coordinates that rely on the geometric toroidal angle. A means of switching between guiding-centre and full-orbit equations at first order in Larmor radius with minimal discrepancy is shown. Techniques are applied to a MAST (Mega Amp Spherical Tokamak) helical core equilibrium in which the inner kinked flux-surfaces are tightly compressed against the outer axisymmetric mantle and where the parallel current peaks at the nearly rational surface. This is put in relation with the simpler situation $\vec{B}(x,y,z) = B_0 [\sin(kx) \vec{e_y} + \cos(kx)\vec{e_z}]$, for which full orbits and lowest order drifts are obtained analytically. In the kinked equilibrium, the full orbits of NBI fast ions are solved numerically and shown to follow helical drift surfaces. This result partially explains the off-axis redistribution of NBI fast particles in the presence of MAST Long-Lived Modes (LLM).
△ Less
Submitted 26 February, 2015; v1 submitted 17 December, 2014;
originally announced December 2014.
-
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.