Turbulence and flows in the plasma boundary of snowflake magnetic configurations

The power exhaust through the scrape-off layer (SOL) in fusion reactors is expected to be significantly higher than in ITER, thus questioning the extrapolation of the ITER exhaust solution to these devices. The snowflake (SF) magnetic configuration is one of the alternative exhaust configurations being considered to mitigate the heat vessel loads in fusion reactors. The SF configuration features a second-order null of the poloidal magnetic field, i.e. a point where the poloidal magnetic field and its spatial derivatives vanish. As a consequence, the null-point is connected to the wall through four legs, which define four strike points. The presence of the four strike points allows for a heat flux distribution on a larger area compared to standard divertor configurations that feature two strike points. SF configurations are obtained experimentally by generating two first-order X-points close to each other. When the two X-points coincide, a second-order null point is obtained. However, in practice, the two X-points never coincide perfectly and, according to their relative position, we distinguish between the snowflake plus (SF+) and snowflake minus (SF-). The configuration with the two X-points coinciding is usually referred to as ideal SF. All these configurations have been experimentally investigated in the TCV, NSTX, EAST, and DIII-D tokamaks and are considered for the DTT tokamak. Numerical simulations of SF configurations, carried out by means of the EMC3-Eirene and the SOLPS codes, are unable to reproduce the heat flux distribution observed experimentally, calling for detailed investigations of the turbulence and flows in the SF plasma boundary. We present the first global turbulent simulations of the plasma dynamics in SF configurations, including the ideal SF, the SF+ and SF- configurations [1]. These simulations carried out by using the GBS code [2], evolve self-consistently the fluctuating and equilibrium quantities, as they result from the interplay of a heat and particle source in the core, turbulent transport, and parallel losses to the vessel wall. The simulations allow us to disentangle the mechanisms behind the heat flux distribution among the different strike points. As pointed out by our simulations, the heat flux can be reduced by a factor of two in the SF configurations, with respect to single-null configurations. The activation of the unconnected strike points in the ideal SF and in the SF+ configurations, also observed in the experiments, is due to the presence of a steady ExB equilibrium flow in the null region that provides a cross-field transport mechanism towards the private flux region. The origin, the properties and the effects of this steady ExB equilibrium flow are carefully analyzed and described as well as its dependence on the direction of the toroidal magnetic field and on the distance between the two X-points.

Turbulent transport regimes in the tokamak boundary

Maurizio Giacomin

The overall performance of a tokamak strongly depends on phenomena that take place in a thin region between the main plasma and the vessel wall, which is denoted as tokamak boundary. In fact, the formation of transport barriers in this region can significantly improve plasma confinement and, therefore, the tokamak fusion performance. In addition, the tokamak boundary determines the peak heat flux to the wall, an essential quantity for the design and the operation of fusion power plants, as well as the level of impurities in the core, the removal of fusion ash and the dynamics of neutral particles.The dynamics in the plasma boundary is strongly nonlinear and characterized by a wide range of length and time scales as well as by a complex magnetic field geometry that mayfeature one or more nulls of the poloidal magnetic field. Large-scale, three-dimensional turbulence simulations are therefore often required to disentangle the complex physicalmechanisms that govern this region.The thesis is focused on the analysis of the different turbulent transport regimes present in the plasma boundary as they appear from three-dimensional, flux-driven, global, two-fluid turbulence simulations carried out by using the GBS code, which is significantly extended here to allow the self-consistent simulation of the plasma dynamics coupled to a kinetic single-species neutral model in arbitrarily complex magnetic geometries. Considering single-null magnetic configurations, three turbulent transport regimes are identified: (i) a regime of suppressed turbulent transport at low values of collisionality and large values of heat source, (ii) a regime of developed turbulent transport at intermediate values of collisionality and heat source, and (iii) a regime of very large turbulent transport at high value of collisionality and density, which can be associated to the crossing of the density limit. By leveraging the results of GBS simulations, theory-based scaling laws of the pressure and density decay lengths in the near and far scrape-off layer are derived in the developed transport regime from a balance among heat source, turbulent transport across the separatrix and parallel losses at the vessel wall. The theoretical scaling of the pressure decay length in the near scrape-off layer is successfully validated against a multi-machine database of SOL width measurements at the outer target.By carefully analysing the transition to the regime of large turbulent transport, we show that the density limit can be explained by an enhancement of turbulent transport at thetokamak boundary when the density increases. This analysis leads to a theory-based scaling law of the maximum edge density achievable in tokamaks, which is in better agreement with a multi-machine database than the widely used Greenwald empirical scaling, thus significantly improving our understanding and predictive capability of thedensity limit, with important implications for the design and the operation of future fusion power plants.The thesis concludes by presenting the first turbulent simulations carried out in various snowflake magnetic configurations, which are used to investigate the effect of turbulenceand equilibrium flow on the heat flux distribution among the different strike points.

