Search NASA⌕ Search

SEARCH · Search NASA

Results for “Toroidal plasma confinement”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An extended variational method for the resistive wall mode in toroidal plasma confinement devices

The external-kink stability of a toroidal plasma surrounded by a rigid resistive wall is investigated. The well-known analysis of Haney and Freidberg is rigorously extended to allow for a wall that is sufficiently thick that the thin-shell approximation does not necessarily hold. A generalized Haney–Freidberg formula for the growth-rate of the resistive wall mode is obtained. Thick-wall effects do not change the marginal stability point of the mode but introduce an interesting asymmetry between growing and decaying modes. Growing modes have growth-rates that exceed those predicted by the original Haney–Freidberg formula. On the other hand, decaying modes have decay-rates that are less than those predicted by the original formula. The well-known Hu–Betti formula for the rotational stabilization of the resistive wall mode is also generalized to take thick-wall effects into account. Increasing wall thickness facilitates the rotational stabilization of the mode, because it decreases the critical toroidal electromagnetic torque that the wall must exert on the plasma. On the other hand, the real frequency of the mode at the marginal stability point increases with increasing wall thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Edge state selection by modulating E × B shearing profile in toroidally confined plasmas

We investigate how the E × B shearing profile impacts the energy transport induced by the MHD mode in the edge of H-mode plasmas. By enhancing the outer layer (ψ nor ≈ 0.9–1) shear, the growth rate spectrum γ(n) gets narrowed with the fastest growth rate having a lower toroidal mode number. In the nonlinear stage, the MHD fluctuations tend to reach a coherent, oscillating state, which is analogous to a quiescent (Q) H mode state with edge-harmonic-oscillations. By enhancing the inner layer (ψ nor ≈ 0.7–0.8) shear, γ(n) shifts to lower n, and its width does not change much. In the nonlinear stage, the edge fluctuations tend to stay in an incoherent, oscillating state, which is analogous to a turbulent QH-mode state. Here, we also explore the underlying physics of the results and show that they are closely related to the different cross-phase dynamics of the edge MHD fluctuations induced by the variation of the edge E × B shearing profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
DOE OSTI2020

Curvature of Radial Electric Field Aggravates Edge Magnetohydrodynamics Mode in Toroidally Confined Plasmas

We show that the radial electric field (E r ) plays a dual role in edge magnetohydrodynamics (MHD) activity. While E r shear (first spatial derivative of E r ) dephases radial velocity and displacement, and so is stabilizing, a new finding here is that E r curvature (second spatial derivative of E r ) tends to synchronize the radial velocity and displacement, and so destabilizes MHD. As a highlighted result, we analytically demonstrate that E r curvature can destabilize an otherwise stable kink mode, and so form a joint vortexkink mode. The synergetic effects of E r shear and E r curvature in edge MHD extend the familiar E × B shearing paradigm. This theory thus explains the experimental findings that a deeper E × B well may aggravate edge MHD, and so trigger the formation of the edge harmonic oscillation. Here, a simple criterion linking E r structure and the edge MHD activity is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Gravitational steady states of solar coronal loops

Coronal loops on the surface of the sun appear to consist of curved, plasma-confining magnetic flux tubes or “ropes,” anchored at both ends in the photosphere. Toroidal loops carrying current are inherently unstable to expansion in the major radius due to toroidal-curvature-induced imbalances in the magnetic and plasma pressures. An ideal MHD analysis of a simple isolated loop with density and pressure higher than the surrounding corona, based on the theory of magnetically confined toroidal plasmas, shows that the radial force balance depends on the loop internal structure and varies over parameter space. It provides a unified picture of simple loop steady states in terms of the plasma beta β o , the inverse aspect ratio ϵ = a/R o , and the MHD gravitational parameter $\hat{G}$ ≡ ga/v$2\atop{A}$, all at the top of the loop, where g is the acceleration due to gravity, a the average minor radius, and v A the shear Alfvén velocity. In the high and low beta tokamak orderings, β o = 2n o T/(B$2\atop{0}$/2μ o ) ~ ϵ 1 and ϵ 2 , that fit many loops, the solar gravity can sustain nonaxisymmetric steady states at $\hat{G}$ ~ ϵβ o that represent the maximum stable height. At smaller \hat{G} ≤ ϵ 2 β o , the loop is axisymmetric to leading order and stabilized primarily by the two fixed loop ends. Very low beta, nearly force-free, steady states with β o ~ ϵ 3 may also exist, with or without gravity, depending on higher order effects. The thin coronal loops commonly observed in solar active regions have ϵ ≃ 0.02 and fit the high beta steady states. $\hat{G}$ increases with loop height. Fatter loops in active regions that form along magnetic neutral lines and may lead to solar flares and Coronal Mass Ejections have ϵ ≃ 0.1–0.2 and may fit the low beta ordering. Larger loops tend to have $\hat{G}$ > ϵβ o and be unstable to radial expansion because the exponential hydrostatic reduction in the density at the loop-top reduces the gravitational force -ρ$\hat{G}\hat{R}$ below the level that balances expansion, also in agreement with the observation that most sufficiently large loops grow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Full-wave simulations of ICRF heating regimes in toroidal plasma with non-Maxwellian distribution functions

