Search NASA⌕ Search

SEARCH · Search NASA

Results for “burning plasmas”

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

Experiments to explore the influence of pulse shaping at the National Ignition Facility

The shaping of the drive pulse in time is a key tool in the design of fusion experiments that use inertia to confine burning plasmas. Here, it is directly related to the adiabat and compressibility of the DT fuel, and the characteristics of the laser and target that are needed to ignite. With this in mind, we have performed experiments at the National Ignition Facility that test small changes in the shape of the pulse. In contrast to theory, we find implosions at lower adiabats can have reduced yield and areal density. We discuss implications to performance and the mechanism(s) that could be responsible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A comparison of past and present computational methods for shape analysis of double-shell x-ray radiographs

Implosion symmetry is a key requirement in achieving a robust burning plasma in inertial confinement fusion experiments. In double-shell capsule implosions, we are interested in the shape of the inner shell as it pushes on the fuel. Shape analysis is a popular technique for studying said symmetry during implosion. Combinations of filtering and contour-finding algorithms are studied for their promise in reliably recovering Legendre shape coefficients from synthetic radiographs of double-shell capsules with applied levels of noise. A radial lineout max(slope) method when used on an image pre-filtered with non-local means and a variant of the marching squares algorithm are able to recover p 0 , p 2 , and p 4 maxslope Legendre shape coefficients with mean pixel discrepancy errors of 2.81 and 3.06, respectively, for the noisy synthetic radiographs we consider. Here, this improves upon prior radial lineout methods paired with Gaussian filtering, which we show to be unreliable and whose performance is dependent on input parameters that are difficult to estimate.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Progress toward hydro-equivalent ignition in OMEGA direct-drive DT-layered implosions

Considerable progress has been made in deuterium-tritium-layered implosion experiments on the OMEGA Laser System, bringing the prospects for thermonuclear ignition in direct-drive configurations with megajoule-class lasers closer to reality. Doing so has required navigating the balance between improved 1D performance and multidimensional stability. Using statistical modeling based on over 350 cryogenic implosions to identify various degradation mechanisms, and combined with multidimensional simulations and experimental techniques such as target offsets to combat residual flows, core conditions have repeatably been achieved that extrapolate to the burning-plasma state when scaled to 2.15 MJ of symmetric laser illumination. Using high implosion velocities (⁠> 450 km/s) and moderately high adiabats (⁠~5⁠), these experiments produced record-high scaled Lawson parameters in direct drive equal to 89 ± 2% of that required for ignition with expected yields of up to 1.5 ± 0.2 MJ. To improve these results still further, focused physics studies are performed to improve physics understanding and identify routes to even greater performance. Recent studies include investigations into the impact of mounting features, laser imprint, reduced fuel temperatures, and greater on-target intensities through subscale experiments. This manuscript gives a summary of the cryogenic direct-drive program on the OMEGA laser, including routes taken to achieve the current best performance, the status of recent focused physics investigations, and future designs—such as target solutions to laser imprint and reducing vapor density to increase convergence—that are expected lead to the demonstration of hydro-equivalent ignition on OMEGA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗
DOE OSTI2022

A New Integrated Analysis Suite for Fast-Ion Study in KSTAR

Here, an integrated workflow for fast-ion analysis was developed by adapting the One Modeling Framework for Integrated Task (OMFIT) workflow manager to support a standard and unified analysis platform for KSTAR users. The newly established analysis suite offers a graphical user interface–based workflow to enable users to readily access and handle experimental data archived in various data formats and servers. Further, users can analyze the data by importing modules designed for conducting certain tasks, such as profile fitting, equilibrium reconstruction, and postprocessing of tokamak data. The procedures for preparing the inputs for fast-ion simulations are streamlined by a common workflow manager, which enables the parallel processing of various tasks to efficiently analyze large fast-ion datasets. The OMFIT platform comprises a flexible Python-based application that enables users to freely manipulate the Python scripts for applications that are unavailable in the standard workflow. The framework also offers mapping tools to translate the output data into the Integrated Modeling and Analysis Suite format to maintain application compatibility for future ITER burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Systematic Source Determination and X-Ray Radiography Detection of Nonuniformities on High-Density Carbon Ablators

Recent deuterium-tritium (D-T)–layered implosion experiments at the National Ignition Facility have achieved a burning plasma and >1-MJ neutron yield. A series of repeat experiments have shown that the degree of performance is very likely dependent on capsule quality, including the quantity of what are collectively termed “high- Z particles.” These particles are detected on a custom-built radiography system, known as the Sagometer, during the final target qualification process. The term particles is misleading, as the source of these nonuniformities in the capsule images is uncertain; the term detection will be used instead. An increased number of D-T targets have been rejected at the final stages of production due to Sagometer detections. Late detections are deleterious in terms of loss of production parts, effort, and overall operating efficiency. In response, we undertook an effort to determine the origin of these detections and to ultimately mitigate target losses caused by them. Through careful testing and analysis, we have determined neither insufficient production cleanliness nor hohlraum shedding is responsible for the detections on the capsule. Here, we determined that the detections are inherent to the capsule and have made efforts to use the Zeiss Xradia to identify them earlier in the production process. While testing revealed the Xradia is not currently sufficient for identifying such particles using radiography images, we continue to look to other forms of metrology to down select the capsules early in the process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

