Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma shaping”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Generation of X-rays and neutrons with a RF-discharge

An experimental study concerning disk shaped plasma structures was performed. Such disk-shaped structures can be obtained using an rf discharge in hydrogen. The applied frequency was 1-2 Mhz. In case of operation in deuterium it was found that the discharge emits neutrons and X-rays, although the applied voltage is only 2 kV. This phenomenon was explained by assuming formation of plasma cavitons which are surrounded by high electric fields. The condition for formation of these cavitons is that the applied rf frequency is equal to the plasma frequency. The ions trapped in these resonance structures acquire sufficient energy that they can undergo fusion reactions with the ions in the surrounding gas.

Schneider, R. T.↗
DOE OSTI2022

Design of a Novel Variable Geometry Divertor for Tokamaks

The divertor is a key component of fusion reactors, allowing exhaust of gas, impurities, and helium ash to preserve plasma purity. The divertor geometry strongly affects plasma performance, and it is designed to be compatible with different plasma shapes in present-day fusion experiments. Here, we present a novel concept for a variable geometry divertor, in which the divertor baffle tiles are reorientable by external actuation. Implementation of this concept in a medium-sized research tokamak would uniquely provide the flexibility to tailor divertor geometry to the plasma configuration and also enable study of the effect of divertor closure on plasma performance. To ensure compatibility with typical tokamak operations, the adjustable divertor must withstand the effects of significant mechanical and thermal stresses such as MW/m 2 -scale heat fluxes and large electromagnetic fields, e.g., disruption forces. The technological solutions for actuation mechanisms, cooling system, gas baffling and plasma-facing components are assessed. A functional reduced-scale model with movable outer divertor target baffle tiles is developed and the actuation mechanism is tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental characterization of turbulence properties in negative triangularity DIII-D plasmas

Negative Triangularity (NT) plasmas have demonstrated robustly ELM-free high-performance operation and provide a unique testbed to study how plasma shaping affects turbulent transport. The Beam Emission Spectroscopy (BES) diagnostic provides localized 2D measurements of low-k ($k_{\theta} \rho_s < 1$) density fluctuations. In a sweep of upper triangularity at fixed power, H-mode access is suppressed and an NT-edge is observed. The turbulence amplitude ($\tilde{n}/n$) is shown to decrease by $\sim50\%$ for $\rho<0.9$ in the NT-edge phase as compared to the H-mode phase. Additionally, low-velocity edge modes below 70 kHz are suppressed by triangularity and the dominant mode propagating in the electron diamagnetic direction is seen to broaden in wavenumber space in the NT-edge phase. In a strong NT plasma ($\delta_{avg}\sim-0.5$), low amplitude modes ($\tilde{n}/n <0.5\%$) propagating in the ion-diamagnetic direction with radially-poloidally symmetric eddy structure are identified for $\rho \sim 0.65-0.83$ consistent with Ion Temperature Gradient (ITG) turbulence. Modes consistent with Trapped Electron Mode (TEM) turbulence are observed propagating in the electron-diamagnetic direction for $\rho\gtrsim0.83$ with poloidally extended eddy structure and reduced amplitudes observed at the separatrix. The turbulence properties presented in this paper help validate our understanding of NT turbulence and help explain the improved confinement and unique edge features of NT plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of elongation and triangularity on plasma response to resonant magnetic perturbations

Abstract Systematic toroidal modeling of the plasma response to the n = 1–4 ( n is the toroidal mode number) resonant magnetic perturbation (RMP) field is carried out in order to understand the plasma-shaping effect on controlling the type-I edge-localized modes (ELMs) in tokamak experiments. Considered are large variations of the plasma elongation and triangularity at a fixed edge safety factor q a , for limiter plasmas with both single-null (SN) and double-null (DN) divertor-like boundary shapes. Numerical results assuming conformal 3D RMP coils show that (i) the optimum coil phasing between the upper and lower rows for ELM control becomes increasingly sensitive to the plasma elongation with higher- n toroidal spectra, (ii) the optimum coil phasing is however essentially independent of the plasma triangularity for all n = 1–4 RMP fields, (iii) with the same coil current and the optimum coil phasing, high elongation generally favors ELM control but it may be more challenging for plasmas with intermediate elongation and with lower- n ( n = 1–2) RMPs, and (iv) higher triangularity is generally always better for ELM control with all n = 1–4 fields for both DN and SN divertor-like plasma boundary shapes.

