Halo Currents Research Articles

Rotation of tokamak halo currents

During tokamak disruptions, halo currents, which can be tenths of the total plasma current, can flow at the plasma edge along the magnetic field lines that intercept the chamber walls. Non-axisymmetric halo currents are required to maintain force balance as the plasma kinks when the edge safety factor drops to about two in a vertical displacement event. The plasma quickly assumes a definite toroidal velocity va(r) with respect to that of the magnetic kink, vk, where va(r) is set by the radial electric field required for ambipolarity. The plasma velocity, vpl=va+vk, near the edge is influenced by the interaction with neutrals and with the potential in the halo required for quasi-neutrality on open magnetic field lines, and the plasma velocity in the core is influenced by external error fields. When plasma effects dominate magnetic locking, the magnetic kink should rotate at a diamagnetic speed of either the edge or the core. If the magnetic field lines of the halo plasma intercept the wall at locations of very different electrical conductivity, the toroidal rotation of the halo currents can intermittently stall at wall locations of high conductivity. Such stalling is seen in experiments. The toroidal phase difference between the stalled halo currents and the kink, which is expected to rotate smoothly, must satisfy δϕ<±π/2. A concern cited by ITER engineers is that the time varying force of the rotating halo could substantially increase the disruption loads on in-vessel components.

Characterization of disruption halo currents in the National Spherical Torus Experiment

This paper describes the general characteristics of disruption halo currents in the National Spherical Torus Experiment (Ono et al 2000 Nucl. Fusion 40 557). The commonly observed types of vertical motion and resulting halo current patterns are described, and it is shown that plasma discharges developing between components can facilitate halo current flow. The halo current fractions and toroidal peaking factors at various locations in the device are presented. The maximum product of these two metrics for localized halo current measurements is always significantly less than the worst-case expectations from conventional aspect ratio tokamaks (which are typically written in terms of the total halo current). The halo current fraction and impulse is often largest in cases with the fastest plasma current quenches and highest quench rates. The effective duration of the halo current pulse is comparable to or shorter than the plasma current quench time. The largest halo currents have tended to occur in lower β and lower elongation plasmas. The sign of the poloidal halo current is reversed when the toroidal field direction is reversed.

Synergetic Effects of Runaway and Disruption Induced by VDE on the First Wall Damage in HL-2A

The plasma facing component in HL-2A has been damaged seriously after disruption, and for this reason its operation is suspended for maintenance. The experimental phenomena and plasma configurations, calculated by the current filament code (CF-code) using the plasma parameters measured by diagnostics and the signals of the magnetic probes, confirm that the first wall is damaged by the synergetic effects of runaway electrons and disruption induced by a vertical displacement event (VDE). When the plasma column is displaced upward/downward, the strong runaway electrons normally hit the baffle plate of the MP3 or MP1 coil in the upper and lower divertor during the disruption, causing the baffle plates to be holed and wrinkled by the energetic runaway current, and water (for cooling or heating the baffle plates) to leak into the vacuum vessel. Another disastrous consequence is that bellows underlying the baffle plate and outside the coil of MP3 for connecting two segments of the jacket casing pipe are punctured by arcing. The arc may be part of the halo current that forms a complete circuit. The experimental phenomena are indirect but compelling evidence for the existence of a halo current during the disruption and VDE, though the halo current has not been measured by the diagnostics in the HL-2A tokamak.

A solid tungsten divertor for ASDEX Upgrade

The conceptual design of a solid tungsten divertor for ASDEX Upgrade (AUG) is presented. The Div-III design is compatible with the existing divertor structure. It re-establishes the energy and heat receiving capability of a graphite divertor and overcomes the limitations of tungsten coatings. In addition, a solid tungsten divertor allows us to investigate erosion and bulk deuterium retention as well as test castellation and target tilting. The design criteria as well as calculations of forces due to halo and eddy currents are presented. The thermal properties of the proposed sandwich structure are calculated with finite element method models. After extensive testing of a target tile in the high heat flux test facility GLADIS, two solid tungsten tiles were installed in AUG for in-situ testing.

