Core-edge 2D fluid modeling of full tokamak discharge with varying magnetic equilibrium: from WEST start-up to ramp-down

In the present work we investigate for the first time the 2D fluid transport of the plasma in WEST during an entire discharge from the start-up to the ramp-down (shot #54487). The evolution of density profile, electron and ion temperatures together with the experimental magnetic equilibrium, total current and gas-puff rate is investigated. Comparisons with the interferometry diagnostic show a remarkable overall qualitative agreement during the discharge that can be quantitative at some locations in the plasma core. If at the onset of the X-points during the ramp-up the electron heat flux is dominant at the target, present results show that the ion heat flux becomes dominant during the stationary phase of the discharge. Using a simple model for erosion, present results assess the tungsten sputtering due to deuterium ions during the start-up and ramp-up phases of the discharge and confirms the need to consider full discharge simulation to accurately treat the W source of contamination. This work also demonstrates the interest of developing magnetic equilibrium free solver including efficient time integration to step toward predictive capabilities in the future for fusion operation.

The model of synchronization between internal reconnections and edge-localized modes

A new model for interaction between the internal reconnections caused by sawtooth and the edge-localized modes (ELM) was presented. The experimental evidence of the coupling between sawtooth crash and ELM events were observed in the Globus-M and Globus-M2 tokamaks. The numerical analysis of magnetic equilibrium showed that internal reconnections can induce the excess current density near the separatrix during the several hundreds of μs. The excess current destabilizes the peeling-ballooning instability. The peeling-ballooning stability analysis showed that the penetration depth of the induced current should be in the range of ψnorm = 0.8 - 0.95 to trigger the instability.

Implementing an X-ray tomography method for fusion devices

In fusion devices, the X-ray plasma emissivity contains essential information on the magnetohydrodynamic activity, the magnetic equilibrium and on the transport of impurities, in particular for tokamaks in the soft X-ray (SXR) energy range of 0.1–20 keV. In this context, tomography diagnostics are a key method to estimate the local plasma emissivity from a given set of line-integrated measurements. Unfortunately, the reconstruction problem is mathematically ill-posed, due to very sparse and noisy measurements, requiring an adequate regularization procedure. The goal of this paper is to introduce, with a didactic approach, some methodology and tools to develop an X-ray tomography algorithm. Based on a simple 1D tomography problem, the Tikhonov regularization is described in detail with a study of the optimal reconstruction parameters, such as the choice of the emissivity spatial resolution and the regularization parameter. A methodology is proposed to perform an in situ sensitivity and position cross-calibration of the detectors with an iterative procedure, by using the information redundancy and data variability in a given set of reconstructed profiles. Finally, the basic steps to build a synthetic tomography diagnostics in a more realistic tokamak environment are introduced, together with some tools to assess the capabilities of the 2D tomography algorithm.

Perturbing an axisymmetric magnetic equilibrium to obtain a quasi-axisymmetric stellarator

It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria are discussed.

Anomalous losses of energetic particles in the presence of an oscillating radial electric field in fusion plasmas

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. It is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

DRIFT-ALFVEN INSTABILITY IN A 2D MAGNETOTAIL CONFIGURATION – THE ADDITION OF BOUNCING ELECTRONS

To explain the possible destabilization of a 2D magnetic equilibrium such as the Near-Earth magnetotail, we developed a kinetic model describing the resonant interaction of electromagnetic fluctuations and bouncing electrons trapped in the magnetosphere, characterized by a high plasma density gradient. A small-β approximation is used in agreement with a small field line curvature. It has been found that for a quasi-dipole configuration, unstable electromagnetic modes may develop in the current sheet in westward direction with a growth rate of the order of a few tenth of seconds provided that the typical scale of density gradient slope responsible for the diamagnetic drift effects is over one Earth radius or less. This instability growth rate is large enough to destabilise the current sheet on time scales often observed during substorm onset.