Gas puff imaging (GPI) is a diagnostic of plasma turbulence which uses
a puff of neutral gas at the plasma edge to increase the local visible
light emission for improved space-time resolution of plasma
fluctuations. This paper reviews gas puff imaging diagnostics of edge
plasma turbulence in magnetic fusion research, with a focus on the
instrumentation, diagnostic cross-checks, and interpretation
issues. The gas puff imaging hardware, optics, and detectors are
described for about 10 GPI systems implemented over the past ~15
years. Comparison of GPI results with other edge turbulence diagnostic
results are described and many common features are observed. Several
issues in the interpretation of GPI measurements are discussed, and
potential improvements in hardware and modeling are suggested.
Petsev, Nikolai D.; Nikoubashman, Arash; Latinwo, Folarin
Abstract:
Source code for our genetic algorithm optimization investigation of conglomerate and racemic chiral crystals. In this work, we address challenges in determining the stable structures formed by chiral molecules by applying the framework of genetic algorithms to predict the ground state crystal lattices formed by a chiral tetramer model. Using this code, we explore the relative stability and structures of the model’s conglomerate and racemic crystals, and extract a structural phase diagram for the stable Bravais crystal types in the zero-temperature limit.
Sharma, A. Y.; Cole, M. D. J.; Görler, T.; Chen, Y.; Hatch, D. R.; Guttenfelder, W.; Hager, R.; Sturdevant, B. J.; Ku, S.; Mishchenko, A.; Chang, C. S.
Boronization is commonly utilized in tokamaks to suppress intrinsic impurities, most notably oxygen from residual water vapor. However, this is a temporary solution, as oxygen levels typically return to pre-boronization levels following repeated plasma exposure. The global impurity migration model WallDYN has been applied to the post-boronization surface impurity evolution in NSTX-U. A “Thin Film Model” has been incorporated into WallDYN to handle spatially inhomogeneous conditioning films of varying thicknesses, together with an empirical boron conditioning model for the NSTX-U glow discharge boronization process. The model qualitatively reproduces the spatial distribution of boron in the NSTX-U vessel, the spatially-resolved divertor emission pattern, and the increase in oxygen levels following boronization. The simulations suggest that oxygen is primarily sourced from wall locations without heavy plasma flux or significant boron deposition, namely the lower and upper passive plates and the lower private flux zone.
The materials include codes and example input / output files for Monte Carlo simulations of lattice chains in the grand canonical ensemble, for determining phase behavior, critical points, and formation of aggregates.
Dielectric tensor for crystalline graphite from X-ray to microwave frequencies, as discussed in the paper "Graphite Revisited" (Draine 2016, Astrophysical Journal, in press). Cross sections for absorption and scattering by graphite spheres and spheroids are also tabulated, as well as Planck-averaged cross sections for absorption and scattering of radiation with a Planck spectrum.
Stoltzfus-Dueck, T; Hornsby, W A; Grosshauser, S R
Abstract:
Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.
Hager, R.; Chang, C. S.; Ferraro, N. M.; Nazikian R.
Abstract:
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement.The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response.With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal.In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.
This dataset is affiliated with the publication https://doi.org/10.1007/s00348-022-03455-0. All of the data provided is necessary to reproduce the results with the aforementioned publication. The data in this repository is for the wake of a wind turbine at high Reynolds numbers. The data is mainly used for reproducing the statistics (deficit and variance profiles) and the phase averaged results.
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.
Kaita, R.; Lucia, M.; Allain, J. P.; Bedoya, F.; Capece, A.; Jaworski, M.; Koel, B. E.; Majeski, R.; Roszell, J.; Schmitt, J.; Scotti, F.; Skinner, C. H.; Soukhanovskii, V.
Abstract:
The application of lithium to plasma-facing components (PFCs) has long been used as a technique for wall conditioning in magnetic confinement devices to improve plasma performance. Determining the characteristics of PFCs at the time of exposure to the plasma, however, is difficult because they can only be analyzed after venting the vacuum vessel and removing them at the end of an operational period. The Materials Analysis and Particle Probe (MAPP) addresses this problem by enabling PFC samples to be exposed to plasmas, and then withdrawn into an analysis chamber without breaking vacuum. The MAPP system was used to introduce samples that matched the metallic PFCs of the Lithium Tokamak Experiment (LTX). Lithium that was subsequently evaporated onto the walls also covered the MAPP samples, which were then subject to LTX discharges. In vacuo extraction and analysis of the samples indicated that lithium oxide formed on the PFCs, but improved plasma performance persisted in LTX. The reduced recycling this suggests is consistent with separate surface science experiments that demonstrated deuterium retention in the presence of lithium oxide films. Since oxygen decreases the thermal stability of the deuterium in the film, the release of deuterium was observed below the lithium deuteride dissociation temperature. This may explain what occurred when lithium was applied to the surface of the NSTX Liquid Lithium Divertor (LLD). The LLD had segments with individual heaters, and the deuterium-alpha emission was clearly lower in the cooler regions. The plan for NSTX-U is to replace the graphite tiles with high-Z PFCs, and apply lithium to their surfaces with lithium evaporation. Experiments with lithium coatings on such PFCs suggest that deuterium could still be retained if lithium compounds form, but limiting their surface temperatures may be necessary.
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.
Using a recently installed impurity powder dropper (IPD), boron powder (< 150 μm) was injected into lower single null (LSN) L-mode discharges in WEST. IPDs possibly enable real-time wall conditioning of the plasma-facing components and may help to facilitate H-mode access in the full-tungsten environment of WEST. The discharges in this experiment featured Ip = 0.5 MA, BT = 3.7 T, q95 = 4.3, tpulse = 12–30 s, ne,0 ~ 4×1019 m-2, and PLHCD ~ 4.5 MW. Estimates of the deuterium and impurity particle fluxes, derived from a combination of visible spectroscopy measurements and their corresponding S/XB coefficients, showed decreases of ~ 50% in O+, N+, and C+ populations during powder injection and a moderate reduction of these low-Z impurities (~ 50%) and W (~ 10%) in the discharges that followed powder injection. Along with the improved wall conditions, WEST discharges with B powder injection observed improved confinement, as the stored energy WMHD, neutron rate, and electron temperature Te increased significantly (10–25% for WMHD and 60–200% for the neutron rate) at constant input power. These increases in confinement scale up with the powder drop rate and are likely due to the suppression of ion temperature gradient (ITG) turbulence from changes in Zeff and/or modifications to the electron density profile.
