Professor Felix Parra Diaz

Abstract:

Conventional radially-local neoclassical calculations become inadequate if the radial gradient scale lengths of the H-mode pedestal become as small as the poloidal ion gyroradius. Here, we describe a radially global $\delta f$ continuum code that generalizes neoclassical calculations to allow stronger gradients. As with conventional neoclassical calculations, the formulation is time-independent and requires only the solution of a single sparse linear system. We demonstrate precise agreement with an asymptotic analytic solution of the radially global kinetic equation in the appropriate limits of aspect ratio and collisionality. This agreement depends crucially on accurate treatment of finite orbit width effects.

Abstract:

New diagnostic, modelling and plant capability on the Mega Ampère Spherical Tokamak (MAST) have delivered important results in key areas for ITER/DEMO and the upcoming MAST Upgrade, a step towards future ST devices on the path to fusion currently under procurement. Micro-stability analysis of the pedestal highlights the potential roles of micro-tearing modes and kinetic ballooning modes for the pedestal formation. Mitigation of edge localized modes (ELM) using resonant magnetic perturbation has been demonstrated for toroidal mode numbers n = 3, 4, 6 with an ELM frequency increase by up to a factor of 9, compatible with pellet fuelling. The peak heat flux of mitigated and natural ELMs follows the same linear trend with ELM energy loss and the first ELM-resolved Ti measurements in the divertor region are shown. Measurements of flow shear and turbulence dynamics during L-H transitions show filaments erupting from the plasma edge whilst the full flow shear is still present. Off-axis neutral beam injection helps to strongly reduce the redistribution of fast-ions due to fishbone modes when compared to on-axis injection. Low-k ion-scale turbulence has been measured in L-mode and compared to global gyro-kinetic simulations. A statistical analysis of principal turbulence time scales shows them to be of comparable magnitude and reasonably correlated with turbulence decorrelation time. Te inside the island of a neoclassical tearing mode allow the analysis of the island evolution without assuming specific models for the heat flux. Other results include the discrepancy of the current profile evolution during the current ramp-up with solutions of the poloidal field diffusion equation, studies of the anomalous Doppler resonance compressional Alfvén eigenmodes, disruption mitigation studies and modelling of the new divertor design for MAST Upgrade. The novel 3D electron Bernstein synthetic imaging shows promising first data sensitive to the edge current profile and flows. © 2013 IAEA, Vienna.

Abstract:

Gyrokinetic simulations have greatly improved our theoretical understanding of turbulent transport in fusion devices. Most gyrokinetic models in use are delta-f simulations in which the slowly varying radial profiles of density and temperature are assumed to be constant for turbulence saturation times, and only the turbulent electromagnetic fluctuations are calculated. New massive simulations are being built to self-consistently determine the radial profiles of density and temperature. However, these new codes have failed to realize that modern gyrokinetic formulations, composed of a gyrokinetic Fokker-Planck equation and a gyrokinetic quasineutrality equation, are only valid for delta-f simulations that do not reach the longer transport time scales necessary to evolve radial profiles. In tokamaks, due to axisymmetry, the evolution of the axisymmetric radial electric field is a challenging problem requiring substantial modifications to gyrokinetic treatments. In this thesis, I study the effect of turbulence on the global electric field and plasma flows. By studying the current conservation equation, or vorticity equation, I prove that the long wavelength, axisymmetric flow must remain neoclassical and I show that the tokamak is intrinsically ambipolar, i.e., the radial current is zero to a very high order for any long wavelength radial electric field. Intrinsic ambipolarity is the origin of the problems with the modern gyrokinetic approach since the lower order gyrokinetic quasineutrality (if properly evaluated) is effectively independent of the radial electric field. I propose a new gyrokinetic formalism to solve for the global radial electric field.

Abstract:

The induced electric field in a tokamak drives a parallel electron current flow. In an inhomogeneous, finite beta plasma, when this electron flow is comparable to the ion thermal speed, the Alfven mode wave solutions of the electromagnetic gyrokinetic equation can become nearly purely growing kink modes. Using the new "low-flow" version of the gyrokinetic code GS2 developed for momentum transport studies [Barnes et al 2013 Phys. Rev. Lett. 111, 055005], we are able to model the effect of the induced parallel electric field on the electron distribution to study the destabilizing influence of current on stability. We identify high mode number kink modes in GS2 simulations and make comparisons to analytical theory in sheared magnetic geometry. We demonstrate reassuring agreement with analytical results both in terms of parametric dependences of mode frequencies and growth rates, and regarding the radial mode structure.

Abstract:

A local magnetic equilibrium solution is sought around the magnetic axis in order to identify the key parameters defining the magnetic-surface's up-down asymmetry in the core of tokamak plasmas. The asymmetry is found to be determined essentially by the ratio of the toroidal current density flowing on axis to the fraction of the external field's odd perturbation that manages to propagate from the plasma boundary into the core. The predictions are tested and illustrated first with an analytical Solovev equilibrium and then using experimentally relevant numerical equilibria. Hollow current-density distributions, and hence reverse magnetic shear, are seen to be crucial to bring into the core asymmetry values that are usually found only near the plasma edge.