Investigation into the formation of the scrape-off layer density shoulder in JET ITER-like wall L-mode and H-mode plasmas

The low temperature boundary layer plasma (Scrape-Off-Layer or SOL) between the hot core and the surrounding vessel determines the level of power-loading, erosion and implantation of material surfaces, and thus the viability of tokamak-based fusion as an energy source. This study explores mechanisms affecting the formation of flattened density profiles, so-called ‘density shoulders’, in the low-field side (LFS) SOL, which modify ion and neutral fluxes to surfaces – and subsequent erosion. There is evidence against local enhancement of ionization inducing shoulder formation. We find that increases in SOL parallel resistivity, Λdiv (=[L||νei Ωi ]/cs Ωe), postulated to lead to shoulder growth through changes in SOL turbulence characteristics, correlates with increases in upstream SOL shoulder amplitude, As only under a subset of conditions (D2-fuelled L-mode density scans with outer strike point on the horizontal target). Λdiv fails to correlate with As for cases of N2 seeding or during sweeping of the strike point across the horizontal target. The limited correlation of Λdiv with As was also found for H-mode discharges. Thus, while Λdiv above a threshold of ~1 may be necessary for shoulder formation and/or growth, another shoulder mechanism is required. More significantly we find that in contrast to parallel resistivity, outer divertor recycling as quantified by the total outer divertor Balmer Dα emission, I-Dα, does scale with shoulder amplitude where Λdiv does and even where Λdiv fails. Divertor recycling could lead to SOL density shoulder formation through: a) reducing the parallel to the field flow (loss) of ions out of the SOL to the divertor; and b) changes in radial electric fields which lead to ExB poloidal flows as well as potentially affecting the SOL turbulence birth characteristics. Thus changes in divertor recycling may be the sole process in bringing about SOL density shoulders or in tandem with parallel resistivity

Neutral pathways and heat flux widths in vertical- and horizontal-target EDGE2D-EIRENE simulations of JET

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

Observations and modelling of ion cyclotron emission observed in JET plasmas using a sub-harmonic arc detection system during ion cyclotron resonance heating

Peer reviewe

Effects of magnetic ripple on 3D equilibrium and alpha particle confinement in the European DEMO

Abstract An assessment of alpha particle confinement is performed in the European DEMO reference design. 3D MHD equilibria with nested flux-surfaces and single magnetic axis are obtained with the VMEC free-boundary code, thereby including the plasma response to the magnetic ripple created by the finite number of TF coils. Populations of fusion alphas that are consistent with the equilibrium profiles are evolved until slowing-down with the VENUS-LEVIS orbit code in the guiding-centre approximation. Fast ion losses through the last-closed flux-surface are numerically evaluated with two ripple models: (1) using the 3D equilibrium and (2) algebraically adding the non-axisymmetric ripple perturbation to the 2D equilibrium. By virtue of the small ripple field and its non-resonant nature, both models quantitatively agree. Differences are however noted in the toroidal location of particles losses on the last-closed flux-surface, which in the first case is 3D and in the second not. Superbanana transport, i.e. ripple-well trapping and separatrix crossing, is expected to be the dominant loss mechanism, the strongest effect on alphas being between 100–200 KeV. Above this, stochastic ripple diffusion is responsible for a rather weak loss rate, as the stochastisation threshold is observed numerically to be higher than analytic estimates. The level of ripple in the current 18 TF coil design of the European DEMO is not found to be detrimental to fusion alpha confinement.</jats:p