Dynamic measurement of impurity ion transport in runaway electron plateaus in DIII-D
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
In this work, the saturation physics of ion temperature gradient (ITG) turbulence is probed by studying how amplitudes and scalings with key parameters vary across a hierarchy of reduced models. The models derive from nonlinear fluid equations for toroidal ITG turbulence under approximations to the mode coupling interactions in wavenumber space and the representation of turbulent decorrelation. Mode coupling approximations include local-in-wavenumber treatments like the spectral density of flux in quasilinear theory, a truncation to three nonlinearly interacting waves, and the interactions in a cascade to high radial wavenumber mediated by a single zonal flow. Turbulent decorrelation treatments are based on the triplet correlation time with and without eddy damping. Model fidelity is assessed by the scalings and magnitudes of the squared amplitudes of unstable mode, stable mode, and zonal flow with respect to the flow-damping rate and temperature gradient. It is shown that all models reproduce fundamental scalings, provided they incorporate the coupling of unstable mode, stable mode, and zonal flow. Accurate amplitude prediction requires eddy damping in the triplet correlation time and proper representation of the zonal-flow drive by interactions associated with the radial wavenumber cascade.
Here. a novel dual-reflection configuration is introduced for the International Tokamak Experimental Reactor (ITER) core x-ray spectrometer to fit the allocated space where it will be placed accompanied by moving the detectors backward to reduce the incident radiation dose. The highly oriented pyrolytic graphite, which has a mosaic structure of microscopic crystallites, is chosen for the front reflector motivated by higher x-ray throughput and stronger misalignment tolerance compared to the perfect crystal reflector. In the ITER core x-ray spectrometer, a combination of several reflector-deflected Lines of Sight (LOSs) and a direct LOS is proposed for the first time named X-Ray Crystal Spectroscopy Core (XRCS-Core). The system is optimized to observe lines from externally seeded xenon and the intrinsic tungsten impurity, meeting both port integration needs and measurement requirements. Its spectral performance is simulated using an analytical-raytracing mixed code——XRSA, showing good imaging quality with a spectral resolution higher than 8000. The XRCS-Core system is thought to be applicable in various ITER scenarios through the assessment taking into account the spectrometers’ specifications and the chosen lines’ emissivity in different plasma parameters.
The analysis of a previous paper [A. Cole and R. Fitzpatrick, Phys. Plasmas 13, 032503 (2006)] that maps out all of the two-fluid response regimes of a linear tearing layer interacting with an externally generated resonant magnetic perturbation, in a large aspect-ratio tokamak plasma, is generalized to take into account realistic levels of perpendicular particle transport. A new response-regime map is obtained that differs substantially, in many respects, from the old one. The improved analysis is first employed to find all of the two-fluid growth-rate regimes of a non-interacting low-mode-number tearing mode in a large aspect-ratio tokamak plasma. The analysis is then used to determine the scaling of the error-field penetration threshold with machine parameters in large aspect-ratio tokamak plasmas.
The microwave electric field at breakdown for argon gas in a microgap (generally, gap size <1 mm) over a wide range of pressures and gap sizes was found to exhibit two behavioral pressure regimes. A high microwave electric field was established in the gap of a re-entrant mode resonator. By slowly raising the microwave electric field level, a breakdown level is reached where impact ionization induces low pressure argon gas to breakdown preferentially outside of the gap. At high pressure, breakdown is found to occur preferentially inside the gap, but only after the pressure is high enough for impact ionization to produce breakdown there. Impact ionization has a lower threshold electric field than the boundary-controlled processes, such as secondary electron emission, that would eventually produce breakdown inside the microgap if impact ionization did not occur first. Thus, the lower pressure regime is characterized by breakdown outside of the gap, while in the high-pressure regime, breakdown occurs inside the gap, in both cases by impact ionization. The transition pressure P T between these two pressure regimes depends on the gap size d gap , roughly as P T ∝ d gap –2 . Consistent with the model developed here, this corresponds to the full-amplitude of electron oscillation in the gap field being commensurate with the gap size.
Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.
The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.
A dense plasma focus (DPF) is a compact coaxial plasma gun, which completes its discharge as a Z-pinch, producing short (<100 ns) pulses of ions, x rays, and/or neutrons. Lawrence Livermore National Laboratory recently constructed and began operating a new device, the MJOLNIR (MegaJOuLe Neutron Imaging Radiography) DPF, which is designed for single-pulse flash neutron radiography. This device has achieved neutron yields of up to 4.1×10 11 neutrons/pulse at 3.3 MA peak current, and higher-current commissioning is under way. Like most DPFs, MJOLNIR exhibits variable yields in some configurations. Here, we present evidence of the role of parasitic current paths within the gun in stochastically influencing the yield. First through “conditioning shots,” where new hardware has been introduced, we show that increased run-down and run-in speeds correlate with higher yields. These observations are consistent with current being delivered to the electrodes but not to the main plasma sheath, degrading the implosion-driving force. Once nominal conditions are established, we correlate low-performing discharges with smaller current dip and associated voltage spike for a fixed machine configuration. A snow-plow model is able to recreate small-magnitude current dips through the introduction of a parasitic current path, and particle-in-cell simulations establish how parasitic current paths lower the ion beam energy available to produce neutrons. Finally, we observe an increased likelihood of shots with low yield and smaller current dip with increasing fill pressure.
Energy resolving the atomic hydrogen outflux created within plasmas by charge exchange (CX) of hot plasma ions with cooler hydrogen neutrals is used to infer the ion energy distribution within the plasma, IED p . In high- β plasma with field nulls and ion gyro-radii comparable to the plasma size, the measured ion energy distribution (IED m ) of the CX outflux will depend on the viewing angle and position of the detector. Here, we describe the physics for this, results from a synthetic diagnostic code that contrasts the IED m to the IED p within relatively small and hot field-reversed-configuration plasmas, and how these data can show the presence of magnetic nulls and different orbit classes in the plasma.
Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.
The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.
Between 15% and 30% of edge-localized-mode (ELM) heat flux can be deposited to regions outside the main divertor including the far scrape-off-layer (SOL), private flux region (PFR), and secondary divertor inner target. Analysis shown here demonstrates that type-I and small ELM plasmas are transported to the PFR and to the secondary inner divertor, which is magnetically isolated from the outer divertor, leading to well-defined heat flux profiles and with peak values that can surpass those at the secondary outer target. Such features are consistently observed for pedestal collisionalities from $ν$$^{*}_{e}$ ~1.5 to $ν$$^{*}_{e}$ ~3.9. Heat flux profiles of the examined ELM types feature rippled structures due to bursts in the outer far-SOL region but not on the secondary inner target, causing long decay lengths in the time-averaged ELM profiles. The contribution of each ELM type to the total time-averaged power deposited to the secondary divertor has been evaluated, showing that grassy ELMs contribute ~8%, small ELMs ~67%, and type-I ELMs ~85%. These findings imply that small ELMs may yet pose a concern for future machines if some regions of the main wall are not designed to withstand significant heat and particle fluxes. Due to the low intra-ELM heat flux contribution, however, the grassy ELM regime is an attractive option for an ELMing scenario in future machines.
The thermal particles contributed neoclassical toroidal viscosity (NTV) have been successfully developed and explored by many impressive works such as the study by Shaing et al. [Phys. Plasmas 10, 1443 (2003)] and Zhu et al. [Phys. Rev. Lett. 96, 225002 (2006)]. In this work, the scope of the NTV study is extended to explore the contribution of energetic particles (EPs) through both theory and experiments. In theory, the existence of the NTV torque due to the precessional drift resonance of trapped EPs is identified based on the equivalence between the NTV torque and the perturbed drift kinetic energy [J. Park, Phys. Plasmas 18, 110702 (2011)]. Toroidal modeling with the Magneto Resistive Spectrum - drift Kinetic code [Y. Liu, Phys. Plasmas 15, 112503 (2008)], based on this equivalence, indicates that trapped EPs can contribute a significant amount of the NTV torque. Meanwhile, this work also focuses on developing the dedicated DIII-D experiments in the presence of the n = 2 external magnetic perturbation to verify the EP induced NTV (EP-NTV) by measuring the change of the NTV torque while varying the angle and the voltage of the neutral beam injection. However, the developed experiments have been unable to create conditions necessary to clearly demonstrate the presence of EP-NTV. The main challenge is separating the resonant and non-resonant momentum transport responses in the plasma. Finally, the experience, gained from this study, can help the further exploration of EP-NTV in the future experiments.
The time dependence of the blob pulse shape and the waiting time between blobs was found using data from the gas puff imaging diagnostic in NSTX. The database used was of 103 shots from 2010 as described in a previous paper (Zweben et al., Phys. Plasmas 29, 012505 (2022).]. The blob pulse shape was well fit by an exponential rise and fall where the average rise time was τ r = 9.0 ± 2.7 μs and the average fall time was τ f = 16.6 ± 5.8 μs. The waiting times between blob pulses above a threshold of three times the mean had a broad distribution with an average of τ w = 1.2 ± 0.85 ms over the database. The blob intermittency parameter γ b = τ d /τ w , where the blob pulse duration was τ d = τ r + τ f , ranged from γ b ~1% to 5% for shots in this database and increased almost linearly with the blob fraction.
A new, four-chord, CO 2 /He–Ne heterodyne interferometer has been designed and built for measuring line-averaged plasma density in the HIT-SI3 and subsequent HIT-SIU sustained spheromak devices. The two-color system successfully eliminates vibration-induced errors caused by mirrors that are secured to the vacuum chamber and is able to resolve electron densities ne in the full operating range of 10 18 –10 20 m –3 in both experiments with an integrated error of 4.7 × 10 17 m –2 . Here, data are presented from high toroidal current plasma discharges, showing the time evolution of electron densities n e and j Φ /n e along multiple chords.
Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.
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.
Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.