Search

An atmospheric Pressure Field Effect (APFE) ionisation source for drift tube ion mobility spectrometry has been developed for operation in positive and negative polarities. The formation of negative and positive ions in synthetic air was studied and compared with the Atmospheric Pressure Corona Discharge (APCD) ionisation source. The APFE ionisation source is of point-to-plane geometry with a 10 μm Pt point electrode, a stainless steel plate electrode and ultra-high resistance (20 GΩ) current limiters. In the case of negative polarity, the ionisation source was able to generate Reactant Ions (RIs) O

In astronomy emission spectroscopy is widely used technique to study various objects. Although emission spectra of individual compounds are quite well known so composition of the studied body can be analysed the information on processes active in distant object can be encoded in the determined spectra as well. To decode such information the individual elementary processes need to be well understood. For example, in case of comets which are simple leftovers from the beginning of the Solar System understanding the connection between the coma and the comet’s nucleus is critical because observations rarely detect the nucleus directly, and its properties must often be inferred from measurements of the coma. As measurements of the coma do not necessarily represent the characteristics of the nucleus due to spatial, temporal and chemical evolution of the emitted material, projecting the coma observations back to the nucleus requires understanding of the processes in the coma. Within the Rosetta various instruments were used for spectroscopy. OSIRIS was equipped with several band filters that are commonly used in ground-based studies (Farnham et al. 2000). The spectroscopic measurements done within the Rosetta mission led to conclusion that electron-molecule processes play much more significant role than it used to be suspected in the past (Bodewits et al. 2016). This fact can be utilized to remotely study the plasma-chemical reactions active in distant objects. For example, the formation of excited OH by electron impact on H2O has an electron temperature threshold of 9.1 eV; H Ly-a has a threshold of 15.4 eV; and the production of OI 130.4 nm emission requires electron temperatures above 23.5 eV (Beenakker et al. 1974; Morgan & Mentall 1974). At electron temperatures below 10 eV, collisions with water molecules result in ro-vibrational excitation and energy is emitted in the IR range (Itikawa & Mason 2005). Electron impact emission thus provides a sensitive diagnostic probe for studying coma (Bodewits et al. 2019). The presence of spectral structures depends on reaching the threshold electron energy for specific transition. Emission cross sections that have thresholds of 10-100 eV – typical for the coma (Broiles et al. 2016) – depend strongly on electron impact energy. Knowledge of emission spectra and cross sections for electron-molecule reactions is critical for evaluation of physical-chemical properties of cometary coma based on spectroscopic observations. For thorough analysis of such studies a reliable database of reference spectra and emission cross sections must be available. Especially the availability of cross section data the Achilles’ heel since most of the values were determined only for a few individual energies but for successful application in models of atmospheres or comas they need to be known for wider range of energies. Such data can be provided by laboratory experiments studying individual electron-molecule interactions in various experimental configurations such as crossed beams of electrons and molecules. The example of the apparatus scheme aimed at studying electron induced excitation reactions is shown in the figure. The electron beam generated by trochoidal electron monochromator collides perpendicularly with molecular beam created by effusive capillary. Photons emitted by excited reaction products are guided out of the vacuum chamber into the spectrometer equipped by CCD or photomultiplier detectors. The details of the experimental system were described elsewhere (Orszagh et al. 2017). Fig. 1. Scheme of the experimental apparatus. The main result of the experiment are the emission spectra (see second figure as an example) and emission cross section curves determined in precisely defined conditions. The shape of the spectrum given by population of individual ro-vibrational states and presence of individual spectral features depends strongly on excitation source and its energy (Bodewits et al. 2019). Fig. 2. Emission spectrum of hydrogen induced by 50 eV electrons impact. From the relative cross section curves it is possible to determine the lowest kinetic energy (threshold energy) of the electrons necessary to create a specific excited product and after calibration to absolute values we get series of cross sections for individual electron impact energies. Below the threshold the spectral line or band will not be present in the spectrum. The third figure shows relative emission cross section of oxygen where two thresholds correspond to different processes generating excited oxygen atom. Fig. 3. Relative emission cross section curve for oxygen atom. By putting spectral and relative cross sections information together by measuring emission spectra at several electron impact energies we can get a surface formed by emission cross section values within selected range of wavelengths and electron energies (see the last figure). Fig. 4. 3D spectral map of emission induced by electron impact on pyridine. Such complex information on the kinetics of electron-molecule excitation reactions is unfortunately rare mostly due to the technically demanding experiment but they provide valuable reference for analysis of astronomical spectroscopic data. Electron-molecule reactions in environment such as cometary comas or atmospheres can serve as a remote probe for characterization of local conditions if reference laboratory data on their kinetics are available. Acknowledgments This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 871149. This research was supported by Slovak grant agency VEGA within the projects nr. 1/0489/21 and 1/0553/22, by Slovak Research and Development Agency within project nr. APVV-19-0386, and by NASA ROSES project nr. NNH18ZDA001N-RDAP. Beenakker, C.I.M. et al., 1974. Chemical Physics, 6(3), pp.445–454. Bodewits, D. et al. 2016. Astronomical J., 152, 130. Bodewits, D. et al., 2019. Astrophys. J., 885(167), p. 16. Broiles, T.W. et al., 2016. Journal of Geophysical Research: Space Physics, 121 (8), pp.7407–7422. Farnham, T.L., Schleicher, D.G. & A'Hearn, M.F., 2000. Icarus, 147, p.180. Itikawa, Y. & Mason, N., 2005. Journal of Physical and Chemical Reference Data, 34(1), pp.1–22. Morgan, H. D. & Mentall, J.E., 1974. Journal of Chemical Physics, Volume 60, Issue 12, pp. 4734-4739. Orszagh, J. et al. 2017, Astrophys. J., 841(17).