Disruption mitigation by massive gas injection in JET

Disruption mitigation is mandatory for ITER in order to reduce forces, to mitigate heat loads during the thermal quench and to avoid runaway electrons (REs). A fast disruption mitigation valve has been installed at JET to study mitigation by massive gas injection. Different gas species and amounts have been investigated with respect to timescales and mitigation efficiency. We discuss the mitigation of halo currents as well as sideways forces during vertical displacement events, the mitigation of heat loads by increased energy dissipation through radiation, the heat loads which could arise by asymmetric radiation and the suppression of REs.

Disruption avoidance by means of electron cyclotron waves

Disruptions are very challenging to ITER operation as they may cause damage to plasma facing components due to direct plasma heating, forces on structural components due to halo and eddy currents and the production of runaway electrons. Electron cyclotron (EC) waves have been demonstrated as a tool for disruption avoidance by a large set of recent experiments performed in ASDEX Upgrade and FTU using various disruption types, plasma operating scenarios and power deposition locations. The technique is based on the stabilization of magnetohydrodynamic (MHD) modes (mainly m/n = 2/1) through the localized injection of EC power on the resonant surface. This paper presents new results obtained in ASDEX Upgrade regarding stable operation above the Greenwald density achieved after avoidance of density limit disruptions by means of ECRH and suitable density feedback control (L-mode ohmic plasmas, Ip = 0.6 MA, Bt = 2.5 T) and NTM-driven disruptions at high-β limit delayed/avoided by means of both co-current drive (co-ECCD) and pure heating (ECRH) with power ⩽1.7 MW (H-mode NBI-heated plasmas, PNBI ∼ 7.5 MW, Ip = 1 MA, Bt = 2.1 T, q95 ∼ 3.6). The localized perpendicular injection of ECRH/ECCD onto a resonant surface leads to the delay and/or complete avoidance of disruptions. The experiments indicate the existence of a power threshold for mode stabilization to occur. An analysis of the MHD mode evolution using the generalized Rutherford equation coupled to the frequency and phase evolution equations shows that control of the modes is due to EC heating close to the resonant surface. The ECRH contribution (Δ′H term) is larger than the co-ECCD one in the initial and more important phase when the discharge is ‘saved’. Future research and developments of the disruption avoidance technique are also discussed.

Techniques for the measurement of disruption halo currents in the National Spherical Torus Experiment

This paper describes techniques for measuring halo currents, and their associated toroidal peaking, in the National Spherical Torus Experiments [M. Ono et al., Nucl. Fusion 40, 557 (2000)]. The measurements are based on three techniques: (1) measurement of the toroidal field created by the poloidal halo current, either with segmented Rogowski coils or discrete toroidal field sensors, (2) the direct measurement of halo currents into specially instrument tiles, and (3) small Rogowski coils placed on the mechanical supports of in-vessel components. For the segmented Rogowski coils and discrete toroidal field detectors, it is shown that the toroidal peaking factor inferred from the data is significantly less than the peaking factor of the underlying halo current distribution, and a simple model is developed to relate the two. For the array of discrete toroidal field detectors and small Rogowski sensors, the compensation steps that are used to isolate the halo current signal are described. The electrical and mechanical design of compact under-tile resistive shunts and mini-Rogowski coils is described. Example data from the various systems are shown.

Overview of JET results

