Dynamical Evolution of Thin Current Sheets in a Three-Dimensional Open System

Dynamical behavior of thin current sheets under the influence of collisionless reconnection in an open system is investigated by using newly developed electromagnetic (EM) particle simulation codes. In a three-dimensional open system collisionless driven reconnection evolves dynamically under the influence of an external driving flow and three different types of plasma instabilities excited in a thin current sheet. Driving electric field imposed at the boundary penetrates into the current sheet in accordance with the propagation of the lower hybrid drift wave excited in the periphery. When the electric field reaches the neutral sheet, collisionless reconnection is triggered. The current sheet is split as a result of collisionless reconnection, and thus small islands appear in the downstream. The accumulation of current density inside the islands excites the kink instability leading to the destruction of the island structure. A low frequency EM instability is triggered in the current sheet after the island structure disappears in the system

Multi-scale interactions between turbulence and magnetohydrodynamic instability driven by energetic particles

In order to realize high performance burning plasmas in magnetic-confinement fusion devices, such as tokamaks, both bulk plasma transport and that of energetic fusion alpha-particles, which result from different scale fluctuations with different free energy sources, have to be reduced simultaneously. Utilizing the advantage of global toroidal non-linear simulations covering a whole torus, here, we found a new coupling mechanism between the low-frequency micro-scale electromagnetic drift-wave fluctuations regulating the former, while the high-frequency macro-scale toroidal Alfven eigenmode (TAE) regulates the latter. This results from the dual spread of micro-scale turbulence due to the macro-scale TAE not only in wavenumber space representing local eddy size but also in configuration space with global profile variations. Consequently, a new class of turbulent state is found to be established, where the turbulence is homogenized on the poloidal cross-section with exhibiting large-scale structure, which increases fluctuation levels and then both transports, leading to deterioration in the fusion performance

Electromagnetic gyrokinetic simulation of turbulence in torus plasmas

Gyrokinetic simulations of electromagnetic turbulence in magnetically confined torus plasmas including tokamak and heliotron/stellarator are reviewed. Numerical simulation of turbulence in finite beta plasmas is an important task for predicting the performance of fusion reactors and a great challenge in computational science due to multiple spatio-temporal scales related to electromagnetic ion and electron dynamics. The simulation becomes further challenging in non-axisymmetric plasmas. In finite beta plasmas, magnetic perturbation appears and influences some key mechanisms of turbulent transport, which include linear instability and zonal flow production. Linear analysis shows that the ion-temperature gradient (ITG) instability, which is essentially an electrostatic instability, is unstable at low beta and its growth rate is reduced by magnetic field line bending at finite beta. On the other hand, the kinetic ballooning mode (KBM), which is an electromagnetic instability, is destabilized at high beta. In addition, trapped electron modes (TEMs), electron temperature gradient (ETG) modes, and micro-tearing modes (MTMs) can be destabilized. These instabilities are classified into two categories: ballooning parity and tearing parity modes. These parities are mixed by nonlinear interactions, so that, for instance, the ITG mode excites tearing parity modes. In the nonlinear evolution, the zonal flow shear acts to regulate the ITG driven turbulence at low beta. On the other hand, at finite beta, interplay between the turbulence and zonal flows becomes complicated because the production of zonal flow is influenced by the finite beta effects. When the zonal flows are too weak, turbulence continues to grow beyond a physically relevant level of saturation in finite-beta tokamaks. Nonlinear mode coupling to stable modes can play a role in the saturation of finite beta ITG mode and KBM. Since there is a quadratic conserved quantity, evaluating nonlinear transfer of the conserved quantity from unstable modes to stable modes is useful for understanding the saturation mechanism of turbulence

Simulation studies on temperature profile stiffness in ITG turbulent transport of helical plasmas for flux-matching technique