Accurate measurement of the nuclear environment within a test tritium breeding-blanket module of a fusion reactor is crucial to determine tritium production rates which are relevant to self-suffici ...

In the last years, it has become apparent that detecting disruptions with sufficient anticipation time is an essential but not exclusive task of predictors. It is also important that the prediction is accompanied by appropriate qualifications of its reliability and it is formulated in mathematical terms appropriate for the task at hand (mitigation, avoidance, classification etc.). In this paper, a wide series of rule-based predictors, of the Classification and Regression Trees (CART) family, have been compared to assess their relative merits. An original refinement of the training, called noise-based ensembles, has allowed not only to obtain significantly better performance but also to increase the interpretability of the results. The final predictors can indeed be represented by a tree or a series of specific and clear rules. Such performance has been proved by analysing large databases of shots on JET with both the carbon wall and the ITER Like Wall. In terms of performance, the developed tools are therefore very competitive with other machine learning techniques, with the specificity of formulating the final models in terms of trees and simple rules.

Numerical scan at constant shows that core and scrape-off layer (SOL) radiations do not depend on the plasma current (I-p). Whereas the SOL radiation increases with seeding, the core radiation, how ...

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.

Activities under the EURO fusion work package JET3 programme have been established to enable the technological exploitation of the planned JET experiments over the next few years, which culminates in a D-T experimental campaign, DTE-2. In the areas of nuclear technology and nuclear safety the programme offers a unique opportunity to provide experimental data that is relevant to ITER. The key purpose of the collected data will be to support bench marking and validation activities relating to neutronics and activation codes, and associated nuclear data, that are used to predict the nuclear behavior of ITER component and materials, during and after operations. This paper details the status and key issues of the ongoing ACT sub-project under work package JET3, which aims to take advantage of the large 14 MeV neutron fluence expected during JET DTE-2 to irradiate samples of real ITER materials used in the manufacturing of the main in-vessel tokamak components. The materials considered, with specified minor elemental impurity levels, include: Nb3Sn, SS316L steels from a range of manufacturers, SS304B, Alloy 660, W, CuCrZr, OF-Cu, XM-19, Al bronze, NbTi and EUROFER. The activities include provision for measurement of nuclide activities for each material and comparison against the predicted quantities through calculation with the FISPACT-II inventory code. Included here are key pre-analysis results for the selected ITER irradiation samples, and corresponding optimization of dosimetry foils (Ti, Mn, Co, Ni, Y, Fe, Co, Sc, Ta) that will be irradiated at selected positions inside JET irradiation stations in order to determine the neutron spectrum. Preliminary experimental activation results through recent JET D–D operations are discussed.

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. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2017.03.031