Since the last IAEA Conference JET has been in operation for one year with a programmatic focus on the qualification of ITER operating scenarios, the consolidation of ITER design choices and preparation for plasma operation with the ITER-like wall presently being installed in JET. Good progress has been achieved, including stationary ELMy H-mode operation at 4.5 MA. The high confinement hybrid scenario has been extended to high triangularity, lower ρ* and to pulse lengths comparable to the resistive time. The steady-state scenario has also been extended to lower ρ* and ν* and optimized to simultaneously achieve, under stationary conditions, ITER-like values of all other relevant normalized parameters. A dedicated helium campaign has allowed key aspects of plasma control and H-mode operation for the ITER non-activated phase to be evaluated. Effective sawtooth control by fast ions has been demonstrated with 3He minority ICRH, a scenario with negligible minority current drive. Edge localized mode (ELM) control studies using external n = 1 and n = 2 perturbation fields have found a resonance effect in ELM frequency for specific q95 values. Complete ELM suppression has, however, not been observed, even with an edge Chirikov parameter larger than 1. Pellet ELM pacing has been demonstrated and the minimum pellet size needed to trigger an ELM has been estimated. For both natural and mitigated ELMs a broadening of the divertor ELM-wetted area with increasing ELM size has been found. In disruption studies with massive gas injection up to 50% of the thermal energy could be radiated before, and 20% during, the thermal quench. Halo currents could be reduced by 60% and, using argon/deuterium and neon/deuterium gas mixtures, runaway electron generation could be avoided. Most objectives of the ITER-like ICRH antenna have been demonstrated; matching with closely packed straps, ELM resilience, scattering matrix arc detection and operation at high power density (6.2 MW m−2) and antenna strap voltages (42 kV). Coupling measurements are in very good agreement with TOPICA modelling.

Contribution of ASDEX Upgrade to disruption studies for ITER

This paper describes the most recent contributions of ASDEX Upgrade to ITER in the field of disruption studies. (1) The ITER specifications for the halo current magnitude are based on data collected from several tokamaks and summarized in the plot of the toroidal peaking factor versus the maximum halo current fraction. Even if the maximum halo current in ASDEX Upgrade reaches 50% of the plasma current, the duration of this maximum lasts a fraction of a ms. (2) Long-lasting asymmetries of the halo current are rare and do not give rise to a large asymmetric component of the mechanical forces on the machine. Differently from JET, these asymmetries are neither locked nor exhibit a stationary harmonic structure. (3) Recent work on disruption prediction has concentrated on the search for a simple function of the most relevant plasma parameters, which is able to discriminate between the safe and pre-disruption phases of a discharge. For this purpose, the disruptions of the last four years have been classified into groups and then discriminant analysis is used to select the most significant variables and to derive the discriminant function. (4) The attainment of the critical density for the collisional suppression of the runaway electrons seems to be technically and physically possible on our medium size tokamak. The CO2 interferometer and the AXUV diagnostic provide information on the highly 3D impurity transport process during the whole plasma quench.

A diffusive model for halo width growth during vertical displacement events

The electromagnetic loads produced by halo currents during vertical displacement events (VDEs) impose stringent requirements on the strength of ITER in-vessel components. A predictive understanding of halo current evolution is essential for ensuring the robust design of these components. A significant factor determining that evolution is the plasma resistance, which is a function of three quantities: the resistivities of the core and halo regions, and the halo region width. A diffusive model of halo width growth during VDEs has been developed, which provides one part of a physics basis for predictive halo current simulations. The diffusive model was motivated by DIII-D observations that VDEs with cold post-thermal quench plasma and a current decay time much faster than the vertical motion (type I VDE) possess much wider halo region widths than warmer plasma VDEs, where the current decay is much slower than the vertical motion (type II). A 2D finite element code is used to model the diffusion of toroidal halo current during selected type I and type II DIII-D VDEs. The model assumes a core plasma region within the last closed flux surface (LCFS) diffusing current into a halo plasma filling the vessel outside the LCFS. LCFS motion and plasma temperature are prescribed from experimental observations. The halo width evolution produced by this model compares favourably with experimental measurements of type I and type II toroidal halo current width evolution.

Using mixed gases for massive gas injection disruption mitigation on Alcator C-Mod

Mixed gases are used for massive gas injection disruption mitigation on Alcator C-Mod in order to optimize radiation efficiency, halo current reduction and response time. Gas mixtures of helium and argon (argon fraction 0–50%) are investigated in detail, as well as mixtures of deuterium, argon, krypton and helium. Experiments show that injecting He/Ar mixtures leads to faster thermal and current quenches than with pure helium or argon injection, thus improving the time response of the disruption mitigation system and reducing the halo current. Small fractions of argon (∼5–10%) in helium also lead to optimized radiation fractions with large electron density increases in the core plasma. These results are consistent with the expectation that small fractions of argon will be entrained with the faster helium in the early phases of gas flow. The gas mixing allows one to simultaneously exploit the fast particle delivery rate of light helium gas and the large radiation capability of argon.

