The divertor solutions obtained with the plasma edge modelling tool SOLPS5.0 are discussed. The code results are benchmarked against carefully analysed L-mode discharges at various density levels with and without impurity seeding in the full-metal tokamaks ASDEX Upgrade and JET. The role of the cross-field drifts and currents in the solutions is analysed in detail, and the improvements achieved by fully activating the drift and current terms in view of matching the experimental signals are addressed. The persisting discrepancies are also discussed.
The stability with respect to a peeling-ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (omega(*i)), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and omega(*i) effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in omega(*i). The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and w*i effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.
Stability to the type-I edge localized mode (ELM) in JET plasmas was investigated numerically by analyzing the stability to a peeling-ballooning mode with the effects of plasma rotation and ion diamagnetic drift. The numerical analysis was performed by solving the extended Frieman-Rotenberg equation with the MINERVA-DI code. To take into account these effects in the stability analysis self-consistently, the procedure of JET equilibrium reconstruction was updated to include the profiles of ion temperature and toroidal rotation, which are determined based on the measurement data in experiments. With the new procedure and MINERVA-DI, it was identified that the stability analysis including the rotation effect can explain the ELM trigger condition in JET with ITER like wall (JET-ILW), though the stability in JET with carbon wall (JET-C) is hardly affected by rotation. The key difference is that the rotation shear in JET-ILW plasmas analyzed in this study is larger than that in JET-C ones, the shear which enhances the dynamic pressure destabilizing a peeling-ballooning mode. In addition, the increase of the toroidal mode number of the unstable MHD mode determining the ELM trigger condition is also important when the plasma density is high in JET-ILW. Though such modes with high toroidal mode number are strongly stabilized by the ion diamagnetic drift effect, it was found that plasma rotation can sometimes overcome this stabilizing effect and destabilizes the peeling-ballooning modes in JET-ILW.
In tokamaks, the role of turbulent transport of heavy impurities, relative to that of neoclassical transport, increases with increasing size of the plasma, as clarified by means of general scalings, which use the ITER standard scenario parameters as reference, and by actual results from a selection of discharges from ASDEX Upgrade and JET. This motivates the theoretical investigation of the properties of the turbulent convection of heavy impurities by nonlinear gyrokinetic simulations in the experimentally relevant conditions of comparable ion and electron heat fluxes. These conditions also correspond to an intermediate regime between dominant ion temperature gradient turbulence and trapped electron mode turbulence. At moderate plasma toroidal rotation, the turbulent convection of heavy impurities, computed with nonlinear gyrokinetic simulations, is found to be directed outward, in contrast to that obtained by quasi-linear calculations based on the most unstable linear mode, which is directed inward. In this mixed turbulence regime, with comparable electron and ion heat fluxes, the nonlinear results of the impurity transport can be explained by the coexistence of both ion temperature gradient and trapped electron modes in the turbulent state, both contributing to the turbulent convection and diffusion of the impurity. The impact of toroidal rotation on the turbulent convection is also clarified.
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.
Recent progress in the understanding and prediction of the tungsten behaviour in the core of JET H-mode plasmas with ITER-like wall is presented. Particular emphasis is given to the impact of poloidal asymmetries of the impurity density. In particular, it is shown that the predicted reduction of temperature screening induced by the presence of low field side localization of the tungsten density produced by the centrifugal force is consistent with the observed tungsten behaviour in a JET discharge in H-mode baseline scenario. This provides first evidence of the role of poloidal asymmetries in reducing the strength of temperature screening. The main differences between plasma parameters in JET baseline and hybrid scenario discharges which affect the impact of poloidally asymmetric density on the tungsten radial transport are identified. This allows the conditions by which tungsten accumulation can be avoided to be more precisely defined.
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 in proximity to the plasma changes the magnetic topology with consequent effects on the magnetic field structure and the plasma boundary. This paper considers the problem of obtaining the free-boundary plasma equilibrium solution in the presence of ferromagnetic material based on measured constraints. The current approach employs, a model described by O'Brien et al. (1992) in which the magnetisation currents at the iron-air boundary are represented by a set of free parameters and appropriate boundary conditions are enforced via a set of quasi-measurements on the material boundary. This can lead to the possibility of overfitting the data and hiding underlying issues with the measured signals. Although the model typically achieves good fits to measured magnetic signals there are significant discrepancies in the inferred magnetic topology compared with other plasma diagnostic measurements that are independent of the magnetic field. An alternative approach for equilibrium reconstruction in iron-core tokamaks, termed the deterministic magnetisation model is developed and implemented in EFIT++. The iron is represented by a boundary current with the gradients in the magnetisation dipole state generating macroscopic internal magnetisation currents. A model for the boundary magnetisation currents at the iron-air interface is developed using B-Splines enabling continuity to arbitrary order; internal magnetisation currents are allocated to triangulated regions within the iron, and a method to enable adaptive refinement is implemented. The deterministic model has been validated by comparing it with a synthetic 2-D electromagnetic model of JET. It is established that the maximum field discrepancy is less than 1.5 mT throughout the vacuum region enclosing the plasma. The discrepancies of simulated magnetic probe signals are accurate to within 1% for signals with absolute magnitude greater than 100 mT; in all other cases agreement is to within 1 mT. The effect of neglecting the internal magnetisation currents increases the maximum discrepancy in the vacuum region to >20 mT, resulting in errors of 5%-10% in the simulated probe signals. The fact that the previous model neglects the internal magnetisation currents (and also has additional free parameters when fitting the measured data) makes it unsuitable for analysing data in the absence of plasma current. The discrepancy of the poloidal magnetic flux within the vacuum vessel is to within 0.1 Wb. Finally the deterministic model is applied to an equilibrium force-balance solution of a JET discharge using experimental data. It is shown that the discrepancies of the outboard separatrix position, and the outer strike-point position inferred from Thomson Scattering and Infrared camera data are much improved beyond the routine equilibrium reconstruction, whereas the discrepancy of the inner strike-point position is similar. (C) 2017 Published by Elsevier B.V.
Quasi-coherent (QC) modes have been reported for more than 10 years in reflectometry fluctuations spectra in the core region of fusion plasmas. They have characteristics in-between coherent and broadband fluctuations as they oscillate at a marked frequency but have a wide spectrum. This work presents further evidences of the link recently established between QC modes and the trapped electron modes (TEM) instabilities (Arnichand et al 2014 Nucl. Fusion 54 123017). In electron cyclotron resonance heated discharges of Tore Supra, an enhancement of QC modes amplitude is observed in a region where TEM cause impurity transport and turbulence. In JET Ohmic plasmas, QC modes disappear during density ramp-up and current ramp-down. This is reminiscent of Tore Supra and TEXTOR observations during transitions from the linear Ohmic confinement (LOC) to the saturated Ohmic confinement (SOC) regimes. Evidencing TEM activity then becomes experimentally possible via analysis of fluctuation spectra.
The gyrokinetic toroidal code (GTC) has been used to study toroidal Alfven eigenmodes (TAEs) in high-performance plasmas. Experiments performed at the Joint European Torus (JET), where TAEs were driven by energetic particles arising from neutral beams, ion cyclotron resonant heating, and resonantly excited by dedicated external antennas, have been simulated. Modes driven by populations of energetic particles are observed, matching the TAE frequency seen with magnetic probes in JET experiments. A synthetic antenna, composed of one toroidal and two neighboring poloidal harmonics has been used to probe the modes' damping rates and quantify mechanisms for this damping in GTC simulations. This method was also applied to frequency and damping rate measurements of stable TAEs made by the Alfven eigenmode active diagnostic in these discharges.
The physical transport processes at the basis of JET typical inductive H-mode scenarios and advanced hybrid regimes, with improved thermal confinement, are analyzed by means of some of the newest and more sophisticated quasi-linear transport models: trapped gyro Landau fluid (TGLF) and QuaLiKiz. The temporal evolution of JET pulses is modelled by CRONOS where the turbulent transport is modelled by either QuaLiKiz or TGLF. Both are first principle models with a more comprehensive physics than the models previously developed and therefore allow the analysis of the physics at the basis of the investigated scenarios. For H-modes, ion temperature gradient (ITG) modes are found to be dominant and the transport models are able to properly reproduce temperature profiles in self-consistent simulations. However, for hybrid regimes, in addition to ITG trapped electron modes (TEM) are also found to be important and different physical mechanisms for turbulence reduction play a decisive role. Whereas E x B flow shear and plasma geometry have a limited impact on turbulence, the presence of a large population of fast ions, quite important in low density regimes, can stabilize core turbulence mainly when the electromagnetic effects are taken into account. The TGLF transport model properly captures these mechanisms and correctly reproduces temperatures.
The density evolution during the transient phase just after the L-H transition is investigated using theoretical transport models. Cases characterized by core densities which evolve in longer timescales than the edge densities, leading to hollow density profiles (R/L-n = -R del n/n < 0) are modelled. This density evolution is particularly interesting because it has been shown to be beneficial in the view of the access to burning plasma conditions in ITER (Loarte et al 2013 Nucl. Fusion 53 083031). Self-consistent simulations of the JET discharge 79676 of the density-only, and of the density and the temperatures are carried out using a quasilinear gyrokinetic code, QuaLiKiz (Bourdelle et al 2007 Phys. Plasmas 14 112501), coupled with a transport code CRONOS (Artaud et al 2010 Nucl. Fusion 50 043001). The slow evolution of the hollow density, associated with the self-consistently calculated hollow NBI particle deposition, is well reproduced in the plasma core. Indeed, QuaLiKiz is shown to reproduce nonlinear gyrokinetic heat and particle fluxes well for both positive and negative R/L-n. That gives a theoretical and general basis for the persistence of the hollowness, laying the groundwork for the extrapolation to ITER.
The phenomena of retention and de-trapping of deuterium (D) and tritium (T) in plasma facing components (PFC) and supporting structures must be understood in order to limit or control total T inventory in larger future fusion devices such as ITER, DEMO and commercial machines. The goal of this paper is to present details of the thermal desorption spectrometry (TDS) system applied in total fuel retention assessment of PFC at the Joint European Torus (JET). Examples of TDS results from beryllium (Be) wall tile samples exposed to JET plasma in PFC configuration mirroring the planned ITER PFC is shown for the first time. The method for quantifying D by comparison of results from a sample of known D content was confirmed acceptable. The D inventory calculations obtained from Ion Beam Analysis (IBA) and TDS agree well within an error associated with the extrapolation from very few data points to a large surface area.
Results of the first dust survey in JET with the ITER-Like Wall (JET-ILW) are presented. The sampling was performed using adhesive stickers from the divertor tiles where the greatest material deposition was detected after the first JET-ILW campaign in 2011-2012. The emphasis was especially on sampling and analysis of metal particles (Be and W) with the aim to determine the composition, size, surface topography and internal dust structure using a large set of methods: high-resolution scanning and transmission electron microscopy, focused ion beam, electron diffraction and also wavelength and energy dispersive x-ray spectroscopy. The most important was the identification of beryllium dust both in the form of flakes and droplets with dimensions in the micrometer range. Tungsten, molybdenum, inconel constituents were identified along with many impurity species. The particles are categorised and the origin of the various constituents discussed.
The neoclassical tearing modes (NTM) increase the effective heat and particle radial transport inside the plasma, leading to a flattening of the electron and ion temperature and density profiles at a given location depending on the safety factor q rational surface (Hegna and Callen 1997 Phys. Plasmas 4 2940). In burning plasma such as in ITER, this NTM-induced increased transport could reduce significantly the fusion performance and even lead to a disruption. Validating models describing the NTM-induced transport in present experiment is thus important to help quantifying this effect on future devices. In this work, we apply an NTM model to an integrated simulation of current, heat and particle transport on JET discharges using the European transport simulator. In this model, the heat and particle radial transport coefficients are modified by a Gaussian function locally centered at the NTM position and characterized by a full width proportional to the island size through a constant parameter adapted to obtain the best simulations of experimental profiles. In the simulation, the NTM model is turned on at the same time as the mode is triggered in the experiment. The island evolution is itself determined by the modified Rutherford equation, using self-consistent plasma parameters determined by the transport evolution. The achieved simulation reproduces the experimental measurements within the error bars, before and during the NTM. A small discrepancy is observed on the radial location of the island due to a shift of the position of the computed q = 3/2 surface compared to the experimental one. To explain such small shift (up to about 12% with respect to the position observed from the experimental electron temperature profiles), sensitivity studies of the NTM location as a function of the initialization parameters are presented. First results validate both the transport model and the transport modification calculated by the NTM model.
The power output of fusion devices is measured from their neutron yields which relate directly to the fusion yield. In this paper we describe the devices and methods that have been prepared to perform a new in situ 14 MeV neutron calibration at JET in view of the new DT campaign planned at JET in the next years. The target accuracy of this calibration is 10% as required for ITER, where a precise neutron yield measurement is important, e.g., for tritium accountancy. In this paper, the constraints and early decisions which defined the main calibration approach are discussed, e.g., the choice of 14 MeV neutron source and the deployment method. The physics preparations, source issues, safety and engineering aspects required to calibrate directly the JET neutron detectors are also discussed. The existing JET remote-handling system will be used to deploy the neutron source inside the JET vessel. For this purpose, compatible tooling and systems necessary to ensure safe and efficient deployment have been developed. The scientific programme of the preparatory phase is devoted to fully characterizing the selected 14 MeV neutron generator to be used as the calibrating source, obtain a better understanding of the limitations of the calibration, optimise the measurements and other provisions, and to provide corrections for perturbing factors (e.g., anisotropy of the neutron generator, neutron energy spectrum dependence on emission angle). Much of this work has been based on an extensive programme of Monte-Carlo calculations which provide support and guidance in developing the calibration strategy. i3/4 (C) 2017 EURATOM.
The design and operation of ITER experimental fusion reactor requires the development of neutron measurement techniques and numerical tools to derive the fusion power and the radiation field in the device and in the surrounding areas. Nuclear analyses provide essential input to the conceptual design, optimisation, engineering and safety case in ITER and power plant studies. The required radiation transport calculations are extremely challenging because of the large physical extent of the reactor plant, the complexity of the geometry, and the combination of deep penetration and streaming paths. This article reports the experimental activities which are carried-out at JET to validate the neutronics measurements methods and numerical tools used in ITER and power plant design. A new deuterium-tritium campaign is proposed in 2019 at JET: the unique 14 MeV neutron yields produced will be exploited as much as possible to validate measurement techniques, codes, procedures and data currently used in ITER design thus reducing the related uncertainties and the associated risks in the machine operation.
In view of the planned DT operations at JET, a calibration of the JET neutron monitors at 14 MeV neutron energy is needed using a 14 MeV neutron generator deployed inside the vacuum vessel by the JET remote handling system. The target accuracy of this calibration is +/- 10% as also required by ITER, where a precise neutron yield measurement is important, e.g. for tritium accountancy. To achieve this accuracy, the 14 MeV neutron generator selected as the calibration source has been fully characterised and calibrated prior to the in-vessel calibration of the JET monitors. This paper describes the measurements performed using different types of neutron detectors, spectrometers, calibrated long counters and activation foils which allowed us to obtain the neutron emission rate and the anisotropy of the neutron generator, i.e. the neutron flux and energy spectrum dependence on emission angle, and to derive the absolute emission rate in 4 pi sr. The use of high resolution diamond spectrometers made it possible to resolve the complex features of the neutron energy spectra resulting from the mixed D/T beam ions reacting with the D/T nuclei present in the neutron generator target. As the neutron generator is not a stable neutron source, several monitoring detectors were attached to it by means of an ad hoc mechanical structure to continuously monitor the neutron emission rate during the in-vessel calibration. These monitoring detectors, two diamond diodes and activation foils, have been calibrated in terms of neutrons/counts within +/- 5% total uncertainty. A neutron source routine has been developed, able to produce the neutron spectra resulting from all possible reactions occurring with the D/T ions in the beam impinging on the Ti D/T target. The neutron energy spectra calculated by combining the source routine with a MCNP model of the neutron generator have been validated by the measurements. These numerical tools will be key in analysing the results from the in-vessel calibration and to derive the response of the JET neutron detectors to DT plasma neutrons starting from the response to the generator neutrons, and taking into account all the calibration circumstances.
Rotating collectors and quartz microbalances (QMBs) are used in JET to provide time-dependent measurements of erosion and deposition. Rotation of collector discs behind apertures allows recording of the long term evolution of deposition. QMBs measure mass change via the frequency deviations of vibrating quartz crystals. These diagnostics are used to investigate erosion/deposition during JET-C carbon operation and JET-ILW (ITER-like wall) beryllium/tungsten operation. A simple geometrical model utilising experimental data is used to model the time-dependent collector deposition profiles, demonstrating good qualitative agreement with experimental results. Overall, the JET-ILW collector deposition is reduced by an order of magnitude relative to JET-C, with beryllium replacing carbon as the dominant deposit. However, contrary to JET-C, in JET-ILW there is more deposition on the outer collector than the inner. This reversal of deposition asymmetry is investigated using an analysis of QMB data and is attributed to the different chemical properties of carbon and beryllium.
This paper reports on the procedure developed as the best method to provide an accurate and reliable estimation of the ion temperature T-i and the toroidal velocity v(phi) from Charge-eXchange Recombination Spectroscopy (CXRS) data from intrinsic rotation experiments at the Joint European Torus with the carbon wall. The low impurity content observed in such plasmas, resulting in low active CXRS signal, alongside low Doppler shifts makes the determination of Ti and v(phi) particularly difficult. The beam modulation method will be discussed along with the measures taken to increase photon statistics and minimise errors from the absolute calibration and magneto-hydro-dynamics effects that may impact the CXRS passive emission.
Future fusion reactors require a safe, steady state divertor operation. A possible solution for the power exhaust challenge is the detached divertor operation in scenarios with high radiated power fractions. The radiation can be increased by seeding impurities, such as N for dominant scrape-off-layer radiation, Ne or Ar for SOL and pedestal radiation and Kr for dominant core radiation. Recent experiments on two of the all-metal tokamaks, ASDEX Upgrade (AUG) and JET, demonstrate operation with high radiated power fractions and a fully-detached divertor by N, Ne or Kr seeding with a conventional divertor in a vertical target geometry. For both devices similar observations can be made. In the scenarios with the highest radiated power fraction, the dominant radiation originates from the confined region, in the case of N and Ne seeding concentrated in a region close to the X-point. Applying these seed impurities for highly radiative scenarios impacts local plasma parameters and alters the impurity transport in the pedestal region. Thus, plasma confinement and stability can be affected. A proper understanding of the effects by these impurities is required in order to predict the applicability of such scenarios for future devices. (C) 2017 Elsevier Ltd.
In this paper, we present the results obtained from the data analysis of neutron spectra measured with a NE213 liquid scintillator at JET. We calculated the neutron response matrix of the instrument combining MCNPX simulations, a generic proton light output function measured with another detector and the fit of data from ohmic pulses. For the analysis, we selected a set of pulses with neutral beam injection heating (NBI) only and we applied a forward fitting procedure of modeled spectral components to extract the fraction of thermal neutron emission. The results showed the same trend of the ones obtained with the dedicated spectrometer TOFOR, even though the values from the NE213 analysis were systematically higher. This discrepancy is probably due to the different lines of sight of the two spectrometers (tangential for the NE213, vertical for TOFOR). The uncertainties on the thermal fraction estimates were from 4 to 7 times higher than the ones from the TOFOR analysis.
The good performance of a detector can be strongly affected by the instrumentation used to acquire the data. The possibility of anticipating how the acquisition chain will affect the signal can help in finding the best solution among different set-ups. In this work we developed a Monte Carlo code that aims to simulate the effect of the various components of a digital Data Acquisition system (DAQ) applied to scintillation detectors. The components included in the model are: the scintillator, the photomultiplier tube (PMT), the signal cable and the digitizer. We benchmarked the code against real data acquired with a NE213 scintillator, comparing simulated and real signal pulses induced by gamma-ray interaction. Then we studied the dependence of the energy resolution of a pulse height spectrum (PHS) on the sampling frequency and the bit resolution of the digitizer. We found that exceeding some values of the sampling frequency and the bit resolution improves only marginally the performance of the system. The method can be applied for the study of various detector systems relevant for nuclear techniques, such as in fusion diagnostics.
In this paper we investigate the possibility of using a NE213 liquid scintillator as aneutron spectrometer to diagnose the fast ion tails produced in experiments with 3rd harmonicradio-frequency heating.We discuss mainly the instrumental effects that need to be considered and corrected for in orderto obtain a good agreement between measured data and models: gain drift, pile-up, impact of theassumption of a standard proton light yield function. We also address problems related to thepresence of triton burn-up events in the spectrum.The expected ion distribution is obtained from a simple 1D Fokker-Planck model. The parametersof the model are estimated using the data collected by the TOFOR neutron spectrometer.The agreement between the data and the model is good and it is possible to make a clear distinctionbetween discharges that had different electron densities and thus different cut-off energies. Wecan conclude that NE213 scintillators can provide useful spectroscopic information for this kind ofexperiments.
In this work we present a method to evaluate the neutron response function of an NE213 liquid scintillator. This method is particularly useful when the proton light yield function of the detector has not been measured, since it is based on a proton light yield function taken from literature, MCNPX simulations, measurements of gammarays from a calibration source and measurements of neutrons from fusion experiments with ohmic plasmas. The inclusion of the latter improves the description of the proton light yield function in the energy range of interest (around 2.46 MeV). We apply this method to an NE213 detector installed at JET, inside the radiation shielding of the magnetic proton recoil (MPRu) spectrometer, and present the results from the calibration along with some examples of application of the response function to perform neutron emission spectroscopy (NES) of fusion plasmas. We also investigate how the choice of the proton light yield function affects the NES analysis, finding that the result does not change significantly. This points to the fact that the method for the evaluation of the neutron response function is robust and gives reliable results.
In AC/DC converters, a peculiar periodic nonsmooth waveform arises, the so-called ripple. In this paper we propose a novel model that captures this nonsmoothness by means of a hybrid dynamical system performing state jumps at certain switching instants, and we illustrate its properties with reference to a three phase diode bridge rectifier. As the ripple corrupts an underlying desirable signal, we propound two observer schemes ensuring asymptotic estimation of the ripple, the first with and the second without knowledge of the switching instants. Our theoretical developments are well placed in the context of recent techniques for hybrid regulation and constitute a contribution especially for our second observer, where the switching instants are estimated. Once asymptotic estimation of the ripple is achieved, the ripple can be conveniently canceled from the desirable signal, and thanks to the inherent robustness properties of the proposed hybrid formulation, the two observer schemes require only that the desirable signal is slowly time varying compared to the ripple. Exploiting this fact, we illustrate the effectiveness of our second hybrid observation law on experimental data collected from the Joint European Torus tokamak.
Use of new 3-strap ICRF antennas with all-tungsten (W) limiters in ASDEX Upgrade results in a reduction of the W sources at the antenna limiters and of the W content in the confined plasma by at least a factor of 2 compared to the W-limiter 2-strap antennas used in the past. The reduction is observed with a broad range of plasma shapes. In multiple locations of antenna frame, the limiter W source has a minimum when RF image currents are decreased by cancellation of the RF current contributions of the central and the outer straps. In JET with ITER-like wall, ITER-like antenna produces about 20% less of main chamber radiation and of W content compared to the old A2 antennas. However the effect of the A2 antennas on W content is scattered depending on which antennas are powered. Experiments in JET with trace nitrogen (N-2) injection show that a presence of active ICRF antenna close to the midplane injection valve has little effect on the core N content, both in dipole and in -90 degrees phasing. This indicates that the effect of ICRF on impurity transport across the scape-off-layer is small in JET compared to the dominant effect on impurity sources leading to increased impurity levels during ICRF operation. (C) 2016 The Authors. Published by Elsevier Ltd.
Recent achievements in studies of the effects of ICRF (Ion Cyclotron Range of Frequencies) power on the SOL (Scrape-Off Layer) and PWI (Plasma Wall Interactions) in ASDEX Upgrade (AUG), Alcator C-Mod, and JET-ILW are reviewed. Capabilities to diagnose and model the effect of DC biasing and associated impurity production at active antennas and on magnetic field connections to antennas are described. The experiments show that ICRF near-fields can lead not only to E x B convection, but also to modifications of the SOL density, which for Alcator C-Mod are limited to a narrow region near antenna. On the other hand, the SOL density distribution along with impurity sources can be tailored using local gas injection in AUG and JET-ILW with a positive effect on reduction of impurity sources. The technique of RF image current cancellation at antenna limiters was successfully applied in AUG using the 3-strap AUG antenna and extended to the 4-strap Alcator C-Mod field-aligned antenna. Multiple observations confirmed the reduction of the impact of ICRF on the SOL and on total impurity production when the ratio of the power of the central straps to the total antenna power is in the range 0.6 < P-cen / P-total < 0.8. Near-field calculations indicate that this fairly robust technique can be applied to the ITER ICRF antenna, enabling the mode of operation with reduced PWI. On the contrary, for the A2 antenna in JET-ILW the technique is hindered by RF sheaths excited at the antenna septum. Thus, in order to reduce the effect of ICRF power on PWI in a future fusion reactor, the antenna design has to be optimized along with design of plasmafacing components.
The paper deals with radiation resistant sensors and their associated measuring instrumentation developed in the course of R and D activities carried out in the framework of an international collaboration. The first trial tests of three-dimensional (3D) probes with Hall sensors have been performed in European tokamaks TORE SUPRA (2004) and JET (2005). Later in 2009 six sets of 3D probes were installed in JET and now continue to operate. The statistical analysis performed in 2014 on the basis of the JET database have demonstrated stable long term operation of all 18 sensors of 3D probes. The results of measurements conducted at the neutron fluxes of nuclear reactors have demonstrated the operability of the sensors up to high neutron fluences of F > 10(18)n , cm(-2) that exceeds the maximum one for the locations of steady state sensors in ITER over its total lifetime.
The products of fusion reactions at JET are measured using different diagnostic techniques. One of the methods is based on measurements of gamma-rays, originating from reactions between fast ions and plasma impurities. During the forthcoming deuterium-tritium (DT) campaign a particular attention will be paid to 4.44 MeV gamma-rays emitted in the Be-9(alpha,n gamma)C-12 reaction. Gamma-ray detectors foreseen for measurements in DT campaigns have to be able to register spectra at high count rates, up to approximately 500 kHz. For the Gamma-ray Camera at JET a new setup will be based on scintillators with a short decay time, e.g., CeBr3, and a multi-pixel photon counter (MPPC). We present two methods of output signal shortening in modules based on MPPC. A short detector output signal is necessary in order to minimize the number of pile up events at high count rates. One method uses a passive RC circuit with a pole zero cancellation, whereas an active transimpedance amplifier is used in the other one. Due to the strong dependence of MPPC properties on temperature variation, a special device MTCD@NCBJ was designed and produced to stabilize the gain in MPPC-based scintillation detectors. We show that this device guarantees stable working conditions.
A set of experiments was carried out in JET ILW (Joint European Torus with ITER-Like Wall) L-mode plasmas in order to study the effects of light impurities on core ion thermal transport. N was puffed into some discharges and its profile was measured by active Charge Exchange diagnostics, while ICRH power was deposited on- and off-axis in (He-3) - D minority scheme in order to have a scan of local heat flux at constant total power with and without N injection. Experimentally, the ion temperature profiles are more peaked for similar heat fluxes when N is injected in the plasma. Gyro-kinetic simulations using the GENE code indicate that a stabilization of Ion Temperature Gradient driven turbulent transport due to main ion dilution and to changes in T-e/T-i and s/q is responsible of the enhanced peaking. The quasi-linear models TGLF and QuaLiKiz are tested against the experimental and the gyro-kinetic results.
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.
A series of experimental observations of light impurity profiles was carried out in JET (Joint European Torus) ITER-like wall (ILW) L-mode plasmas in order to investigate their transport mechanisms. These discharges feature the presence of He-3, Be, C, N, Ne, whose profiles measured by active Charge Exchange diagnostics are compared with quasi-linear and non-linear gyro-kinetic simulations. The peaking of He-3 density follows the electron density peaking, Be and Ne are also peaked, while the density profiles of C and N are flat in the mid plasma region. Gyro-kinetic simulations predict peaked density profiles for all the light impurities studied and at all the radial positions considered, and fail predicting the flat or hollow profiles observed for C and N at mid radius in our cases.
This paper presents results of JET ITER-like wall L-mode experiments in hydrogen and deuterium (D) plasmas, dedicated to the study of the isotope dependence of ion heat transport by determination of the ion critical gradient and stiffness by varying the ion cyclotron resonance heating power deposition. When no strong role of fast ions in the plasma core is expected, the main difference between the two isotope plasmas is determined by the plasma edge and the core behavior is consistent with a gyro-Bohm scaling. When the heating power (and the fast ion pressure) is increased, in addition to the difference in the edge region, also the plasma core shows substantial changes. The stabilization of ion heat transport by fast ions, clearly visible in D plasmas, appears to be weaker in H plasmas, resulting in a higher ion heat flux in H with apparent anti-gyro-Bohm mass scaling. The difference is found to be caused by the different fast ion pressure between H and D plasmas, related to the heating power settings and to the different fast ion slowing down time, and is completely accounted for in non-linear gyrokinetic simulations. The application of the TGLF quasi-linear model to this set of data is also discussed.
The main purpose of this work is to study the dependence of trapped electron modes (TEM) threshold and of electron stiffness on the most relevant plasma parameters. Dedicated transport experiments based on heat flux scans and T-e modulation have been performed in JET in TEM dominated plasmas with pure ICRH electron heating and a numerical study using gyrokinetic simulations has been performed with the code GKW. Using multilinear regressions on the experimental data, the stabilizing effect of magnetic shear predicted by theory for our plasma parameters is confirmed while no significant effect of safety factor was found. Good quantitative agreement is found between the TEM thresholds found in the experiments and calculated with linear GKW simulations. Non-linear simulations have given further confirmation of the threshold values and allowed comparison with the values of stiffness found experimentally. Perturbative studies using RF power modulation indicate the existence of an inward convective term for the electron heat flux. Adding NBI power, ion temperature gradient (ITG) modes become dominant and a reduction of vertical bar del T-e vertical bar/T-e with respect to pure ICRH, TEM dominant discharges has been experimentally observed, in spite of increased total electron power. Possible explanations 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.
Physical and chemical assisted physical sputtering were characterised by the Be I and Be II line and BeD band emission in the observation chord measuring the sightline integrated emission in front of the inner beryllium limiter at the torus midplane. The 3D local transport and plasma-surface interaction Monte-Carlo modelling (ERO code [18]) is a key for the interpretation of the observations in the vicinity of the shaped solid Be limiter. The plasma parameter variation (density scan) in limiter regime has provided a useful material for the simulation benchmark. The improved background plasma parameters input, the new analytical expression for particle tracking in the sheath region and implementation of the BeD release into ERO has helped to clarify some deviations between modelling and experiments encountered in the previous studies [4,5]. Reproducing the observations provides additional confidence in our 'ERO-min' fit for the physical sputtering yields for the plasma-wetted areas based on simulated data.
ERO is a 3D Monte-Carlo impurity transport and plasma-surface interaction code. In 2011 it was applied for the ITER first wall (FW) life time predictions [1] (critical blanket module BM11). After that the same code was significantly improved during its application to existing fusion-relevant plasma devices: the tokamak JET equipped with an ITER-like wall and linear plasma device PISCES-B. This has allowed testing the sputtering data for beryllium (Be) and showing that the "ERO-min" fit based on the large (50%) deuterium (D) surface content is well suitable for plasma-wetted areas (D plasma). The improved procedure for calculating of the effective sputtering yields for each location along the plasma-facing surface using the recently developed semi-analytical sheath approach was validated. The re-evaluation of the effective yields for BM11 following the similar revisit of the JET data has indicated significant increase of erosion and motivated the current re-visit of ERO simulations.
For simulation of plasma-facing component erosion in fusion experiments, an analytical expression for the ion velocity just before the surface impact including the local electric field and an optional surface biasing effect is suggested. Energy and angular impact distributions and the resulting effective sputtering yields were produced for several experimental scenarios at JET ILW mostly involving PFCs exposed to an oblique magnetic field. The analytic solution has been applied as an improvement to earlier ERO modelling of localized, Be outer limiter, RF-enhanced erosion, modulated by toggling of a remote, however magnetically connected ICRH antenna. The effective W sputtering yields due to D and Be ion impact in Type-I and Type-III ELMs and inter-ELM conditions were also estimated using the analytical approach and benchmarked by spectroscopy. The intra-ELM W sputtering flux increases almost 10 times in comparison to the inter-ELM flux. (C) 2017 Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license.
A new analytical approximation for the electric potential profile in the presence of an oblique magnetic field and the analytical solution for the particle motion just before the impact with a plasma-facing surface are presented. These approximations are in good agreement with fluid solutions and the corresponding PIC simulations. These expressions were applied to provide effective physical erosion yields for Be, which have in a second step been used in ERO code simulations of spectroscopy at Be limiters of the JET ITER-like wall. These new analytical expressions lead to an increase of the effective physical sputtering yields of Be by deuteron impact up to 30% in comparison with earlier pure numerical simulations. ((c) 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)
The W source strength in JET H-mode discharges depends on the W sputtering in the inter and the intra-ELM phase due to impinging hydrogenic ions (D or H) and impurities (mainly Be). The analytical approach for interpretation of the Langmuir probe measurements is applied to model the ELM ion parallel transport and the W sputtering flux in intra-ELM and inter-ELM conditions in JET-ITER-like wall (ILW) hydrogen and deuterium plasmas. The impact of the Be ion charge and the Be concentration in the impinging ion flux on the W sputtering was estimated. Be2+ concentrations of 0.5% and 1% in the impinging ion flux increases the W sputtering fluence per ELM by 20%-30% and 35%-55% correspondingly with respect to pure deuterium plasma; the charge state of Be ions has no substantial effect on W sputtering in the intra-ELM phase. The analysis of JET ELMy H-mode discharges in hydrogen with different types of ELMs is presented. The W sputtering source under inter-and intra-ELM conditions is estimated using the analytical approach and validated by optical emission spectroscopy in these discharges. The intra-ELM W sputtering flux increases 2-4 times in comparison to the inter-ELM flux.
Nonlinear gyrokinetic codes allow for detailed understanding of tokamak core turbulent transport. However, their computational demand precludes their use for predictive profile modeling. An alternative approach is required to bridge the gap between theoretical understanding and prediction of experiments. A quasilinear gyrokinetic model, QuaLiKiz (Bourdelle et al 2007 Phys. Plasmas 14 112501), is demonstrated to be rapid enough to ease systematic interface with experiments. The derivation and approximation of this approach are reviewed. The quasilinear approximation is proven valid over a wide range of core plasma parameters. Examples of profile prediction using QuaLiKiz coupled to the CRONOS integrated modeling code (Artaud et al 2010 Nucl. Fusion 50 043001) are presented. QuaLiKiz is being coupled to other integrated modeling platforms such as ETS and JETTO. QuaLiKiz quasilinear gyrokinetic turbulent heat, particle and angular momentum fluxes are available to all users. It allows for extensive stand-alone interpretative analysis and for first principle based integrated predictive modeling.
In ion-temperature-gradient (ITG) driven turbulence, the resonance condition leads to ion particle turbulent transport coefficients significantly larger than electron particle turbulent transport coefficients. This is shown in nonlinear gyrokinetic simulations and explained by an analytical quasilinear model. It is then illustrated by JETTO-QuaLiKiz integrated modelling. Large ion particle transport coefficients implies that the ion density profiles are uncorrelated to the corresponding ion source, allowing peaked isotope density profiles even in the absence of core source. This also implies no strong core accumulation of He ash. Furthermore, the relaxation time of the individual ion profiles in a multi-species plasma can be significantly faster than the total density profile relaxation time which is constrained by the electrons. This leads to fast isotope mixing and fast impurity transport in FM regimes. In trapped-electron- mode (TEM) turbulence, in presence of electron heating about twice the ion heating, the situation is the inverse: ion particle turbulent transport coefficients are smaller than their electron counterpart.
The L to H mode transition occurs at a critical power which depends on various parameters, such as the magnetic field, the density, etc. Experimental evidence on various tokamaks (JET, ASDEX-Upgrade, DIII-D, Alcator C-Mod) points towards the existence of a critical temperature characterizing the transition. This criterion for the L-H transition is local and is therefore easier to be compared to theoretical approaches. In order to shed light on the mechanisms of the transition, simple theoretical ideas are used to derive a temperature threshold (T-th). They are based on the stabilization of the underlying turbulence by a mean radial electric field shear. The nature of the turbulence varies as the collisionality decreases, from resistive ballooning modes to ion temperature gradient and trapped electron modes. The obtained parametric dependencies of the derived T-th are tested versus magnetic field, density, effective charge. Various robust experimental observations are reproduced, in particular T-th increases with magnetic field B and increases with density below the density roll-over observed on the power threshold.
With WEST (Tungsten Environment in Steady State Tokamak) (Bucalossi et al 2014 Fusion Eng. Des. 89 907-12), the Tore Supra facility and team expertise (Dumont et al 2014 Plasma Phys. Control. Fusion 56 075020) is used to pave the way towards ITER divertor procurement and operation. It consists in implementing a divertor configuration and installing ITER-like actively cooled tungsten monoblocks in the Tore Supra tokamak, taking full benefit of its unique long-pulse capability. WEST is a user facility platform, open to all ITER partners. This paper describes the physics basis of WEST: the estimated heat flux on the divertor target, the planned heating schemes, the expected behaviour of the L-H threshold and of the pedestal and the potential W sources. A series of operating scenarios has been modelled, showing that ITER-relevant heat fluxes on the divertor can be achieved in WEST long pulse H-mode plasmas.
The pressure gradient of the high confinement pedestal region at the edge of tokamak plasmas rapidly collapses during plasma eruptions called edge localised modes (ELMs), and then re-builds over a longer time scale before the next ELM. The physics that controls the evolution of the JET pedestal between ELMs is analysed for 1.4 MA, 1.7 T, low triangularity, delta = 0.2, discharges with the ITER-like wall, finding that the pressure gradient typically tracks the ideal magneto-hydrodynamic ballooning limit, consistent with a role for the kinetic ballooning mode. Furthermore, the pedestal width is often influenced by the region of plasma that has second stability access to the ballooning mode, which can explain its sometimes complex evolution between ELMs. A local gyrokinetic analysis of a second stable flux surface reveals stability to kinetic ballooning modes; global effects are expected to provide a destabilising mechanism and need to be retained in such second stable situations. As well as an electronscale electron temperature gradient mode, ion scale instabilities associated with this flux surface include an electro-magnetic trapped electron branch and two electrostatic branches propagating in the ion direction, one with high radial wavenumber. In these second stability situations, the ELM is triggered by a peeling-ballooning mode; otherwise the pedestal is somewhat below the peeling-ballooning mode marginal stability boundary at ELM onset. In this latter situation, there is evidence that higher frequency ELMs are paced by an oscillation in the plasma, causing a crash in the pedestal before the peeling-ballooning boundary is reached. A model is proposed in which the oscillation is associated with hot plasma filaments that are pushed out towards the plasma edge by a ballooning mode, draining their free energy into the cooler plasma there, and then relaxing back to repeat the process. The results suggest that avoiding the oscillation and maximising the region of plasma that has second stability access will lead to the highest pedestal heights and, therefore, best confinement-a key result for optimising the fusion performance of JET and future tokamaks, such as ITER.
A new Safety Case is required to permit tritium operations on JET during the forthcoming DTE2 campaign. The outputs, benefits and lessons learned associated with the production of this Safety Case are presented. The changes that have occurred to the Safety Case methodology since the last JET tritium Safety Case are reviewed. Consideration is given to the effects of modifications, particularly ITER related changes, made to the JET and the impact these have on the hazard assessments as well as normal operations. Several specialized assessments, including recent MELCOR modelling, have been undertaken to support the production of this Safety Case and the impact of these assessments is outlined. Discussion of the preliminary actions being taken to progress implementation of this Safety Case is provided, highlighting new methods to improve the dissemination of the key Safety Case results to the plant operators. Finally, the work required to complete this Safety Case, before the next tritium campaign, is summarized. (C) 2016 EURATOM. Published by Elsevier B.V. All rights reserved.
Nonlinear gyrokinetic simulations using the GENE code have previously predicted a significant nonlinear enhanced electromagnetic stabilization in certain JET discharges with high neutral-beam power and low core magnetic shear (Citrin et al 2013 Phys. Rev. Lett. 111 155001, 2015 Plasma Phys. Control. Fusion 57 014032). This dominates over the impact of E x B flow shear in these discharges. Furthermore, fast ions were shown to be a major contributor to the electromagnetic stabilization. These conclusions were based on results from the GENE gyrokinetic turbulence code. In this work we verify these results using the GYRO code. Comparing results (linear frequencies, eigenfunctions, and nonlinear fluxes) from different gyrokinetic codes as a means of verification (benchmarking) is only convincing if the codes agree for more than one discharge. Otherwise, agreement may simply be fortuitous. Therefore, we analyze three discharges, all with a carbon wall: a simplified, two-species, circular geometry case based on an actual JET discharge; an L-mode discharge with a significant fast-ion pressure fraction; and a low-triangularity high-beta hybrid discharge. All discharges were analyzed at normalized toroidal flux coordinate rho = 0.33 where significant ion temperature peaking is observed. The GYRO simulations support the conclusion that electromagnetic stabilization is strong, and dominates E x B shear stabilization.