Beryllium film deposition was studied with cavity samples in remote areas of the inner and outer JET divertor and below divertor tile 5 during the 2011-2012 campaign with the ITER-like wall. Predominantly beryllium films were formed inside the cavities with some additional carbon, the ratio Be/C was > 2. These deposited layers had high D/(Be+C) ratios of about 0.3. The formation of these films is mainly due to sticking of beryllium-containing particles with low sticking coefficients < 0.5. The observed surface loss probabilities depend on the position in the divertor. The particles responsible for film deposition originated from the location of in the divertor strike points. (C) 2016 Elsevier Ltd. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2016.12.005

Tungsten coatings deposited on carbon materials such as carbon fiber composite (CFC) or fine grain graphite are currently used in fusion devices as amour for plasma facing components (PFC). More than 4000 carbon tiles were W-coated by Combined Magnetron Sputtering and Ion Implantation technology for the ITER-like Wall at JET, ASDEX Upgrade and WEST tokamaks. The emissivity of W coatings is a key parameter required by protection systems of the W-coated PFC and also by the diagnostic tools in order to get correct values of temperature and heat loading. The emissivity of tungsten is rather well known, but the literature data refer to bulk tungsten or tungsten foils and not to coatings deposited on carbon materials. The emissivity was measured at the wavelengths of 1.064 μm, 1.75 μm, 3.75 μm and 4.0 μm. It was found that the structure of the substrate has a significant influence on the emissivity values. The temperature dependence of the emissivity in the range of 400 °C–1200 °C and the influence of the viewing angle were investigated as well. The results are given in a table for W coatings and for other materials of interest for fusion such as bulk W and bulk Be.

In the frame of the WPJET3-DT Technology project within the EUROfusion Consortium program, neutronics experiments are in preparation for the future deuterium–tritium campaign on JET (DTE2). The experiments will be conducted with the purpose to validate the neutronics codes and tools used in ITER, thus reducing the related uncertainties and the associated risks in the machine operation. This paper summarizes the status of previous shutdown dose rate benchmarks experiments and analyses performed at JET and focuses on the computational and experimental efforts conducted in preparation of the future DTE2 experiments. In particular, preliminary calculations and studies to select detectors and positions aimed to reduce uncertainties in the shutdown dose rate experiment are presented and discussed.

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.

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

Ionization of 2-nonanone, cyclopentanone, acetophenone, pyridine, and di- tert-butylpyridine (DTBP) in a corona discharge (CD) atmospheric pressure chemical ionization (APCI) ion source was studied using ion mobility (IMS) and time-of-flight mass spectrometry (TOF-MS). The IMS and MS spectra were recorded in the absence and presence of ammonia dopant. Without NH

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 MHD response and the penetration of a deuterium shattered pellet into a JET plasma is investigated via the non-linear reduced MHD code JOREK with the neutral gas shielding (NGS) ablation model. The dominant MHD destabilizing mechanism by the injection is identified as the local helical cooling at each rational surface, as opposed to the global current profile contraction. Thus the injected fragments destabilize each rational surface as they pass through them. The injection penetration is found to be much better compared to MGI, with the convective transport caused by core MHD instabilities (e.g. 1/1 kink) contributing significantly to the core penetration. Moreover, the injection with realistic JET SPI system configurations is simulated in order to provide some insights into future operations, and the impact on the total assimilation and penetration depth of varying injection parameters such as the injection velocity or fineness of shattering is assessed. Further, the effect of changing the target equilibrium temperature or q profile on the assimilation and penetration is also investigated. Such analysis will form the basis of further investigation into a desirable configuration for the future SPI system in ITER. EURATOM 633053