EPFL2022

A new generation plasma turbulence fluid code for tokamak scrape-off layer dynamics

Federico David Halpern, Sébastien Jolliet, Joaquim Loizu Cisquella, Annamaria Mosetto, Félix Benedito Clément Musil, Paolo Ricci, Fabio Riva, Minh Quang Tran, Christoph Wersal

Understanding the turbulent dynamics of the tokamak scrape-off layer (SOL) is critical for the safe operation of future fusion reactors. In this region of the tokamak device, the plasma exhausts particles and heat onto plasma facing hardware components, which may limit their lifetime if material constraints are not met. Recently, we have gained deep insights into SOL dynamics by means of electromagnetic fluid turbulence simulations carried out with the Global Braginskii Solver (GBS). GBS is currently capable of carrying out simulations of the SOL of medium-size tokamaks such as TCV or Alcator C-Mod. We present the new version of the GBS code, which is intended to address larger SOL plasmas using advanced two-fluid models. GBS has undergone major transformations, including (a) porting of the code kernel to F2008, (b) upgrading pure MPI parallelism to employ a 3D cartesian communicator, and (c) including a matrix-free parallel multigrid solver to address generalized Poisson and Ampere operators. There is also an ongoing effort to port this code to MIC/GPU architecture. The new GBS code achieves excellent parallel scalability up to 8192 cores for a C-Mod size plasma (e.g. weak/strong scalings), while at the same time it is possible to achieve a fast and scalable solution of time dependent operators. As an example, we discuss two new physics capabilities of GBS allowed by numerical advances. The first is the inclusion of electromagnetic fluctuations at realistic plasma sizes, which is now possible due to the reduced numerical cost and improved parallel scalability. Second, the Boussinesq approximation in the evaluation of the divergence of the polarization drift is not justified in the tokamak scrape-off layer, where dn/n ~ 1. The drift-reduced equations have been rederived, and we have implemented and benchmarked a non-Boussinesq Poisson operator in GBS. Finally, we present the latest GBS simulation results using realistic size and parameters from inner-wall limited discharges, which is the chosen scenario for ITER start-up. Our simulations reveal the presence of steep plasma gradients near the last closed flux surface, which implies that a large portion of the exhaust heat can be deposited in a very small area. The turbulent dynamics recovered show low frequency, large amplitude drift and ballooning modes with a 1cm radial and poloidal correlation lengths. Intermittent transport events are present, leading to strongly skewed fluctuation probability distribution functions. The simulations have been benchmarked against measurements from the same experiments, showing very good agreement in many observables.

2015

Simulation of plasma turbulence in the periphery of diverted tokamaks

Paola Paruta

The turbulent plasma dynamics in the periphery of a fusion device plays a key role in determining its overall performance. In fact, the periphery controls the heat load on the vessel walls, the plasma confinement, the level of impurities in the core, the plasma fuelling and the removal of fusion ashes. Hence, understanding and predicting the plasma turbulence in this region is of crucial importance for the success of the fusion program. The GBS code has been developed in past years to simulate plasma turbulence in the periphery of limited tokamaks. The goal of the present thesis is to extend GBS to the treatment of diverted scenarios. Such configurations are of interest for present state-of-the-art experiments and future fusion reactors. For the implementation of this geometry, we express the model in toroidal coordinates, abandoning the flux coordinates previously used in limited configuration, and overcoming the singularity that this coordinate system presents at the X-point of diverted configurations. The accuracy of the numerical scheme is improved by upgrading the second order finite differences scheme to fourth order on staggered grids. The resulting version of GBS is carefully verified through a series of tests (i.e., a benchmark with the previous version of GBS in limited configuration, a rigorous check of the correctness of the code implementation with the method of manufactured solutions, and a convergence study on a relatively simple diverted configuration). The results of a GBS simulation is then used to investigate the dynamics of coherent turbulent structures, called blobs, that characterise plasma turbulence in the periphery of fusion devices. A diverted double-null configuration is considered, and the blob motion is studied using a pattern recognition algorithm. The velocity of the blobs in the presence of an X-point matches the analytical scaling that we derived by considering the different blob properties in the divertor and main SOL regions, retaining the correction terms that account for blob density and ellipticity. In addition, we show that the blob current pattern observed in the simulation results match the theoretical expectations. Finally, the new version of GBS is run with a realistic diverted magnetic equilibrium, taken from an experiment carried out on the TCV tokamak. First insights of the turbulence properties are in good agreement with the current physical understanding of plasma dynamics in the periphery of diverted tokamaks.

EPFL2018