Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal magnetic fields”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonlinear coupling of tearing modes in reversed field pinch plasmas with stepped pressure profiles

A theory of the three-wave coupling of triplets of tearing modes in toroidal pinches [i.e., either reversed field pinches (RFPs) or tokamaks] was proposed by R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)]. However, this theory only applies to toroidal pinches with negligible equilibrium plasma pressure gradients. Such a limitation is particularly inappropriate to RFPs. This paper generalizes the analysis of R. Fitzpatrick [Phys. Plasmas 6, 1168 (1999)] in order to take the equilibrium pressure gradient into account. However, for the sake of simplicity, a stepped pressure profile, rather than a continuous profile, is employed. In the limit in which the number of steps becomes very large, the results obtained from the generalized theory are presumably equivalent to those that would have been achieved using a continuous pressure profile. The generalized theory is used to investigate the formation of the characteristic toroidally localized pattern of phase-locked m = 1 and m = 0 tearing modes in RFP plasmas that is known as the “slinky” pattern. The incorporation of the equilibrium plasma pressure into the analysis is found to be of crucial importance when determining the properties of the pattern. This is the case because the plasma pressure controls the number of unstable m = 1 and m = 0 tearing modes, and also significantly affects the strength of three-wave coupling, as well as the phase relation between the phase-locked m = 1 and m = 0 modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2021

Upgrades to the ion cyclotron emission diagnostic on the DIII-D tokamak

The ion cyclotron emission diagnostic on the DIII-D tokamak comprises seven single-turn loops that measure high-frequency (1–100 MHz) magnetic field fluctuations that are often excited by energetic particles in the plasma. The raw voltage signals induced in the loops in response to these fluctuations travel through a series of cables, isolation transformer DC blocks, low-pass filters, and finally a digitizer before being analyzed in frequency space. The diagnostic has been recently upgraded, most notably to include four additional graphite tile loops and a new eight-channel digitizer. The previous three loops are all on the low-field side of the tokamak. The measurement capabilities of the system have been expanded by the addition of a new horizontally oriented loop on the low-field side, an additional toroidal loop on the low-field side, and two toroidal loops on the high-field side. These loops will be used to provide approximate mode polarization, improved toroidal mode number calculations, and information on modes in inward-shifted plasmas, respectively.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Solution of the linear wave-particle kinetic equation for global modes of arbitrary frequency in a tokamak

The linear response of a plasma to perturbations of arbitrary frequency and wavelength is derived for any axisymmetric magnetized toroidal plasma. An explicit transformation to action-angle coordinates is achieved using orthogonal magnetic coordinates and the Littlejohn Lagrangian, establishing the validity of this result to arbitrary order in normalized Larmor radius. The global resonance condition for compressional modes is clarified in more detail than in previous works, confirming that the poloidal orbit-average of the cyclotron frequency gives the desired result at lowest order in Larmor radius. The global plasma response to the perturbation at each resonance is captured by a poloidal and gyroaverage of the perturbing potential. A “global gyroaveraging” of the potential is a natural by-product of this analysis which takes into account the changing of the magnetic field over an orbit. The resonance condition depends on two arbitrary integers which completely separately capture the effects poloidal non-uniformity and finite Larmor radius in generating sidebands. We learn that poloidal sidebands generated for compressional modes are dominated by the change in gyrofrequency over the orbit, which is very different to shear modes where the gyrofrequency only contributes via a finite Larmor radius effect. This increases the number of bounce harmonics required to compute the linear drive, giving a more complicated resonance map. An example calculation is given comparing resonance of shear and compressional modes in a published DIII-D case.

Compressional↗
DOE OSTI2025

Microtearing stability and turbulence in the pedestal: Linear gyrokinetics, reduced models, and nonlinear turbulent transport

Microtearing modes can play a crucial role in electron heat transport in tokamak plasmas, affecting both energy confinement and overall performance. This study investigates microtearing modes (MTM) stability and turbulence in a JET pedestal through gyrokinetic simulations using the Gene code, complemented by a reduced eigenvalue model. The focus is on how MTM properties depend on key plasma parameters, including collisionality and plasma beta β—the ratio of plasma pressure to magnetic pressure—the normalized toroidal wavenumber k y ρ s ⁠, where ρ s denotes the ion sound gyroradius (typically a few millimeters in edge plasmas) and isotope mass. Collisionality enhances MT growth rates, while increasing β leads to a shift from MTMs to kinetic-ballooning modes, typically for k y ρ s ⁠, where ρ s ≲ 0.2⁠. A purely collisionless branch of MTMs persists at low k y ρ s ⁠, where ρ s with distinctive properties including non-negligible particle flux and ion thermal transport. Isotope mass scans reveal modest reduction of MTM growth rates as ion mass decreases. Nonlinear simulations produce experimentally relevant transport levels. Numerical experiments turning off zonal flows and fields identify the critical role of zonal flows and zonal fields in regulating MTM turbulence. Their removal leads to a significant increase in electron heat flux. These findings provide new insight into MTM-driven transport and its impact on tokamak confinement and lay a foundation for reduced modeling and predictive capabilities.