Tungsten (W) melting due to transient power loads, for example those delivered by edge localised modes (ELMs), is a major concern for next step fusion devices. A series of experiments has been performed on JET to investigate the dynamics of Type-I ELM-induced transient melting. Following initial exposures in 2013 of a W-lamella with sharp leading edge in the bulk W outer divertor, new experiments have been performed in 2016–2017 on a protruding W-lamella with a 15° slope, allowing direct and spatially resolved (0.85 mm/pixel) observation of the top surface using the IR thermography system viewing from the top of the poloidal cross-section. Thermal and IR analysis have already been conducted assuming the geometrical projection of the parallel heat flux on the W-lamellas, thus ignoring the gyro-radius orbit of plasma particles. Although it is well justified during L-mode or inter-ELM period, the hypothesis becomes questionable during ELM when the ion Larmor radius is larger. The goal of this paper is to extend the previous analysis based on the forward approach to the H-mode discharges and investigate in particular the gyro-radius effect during the Type-I ELMs, those used to achieve transient melting on the slope of the protruding W-lamella. Surface temperatures measured by the IR camera are compared with reconstructed synthetic data from 3D thermal modelling using heat loads derived from optical projection of the parallel heat flux (ignoring the gyro-radius orbit), 2D gyro-radius orbit and particle-in-cell (PIC) simulations describing the influence of finite Larmor-radius effects and electrical potential on the deposited power flux. Results show that the ELM power deposition behaves differently than the optical projection of the parallel heat flux, contrary to the L-mode observations, and may thus be due to the much larger gyro-orbits of the energetic ELM ions in comparison to L-mode or inter-ELM conditions. Keywords: ELM heat flux analysis, Optical projection, PIC modelling, Ion orbit modelling, IR data, Spatial resolution effect

Micro/nanoscopic observations on the surface of the divertor tiles used in the first campaign (2011-2012) of the JET tokamak with ITER-like Wall (JET ILW) have been carried out by means of several material analysis techniques. Previous results from the inner divertor were reported for a single poloidal section of the tile numbers 1, 3 and 4, i.e., upper, vertical and horizontal targets, respectively. The formation of the thick stratified mixed-material deposition layer on tiles 1 and 4, and erosion on tile 3 were identified. This study is mostly focused on the outer divertor: tiles 6, 7 and 8. In contrast to the inner tile, remarkable surface modifications have not been observed on the vertical target (tiles 7 and 8) where sputtering erosion and impurity deposition would have been almost balanced. Only a specific part of tile 6 (horizontal target) located near the exhaust channel was covered with a stratified ("geological-like") mixed-material deposition layer which mainly included Be and Ni with the thickness of similar to 2 mu m. Special feature of this mixed layer was that a certain amount of nitrogen (N) was clearly detected in the layer. Since the concentration of N varied with the depth position, it could be depended on the amount of that gas puffed for plasma edge cooling during the JET experimental campaign. In addition to the outer divertor tiles, a very interesting feature of the local erosion and deposition effects is reported in this paper.

There is experimental evidence that the pedestal dynamics in type-I ELMy H-mode discharges is significantly affected by a change in the recycling conditions at the tungsten plasma-facing components ...

In this paper we present the comparison of simulations with the numerical codes COREDIV and SOLPS5.0 for JET L-mode discharges with ITER like wall (ILW). The simulations have been performed for L-m ...

In this paper numerical simulations with the self-consistent COREDIV code of the planned JET DT experiments have been performed. First, record shot from the 1997 experiments was simulated and good ...

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 plas ...

A new Deuterium-Tritium (DT) campaign is planned at JET. An accurate calibration for the 14MeV neutron yield monitors is necessary. In order to perform the calibration a 14MeV Neutron Generator with suitable intensity (∼108n/s) will be used. Due to the intensity change during the Neutron Generator lifetime it would be necessary to monitor continuously the neutron emission intensity during the calibration using a compact detector attached to it. A high quality diamond detector has been chosen as one of the monitors. This detector has been fully characterized at the 14MeV Frascati Neutron Generator facility. The characterization procedure and the resulting 14MeV neutron response of the detector are described in this paper together with the obtained uncertainties.

