Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Toroidal Field”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Characterisation of ion temperature and toroidal rotation on the ST40 tokamak

ST40 is a high field low-aspect ratio spherical tokamak built and operated by Tokamak Energy Ltd. Recent plasma operations were aimed at exploring operational scenarios to maximise the central plasma temperature and have culminated in the achievement of thermal ion temperatures of over 9 keV. Here this paper presents ion temperature and toroidal rotation measurements performed on ST40 during the 2021–22 campaign for a range of different scenarios. Several independent diagnostic systems are used, analysing their correlation and interpreting their differences using new diagnostic forward models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

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 OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2020

Role of 3D neoclassical particle flux in density pump-out during ELM control by RMP in DIII-D

Application of three-dimensional (3-D) resonant magnetic perturbations, for the purpose of controlling edge localized modes (ELMs) in high confinement (H-mode) tokamak fusion operations, routinely leads to significantly reduction of the plasma density during the discharge, a phenomenon termed 'density pump-out'. Understanding the density pump-out phenomenon due to 3-D fields is a longstanding challenging issue. This work discusses a self-consistent toroidal model that allows initial value simulation of the density pump-out effect. The model successfully reproduces pump-out phenomenology in DIII-D ELM suppression experiments. The radial particle flux, associated with the 3-D neoclassical effects (the neoclassical toroidal viscosity, NTV), is discovered to play an important role in producing the observed large density pump- out. The key physics is associated with the drift kinetic Landau resonances between the toroidal precession of trapped thermal ions and electrons from one side, and the magnetic field perturbation from the other side. Such resonances, occurring near the top of the pedestal region, substantially enhances the radial particle flux associated with NTV. As a general consequence of the resonant NTV flux, larger density pump- out occurs in plasmas with lower collisionality and slower toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Design and simulation of multilayer coatings for a multi-channel Wolter-like x-ray imager with large field of view and high resolution

X-ray diagnostics are key instruments for understanding the physics behind inertial confinement fusion experiments. We report on the multilayer design optimization for the Toroidal X-ray Imager (TXI), a hard x-rays microscope instrument designed by Commissariat à l’énergie atomique (CEA) and Laboratoire Charles Fabry (LCF) to be installed on the National Ignition Facility. TXI includes six channels designed for three different energy bands centered on 8.7, 13, and 17.5 keV. Each channel is made up of two toroidal mirrors arranged in a Wolter-like configuration. The required field of view is 800 × 400 µm2, and the resolution should be better than 5 µm. In addition, we seek to estimate the spatial distribution of the temperature, which requires no spectral overlap of the different energy bands and a good spectral homogeneity of the image produced. The development of the multilayer coatings was performed in a two-step method. First, the coatings were optimized to obtain proper energy bands. Then, an x-ray tracing code was used to calculate the integrated optical response of each channel and adjust the response of the mirror to fulfill the requirements. To fulfill all the specifications, we propose an original design using a combination of two aperiodic coatings, one with a narrow bandwidth and the other one with a larger bandwidth.

hard x-rays↗
DOE OSTI2026

Experimental study of ECH pre-ionization on J-TEXT

An experimental study on electron cyclotron heating (ECH) pre-ionization has been conducted on J-TEXT in support of the joint experiment research for ITER plasma initiation. In this experiment, ECH power was injected to the vessel before the application of loop voltage ionize the neutral gas and form the initial plasma or so-called pre-plasma. The impact of several significant factors, such as magnetic field configuration, pre-fill gas pressure, ECH toroidal injection angle and ECH power on the evolution of pre-plasma are systematically studied, aiming to identify shared features, clarify their potential relationship and optimize the discharge parameters to generate a rather high pre-plasma density. By separating ECH power from the inductive start-up, the effect of pre-plasma on tokamak start-up can be observed. A dynamic magnetic configuration facilitates the transition of pre-plasma to tokamak plasma. To assess the influence of pre-plasma density on tokamak start-up, two kinds of magnetic field configurations are examined. While the effect of different pre-plasma densities on tokamak start-up is negligible, a significant difference is observed between pure ohmic start-up and start-up with pre-ionization. The studies presented here show evolutionary trends and threshold values needed to optimize ECH pre-ionization and also a feasible way to improve pre-plasma density and a configuration to stabilize pre-plasma for transition. Eventually, these results may contribute to multi-machine research and physics analysis that can assist ITER with its optimal preparation for first plasma operation.

ECH↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2022

Momentum transfer and plasma rotation caused by destabilized eigenmodes in tokamaks

