Search

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n ≥ 3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The "double tearing mode" (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature). [ABSTRACT FROM AUTHOR]

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. We demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. We generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values. [ABSTRACT FROM AUTHOR]

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). The simulations show important 3D features for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here. [ABSTRACT FROM AUTHOR]

Recently, it has been shown that a vertical displacement event (VDE) can occur in ITER even when the walls are perfect conductors, as a consequence of the current quench (CQ) [A. H. Boozer, Phys. Plasmas 26, 114501 (2019)]. We used the extended-MHD code M3D-C1 with an ITER-like equilibrium and induced a CQ to explore cold VDEs in the limit of perfectly conducting walls, using different wall geometries. In the case of a rectangular first wall with the side walls far away from the plasma, we obtained very good agreement with the analytical model developed by Boozer that considers a top/bottom flat-plates wall. We show that the solution in which the plasma remains at the initial equilibrium position is improved when bringing the side walls closer to the plasma. When approximating the ITER first wall as a perfect conductor, the plasma remains stable at the initial equilibrium position far beyond the value predicted by the flat-plates wall limit. When considering an opposite limit in which only the inner shell of the ITER vacuum vessel acts as a perfect conductor, the plasma is displaced during the CQ, but the edge safety factor remains above 2 longer in the current decay compared to the flat-plates wall limit. In all the simulated cases, the VDE is found to be strongly dependent on the plasma current, in agreement with a similar finding in the flat-plates wall limit, showing an important difference with hot VDEs in which the CQ is not a necessary condition. [ABSTRACT FROM AUTHOR]

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C1, JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a 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-C1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall. [ABSTRACT FROM AUTHOR]

Employment of non-inductive plasma start-up techniques would considerably simplify the design of a spherical tokamak fusion reactor. Transient coaxial helicity injection (CHI) is a promising method, expected to scale favorably to next-step reactors. However, the implications of reactor-relevant parameters on the initial breakdown phase for CHI have not yet been considered. Here, we evaluate CHI breakdown in reactor-like configurations using an extension of the Townsend avalanche theory. We find that a CHI electrode concept in which the outer vessel wall is biased to achieve breakdown, while previously successful on NSTX and HIT-II, may exhibit a severe weakness when scaled up to a reactor. On the other hand, concepts which employ localized biasing electrodes such as those used in QUEST would avoid this issue. Assuming that breakdown can be successfully attained, we then apply scaling relationships to predict plasma parameters attainable in the transient CHI discharge. Assuming the use of 1 Wb of injector flux, we find that plasma currents of 1 MA should be achievable. Furthermore, these plasmas are expected to Ohmically self-heat with more than 1 MW of power as they decay, facilitating efficient hand-off to steady-state heating sources. These optimistic scalings are supported by Tokamak Simulation Code simulations. [ABSTRACT FROM AUTHOR]

We demonstrate that in a 3D resistive magnetohydrodynamic simulation, for some parameters it is possible to form a stationary state in a tokamak where a saturated interchange mode in the center of the discharge drives a near helical flow pattern that acts to nonlinearly sustain the configuration by adjusting the central loop voltage through a dynamo action. This could explain the physical mechanism for maintaining stationary nonsawtoothing "hybrid" discharges, often referred to as "flux pumping." [ABSTRACT FROM AUTHOR]

The National Spherical Torus Experiment (NSTX) at Princeton University, Princeton, NJ, is designed to assess the potential of the low-aspect-ratio spherical torus concept for magnetic plasma confinement. The plasma has been heated by up to 7 MW of neutral beam injection (NBI) at an injection energy of 100 keV and up to 6 MW of high harmonic fast wave (HHFW) at 30 MHz. NSTX has achieved βT of 32%. A variety of MHD phenomena have been observed to limit β. NSTX has now begun addressing TE scaling, β limits, and current drive issues. During the NBI heating experiments, a broad Ti profile with Ti up to 2 keV, Ti > Te and a large toroidal rotation were observed. Transport analysis suggests that the impurity ions have diffusivities approaching neoclassical. For L-Mode plasmas, TE is up to two times the ITER97L L-Mode scaling and exceeds the ITER98pby2 H-Mode scaling in some cases. Transitions to H-Mode have been observed which result in an approximate doubling of TE after the transition in some conditions. During HHFW heating, Te > Ti and Te up to 3.5 keV were observed. Current drive has been studied using both coaxial helicity injection with up to 390 kA of toroidal current and HHFW. HHFW has produced H-modes with significant bootstrap current fraction at low Ip, high q, and high βp. [ABSTRACT FROM AUTHOR]

Recent operation of the Tokamak Fusion Test Reactor (TFTR) [Plasma Phys. Controlled Nucl. Fusion Research 1, 51 (1986) ] has produced plasma equilibria with values of Γ ≡ β[sub p eq] + l[sub i]/2 as large as 7, ∈β[sub p dia] ≡ 2μ[sub 0] ∈ 〈p[sub ⊥]〉 / 〈〈B[sub p]〉〉² as large as 1.6, and Troyon normalized diamagnetic beta [Plasma Phys. Controlled Fusion 26, 209 ( 1984); Phys. Lett. ll0A, 29 (1985)], β[sub N dia] ≡ 10[sup 8] 〈β[sub t⊥]〉aB[sub o]/I[sub p] as large as 4.7. When ∈β[sub p dia] ... 1.25, a separatrix entered the vacuum chamber, producing a naturally diverted discharge that was sustained for many energy confinement times, τ[sub E]. The largest values of ∈β[sub p] and plasma stored energy were obtained when the plasma current was ramped down prior to neutral beam injection. The measured peak ion and electron temperatures were as large as 24 and 8.5 keV, respectively. Plasma stored energy in excess of 2.5 MJ and τ[sub E] greater than 130 msec were obtained. Confinement times of greater than 3 times that expected from L-mode predictions have been achieved. The fusion power gain Q[sub DD] reached a value of 1.3 × 10[sup -3] in a discharge with I[sub p] = 1 MA and ∈β[sub p dia] = 0.85. A large, sustained negative loop voltage during the steady-state portion of the discharge indicates that a substantial noninductive component of I[sub p] exists in these plasmas. Transport code analysis indicates that the bootstrap current constitutes up to 65% of I[sub p]. Magnetohydrodynamic (MHD) ballooning stability analysis shows that, while these plasmas are near, or at the β[sub p], limit, the pressure gradient in the plasma core is in the first region of stability to high-n modes. [ABSTRACT FROM AUTHOR]