The JET tokamak is operated with beryllium limiter tiles in the main chamber and tungsten coated carbon fiber composite tiles and solid W tiles in the divertor. One important issue is how wall materials are migrating during plasma operation. To study beryllium redistribution in the main chamber and in the divertor, a 10Be enriched limiter tile was installed prior to plasma operations in 2011-2012. Methods to take surface samples have been developed, an abrasive method for bulk Be tiles in the main chamber, which permits reuse of the tiles, and leaching with hot HCl to remove all Be deposited at W coated surfaces in the divertor. Quantitative analysis of the total amount of Be in cm2 sized samples was made with inductively coupled plasma atomic emission spectroscopy (ICP-AES). The 10Be/9Be ratio in the samples was measured with accelerator mass spectrometry (AMS). The experimental setup and methods are described in detail, including sample preparation, measures to eliminate contributions in AMS from the 10B isobar, possible activation due to plasma generated neutrons and effects of diffusive isotope mixing. For the first time marker concentrations are measured in the divertor deposits. They are in the range 0.4-1.2% of the source concentration, with moderate poloidal variation.

First wall materials in controlled fusion devices undergo serious modification by several physical and chemical processes arising from plasma-wall interactions. Detailed information is required for ...

International audience; 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 non-linear 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 ITG 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 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

The EUROfusion Work Package JET3 programme, established to enable the technological exploitation of the JET experiments over the next years, includes, within the NEXP subproject, a novel Shutdown Dose Rate (SDR) experiment. Considering its ITER-relevance, SDR experiments at JET represent a unique opportunity to validate the numerical tools for ITER nuclear analysis, through the comparison between numerical predictions and measured quantities (C/E). Within this framework, two active gamma dosimeters based on spherical air-vented ionization chambers (ICs) have been installed in ex-vessel positions close to the horizontal ports of the tokamak in Octants 1 and 2. The first JET campaign exploited in the novel SDR experiment is the latest 5-week Deuterium-Deuterium campaign (c36b), which achieved the best results in recent years in terms of high power operation. The present work is dedicated to the analysis of dose rate measurements carried out during this campaign and after shutdown. Proper correction factors are evaluated and applied to the instrument reading, while influence quantities and error sources are analyzed in order to calculate the overall experimental uncertainty.

Plasma tomography consists of reconstructing a two-dimensional radiation profile of a poloidal cross section of a fusion device based on line-integrated measurements along several lines of sight. The reconstruction process is computationally intensive, and in practice, only a few reconstructions are usually computed per pulse. In this work, we trained a deep neural network based on a large collection of sample tomograms that have been produced at JET over several years. Once trained, the network is able to reproduce those results with high accuracy. More importantly, it can compute all the tomographic reconstructions for a given pulse in just a few seconds. This makes it possible to visualize several phenomena—such as plasma heating, disruptions, and impurity transport—over the course of the entire pulse.

Beryllium oxide (BeO) and deuteroxide (BeOxDy) have been found on the melted zone of a beryllium tile extracted from the upper dump plate of JET-ILW (2011-2012 campaign). Results have been obtained using Raman microscopy, which is sensitive to both the chemical bond and crystal structure, with a micrometric lateral resolution. BeO is found with a wurtzite crystal structure. BeOxDy is found as three different types which are not the beta-phase but behaves as molecular species like Be(OD)(2), O(Be-D)(2) and DBeOD. The presence of a small amount of trapped D2O is also suspected. Our results therefore strongly suggest that D trapping occurs after melting through the formation of deuteroxides. The temperature increase favors the formation of crystal BeO which favors deuterium trapping through OD bonding. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2018.11.008

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

