Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOL”

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.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

In situ saxs characterization of thermoresponsive behavior of a poly(ethylene glycol)-graft-(poly(vinyl caprolactam)-co-poly(vinyl acetate)) amphiphilic graft copolymer

In this work, we report the thermoresponsive assembly and rheology of an amphiphilic thermosensitive graft copolymer, poly(ethylene glycol)-graft-(poly(vinyl caprolactam)-co-poly(vinyl acetate)) (commercial name Soluplus ® ), which has been investigated for potential biomedical applications. It has received attention due to is ability to solubilize hydrophobic drugs and for its thickening behavior close to body temperature. Through use of the synchrotron at Brookhaven National Lab, and collaboration with the department of energy, the nanoscale structure and properties can be probed in greater detail. Soluplus ® undergoes two structural changes as temperature is increased; the first, a concentration independent change where samples become turbid at 32 °C. Increasing the temperature further causes the formation of physically associated hydrogels. This sol-gel transition is concentration dependent and occurs at 32 °C for 40 wt% samples, and increases to 42 °C for 10 wt% samples. From variable temperature SAXS characterization micelles of 20–25 nm in radius can be seen and maintain their size and packing below 32 °C. A gradual increase in the aggregation of micelles corresponding to a thickening of the material is also observed. Close to and above the gelation temperature, micelles collapse and form a physically associated 3D network. A model is proposed to explain these physical effects, where the poly(vinyl caprolactam) group transitions from the hydrophilic corona at room temperature to the hydrophobic core as temperature is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Plasma flux expansion control on the DIII-D tokamak

A new controller has been developed with help of the flexible divertor poloidal-field coil set of the DIII-D tokamak, to aid in the precise control of the flux expansion in the scrape-off layer (SOL). Here, the single-input multiple-output (SIMO) architecture ensures flexibility through a complementary set of orthogonal actuator direction to guarantee minimum effect on existing controlled variables, e.g. radial and vertical position control of the X-point. A non-linear free-boundary simulation code (GSevolve) is used for simulating the closed-loop response and for verifying the implementation of the control algorithm on the DIII-D plasma control system (PCS). First results of the experimental commissioning of the new controller during 2020 DIII-D campaign are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reversed-direction 2-point modelling applied to divertor conditions in DIII-D *

Abstract A predictive form of the extended 2-point model known as the ‘reverse 2-point model’, Rev2PM, is applied to a range of detachment levels in the open lower divertor of DIII-D, showing that the experimentally measured electron temperature ( T e ) and pressure ( p e ) at the divertor entrance can be calculated within 50% from target measurements, if and only if a posteriori corrections for convective heat flux are included in the model. Unlike the standard 2-point model, the Rev2PM calculates upstream scrape-off layer (SOL) quantities (such as separatrix T e and p e ) from target conditions (such as T e and parallel heat flux), with volumetric power and momentum losses depending solely on target T e . The Rev2PM is tested against a database of DIII-D inter-ELM divertor Thomson scattering measurements, built from a series of 6 MW, 1.3 MA, LSN H-mode discharges with varied main ion density, drift direction, and nitrogen puffing rate. Measured target T e ranged from 0.4–25 eV over this database, and upstream T e ranged from 5–60 eV. Poor agreement is found between upstream measurements and Rev2PM calculations that assume purely conductive parallel heat transport. However, introducing a posteriori corrections to account for convective heat transport brings the Rev2PM calculations within 50% of the measured upstream values across the dataset. These corrections imply that up to 99% of the parallel heat flux is carried by convection in detached conditions in the DIII-D open lower divertor, though further work is required to assess any potential dependencies on device size or divertor closure.

2-point model↗
DOE OSTI2025

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2020

Identity of the JET M-mode and the ASDEX Upgrade I-phase phenomena

An H-mode plasma state free of edge-localized mode (ELM), close to the L-H transition with clear density and temperature pedestal has been observed both at the Joint European Torus (JET) and at the ASDEX Upgrade (AUG) tokamaks usually identified by a low frequency (LFO, 1–2 kHz), m = 1, n = 0 oscillation of the magnetics and the modulation of pedestal profiles. The regime at JET is referred to as M-mode while at AUG as intermediate phase or I-phase. This contribution aims at a comparative analysis of these phenomena in terms of the density and temperature pedestal properties, the magnetic oscillations and symmetries. Lithium beam emission spectroscopy (Li-BES) and reflectometer measurements at JET and AUG show that the M-mode and the I-phase modulates the plasma edge density. A high frequency oscillation of the magnetics and the density at the pedestal is also associated with both the M-mode and the I-phase, and its power is modulated with the LFO frequency. The power modulation happens simultaneously in every Mirnov coil signal where it can be detected. The bursts of the magnetic signals and the density at the pedestal region are followed by the flattening of the density profile, and by a radially outward propagating density pulse in the scrape-off layer (SOL). The analysis of the helium line ratio spectroscopy (He-BES) signals at AUG revealed that the electron temperature is modulated in phase with the density, thus the I-phase modulates the pressure profile gradient. Overall, this analysis gave opportunity to compare Li-BES and He-BES density profiles at different locations suggesting a toroidal and poloidal symmetry of the density modulation. The presented results indicate that the regimes, the AUG I-phase and the JET M-mode, exhibit similar characteristics, which leads to the conclusion that the regimes are likely the same.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

