Final Technical Report Transport Task Force

The Transport Task Force has functioned as the primary scientific organization in the area of magnetic-fusion confinement and transport since its inception in 1988. It has defined and set research directions, coordinated broad research efforts, advocated new funding initiatives, and created a highly successful and widely admired interactive culture between experiment, theory and modeling. The Transport Task Force carries out its activities under the direction of its chair and the Executive Committee. The Executive Committee is comprised of the leaders and deputy leaders of the scientific working groups. The working groups are structured and organized according to research needs and priorities and have been organized around the areas of Core Transport, H Mode and Pedestal, Fast Particle Transport, Transient Transport Phenomena, and Modeling and Simulation. A steering committee provides advise on TTF activities. Further information on the working groups and the structure and management of the TTF can be found at http://psfcwww2.psfc.mit.edu/ttf/index.html. The TTF holds an annual workshop. A summary of the workshops held during the period of this report is given in Appendix I. During the period of this report the Transport Task Force was involved in several significant activities. Foremost of these was a sweeping review of the status of transport science, the key research tasks for progress during the next 5-10 years, and a proposal for a funding initiative to ensure application of adequate resources to these problems. The conclusions of this study were incorporated into a white paper, which is copied below in Appendix II. Other significant activities have included the introduction of an extended, ongoing discussion on verification and validation as a requisite for defining and codifying the path toward predictive capability, the orchestration of a gradual shift of focus from ion thermal confinement to electron thermal confinement, and a joining of efforts on edge physics by coordinating and uniting efforts of the Transport Task Force and Edge Coordinating Committee. During the next biennium the TTF is chaired by Keith Burrell

Threshold Heat-Flux Reduction by Near-Resonant Energy Transfer

Near-resonant energy transfer to large-scale stable modes is shown to reduce transport above the linear critical gradient, contributing to the onset of transport at higher gradients. This is demonstrated for a threshold fluid theory of ion temperature gradient turbulence based on zonal-flow-catalyzed transfer. The heat flux is suppressed above the critical gradient by resonance in the triplet correlation time, a condition enforced by the wave numbers of the interaction of the unstable mode, zonal flow, and stable mode.</p

Saturation Physics of Threshold Heat-Flux Reduction

The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. The heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.</p

Three-dimensional shear-flow instability saturation via stable modes

Turbulence in three dimensions (3D) supports vortex stretching that has long been known to accomplish energy transfer to small scales. Moreover, net energy transfer from large-scale, forced, unstable flow-gradients to smaller scales is achieved by gradient-flattening instability. Despite such enforcement of energy transfer to small scales, it is shown here that the shear-flow-instability-supplied 3D-fluctuation energy is largely inverse-transferred from the fluctuation to the mean-flow gradient, and such inverse transfer is more efficient for turbulent fluctuations in 3D than in two dimensions (2D). The transfer is due to linearly stable eigenmodes that are excited nonlinearly. The stable modes, thus, reduce both the nonlinear energy cascade to small scales and the viscous dissipation rate. The vortex-tube stretching is also suppressed. Up-gradient momentum transport by the stable modes counters the instability-driven down-gradient transport, which also is more effective in 3D than in 2D ( ≈ 70 % vs ≈ 50 % ). From unstable modes, these stable modes nonlinearly receive energy via zero-frequency fluctuations that vary only in the direction orthogonal to the plane of 2D shear flow. The more widely occurring 3D turbulence is thus inherently different from the commonly studied 2D turbulence, despite both saturating via stable modes.</p

Effect of Triangularity on Ion-Temperature-Gradient-Driven Turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven (ITG) turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. The effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ| ≤ 0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.</p

Global linear and nonlinear gyrokinetic simulations of tearing modes

To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes (n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.</p

Kinetic Turbulence

The weak collisionality typical of turbulence in many diffuse astrophysical plasmas invalidates an MHD description of the turbulent dynamics, motivating the development of a more comprehensive theory of kinetic turbulence. In particular, a kinetic approach is essential for the investigation of the physical mechanisms responsible for the dissipation of astrophysical turbulence and the resulting heating of the plasma. This chapter reviews the limitations of MHD turbulence theory and explains how kinetic considerations may be incorporated to obtain a kinetic theory for astrophysical plasma turbulence. Key questions about the nature of kinetic turbulence that drive current research efforts are identified. A comprehensive model of the kinetic turbulent cascade is presented, with a detailed discussion of each component of the model and a review of supporting and conflicting theoretical, numerical, and observational evidence.Comment: 31 pages, 3 figures, 99 references, Chapter 6 in A. Lazarian et al. (eds.), Magnetic Fields in Diffuse Media, Astrophysics and Space Science Library 407, Springer-Verlag Berlin Heidelberg (2015

On the effect of flux-surface shaping on trapped-electron modes in quasi-helically symmetric stellarators

Using a novel optimization procedure, it has been shown that the Helically Symmetric eXperiment stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [Gerard et al., Nucl. Fusion 63, (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available energy E A , quasi-helical symmetry, and flux-surface shaping parameters as metrics for TEM stabilization. It is found that lower values of E A correlate with reduced growth rates, but only when separate flux-surface shaping regimes are considered. Moreover, configurations with improved quasi-helical symmetry demonstrate a similar reduction in growth rates and less scatter compared to E A . Regarding flux-surface shaping, a set of helical shaping parameters is introduced that show increased elongation is strongly correlated with reduced TEM growth rates, however, only when the quasi-helical symmetry is preserved. Using a newly derived velocity-space-averaged TEM resonance operator, these trends are analyzed to provide insights into the physical mechanism of the observed stabilization. For elongation, stabilization is attributed to geometric effects that reduce the destabilizing particle drifts across the magnetic field. Regarding quasi-helical symmetry, the TEM resonance in the maximally resonant trapping well is shown to increase as the quasi-helical symmetry is broken, and breaking quasi-helical symmetry increases the prevalence of highly resonant trapping wells. While these results demonstrate the limitations of using any single metric as a linear TEM proxy, it is shown that quasi-helical symmetry and plasma elongation are highly effective metrics for reducing TEM growth rates in helical equilibria.</p