In the framework of the flux-matching method, which is a useful way for the validation of the gyrokinetic turbulence simulations, it is strongly demanded to evaluate the plasma profile sensitivity of the transport coefficients obtained in the employed simulation model within the profile gradient ranges estimated from the experimental observations. The sensitivity causes the plasma profile stiffness for wide ranges of the transport fluxes. In the nonlinear gyrokinetic simulations for the ion temperature gradient (ITG) turbulence in the Large Helical Device (LHD) [Takeiri et al., Nucl. Fusion 57, 102023 (2017)], it is found that the temperature gradients around the experimental nominal observations are slightly larger than the threshold of the instability, and the ion heat diffusivities are quite sensitive to the temperature gradient. The growth rates of the instability, the generations of the zonal flows, and the sensitivities of the transport coefficients to the temperature profiles depend on the radial locations, the employed simulation models, and the field configurations. Specifically, in the optimized LHD field configuration, the sensitivities are relaxed in the outer radial region due to the enhancement of the zonal flows and the reduction of the ITG instability. In order to estimate the range of the temperature gradients possible given the experimentally obtained data of the temperature with errorbars, the statistical technique, Akaike\u27s Information Criterion [H. Akaike, in Proceedings of the 2nd International Symposium on Information Theory, edited by B. N. Petrov and F. Caski (Akadimiai Kiado, Budapest, 1973), pp. 267–281] is applied. Against the range of the temperature gradients, the flux-matching method to predict the temperature gradient in helical plasmas is demonstrated

Modeling of turbulent particle and heat transport in helical plasmas based on gyrokinetic analysis

The particle and heat transport driven by the ion temperature gradient instability in helical plasmas is investigated by the gyrokinetic analysis taking into account the kinetic electron response. High and low ion temperature plasma cases for the discharge in the Large Helical Device (LHD) are studied. Two types of transport models with a lower computational cost to reproduce the nonlinear gyrokinetic simulation results within allowable errors are presented for application in quick transport analyses. The turbulent electron and ion heat diffusivity models are given in terms of the linear growth rate and the characteristic quantity for the linear response of zonal flows, while the model of the effective particle diffusivity is not obtained for the flattened density profile observed in the LHD. The quasilinear flux model is also shown for the heat transport. The quasilinear flux models for the energy fluxes are found to reproduce the nonlinear simulation results at the accuracy similar to that of the heat diffusivity models. In addition, the quasilinear particle flux model, which is applicable to the transport analysis for LHD plasmas, is constructed. These turbulent reduced models enable coupling to the other simulation in the integrated codes for the LHD

Spatially resolved measurement of helium atom emission line spectrum in scrape-off layer of Heliotron J by near-infrared Stokes spectropolarimetry

1視線の観測のみで核融合プラズマ中のヘリウム近赤外輝線の発光分布を推定. 京都大学プレスリリース. 2022-09-26.For plasma spectroscopy, Stokes spectropolarimetry is used as a method to spatially invert the viewing-chord-integrated spectrum on the basis of the correspondence between the given magnetic field profile along the viewing chord and the Zeeman effect appearing on the spectrum. Its application to fusion-related toroidal plasmas is, however, limited owing to the low spatial resolution as a result of the difficulty in distinguishing between the Zeeman and Doppler effects. To resolve this issue, we increased the relative magnitude of the Zeeman effect by observing a near-infrared emission line on the basis of the greater wavelength dependence of the Zeeman effect than of the Doppler effect. By utilizing the increased Zeeman effect, we are able to invert the measured spectrum with a high spatial resolution by Monte Carlo particle transport simulation and by reproducing the measured spectra with the semiempirical adjustment of the recycling condition at the first walls. The inversion result revealed that when the momentum exchange collisions of atoms are negligible, the velocity distribution of core-fueling atoms is mainly determined by the initial distribution at the time of recycling. The inversion result was compared with that obtained using a two-point emission model used in previous studies. The latter approximately reflects the parameters of atoms near the emissivity peak

Turbulent transport of heat and particles in a high ion temperature discharge of the Large Helical Device

Turbulent transport in a high ion temperature discharge of the Large Helical Device (LHD) is investigated by means of electromagnetic gyrokinetic simulations, which include kinetic electrons, magnetic perturbations, and full geometrical effects. Including kinetic electrons enables us to firstly evaluate the particle and the electron heat fluxes caused by turbulence in LHD plasmas. It is found that the electron energy transport reproduces the experimental result, and that the particle flux is negative. The contribution of magnetic perturbation to the transport is small because of very low beta. The turbulence is driven by the ion temperature gradient instability, and the effect of kinetic electrons enhances the growth rate larger than that from the adiabatic electron calculation. The ion energy flux is larger than that observed in the experiment, while the flux is close to the experimental observation when the temperature gradient is reduced 20% in the simulation. This significant sensitivity of the energy flux implies that the profile in the experiment is close to the critical temperature gradient. The critical gradient for turbulent energy flux is similar to that for the linear instability, i.e., the Dimits shift is small. This is because the zonal flow in the LHD is weaker than that in tokamaks