Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical tokamak”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

Editorial for the 28th Workshop on Magnetohydrodynamic (MHD) Stability Control

This special issue contains four papers which were presented at the 28th Workshop on Magnetohydrodymanic (MHD) Stability Control. The workshop was held at Auburn University and on Zoom from 3–5 October 2024. Furthermore, this annual workshop covers the active and passive control of MHD instabilities such as kink, tearing and edge-localized modes in tokamaks, reversed field pinches, spherical tori or other magnetic confinement devices.

Yang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2021

Effect of boronization on plasma-facing graphite surfaces and its correlation with the plasma behavior in NSTX-U

Boronization is a Plasma Facing Component (PFC) conditioning technique widely used in tokamak machines. The National Spherical Torus Experiment-Upgrade (NSTX-U) applied this conditioning, using a plasma glow with a deuterated Trimethyl-boron (d-TMB) and He mixture. The use of boronization during the campaign improved the plasma performance, allowing longer plasma discharges and H-mode access. The chemical state of an ATJ graphite sample, used as a proxy for the NSTX-U PFCs, was monitored in-situ using the Materials Analysis Particle Probe (MAPP) diagnostic and X-ray Photoelectron Spectroscopy (XPS). The XPS data showed a progressive rise (from + fluence increased. Filterscopes were used to measure the light emitted by oxygen impurities in the plasma near the surface of the PFC. An increase in the registered magnitude of the OII line, normalized to the Dγ intensity, was observed as the concentration of O on the ATJ surface increased. The plasma performance was found to be strongly correlated to oxygen impurity concentrations at the plasma edge and on the PFC surface, as measured by the discharge length and access to the H-mode regime. In this work, we present a quantitative analysis of the evolution of the chemistry of the ATJ surface, and the oxygen presence in the plasma-material interface, and report relevant plasma parameters observed during the same period of time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Recycling Measurements for the Lithium Tokamak Experiment

The lithium tokamak experiment (LTX) is a spherical torus designed to investigate a low-recycling wall regime for magnetically confined plasmas. Traditional measurements of recycling using Dα emission are not feasible in LTX due to reflections off the plasma-facing-components: stainless-steel shell, liquid lithium layer, and molybdenum limiters. Since lithium has low VUV reflectivity, LTX will use Lyman-α emission to measure particle recycling. LTX will use 10- and 20-channel AXUV diode arrays with Lyman-α filters to view recycling off the inboard and outboard shells. A single detector views recycling off a limiter. Measurements will be used with neutral modeling to calculate recycling and neutral deposition. To investigate the performance of an AXUV diode with directly-deposited filter, measurements were obtained with this detector and a visible detector viewing a calibration plasma then compared with emission calculated using a collisional-radiative model.

Physics - Plasma physics↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2021

Modeling of small tungsten dust grains in EAST tokamak with NDS-BOUT ++

In order to investigate the transport of small dusts as well as their evolution property along their trajectories, the NDS module is developed under the BOUT++ framework, a highly desirable C++ code package to perform parallel plasma fluid simulations with an arbitrary number of equations in three-dimensional curvilinear coordinates. Due to the severe dust ablation in fusion plasmas, the dust size would decrease from micrometer to nanometer, resulting in impurities. Small dusts in the simulations here are specified as tungsten spheres with the radii on or below the order of submicrometer. The Rayleigh limit is included in the charging process when the dust is ablated to the droplet phase. The simulation results from the NDS module show that a 200 nm radius spherical tungsten dust originated from upper divertor region of EAST Tokamak is ablated completely due to the intense heating from the incoming plasma inside the core region, well consistent with the CCD footage of EAST shot # 81459. Furthermore it is found that the magnetic field dominates the dust transport when the dust radius is below 100 nm during the ablation along the trajectory. Our simulations predict that a 10 nm radius spherical tungsten dust injected from the inner midplane is well constrained by the magnetic field, and it reaches the inner divertor target with a velocity on the order of km/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multi-machine validation of plasma initiation modelling and prospects for future devices: Predicting plasma initiation using only hardware design and control room input data

This paper reports on the generic prediction capability of full electromagnetic plasma initiation modelling with DYON, which was carried out for the first time in fusion research by the joint modelling of the International Tokamak Physics Activity—Integrating Operation Scenario group. The following devices were included in the experiment database: VEST (spherical torus, copper coils, Stainless steel wall, R/a = 0.3 m/0.2 m, V v = 3.7 m 3 ), MAST-U (spherical torus, copper coils, C wall, R/a = 0.7 m/0.5 m, V v = 55 m 3 ), EAST (conventional tokamak, superconducting coils, metallic wall, R/a = 1.85 m/0.5 m, V v = 38 m 3 ), DIII-D (conventional tokamak, copper coils, C wall, R/a = 1.67 m/0.65 m, V v = 35 m 3 ), and KSTAR (conventional tokamak, superconducting coils, C wall, R/a = 1.8 m/0.5 m, V v = 55 m 3 ). Despite the different hardware features of the devices, the required operating spaces of the loop voltage induction and prefill gas pressure for inductive plasma initiation in each device were successfully reproduced by the predictive simulations with DYON using only the individual hardware design and the control room input data for each discharge. This successful validation across multiple machines demonstrates that the full electromagnetic DYON modelling can capture the essential physics of inductive plasma initiation. The simulation settings commonly employed for all modelling and the modifications necessary to account for the discrepancies between individual devices are reported. Predictions for ITER based on the multi-machine validation indicate that a wide range of prefill gas pressures exists for the Townsend breakdown and the plasma burn-through (0.01–1.5 mPa).

DYON↗
DOE OSTI2020

Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods

Here, we present an energy-conserving discontinuous Galerkin scheme for the full-f electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for ∂A II /∂t, the inductive component of the parallel electric field, using a generalized Ohm’s law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. Here, an algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Basic Physical Processes Involving Dust in Fusion Plasmas

This report presents the main outcomes of our research focused on understanding the behavior and effects of dust particles in fusion plasmas, particularly in the edge regions of tokamaks and stellarators. We conducted computer modeling studies of burst injections of carbon and tungsten dust particles in DIII-D like divertor plasmas. The studies investigated effects of transient influx of the low- and high-Z material plasma contaminants in form of dust on the edge plasma dynamics in a modern mid-size tokamak. We also performed theoretical and computational studies of the forces acting on non-spherical dust grains in magnetized and nonmagnetized plasmas. In addition, we cooperated with General Atomics slag management group on development of DIII-D dust injection experiments to measure trajectories of dust grains in the divertor plasma and compare them with DUSTT code predictions for validation of the dust modeling capabilities. We collaborated with JET experimentalists on evaluation of radiative losses induced by injection of mixed neon-deuterium ice pellets for disruption and run-away electron generation mitigation during thermal and current quench phases. We also cooperated with experimentalists at LHD fusion device (Japan) on modeling support of experiments on injection of boron granules in fusion plasma discharges to assess their dynamics and effectiveness for in situ dynamic wall conditioning.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Optimization and active stabilization of a far-infrared laser for NSTX-U high poloidal wavenumber scattering diagnostics

The far-infrared (FIR) laser output beam power and profile are important parameters in the laser-aided diagnostics, directly influencing the spatial resolution and signal-to-noise ratio of measurements. Here, this work focuses on developing a systematic control method to enhance FIR laser beam quality through optimized mirror alignment and real-time feedback-based precision cavity length tuning. A 150 W CO 2 laser, aligned with the waveguide axis using a HeNe reference laser, serves as the pump source. The sensitivity of FIR beam intensity to pump gas pressure and thermal expansion is investigated, revealing that even a 1 µm cavity expansion can significantly degrade output power stability to about two-thirds of its original value. To address this, a feedback control module has been designed and implemented for active cavity length adjustment, stabilizing the output power at ∼30 mW. In addition, maintaining a high formic acid gas pressure ($>$190 mTorr) within the cavity ensures reliable operation. The optimized FIR laser will be deployed on the National Spherical Torus eXperiment-U high poloidal wavenumber scattering system for studying electron-scale turbulence in tokamak plasmas.

Xu, Xinhang [Univ. of California, Davis, CA (Unite↗
DOE OSTI2022

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high- β are reduced to the level of that of collisionless trapped electron modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2021

Two-color scattering for the measurement of neutrals at the edge of fusion devices

Laser two-color scattering (TCS) is proposed to detect the neutral species in the edge of fusion devices, namely, tokamaks. Here, TCS uses two wavelengths to probe both the laser Rayleigh scattering and Thomson scattering of the neutral-electron bath, with emphasis on neutral density measurements such as that of hydrogen and deuterium. Modeling of the Rayleigh scattering of tokamak neutral species under various plasma conditions (electron density and temperature) shows that, with an appropriate filtering of the Thomson signal and by going to ultraviolet-region wavelengths, identification of the Rayleigh signal can be achieved. Photon count and signal fractions are calculated in two test cases, one in the midplane region of the National Spherical Torus Experiment and one in the divertor region of DIII-D. An uncertainty analysis and discussion of the feasibility of the TCS diagnostic is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

HTS Cable Conductor for Compact Fusion Tokamak Solenoids

Significant progress has been made recently in the US fusion community to develop a strategic plan to enable engineering design and construction of a fusion pilot plant (FPP). Princeton Plasma Physics Laboratory (PPPL) is working on developing high current density HTS conductors for next fusion experiments. Partnering with the US industry, we are evaluating feasibility and affordability of Conductor on Round Core (CORC ® ) developed by Advanced Conductor Technologies (ACT) for the next compact fusion tokamak facility. High current density achieved by a CORC ® cable based on a four-layer model coil recently tested at National High Magnetic Field Laboratory (NHMFL) motivated its consideration for low cost, reduced size fusion magnet application. This is of interest to PPPL because of its scalability to tokamak central solenoid (CS) coils in terms of required flux swings for plasma startup operations. Furthermore, partnering with ACT, we designed and built a two-layer model coil solenoid directly wound with CORC to demonstrate its applicability for compact solenoids such as that used in the national spherical torus experiment (NSTX), NSTX-upgrade (NSTX-U) and the US sustained high-power density test facility (SHPD). The ~160 mm diameter solenoid wound by a two-layer CORC is being tested in early 2022 under electromagnetic cyclic loading at NHMFL in a unique 160 mm bore, 14 T background field magnet facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Observation of stationary filaments with resonant magnetic perturbations in NSTX

The application of resonant magnetic perturbations (RMPs) induces striated heat fluxes on the divertor plate in tokamaks. These striated heat fluxes were previously considered strike point splitting by the intersection between the divertor and predicted magnetic lobe structures. Recently, the existence of stationary 3D filaments has been observed for the first time in National Spherical Torus Experiment (NSTX) plasmas with toroidal mode number n = 3 RMPs. Magnetic field line tracing of these filaments from the midplane to the divertor is consistent with the quantity and the radial and toroidal location of striated heat fluxes observed on the divertor. This study suggests that divertor striated heat flux patterns are linked to these stationary filaments rather than to strike point splitting or magnetic lobe structures for NSTX. This finding will help guide future efforts to understand the heat flux striations in NSTX with RMPs in order to control and reduce the divertor peak heat flux.

RMPs↗