The maximum normalized beta achieved in long-pulse tokamak discharges at low collisionality falls significantly below both that observed in short pulse discharges and that predicted by the ideal MHD theory. Recent long-pulse experiments, in particular those simulating the International Thermonuclear Experimental Reactor ͑ITER͒ ͓M. Rosenbluth et al., Plasma Physics and Controlled Nuclear Fusion ͑International Atomic Energy Agency, Vienna, 1995͒, Vol. 2, p. 517͔ scenarios with low collisionality e * , are often limited by low-m/n nonideal magnetohydrodynamic ͑MHD͒ modes. The effect of saturated MHD modes is a reduction of the confinement time by 10%-20%, depending on the island size and location, and can lead to a disruption. Recent theories on neoclassical destabilization of tearing modes, including the effects of a perturbed helical bootstrap current, are successful in explaining the qualitative behavior of the resistive modes and recent results are consistent with the size of the saturated islands. Also, a strong correlation is observed between the onset of these low-m/n modes with sawteeth, edge localized modes ͑ELM͒, or fishbone events, consistent with the seed island required by the theory. We will focus on a quantitative comparison between both the conventional resistive and neoclassical theories, and the experimental results of several machines, which have all observed these low-m/n nonideal modes. This enables us to single out the key issues in projecting the long-pulse beta limits of ITER-size tokamaks and also to discuss possible plasma control methods that can increase the soft ␤ limit, decrease the seed perturbations, and/or diminish the effects on confinement.

No abstract

Demonstrating improved confinement of energetic ions is one of the key goals of the Wendelstein 7-X (W7-X) stellarator. In the past campaigns, measuring confined fast ions has proven to be challenging. Future deuterium campaigns would open up the option of using fusion-produced neutrons to indirectly observe confined fast ions. There are two neutron populations: 2.45 MeV neutrons from thermonuclear and beam-target fusion, and 14.1 MeV neutrons from DT reactions between tritium fusion products and bulk deuterium. The 14.1 MeV neutron signal can be measured using a scintillating fiber neutron detector, whereas the overall neutron rate is monitored by common radiation safety detectors, for instance fission chambers. The fusion rates are dependent on the slowing-down distribution of the deuterium and tritium ions, which in turn depend on the magnetic configuration via fast ion orbits. In this work, we investigate the effect of magnetic configuration on neutron production rates in W7-X. The neutral beam injection, beam and triton slowing-down distributions, and the fusion reactivity are simulated with the ASCOT suite of codes. The results indicate that the magnetic configuration has only a small effect on the production of 2.45 MeV neutrons from DD fusion and, particularly, on the 14.1 MeV neutron production rates. Despite triton losses of up to 50 %, the amount of 14.1 MeV neutrons produced might be sufficient for a time-resolved detection using a scintillating fiber detector, although only in high-performance discharges.

After completing the main construction phase of Wendelstein 7-X (W7-X) and successfully commissioning the device, first plasma operation started at the end of 2015. Integral commissioning of plasma start-up and operation using electron cyclotron resonance heating (ECRH) and an extensive set of plasma diagnostics have been completed, allowing initial physics studies during the first operational campaign. Both in helium and hydrogen, plasma breakdown was easily achieved. Gaining experience with plasma vessel conditioning, discharge lengths could be extended gradually. Eventually, discharges lasted up to 6 s, reaching an injected energy of 4 MJ, which is twice the limit originally agreed for the limiter configuration employed during the first operational campaign. At power levels of 4 MW central electron densities reached 3  ×  1019 m−3, central electron temperatures reached values of 7 keV and ion temperatures reached just above 2 keV. Important physics studies during this first operational phase include a first assessment of power balance and energy confinement, ECRH power deposition experiments, 2nd harmonic O-mode ECRH using multi-pass absorption, and current drive experiments using electron cyclotron current drive. As in many plasma discharges the electron temperature exceeds the ion temperature significantly, these plasmas are governed by core electron root confinement showing a strong positive electric field in the plasma centre.

Compact optimized stellarators offer novel solutions for confining high-β plasmas and developing magnetic confinement fusion. The three-dimensional plasma shape can be designed to enhance the magnetohydrodynamic (MHD) stability without feedback or nearby conducting structures and provide driftorbit confinement similar to tokamaks. These configurations offer the possibility of combining the steady-state low-recirculating power, external control, and disruption resilience of previous stellarators with the low aspect ratio, high β limit, and good confinement of advanced tokamaks. Quasiaxisymmetric equilibria have been developed for the proposed National Compact Stellarator Experiment (NCSX) with average aspect ratio 4-4.4 and average elongation ∼1.8. Even with bootstrap-current consistent profiles, they are passively stable to the ballooning, kink, vertical, Mercier, and neoclassicaltearing modes for β > 4%, without the need for external feedback or conducting walls. The bootstrap current generates only 1/4 of the magnetic rotational transform at β = 4% (the rest is from the coils); thus the equilibrium is much less non-linear and is more controllable than similar advanced tokamaks. The enhanced stability is a result of 'reversed' global shear, the spatial distribution of local shear, and the large fraction of externally generated transform. Transport simulations show adequate fast-ion confinement and thermal neoclassical transport similar to equivalent tokamaks. Modular coils have been designed which reproduce the physics properties, provide good flux surfaces, and allow flexible variation of the plasma shape to control the predicted MHD stability and transport properties.

