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 91 records · Page 5

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hotspot parameter scaling with velocity and yield for high-adiabat layered implosions at the National Ignition Facility

Here, in this paper, we present a study on hotspot parameters in indirect-drive, inertially confined fusion implosions as they proceed through the self-heating regime. The implosions with increasing nuclear yield reach the burning-plasma regime, hotspot ignition, and finally propagating burn and ignition. These implosions span a wide range of alpha heating from a yield amplification of 1.7–2.5. We show that the hotspot parameters are explicitly dependent on both yield and velocity and that by fitting to both of these quantities the hotspot parameters can be fit with a single power law in velocity. The yield scaling also enables the hotspot parameters extrapolation to higher yields. This is important as various degradation mechanisms can occur on a given implosion at fixed implosion velocity which can have a large impact on both yield and the hotspot parameters. The yield scaling also enables the experimental dependence of the hotspot parameters on yield amplification to be determined. The implosions reported have resulted in the highest yield (1.73×10 16 ±2.6%), yield amplification, pressure, and implosion velocity yet reported at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗

A laser inverse Compton scattering diagnostic to study runaway electron dynamics during tokamak disruptions

The goal is to design a new diagnostic based on laser inverse Compton scattering (LICS) to study the dynamics of runaway electron formation and their subsequent loss during triggered disruptions in DIII-D. It is known that better measurements of the runaway electron population, throughout the evolution of a disruption, is critical for mitigating the impact of energetic electron beams especially in high plasma current discharges such as those in reactor relevant burning plasmas. The goal of the LICS approach being designed for DIII-D is to measure the timeresolved spatial profile and energy distribution of the runaway electrons while they are in the core of the tokamak plasma. During the last year, a code was developed to simulate the signals from a new laser scattering diagnostic for DIII-D, which will measure the relativistic electron population expected in DIII-D after disruptions triggered by argon pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inertial Confinement Fusion—Experimental Physics: Laser Drive

Laser-driven inertial confinement fusion creates a central hot-spot plasma by irradiating a spherical target containing a thermonuclear fuel layer of deuterium (D) and tritium (T) with temporally-shaped, high-intensity lasers or X rays. Inertial confinement fusion relies on the DT-fusion alpha particles depositing their energy in a hot-spot plasma, causing its temperature to rise sharply and a thermonuclear burn wave to propagate out through the surrounding cold, dense DT fuel, producing significantly more energy than was used to heat and compress the fuel. Alpha heating has been demonstrated and burning plasma (i.e., yield amplification due to alpha heating >3.5) is being explored with megajoule-scale lasers on the path to ignition (i.e., yield amplification due to alpha heating >20)

alpha particles↗

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of compact tokamak fusion reactor use cases to inform future transport studies

The OMFIT STEP (Meneghini et al. , Nucl. Fusion , vol. 10, 2020, p. 1088) workflow has been used to develop inductive and steady-state H-mode core plasma scenario use cases for a $B_0 = 8 \, {\rm T}$ , $R_0 = 4 \, {\rm m}$ machine to help guide and inform future higher-fidelity studies of core transport and confinement in compact tokamak reactors. Both use cases are designed to produce 200 MW or more of net electric power in an up-down symmetric plasma with minor radius $a = 1.4 \, {\rm m}$ , elongation $\kappa = 2.0$ , triangularity $\delta = 0.5$ and effective charge $Z_{{\rm eff}} \simeq 2$ . Additional considerations based on the need for compatibility of the core with reactor-relevant power exhaust solutions and external actuators were used to guide and constrain the use case development. An extensive characterization of core transport in both scenarios is presented, the most important feature of which is the extreme sensitivity of the results to the quantitative stiffness level of the transport model used as well as the predicted critical gradients. This sensitivity is shown to arise from different levels of transport stiffness exhibited by the models, combined with the gyroBohm-normalized fluxes of the predictions being an order of magnitude larger than other H-mode plasmas. Additionally, it is shown that although heating in both plasmas is predominantly to the electrons and collisionality is low, the plasmas remain sufficiently well coupled for the ions to carry a significant fraction of the thermal transport. As neoclassical transport is negligible in these conditions, this situation inherently requires long-wavelength ion gyroradius-scale turbulence to be the dominant transport mechanism in both plasmas. These results are combined with other basic considerations to propose a simple heuristic model of transport in reactor-relevant plasmas, along with simple metrics to quantify coupling and core transport properties across burning and non-burning plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polarized fusion and potential in situ tests of fuel polarization survival in a tokamak plasma

Abstract The use of spin-polarized fusion fuels would provide a significant boost towards the ignition of a burning plasma. The cross section for D + T → α + n, would be increased by 1.5 if the fuels were injected with parallel polarization. Furthermore, our simulations demonstrate additional non-linear power gains in large-scale machines such as ITER, due to increased alpha heating. Such benefits require the survival of spin polarizations for periods comparable to the particle confinement time. During the 1980s, calculations predicted that polarizations could survive a plasma environment, although concerns persisted regarding the cumulative impacts of wall recycling. In that era, technical challenges prevented direct tests and left the large scale fueling of a power reactor beyond reach. Over the last decades, this situation has changed dramatically. Detailed simulations of ITER have predicted negligible wall recycling in a high-power reactor, and recent advances in laser-driven sources project the capability of producing large quantities of ∼100% polarized D and T. The remaining crucial step is an in-situ demonstration of polarization survival in a plasma. For this, we outline a measurement strategy using the isospin-mirror reaction, D + 3 He → α + p. Polarized 3 He avoids the complexities of handling tritium, while encompassing the same spin-physics. We evaluate two methods of delivering deuterium, using dynamically polarized Lithium-Deuteride (with vector polarization P V D of 70%) or frozen-spin Hydrogen-Deuteride (with P V D of 40%), together with a method of injecting optically-pumped 3 He (with 65% polarization). Pellets of these materials all have long polarization decay times (∼6 min for LiD at 2 K, ∼2 months for HD at 2 K, and ∼3 d for 3 He at 77 K), all far greater than a plasma shot in a research tokamak such as DIII-D (∼20 s). Both species can be propelled from a single cryogenic injection gun. We review plasma requirements and strategies for detecting polarization survival. Polarization alters both fusion yields and the angular distribution of fusion products, and each of these provides a potential signal. In this paper we simulate a selection of shots with similar characteristics in a future high-T ion H plasma, and find ratios of yields from shots with fuel spins parallel and antiparallel reaching 1.3 (HD + 3 He) to 1.6 (LiD + 3 He) over a wide range of poloidal angles. (A companion paper finds sensitivity to fusion product angular distributions as reflected in the pitch angles of protons and alphas reaching the plasma facing wall.)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗

