Search

In preparation for the upcoming deuterium-tritium campaign at the JET tokamak, the ex-vessel fission chamber neutron diagnostics and the neutron activation system will be calibrated in absolute terms at 14 MeV neutron energy, to a required accuracy of less than 10%. Two nominally identical DT neutron generators were chosen as the calibration sources, both of which were fully calibrated and characterized at the UK's National Physical Laboratory. The neutron activation method was adopted as a complementary method for the purpose of determining the absolute value of the neutron yield from the neutron generators and to provide a means of cross check for the active detection methods being employed. The work being presented here shows the derivation of the neutron emission rate from the neutron generators based upon experimental activation foil measurements. © 2017 Elsevier B.V.

The first divertor was installed in the JET machine between 1992 and 1994 and was operated with carbon tiles and then beryllium tiles in 1994-5. Post-mortem studies after these first experiments demonstrated that most of the impurities deposited in the divertor originate in the main chamber, and that asymmetric deposition patterns generally favouring the inner divertor region result from drift in the scrape-off layer. A new monolithic divertor structure was installed in 1996 which produced heavy deposition at shadowed areas in the inner divertor corner, which is where the majority of the tritium was trapped by co-deposition during the deuterium-tritium experiment in 1997. Different divertor geometries have been tested since such as the Gas-Box and High-Delta divertors; a principle objective has been to predict plasma behaviour, transport and tritium retention in ITER. Transport modelling experiments were carried out at the end of four campaigns by puffing C-13-labelled methane, and a range of diagnostics such as quartz-microbalance and rotating collectors have been installed to add time resolution to the post-mortem analyses. The study of material migration after D-D and D-T campaigns clearly revealed important consequences of fuel retention in the presence of carbon walls. They gave a strong impulse to make a fundamental change of wall materials. In 2010 the carbon divertor and wall tiles were removed and replaced with tiles with Be or W surfaces for the ITER-Like Wall Project. For complete list of authors see http://dx.doi.org/10.1016/j.fusengdes.2018.10.002

In this paper we validate the finite element code RAPLICASOL, which models radiofrequency wave propagation in edge plasmas near ICRF antennas, against calculations with the TOPICA code. We compare the output of both codes for the ASDEX Upgrade 2-strap antenna, and for a 4-strap WEST-like antenna. Although RAPLICASOL requires considerably fewer computational resources than TOPICA, we find that the predicted quantities of experimental interest (including reflection coefficients, coupling resistances, S- and Z-matrix entries, optimal matching settings, and even radiofrequency electric fields) are in good agreement provided we are careful to use the same geometry in both codes.

The experimental method developed at ASDEX Upgrade for the determination of the intrinsic tungsten (W) density profile coupling data from the soft X-ray (SXR) diagnostic and vacuum-ultra-violet (VUV) spectroscopy has been upgraded for application to JET plasmas. The strong poloidal asymmetries in the SXR emission are modeled assuming a ln(ϵ(ρ,R)/ϵ(ρ,R0))=λ(ρ)(R2−R02) distribution, where ρ is the flux coordinate, R is the major radius, and λ is the fit parameter. The W density is calculated from the resulting 2D SXR emissivity maps accounting for contributions from a low-Z impurity (typically beryllium) and main ion with the assumption that their contributions are poloidally symmetric. Comparing the result with the independent W concentration measurement of VUV spectroscopy, a recalibration factor for the SXR emissivity is calculated making the method robust against the decrease in the sensitivity of the SXR diodes which has been observed across multiple campaigns. The final 2D W density map is checked for consistency versus the time-evolution of the W concentration measurement from VUV spectroscopy, toroidal rotation measurements from charge exchange recombination spectroscopy, and tomographic reconstructions of bolometry data. The method has been found to be robust for W concentrations above a few 10−5 and in cases where the contributions from other medium-Z impurities such as Ni are negligible.The experimental method developed at ASDEX Upgrade for the determination of the intrinsic tungsten (W) density profile coupling data from the soft X-ray (SXR) diagnostic and vacuum-ultra-violet (VUV) spectroscopy has been upgraded for application to JET plasmas. The strong poloidal asymmetries in the SXR emission are modeled assuming a ln(ϵ(ρ,R)/ϵ(ρ,R0))=λ(ρ)(R2−R02) distribution, where ρ is the flux coordinate, R is the major radius, and λ is the fit parameter. The W density is calculated from the resulting 2D SXR emissivity maps accounting for contributions from a low-Z impurity (typically beryllium) and main ion with the assumption that their contributions are poloidally symmetric. Comparing the result with the independent W concentration measurement of VUV spectroscopy, a recalibration factor for the SXR emissivity is calculated making the method robust against the decrease in the sensitivity of the SXR diodes which has been observed across multiple campaigns. The final 2D W density map is checked for...

