No Arabic abstract
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
Turbulence is a major factor limiting the achievement of better tokamak performance as it enhances the transport of particles, momentum and heat which hinders the foremost objective of tokamaks. Hence, understanding and possibly being able to control turbulence in tokamaks is of paramount importance, not to mention our intellectual curiosity of it.
Results of the first validation of large guide field, $B_g / delta B_0 gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $beta_i = 0.01$ and solar wind relevant $beta_i = 1$ are presented, where $delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $beta_i B_g / delta B_0 gtrsim 1$ and $B_g / delta B_0 gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.
Recent efforts to include kinetic effects in fluid simulations of plasmas have been very promising. Concerning collisionless magnetic reconnection, it has been found before that damping of the pressure tensor to isotropy leads to good agreement with kinetic runs in certain scenarios. An accurate representation of kinetic effects in reconnection was achieved in a study by Wang et al. (Phys. Plasmas, volume 22, 2015, 012108) with a closure derived from earlier work by Hammett and Perkins (PRL, volume 64, 1990, 3019). Here, their approach is analyzed on the basis of heat flux data from a Vlasov simulation. As a result, we propose a new local closure in which heat flux is driven by temperature gradients. That way, a more realistic approximation of Landau damping in the collisionless regime is achieved. Previous issues are addressed and the agreement with kinetic simulations in different reconnection setups is improved significantly. To the authors knowledge, the new fluid model is the first to perform well in simulations of the coalescence of large magnetic islands.
The results of flux-driven, two-fluid simulations in single-null configurations are used to investigate the processes determining the turbulent transport in the tokamak edge. Three turbulent transport regimes are identified: (i) a developed transport regime with turbulence driven by an interchange instability, which shares a number of features with the standard L-mode of tokamak operation, (ii) a suppressed transport regime, characterized by a higher value of the energy confinement time, low-amplitude relative fluctuations driven by a Kelvin-Helmholtz instability, a strong E x B sheared flow, and the formation of a transport barrier, which recalls the H-mode, and (iii) a degraded confinement regime, characterized by a catastrophically large interchange-driven turbulent transport, which reminds the crossing of the Greenwald density limit.We derive an analytical expression of the pressure gradient length in the three regimes. The transition from the developed to the suppressed transport regime is obtained by increasing the heat source or decreasing the collisionality and vice versa for the transition from the developed transport regime to the degraded confinement regime. An analytical expression of the power threshold to access the suppressed transport regime, linked to the power threshold for H-mode access, as well as the maximum density achievable before entering the degraded confinement regime, related to the Greenwald density, are also derived. The experimental dependencies of the power threshold for H-mode access on density, tokamak major radius, and isotope mass are retrieved. The analytical estimate of the density limit contains the correct dependence on the plasma current and on the tokamak minor radius.
A minimal model for magnetic reconnection and, generally, low-frequency dynamics in low-beta plasmas is proposed. The model combines analytical and computational simplicity with physical realizability: it is a rigorous limit of gyrokinetics for plasma beta of order the electron-ion mass ratio. The model contains collisions and can be used both in the collisional and collisionless reconnection regimes. It includes gyrokinetic ions (not assumed cold) and allows for the topological rearrangement of the magnetic field lines by either resistivity or electron inertia, whichever predominates. The two-fluid dynamics are coupled to electron kinetics --- electrons are not assumed isothermal and are described by a reduced drift-kinetic equation. The model therefore allows for irreversibility and conversion of magnetic energy into electron heat via parallel phase mixing in velocity space. An analysis of the exchanges between various forms of free energy and its conversion into electron heat is provided. It is shown how all relevant linear waves and regimes of the tearing instability (collisionless, semicollisional and fully resistive) are recovered in various limits of our model. An efficient way to simulate our equations numerically is proposed, via the Hermite representation of the velocity space. It is shown that small scales in velocity space will form, giving rise to a shallow Hermite-space spectrum, whence it is inferred that, for steady-state or sufficiently slow dynamics, the electron heating rate will remain finite in the limit of vanishing collisionality.