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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗
DOE OSTI2024

The formation of an radial edge electric field due to finite ion orbit width effects is the possible root cause of the H-mode edge

Abstract Full orbit-following simulations of thermal ions show that finite ion-orbit width effects create charge separation near the last closed flux surface (LCFS) which generates a localized radial electric field. Experimentally, edge electric fields are observed in H-mode plasmas and they are necessary for the edge turbulence suppression via the E × B flow shear mechanism. Confined trapped (and to a lesser extent co-passing) ions near the plasma edge form a positive charge distribution outside the LCFS, while thermal electrons are tied more tightly to field lines owing to their small mass and are poorly confined outside the LCFS, hence charge neutrality is violated outside the LCFS. A large number of reported observations from spherical and conventional tokamaks support the results from the simulations although the simulations were not performed fully self consistently. The results suggest ways to lower the H-mode power threshold and optimize the H-mode plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-dependent modeling of dust outburst into tokamak divertor plasma

The first self-consistent simulation studies of the impact of transient injection of tungsten dust on ITER-like edge plasma are presented. The simulations are performed using DUSTT and UEDGE codes coupled in a fully time-dependent manner. The evolution of the divertor and upstream plasma conditions, following the injection of tungsten dust of various quantities and sizes at the outer strike point, is investigated. The estimates of the core plasma impurity fraction and the core impurity screening factor associated with the dust injection are also obtained. Here, the simulations show that the bursts of even modest amounts of tungsten dust injected into the divertor plasmas can result in large perturbations of the divertor power load and unacceptable levels of the core plasma contamination with tungsten impurities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗
DOE OSTI2024

How turbulence spreading improves power handling in quiescent high confinement fusion plasmas

Viable magnetic fusion devices necessitate combining good confinement with effective power flux handling. A major concern for ITER, and devices beyond, is the divertor heat load width, which sets peak boundary heat loads on the plasma-facing materials. Current estimates of the heat flux width are narrow for future reactors. Here, we demonstrate how pedestal turbulence can expand into, or entrain, the stable scrape-off-layer and so broaden the heat flux width beyond these neoclassical predictions. Employing combined theoretical, computational, and experimental approaches, we focus on quiescent high confinement discharges on the DIII-D tokamak, but the results are of broader significance. Our findings uncover common trends in the edge turbulence intensity flux, the pressure perturbation skewness, and the turbulence mixing length, which together determine the heat flux width. This research demonstrates the physics of scrape-off-layer broadening by turbulence and highlights the promise of a turbulent pedestal for successful core-edge integration in ITER and future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First observation of RMP ELM mitigation on MAST Upgrade

Abstract The first experimental attempts at controlling edge localised modes (ELMs) via the application of resonant magnetic perturbations on the MAST Upgrade tokamak are reported. Using the linear MHD model MARS-F, the phase shift between the upper and lower coil rows ΔΦ was optimised for toroidal mode number n = 1 and n = 2 fields, to provide forward guidance to experiments. In low β N discharges, the application of n = 1 3D fields caused the ELM frequency f E L M to increase by over a factor 20 relative to the reference, and also induced a locked mode, which did not cause a plasma termination nor an H-L back transition. However when β N was raised, this induced locked mode caused plasma termination which precluded mitigation access. Initially, applying a numerically optimised n = 2 field had no effect. However applying a rigid toroidal shift to this field caused a locked mode disruption, demonstrating the presence of a substantial n = 2 error field. Coil current ramps were conducted with ΔΦ set at 6 different values, resulting in either locked mode disruptions or no effect, but mitigation with n = 2 fields was not established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Plasma Facing Component Characterization and Correlation With Plasma Conditions in Lithium Tokamak Experiment-β

