Engineering Papers⌕ Search

Engineering topics

Kramer, G. J.

Publications and source records attributed to Kramer, G. J..

Plasma shape effects on the Alfvén eigenmode spectrum through Alfvén slow-magnetosonic wave coupling

In this report the effect of the plasma shaping (triangularity and elongation) on the continuous spectrum, frequency gaps, and number of eigenmodes is studied. The plasma shaping affects the coupling between Alfvén and slow magnetosonic waves via the geodesic magnetic curvature and the plasma pressure. The Alfvén-slow coupling creates a large number of new gaps in the continuous spectrum below the Alfvén frequency where discrete modes can reside. In ideal magnetohydrodynamic simulations a large number (a few hundred or more) of potential eigenmodes are found. The number of eigenmodes is correlated with the maximum geodesic curvature and a minimum number of possible discrete eigenmodes is found at a negative triangularity of -0.3. It is hypothesized that these possible eigenmodes form a low amplitude and dense discrete spectrum, which can be studied experimentally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Integrated package of electron cyclotron emission imaging data processing and forward modeling in OMFIT

An Electron Cyclotron Emission Imaging (ECEI) data analysis module has been developed for the OMFIT platform to accommodate the needs of users at the DIII-D tokamak for physics applications. The user can easily access the ECEI spatial observation windows in the plasma that are calculated based on the automatically retrieved hardware setup and available DIII-D equilibria, perform spectral analysis, and make 2D electron temperature fluctuation images. The module provides a powerful data post processing package for extracting important physics parameters from the 2D measurements, including the radial structure and poloidal mode number of Alfven Eigenmodes, as well as the frequency-vs-wavenumber dispersion relationship of broadband MHD. Lastly, the module propagates characterized synthetic fluctuations for the user, so one can perform forward modeling tasks with simple analytical fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Noise suppression for MHD characterization with electron cyclotron emission imaging 1D technique

Significant noise suppression for magnetohydrodynamics (MHD) mode characterization in the spatial and spectral domain is achieved by processing two-dimensional (2D) electron cyclotron emission imaging (ECEI) data with a one-dimensional (1D) ECEI technique using a short time window ($1\;{\text{ms}}$). The technique is applied to detect toroidal Alfven eigenmodes (TAEs) in the temporal spectrum and fit their radial envelope using the data from the DIII-D tokamak W-band 2D ECEI system. Using the data length (time window) of only 1 ${\text{ms}}$, the 1D ECEI can clearly detect the TAEs (~100 ${\text{kHz}}$) on the spectrum, while similar spectrum quality requires ~10 ${\text{ms}}$ data length with the cross power spectrum between two midplane ECEI channels. The 1D ECEI technique also effectively avoids biased fitting when resolving the fine structure of the TAE's radial envelope. The radially spatial resolution of 1D ECEI is constrained by the finite ECE radiation volume of the ECEI receiver. Here, with forward radiation modeling, we find the DIII-D ECEI system can sensitively measure the even parity MHD activities, for which the mode width is >15 mm, and tearing modes (odd parity MHD activities), for which the island full width is >30 mm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance demonstration of vacuum microwave components critical for the operation of the ITER low-field side reflectometer

Final design studies in preparation for manufacturing have been performed for functional components of the vacuum portion of the ITER Low-Field Side Reflectometer (LFSR). These components consist of an antenna array, electron cyclotron heating (ECH) protection mirrors, phase calibration mirrors, and vacuum windows. Evaluation of these components was conducted at the LFSR test facility and DIII-D. The antenna array consists of six corrugated-waveguide antennas for simultaneous profile, fluctuation, and Doppler measurements. A diffraction grating, incorporated into the plasma-facing miter bend, provides protection of sensitive components from stray ECH at 170 GHz. For in situ phase calibration of the LFSR profile reflectometer, an embossed mirror is incorporated into the adjacent miter bend. Measurements of the radiated beam profile indicate that these components have a small, acceptable effect on mode conversion and beam quality. Baseline transmission characteristics of the dual-disk vacuum window are obtained and are used to guide ongoing developments. Preliminary simulations indicate that a surface-relief structure on the window surfaces can greatly improve transmission. The workability of real-time phase measurements was demonstrated on the DIII-D profile reflectometer. The new automated real-time analysis agrees well with the standard post-processing routine.

47 OTHER INSTRUMENTATION↗

Response of thermal and fast-ion transport to beam ion population, rotation and T e / T i in the DIII-D steady state hybrid scenario

The thermal and fast-ion transport properties of DIII-D steady-state hybrid discharges with normalized beta β N ≳ 3 are studied at low injected torque and an increased electron to ion temperature ratio T e / T i . Linear stability analysis performed with the TGLF turbulent code indicates that a high-k mode is usually dominant at smaller radii, whereas a low-k mode is usually dominant at larger radii in these plasmas. A reduction in the net injected torque from 8.6 to 4.3 N-m leads to reduced E × B shear and hence, an enhanced turbulence that was observed on the Doppler backscattering diagnostic and was also computed with TGLF. As T e / T i in the core was increased from 0.57 to 0.66 by adding electron cyclotron current drive (ECCD) to these plasmas, higher levels of transport are observed with increased high-k modes indicated by TGLF. The fast-ion transport level varied over an order of magnitude in these discharges depending on whether Alfvén eigenmodes, fishbones, or no instabilities were observed. Hybrid plasmas with fishbones have decreased fast-ion transport, compared to plasmas with Alfvén eigenmodes, since they are resonant with a smaller portion of phase space and their resonance is farther from the wall. This reduction in fast-ion transport with ECCD mitigates the increase in turbulent transport, resulting in higher performance than expected during strong electron heating. Similarly, the lowest fast-ion transport was observed in the low torque plasma, which also led to better than expected performance at this torque value. Finally, the thermal and fast-ion transport changes observed as the torque/rotation and T e / T i are varied indicate possible methods for transferring this scenario to a reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Role of fast-ion transport manipulating safety factor profile in KSTAR early diverting discharges

Magnetohydrodynamics (MHD) activity is analyzed for a particular KSTAR discharge category, the main features of which are a high-$q_\mathrm{min}$ and a broad q profile. An interesting finding in this scenario is that a steady-state high-$q_\mathrm{min}$ and broad q profile has been constructed even without a strong off-axis current drive scheme. Time traces of magnetics spectrograms and electron cyclotron emission imaging (ECEI) clearly show that Alfvénic type MHD activity appears as a high-$q_\mathrm{min}$ and broad q profile is formed. Conversely, this activity disappears as the current profiles evolve toward a low-$q_\mathrm{min}$ profile. In this study, the effect of MHD modes and accompanying fast-particle transport on the formation of various current profiles is explored. The frequency patterns of the magnetics spectrograms and ECEI coherence are compared using a NOVA analysis to determine which series of Alfvénic modes are active. Then, the kick-model is applied to potential modes, and used to estimate how much fast-ions are transported and beam driven current profiles are varied. Here, by comparing the reconstructed total current profile from the kinetic EFIT with the beam-driven current profile estimated from the kick model, it is determined how energetic particle transport is responsible for maintaining the safety factor profile and which q profile can provide high β N operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