Understanding the details of the correlation between time series is an essential step on the route to assessing the causal relation between systems. Traditional statistical indicators, such as the Pearson correlation coefficient and the mutual information, have some significant limitations. More recently, transfer entropy has been proposed as a powerful tool to understand the flow of information between signals. In this paper, the comparative advantages of transfer entropy, for determining the time horizon of causal influence, are illustrated with the help of synthetic data. The technique has been specifically revised for the analysis of synchronization experiments. The investigation of experimental data from thermonuclear plasma diagnostics proves the potential and limitations of the developed approach. For complete list of authors see http://dx.doi.org/10.3390/e20090627

This paper deals with a divertor coil currents optimized procedure to design High Flux Expansion (HFE) configurations in the JET tokamak aimed to study the effects of flux expansion variation on the radiation fraction and radiated power re-distribution. A number of benefits of HFE configuration have been experimentally demonstrated on TCV, EAST, NSTX and DIII-D tokamaks and are under investigation for next generation devices, as DEMO and DTT. The procedure proposed here exploits the linearized relation between the plasma-wall gaps and the Poloidal Field (PF) coil currents. Once the linearized model is provided by means of CREATE-NL code, the divertor coils currents are calculated using a constrained quadratic programming optimization procedure, in order to achieve HFE configuration. Flux expanded configurations have been experimentally realized both in ohmic and heated plasma with and without nitrogen seeding. Preliminary results on the effects of the flux expansion variation on total power radiation increase will be also briefly discussed. © 2018 Elsevier B.V.

This paper analyses the final phase of hybrid scenario discharges at JET, the reduction of auxiliary heating towards finally the Ohmic phase. The here considered Ohmic phase is mostly still in the current flattop but may also be in the current ramp down. For this purpose a database is created of 54 parameters in 7 phases distributed in time of the discharge. It is found that the occurrence of a locked mode is in most cases preceded by a radiation peaking after the main heating phase either in a low power phase and/or in the Ohmic phase. To gain insight on the importance of different parameters in this process a correlation analysis to the radiation peaking in the Ohmic phase is done. The first finding is that the further away in time the analysed phases are the less the correlation is. This means in the end that a good termination scenario might also be able to terminate unhealthy plasmas safely. The second finding is that remaining impurities in the plasma after reducing the heating power in the termination phase are the most important reason for generating a locked mode which can lead to a disruption. © EURATOM 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

A novel Shutdown dose rate benchmark experiment has been performed at Joint European Torus (JET) machine during the last high performance Deuterium-Deuterium (DD) campaign in preparation of future ...

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.

Fusion energy systems pose unique challenges to the modeling and simulation community. These challenges must be met to ensure the success of the ITER experimental fusion reactor. ITER’s complex systems require detailed modeling that goes beyond the scale of comparable simulations to date. In this work, the Denovo radiation transport code was used to calculate neutron fluence and kerma for the JET streaming benchmark. This work was performed on the Titan supercomputer at the Oak Ridge Leadership Computing Facility. Denovo is a novel three-dimensional discrete ordinates transport code designed to be highly scalable. Sensitivity studies have been completed to examine the impact of several deterministic parameters. Results were compared against experiment as well as the MCNP and Shift Monte Carlo codes.

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.

The total emission of radiation is a crucial quantity to calculate the power balances and to understand the physics of any Tokamak. Bolometric systems are the main tool to measure this important physical quantity through quite sophisticated tomographic inversion methods. On the Joint European Torus, the coverage of the bolometric diagnostic, due to the availability of basically only two projection angles, is quite limited, rendering the inversion a very ill-posed mathematical problem. A new approach, based on the maximum likelihood, has therefore been developed and implemented to alleviate one of the major weaknesses of traditional tomographic techniques: the difficulty to determine routinely the confidence intervals in the results. The method has been validated by numerical simulations with phantoms to assess the quality of the results and to optimise the configuration of the parameters for the main types of emissivity encountered experimentally. The typical levels of statistical errors, which may significantly influence the quality of the reconstructions, have been identified. The systematic tests with phantoms indicate that the errors in the reconstructions are quite limited and their effect on the total radiated power remains well below 10%. A comparison with other approaches to the inversion and to the regularization has also been performed.

