Professor Felix Parra Diaz

Abstract:

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

Abstract:

Ballooning mode saturation is investigated in realistic stellarator configurations using the flux tube approach of Ham et al (2018 Plasma Phys. Control. Fusion 60 075017), Ham et al (2016 Phys. Rev. Lett. 116 235001). The method is adapted to account for the lack of exact force balance in stellarator equilibrium solvers that assume existence of nested flux surfaces. A variational approach for calculating flux tube energy is developed to overcome this force error problem in stellarator numerical equilibria. Saturated (equilibrium) flux tube states that cross 10%–20% of the plasma minor radius are shown to exist for linearly ballooning unstable profiles. It is shown that several features of the displaced flux tube structure in a full nonlinear MHD simulation of Wendelstein 7X are reproduced by our model. Saturated states are found in a compact stellarator equilibrium close but below the marginal ballooning linear instability, i.e. the unperturbed equilibrium is metastable. This suggests that edge-localized-mode-like explosive MHD behavior may be possible in stellarators.

Abstract:

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Abstract:

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Abstract:

The toroidal geometry of tokamaks and stellarators is known to play a crucial role in the linear physics of zonal flows (ZFs), leading to e.g. the Rosenbluth–Hinton residual and geodesic acoustic modes. However, descriptions of the nonlinear ZF growth from a turbulent background typically resort to simplified models of the geometry. We present a generalised theory of the secondary instability to model the ZF growth from turbulent fluctuations in toroidal geometry, demonstrating that the radial magnetic drift substantially affects the nonlinear ZF dynamics. In particular, the toroidicity gives rise to a new branch of propagating ZFs, the toroidal secondary mode, which is nonlinearly supported by the turbulence. We present a theory of this mode and compare the theory against gyrokinetic simulations of the secondary mode. The connection with other secondary modes—the ion-temperature-gradient and Rogers–Dorland–Kotschenreuther secondary modes—is also examined.