Electrostatics↗
DOE OSTI2024

Enhanced plasma performance in C-2W advanced beam-driven field-reversed configuration experiments

TAE Technologies’ fifth-generation fusion device, C-2W (also called ‘Norman’), is the world’s largest compact-toroid device and has made significant progress in field-reversed configuration (FRC) plasma performance. C-2W produces record breaking, macroscopically stable, high-temperature advanced beam-driven FRC plasmas, dominated by injected fast particles and sustained in steady state, which is primarily limited by neutral-beam (NB) pulse duration. The NB power supply system has recently been upgraded to extend the pulse length from 30 ms to 40 ms, which allows for a longer plasma lifetime and thus better characterization and further enhancement of FRC performance. An active plasma control system is routinely used in C-2W to produce consistent FRC performance as well as for reliable machine operations using magnet coils, edge-biasing electrodes, gas injection and tunable-energy NBs. Google’s machine learning framework for experimental optimization has also been routinely used to enhance plasma performance. D

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electronic scattering off a magnetic hopfion

We study scattering of itinerant electrons off a magnetic hopfion in a three-dimensional metallic magnet described by a magnetization vector S(r). A hopfion is a confined topological soliton of S(r) characterized by an emergent magnetic field Bγ(r)≡ε αβγ S∙(∇ α S×∇ β S)/4≠0 with vanishing average value [B(r)]=0. We evaluate the scattering amplitude in the opposite limits of large and small hopfion radius R using the eikonal and Born approximations, respectively. In both limits, we find that the scattering cross section contains a skew-scattering component giving rise to the Hall effect within a hopfion plane. That conclusion contests the popular notion that the Hall effect in noncollinear magnetic structures necessarily implies [B(r)]≠0. In the limit of small hopfion radius pR<<1, we expand the Born series in powers of momentum p and identify different expansion terms corresponding to the hopfion anisotropy, toroidal moment, and skew scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Optimization of confinement in a toroidal plasma subject to strong radial electric fields

The toroidal ring of plasma contained in the NASA Lewis Bumpy Torus facility may be biased to positive or negative potentials approaching 50 kV by applying dc voltages of the respective polarity to 12 or fewer midplane electrode rings. The radial electric fields, which are responsible for raising the ions to high energies by E x B/B-squared drift, then point out of or into the plasma. A preliminary report is given on the identification and optimization of independent variables which affect the ion density and confinement time in the Bumpy Torus plasma. The independent variables include the polarity, position, and number of the midplane electrode rings, the method of gas injection, and the polarity and strength of a weak vertical magnetic field. Some characteristic data taken under conditions where most of the independent variables were optimized are presented. The highest value of the electron number density on the plasma axis is 3.2 trillion per cu cm, the highest ion heating efficiency is 47%, and the longest particle containment time is 2.0 msec.

Roth, J. R.↗
NASA NTRS1990

Observational features of field line resonances excited by solar wind pressure variations on 4 September 1984

The purpose of the study is to establish the most probable excitation mechanism of the magnetic storm occurred after an inverse sudden impulse on September 4, 1984. Geomagnetic pulsations in the Pc5-frequency range observed at magnetometer stations are evaluated. Attention is focused on two events of the enhanced activity: for the first one, conjugate observations on the ground are assessed and then compared with satellite-based observations on adjacent field lines; for the second event two hours later, data from an extended azimuthal range is employed. It is pointed out that the observations are consistent with the theory of filed-line resonance, and may be interpreted as excitations caused by pressure variations in the solar wind. Both magnetopause-surface waves and cavity resonances are excited; the cavity mode drives toroidal field-line oscillations at locations where its frequency matches the resonance frequency of the field lines.

Warnecke, J.↗
NASA NTRS2006

Constraining Numerical Geodynamo Modeling with Surface Observations