Physics↗
DOE OSTI2024

Inverse aspect-ratio expanded tokamak equilibria

Following Greene et al. [Phys. Fluids 14, 671 (1971)] and Connor et al. [Phys. Plasmas 31, 577 (1988); Plasma Phys. Control. Fusion 34, 161 (1992); and Nucl. Fusion 33, 1533 (1993)], the Grad-Shafranov equation for an axisymmetric tokamak plasma equilibrium is solved via an expansion in the, supposedly small, inverse aspect-ratio of the plasma, ϵ. The displacements of equilibrium magnetic flux-surfaces due to plasma shaping are assumed to be $\mathcal{O}$(ϵ) smaller than the minor radii of the surfaces, but no other restriction is placed on the nature of the shaping. The solution of the Grad-Shafranov equation is matched to a vacuum solution that extends to infinity, and consists of an expansion in toroidal functions. The external poloidal magnetic field generated by a finite set of discrete external poloidal magnetic field-coils is calculated, and incorporated into the toroidal function expansion. In this manner, the shape of a large aspect-ratio tokamak plasma is directly related to the currents flowing in the external poloidal field-coils. Finally, a pedestal in the plasma pressure, and the associated spike in the bootstrap current, are incorporated into the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of plasma squareness on RMP response in MAST-U

The effect of the squareness of the plasma boundary shape on the plasma response to the applied n = 2 (n is the toroidal mode number) resonant magnetic perturbation, and consequently on controlling the edge-localized mode (ELM) in MAST-U in the double-null configuration, is numerically systematically investigated. The modeling is performed utilizing three different plasma response models: the single-fluid model, the two-fluid model and the magnetohydrodynamic-kinetic hybrid model. The effect of the plasma squareness on ELM control is found to depend on the equilbrium kinetic profiles. Fixing the equilibrium radial profiles to that from a recent MAST-U discharge, the experimental plasma shape is found to be close to optimal for ELM control in terms of squareness. The optimal coil phasing does not significantly vary with varying plasma squareness. These findings robustly hold with (i) two figures of merit (FoMs) defined for predicting the ELM control feasibility, and (ii) three plasma response models employed for evaluating these FoMs.

MAST-U↗
DOE OSTI2022

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional modeling of a plasma in a strong azimuthal magnetic field

Three-dimensional magnetohydrodynamic (MHD) simulations are able to model the generation of disc-shaped plasma, driven by laser ablation from a current carrying rod in a pulsed-power machine producing azimuthal magnetic fields of 2–3 MG. The plasma at such extreme conditions is unique in that the parameter space for the plasma β and Hall parameter χ transition from below unity to greater than unity at different stages of the plasma generation. In simulations, the formation of the plasma disc in the azimuthal direction, is driven by heat flux from the laser spot and depends on the set of transport coefficients used in simulations. The most recent set of transport coefficients leads to the formation of plasma ejecta at the back end of the rod, that qualitatively matches experiments. Specifically, the cross-gradient Nernst effect, which twists magnetic field, is shown to have a large effect on the shape of the back end ejecta. In the direction along the axis of the rod there is propagation of perturbations from the disc as observed in experiments. In simulations the period of temperature perturbations is in good agreement with experimental results. Furthermore, an instability due to coupling of heat flux and the magnetic field advection provides a possible explanation for perturbation growth along the axis of the rod, and the instability growth rate is consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of anisotropic thermal transport on tearing mode stability in negative versus positive triangularity plasmas

The combined effects of anisotropic thermal transport and the plasma shaping, including negative triangularity, on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability are numerically investigated utilizing the MARS-F code [Liu et al., Phys. Plasmas 7, 3681–3690 (2000)]. While varying the plasma boundary triangularity, the TM stability is found to be dictated by the competing effects of the Shafranov shift induced stabilization and the bad-curvature induced destabilization. The negative triangularity shape increases the Shafranov shift (stabilizing) in the plasma core but also enlarges bad-curvature regions (destabilizing) near the plasma edge, with the net effect being largely destabilizing for the TM as compared to the positive triangularity counter-part. Large negative triangularity however can also lead to more stabilization for the plasma core-localized TM. Anisotropic thermal transport reduces the stabilizing effect on the TM associated with the favorable averaged curvature, resulting in more unstable core-localized TMs in both negative and positive triangularity plasmas. But the opposite effect can also take place for the edge-localized TM in finite-pressure plasmas with negative triangularity.

