Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tokamak edge”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Millimeter-Wave Imaging Technology Development for Real-Time 2D and 3D Fusion Plasma Diagnostics

Including a one year No-Cost Extension, The University of California, Davis has been funded since 2020 to advance innovative millimeter-wave diagnostics for fusion plasma research. This work includes: (1) Developing and applying cutting-edge System-on-Chip millimeter wave Integrated Circuit (IC) technology for Electron Cyclotron Emission Imaging (ECEI) and Microwave Imaging Reflectometry (MIR) on the DIII-D tokamak and the NSTX-U spherical tokamak. (2) Development of radiation-hardened millimeter-wave chip for diagnostics in harsh environments, utilizing wide bandgap Gallium Nitride (GaN) semiconductors. (3) Integrating forward-modeling synthetic diagnostics to enable advanced, real-time plasma prediction and feedback control. Originally planned as a three-year program (2020–2023), the project was granted a one-year no-cost extension, extending its duration through the end of 2024.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Prediction of DIII-D Pedestal Structure from Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. An experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (ne) and electron temperature (Te) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (Ip), toroidal magnetic field (Bφ), neutral beam heating power (PNBI) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 2. Pedestals and flow

Part 1 described a wide range of analytic tokamak equilibria modelling smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to further include edge pedestals and toroidal flow. Specifically, edge pedestals are allowed in the pressure, pressure gradient and toroidal current density. Also, an edge-localized contribution to the bootstrap current is treated. In terms of flow, analytic solutions are obtained for two cases: a $\gamma = 2$ adiabatic and a $\gamma = \infty $ incompressible energy conservation relation.

Physics↗
DOE OSTI2022

Achievements of actively controlled divertor detachment compatible with sustained high confinement core in DIII-D and EAST

The compatibility of efficient divertor detachment with high-performance core plasma is vital to the development of magnetically controlled fusion energy. The joint research on the EAST and DIII-D tokamaks demonstrates successful integration of divertor detachment with excellent core plasma confinement quality, a milestone towards solving the critical plasma–wall-interaction (PWI) issue and core-edge integration for ITER and future reactors. In EAST, actively controlled partial detachment with T et,div ~ 5 eV around the strike point and H 98 > 1 in different H-mode scenarios including the high β P H-mode scenario have been achieved with ITER-like tungsten divertor, by optimizing the detachment access condition and performing detailed experiments for core-edge integration. For active long-pulse detachment feedback control, a 30 s H-mode operation with detachment-control duration being 25 s has been successfully achieved in EAST. DIII-D has achieved actively controlled fully detached divertor with low plasma electron temperature (T et,div ≤ 5 eV across the entire divertor target) and low particle flux (degree of detachment, DoD > 3), simultaneously with very high core performance (β N ~3, β P > 2 and H 98 ~ 1.5) in the high βP scenario being developed for ITER and future reactors. The high-β P high confinement scenario is characterized by an internal transport barrier (ITB) at large radius and a weak edge transport barrier (ETB, or pedestal), which are synergistically self-organized. Both the high-β P scenario and impurity seeding facilitate divertor detachment. The detachment access leads to the reduction of ETB, which facilitates the development of an even stronger ITB at large radius in the high β P scenario. Thus, this strong large radius ITB enables the core confinement improvement during detachment. In conclusion, these significant joint DIII-D and EAST advances on the compatibility of high confinement core and detached divertor show a great potential for achieving a high-performance core plasma suitable for long-pulse operation of fusion reactors with controllable steady-state PWIs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C +6 and D +1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM'ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. Here, the ELMs result in a ~25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comprehensive new insights on the potential use of SiC as plasma-facing materials in future fusion reactors

Abstract The performance of silicon carbide as an alternative plasma facing material (PFM) was studied at various irradiation conditions relevant to ion energies and fluxes of a fusion reactor. This analysis involves detailed modeling of subsurface plasma/material interactions, sputtered particle transport above the surface and redeposition, and related changes in material composition and microstructure induced by steady-state and Edge Localized Mode ion fluxes. Transition of a crystalline SiC surface to semi-crystalline and amorphous phases was analyzed based on advanced modeling of DIII-D tokamak experiments where SiC was irradiated in single- and multiple- L-mode and H-mode discharges. This analysis shows that displacement damage, particle deposition/redeposition, and D accumulation on the SiC divertor surface can lead to significant microstructural changes that result in enhanced sputtering erosion in comparison with the original crystalline material. However, the resulting total net erosion rate for a full-coverage, advanced tokamak, SiC coated divertor may well be acceptably low. Moreover, the C sputtering yield from the evolved SiC surface can be seven times lower than from a pure graphite surface; this would imply significantly reduced tritium co-deposition rates in a D-T tokamak reactor, compared with a pure carbon surface. It was also determined that chemical sputtering of both C and Si should not result in any noticeable effect on the net erosion, for attached plasma regimes. Our results thus show encouraging results overall for use of SiC as a PFM in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Control of Alfvén Eigenmodes by Externally Applied 3D Magnetic Perturbations

