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.
Explanation of ST40 Langmuir probe and infrared thermography diagnostics and their planned upgrades as part of the Lithium Evaporations to Advance PFCs in ST40 (LEAPS) program.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.
The Tore Supra tokamak was transformed into an X-point divertor fusion device in the frame of the WEST (W-for tungsten-Environment in Steady-state Tokamak) project, launched in support of the ITER tungsten divertor strategy. WEST began operation in 2017. A new reciprocating probe head, equipped with a Titanium-Zirconium-Molybdenum (TZM) heat shield, was recently built and operated on WEST during the C9 Campaign in early 2024. TZM was chosen as armour material for its very good mechanical properties between 700 °C and 1400 °C. This material was already extensively used for probe heads and other plasma facing components on the Alcator C-Mod tokamak. Here, this paper presents the design of the new WEST probe head. It is equipped with tungsten filaments that can be heated to emit electrons and provide a direct measurement of the plasma potential in the scrape-off layer, and it is also equipped with a pair of Langmuir probes. The probe plunges vertically into the plasma, reaches its maximum dive depth, then returns to its protected resting position in about 0.25 s. Thermal loads are calculated using as input the thermal properties of TZM and radial profiles of the heat flux measured in the past by other reciprocating probes, and scaled upwards to the expected heat loads in the WEST scrape-off layer. An analytical estimation of Electromagnetic (EM) loads will also be presented, as well as a thermal and structural analysis of the probe for 15 cycles at maximum plunge depth. The probe is observed by an infrared camera. Comparisons between the heat flux deduced by the Langmuir probe and the infrared radiance of the probe housing during plasma exposure are finally presented.
Measurements of Balmer-alpha (B-α) and Balmer-gamma (B-γ), where B = H, D, or T isotopes of hydrogen, emission across the low-field side (LFS) divertor in JET ITER-like Wall (ILW), low-confinement mode (L-mode) plasma indicate that the divertor plasma was more strongly detached in deuterium, tritium and deuterium–tritium than in hydrogen plasmas. These measurements are consistent with the lower degree of detachment in hydrogen versus tritium plasmas inferred from Langmuir probes (Groth et al., 2023). The ratio of B-γ to B-α is a suitable indicator of high density/low temperature conditions because B-α emission is dominated by electron impact excitation with a small electron–ion recombination component regardless of electron temperature and density (for conditions in tokamak plasmas) whilst for B-γ the recombination component increases and the excitation component decreases as the electron temperature decreases and the electron density increases. Within the uncertainty of the measurements, the electron density at the low-field side separatrix and the power crossing the separatrix are similar to permit these comparisons. The B-©/B-α ratio is well correlated with the onset of detachment as measured by Langmuir probes (McCracken et al., 1998; Groth et al., 2023). The ratios of B-γ and B-α varied between 0.005 and 0.035 for low-recycling, high-recycling and partially detached conditions, and were 60–70 % higher in partially detached hydrogen than in tritium plasmas. In contrast, the difference between deuterium, tritium and 40–60 % deuterium–tritium is within 10 % and thus within the uncertainties of the measurements. Profiles of the B-α and B-γ emission was measured with two filtered CCDs connected to the visible lost light leg of JET’s mirror link spectroscopy system (Meigs et al., 2010; Lomanowski et al., 2014). Spatially (radial) narrower B-α profiles were measured with increasing isotope mass, while similarly broadened B-γ profiles were measured.
Abstract We describe the results of a multi-scale, multi-physics modeling assessment of SOLPS-ITER, hPIC2, RustBCA and Xolotl, in which five single-crystal tungsten (W) samples were placed in a reciprocating collector probe and exposed to helium (He) plasma in the WEST fusion device. In our models, we considered a pure (100 %) He plasma, as well as one with oxygen (O) present (95% He 5% O) corresponding to the impurity concentration estimated during the C4 He campaign in WEST. Our SOLPS simulations approximately match experimental reciprocating Langmuir probe plasma measurements of plasma density and temperature. Using these plasma parameters as input, hPIC2 and RustBCA predict that the presence of oxygen impurities lead to a 15%–20% decrease in ion and heat fluxes to the surface, and an order of magnitude higher sputtering yields (compared with a pure He plasma). Xolotl predictions for the response of tungsten to plasma surface interactions (PSIs) agree with experimental LAMS analysis, and indicate large near-surface He concentrations, which quickly decay with depth. Our model also shows an increasing role of erosion—in removing the near-surface He—with time. Overall, slightly higher retention is predicted for tungsten exposed to a pure He plasma, with the largest differences in the near-surface gas content caused by the large oxygen-induced erosion. This highlights the important role that impurities play in PSI. Therefore, future work will focus on providing a fully self-consistent description of oxygen (and oxides, etc.) in our models, through multi-species implementation in GITR and inclusion of oxygen and tungsten oxide formation in Xolotl.
Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.
Gold-tantalum (Au—Ta) alloys are promising for next-generation hohlraums for magnetically-assisted indirect-drive inertial confinement fusion. Hohlraum fabrication involves sputter deposition of ultrathick coatings on sphero-cylindrical substrates. The control of physical properties of such coatings remains a challenge. Here, in this work, we systematically study effects of the growth rate and dynamic substrate tilt on properties of AuTa 4 films deposited onto rotating planar substrates mounted at different tilt angles to mimic different regions of the sphero-cylindrical surface of a hohlraum. We use high-power impulse magnetron sputtering (HiPIMS) in the constant charge per pulse mode with a pulse duration of 100 μs, a charge per pulse of 190 μC, a peak target current of about 3 A, and the pulse frequency varied in the range of 400–1400 Hz. The Langmuir probe, mass-resolved ion energy spectrometry, and optical emission spectroscopy are used to monitor plasma discharge characteristics in order to isolate and study effects of the growth rate. The deposition rate and the ballistics and energetics of depositing species are estimated by Monte Carlo simulations. Results show that the film microstructure, crystallographic phase, residual stress, and electrical resistivity strongly depend on both the deposition rate and substrate tilt, highlighting their critical role in tailoring properties of Au—Ta films during hohlraum fabrication.
Validated and accurate edge profiles (temperature, density, etc.) are vitally important to the Mega Ampere Spherical Tokamak Upgrade (MAST-U) divertor and confinement effort. Density profile reflectometry has the potential to significantly add to the measurement capabilities currently available on MAST-U (e.g., Thomson scattering and Langmuir probes). This work presents the diagnostic requirements, problems, and solutions facing profile reflectometry in spherical tokamaks and MAST-U in particular. Requirements include density measurements near zero electron density in the scrape off layer region, coverage for a broad range of MAST-U plasma parameters, high time (≤10 microseconds) and spatial resolutions (≤1 cm), reliability, and identification of the plasma start frequency.
A particle balance analysis was conducted during a deuterium (D 2 ) shattered pellet injection-induced plasma shutdown on the DIII-D tokamak to determine why less than 20% of the pellet material is assimilated into the core plasma by the mid-current quench (CQ). Initially, most of the D 2 is injected as frozen shards and ionized upon entering the vessel. During the thermal quench, ionized particles move to the divertors and subsequently to the center post (CP) walls, where they rapidly recycle and partially accumulate as neutrals without assimilating into the core plasma. In contrast, the particle flux to the outer midplane walls is negligible, despite being accompanied by hot plasma with electron temperatures exceeding 100 eV. During mid-CQ, volume recombination effects, although not large enough to impact overall particle balance, were significant enough to require accounting for accurate interpretation of fast-framing camera D-alpha signals and the estimation of the CP wall particle flux. In addition, toroidal asymmetries, observed in measurements of toroidal electron density perturbations and the phase of magnetohydrodynamic modes, are present throughout the shutdown and can account for a discrepancy in the assimilation rate for up to 50% of the observed D 2 particle inventory. These sources and sinks of particles and fluxes were identified using absolutely calibrated D-alpha brightness and Langmuir probes.
While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.
Determining nanoparticle charge is more challenging than that for microparticles due to change in the particle size during the synthesis substantial plasma property variations, and difficulties in visualizing individual particles, rendering conventional microparticle charge diagnostics ineffective in dusty plasma. In this work, we utilized laser-stimulated photodetachment (LSPD) to deduce the mean charge of nanoparticles. Nanoparticles were grown in an Ar/C 2 H 2 mixture using a capacitively coupled RF discharge and the LSPD induced changes in the electron current monitored by a cylindrical Langmuir probe. LSPD signals were obtained and analyzed across different dust growth phases. The prolonged decay of electron current pulses was attributed to the presence of residual negative ions, caused by the effective electrostatic trapping of these ions and the potential post—LSPD re-formation of new ones. The particle charge was estimated by combining the laser-stimulated photodetachment signal from the probe with the dust density obtained from laser-light extinction using the measured nanoparticle size distribution. For a nanoparticles size range of approximately 100–250 nm and mean diameter of $d$ p ∼154.37 nm, the effective mean charge was estimated to be $\langle$$Q$$\rangle$ d $≈$ 37 elementary charge units. The measured charge values are lower than those predicted by orbital motion limited theory, which may be attributed to significant electron depletion in the nanodusty plasma. LSPD results in Ar/C 2 H 2 nano-dusty plasma confirm the applicability of this method for estimating individual nanoparticle charges. However, it has also been demonstrated that electron detachment from residual background negative ions can influence the detachment current decay and must be carefully considered.
The low frequency edge oscillation (LFEO) is a low frequency fluctuation in many plasma quantities in the pedestal region of the I-Mode confinement regime. It is observed on Alcator C-Mod between 10–30 kHz and on ASDEX Upgrade between 5–10 kHz. On both tokamaks it has been previously identified as a geodesic acoustic mode (GAM), however the recent discovery of the edge temperature ring oscillation (ETRO) in a similar frequency and spatial location as the LFEO in I-Modes on EAST has called this identification into question. In this paper we investigate the LFEO on C-Mod and AUG using a variety of different experimental techniques including spectral analysis, magnetic mode number analysis, localization, and direct measurement of the LFEO zonal structure and propagation using a mirror Langmuir probe. This investigation has reconfirmed the identification of the LFEO as a GAM and determined that it has several key differences from the ETRO.
A new self-consistent 1D scrape-off layer model has been recently developed in BOUT++ framework, named SD1D, which includes equations for various particle species (e.g. main plasma, neutrals and impurities) and couples open databases like ADAS and AMJUEL. It is able to quickly and effectively simulate divertor detachment experiments. In this work, a typical detachment experiment (shot #6270) on HL-3 with neon seeding is simulated using the SD1D code. It is found that the target electron temperature and the target ion saturation current in the simulations are consistent with experimental results measured by Langmuir probes on the target plate. The variation of D α radiation intensity in the divertor is qualitatively similar to the measured D α signal. Following the experimental validations, different upstream densities are set in the simulations to study the impurity distribution under different plasma density conditions. It is found that increasing upstream density can be helpful for the control of the neon radiation front (closer to the target). In this work we also compare two detachment regimes in simulations. Based on the same initial experimental parameters (shot #6270) on HL-3, a scan of upstream density and a scan of neon seeding rate are carried out respectively. It is found that the role of atomic and molecular processes is different in the two detachment regimes. The current density roll-over is ascribed to a drop in the divertor ion source, and the variation of D α radiation intensity via different excitation channels is associated with the relevant collisional reaction sources.
KSTAR has recently undergone an upgrade to use a new tungsten divertor to run experiments in ITER-relevant scenarios. Even with a high melting point of tungsten, it is important to control the heat flux impinging on tungsten divertor targets to minimize sputtering and contamination of the core plasma. Heat flux on the divertor is often controlled by increasing the degree of detachment of scrape-off layer plasma from the target plates. In this work, we have demonstrated successful divertor detachment and heat exhaust dissipation control experiments using two different methods. The first method uses attachment fraction as a control variable which is estimated using ion saturation current measurements from embedded Langmuir probes in the divertor. The second method uses a novel machine-learning-based surrogate model of 2D UEDGE simulation database, DivControlNN. We demonstrated running inference operation of DivControlNN in realtime to estimate heat flux at the divertor and use it as the control variable in a feedback loop with impurity gas flow. We present interesting insights from these experiments including a systematic approach to tuning controllers and discuss future improvements in the control infrastructure and control variables for future burning plasma experiments.
Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.