Lithium coatings in the Lithium Tokamak eXperiment (LTX) led to flat temperature profiles. The flat temperature profiles were observed along with a hot, low density edge, implying a broad, collisionless scrape-off layer (SOL). Additionally, in vacuo X-ray photoelectron spectroscopy (XPS) measurements established that lithium coatings evaporatively deposited onto high-Z plasma facing components (PFCs) became oxidized while retaining the ability to achieve good plasma performance long after lithium was applied to the PFCs. Longstanding theory predicted flat temperature profiles with low recycling walls, which was presumed to be due to hydrogen binding with elemental lithium to form lithium hydride. The presence of oxidized lithium, however, raised questions regarding the exact mechanism of hydrogen retention in LTX. To investigate these questions, the upgraded facility LTX- β includes a new sample exposure probe (SEP) for more detailed in vacuo analysis of PFC samples. The SEP is equipped with a vacuum suitcase capable of transporting samples representative of the LTX- β outer midplane PFCs to a stand-alone XPS system while maintaining pressures lower than the LTX- β base vacuum to limit the contamination between sample exposure and analysis. The low-energy resolution XPS system used in past experiments could only enable the determination of elemental percentages on the PFC sample surfaces. Because the new XPS system has higher energy resolution, it is more direct to assign chemical compounds to the measured binding energies. This capability has been confirmed by comparing XPS data from PFC test samples with measurements using a commercial high-resolution XPS system. Quartz crystal microbalances (QCMs) were used to quantify the thickness of the deposited lithium on the LTX- β PFCs. This article describes the application of the SEP to characterize the PFC surfaces using XPS and their relationship to plasma conditions.

Maan, Anurag↗
DOE OSTI2024

Non-linear dependence of ion heat flux on plasma density at the L–H transition of JET NBI-heated deuterium–tritium plasmas

Abstract Recent JET D–T campaigns opened the possibility of unique isotope studies to investigate the L–H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L–H transitions in the low and high-density branches of JET NBI-heated D–T plasmas. As discussed in the paper, L–H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux ( Q i ). We present the first power balance analysis of JET NBI-heated D–T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D–T Q i , with a lower Q i in D–T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L–H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Q i deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Q i and the power threshold minimum at JET. At low plasma density, NBI power dominates Q i , while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D–T NBI-heated plasmas at JET.

Vincenzi, P. (ORCID:0000000251564354)↗
DOE OSTI2024

Design and analysis of actively-cooled, edge-transport diagnostic for long-pulsed operation in WEST

Next step fusion devices that will operate in steady-state will require complex plasma-facing components (PFCs) that can survive the harsh environment over long timescales not common in current devices. This will require robust plasma facing surfaces that are integrated with active cooling systems. In a collaboration between CEA and ORNL, a plasma-interacting diagnostic is being designed for the W Environment in Steady-state Tokamak (WEST) in Cadarache, France which requires plasma-facing protection like those needed for steady-state PFCs. This integrated diagnostic studies edge transport and impurity migration within WEST, and is planned to include imbedded temperature Langmuir probe sensors, as well as removeable sample slots for ex-situ surface analysis of plasma-material interactions. The entire assembly is expected to move into the plasma edge for periods up to 1000 s having an energy removeable capability of ∼6 kW and seeing a peak heat flux of more than 7 MW/m2. The assembly compliments WEST high fluence campaigns that plan for multiple 1000 s pulses. Within these specifications, the assembly will require a refractory metal plasma facing surface and integral cooling in order to function within the limited space allotted for such diagnostics. Because of the limited space and linear actuator needs, additive manufacturing of the high heat flux working end of the assembly is being considered which could allow for precision cooling-channels and lighter weight designs. Conceptual design along with simulation and analysis results will be presented for this complex diagnostic with novel PFCs.

Lumsdaine, Arnold↗
DOE OSTI2022

Two-dimensional plasma density evolution local to the inversion layer during sawtooth crash events using Beam Emission Spectroscopy

Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2022

Local transport dynamics of cold pulses in tokamak plasmas