The Impact of Disruptions on the Economics of a Tokamak Power Plant

Tokamaks are often considered to be a leading candidate for near-term, cost-effective fusion energy, but these devices are susceptible to sudden loss of confinement events called disruptions. The threat of disruptions has garnered serious attention in research for the next generation of burning plasma experiments, such as ITER, but has received little treatment in economic studies of magnetic fusion energy. In this paper, we present a model for quantifying the effect of disruptions on the cost of electricity produced by a tokamak power plant (TPP). We outline the various ways disruptions increase costs and decrease revenues, introduce metrics to quantify these effects, and add them to a levelized cost of electricity (LCOE) model. Additionally, we identify several rate-limiting repair steps and introduce a classification system of disruption types based on the time to return to operations. We demonstrate how the LCOE model can be used to find the cost of electricity and the requirements for disruption handling of a TPP, and we further highlight where future research can have a strong impact in neutralizing the “showstopping” potential of disruptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Overview of density pedestal structure: Role fueling versus transport

How much of the density pedestal structure profile is determined by fueling versus transport is still a valid question. In most current tokamaks, the Scrape-Off Layer is not opaque enough to screen the neutrals and limit their penetration past the separatrix. As such recycled neutrals as well as additional gas puffing can directly affect the pedestal structure through ionization. Yet, in ITER and other future burning plasma magnetic confinement devices, this ionization source will be pushed further out into the Scrape-Off Layer. This poses the question, how much of the pedestal density structure is governed by this edge ionization source versus plasma transport effects. Theoretically, several turbulent modes have been identified which could provide a 'pinch' like up-gradient transport mechanism. Up-gradient transport is necessary in a source-free region to obtain a peaked density profile, which is often observed in the core of tokamaks. In an opaque Scrape-Off Layer, without the existence of such a pinch, the pedestal density structure would eventually disappear. Not only would this limit our ability to fuel a fusion reactor, it would also affect the pedestal stability. In this paper we will give an overview of the state of the art on what sets the pedestal density structure based on experimental observations as well as theoretical models. We will complement these studies with new results from Alcator C-Mod and DIII-D where the opaqueness is increased to reach values similar to those in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Follow the power—pathways to steady-state tokamak reactors

Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical investigation of alpha particle confinement under the perturbation of neoclassical tearing modes and toroidal field ripple in CFETR

The confinement of alpha particles in burning plasma is a key issue in fusion reactor design, including particle interaction with instabilities. Here, we include two topics: the effect of neoclassical tearing modes (NTMs) and toroidal field ripple on alpha particle loss, and the assessment of particle redistribution under an NTM with a reduced model. We consider Chinese fusion engineering test reactor parameters, the alpha particle distribution given by TRANSP/NUBEAM and the NTM perturbation function given by the initial value code TM1. We show that the synergistic effect of the NTM and ripple is negligible; the particle loss fraction does not change with increasing NTM amplitude. Only passing particles are affected by the mode particle resonance, producing profile flattening but no increased loss because only trapped particles are influenced by ripple. To study alpha particle profile flattening, the work adopts an innovative method of phase vector rotation to determine regions of good and broken Kolmogorov–Arnold–Moser surfaces and equilibrates the particle density according to local stochasticity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of machine learning systems trained to detect Alfvén eigenmodes using the CO 2 interferometer on DIII-D

Abstract A Machine-Learning (ML) based detection scheme that automatically detects Alfvén Eigenmodes (AE) in a labelled DIII-D database is presented here. Controlling AEs is important for the success of planned burning plasma devices such as ITER, since resonant fast ions can drive AEs unstable and degrade the performance of the plasma or damage the first walls of the machine vessel. Artificial Intelligence could be useful for real-time detection and control of AEs in steady-state plasma scenarios by implementing ML-based models into control algorithms that drive actuators for mitigation of AE impacts. Thus, the objective is to compare differences in performance between using two different recurrent neural network systems (Reservoir Computing Network and Long Short Term Memory Network) and two different representations of the C O 2 phase data (simple and crosspower spectrograms). All C O 2 interferometer chords are used to train both models, but only one is processed during each training step. The results from the model and data comparison show higher performance for the RCN model (True Positive Rate = 90% and False Positive Rate = 14%), and that using simple magnitude spectrograms is sufficient to detect AEs. Also, the vertical C O 2 interferometer chord passing near the center is better for ML-based detection of AEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