SOLPS analysis of the necessary conditions for detachment cliff

SOLPS modeling including full cross-filed drifts is used, for the first time, to reproduce the step-like detachment onset observed in the DIII-D tokamak previously simulated with UEDGE [Jaervinen A. E. et al., 2018 Phys. Rev. Lett.121 075001], and to further investigate the underlying physics processes. The SOLPS analysis reveals a strong interplay between the E×B drift-driven flows and carbon radiation in DIII-D. The modeling shows that the high confinement H-mode facilitates the formation of such a detachment cliff due to the narrow SOL width and hence steep radial gradients near the separatrix, which enhances the poloidal E×B drift. Furthermore, the SOLPS modeling indicates that carbon radiation also plays a key role in the formation of the detachment cliff in DIII-D without additional impurity seeding. These findings may have significant implications for the control and maintenance of divertor detachment for the steady-state operation of high-power fusion devices.

Physics - Plasma physics↗
DOE OSTI2020

New heat flux model for non-axisymmetric divertor infrared structures

A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First attempt to quantify W7-X island divertor plasma by local experiment-model comparison

As a complement to our recent work, which focused on understanding the basic detachment physics and general experimental and numerical trends observed in the W7-X island divertor, this paper compares EMC3-Eirene simulation results with different local diagnostics, including IR cameras, Langmuir probes, H α -cameras, and Thomson scattering. The main purposes are to (1) justify the simulation setup in the previous work, (2) identify the application limitations of the current EMC3-Eirene model, (3) verify the consistency of different diagnostics, and (4) isolate the main geometric and physical effects that need to be prioritized in further developing the EMC3-Eirene code and improving diagnostic capabilities. It turns out that the current version of the EMC3-Eirene code (without drifts) is not yet able to quantitatively reproduce all selected local measurements simultaneously under the current experimental conditions (in particular, the existence of error fields). Nevertheless, it can be shown that within a reasonable range of variation in magnetic configuration, cross-field transport, and SOL plasma state in the modeling, a region of overlap between the numerical results and the local measurements can be established. More accurate model-experiment comparisons will require clarification of error fields, implementation of drifts in EMC3, and improvement of diagnostic capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High gas throughput SOLPS-ITER simulations extending the ITER database to strong detachment

Abstract SOLPS-ITER simulations performed for Q DT = 10, P SOL = 100 MW burning plasmas on ITER extend the existing database to high values of separatrix averaged neon impurity concentration (⟨ c Ne ⟩ ≈ 6%) and divertor neutral pressure (⟨ p div ⟩ > 25 Pa) in order to determine the heat flux mitigation capability of these scenarios and whether strongly detached states are accessible. In the existing database of ITER simulations, the level of detachment was limited to cases where the integral ion flux to the outer target was greater than 80% of the value at rollover, with the impurity radiation localized near the target. With the possibility of narrow heat flux channels and increased deposited power due to tile shaping, it is important to explore operation at a higher degree of detachment. Two series of simulations were explored to extend the database of SOLPS simulations. By increasing the deuterium and neon puff rates proportionally, the peak divertor energy flux ( q ⊥,max ) is decreased from 5 to 3 MW m −2 while ⟨ p div ⟩ increased from 11 to 27 Pa. By increasing only the neon puff, q ⊥, max can be reduced to <1MW m −2 while ⟨ p div ⟩ is maintained at ∼ 11 Pa. As the neon puff level is increased, the position of the impurity radiation peak is shifted towards the X-point. At the highest neon puff levels with steady-state solutions, the electron temperature is reduced below 1 eV across 50 cm of each divertor target. The new cases extend previously observed tight relationships in power and momentum loss factors to low electron temperature improving their utility for highly detached regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: A. identifying options for the control parameters*

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter-erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding. This part A of the present study reports new results for largely the same ITER cases, confirming the strong correlations reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne [Ne/s], and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In the companion report, part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction 2 point modelling, Rev2PM. The input to the reduced model is a value of the variable pair (T et @q ⟂,pk ,Inj Ne ) and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , n e,Xpt (the electron density at the poloidal location of the X-point), etc. In part B the reduced model is quantitatively characterized using one half of the code cases; it is then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half of the cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation and analysis of lower-hybrid-current-drive density limit in EAST