At the power levels required for significant heating and current drive in magnetically-confined toroidal plasma, modification of the particle distribution function from a Maxwellian shape is likely (Stix 1975 Nucl. Fusion 15 737), with consequent changes in wave propagation and in the location and amount of absorption. In order to study these effects computationally, both the finite-Larmor-radius and the high-harmonic fast wave (HHFW), versions of the full-wave, hot-plasma toroidal simulation code TORIC (Brambilla 1999 Plasma Phys. Control. Fusion 41 1 and Brambilla 2002 Plasma Phys. Control. Fusion 44 2423), have been extended to allow the prescription of arbitrary velocity distributions of the form f(v(parallel to), v(perpendicular to) , psi, theta). For hydrogen (H) minority heating of a deuterium (D) plasma with anisotropic Maxwellian H distributions, the fractional H absorption varies significantly with changes in parallel temperature but is essentially independent of perpendicular temperature. On the other hand, for HHFW regime with anisotropic Maxwellian fast ion distribution, the fractional beam ion absorption varies mainly with changes in the perpendicular temperature. The evaluation of the wave-field and power absorption, through the full wave solver, with the ion distribution function provided by either a Monte-Carlo particle and Fokker-Planck codes is also examined for Alcator C-Mod and NSTX plasmas. Non-Maxwellian effects generally tend to increase the absorption with respect to the equivalent Maxwellian distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Full-wave simulations of ICRF heating regimes in toroidal plasmas with non-Maxwellian distribution functions

At the power levels required for significant heating and current drive in magnetically-confined toroidal plasma, modification of the particle distribution function from a Maxwellian shape is likely [T. H. Stix, Nucl. Fusion, 15 737 (1975)], with consequent changes in wave propagation and in the location and amount of absorption. In order to study these effects computationally, both the finite-Larmor-radius and the high-harmonic fast wave (HHFW), versions of the full-wave, hot-plasma toroidal simulation code TORIC [M. Brambilla, Plasma Phys. Control. Fusion 41, 1 (1999) and M. Brambilla, Plasma Phys. Control. Fusion 44, 2423 (2002)], have been extended to allow the prescription of arbitrary velocity distributions of the form f(v||, v_perp, psi , theta). For hydrogen (H) minority heating of a deuterium (D) plasma with anisotropic Maxwellian H distributions, the fractional H absorption varies significantly with changes in parallel temperature but is essentially independent of perpendicular temperature. On the other hand, for HHFW regime with anisotropic Maxwellian fast ion distribution, the fractional beam ion absorption varies mainly with changes in the perpendicular temperature. The evaluation of the wave-field and power absorption, through the full wave solver, with the ion distribution function provided by either aMonte-Carlo particle and Fokker-Planck codes is also examined for Alcator C-Mod and NSTX plasmas. Non-Maxwellian effects generally tends to increase the absorption with respect to the equivalent Maxwellian distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Noninductively Driven Tokamak Plasmas at Near-Unity Toroidal Beta

Access to and characterization of sustained, toroidally confined plasmas with a very high plasma-to-magnetic pressure ratio (β t ), low internal inductance, high elongation, and nonsolenoidal current drive is a central goal of present tokamak plasma research. Stable access to this desirable parameter space is demonstrated in plasmas with ultralow aspect ratio and high elongation. Local helicity injection provides nonsolenoidal sustainment, low internal inductance, and ion heating. Equilibrium analyses indicate β t up to ~100% with a minimum |B| well spanning up to ~50% of the plasma volume.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

What Can and Cannot Be Expected from Tokamak Fusion

