SEARCH · Search NASA

Results for “Tokamak plasma”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2022

To dee or not to dee: costs and benefits of altering the triangularity of a steady-state DEMO-like reactor

Shaping a tokamak plasma to have a negative triangularity may allow operation in an ELM-free L-mode regime and with a larger strike-point radius, ameliorating divertor power-handling requirements. However, the shaping has a potential drawback in the form of a lower no-wall ideal beta limit, found using the MHD codes CHEASE and DCON. Using the new fusion systems code FAROES, we construct a steady-state DEMO2 reactor model. This model is essentially zero-dimensional and neglects variations in physical mechanisms like turbulence, confinement, and radiative power limits, which could have a substantial impact on the conclusions deduced herein. Keeping its shape otherwise constant, we alter the triangularity and compute the effects on the levelized cost of energy (LCOE). If the tokamak is limited to a fixed B field, then unless other means to increase performance (such as reduced turbulence, improved current drive efficiency or higher density operation) can be leveraged, a negative-triangularity reactor is strongly disfavored in the model due to lower \beta_N limits at negative triangularity, which leads to tripling of the LCOE. However, if the reactor is constrained by divertor heat fluxes and not by magnet engineering, then a negative-triangularity reactor with higher B0 could be favorable: we find a class of solutions at negative triangularity with lower peak heat flux and lower LCOE than those of the equivalent positive triangularity reactors.

CHEASE↗

DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Kinetic equilibrium reconstruction and the impact on stability analysis of KSTAR plasmas

High fidelity kinetic plasma equilibrium reconstructions are an essential requirement for accurate stability and disruption prediction analyses to support continuous operation of high beta tokamak plasmas. The present kinetic equilibrium reconstructions of plasmas in the KSTAR device include plasma density and temperature profiles from Thomson scattering and ion temperature from charge exchange spectroscopy diagnostics, and allowance for fast particle pressure. In addition, up to 25 channels of motional Stark effect diagnostic data are used to constrain the magnetic field pitch angle profile in the plasma to produce a reliable computation of the safety factor, q, profile. In this work, H-mode plasmas exhibit clear pedestal characteristics in the reconstructed pressure profile compared to internal transport barrier or L-mode plasmas. The plasma configuration and vertical position of inner strike points are validated by CCD and infrared camera images. Ideal and resistive magnetohydrodynamic (MHD) stability analyses using the DCON and resistive DCON codes utilize these kinetic equilibrium reconstructions to compare to experimental plasma stability. Equilibria with sufficiently low convergence error can provide reliable computation of ideal and resistive MHD stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

ARC disruption physics and strategy

Commonwealth Fusion Systems (CFS) plans to operate a tokamak power plant called ARC in the early 2030s. Tokamak plasmas have stability limits that, if crossed, lead to a rapid termination of the plasma, referred to as a disruption. Disruptions pose a melt risk to the first wall resulting from thermal and non-thermal particle heat fluxes, and an electromagnetic loading risk on all metal components within the equilibrium coils. A comprehensive set of models is used herein to provide an assessment of both mitigated and unmitigated ARC disruption loads. A preliminary massive gas injection system is baselined and a runaway electron mitigation coil option is proposed to close possible gaps in the baseline. It is predicted that all ARC disruption loads are within a factor of 2 of the disruption loads in SPARC, a tokamak presently under construction by CFS, and therefore SPARC provides an opportunity to calibrate models, test solutions and inform the design of ARC. The goal for ARC is disruption-free operation, however, the pragmatic design target is to withstand one mitigated disruption per day, and to restart the plasma following mitigation in tens of seconds without interrupting the power output. Unmitigated disruptions must be rare, and experience with unmitigated disruption impacts in SPARC will better define what rare means. The implications of this strategy for plasma disruptivity and disruption prediction are discussed, and operating the ARC scenario on SPARC is expected to refine the ARC final design and operational plan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

47 OTHER INSTRUMENTATION↗

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 OSTI2023

Species mix & distribution function dependence of energetic-ion driven magnetospheric instabilities on DIII-D (Final Report)

Energetic ions drive a rich variety of instabilities in the complex multi-species plasma of Earth’s magnetosphere. This experiment exploited the highly flexible, high power particle beam injection capability and sophisticated diagnostic suite of the DIII-D tokamak to advance understanding of two particular instabilities of significant interest to the magnetospheric research community that have closely-related counterparts in beam-heated tokamak plasmas: electromagnetic ion cyclotron waves (EMIC) [1], observed at frequencies somewhat below the ion cyclotron frequencies of bulk ion species, and equatorial noise [2], observed at energetic-ion cyclotron harmonics. In the magnetic fusion energy domain, such waves are called ICE (Ion Cyclotron Emission). Recent substantial advances in diagnostic capability in the ion cyclotron range of frequencies in DIII-D, including measurement of wave polarization, spa)al structure and amplitude created an exciting opportunity to challenge and drive advances in theoretical understanding of these phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Experimental Observation of Magnetic Island Heteroclinic Bifurcation in Tokamaks