Diagnostic alkali atom (e.g., lithium) beams are routinely used to diagnose magnetically confined plasmas, namely, to measure the plasma electron density profile in the edge and the scrape off layer region. A light splitting optics system was installed into the observation system of the lithium beam emission spectroscopy diagnostic at the Joint European Torus (JET) tokamak, which allows simultaneous measurement of the beam light emission with a spectrometer and a fast avalanche photodiode (APD) camera. The spectrometer measurement allows density profile reconstruction with ∼10 ms time resolution, absolute position calculation from the Doppler shift, spectral background subtraction as well as relative intensity calibration of the channels for each discharge. The APD system is capable of measuring light intensities on the microsecond time scale. However ∼100 μs integration is needed to have an acceptable signal to noise ratio due to moderate light levels. Fast modulation of the beam up to 30 kHz is implemented which allows background subtraction on the 100 μs time scale. The measurement covers the 0.9 < ρpol < 1.1 range with 6-10 mm optical resolution at the measurement location which translates to 3-5 mm radial resolution at the midplane due to flux expansion. An automated routine has been developed which performs the background subtraction, the relative calibration, and the comprehensive error calculation, runs a Bayesian density reconstruction code, and loads results to the JET database. The paper demonstrates the capability of the APD system by analyzing fast phenomena like pellet injection and edge localized modes.

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 + ...

In this paper, the dynamic characteristics of type-I ELM time-series from the JET tokamak, the world's largest magnetic confinement plasma physics experiment, have been investigated. The dynamic analysis has been focused on the detection of nonlinear structure in D a radiation time series. Firstly, the method of surrogate data has been applied to evaluate the statistical significance of the null hypothesis of static nonlinear distortion of an underlying Gaussian linear process. Several nonlinear statistics have been evaluated, such us the time delayed mutual information, the correlation dimension and the maximal Lyapunov exponent. The obtained results allow us to reject the null hypothesis, giving evidence of underlying nonlinear dynamics. Moreover, no evidence of low-dimensional chaos has been found; indeed, the analysed time series are better characterized by the power law sensitivity to initial conditions which can suggest a motion at the 'edge of chaos', at the border between chaotic and regular non-chaotic dynamics. This uncertainty makes it necessary to further investigate about the nature of the nonlinear dynamics. For this purpose, a second surrogate test to distinguish chaotic orbits from pseudoperiodic orbits has been applied. In this case, we cannot reject the null hypothesis which means that the ELM time series is possibly pseudo-periodic. In order to reproduce pseudo-periodic dynamical properties, a periodic state-of-the-art model, proposed to reproduce the ELM cycle, has been corrupted by a dynamical noise, obtaining time series qualitatively in agreement with experimental time series.

The impact of the finite grid size in SOLPS-ITER edge plasma simulations is assessed for JET H-mode discharges with a metal wall. For a semi-horizontal divertor configuration it is shown that the separatrix density is at least 30% higher when a narrow scrape-off layer (SOL) grid width is chosen in SOLPS-ITER compared to the case for which the SOL grid width is maximised. The density increase is caused by kinetic neutrals being not confined inside the divertor region because of the reduced extent of the plasma grid. In this case, an enhanced level of reflections of energetic neutrals at the low-field side (LFS) metal divertor wall is observed. This leads to a shift of the ionisation source further upstream which must be accounted for as a numerical artefact. An overestimate in the cooling at the divertor entrance is observed in this case, identified by a reduced heat flux decay parameters lambda(div)(q). Otherwise and further upstream the mid-plane heat decay length lambda(q) parameter is not affected by any change in divertor dissipation. This confirms the assumptions made for the ITER divertor design studies, i.e. that lambda(q) upstream is essentially set by the assumptions for the ratio radial to parallel heat conductivity. It is also shown that even for attached conditions the decay length relations lambda(ne)>lambda(Te)>lambda(q) hold in the near-SOL upstream. Thus for interpretative edge plasma simulations one must take the (experimental) value of lambda(ne) into account, rather than lambda(q), as the former actually defines the required minimum upstream SOL grid extent. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2018.10.013

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.