David Pfefferlé

NBI fast ion confinement in the helical core of MAST hybrid-like plasmas

Energetic ions are found to be transported strongly from the core of MAST hybrid-like plasmas during long-lived mode (LLM) magnetohydrodynamic activity. The resulting impact on the neutral beam ion deposition and concurrent current drive is modelled using the guiding-centre approximation in the internal kinked magnetic topology. General coordinate guiding-centre equations are extended for this purpose. It is found that the kinked core spirals around the position of strongest ionization, which remains geometrically centred, so that a large fraction of the population is deposited in the high shear external region where the plasma is almost axisymmetric. Those particles ionized in the low shear region exhibit exotic drift motion due to the strongly non-axisymmetric equilibrium, periodically passing near the magnetic axis and then reflected by the boundary of the kinked equilibrium, which in this respect acts as a confining pinch. Broad agreement is found against experimental measurement of fast ion particle confinement degradation as the MAST LLM amplitude varies.

Impact of RMP magnetic field simulation models on fast ion losses

Two opposing approaches to include resonant magnetic perturbations (RMPs) in fast ion simulations are compared, one where the vacuum field caused by the RMP current coils is added to the axisymmetric MHD equilibrium, the other where the MHD equilibrium includes the plasma response within the 3D deformation of its flux-surfaces. The first model admits large regions of stochastic field-lines that penetrate the plasma without alteration. The second assumes nested flux-surfaces with a single magnetic axis, which excludes stochastic field-lines, and embeds the RMPs within a 3D saturated ideal MHD state. The two descriptions of RMPs have been implemented in the VENUS-LEVIS guiding-centre orbit code. Simulations of fast ion populations resulting from MAST neutral beam injection have been applied to MAST n = 3 RMP coil configuration. At low beam energies, particle losses are dominated by parallel transport due to the stochasticity of the field-lines (vacuum-RMP model), whereas at higher energies, losses are accredited to the 3D structure of the perturbed plasma and the resulting drifts (equilibrium-RMP model).

Hybrid guiding-centre/full-orbit simulations in non-axisymmetric magnetic geometry exploiting general criterion for guiding-centre accuracy

To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details are provided on how to implement it in cylindrical coordinates and in flux coordinates that rely on the geometric toroidal angle. A means of switching between guiding-centre and full-orbit equations at first order in Larmor radius with minimal discrepancy is shown. Techniques are applied to a MAST (mega amp spherical tokamak) helical core equilibrium in which the inner kinked flux-surfaces are tightly compressed against the outer axisymmetric mantle and where the parallel current peaks at the nearly rational surface. This is put in relation with the simpler situation B(x, y, z) = B0[sin(kx)ey +cos(kx)ez], for which full orbits and lowest order drifts are obtained analytically. In the kinked equilibrium, the full orbits of NBI fast ions are solved numerically and shown to follow helical drift surfaces. This result partially explains the off-axis redistribution of neutral beam injection fast particles in the presence of MAST long-lived modes (LLM).

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

Curvilinear guiding-centre drift and full-orbit equations of motion are presented as implemented in the code{VENUS-LEVIS} code. A dedicated interpolation scheme based on Fourier reconstruction in the toroidal and poloidal direction and cubic spline in the radial direction of flux coordinate systems is detailed. This interpolation method exactly preserves the order of the RK4 integrating scheme which is crucial for the investigation of fast particle trajectories in 3D magnetic structures such as helical saturated tokamak plasma states, stellarator geometry and resonant magnetic perturbations (RMP). The initialisation of particles with respect to the guiding-centre is discussed. Two approaches to implement RMPs in orbit simulations are presented, one where the vacuum field is added to the 2D equilibrium, creating islands and stochastic regions, the other considering 3D nested flux-surfaces equilibrium including the RMPs.

Exact collisional moments for plasma fluid theories