We report empirical observations of magnetic island heteroclinic bifurcation for the first time. This behavior is observed in interacting coupled 2/1 tearing modes in the core of a DIII-D tokamak plasma. Poincaré maps constrained by measured magnetic amplitudes and phasing show bifurcation from heteroclinic to homoclinic topology in the 2/1 island as the 4/2 relative amplitude (R 4/2 ) decreases. Initially, the local electron temperature peak in the 2/1 island splits, consistent with two O points. As R 4/2 decreases a single peak forms, consistent with one O point. Lastly, these call for developing tearing stability theory and control solutions for heteroclinic islands in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Experiments and Simulations of the Dusty Properties of Hypervelocity Impact Plasmas

Hypervelocity impacts, referring to impacts occurring at several times the speed of sound in the target, result in inertial stresses greatly exceeding the material strength and produce behavior that is not fully understood. These impact phenomena are seen in space through impacts of orbital debris and meteoroids on spacecraft and airless natural bodies but can also be produced in laboratory environments using light-gas guns and electrostatic dust accelerators. When a hypervelocity impact occurs, plasma is generated both from thermal ionization and from pressure ionization, producing a wide range of plasma densities depending on the impactor’s velocity. This plasma is initially formed in the warm dense matter (WDM) regime and governed by high energy density (HED) hydrodynamics before it rapidly expands and dissipates, spanning many orders of magnitude in density and length scale. The plasma can also contain a dust component, which is particularly relevant for the lower-velocity range of hypervelocity impacts with parameters that are particularly relevant to Tokamak plasmas. Therefore, the properties and dynamics of these plasmas cannot be fully characterized without considering dusty plasma effects, which are ubiquitous but poorly understood.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Spatial heterodyne spectroscopy for fast local magnetic field measurements of magnetized fusion plasmas

A novel spectroscopy diagnostic for measuring internal magnetic fields in high temperature magnetized plasmas has been developed. It involves spectrally resolving the Balmer- α (656 nm) neutral beam radiation split by the motional Stark effect with a spatial heterodyne spectrometer (SHS). The unique combination of high optical throughput (3.7 mm 2 sr) and spectral resolution (δλ ~ 0.1 nm) allows these measurements to be made with time resolution $\ll$1 ms. The high throughput is effectively utilized by incorporating a novel geometric Doppler broadening compensation technique in the spectrometer. The technique significantly reduces the spectral resolution penalty inherent to using large area, high-throughput optics while still collecting the large photon flux provided by such optics. In this work, fluxes of order 10 10 s –1 support the measurement of deviations of <5 mT (Δ λ Stark ~ 10 –4 nm) in the local magnetic field with 50 µs time resolution. Example high time resolution measurements of the pedestal magnetic field throughout the ELM cycle of a DIII-D tokamak plasma are presented. Local magnetic field measurements give access to the dynamics of the edge current density, which is essential to understanding stability limits, edge localized mode generation and suppression, and predicting performance of H-mode tokamaks.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Role of electrostatic perturbation on kinetic resistive wall mode with application to spherical tokamak

Abstract A more complete non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid formulation is developed by including the perturbed electrostatic potential δφ in the particle Lagrangian. The fluid-like counter-parts of the hybrid equations, in the Chew-Goldberger-Low high-frequency limit, are also derived and utilized to test the new toroidal implementation in the MARS-K code. Application of the updated non-perturbative hybrid model for a high- β spherical tokamak plasma in MAST finds that the perturbed electrostatic potential generally plays a minor role in the n = 1 ( n is the toroidal mode number) resistive wall mode instability. The effect of δφ is largely destabilizing, with the growth rate of the instability increased by several (up to 20) percent as compared to the case without including δφ . A similar relative change is also obtained for the kinetic-induced resonant field amplification effect at high- β in the MAST plasma considered. The updated capability of the MARS-K code allows quantitative exploration of drift kinetic effects on various MHD instabilities and the antenna-driven plasma response where the electrostatic perturbation, coupled to magnetic perturbations, may play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Toroidal modeling of Alfvén eigenmodes excited by runaway electrons in DIII-D and ITER