In this work, the influence of magnetohydrodynamic eigenmodes destabilized by energetic ions on the momentum of these ions and concomitant sheared plasma rotation are studied. Two mechanisms affecting rotation are revealed: (i) spatial channelling (SC) – radially separated emission and absorption of the momentum; (ii) mode induced redistribution (MIR) across the magnetic field of the momentum of energetic ions. Forces arising during SC and MIR produce both toroidal and poloidal rotations. In addition, the momentum emission during SC leads to a radial flux of fast ions and generation of a radial electric field. Using the developed theory, estimates were made for the ITER (International Thermonuclear Experimental Reactor) 15 MA baseline scenario. They show that a global toroidicity-induced Alfvén eigenmode destabilized by alpha particles and neutral beam injection can result in significant radial electric field and forces applied to plasma. However, available data are not sufficient for a reliable prediction of the effects of SC and MIR in ITER. In general, one can expect that sheared rotation arising after destabilization of Alfvén modes and fast magnetoacoustic modes by energetic ions will tend to suppress the turbulence and improve plasma performance. The importance of plasma rotation is supported, in particular, by the fact that during the JET DTE1 experimental campaign the best parameters were achieved in a deuterium–tritium discharge where the rotation frequency was largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Islands and current singularities in quasisymmetric toroidal plasmas

The presence of current singularities in a quasisymmetric magnetic field is explored. Quasisymmetry is shown effective in isolating Pfirsch–Schlüter singularities, to leading order, to a single magnetic surface resonant with the helicity of the symmetry. The effects of pressure driven currents are analyzed, indicating that exclusion of this surface from the plasma volume reduces the potential opening of islands but does not generally eliminate them completely due to higher order asymmetric geometric effects. These three-dimensional effects are contained in quasisymmetry and indicate the complexity of finding consistent solutions and their potential sensitivity. The δ-function current singularities show a distinct quasisymmetric behavior only when the higher-order Fourier content of B is relevant for the resistive stability parameter D R (not included in leading-order near-axis expansions). In such scenarios, quasisymmetry proves advantageous, both in simplicity and avoidance of amplification by low-order rational surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Charge exchange recombination spectroscopy measurements of DIII-D poloidal rotation with poloidal asymmetry in angular rotation

Sixteen new tangential views for the charge exchange recombination (CER) spectroscopy diagnostic at DIII-D were installed in 2019 on the high-field side (HFS) of the tokamak with the main goal being the measurement of main-ion (deuterium) poloidal rotation. Eight of the new views are connected to spectrometers, which view the main-ion spectrum, adding main-ion measurements where there were previously none, and another eight new views increased the spatial resolution of existing impurity (carbon) measurements on the HFS. When combined with the existing low-field side measurements, measurements at two locations on flux surfaces out to a normalized minor radius of ≈0.6 are possible. The new tangential views have been used to measure the deuterium poloidal rotation directly for the first time using the Poloidal Asymmetry in Angular Rotation (PAAR) method. Furthermore, these new measurements enable further testing of the validity of neoclassical poloidal rotation predictions. Separate measurements of the radial electric field can be made for an impurity ion and the main-ion by combining the PAAR measurements with additional CER measurements of toroidal rotation, temperature, and density. These independent measurements of the radial electric field agree reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

$\mathrm{T}$ransient $\mathrm{CHI}$ System Design Studies for $\mathrm{P}$ EGASUS -$\mathrm{III}$

We report transient coaxial helicity injection (transient CHI), first developed on the Helicity Injected Torus-II (HIT-II) and later on the National Spherical Torus Experiment (NSTX) for implementing solenoid-free plasma current startup capability in a spherical tokamak (ST), is now planned to be tested on the Pegasus-III ST using a novel double-biased configuration. Such a configuration is likely needed for transient CHI deployment in a reactor. The transient CHI system optimization will be studied on Pegasus-III to enable startup toroidal persisting currents at the limits permitted by the external poloidal field coils. A transient CHI discharge is generated by driving injector current along magnetic field lines that connect the inner and outer divertor plates on one end of the ST. Simulations using the Tokamak Simulation Code are used to assess the transient CHI toroidal current generation potential and electrode gap location on the Pegasus-III. While past transient CHI systems have used high-voltage, oil-filled capacitors for driving the injector current, for improved safety, Pegasus-III will use a high-current capacitor bank based on low-voltage electrolytic capacitors. The designed and fabricated system is capable of over 32 kA. The modular design features permit the system to be upgraded to higher currents, as needed, to meet the future needs of the Pegasus-III facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