Specification of asymmetric VDE loads of the ITER tokamak

During asymmetric vertical displacement events (AVDEs) associated with the kink mode of the plasma two asymmetry phenomena were observed in existing tokamaks, in particular in JET [1]. The related halo currents flowing in the passive structure were identified as the cause of asymmetric EM loads on tokamak components. The first phenomenon is a toroidal peak of the poloidal halo current that flows in the passive structure. The second phenomenon is that the toroidal plasma current is not uniform toroidally, so a toroidally non-uniform current flows in the vessel [2]. The specification of the expected characteristics of both phenomena as well as of the consequent asymmetric loads in ITER are summarized here. The related loads are specified for likely, unlikely and extremely unlikely AVDEs.

The halo current in ASDEX Upgrade

Due to the complexity of the phenomena involved, a self-consistent physical model for the prediction of the halo current is not available. Therefore the ITER specifications of the spatial distribution and evolution of the halo current rely on empirical assumptions. This paper presents the results of an extensive analysis of the halo current measured in ASDEX Upgrade with particular emphasis on the evolution of the halo region, on the magnitude and time history of the halo current, and on the structure and duration of its toroidal and poloidal asymmetries. The effective length of the poloidal path of the halo current in the vessel is found to be rather insensitive to plasma parameters. Large values of the toroidally averaged halo current are observed in both vertical displacement events and centred disruptions but last a small fraction of the current quench; they coincide typically with a large but short-lived MHD event.

Evaluation of electromagnetic loads on various design options of the ITER diagnostic upper port plug during plasma disruptions

Electromagnetic (EM) loads due to eddy current and halo current during plasma disruptions are evaluated for the ITER diagnostic upper port plug. To reduce strong EM loads acting on the port plug fixed to the vacuum vessel like a cantilever beam, three design options have been considered: removal of the diagnostic first wall, slitting of the diagnostic shield module and recess of the port plug. The main focus of the present study is to examine the efficacy of these options in terms of EM loads on the upper port plug. It is found that making slits is more effective than removing the first wall. It is also shown that the upper port plug needs to be recessed to reduce the EM load induced by halo current.

Structure of the controlled halo-current magnetic field in the T-10 tokamak

In the T-10 tokamak, the magnetic field spatially resonant with a helical MHD perturbation is generated using the controlled halo current supplied using a contact method in the scrape-off-layer plasma. This paper is concerned with studying the spatial structure of the halo current and its magnetic field. For this purpose, the magnetic field of the halo current was measured in one of the cross sections of the torus near the tokamak vacuum vessel wall. The spatial distribution of the magnetic field as a function of the halo current configuration was calculated in the cylindrical approximation. The terms proportional to the plasma pressure were disregarded. The configuration of the halo current and the spatial structure of its magnetic field were determined by comparing the calculated and experimental results.

JET disruption studies in support of ITER

Plasma disruptions affect plasma-facing and structural components of tokamaks due to electromechanical forces, thermal loads and generation of high energy runaway electrons (REs). Asymmetries in poloidal halo and toroidal plasma current can now be routinely measured in four positions 90° apart. Their assessment is used to validate the design of the ITER vessel support system and its in-vessel components. The challenge of disruption thermal loads comes from both the short duration over which a large energy has to be lost and the potential for asymmetries. The focus of this paper will be on localized heat loads. Resonant magnetic perturbations failed to reduce the generation of REs in JET. An explanation of the limitations applying to these attempts is offered together with a minimum guideline. The REs generated by a moderate, but fast, Ar injection in limiter plasmas show evidence of milder and more efficient losses due to the high Ar background density.

Plasma–surface interactions during tokamak disruptions and rapid shutdowns

