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.
A neutron spectrometer at ITER is expected to provide estimates of plasma parameters such as ion temperature, Ti, fuel ion ratio, nt/nd, and Qthermal/Qtot, with 10-20% precision at a time resolution, Δt, of at least 100 ms. The present paper describes a method for evaluating different neutron spectroscopy techniques based on their instrumental response functions and synthetic measurement data. We include five different neutron spectrometric techniques with realistic response functions, based on simulations and measurements where available. The techniques are magnetic proton recoil, thin-foil proton recoil, gamma discriminating organic scintillator, diamond and time-of-flight. The reference position and line of sight of a high resolution neutron spectrometer on ITER are used in the study. ITER plasma conditions are simulated for realistic operating scenarios. The ITER conditions evaluated are beam and radio frequency heated and thermal deuterium-tritium plasmas. Results are given for each technique in terms of the estimated time resolution at which the parameter determination can be made within the required precision (here 10% for Ti and the relative intensities of NB and RF emission components). It is shown that under the assumptions made, the thin-foil techniques out-perform the other spectroscopy techniques in practically all measurement situations. For thermal conditions, the range of achieved Δt in the determination of Ti varies in time scales from ms (for the magnetic and thin-foil proton recoil) to s (for gamma discriminating organic scintillator).
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 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.
The plasma diagnostic and control (D&C) system for a future tokamak demonstration fusion reactor (DEMO) will have to provide reliable operation near technical and physics limits, while its front-end components will be subject to strong adverse effects within the nuclear and high temperature plasma environment. The ongoing developments for the ITER D&C system represent an important starting point for progressing towards DEMO. Requirements for detailed exploration of physics are however pushing the ITER diagnostic design towards using sophisticated methods and aiming for large spatial coverage and high signal intensities, so that many front-end components have to be mounted in forward positions. In many cases this results in a rapid aging of diagnostic components, so that additional measures like protection shutters, plasma based mirror cleaning or modular approaches for frequent maintenance and exchange are being developed. Under the even stronger fluences of plasma particles, neutron/gamma and radiation loads on DEMO, durable and reliable signals for plasma control can only be obtained by selecting diagnostic methods with regard to their robustness, and retracting vulnerable front-end components into protected locations. Based on this approach, an initial DEMO D&C concept is presented, which covers all major control issues by signals to be derived from at least two different diagnostic methods (risk mitigation).
An initial concept for the plasma diagnostic and control (D&C) system has been developed as part of European studies towards the development of a demonstration tokamak fusion reactor (DEMO). The main objective is to develop a feasible, integrated concept design of the DEMO D&C system that can provide reliable plasma control and high performance (electricity output) over extended periods of operation. While the fusion power is maximized when operating near to the operational limits of the tokamak, the reliability of operation typically improves when choosing parameters significantly distant from these limits. In addition to these conflicting requirements, the D&C development has to cope with strong adverse effects acting on all in vessel components on DEMO (harsh neutron environment, particle fluxes, temperatures, electromagnetic forces, etc.). Moreover, space allocation and plasma access are constrained by the needs for first wall integrity and optimization of tritium breeding. Taking into account these boundary conditions, the main DEMO plasma control issues have been formulated, and a list of diagnostic systems and channels needed for plasma control has been developed, which were selected for their robustness and the required coverage of control issues. For a validation and refinement of this concept, simulation tools are being refined and applied for equilibrium, kinetic and mode control studies.
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 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.
In ITER and DEMO, various control objectives related to plasma control must be simultaneously achieved by the plasma control system (PCS), in both normal operation as well as off-normal conditions. The PCS must act on off-normal events and deviations from the target scenario, since certain sequences (chains) of events can precede disruptions. It is important that these decisions are made while maintaining a coherent prioritization between the real-time control tasks to ensure high-performance operation. In this paper, a generic architecture for task-based integrated plasma control is proposed. The architecture is characterized by the separation of state estimation, event detection, decisions and task execution among different algorithms, with standardized signal interfaces. Central to the architecture are a plasma state monitor and supervisory controller. In the plasma state monitor, discrete events in the continuous-valued plasma state arc modeled using finite state machines. This provides a high-level representation of the plasma state. The supervisory controller coordinates the execution of multiple plasma control tasks by assigning task priorities, based on the finite states of the plasma and the pulse schedule. These algorithms were implemented on the TCV digital control system and integrated with actuator resource management and existing state estimation algorithms and controllers. The plasma state monitor on TCV can track a multitude of plasma events, related to plasma current, rotating and locked neoclassical tearing modes, and position displacements. In TCV experiments on simultaneous control of plasma pressure, safety factor profile and NTMs using electron cyclotron heating (ECI I) and current drive (ECCD), the supervisory controller assigns priorities to the relevant control tasks. The tasks are then executed by feedback controllers and actuator allocation management. This work forms a significant step forward in the ongoing integration of control capabilities in experiments on TCV, in support of tokamak reactor operation.
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.
For the first time, over five confinement times, the self-consistent flux driven time evolution of heat, momentum transport and particle fluxes of electrons and multiple ions including Tungsten (W) is modeled within the integrated modeling platform JFTTO (Romanelli et al 2014 Plasma Fusion Res. 9 1-4), using first principle-based codes: namely, QuaLiKiz (Bourdelle et al 2016 Plasma Phys. Control. Fusion 58 014036) for turbulent transport and NEO (Belli and Candy 2008 Plasma Phys. Control. Fusion 50 95010) for neoclassical transport. For a JET-ILW pulse, the evolution of measured temperatures, rotation and density profiles are successfully predicted and the observed W central core accumulation is obtained. The poloidal asymmetries of the W density modifying its neoclassical and turbulent transport are accounted for. Actuators of the W core accumulation are studied: removing the central particle source annihilates the central W accumulation whereas the suppression of the torque reduces significantly the W central accumulation. Finally, the presence of W slightly reduces main ion heat turbulent transport through complex nonlinear interplays involving radiation, effective charge impact on ITG and collisionality.
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-Schluter 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.