Measuring Fluxes Of Heat To A Plasma-Arc Anode

Three probes constructed to provide measurements indicative of conductive, convective, and radiative transfer of heat from free-burning plasma arc to water-cooled copper anode used in generating arc. Each probe consists mainly of copper body with two thermocouples embedded at locations 4 mm apart along length. Thermocouples provide measure of rate of conduction of heat along probe and transfers of heat from plasma to sensing surface at tip of probe. Probes identical except sensing surface of one uncoated and other two coated with different materials to make them sensitive to different components of overall flux of heat.

Sankovic, John M.↗

Evaluation of ITER divertor shunts as a synthetic diagnostic for detachment control

Reliable diagnostics that measure the detached state of the ITER divertor plasma will be necessary to control heat flux to the divertor targets during steady state, burning plasma operation. This paper conducts an initial exploration into the feasibility of the divertor shunt diagnostic as a lightweight, robust, and real-time detachment sensor. This diagnostic is a set of shunt lead pairs that measure the voltage drop along the divertor cassette body, from which the plasma scrape-off layer (SOL) current is calculated. Using SOLPS-ITER simulations for control-relevant ITER plasma scenarios, the thermoelectric current magnitude along the SOL is shown to decrease significantly with the onset of partial detachment at the outer divertor target. Electromagnetic modelling of a simplified divertor cassette is used to develop a control-oriented inductance-resistance circuit model, from which SOL currents can be calculated from shunt pair voltage measurements. The sensitivity and frequency-response of the resulting system indicates that the diagnostic will accurately measure SOL thermoelectric currents during ITER operation. These currents will be a good measure of the detached state of the divertor plasma, making the divertor shunt diagnostic a potentially extremely valuable and physically robust sensor for real-time detachment control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding NBI heating and fueling in LTX-β tokamak

My work on the contract DE-SC=0023274 (July, 2022 - July, 2025) was focused on understanding the plasma fueling by NBI. On Dec. 1998 with Sergei Krasheninnikov (UCSD, San Diego, CA) we initiated the Li Wall Fusion (LiWFusion) as a new concept of magnetic fusion [1]. The new 1 2 plasma 3 4 5 6 2.1 Overview of rate coefficients, mean-free paths λH0 , and diffusions coefficients 6 2. GSV code for ⟨σv⟩ analysis of NBI fueling concept was a reaction of the lack of luck of tokamak fusion with QDT = 1 on TFTR and JET. We recognized the edge plasma cooling by recycling to be the route reason of the tokamak problems on the way to burning plasma. LiWFusion relies on plasma pumping by a lithium layer on the inner walls of the plasma chamber, combined with the plasma heating and fueling by the Neutral Beam Injection (NBI). These two innovations eliminate the route problem of the tokmak fusion. The concept became theoretically mature in 2006. In 2012, the technology of continuously Flowing Liquid Lithium (24/7-FLiLi) was invented by me for the future implementations of the LWFusion concept. At the same time, my every presentation on LiWFusion to International Symposium on Lithium Applications since 2010 was objected by some person with words “Everybody since the 1970s knows that plasma fueling by NBI is impossible”. I dropped the name of the author of this objection and ignored his views. My work on the current grant on NBI fueling of LTX-β not only clarified the issue but resolved it in an astonishing way. The NBI fueling was fully understood (thus humiliating the current dogmatic fusion community). The new future of LTX-β dedicated for decades to Li in tokamaks, as well as of the entire magnetic fusion program was envisioned, in sharp contrast with the fallure of OFES in the post-TFTR era of 21st century.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmark between antenna code TOPICA, RAPLICASOL and Petra-M for the ICRH ITER antenna

ITER will be equipped with three plasma heating systems: neutral beam (NB), electron cyclotron (EC), and ion cy-clotron resonance heating (ICRH). The latter consists of two identical ICRH antennas to deliver 20 MW to the plasma (baseline, upgradable to 40 MW). ICRH will play a crucial role in the ignition and sustainment of burning plasmas in ITER. A high fidelity and robust modeling effort to understand the interaction of the IC waves with the scrape-off-layer (SOL) plasma is a very important aspect. Among the main important research topics, we have the assessment of the antenna loading for different plasma scenarios, the role of the lower hybrid resonance in front of the antenna and how to include it in our models, and the RF sheath boundary conditions to evaluate the antenna impurity generation. In this work, we tackle the first of these by reporting on ICRF simulations employing the Petra-M code, which is an electromagnetic simulation tool for modeling RF wave propagation based on MFEM [http://mfem.org] for the ITER ICRH antenna. Moreover, a benchmark between the well tested antenna codes TOPICA, RAPLI-CASOL, which is based on COMSOL [www.comsol.com], and the Petra-M code is also presented. In conclusion, S- and Z-matrices and wave electric field are compared showing an excellent agreement among these codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