Recent progress in understanding of disruptions and in developing methods to avoid disruption damage is presented. Nearly complete mitigation of conducted heat loads has been achieved with high-Z gas jet shutdown. The resulting local radiation heat flash melting in the main chamber might be a concern in ITER, especially with beryllium walls. During the current quench, significant vessel forces can occur due to halo currents Ihalo; however, these are found to fall reliably below a boundary of (halo current fraction times halo current peaking factor) <0.7 both experimentally and numerically. Numerical simulations indicate that runaway electrons (REs) could cause serious damage to hard-to reach components in ITER, making their suppression a high priority. During the current quench, less than 20% of the density required for collisional suppression of REs appears to have been achieved. Collisional suppression of REs may have been achieved, however, in full-current RE beams with gas injection.

Eddy current induced electromagnetic loads on shield blankets during plasma disruptions in ITER: A benchmark exercise

According to recent updates of ITER shield blanket design, electromagnetic loads during the plasma disruption are being evaluated to verify the mechanical confidence and reliability. As a course of such evaluations, a benchmark activity for the electromagnetic analysis, coordinated by ITER Organization, is underway between ITER parties to compare the calculation results for disruption loads on the blankets. In this paper, we present calculation results for the electromagnetic loads on the simplified but practical model of ITER shield blankets with respect to six representative disruption scenarios of which ITER distributes simulation results based on the DINA code as a reference of the design and analysis. Commercial finite element method software, ANSYS/Emag™, was employed to evaluate the eddy current on the blanket modules with the 40° sector model for major conducting structure of the tokamak including double-walled vacuum vessel, triangular support, and vertical targets of divertors. An interface between ANSYS/Emag™ and plasma simulator was implemented with a conversion tool assigning the plasma current density on the ANSYS elements corresponding to the current filaments in DINA outputs. Discussions are made of the possible improvement of the blanket model taking more realistic blanket configuration into account at the cost of the moderate increase in computational time. A final remark is given of the possibility of incorporating halo currents into ANSYS disruption simulations, which are major sources of electromagnetic loads on in-vessel components including blankets.

Disruption-Mitigation-Technology Concepts and Implications for ITER

Disruptions on ITER present challenges to handle the intense heat flux, the large forces from halo currents, and the potential first wall damage from energetic runaway electrons. Injecting large quantities of material into the plasma during the disruption can reduce the plasma energy and increase its resistivity to mitigate these effects. Assessments of the amount of various mixtures and quantities of the material required have been made to provide collision mitigation of runaway-electron conversion, which is the most difficult challenge. The quantities of the material required (~0.5 MPa·m3 for deuterium or helium gas) are large enough to have implications on the design and operation of the vacuum system and tokamak exhaust processing system.

Overview of the Alcator C-Mod Research Program

This paper summarizes highlights of research results from the Alcator C-Mod tokamak covering the period 2006–2008. Active flow drive, using mode converted ion cyclotron waves, has been observed for the first time in a tokamak plasma, using a mix of D and 3He ion species; toroidal and poloidal flows are driven near the location of the mode conversion layer. ICRF induced edge sheaths are implicated in both the erosion of thin boron coatings and the generation of metallic impurities. Lower hybrid range of frequencies (LHRF) microwaves have been used for efficient current drive, current profile modification and toroidal flow drive. In addition, LHRF has been used to modify the H-mode pedestal, increasing temperature, decreasing density and lowering the pedestal collisionality. Studies of hydrogen isotope retention in solid metallic plasma facing components reveal significantly higher retention than expected from ex situ laboratory studies; a model to explain the results, based on plasma/neutral induced lattice damage, has been developed and tested. During gas-puff mitigation of disruptions, induced MHD instabilities cause the magnetic field to become stochastic, resulting in reduction of halo currents, spreading of plasma power loading and loss of runaway electrons before they cause damage. Detailed pedestal rotation profile measurements have been used to infer Er profiles, and correlation with global H-mode confinement. An improved L-mode regime, obtained at q95 ⩽ 3 with ion drift away from the active X-point, shows very good energy confinement with a strong temperature pedestal, a weak density pedestal, and no evidence of particle or impurity accumulation, without the need for ELMs or any additional edge density regulation mechanism.