Engineering topics

Belli, Emily A.

Publications and source records attributed to Belli, Emily A..

5 records

DOE OSTI2023

Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations

The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Optimization and Portability of a Fusion OpenACC-based FORTRAN HPC Code from NVIDIA to AMD GPUs

NVIDIA has been the main provider of GPU hardware in HPC systems for over a decade. Most applications that benefit from GPUs have thus been developed and optimized for the NVIDIA software stack. Recent exascale HPC systems are, however, introducing GPUs from other vendors, e.g. with the AMD GPU-based OLCF Frontier system just becoming available. AMD GPUs cannot be directly accessed using the NVIDIA software stack, and require a porting effort by the application developers. This paper provides an overview of our experience porting and optimizing the CGYRO code, a widely-used fusion simulation tool based on FORTRAN with OpenACC-based GPU acceleration. While the porting from the NVIDIA compilers was relatively straightforward using the CRAY compilers on the AMD systems, the performance optimization required more fine-tuning. In the optimization effort, we uncovered code sections that had performed well on NVIDIA GPUs, but were unexpectedly slow on AMD GPUs. After AMD-targeted code optimizations, performance on AMD GPUs has increased to meet our expectations. Modest speed improvements were also seen on NVIDIA GPUs, which was an unexpected benefit of this exercise.

Sfiligoi, Igor↗

DOE OSTI2022

Neoclassical transport due to resonant magnetic perturbations in DIII-D

In this work, the role of neoclassical physics in the particle and energy transport during the application of resonant magnetic perturbations (RMPs) to suppress the edge localised modes in a tokamak is analysed. The neoclassical fluxes in non-axisymmetric DIII-D equilibria with applied RMPs are calculated using the NEO code. The magnetic field provided to NEO as an input is calculated using M3D-C1 and includes the nonlinear one-fluid plasma response. Neoclassical fluxes obtained in this study are found to dramatically increase in the presence of applied resonant magnetic perturbations, and are in same range as the total radial particle fluxes calculated in comparable RMP discharges in DIII-D. This suggests that neoclassical transport plays a significant role in edge transport when RMPs are present. An increase in neoclassical fluxes during the ELM suppressed phase in DIII-D plasmas is calculated and is strongly correlated with the observation of density pump-out in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Geometry dependence of the fluctuation intensity in gyrokinetic turbulence

The findings of an investigation into the properties of the three dimensional (3D) saturated fluctuation intensity of the electric potential in gyrokinetic turbulence simulations is presented. Scans in flux surface elongation and Shafranov shift are used to isolate the tokamak geometric dependencies. The potential intensity required in order to compute exact fluxes by a quasilinear method is determined using linear eigenmodes computed with the gyrokinetic code. A model of this non-linear intensity is constructed using the linear eigenmode properties and the geometry shape functions obtained from the 3D intensity spectrum. The model computes the poloidal wavenumber spectrum of the electron and ion energy fluxes with unprecedented accuracy. In conclusion, new insights are gained into the way zonal flow mixing saturates ion-scale turbulence by controlling the radial wavenumber width of the turbulence spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Rapidly-convergent flux-surface shape parameterization

We propose a novel flux-surface parameterization suitable for local MHD equilibrium calculations with strongly-shaped flux surfaces. The method is based on a systematic expansion in a small number of intuitive shape parameters, and reduces to the well-known Miller D-shaped parameterization in the limit where some of the coefficients are set to zero. Here, the new parameterization is valid for up-down asymmetric plasmas and provides an improvement to the Miller form. Simultaneously, the method is rapidly convergent and requires only about half the number of shape parameters as a general Fourier representation in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