This article is a response to the analysis of nuclear fusion by V. V. Orlov and L. I. Ponomarev in their article “Nuclear problems of thermonuclear power generation” [Atomnaya Energiya, 124, No. 2, 105–110 (2018); Atomic Energy, 124, No. 2, 129–138 (2018)], pointing to the decades-long lack-of-progress in thermonuclear energy research. Here, the authors of the article attribute the lack-of-progress to fundamental problems in the nuclear aspects of thermonuclear fusion. Here, the reason for the current stagnation – essentially, the degeneration of thermonuclear fusion that commenced in the 1950s–1970s with rapid development – is presented. The scientific investigations have been and remain within the scope of the Tamm–Sakharov concept of toroidal plasma confinement, which provided the first impetus to progress but soon showed its insufficiency. An additional element is considered – the suppression of plasma cooling by using moving liquid lithium to pump out plasma. For tokamaks this opens up new prospects for obtaining burning plasma suitable for use in hybrid (fission-fusion) reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Transport dynamics of self-consistent, near-marginal drift-wave turbulence. I. Investigation of the ability of external flows to tune the non-diffusive dynamics

The reduction of turbulent transport across sheared flow regions has been known for a long time in magnetically confined toroidal plasmas. However, details of the dynamics are still unclear, in particular, in what refers to the changes caused by the flow on the nature of radial transport itself. In Paper II, we have shown in a simplified model of drift wave turbulence that, when the background profile is allowed to evolve self-consistently with fluctuations, a variety of transport regimes ranging from superdiffusive to subdiffusive open up depending on the properties of the underlying turbulence. In this study, we show that externally applied sheared flows can, under the proper conditions, cause the transport dynamics to be diffusive or subdiffusive

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗
DOE OSTI2025

Efficient analysis of magnetic field line behavior in toroidal plasmas

The confinement of plasmas in tokamaks and stellarators depends on magnetic field lines lying in nested toroidal surfaces. The transition near the plasma edge away from the lines lying in magnetic surfaces defines properties of divertors. The transition in time defines properties of disruptions. Divertor design and disruption analyses require a detailed understanding of these transitions. The use of a Fourier transform coupled with a Gaussian window function allows far more information to be extracted about these transitions using far shorter field line integrations than can be obtained using traditional methods based on Poincaré plots. The physics of divertors and disruptions is reviewed to clarify why the type of information that can be gained from more efficient methods of analysis is of central importance to the fusion program based on magnetic confinement.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2024

Transport from electron-scale turbulence in toroidal magnetic confinement devices

Plasma transport driven by turbulence ultimately determines the energy confinement performance of controlled fusion devices regardless of their confinement schemes and configurations. A large variety of plasma instabilities have been proposed for driving turbulence responsible for anomalous plasma transport beyond classical/neoclassical transport due to collisions. Although ion-scale turbulence usually dominates due to its large eddy size and saturation level, electron-scale turbulence has been recognized to be important in regions where ion-scale turbulence is suppressed (e.g., in internal transport barrier and in spherical tokamak H-mode plasmas) or is close to marginality. Electron-scale turbulence has been shown to nonlinearly interact with ion-scale turbulence, which modifies the dynamics of both and affects the resulting plasma transport, particularly when ion-scale instability is weakly driven. In this review paper, we focus on electron-scale turbulence that is believed to operate in magnetically confinement fusion devices and aim to provide a review of theoretical, numerical, and experimental developments in understanding electron-scale turbulence and its role in driving anomalous plasma turbulence. In particular, we focus on the electrostatic electron temperature gradient (ETG) mode which is the most widely recognized plasma instability underlying electron-scale turbulence observed in magnetically confined plasmas. Finally, we note that there are other less studied instabilities that might be responsible for observed electron-scale turbulence, most notably ubiquitous mode, and short-wavelength ion temperature gradient (SWITG) mode, which will be briefly touched on in this review.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Is it possible that MHD instability triggers a transition into the improved confinement regime of toroidal plasmas?

Analyzing the experimental data from L-2M stellarator and performing appropriate theoretical analysis, the authors affirmatively answer the question raised in the title of the article. The external peeling mode is proposed as the instability triggering the transition. Peeling mode stability in magnetic-hill Mercierstable plasmas is analyzed analytically. It is shown that correctly taking into account vacuum region forbids internal peeling modes, while the external peeling mode has a threshold with respect to the plasma pressure gradient. Calculations and experimental data are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