The velocity-space moments of the often troublesome nonlinear Landau collision operator are expressed exactly in terms of multi-index Hermite-polynomial moments of distribution functions. The collisional moments are shown to be generated by derivatives of two well-known functions, namely, the Rosenbluth-MacDonald-Judd-Trubnikov potentials for a Gaussian distribution. The resulting formula has a nonlinear dependency on the relative mean flow of the colliding species normalised to the root-mean-square of the corresponding thermal velocities and a bilinear dependency on densities and higher-order velocity moments of the distribution functions, with no restriction on temperature, flow, or mass ratio of the species. The result can be applied to both the classic transport theory of plasmas that relies on the Chapman-Enskog method, as well as to derive collisional fluid equations that follow Grads moment approach. As an illustrative example, we provide the collisional ten-moment equations with exact conservatio...

Fluid moments of the nonlinear Landau collision operator

An important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multivariate Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The proposed methodology provides a systematic and rigorous means of extending the validity of fluid models that have an underlying inverse-square force particle dynamics to arbitrary collisionality and flow.

Nonlinear resistivity for magnetohydrodynamical models

A new formulation of the plasma resistivity that stems from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The full expression for the collisional momentum-transfer rate is determined and is used to formulate the nonlinear anisotropic resistivity. The new formalism recovers the Spitzer resistivity, as well as the concept of thermal force if the heat flux is assumed to be proportional to a temperature gradient. Furthermore, if the pressure tensor is related to viscous stress, the latter enters the expression for the resistivity. The relative importance of the nonlinear term(s) with respect to the well-established electron inertia and Hall terms is also examined. The subtle implications of the nonlinear resistivity, and its dependence on the fluid vari...

Applications of the SCENIC code package to the minority ion-cyclotron heating in Wendelstein 7-X plasmas

We present SCENIC simulations of a W7X 4He plasma with 1% H minority and with an antenna model close to the design foreseen for the W7X ICRF antenna [1, 2]. A high mirror and a standard equilibrium are considered. The injected wave frequency is fixed at 33.8 MHz and 39.6MHz respectively and only fundamental minority heating is considered. Included in this calculation is a new realistic model of the antenna, where it is found that the localization of the antenna geometry tends to break the five-fold periodicity of the system. We assess the heat transfer through the toroidal periods via Coulomb collisions.

An approximate single fluid 3-dimensional magnetohydrodynamic equilibrium model with toroidal flow

An approximate model for a single fluid three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium with pure isothermal toroidal flow with imposed nested magnetic flux surfaces is proposed. It recovers the rigorous toroidal rotation equilibrium description in the axisymmetric limit. The approximation is valid under conditions of nearly rigid or vanishing toroidal rotation in regions with significant 3D deformation of the equilibrium flux surfaces. Bifurcated helical core equilibrium simulations of long-lived modes in the MAST device demonstrate that the magnetic structure is only weakly affected by the flow but that the 3D pressure distortion is important. The pressure is displaced away from the major axis and therefore is not as noticeably helically deformed as the toroidal magnetic flux under the subsonic flow conditions measured in the experiment. The model invoked fails to predict any significant screening by toroidal plasma rotation of resonant magnetic perturbations in MAST free boundary computations.

Non-planar elasticae as optimal curves for the magnetic axis of stellarators

The Euler-Lagrange equations are derived for finite length three-dimensional curves that optimize their bending energy while yielding fixed integrated torsion. The obvious translational and rotational symmetry is exploited to express solutions in a preferred cylindrical coordinate system in terms of elliptic Jacobi functions. These solution curves, which, up to similarity transformations, depend on three dimensionless parameters, do not necessarily close. Two closure conditions are obtained for the vertical and toroidal displacement (the radial coordinate being trivially periodic) to yield an infinite countable set of one-parameter families of closed non-planar curves. The behaviour of the integrated torsion (Twist), the Linking of the Frenet frame and the Writhe of the solution curves is studied in light of the Calugareanu theorem. A refreshed interpretation of Merciers formula for the on-axis rotational transform of stellarator magnetic field-lines is proposed.