The suppression and excitation of Alfvén eigenmodes have been experimentally obtained, for the first time, by means of externally applied 3D perturbative fields with different spatial spectra in a tokamak plasma. The applied perturbation causes an internal fast-ion redistribution that modifies the phase-space gradients responsible for driving the modes, determining, ultimately their existence. Hybrid kinetic-magnetohydrodynamic simulations reveal an edge resonant transport layer activated by the 3D perturbative field as the responsible mechanism for the fast-ion redistribution. Finally, the results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design elements and first data from a new Doppler backscattering system on the MAST-U spherical tokamak

A new Doppler backscattering (DBS) system has been installed and tested on the MAST-U spherical tokamak. It utilizes eight simultaneous fixed frequency probe beams (32.5, 35, 37.5, 40, 42.5, 45, 47.5, and 50 GHz). These frequencies provide a range of radial positions from the edge plasma to the core depending on plasma conditions. The system utilizes a combination of novel features to provide remote control of the probed density wavenumber, the launched polarization (X vs O-mode), and the angle of the launched DBS to match the magnetic field pitch angle. The range of accessible density turbulence wavenumbers (kθ) is reasonably large with normalized wavenumbers kθρs ranging from ≤0.5 to 9 (ion sound gyroradius ρs = 1 cm). This wavenumber range is relevant to a variety of instabilities believed to be important in establishing plasma transport (e.g., ion temperature gradient, trapped electron, electron temperature gradient, micro-tearing, kinetic ballooning modes). The system is specifically designed to address the requirement of density fluctuation wavevector alignment which can significantly reduce the SNR if not accounted for.

Instruments & Instrumentation↗
DOE OSTI2025

Optimization of simulated high-field side lower hybrid current drive coupling using machine learning predictions of scrape-off layer density

Lower hybrid current drive (LHCD) is a potential source of non-inductive off-axis current drive (CD) for tokamaks. Although LHCD has been successfully deployed on a number of tokamaks, it is highly sensitive to the scrape-off layer (SOL) conditions local to the LHCD launcher. Large gaps between the launcher and plasma core, SOL turbulence, or edge density perturbations due to edge-localized modes can hamper CD or cause large reflected power. These coupling issues in part motivated the installation of an LHCD launcher on the high-field side (HFS) of DIII-D. On the HFS, the SOL is less turbulent and more controllable compared to the low-field side. This quiescence may result in more predictable edge conditions and thus a more predictable CD. Here, in this work, HFS SOL reflectometry measurements are predicted from global plasma parameters using machine learning models. The SOL predictions coupled with the full-wave simulation of the LHCD launcher allow for the prediction of reflected power, directivity, and arcing risk before the discharge. Launcher performance is then optimized using multi-objective Bayesian optimization, finding the shot parameters that result in an optimal SOL density that maximizes CD while minimizing the risk of arcing. The predictions and optimizations of LHCD performance are then accelerated using a surrogate model of the full-wave LHCD simulation.

Bayesian optimization↗
DOE OSTI2020

Oxidation of lithium plasma facing components and its effect on plasma performance in the lithium tokamak experiment-β

The characteristics of lithium-coated plasma-facing components (PFCs) have been correlated with plasma performance on the lithium tokamak experiment-β (LTX-β). Previous experiments on LTX showed that the application of lithium to PFCs was needed to achieve higher performance discharges with flat electron temperature profiles and high edge temperatures. Samples that match the LTX-β PFCs were exposed to plasmas before and after PFCs were coated with lithium and transferred under vacuum to a surface analysis station. Measurements using x-ray photoelectron spectroscopy (XPS) revealed that the primary surface constituent was lithium oxide. Earlier XPS analysis of lithium-coated PFCs on LTX was only able to show the presence of surface oxygen. The new XPS data from LTX-β have sufficient resolution to clearly identify lithium compounds for the first time, and enable them to be correlated with how lithium-coated PFCs can reduce impurities and retain hydrogen to reduce recycling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Solving the Grad–Shafranov equation using spectral elements for tokamak equilibrium with toroidal rotation

The Grad–Shafranov equation is solved using spectral elements for tokamak equilibrium with toroidal rotation. The Grad–Shafranov solver builds upon and extends the NIMEQ code (Howell and Sovinec, 2014) previously developed for static tokamak equilibria. Both geometric and algebraic convergence are achieved as the polynomial degree of the spectral-element basis increases. A new analytical solution to the Grad–Shafranov equation is obtained for Solov’ev equilibrium in presence of rigid toroidal rotation, in addition to a previously obtained analytical solution for a different set of equilibrium and rotation profiles. The numerical solutions from the extended NIMEQ are benchmarked with the analytical solutions, with good agreements. Besides, the extended NIMEQ code is benchmarked with the FLOW code (Guazzotto et al., 2004). The modification of pressure profile induced by toroidal flow is investigated. The relative change of pressure profile is found significant around the edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Active Divertor Heat Flux Control using Impurity Powder Dropper

