Hill, K. W.; Gao, L.; Kraus, B. F.; Bitter, M.; Efthimion, P. C.; Pablant, N. A.; Schneider, M. B.; Thorn, D. B.; Chen, H.; Kauffman, R. L.; Liedahl, D. A.; MacDonald, M. J.; MacPhee, A. G.; Scott, H. A.; Stoupin, S.; Doron, R.; Stambulchik, E.; Maron, Y.; Lahmann, B.
Abstract:
Numerical data used to draw the figures in the manuscript
Schwartz, Jacob A.; Nelson, A. O.; Kolemen, Egemen
Abstract:
Shaping a tokamak plasma to have a negative triangularity may allow operation in an ELM-free L-mode regime and with a larger strike-point radius, ameliorating divertor power-handling requirements. However, the shaping has a potential drawback in the form of a lower no-wall ideal beta limit, found using the MHD codes CHEASE and DCON. Using the new fusion systems code FAROES, we construct a steady-state DEMO2 reactor model. This model is essentially zero-dimensional and neglects variations in physical mechanisms like turbulence, confinement, and radiative power limits, which could have a substantial impact on the conclusions deduced herein. Keeping its shape otherwise constant, we alter the triangularity and compute the effects on the levelized cost of energy (LCOE). If the tokamak is limited to a fixed B field, then unless other means to increase performance (such as reduced turbulence, improved current drive efficiency or higher density operation) can be leveraged, a negative-triangularity reactor is strongly disfavored in the model due to lower \beta_N limits at negative triangularity, which leads to tripling of the LCOE. However, if the reactor is constrained by divertor heat fluxes and not by magnet engineering, then a negative-triangularity reactor with higher B0 could be favorable: we find a class of solutions at negative triangularity with lower peak heat flux and lower LCOE than those of the equivalent positive triangularity reactors.
Gilson, Erik; Lee, H; Bortolon, A; Choe, W; Diallo, A; Hong, SH; Lee, HM; Maingi, R; Mansfield, DK; Nagy, A; Park, SH; Song, IW; Song, JI; Yun, SW; Nazikian, R
Abstract:
Results from KSTAR powder injection experiments, in which tens of milligrams of boron nitride (BN) were dropped into low-power H-mode plasmas, show an improvement in wall conditions in subsequent discharges and, in some cases, a reduction or elimination of edge-localized modes (ELMs). Injected powder is distributed by the plasma flow and is deposited on the wall and, over the course of several discharges, was observed to gradually reduce recycling by 33%, and decrease both the ELM amplitude and frequency. This is the first demonstration of the use of BN for ELM mitigation. In all of these experiments, an Impurity Powder Dropper (IPD) was used to introduce precise, controllable amounts of the materials into ELMy H-mode KSTAR discharges. The plasma duration was between 10 s and 15 s, 𝐼𝑝 = 500 kA, 𝐵𝑇 = 1.8 T, 𝑃NBI = 1.6 MW, and 𝑃ECH = 0.6 MW. Plasma densities were between 2 and 3 × 1019 m−3. In all cases, the pre-fill and startup gas-fueling was kept constant, suggesting that the decrease in baseline D𝛼 emission is in fact due to a reduction in recycling. The results presented herein highlight the viability of powder injection for intra-shot and between-shot wall conditioning.
Lampert,Mate; Diallo,Ahmed; Myra,James R.; Zweben, Stewart J.
Abstract:
Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging (GPI) at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g. ELM filaments) in the scrape-off layer and the at the plasma edge with 2.5 us temporal and 10 mm optical resolution.A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics, and then extended to a database of 169 ELM events.Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average an ELM crash consists of a filament with a circular cross-section which is propelled outwards with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100us until the radial propulsion settles back to the pre-ELM level. The experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current-filament interaction model.
Data supporting the manuscript "Enhancement of edge turbulence concomitant with ELM suppression during boron powder injection in EAST" published in Plasma of Physics, 2021.
A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron frequency compressional (CAE) and global (GAE) Alfven Eigenmodes in low aspect ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.
Mollen Albert; Adams Mark F.; Knepley Matthew G.; Hager Robert; Chang C. S.
Abstract:
The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.
Hammond, K. C.; Laggner, F. M.; Diallo, A.; Doskoczynski, S.; Freeman, C.; Funaba, H.; Gates, D.A.; Rozenblat, R.; Tchilinguirian, G.; Xing, Z.; Yamada, I.; Yasuhara, R.; Zimmer, G.; Kolemen, E.
Abstract:
A scalable system for real-time analysis of electron temperature and density based on signals from the Thomson scattering diagnostic, initially developed for and installed on the NSTX-U experiment, was recently adapted for the Large Helical Device (LHD) and operated for the first time during plasma discharges. During its initial operation run, it routinely recorded and processed signals for four spatial points at the laser repetition rate of 30 Hz, well within the system's rated capability for 60 Hz. We present examples of data collected from this initial run and describe subsequent adaptations to the analysis code to improve the fidelity of the temperature calculations.
A matrix inversion technique is derived to calculate local ion temperature from line-integrated measurements of an extended emission source in an axisymmetric plasma which exactly corrects for both toroidal velocity and radial velocity components. Local emissivity and toroidal velocity can be directly recovered from line-integrated spectroscopic measurements, but an independent measurement of the radial velocity is necessary to complete the temperature inversion. The extension of this technique to handle the radial velocity is relevant for magnetic reconnection and merging compression devices where temperature inversion from spectroscopic measurements is desired. A simulation demonstrates the effects of radial velocity on the determination of ion temperature.
Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet-Parker scaling, where neutrals are unimportant, or a fully coupled Sweet-Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfv\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{ni}$ is the neutral-ion collision frequency and $L$ is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.
Choi, W.; Poli, F. M.; Li, M. H.; Baek, S. G.; Gorenlenkova, M.; Ding, B. J.; Gong, X. Z.; Chan, A.; Duan, Y. M.; Hu, J. H.; Lian, H.; Lin, S. Y.; Liu, H. Q.; Qian, J. P.; Wallace, G.; Wang, Y. M.; Zang, Q.; Zhao, H. L.
Abstract:
Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (n_e,lin ≈ 2.0 × 10^19m^−3), simultaneous injection of a 65/35 mix of 2.45 GHz/4.6 GHz power achieved an LHCD efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker-Planck solver CQL3D. For each discharge, profiles of current and hard x-ray from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.