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

Stable compression of a spherical tokamak plasma

In this study, an analysis of toroidal plasma compression by a collapsing flux conserver reveals linearly stable scenarios of operation to high compression ratios. The resistive and ideal MHD stability is calculated in full toroidal geometry using the asymptotic matching method in realistic conditions. A time dependent MHD simulation of the compression is conducted to confirm the conservation principles used to calculate the equilibrium states for the analysis. The near edge current profile, controlled by toroidal field ramping during compression, is shown to be critical to stability due to coupling between poloidal components of the least stable mode. Two cases with slight differences in edge current are examined, one with a stable corridor to high compression and one without.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Microtearing stability and turbulence in the pedestal: Linear gyrokinetics, reduced models, and nonlinear turbulent transport

Microtearing modes can play a crucial role in electron heat transport in tokamak plasmas, affecting both energy confinement and overall performance. This study investigates microtearing modes (MTM) stability and turbulence in a JET pedestal through gyrokinetic simulations using the Gene code, complemented by a reduced eigenvalue model. The focus is on how MTM properties depend on key plasma parameters, including collisionality and plasma beta β—the ratio of plasma pressure to magnetic pressure—the normalized toroidal wavenumber k y ρ s ⁠, where ρ s denotes the ion sound gyroradius (typically a few millimeters in edge plasmas) and isotope mass. Collisionality enhances MT growth rates, while increasing β leads to a shift from MTMs to kinetic-ballooning modes, typically for k y ρ s ⁠, where ρ s ≲ 0.2⁠. A purely collisionless branch of MTMs persists at low k y ρ s ⁠, where ρ s with distinctive properties including non-negligible particle flux and ion thermal transport. Isotope mass scans reveal modest reduction of MTM growth rates as ion mass decreases. Nonlinear simulations produce experimentally relevant transport levels. Numerical experiments turning off zonal flows and fields identify the critical role of zonal flows and zonal fields in regulating MTM turbulence. Their removal leads to a significant increase in electron heat flux. These findings provide new insight into MTM-driven transport and its impact on tokamak confinement and lay a foundation for reduced modeling and predictive capabilities.

Electrostatics↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2021

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 1. Limiter and divertor tokamaks

Tokamak equilibria have been derived that are analytic solutions to the Grad–Shafranov equation. Here in this paper, Part 1, describes a wide range of such equilibria including smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to include edge pedestals and toroidal flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Direct Measurement of ICRF-Enhanced Plasma Potentials on WEST Using Reciprocating Emissive Probes

An extensive documentation of ICRF-enhanced plasma potentials has been conducted over two experimental campaigns on the WEST tokamak using reciprocating emissive probes magnetically connected to two ICRF antennas. The collected data spans a wide range of antenna electrical settings (coupled power, toroidal phasing, left–right power balance) and plasma parameters (density at the antenna limiter above and below the lower hybrid resonance, plasma current, minority fraction). By scanning the edge safety factor across multiple probe plunges, the magnetic connection between the probe and the antenna varied, enabling the construction of a 2D map of the plasma and floating potentials around an active ICRF antenna. This dataset will be used to validate RF simulation tools equipped with the sheath boundary condition and used to predict RF rectified potentials and ICRF-induced impurity sputtering in future machines. This paper presents the diagnostic and some initial measurements, while the rest will be reported elsewhere.