Divertor plasma-facing components (PFCs) in a tokamak are typically designed to withstand average steady-state heat loads of about 5–10 MW/m 2 , a limit that applies to both solid and liquid lithium (LL) PFCs. Exceeding these design values can result in surface damage to tungsten PFCs or excessive lithium (Li) evaporation in liquid lithium divertor (LLD) PFCs. Since exceeding the divertor heat load limits has serious consequences, it is therefore prudent to develop a tool to reduce the divertor heat load and bring the heat load to within the design limit without affecting the plasma performance. Active low Z impurity injection such as Li has been suggested as a potential solution to mitigate excess heat flux as suggested previously, given that non-coronal radiation can be quite large ~ 20–30 MJ per mole of injected Li. Li is considered desirable for reducing the edge neutral recycling helping to improve plasma energy confinement. In this paper, we model the Impurity Power Dropper (IPD) to investigate its potential of divertor heat flux control. The IPD is typically located at the top of the tokamak device and uses a vertical drift tube of a few meters. In the 2 m drift tube case, the IPD powder is accelerated to ~ 6 m/sec before reaching the plasma with the upper divertor configuration, matching the condition for the in-board side pellet injection case. By modeling the IPD geometry we determined the IPD powder deposition profile, and thus the non-coronal radiation and ionization profiles in time as well. From the enhanced radiation power loss, it is therefore possible to reduce the divertor heat load using the divertor simulation code. In conclusion, the IPD divertor heat flux control can be tested in the facilities with IPD including ST-40, DIII-D, EAST, WEST and NSTX-U.

Active lithium injection↗
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
DOE OSTI2021

Observation of fully detached divertor integrated with improved core confinement for tokamak fusion plasmas

Integration of divertor detachment with a high-performance (β N ~ 3, β p > 2, H 98 ~ 1.5) core plasma has been demonstrated in DIII-D high-β p (poloidal beta) plasmas associated with a sustained core internal transport barrier (ITB) and an H-mode edge transport barrier (ETB). Such good core-edge integration has been achieved for both neon and nitrogen seeding, for both favorable and unfavorable B-field directions, independently from the impurity puffing locations, though these variations play important roles on divertor characteristics. Compared to the standard H-mode plasmas, the high-β p plasma exhibits a much wider window of detachment compatible with high confinement core. Fully detached divertor plasmas with low plasma temperature (Te < 5 eV), low particle flux, and low heat flux across the entire divertor target plate were obtained by using nitrogen seeding. This detached high-β p plasma is compatible with a newly developed detachment control system which can help optimize the nitrogen gas flow rate. Several features, i.e., the high edge safety factor in the high-β p scenario, impurity injection, closed divertor and reduced heating power requirement due to the high confinement, facilitate the achievement of full divertor detachment at lower density. Instead of degrading global performance, the divertor detachment facilitates the access to an even stronger ITB at large radius with a relatively weak ETB through self-organized synergy between ITB and ETB, leading to sustained high confinement. The strengthening of the large-radius ITB compensates for the ETB degradation associated with divertor detachment. In addition, a weak ETB naturally has smaller edge localized modes (ELMs). In particular, with neon injection, a long-period no-ELM H-mode phase has been achieved simultaneously with high-performance core and partially detached divertor plasmas. Furthermore, these results demonstrate the possibility of integrating excellent core plasma performance with an effective divertor solution, an essential step toward steady-state operation of reactor-grade plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observation and characterisation of trapped electron modes in Wendelstein 7-X

In the past, quasi coherent (QC) modes were reported for nearly all tokamaks. The general definition describes modes as QC when the magnitude squared coherence is in the range of 0.3 to 0.6. QC modes are observed in the plasma core as well as in the plasma edge and can have quite different physical origins. The one in the core are observed in plasmas with low collisionality, where the electron temperature exceeds the ion temperature in the plasma core. This is the case for electron cyclotron heating in general. The origin of these modes are electrons trapped within a magnetic mirror, as reported in the past from various fusion devices. The so-called trapped-electron modes (TEMs) belong to drift wave instabilities and can be destabilized by electron-temperature gradients in the plasma core. From the diagnostic point of view, QC modes appear as fluctuations in electron density and temperature. Therefore, the microwave reflectometer is very well suited to monitor these modes. This paper describes experiments, conducted at the Wendelstein 7-X stellarator (W7-X), which aim at detecting QC modes at low wave numbers. A poloidal correlation reflectometer installed at W7-X, is able to measure low wave numbers ($k_\perp\unicode{x2A7D}{3.5}\,\textrm{cm}^{-1}$). For medium line-averaged densities ($\int n_\textrm{e}\unicode{x2A7D}{6 \times10^{19}}\,\textrm{m}^{-2}$) the plasma core is accessible for this diagnostic. For different magnetic configurations and plasma parameters, broad QC structures are observed in the coherence spectra. From the analysis of the rotation and the poloidal structure, these QC modes show the properties of electron-temperature-gradient driven TEMs. A linear relation between the mode velocity and the rotation frequency is found. The relation is uniform and confirms the nature of QC-mode observation as TEM in tokamaks, too.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