Engineering topics

Wilcox, Robert S.

Publications and source records attributed to Wilcox, Robert S..

8 records

DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Modeling a divertor with mid-leg pumping for high-power H-mode scenarios in DIII-D considering E x B drift flows

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E x B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial E θ x B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Tungsten erosion and divertor leakage from the DIII-D SAS-VW tungsten-coated divertor in experiments with neon gas seeding

Collector probes have been used to examine tungsten divertor leakage in a variety of scenarios with low-Z impurity seeding during operation with the new tungsten-coated SAS-VW divertor in DIII-D. Measurements of tungsten deposition on collector probes inserted into the far Scrape-off-Layer (SOL) are used to deduce how efficiently tungsten leaks out of the closed, V-shaped divertor after it is eroded from the target surfaces. Qualitative differences in the tungsten deposition patterns across the collector probes provide clear experimental evidence that the SOL conditions depend on the low-Z impurity seeding conditions. These measurements show that in scenarios where neon gas is injected into the plasma, the tungsten divertor leakage and SOL transport depend on the poloidal location from which the neon is injected. In particular, neon injection from the Inner Midplane and Outer Midplane appear to each result in higher divertor leakage by a factor of 2 to 3 compared to cases with neon injection from either the SOL Crown or from the SAS-VW divertor itself.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Reducing the L-H transition power threshold in ITER-similar-shape DIII-D hydrogen plasmas

Recent dedicated DIII-D experiments in low-torque, ITER-similar-shape (ISS) hydrogen plasmas (at a plasma current I p ~ 1.5 MA and ITER-similar edge safety factor q 95 ~ 3.6) show that the L-H transition power threshold P LH can be reduced substantially (~30%) with L-mode helium admixtures n He /n e $\leqslant$ 25%. In the ensuing H-mode, helium ion fractions n He /n H remain below 25%. H-mode normalized pressure and confinement quality are only slightly affected by helium seeding, and Z eff $\leqslant$ 2.15 (including helium and carbon content). The plasmas investigated here are electron-heat dominated, with temperatures T e (0)/T i (0) $\geqslant$ 1 and edge heat flux ratio Q e /Q i (ρ = 0.95) ~ 1.2–1.5. Without mitigation, P LH is higher by a factor of 2–3 in comparison to similar ISS deuterium plasmas. ISS hydrogen plasmas with lower plasma current I p ~ 1 MA (increased edge safety factor q 95 ~ 5.1) exhibit a substantially lower power threshold. This plasma current dependence, also observed previously on ASDEX-U and in JET, is not accounted for by the commonly used 2008 ITPA multi-machine threshold scaling, but could potentially allow H-mode access at marginal heating power during the initial plasma current ramp-up. Attempts to reduce P LH with low-field- and high-field-side hydrogen pellet injection, using 1.7 mm diameter pellets, have not demonstrated a robust threshold reduction, in contrast to successful earlier experiments with larger 2.7 mm pellets. Finally, techniques for reducing P LH are very important for ITER, in particular for accessing H-mode in hydrogen plasmas during the Pre-Fusion Power Operation-1 (PFPO-1) campaign with marginal auxiliary heating power (20–30 MW of ECH).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

E × B flow driven electron temperature bifurcation in a closed slot divertor with ion B × ∇ B away from the X-point in the DIII-D tokamak

An electron temperature bifurcation is observed in the small angle slot divertor, which has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Experiments in the DIII-D tokamak and associated SOLPS-ITER modeling with full drifts find a strong interplay between drifts and divertor geometry on divertor dissipation. The coupling of divertor geometry and drift flows can strongly affect the path towards divertor detachment onset as the plasma density is raised. With the strike point on the inner slanted surface and ion B × ∇B away from the magnetic X-point, bifurcative transitions were observed with sharp decrease of T e towards detachment onset both experimentally and computationally. This differs from the situation for the open divertor where the T e cliff was only observed for ion B × ∇B towards the X-point. SOLPS-ITER modeling with full drifts demonstrates that the magnitude of the E × B drift flow is comparable with the main plasma flow. The reversal of both the poloidal and radial E × B flows near the strike point leads to rapid density accumulation right near the separatrix, which results in bifurcative step transition of divertor conditions with cold plasma across the entire divertor target plate. Furthermore, these results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D

Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. Cryogenically frozen deuterium pellets are injected from the low-field side of the DIII-D tokamak at frequencies lower than the natural ELM frequency, and heat flux is measured by infrared cameras. Ideal MHD pedestal stability calculations show that without pellet injection, these low collisionality pedestals were limited by their current density (peeling-limited) rather than their pressure gradient (ballooning-limited). ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. Cameras observing both heat flux and D-alpha emission often find significant toroidally asymmetric striations in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets. Toroidal asymmetries at the inner strike point are similar between natural and pellet-triggered ELMs, indicating that the reduction in peak heat flux and total fluence at that location is robust for the conditions reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Evaluation of silicon carbide as a divertor armor material in DIII-D H-mode discharges

Silicon carbide (SiC) represents a promising but largely untested plasma-facing material (PFM) for next-step fusion devices. In this work, an analytic mixed-material erosion model is developed by calculating the physical (via SDTrimSP) and chemical (via empirical scalings) sputtering yield from SiC, Si, and C. The Si content in the near-surface SiC layer is predicted to increase during D plasma bombardment due to more efficient physical and chemical sputtering of C relative to Si. Silicon erosion from SiC thereby occurs primarily from sputtering of the enriched Si layer, rather than directly from the SiC itself. SiC coatings on ATJ graphite, manufactured via chemical vapor deposition, were exposed to repeated H-mode plasma discharges in the DIII-D tokamak to test this model. The qualitative trends from analytic modeling are reproduced by the experimental measurements, obtained via spectroscopic inference using the S/XB method. Quantitatively the model slightly under-predicts measured erosion rates, which is attributed to uncertainties in the ion impact angle distribution, as well as the effect of edge-localized modes. After exposure, minimal changes to the macroscopic or microscopic surface morphology of the SiC coatings were observed. Compositional analysis reveals Si enrichment of about 10%, in line with expectations from the erosion model. Extrapolating to a DEMO-type device, an order-of-magnitude decrease in impurity sourcing, and up to a factor of 2 decrease in impurity radiation, is expected with SiC walls, relative to graphite, if low C plasma impurity content can be achieved. Furthermore, these favorable erosion properties motivate further investigations of SiC as a low-Z, non-metallic PFM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