Plasma Physics (physics.plasm-ph)

Abstract: In this work, we present the results of first-of-their-kind nonlinear local gyrokinetic simulations of electromagnetic turbulence at mid-radius in the burning plasma phase of the conceptual high-$\beta$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production). A prior linear analysis in D.Kennedy et al. submitted to Nucl. Fusion [1] reveals the presence of unstable hybrid kinetic ballooning modes and subdominant microtearing modes at binormal scales approaching the ion- Larmor radius. Local nonlinear gyrokinetic simulations, using three different codes, are in qualitative and quantitative agreement and suggest that hybrid kinetic ballooning modes drive very large turbulent transport in the absence of equilibrium flow shear. The heat flux rises to values that exceed the available heating power by orders of magnitude and the turbulent eddies are highly extended radially so that they may not be well described by the local gyrokinetic model. The saturated transport fluxes are extremely sensitive to equilibrium flow shear, and diamagnetic levels of flow shear can suppress the fluxes to more reasonable values on the chosen surface. Given this sensitivity there is a large uncertainty in the saturated fluxes. The possible transport impact of the subdominant microtearing modes is also analysed in isolation by artificially and unphysically removing compressional magnetic perturbations from nonlinear calculations, to suppress the dominant hybrid kinetic ballooning mode. The microtearing heat flux is found to saturate at negligible values, though we cannot exclude the possibility that microtearing turbulence may be more transport relevant in other regions of parameter space.

Abstract: We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$\beta$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production). We examine a range of flux surfaces between the deep core and the pedestal top for the two candidate flat-top operating points of the prototype device (EC and EBW operating points). Local linear gyrokinetic analysis is performed to determine the type of microinstabilities that arise under these reactor-relevant conditions. We find that the equilibria are dominated by a hybrid version of the Kinetic Ballooning Mode (KBM) instability at ion binormal and radial scales, with collisional Microtearing Modes (MTMs) sub-dominantly unstable at very similar binormal scales but different radial scales. We study the sensitivity of these instabilities to physics parameters, and discuss potential mechanisms for mitigating them. The results of this investigation are compared to a small set of similar conceptual reactor designs in the literature. A detailed benchmark of the linear results is performed using three gyrokinetic codes; alongside extensive resolution testing and sensitivity to numerical parameters providing confidence in the results of our calculations, and paving the way for detailed nonlinear studies in a companion article.

Abstract: This paper contains the foundation for a new Particle-In-Cell model for gas discharges, based on Ito diffusion and Kinetic Monte Carlo (KMC). In the new model the electrons are described with a microscopic drift-diffusion model rather than a macroscopic one. We discuss the connection of the Ito-KMC model to the equations of fluctuating hydrodynamics and the advection-diffusion-reaction equation which is conventionally used for simulating streamer discharges. The new model is coupled to a particle description of photoionization, providing a non-kinetic all-particle method with several attractive properties, such as: 1) Taking the same input as a fluid model, e.g. mobility coefficients, diffusion coefficients, and reaction rates. 2) Guaranteed non-negative densities. 3) Intrinsic support for reactive and diffusive fluctuations. 4) Exceptional stability properties. The model is implemented as a particle-mesh model on cut-cell grids with Cartesian adaptive mesh refinement. Positive streamer discharges in atmospheric air are considered as the primary application example, and we demonstrate that we can self-consistently simulate large discharge trees.