Abstract For over two decades, our fundamental understanding of energy transport dynamics in the core of tokamak plasmas had been challenged by the striking observation of temperature perturbation reversals following the injection of cold pulses at the plasma edge. These phenomena were first discovered by Gentle et al. (Phys. Rev. Lett. 74(18):3620–3623, 1995) in 1995 and had long been suggested to be evidence of nonlocal transport effects. In recent years, a new explanation to these phenomena has emerged, fully consistent with the theory of turbulent transport in magnetized plasmas and in remarkable agreement with experiment. This article reviews the experimental observation of temperature reversals in tokamak plasmas and presents the explanation based on local transport physics.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗
DOE OSTI2021

A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments

The interaction of radio-frequency (RF) waves with edge turbulence modifies the incident wave spectrum, and can significantly affect RF heating and current drive in tokamaks. Previous lower hybrid (LH) scattering models have either used the weak-turbulence approximation, or treated more realistic, filamentary turbulence in the ray tracing limit. In this study, a new model is introduced which retains full-wave effects of RF scattering in filamentary turbulence. First, a Mie-scattering technique models the interaction of an incident wave with a single Gaussian filament. Next, an effective differential scattering width is derived for a statistical ensemble of filaments. Lastly, a Markov chain solves for the transmitted wave spectrum in slab geometry. This model is applied to LH launching for current drive. The resulting wave spectrum is asymmetrically broadened in angular wavenumber space. This asymmetry is not accounted for in previous LH scattering models. The modified wave spectrum is coupled to a ray tracing/Fokker–Planck solver (GENRAY/CQL3D) to study its impact on current drive. The resulting current profile is greatly altered, and there is significant increase in the on-axis current and decrease in the off-axis peaks. This is attributed to a portion of the modified wave spectrum that is strongly dampened on-axis during the first pass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the experimental advanced superconducting tokamak (EAST) device with tungsten plasma-facing components. The maximum ELM size was reduced by ~70% in high β N H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet Futatani et al (2014 Nucl. Fusion 54 073008). This study extends previous ELM elimination with Li powder injection Maingi et al (2018 Nucl. Fusion 58 024003) in EAST because (1) use of small, dust-like powder and the related potential health hazards were eliminated, and (2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~10 m s -1 in EAST.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma performance and operational space with an RMP-ELM suppressed edge

Abstract The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3 × 10 19 m −3 . This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The advanced tokamak path to a compact net electric fusion pilot plant

Abstract Physics-based simulations project a compact net electric fusion pilot plant with a nuclear testing mission is possible at modest scale based on the advanced tokamak concept, and identify key parameters for its optimization. These utilize a new integrated 1.5D core-edge approach for whole device modeling to predict performance by self-consistently applying transport, pedestal and current drive models to converge fully non-inductive stationary solutions, predicting profiles and energy confinement for a given density. This physics-based approach leads to new insights and understanding of reactor optimization. In particular, the levering role of high plasma density is identified, which raises fusion performance and self-driven ‘bootstrap currents’, to reduce current drive demands and enable high pressure with net electricity at a compact scale. Solutions at 6–7 T, ∼4 m radius and 200 MW net electricity are identified with margins and trade-offs possible between parameters. Current drive comes from neutral beam and ultra-high harmonic (helicon) fast wave, though other advanced approaches are not ruled out. The resulting low recirculating power in a double null configuration leads to a divertor heat flux challenge that is comparable to ITER, though reactor solutions may require more dissipation. Strong H-mode access (x2 margin over L–H transition scalings) and ITER-like heat fluxes are maintained with ∼20%–60% core radiation, though effects on confinement need further analysis. Neutron wall loadings appear tolerable. The approach would benefit from high temperature superconductors, as higher fields would increase performance margins while potential for demountability may facilitate nuclear testing. However, solutions are possible with conventional superconductors. An advanced load sharing and reactive bucking approach in the device centerpost region provides improved mechanical stress handling. The prospect of an affordable test device which could close the loop on net-electric production and conduct essential nuclear materials and breeding research is compelling, motivating research to validate the techniques and models employed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