No Arabic abstract
To faithfully simulate ITER and other modern fusion devices, one must resolve electron and ion fluctuation scales in a five-dimensional phase space and time. Simultaneously, one must account for the interaction of this turbulence with the slow evolution of the large-scale plasma profiles. Because of the enormous range of scales involved and the high dimensionality of the problem, resolved first-principles global simulations are very challenging using conventional (brute force) techniques. In this thesis, the problem of resolving turbulence is addressed by developing velocity space resolution diagnostics and an adaptive collisionality that allow for the confident simulation of velocity space dynamics using the approximate minimal necessary dissipation. With regard to the wide range of scales, a new approach has been developed in which turbulence calculations from multiple gyrokinetic flux tube simulations are coupled together using transport equations to obtain self-consistent, steady-state background profiles and corresponding turbulent fluxes and heating. This approach is embodied in a new code, Trinity, which is capable of evolving equilibrium profiles for multiple species, including electromagnetic effects and realistic magnetic geometry, at a fraction of the cost of conventional global simulations. Furthermore, an advanced model physical collision operator for gyrokinetics has been derived and implemented, allowing for the study of collisional turbulent heating, which has not been extensively studied. To demonstrate the utility of the coupled flux tube approach, preliminary results from Trinity simulations of the core of an ITER plasma are presented.
Scaling laws for the transport and heating of trace heavy ions in low-frequency, magnetized plasma turbulence are derived and compared with direct numerical simulations. The predicted dependences of turbulent fluxes and heating on ion charge and mass number are found to agree with numerical results for both stationary and differentially rotating plasmas. Heavy ion momentum transport is found to increase with mass, and heavy ions are found to be preferentially heated, implying a mass-dependent ion temperature for very weakly collisional plasmas and for partially-ionized heavy ions in strongly rotating plasmas.
Magnetic reconnection in strongly magnetized (low-beta), weakly collisional plasmas is investigated using a novel fluid-kinetic model [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] which retains non-isothermal electron kinetics. It is shown that electron heating via Landau damping (linear phase mixing) is the dominant dissipation mechanism. In time, electron heating occurs after the peak of the reconnection rate; in space, it is concentrated along the separatrices of the magnetic island. For sufficiently large systems, the peak reconnection rate is $cE_{max}approx 0.2v_AB_{y,0}$, where $v_A$ is the Alfven speed based on the reconnecting field $B_{y,0}$. The island saturation width is the same as in MHD models except for small systems, when it becomes comparable to the kinetic scales.
We analyze how the turbulent transport of $mathbf{E}times mathbf{B}$ type in magnetically confined plasmas is affected by intermittent features of turbulence. The latter are captured by the non-Gaussian distribution $P(phi)$ of the turbulent electric potential $phi$. Our analysis is performed at an analytical level and confirmed numerically using two statistical approaches. We have found that the diffusion is inhibited linearly by intermittency, mainly via the kurtosis of the distribution $P(phi)$. The associated susceptibility for this linear process is shown to be dependent on the poloidal velocity $V_p$ and on the correlation time $tau_c$ with a maxima at the time-of-flight $tau_{fl}$. Intermittency does not affect the Kubo number scaling in the strong regime.
A limit for the edge density, ruled by radiation losses from light impurities, is established by a minimal cylindrical magneto-thermal equilibrium model. For ohmic tokamak and reversed field pinch the limit scales linearly with the plasma current, as the empirical Greenwald limit. The auxiliary heating adds a further dependence, scaling with the 0.4 power, in agreement with L-mode tokamak experiments. For a purely externally heated configuration the limit takes on a Sudo-like form, depending mainly on the input power, and is compatible with recent Stellarator scalings.
Nonlinear axisymmetric cylindrical plasma oscillations in magnetized collisionless plasmas are a model for the electron fluid collapse on the axis behind an ultrashort relativisically intense laser pulse exciting a plasma wake wave. We present an analytical description of the strongly nonlinear oscillations showing that the magnetic field prevents closing of the cavity formed behind the laser pulse. This effect is demonstrated with 3D PIC simulations of the laser-plasma interaction. An analysis of the betatron oscillations of fast electrons in the presence of the magnetic field reveals a characteristic Four-Ray Star pattern which has been observed in the image of the electron bunch in experiments [T. Hosokai, et al., Phys. Rev. Lett. 97, 075004 (2006)].