Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Toroid”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Observations of a plectonemic configuration in a stable magnetized plasma jet

Astrophysical jets are collimated high-speed outflows emerging from spinning and accreting matter around celestial objects and may spontaneously result from self-organized processes. Magnetic self-organization is commonly observed in laboratory plasma physics experiments; however, they require close-fitting flux conservers to constrain and stabilize the toroidal or cylindrical structures. Here we report the first observations of a long, stable, free-boundary plasma jet far from chamber walls, embedding a double-helix magnetic structure resembling a force-free plectonemic Taylor state. The jets arise from an experimental setup that mimics an accretion disk and has no close-fitting solid flux conserver. The results support the hypothesis that self-organization could be a universal, intrinsic explanation for jet formation, collimation, and stability and may help explain double-helix features in celestial observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-speed feedback control of an oscillating magnetic helicity injector using a graphics processing unit

A real-time control system has been developed to control the amplitude, phase, and offset of bulk plasma parameters inside an oscillating magnetic helicity injector. Control software running entirely on an Nvidia Tesla P40 graphical processing unit is able to receive digitizer inputs and send response patterns to a Pulse Width Modulation (PWM) controller with a minimum control loop period of 12.8μs. With an input digitization rate of 10 MS/s, a three-parameter proportional integral differential controller is shown to be sufficient to inform the PWM controller to drive the desired oscillating plasma waveform with a frequency of 16.6 kHz that is located near the resonance of a coupled RLC circuit. In particular, the temporal phase of the injector waveform is held within 10 degrees of the target value. Control is demonstrated over the toroidal modal structure of the imposed magnetic perturbations of the helicity injection system, allowing a new class of discharges to be studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of divertor electrical drifts on particle distribution and detachment near the divertor target plate in DIII-D

Strong impacts of drifts on the divertor plasma in–out asymmetry and detachment are demonstrated in DIII-D with an open divertor configuration. For forward toroidal field, B T , i.e., with the ion B × ∇B drift toward the divertor, the particle flux to the inner divertor, as represented by the Langmuir probe measured ion saturation current (J sat ), exhibits a double peak structure, with electron temperature, lower at the inner target. Reversing the B T direction reverses both the radial and poloidal E × B flows, leading to a broad particle flux profile in the outboard scrape-off layer (SOL) with a similar double-peak structure to that observed at the inner target with forward B T . The correlation of a double peak structure with divertor temperature profiles confirms physical coupling between the drift flow and sheath boundary condition and their strong impact on divertor profiles. In addition, under reversed B T conditions, increasing the density flattens the target temperature profile. However, J sat remains high away from the strike point, rendering it difficult to achieve an “effective” detached plasma, i.e., with effective reduction in both peak heat flux and peak temperature (in the far SOL). In contrast, divertor detachment with a cold and flat temperature profile can be achieved at both target plates with the forward B T .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Mode spectrum characteristics and onset of the low-shear MHD stability regime

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q 0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q 0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q 0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Saturation physics of threshold heat-flux reduction

The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. Here, the heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Generalized Boozer coordinates: A natural coordinate system for quasisymmetry

We prove the existence of a straight-field-line coordinate system we call generalized Boozer coordinates. This coordinate system exists for magnetic fields with nested toroidal flux surfaces provided ∮ d l / B   ( j · ∇ ψ ) = 0 , where symbols have their usual meaning, and the integral is taken along closed magnetic field lines. All quasisymmetric fields, regardless of their associated form of equilibria, must satisfy this condition. This coordinate system presents itself as a convenient form to describe general quasisymmetric configurations and their properties. Insight can be gained analytically into the difference between strong and weak forms of quasisymmetry, as well as axisymmetry, and the interaction of quasisymmetry with different forms of equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Weakly quasisymmetric near-axis solutions to all orders

Here, we show that the equations satisfied by weakly quasisymmetric magnetic fields can be solved to arbitrarily high order in powers of the distance from the magnetic axis. This demonstration does not consider force balance. The existence of solutions requires an appropriate choice of parameters, most notably the toroidal current or rotational transform profiles. We do not prove that the expansion converges (it is likely divergent but asymptotic), and thus, the demonstration here should not be taken as definitive proof of the existence of global solutions. Instead, we provide a systematic construction of solutions to an arbitrarily high order.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Assessing physics of ion temperature gradient turbulence via hierarchical reduced-model representations

In this work, the saturation physics of ion temperature gradient (ITG) turbulence is probed by studying how amplitudes and scalings with key parameters vary across a hierarchy of reduced models. The models derive from nonlinear fluid equations for toroidal ITG turbulence under approximations to the mode coupling interactions in wavenumber space and the representation of turbulent decorrelation. Mode coupling approximations include local-in-wavenumber treatments like the spectral density of flux in quasilinear theory, a truncation to three nonlinearly interacting waves, and the interactions in a cascade to high radial wavenumber mediated by a single zonal flow. Turbulent decorrelation treatments are based on the triplet correlation time with and without eddy damping. Model fidelity is assessed by the scalings and magnitudes of the squared amplitudes of unstable mode, stable mode, and zonal flow with respect to the flow-damping rate and temperature gradient. It is shown that all models reproduce fundamental scalings, provided they incorporate the coupling of unstable mode, stable mode, and zonal flow. Accurate amplitude prediction requires eddy damping in the triplet correlation time and proper representation of the zonal-flow drive by interactions associated with the radial wavenumber cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A survey of pedestal magnetic fluctuations using gyrokinetics and a global reduced model for microtearing stability

This article presents a global reduced model for slab-like microtearing modes (MTMs) in the H-mode pedestal, which reproduces distinctive features of experimentally observed magnetic fluctuations, such as chirping and discrete frequency bands at noncontiguous mode numbers. Our model, importantly, includes the global variation of the diamagnetic frequencies, which is necessary to reproduce the experimental observations. The key insight underlying this model is that MTM instability is enabled by the alignment of a rational surface with the peak in the profile of the diamagnetic frequency. Conversely, MTMs are strongly stabilized for toroidal mode numbers for which these quantities are misaligned. This property explains the discrete fluctuation bands in several DIII-D and JET discharges, which we survey using our reduced model in conjunction with global gyrokinetic simulations. A fast yet accurate reduced model for MTMs enables rapid interpretation of magnetic fluctuation data from a wide range of experimental conditions to help assess the role of MTM in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation

A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