A three-dimensional MHD equilibrium code is described that does not assume the existence of good flux surfaces. Given an initial guess for the magnetic field, the code proceeds by calculating the pressure-driven current and then by updating the field using Ampere's law. The numerical algorithm to solve the magnetic differential equation for the pressure-driven current is described, and demonstrated for model fields having islands and stochastic regions. The numerical algorithm which solves Ampere's law In three dimensions is also described. Finally, the convergence of the code is illustrated for a particular steilarator equilibrium with no large islands."Presented as an invited talk at the 8th Europhysics Conference on Computational Physics, Elbsee, Fed. Rep. of Germany, May 13-16, 1986. To appear in a special issue of Computer Physics Communications. In this paper we describe a three-dimensional MHD equilibrium code that does not assume the existence of good flux surfaces. The code uses an algorithm for solving the equilibrium equations that is different from the energy-minimizing algorithms which have been developed extensively and successfully for 3D equilibrium codes that assume good surfaces [1][2][3]. We describe the algorithm, our reasons for choosing It, and the numerical methods we have developed to implement it. We also describe the numerical issues raised by the presence of Islands and stochastic regions, and how we have dealt with those issues. In particular, we describe how the code solves magnetic differential equations, using as examples model fields with islands and stochastic regions. Finally, we describe the convergence of the code for a particular stellarator equilibrium with no large inlands.We specialize In this paper to the case of zero net current within each flux surface, as is appropriate for applications to stellarators. Our code has been written to deal with the more general case of arbitrary net currents, except for some small modifications that need to be made to deal with the currents in the islands and stochastic regions.The applications we have in mind for the code are the study of the breaking of flux surfaces in stellarators, and the study of The magnetic field becomes increasingly distorted, it is not known under what conditions the field produced by the pressure-driven currents will cause the formation of large islands and stochastic regions. In particular, it is important to know whether there is an equilibrium beta limit due to such an effect, and what that beta limit is. The convention In evaluating stellarator designs has been to estimate the equilibrium beta limit to be the value of beta at which the magnetic axis shifts halfway to the wall. An analytical calculation of equilibrium island formation in one proposed helical axis stellarator found a much lower beta limit [4].Nominally axisymmetric devices, such as tokamaks, are in practice generally nonaxisymmetrlc because of the presence of tearing instabilities, magnetic ripple from discrete coils, and field error...

The MFBE procedure developed by Strumberger (1997 Nucl. Fusion 37 19) is used to provide an improved starting point for free boundary equilibrium computations in the case of W7-X (Nührenberg and Zille 1986 Phys. Lett. A 114 129) using the Princeton iterative equilibrium solver (PIES) code (Reiman and Greenside 1986 Comput. Phys. Commun. 43 157). Transferring the consistent field found by the variational moments equilibrium code (VMEC) (Hirshmann and Whitson 1983 Phys. Fluids 26 3553) to an extended coordinate system using the VMORPH code, a safe margin between plasma boundary and PIES domain is established. The new EXTENDER P code implements a generalization of the virtual casing principle, which allows field extension both for VMEC and PIES equilibria. This facilitates analysis of the 5/5 islands of the W7-X standard case without including them in the original PIES computation.

With the installation of non-axisymmetric coil systems on major tokamaks for the purpose of studying the prospects of ELM-free operation, understanding the plasma response to the applied fields is a crucial issue. Application of different response models, using standard tools, to DIII-D discharges with applied non-axisymmetric fields from internal coils, is shown to yield qualitatively different results. The plasma response can be treated as an initial value problem, following the system dynamically from an initial unperturbed state, or from a nearby perturbed equilibrium approach, and using both linear and nonlinear models [A. D. Turnbull, Nucl. Fusion 52, 054016 (2012)]. Criteria are discussed under which each of the approaches can yield a valid response. In the DIII-D cases studied, these criteria show a breakdown in the linear theory despite the small 10−3 relative magnitude of the applied magnetic field perturbations in this case. For nonlinear dynamical evolution simulations to reach a saturated nonlinear steady state, appropriate damping mechanisms need to be provided for each normal mode comprising the response. Other issues arise in the technical construction of perturbed flux surfaces from a displacement and from the presence of near nullspace normal modes. For the nearby equilibrium approach, in the absence of a full 3D equilibrium reconstruction with a controlled comparison, constraints relating the 2D system profiles to the final profiles in the 3D system also need to be imposed to assure accessibility. The magnetic helicity profile has been proposed as an appropriate input to a 3D equilibrium calculation and tests of this show the anticipated qualitative behavior.