The non-perturbative MHD-kinetic hybrid code MARS-K (Liu et al 2014 Phys. Plasmas 21 056105) is updated to include relativistic effects for kinetic fast particles, enabling the code to model excitation of Alfvén eigenmodes (AEs) by runaway electrons (REs) in post-disruption tokamak plasmas. Applying the updated code to RE beams in both DIII-D and ITER, a zoo of AE modes triggered by trapped REs due to precessional drift-kinetic resonances is computed while scanning the RE energy. At fixed RE energy, multiple unstable roots are also excited. These AE modes possess radially different eigenmode structures, ranging from global modes to core-localized ones. The computed mode frequency is in the Alfvén frequency range, increasing with the assumed RE energy in a staircase fashion and quantitatively matching the experimental measurement (in DIII-D). At the (more relevant) high-frequency range (above 1 MHz), the modeled eigenmodes are identified as compressional AEs (CAEs) in DIII-D and a mixture of CAE and shear Alfvén waves in ITER.

Alfvén eigenmodes↗

DOE OSTI2021

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Implementation and assessment of an extended hydrogenic molecular model in UEDGE

Abstract Two sets of pure deuterium plasmas are simulated using UEDGE, one including atoms and molecules and the other including atoms only as separate fluid species. Simulation results in one and two dimensions are reported to assess the role of molecules in tokamak plasmas. It is shown that thermal coupling of the molecules to the plasma can be relevant for the onset of detachment under highly collisional conditions in simple geometries, but play a modest role at low collisionality. Ion‐molecule equipartition presents an additional ion energy loss channel, dissipating ion energy that would otherwise heat the electrons by equipartition. The resulting reduced electron temperature increases the ion‐electron recombination by an order of magnitude as the plasma density is increased, yielding earlier and deeper detachment.

Holm, A.↗

DOE OSTI2020

Microturbulence-mediated route for energetic ion transport and Alfvénic mode amplitude oscillations in tokamaks

The new regimes of Alfvén eigenmode (AE) induced fast ion transport in tokamak plasmas are reported, in which microturbulence plays the role of a mediator of fast ion relaxation. Coulomb collisional scattering alone leads to small AE amplitudes and does not reproduce the steady state regimes observed in experiments. We show that in nonlinear regimes the sufficiently large effective pitch angle scattering due to microturbulence can lead to steady state AE amplitude evolution. This indicates a new route for fast ion losses, which is beyond the scenarios described in “Energetic ion transport by microturbulence is insignificant in tokamaks”. As a result, microturbulence can significantly increase the amplitude of AEs in predictive simulations of burning plasma experiments such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Ion cyclotron emission in Maxwellian plasmas

Destabilization of ion cyclotron waves—waves with frequencies close to ion cyclotron harmonics—in inhomogeneous plasmas with Maxwell velocity distribution is considered. A new mechanism of destabilization of these waves is found, whereas the known interaction between cyclotron waves and drift waves is shown to be hardly able to lead to instabilities in realistic tokamak plasmas. Our finding is that the resonance wave–particle interaction in the presence of temperature gradient can change the diamagnetic drift frequency in such a way that the destabilizing influence of the diamagnetic drift exceeds Landau damping. This occurs when the energy of resonance particles is sufficiently high, which is always the case due to the infinite tail of the Maxwell distribution. Particles with higher energies provide a larger ratio of drive to damping, but the number of resonance particles and, thus, instability growth rates are exponentially small when particle energy is very high. Therefore, only moderately suprathermal particles can lead to observable instabilities. A condition for instabilities driven by these particles is obtained. Destabilization of electrostatic cyclotron waves and ordinary cyclotron waves is studied.

Kolesnichenko, Ya. I. (ORCID:0000000323219681)↗

DOE OSTI2022

Predicting resistive wall mode stability in NSTX through balanced random forests and counterfactual explanations

Abstract Recent progress in the disruption event characterization and forecasting framework has shown that machine learning guided by physics theory can be easily implemented as a supporting tool for fast computations of ideal stability properties of spherical tokamak plasmas. In order to extend that idea, a customized random forest (RF) classifier that takes into account imbalances in the training data is hereby employed to predict resistive wall mode (RWM) stability for a set of high beta discharges from the NSTX spherical tokamak. More specifically, with this approach each tree in the forest is trained on samples that are balanced via a user-defined over/under-sampler. The proposed approach outperforms classical cost-sensitive methods for the problem at hand, in particular when used in conjunction with a random under-sampler, while also resulting in a threefold reduction in the training time. In order to further understand the model’s decisions, a diverse set of counterfactual explanations based on determinantal point processes (DPP) is generated and evaluated. Via the use of DPP, the underlying RF model infers that the presence of hypothetical magnetohydrodynamic activity would have prevented the RWM from concurrently going unstable, which is a counterfactual that is indeed expected by prior physics knowledge. Given that this result emerges from the data-driven RF classifier and the use of counterfactuals without hand-crafted embedding of prior physics intuition, it motivates the usage of counterfactuals to simulate real-time control by generating the β N levels that would have kept the RWM stable for a set of unstable discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