Lower hybrid current drive (LHCD) experiments with line-averaged density up to ~5.1 × 10 19 m –3 were performed in EAST L-mode plasmas. When the line-averaged density rises above a critical value, the hard x-ray (HXR) emission falls to the noise level, indicating that the LHCD density limit is encountered. The experimental results show that the LHCD density limit can be increased with higher wave source frequency (ƒ 0 ) and higher magnetic field (B t ). Although a higher LHCD density limit is obtained by a higher magnetic field for both 2.45 GHz and 4.6 GHz waves, the results show a stronger dependence on the magnetic field for the 4.6 GHz case. Analysis suggests that, for normal operation with a relatively low magnetic field (1.6 T $\leqslant$ B t $\leqslant$ 2.5 T) on EAST, the dominant mechanisms responsible for the LHCD density limit are different between the 2.45 GHz and 4.6 GHz waves. The wave accessibility plays a more significant role during 4.6 GHz LHCD experiments, while parasitic losses due to parametric decay instability (PDIs) dominate the accessibility issue in the 2.45 GHz case. Collisional loss in the scrape-off layer (SOL) may explain the 4.6 GHz result when combined with the accessibility limit at high density and low temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental study on boron distribution and transport at plasma-facing components during impurity powder dropping in the Large Helical Device

Toward real-time wall conditioning, impurity powder dropping experiments with boron powder were performed in the 22nd experimental campaign of the Large Helical Device. To examine the deposition and desorption process of boron, we focus on boron hydride (BH) molecules which presumably populate near plasma-facing components. We performed spatially-resolved spectroscopic measurements of emission by boron ions and BH molecules. From the measurement, we found that BH and B + were concentrated on the divertor viewing chord, which suggest boron deposition in the divertor region. By comparing Hγ emissions with and without boron injection, neutral hydrogen shows uniform reduction in the SOL region, whereas less reduction of neutral hydrogen is confirmed in the divertor region. Although emissions from BH and B + increased linearly, emissions by B0 and B4 + became constant after the middle of the discharge. Continuous reduction of carbon density in the core plasma was confirmed even after B 0 and B 4 + became constant. The results may show reduction of hydrogen recycling and facilitation of impurity gettering by boron in the divertor region and thus effective real-time wall conditioning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bifurcation-like transition of divertor conditions induced by X-point radiation in KSTAR L-mode plasmas *

Abstract Density ramps with ion grad B drift directed into lower single null KSTAR L-mode plasmas are associated with a simultaneous and abrupt reduction of the divertor particle flux on both low- and high-field-side targets when the mid-plane line averaged electron density reaches a given level. Target embedded Langmuir probe signals show a clear ‘cliff edge’ behavior similar to that observed in the divertor target electron temperature in DIII-D H-mode plasmas (Eldon et al 2017 Nucl. Fusion 57 066039; McLean et al 2015 J. Nucl. Mater. 463 533–6). The collapse of the particle flux is observed along the whole divertor target area (from private flux region to the far scrape-off layer (SOL)). The critical upstream density of this target flux cliff is invariant under fuel gas throughput modulation. The transition along the cliff occurs in tens of milliseconds. With the cliff, carbon impurities and deuterium neutrals transported through the X-point to the core produce a strong radiation spot near the X-point, seen on bolometric signals, and increase the upstream density. The experimental observations are consistent with time-dependent SOLPS-ITER simulations, which also demonstrate an abrupt transition of the target flux and upstream density with the increase in X-point radiation. The timescale of the cliff predicted by SOLPS-ITER is consistent with the experiment, although, it is influenced by gas throughput or time-dependent numerical methods. In the L-mode phase space of separatrix electron density and temperature, branches are divided based on target temperature, because the latter is strongly coupled to the radiation front and ionization front due to the monotonic characteristic of the parallel electron temperature distribution. Since the H-mode condition operates at a much higher upstream density and electron temperature in phase space, dissipation from sputtered carbon alone leads to the density limit before reaching the X-point radiation condition. This is therefore consistent with the fact that cliffs have never been observed in H-mode KSTAR experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