Numerical dynamo solutions have traditionally been generated entirely by a set of self-consistent differential equations that govern the spatial-temporal variation of the magnetic field, velocity field and other fields related to dynamo processes. In particular, those solutions are obtained with parameters very different from those appropriate for the Earth s core. Geophysical application of the numerical results therefore depends on correct understanding of the differences (errors) between the model outputs and the true states (truth) in the outer core. Part of the truth can be observed at the surface in the form of poloidal magnetic field. To understand these differences, or errors, we generate new initial model state (analysis) by assimilating sequentially the model outputs with the surface geomagnetic observations using an optimal interpolation scheme. The time evolution of the core state is then controlled by our MoSST core dynamics model. The final outputs (forecasts) are then compared with the surface observations as a means to test the success of the assimilation. We use the surface geomagnetic data back to year 1900 for our studies, with 5-year forecast and 20-year analysis periods. We intend to use the result; to understand time variation of the errors with the assimilation sequences, and the impact of the assimilation on other unobservable quantities, such as the toroidal field and the fluid velocity in the core.

Kuang, Weijia↗
DOE OSTI2023

Required toroidal confinement for fusion and omnigeneity

Deuterium–tritium (DT) burning requires a long energy confinement times compared to collision times, so the particle distribution functions must approximate local Maxwellians. Non-equilibrium thermodynamics is applicable, which gives relations among transport, entropy production, the collision frequency, and the deviation from a Maxwellian. The distribution functions are given by the Fokker–Planck equation, which is an advection–diffusion equation. A large hyperbolic operator, the Vlasov operator with the particle trajectories as its characteristics, equals a small diffusive operator, the collision operator. The collisionless particle trajectories would be chaotic in stellarators without careful optimization. This would lead to rapid entropy production and transport—far beyond what is consistent with a self-sustaining DT burn. Omnigeneity is the weakest general condition that is consistent with a sufficiently small entropy production associated with the thermal particle trajectories. Omnigeneity requires that the contours of constant magnetic field strength be unbounded in at least one of the two angular coordinates in magnetic surfaces and that there be a symmetry in the field-strength wells along the field lines. Even in omnigenous plasmas, fluctuations due to microturbulence can produce chaotic particle trajectories and the gyro-Bohm transport is seen in many stellarator and tokamak experiments. The higher the plasma temperature above 10 keV, the smaller the transport must be compared to gyro-Bohm for a self-sustaining DT burn. The hot alphas of DT fusion heat the electrons. When the ion–electron equilibration time is long compared to the ion energy confinement time, a self-sustaining DT burn is not possible, which sets a limit on the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of aspect ratio on plasma response to resonant magnetic perturbations in tokamak devices

A systematic numerical study is carried out, computing and comparing the plasma response to the resonant magnetic perturbation (RMP) field, applied for controlling edge localized modes (ELMs), in a series of tokamak plasmas with varying aspect ratio and utilizing the MARS-F code. The aspect ratio is scanned either by varying the plasma major radius at a fixed minor radius or by varying the latter while fixing the former. Both approaches yield similar results when compared in terms of quantities with proper normalizations. In general, a non-monotonic dependence of the resonant response field (normalized by the vacuum counterpart) near the plasma edge is found with varying aspect ratio, indicating that a given ELM control coil current configuration strongly favors plasmas with a certain aspect ratio. This optimal aspect ratio, on the other hand, depends on the toroidal as well as poloidal (i.e., coil phasing) spectra of the applied RMP field. The equilibrium (edge) safety factor, the plasma shape, and the plasma toroidal flow are all fixed to ensure that the effects identified here are predominantly due to the plasma aspect ratio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
DOE OSTI2023

Loss of energetic ions due to n = 1 internal kink instability in HL-2M

Effects of three-dimensional (3D) perturbations due to an unstable n = 1 (n is the toroidal mode number) internal kink (IK) on the energetic particles (EPs) are systematically investigated for the HL-2M tokamak, utilizing the MARS-F/K code and a recently developed test particle tracing module. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code, is chosen for this study. In general, the 3D perturbation associated with an unstable IK is found to affect the EP drift orbit, confinement, and loss in HL-2M. The instability-induced EP loss fraction is found to be typically less than 10%, without counting for the prompt orbit loss associated with the 2D equilibrium field for counter-current particles. The latter reaches about 16% in HL-2M. For co-current EPs, a 100 G 3D magnetic field (inside the plasma) due to the IK does not induce any EP loss assuming a static perturbation. A sawtooth-like time-varying perturbation field, with the peak amplitude reaching 1000 G, can however produce about 30% loss for the co-current EPs in HL-2M. The majority of lost EPs tend to strike the lower divertor region, with a small fraction of particles striking the low-field side mid-plane region of the limiting surface.

36 MATERIALS SCIENCE↗