The x-ray transitions in Cu- and Ni-like tungsten ions in the 5.19–5.26 A wavelength range that are relevant ˚ as a high-temperature tokamak diagnostic, in particular for JET in the ITER-like wall configuration, have been studied. Tungsten spectra were measured at the upgraded Shanghai- Electron Beam Ion Trap operated with electron-beam energies from 3.16 to 4.55 keV. High-resolution measurements were performed by means of a flat Si 111 crystal spectrometer equipped by a CCD camera. The experimental wavelengths were determined with an accuracy of 0.3–0.4 mA. The wavelength of the ground-state transition in Cu-like tungsten from the ˚ 3p53d104s4d [(3/2,(1/2,5/2)2]1/2 level was measured. All measured wavelengths were compared with those measured from JET ITER-like wall plasmas and with other experiments and various theoretical predictions including COWAN, RELAC, multiconfigurational Dirac-Fock (MCDF), and FAC calculations. To obtain a higher accuracy from theoretical predictions, the MCDF calculations were extended by taking into account correlation effects (configuration-interaction approach). It was found that such an extension brings the calculations closer to the experimental values in comparison with other calculations. National Magnetic Confinement Fusion Program of China 2015GB117000 National Natural Science Foundation of China 11374061 EURATOM 633053

Measurements of neutron streaming through penetrations in biological shields are being carried out at JET fusion device by means of thermoluminescence (TL) detectors with the objective to validate the neutronics codes and nuclear data widely applied in ITER nuclear analyses in a real fusion environment. TLDs response due to the neutron component of the radiation field is related to the neutron fluence in a well-defined neutron energy spectrum. Therefore, a TLDs calibration in real fusion radiation fields is necessary to allow neutron fluence from TL measurements at JET to be more precisely calculated. Hence, a MCP-N and MCP-7 TLDs produced at the IFJ PAN in Kraków were calibrated at the ENEA facilities of Frascati and Casaccia laboratories. The obtained results have been analysed and new calibration factors are proposed. The detection system based on TLDs developed and calibrated for JET experiments can then be generally applied not only to fusion neutron fields, but also to ITER to monitor the neutron fluence outside the biological shield. © 2018

Fast ions have recently been found to significantly impact and partially suppress plasma turbulence both in experimental and numerical studies in a number of scenarios. Understanding the underlying physics and identifying the range of their beneficial effect is an essential task for future fusion reactors, where highly energetic ions are generated through fusion reactions and external heating schemes. However, in many of the gyrokinetic codes fast ions are, for simplicity, treated as equivalent-Maxwellian-distributed particle species, although it is well known that to rigorously model highly non-thermalised particles, a non-Maxwellian background distribution function is needed. To study the impact of this assumption, the gyrokinetic code GENE has recently been extended to support arbitrary background distribution functions which might be either analytic, e.g. slowing down and bi-Maxwellian, or obtained from numerical fast ion models. A particular JET plasma with strong fast-ion related turbulence suppression is revised with these new code capabilities both with linear and nonlinear gyrokinetic simulations. It appears that the fast ion stabilization tends to be less strong but still substantial with more realistic distributions, and this improves the quantitative power balance agreement with experiments., Submitted to Physics of Plasmas

In order to contribute to the understanding of runaway electron generation mechanisms during tokamak disruptions, a test particle tracker is introduced in the JOREK 3D non-linear MHD code, able to compute both full and guiding center relativistic orbits. Tests of the module show good conservation of the invariants of motion and consistency between full orbit and guiding center solutions. A first application is presented where test electron confinement properties are investigated in a massive gas injection-triggered disruption simulation in JET-like geometry. It is found that electron populations initialised before the thermal quench (TQ) are typically not fully deconfined in spite of the global stochasticity of the magnetic field during the TQ. The fraction of ‘survivors’ decreases from a few tens down to a few tenths of percent as the electron energy varies from 1 keV to 10 MeV. The underlying mechanism for electron ‘survival’ is the prompt reformation of closed magnetic surfaces at the plasma core and, to a smaller extent, the subsequent reappearance of a magnetic surface at the edge. It is also found that electrons are less deconfined at 10 MeV than at 1 MeV, which appears consistent with a phase averaging effect due to orbit shifts at high energy.

Heavy impurities, such as tungsten (W), can exhibit strongly poloidally asymmetric density profiles in rotating or radio frequency heated plasmas. In the metallic environment of JET, the poloidal asymmetry of tungsten enhances its neoclassical transport up to an order of magnitude, so that neoclassical convection dominates over turbulent transport in the core. Accounting for asymmetries in neoclassical transport is hence necessary in the integrated modeling framework. The neoclassical drift kinetic code, NEO [E. Belli and J. Candy, Plasma Phys. Controlled Fusion P50, 095010 (2008)], includes the impact of poloidal asymmetries on W transport. However, the computational cost required to run NEO slows down significantly integrated modeling. A previous analytical formulation to describe heavy impurity neoclassical transport in the presence of poloidal asymmetries in specific collisional regimes [C. Angioni and P. Helander, Plasma Phys. Controlled Fusion 56, 124001 (2014)] is compared in this work to numerical results from NEO. Within the domain of validity of the formula, the factor for reducing the temperature screening due to poloidal asymmetries had to be empirically adjusted. After adjustment, the modified formula can reproduce NEO results outside of its definition domain, with some limitations: When main ions are in the banana regime, the formula reproduces NEO results whatever the collisionality regime of impurities, provided that the poloidal asymmetry is not too large. However, for very strong poloidal asymmetries, agreement requires impurities in the Pfirsch-Schlüter regime. Within the JETTO integrated transport code, the analytical formula combined with the poloidally symmetric neoclassical code NCLASS [W. A. Houlberg et al., Phys. Plasmas 4, 3230 (1997)] predicts the same tungsten profile as NEO in certain cases, while saving a factor of one thousand in computer time, which can be useful in scoping studies. The parametric dependencies of the temperature screening reduction due to poloidal asymmetries would need to be better characterised for this faster model to be extended to a more general applicability.

In many tokamaks ferromagnetic material, usually referred to as an iron-core, is present in order to improve the magnetic coupling between the solenoid and the plasma. The presence of the iron core ...

This paper further analyses the EDGE2D-EIRENE simulations presented by [Chankin et al 2017 Nucl. Mater. Energy 12 273], of L-mode JET plasmas in vertical-vertical (VV) and vertical-horizontal (VH) divertor configurations. As expected, the simulated outer divertor ionisation source peaks near the separatrix in VV and radially further out in VH. We identify the reflections of recycled neutrals from lower divertor tiles as the primary mechanism by which ionisation is concentrated on the outer divertor separatrix in the VV configuration. These lower tile reflection pathways (of neutrals from the outer divertor, and to an even greater extent from the inner divertor) dominate the outer divertor separatrix ionisation. In contrast, the lower-tile-reflection pathways are much weaker in the VH simulation and its outer divertor ionisation is dominated by neutrals which do not reflect from any surfaces. Interestingly, these differences in neutral pathways give rise to strong differences in the heat flux density width λq at the outer divertor entrance: λq = 3.2 mm in VH compared to λq = 11.8 mm in VV. In VH, a narrow channel exists in the near scrape-off-layer (SOL) where the convected heat flux, driven by strong Er × B flow and thermoelectric current, dominates over the conducted heat flux. The width of this channel sets λq and is determined by the radial distance between the separatrix and the ionisation peak in the outer divertor EURATOM 633053 RCUK Energy Programme P012450/1

In the preparation for the Deuterium-Tritium campaign, JET will operate with a tritium plasma. The T + T reaction consists of two notable channels: (1) T + T -> He-4 + 2n, (2) T + T -> He-5 + ...

ASCOT Fusion Source Integrator AFSI, an efficient tool for calculating fusion reaction rates and characterizing the fusion products, based on arbitrary reactant distributions, has been developed and is reported in this paper. Calculation of reactor-relevant D-D, D-T and D-3He fusion reactions has been implemented based on the Bosch-Hale fusion cross sections. The reactions can be calculated between arbitrary particle populations, including Maxwellian thermal particles and minority energetic particles. Reaction rate profiles, energy spectra and full 4D phase space distributions can be calculated for the non-isotropic reaction products. The code is especially suitable for integrated modelling in self-consistent plasma physics simulations as well as in the Serpent neutronics calculation chain. Validation of the model has been performed for neutron measurements at the JET tokamak and the code has been applied to predictive simulations in ITER.

The impact of the change of the mass of hydrogen isotopes on the turbulent particle flux is studied. The trapped electron component of the turbulent particle convection induced by collisionality, which is outward in ion temperature gradient turbulence, increases with decreasing thermal velocity of the isotope. Thereby, the lighter is the isotope, the stronger is the turbulent pinch, and the larger is the predicted density gradient at the null of the particle flux. The passing particle component of the flux increases with decreasing mass of the isotope and can also affect the predicted density gradient. This effect is however subdominant for usual core plasma parameters. The analytical results are confirmed by means of both quasi-linear and nonlinear gyrokinetic simulations, and an estimate of the difference in local density gradient produced by this effect as a function of collisionality has been obtained for typical plasma parameters at mid-radius. Analysis of currently available experimental data from the JET and the ASDEX Upgrade tokamaks does not show any clear and general evidence of inconsistency with this theoretically predicted effect outside the errorbars and also allows the identification of cases providing weak evidence of qualitative consistency.

International audience; The paper is dedicated to the characterisation of multi-strap ICRH antenna coupling to plasma. Relevance of traditional concept of coupling resistance to antennas with mutually coupled straps is revised and the importance of antenna port excitation consistency for application of the concept is highlighted. A method of antenna S-matrix measurement in presence of plasma is discussed allowing deeper insight into the problem of antenna-plasma coupling. The method is based entirely on the RF plant hardware and control facilities available at JET and it involves application of variable phasing between the antenna straps during the RF plant operations at >100 kW. Unlike traditional techniques relying on low-power (~10 mW) network analysers, the applied antenna voltage amplitudes are relevant to practical conditions of ICRH operations; crucially, they are high enough to minimise possible effects of antenna loading non-linearity due to the RF sheath effects and other phenomena which could affect low-power measurements. The method has been successfully applied at JET to conventional 4-port ICRH antennas energised at frequencies of 33 MHz, 42 MHz and 51 MHz during L-mode plasma discharges while different gas injection modules (GIMs) were used to maintain comparable plasma densities during the pulses. The S-matrix assessment and its subsequent processing yielding `global' antenna coupling resistances in conditions of equalised port maximum voltages allowed consistent description of antenna coupling to plasma at different strap phasing, operational frequencies and applied GIMs. Comprehensive experimental characterisation of mutually coupled antenna straps in presence of plasma also provided a unique opportunity for in-depth verification of TOPICA computer simulations.

We report the first experimental observation of the depolarization of the Thomson scattering (TS) radiation, a relativistic effect expected to occur in very high T e plasmas and never observed so far in a fusion machine. A set of unused optical fibers in the collection optics of the high resolution Thomson scattering system of JET has been used to detect the depolarized TS radiation during a JET campaign with T e 8keV . A linear polarizer with the axis perpendicular to the direction of the incident E -field was placed in front of a fiber optic pair observing a region close to the plasma core, while another fiber pair with no polariser simultaneously observed an adjacent plasma region. The measured intensity ratio was found to be consistent with the theory, taking into account sensitivity coefficients of the two measurement channels determined with post-experiment calibrations and Raman scattering. This depolarization effect is at the basis of polarimetric TS, a different and complementary method for the analysis of TS spectra that can provide significant advantages for T e measurements in very hot plasmas such as in ITER ( T e 40keV ).

Experimental observations in JET tokamak plasmas and gyrokinetic simulations point to an important role, for electron heat transport, of electron-scale instabilities and of their interaction with ion-scale instabilities. Since these effects are maximized for strong electron heating and ion-scale modes close to marginal stability, these findings are of high relevance for ITER plasmas, featuring both conditions. Gyrokinetic and quasi-linear transport models accounting for multi-scale effects are assessed against JET experimental results.

Neutron emission spectroscopy is a diagnostic technique that allows for energy measurements of neutrons born in nuclear reactions. The JET tokamak fusion experiment (Culham, UK) has a special role in this respect as advanced spectrometers for 2.5 MeV and 14 MeV neutrons have been developed here for the first time for measurements of the neutron emission spectrum from D and DT plasmas with unprecedented accuracy. Twin liquid scintillating neutron spectrometers were built and calibrated at the Physikalisch-Technische Bundesanstalt (PTB) (Braunschweig, Germany) and installed on JET in the recent years with tangential-equatorial (KM12) and vertical-radial (KM13) view lines, with the latter only recently operational. This article reports on the performance of KM12 and on the development of the data analysis methods in order to extract physics information upon D ions kinematics in JET auxiliary-heated D plasmas from 2.5 MeV neutron measurements. The comparison of these results with the correspondents from other JET neutron spectrometers is also presented: their agreement allows for JET unique capability of multi-lines of sight neutron spectroscopy and for benchmarking other 14 MeV neutron spectrometers installed on the same lines of sight in preparation for the DT experimental campaign at JET. © 2018 EURATOM.

The connection between plasma physics and neutronics is crucial for the understanding of the operation and performance of modern and future tokamak devices. Neutrons are one of the primary carriers ...

Compact DT neutron generators based on accelerators are often built on the principle of a mixed beam operation, meaning that deuterium (D) and tritium (T) are both present in the ion beam and in the target. Moreover, the beam consists of a mixture of ions and ionized molecules (D, T ions, and ionized D-D, T-T and D-T molecules) so the relevant source components come from T(d, n), D(t, n), D(d, n) and T(t, 2n) reactions at different ion energies. The method for assessing the relative amplitudes of different source components (DD, DT, TT) is presented. The assessment relies on the measurement of the neutron spectrum of different DT components (T(d, n) and D(t, n) at different energies) using a high resolution neutron spectrometer, e.g. a diamond detector, fusion reaction cross-sections, and simulations of neutron generation in the target. Through this process a complete description of the neutron source properties of the mixed beam neutron generator can be made and a neutron source description card, in a format suitable for Monte Carlo code MCNP, produced.

The design and operation of future fusion devices relying on H-mode plasmas requires reliable modelling of edge-localized modes (ELMs) for precise prediction of divertor target conditions. An exten ...

The present paper describes the findings of the calibration of the neutron yield monitors on the Joint European Torus (JET) performed in 2013 using a 252Cf source deployed inside the torus by the remote handling system, with particular regard to the calibration of fission chambers which provide the time resolved neutron yield from JET plasmas. The experimental data obtained in toroidal, radial, and vertical scans are presented. These data are first analysed following an analytical approach adopted in the previous neutron calibrations at JET. In this way, a calibration function for the volumetric plasma source is derived which allows us to understand the importance of the different plasma regions and of different spatial profiles of neutron emissivity on fission chamber response. Neutronics analyses have also been performed to calculate the correction factors needed to derive the plasma calibration factors taking into account the different energy spectrum and angular emission distribution of the calibrating (point) 252Cf source, the discrete positions compared to the plasma volumetric source, and the calibration circumstances. All correction factors are presented and discussed. We discuss also the lessons learnt which are the basis for the on-going 14 MeV neutron calibration at JET and for ITER. © 2017 EURATOM.

The ASCOT Fusion Source Integrator (AFSI) has been used to calculate neutron production rates and spectra corresponding to the JET 19-channel neutron camera (KN3) and the time-of-flight spectrometer (TOFOR) as ideal diagnostics, without detector-related effects. AFSI calculates fusion product distributions in 4D, based on Monte Carlo integration from arbitrary reactant distribution functions. The distribution functions were calculated by the ASCOT Monte Carlo particle orbit following code for thermal, NBI and ICRH particle reactions. Fusion cross-sections were defined based on the Bosch-Hale model and both DD and DT reactions have been included. Neutrons generated by AFSI-ASCOT simulations have already been applied as a neutron source of the Serpent neutron transport code in ITER studies. Additionally, AFSI has been selected to be a main tool as the fusion product generator in the complete analysis calculation chain: ASCOT - AFSI - SERPENT (neutron and gamma transport Monte Carlo code) - APROS (system and power plant modelling code), which encompasses the plasma as an energy source, heat deposition in plant structures as well as cooling and balance-of-plant in DEMO applications and other reactor relevant analyses. This conference paper presents the first results and validation of the AFSI DD fusion model for different auxiliary heating scenarios (NBI, ICRH) with very different fast particle distribution functions. Both calculated quantities (production rates and spectra) have been compared with experimental data from KN3 and synthetic spectrometer data from ControlRoom code. No unexplained differences have been observed. In future work, AFSI will be extended for synthetic gamma diagnostics and additionally, AFSI will be used as part of the neutron transport calculation chain to model real diagnostics instead of ideal synthetic diagnostics for quantitative benchmarking.

Future devices like JT-60SA, ITER and DEMO require quantitative predictions of pedestal density and temperature levels, as well as inter-ELM and ELM divertor heat fluxes, in order to improve global confinement capabilities while preventing divertor erosion/melting in the planning of future experiments. Such predictions can be obtained from dedicated pedestal models like EPED, and from non-linear MHD codes like JOREK, for which systematic validation against current experiments is necessary. In this paper, we show progress in the quantitative validation of the JOREK code using JET simulations. Results analyse the impact of diamagnetic terms on the dynamics and size of the ELMs, and evidence is provided that the onset of type-I ELMs is not governed by linear MHD stability alone, but that a nonlinear threshold could be responsible for large MHD events at the plasma edge. This work has been carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014–2018 under grant agreement No 633053, and from the RCUK Energy Programme (grant number EP/I501045). To obtain further information on the data and models underlying this paper please contact PublicationsManagerccfe.ac.uk. This work used the HELIOS supercomputer (IFERC-CSC), Japan, under the Broader Approach collaboration, implemented by Fusion for Energy and JAEA. The views and opinions expressed herein do not necessarily reflect those of the European Commission or the ITER Organization. The HEC ARCHER computer (UK), as part of the Plasma HEC Consortium EPSRC grant EP/L000237/1, and the MARCONI computer at CINECA in Italy, were also used.

In metallic machines ICRH heating is playing an increasingly important role. One of its most recent applications on the Joint Europena Torus (JET) is sawtooth control by ICRH modulation, for avoiding triggering dangerous neo-classical tearing modes (NTMs) and counteracting impurity accumulation. Some of the main difficulties of these experiments are the assessment of the synchronization efficiency and the understanding of the main physical mechanisms at play. In this paper, three independent classes of statistical indicators are introduced to address these issues: Recurrence Plots, Convergent Cross Mapping and Transfer Entropy. The application to JET experiments with the ILW shows that the proposed indicators agree quite well among themselves and provide sound estimates of the efficiency of the synchronisation scheme investigated. They also support, with a shot to shot basis analysis and an estimate of the uncertainties, the interpretation that the fast ions play a fundamental role in the stabilization of the sawteeth, in both L and H mode. Proposals for experiments to be carried out in the future to consolidate the interpretation of the results are discussed.

In this work, an investigation of the neutral gas flow in the JET sub-divertor area is presented, with respect to the interaction between the plasma side and the pumping side. The edge plasma side is simulated with the SOLPS code, while the sub-divertor area is modeled by means of the direct simulation Monte Carlo (DSMC) method, which in the last few years has proved well able to describe rarefied, collisional flows in tokamak sub-divertor structures. Four different plasma scenarios have been selected, and for each of them a user-defined, iterative procedure between SOLPS and DSMC has been established, using the neutral flux as the key communication term between the two codes. The goal is to understand and quantify the mutual influence between the two regions in a self-consistent manner, that is to say, how the particle exhaust pumping system controls the upstream plasma conditions. Parametric studies of the flow conditions in the sub-divertor, including additional flow outlets and variations of the cryopump capture coefficient, have been performed as well, in order to understand their overall impact on the flow field. The DSMC analyses resulted in the calculation of both the macroscopic quantities—i.e. temperature, number density and pressure—and the recirculation fluxes towards the plasma chamber. The consistent values for the recirculation rates were found to be smaller than those according to the initial standard assumption made by SOLPS.

In several tokamaks, non-axisymmetric magnetic field studies show that applied magnetic fields with a toroidal harmonic n = 2 can lead to disruptive n = 1 locked modes. In Ohmic plasmas, n = 2 magnetic reconnection thresholds in otherwise stable discharges are readily accessed at edge safety factors q similar to 3, low density, and low rotation. Similar to previous studies with n = 1 fields, the thresholds are correlated with the "overlap" field computed with the IPEC code. The overlap field quantifies the plasma-mediated coupling of the external field to the resonant field. Remarkably, the "critical overlap fields" at which magnetic islands form are similar for applied n = 1 and 2 fields. The critical overlap field increases with plasma density and edge safety factor but is independent of the toroidal field. Poloidal harmonics m> nq dominate the drive for resonant fields while m < nq harmonics have a negligible impact. This contrasts with previous results in H-mode discharges at high plasma pressure in which the toroidal angular momentum is sensitive to low poloidal harmonics. Together, these results highlight unique requirements for n > 1 field control including the need for multiple rows of coils to control selected plasma parameters for specific functions (e.g., rotation control or ELM suppression).

Nonlinear 3D MHD asymmetric vertical displacement disruption simulations have been performed using JET equilibrium reconstruction initial data. Several experimentally measured quantities are compared with the simulation. These include vertical displacement, halo current, toroidal current asymmetry, and toroidal rotation. The experimental data and the simulations are in reasonable agreement. Also compared was the correlation of the toroidal current asymmetry and the vertical displacement asymmetry. The Noll relation between asymmetric wall force and vertical current moment is verified in the simulations. Also verified is the toroidal flux asymmetry. Although in many ways, JET is a good predictor of ITER disruption behavior, JET and ITER can be in different parameter regimes, and extrapolating from JET data can overestimate the ITER wall force.

Charge-exchange spectroscopy on JET has become particularly challenging with the introduction of the ITER-like wall. The line intensities are weaker and contaminated by many nuisance lines. We have therefore upgraded the instrumentation to improve throughput and allow the simultaneous measurement of impurity and fuel-ion charge exchange by splitting the light between two pairs of imaging spectrometers using dichroic beam splitters. Imaging instruments allow us to stack 11 × 1 mm diameter fibres on the entrance slits without cross talk. CCD cameras were chosen to have 512 × 512 pixels to allow frame transfer times

The observation that fast ions stabilize ion-temperature-gradient-driven microturbulence has profound implications for future fusion reactors. It is also important in optimizing the performance of present-day devices. In this work, we examine in detail the phenomenology of fast ion stabilization and present a reduced model which describes this effect. This model is derived from the high-energy limit of the gyrokinetic equation and extends the existing "dilution" model to account for nontrivial fast ion kinetics. Our model provides a physically-transparent explanation for the observed stabilization and makes several key qualitative predictions. Firstly, that different classes of fast ions, depending on their radial density or temperature variation, have different stabilizing properties. Secondly, that zonal flows are an important ingredient in this effect precisely because the fast ion zonal response is negligible. Finally, that in the limit of highly-energetic fast ions, their response approaches that of the "dilution" model; in particular, alpha particles are expected to have little, if any, stabilizing effect on plasma turbulence. We support these conclusions through detailed linear and nonlinear gyrokinetic simulations., Comment: 29 pages, 10 figures, 3 tables

This work aims at establishing the parameter space for the existence of geodesic acoustic modes (GAMs) on JET as well as investigating their driving and damping mechanisms predicted by different th ...

Future fusion reactors are foreseen to be heated by the energetic alpha particles produced in fusion reactions. For this to happen, it is important that the energetic ions are sufficiently confined. In present day fusion experiments, energetic ions are primarily produced using external heating systems such as neutral beam injection and ion cyclotron resonance heating. In order to diagnose these fast ions, several different fast-ion diagnostics have been developed and implemented in the various experiments around the world. The velocity-space sensitivities of fast-ion diagnostics are given by so-called weight functions. Here instrument-specific weight functions are derived for neutron emission spectrometry detectors at the tokamaks JET and ASDEX Upgrade for the 2.45 MeV neutrons produced in deuterium-deuterium reactions in deuterium plasmas. Using these, it is possible to directly determine which part of velocity space each detector observes.

A calibration of the JET neutron detectors was carried out prior to the upcoming deuterium-tritium experimental campaign. Two Compact DT neutron generators (NGs) were purchased for this purpose from VNIIA, Russia. These generators are capable of producing approximately 2 × 108 neutrons/s with a DT fusion energy spectrum. Preceding the in-vessel calibration measurements, these compact generators were tested and fully characterised at the UK's National Physical Laboratory (NPL). In order to support the characterisation measurements, detailed neutronics models were developed of the NGs, monitoring detectors and remote handling (RH) apparatus. Neutron spectra calculated from these models have been used to help determine NPL long counter efficiencies and effective centres, as well as NPL reference iron and aluminium activation foil reaction rates. The neutron emission rate has been measured for both generators as a function of angle using absolutely calibrated long counters and the relative emission rate by monitoring single crystal diamond detectors. The measured anisotropy profile is shown to be reproducible with a detailed NG MCNP model. Consequently, the neutron source routine and the MCNP model of the NGs can be reliably used for the analysis of the in-vessel calibration at JET.

At the Joint European Torus (JET) the ex-vessel fission chambers and in-vessel activation detectors are used as the neutron production rate and neutron yield monitors respectively. In order to ensure that these detectors produce accurate measurements they need to be experimentally calibrated. A new calibration of neutron detectors to 14 MeV neutrons, resulting from deuterium–tritium (DT) plasmas, is planned at JET using a compact accelerator based neutron generator (NG) in which a D/T beam impinges on a solid target containing T/D, producing neutrons by DT fusion reactions. This paper presents the analysis that was performed to model the neutron source characteristics in terms of energy spectrum, angle–energy distribution and the effect of the neutron generator geometry. Different codes capable of simulating the accelerator based DT neutron sources are compared and sensitivities to uncertainties in the generator's internal structure analysed. The analysis was performed to support preparation to the experimental measurements performed to characterize the NG as a calibration source. Further extensive neutronics analyses, performed with this model of the NG, will be needed to support the neutron calibration experiments and take into account various differences between the calibration experiment and experiments using the plasma as a source of neutrons.

Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machinespecific minimum current quench time, τCQ; (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, fh . These data reveal that the normalized rotation duration, trot/τCQ, and the average rotation velocity, vh , are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above fh = 20 Hz is to be expected. More importantly, depending on the projected value of τCQ in ITER, substantial rotation could also occur in the resonant frequency range of 6–20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out. EURATOM 633053 RCUK Energy Programme EP/ I501045 Princeton University DE-AC02-09CH11466

This paper presents a diagnostic system, gamma-ray alpha-particle monitor (GRAM), for continuous monitoring of deuterium-tritium fusion alpha-particles in the MeV energy range escaped from the plas ...