Physics↗
DOE OSTI2024

Real-time plasma equilibrium reconstruction and shape control for the MAST Upgrade tokamak

Real-time magnetic control has been developed to deliver precise control of multiple plasma shape parameters for advanced divertor configurations, including double-null, Super-X, X-point target and X-divertor for the first time on the MAST Upgrade (MAST-U) spherical tokamak. Successful real-time magnetic equilibrium control of different plasma shape variables has been accomplished in the 2022–2023 MAST-U experimental campaign for the advanced MAST-U divertor configurations. Application of the MAST-U boundary reconstruction algorithm, LEMUR, is described and compared with off-line equilibrium reconstruction and diagnostic measurements. The process of design and verification of the axisymmetric magnetic control schemes using a suite of control analysis tools (known collectively as TokSys) is also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predictive Integrated Modeling of Heated Tokamak Plasmas (Final Technical Report)

The research carried out by the Lehigh University fusion group during this reporting period: (1) Self-consistent predictive core-pedestal ITER scenario modeling; (2) Microtearing instabilities and electron thermal transport in standard and spherical tokamak discharges; (3) Neural network model of the Multi-Mode anomalous transport module; (4) Anomalous transport in KSTAR, JET, EAST, Alcator C-Mod and DIII-D discharges; (5) Theory for transport in magnetized plasmas; (6) Control-oriented modeling of the q-profile dynamic response in preparation for advanced scenario control in KSTAR and in EAST; (7) Modeling of edge harmonic oscillation in DIII-D QH mode discharges; (8) Elucidating plasma dynamics in Hasegawa-Wakatani turbulence by information geometry; (9) Non-linear Langevin and Fractional Fokker–Planck equations for anomalous diffusion by Levy stable processes; and (10) Investigation of the plasma shaping effects on the H-mode pedestal structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic-ballooning-bifurcation in tokamak pedestals across shaping and aspect-ratio

We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of anisotropic fast ions on internal kink stability in DIII-D negative and positive triangularity plasmas

Recent DIII-D experiments show that sawtooth stability is strongly affected by anisotropic fast ions from neutral beam injection (NBI) in both negative and positive triangularity plasmas. Fast ions from co-current NBI are stabilizing for the sawtooth stability, resulting in longer sawtooth periods. On the other hand, fast ions from counter-current NBI are destabilizing, leading to small and frequent sawteeth. The relative change of sawtooth period and amplitude is more than a factor of two. These observations appear to hold in both plasma shapes. Non-perturbative toroidal modeling, utilizing the magnetohydrodynamic-kinetic hybrid stability code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503), reveals an asymmetric dependence of the stability of the n = 1 (n is the toroidal mode number) internal kink mode on the injection direction of NBI, being qualitatively consistent with the experimentally observed sawtooth behavior. The MARS-K modeling results suggest that anisotropic fast ions affect the mode growth rate and frequency through both adiabatic and non-adiabatic contributions. Here, the asymmetry of the internal kink mode instability relative to the NBI direction is mainly due to the non-adiabatic contribution of passing fast ions, which stabilize (destabilize) the internal kink with the co-(counter-) current NBI as compared to the fluid counterpart. However, finite orbit width (FOW) correction to passing particles partially cancels the asymmetry. Trapped particles are always stabilizing due to precessional drift resonance. Modeling also shows that fast ions affect the internal kink in a similar manner in both negative and positive triangularity plasmas, although being slightly more unstable in the negative triangularity configuration already in the fluid limit. The similarity is mainly attributed to the fact that the mode is localized in the plasma core region, with very similar eigenmode structures in both negative and positive configurations. Furthermore, MARS-K modeling indicates that other factors, such as the plasma rotation and the drift kinetic effects of thermal plasmas, weakly modify the mode stability as compared to the drift kinetic resonance effects and FOW correction of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