Diab, Raymond [Massachusetts Inst. of Technology (↗
DOE OSTI2020

Stabilizing effect of enhanced resistivity on peeling-ballooning instabilities on EAST

Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.

EAST tokamak↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗
DOE OSTI2023

Boundary, SOL, and Divertor Physics Studies on TCV (Final Report for Research Funded by DoE Award DE-SC0020327)

Research supported under this grant was organized around two diagnostic systems installed on the TCV tokamak located at EPFL in Lausanne, Switzerland: Multi-Spectral Imaging (MSI) and Gas-Puff Imaging (GPI). MSI monitors 2-d time-history images of up to 10 spectral lines, in the visible, using identical views for all wavelengths. The primary application of MSI is to infer 2-d time histories of electron temperature and density in the Scrape-Off-Layer (SOL) of diverted plasmas, utilizing ratios of emission from multiple neutral helium lines. The atomic models are validated against 1-d Thomson scattering measurements. SOL density and temperature maps are analyzed to infer divertor dynamics, including particle and power flows, for a variety of divertor configurations, and specifically as flux expansion along the long-legged outboard divertor is systematically varied. GPI is used to measure details of edge fluctuations in the SOL, particularly as a function of plasma triangularity (δ), including both positive and negative δ. A threshold is found, whereby for δ<-0.3, edge turbulence in the far SOL is completely suppressed. Dynamics of the filamentary turbulence was also analyzed using advance machine learning techniques. Research in support of four PhD theses was carried out within the scope of this Award, including 2 at MIT and 2 at EPFL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

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↗
DOE OSTI2024

Versatile multi-energy hard x-ray camera to study confined and unconfined fast electron dynamics and anisotropies (Invited)

A powerful and flexible hard x-ray (HXR) camera has been recently installed and tested on the WEST tokamak (CEA, France) in collaboration with the Princeton Plasma Physics Laboratory. The diagnostic is a pinhole camera fielded with a 2D pixel detector equipped with a 1 mm thick CdTe sensor. The novelty of this diagnostic technique is the detector’s capability of adjusting the threshold energy at the pixel level. This innovation provides great flexibility in the energy configuration, allowing simultaneous space, energy, and time resolved x-ray measurements. The novel camera has been used to measure the core radiation from non-Maxwellian (fast) electrons accelerated by Lower Hybrid (LH) waves and also the beam–target emission of tungsten in the divertor region produced by fast electron losses interacting with the target. In addition, anisotropic hard x-ray emission has been detected for the first time at the WEST core and edge plasma, with opposite toroidal intensity trends. Finally, experimental vertical and toroidal HXR profiles have been successfully reproduced with the LH code LUKE.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Stellarators with enhanced tritium confinement and edge radiation control

A stellarator design is described with the purpose of achieving three goals: (1) enhance the confinement time of tritium. (2) Have a sufficient density of high-Z impurities to radiate the thermal power escaping from the core while having an extremely low impurity density in the core. (3) Maintain a large fraction of the plasma in a burning plasma state with an optimal tritium fraction. Some features of this design could be used in tokamaks. Although having three confinement zones is natural for stellarators, it is not for tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Divertor heat flux challenge and mitigation in SPARC

Owing to its high magnetic field, high power, and compact size, the SPARC experiment will operate with divertor conditions at or above those expected in reactor-class tokamaks. Power exhaust at this scale remains one of the key challenges for practical fusion energy. Based on empirical scalings, the peak unmitigated divertor parallel heat flux is projected to be greater than 10 GW m -2 . This is nearly an order of magnitude higher than has been demonstrated to date. Furthermore, the divertor parallel Edge-Localized Mode (ELM) energy fluence projections (~11–34 MJ m -2 ) are comparable with those for ITER. However, the relatively short pulse length (~25 s pulse, with a ~10 s flat top) provides the opportunity to consider mitigation schemes unsuited to long-pulse devices including ITER and reactors. The baseline scenario for SPARC employs a ~1 Hz strike point sweep to spread the heat flux over a large divertor target surface area to keep tile surface temperatures within tolerable levels without the use of active divertor cooling systems. In addition, SPARC operation presents a unique opportunity to study divertor heat exhaust mitigation at reactor-level plasma densities and power fluxes. Not only will SPARC test the limits of current experimental scalings and serve for benchmarking theoretical models in reactor regimes, it is also being designed to enable the assessment of long-legged and X-point target advanced divertor magnetic configurations. Experimental results from SPARC will be crucial to reducing risk for a fusion pilot plant divertor design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of shape on pedestal characteristics in the wide pedestal quiescent H-mode in the DIII-D tokamak

Experiments on the DIII-D tokamak have advanced the operational limits of wide pedestal quiescent H-mode (WPQH) plasmas towards increased ITER relevance by simultaneously demonstrating well-matched plasma shape and net zero injected torque. Wide pedestal QHmodes are a compelling candidate regime for a future power producing device because they maintain a stationary pedestal without ELMs via additional edge transport. The pedestal is wider than what would be predicted from kinetic ballooning mode physics due to enhanced edge transport generated by broadband turbulence, a limit cycle oscillation, or some combination thereof. Here, compared to the double null shape, the lower single null shape is observed to have a lower density, narrower pedestal width, larger density fluctuations over a broad range of wavenumber, and operates closer to the peeling-ballooning instability boundary calculated from the simple pedestal scaling, Δ ψΝ , which is still observed to be wider than the EPED prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