SEARCH · Search NASA

Results for “Fusion plasma physics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Plasma potential shaping using end-electrodes in the Large Plasma Device

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Nonlinear solution of classical three-wave interaction via finite-dimensional quantum model

The quantum three-wave interaction, the lowest-order nonlinear interaction in plasma physics, describes energy–momentum transfer between three resonant waves in the quantum regime. We describe how it may also act as a finite-degree-of-freedom approximation to the classical three-wave interaction in certain circumstances. By promoting the field variables to operators, we quantize the classical system, show that the quantum system has more free parameters than the classical system and explain how these parameters may be selected to optimize either initial or long-term correspondence. We then numerically compare the long-time quantum–classical correspondence far from the fixed point dynamics. We discuss the Poincaré recurrence of the system and the mitigation of quantum scrambling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Stellarator optimization with constraints

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a quasi-axisymmetric stellarator for favourable physics properties while preventing strong shaping of the plasma boundary, which can be difficult to create with external current sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation

Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared with unoptimised configurations. Although stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratios and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared with solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface, three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Modelling of the electron cyclotron emission burst from a laboratory tokamak plasma with loss-cone maser instability

The maser instability associated with the loss-cone distribution has been widely invoked to explain the radio bursts observed in the astrophysical plasma environment, such as aurora and corona. In the laboratory plasma of a tokamak, events reminiscent of these radio bursts have also been frequently observed as an electron cyclotron emission (ECE) burst in the microwave range (~2f ce near the last closed flux surface) during transient magnetohydrodynamic events. These bursts have a short duration of ~10 μs and display a radiation spectrum corresponding to a radiation temperature T e,rad of over 30 keV while the edge thermal electron temperature T e is only in the range of 1 keV. Suprathermal electrons can be generated through magnetic reconnection, and a loss-cone distribution can be generated through open stochastic field lines in the magnetic mirror of the near-edge region of a tokamak plasma. Radiation modelling shows that a sharp distribution gradient ∂f/∂v ⊥ > 0 at the loss-cone boundary can cause a negative absorption of ECE radiation through the maser instability. The negative absorption then amplifies the radiation so that the microwave intensity is significantly stronger than the thermal value. The significant T e,rad from the simulations suggests the potential role of the loss-cone maser instability in generating the ECE burst in a tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

MagNetUS: a magnetized plasma research ecosystem

MagNetUS is a network of scientists and research groups that coordinates and advocates for fundamental magnetized plasma research in the USA. Its primary goal is to bring together a broad community of researchers and the experimental and numerical tools they use in order to facilitate the sharing of ideas, resources and common tasks. Discussed here are the motivation and goals for this network and details of its formation, history and structure. An overview of associated experimental facilities and numerical projects is provided, along with examples of scientific topics investigated therein. Finally, a vision for the future of the organization is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Strong gradient effects on neoclassical electron transport and the bootstrap current in large aspect ratio tokamaks

Standard approaches to neoclassical theory do not extend into regions of strong gradients in tokamaks such as the pedestal and internal transport barriers. Here, we calculate the modifications to neoclassical electron physics inside strong gradient regions of large aspect ratio tokamaks in the banana regime. We show that these modifications are due to the different ion flow and the strong poloidal variation of the potential. We also provide a physical interpretation of the mechanisms that drive poloidal asymmetries and hence a poloidal electric field. We apply our model to two specific example cases of pedestal profiles, calculating the neoclassical electron flux and the bootstrap current. We find that, depending on the ion flow, weak gradient neoclassical theory overestimates or underestimates the neoclassical electron transport and the bootstrap current in regions with strong gradients. We show that the determination of the mean parallel flow is more complex than in weak gradient neoclassical theory. For vanishing turbulence, we can determine the radial electric field for a given flow profile in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modeling study of divertor particle and heat flux asymmetries for EAST H-mode discharges

The BOUT++ transport code is run to study the effects of plasma drifts on the divertor out-in asymmetries (DOIAs) of particle and heat fluxes and their decay widths for EAST lower single null H-mode discharges. The diamagnetic drift seems to have no effects on the DOIAs of total particle and heat fluxes due to its divergence-free nature. However, it could significantly increase the DOIAs of peak particle and heat fluxes and the flux decay widths. The E × B drift is found to induce a large plasma flow to the divertor region, enhancing the DOIAs of both total and peak particle and heat fluxes and the flux decay widths. Both the radial and poloidal components of the E × B drift are necessary in increasing the DOIAs, however, the radial E × B drift seems to play a more important role. The effects on the DOIAs caused by both diamagnetic and E × B drifts are reversed with the reverse of toroidal magnetic field. The heat flux decay width λ q and spreading width S q are important physical and engineering parameters for the divertors and could be obtained by fitting the heat flux profiles at divertor targets. The λ q at the outer target from the simulation case with all drifts could well match with the multi-machine scaling proposed by Eich and the DOIA of λ q is in reasonable agreement with the scaling proposed by Goldston.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Probing strong-field QED in beam-plasma collisions

Ongoing progress in laser and accelerator technology opens new possibilities in high-field science, notably to investigate the largely unexplored strong-field quantum electrodynamics (SFQED) regime where electron-positron pairs can be created directly from light-matter or even light-vacuum interactions. Laserless strategies such as beam-beam collisions have also been proposed to access the nonperturbative limit of SFQED. Here we report on a concept to probe SFQED by harnessing the interaction between a high-charge, ultrarelativistic electron beam and a solid conducting target. When impinging onto the target surface, the beam self fields are reflected, partly or fully, depending on the beam shape; in the rest frame of the beam electrons, these fields can exceed the Schwinger field, thus triggering SFQED effects such as quantum nonlinear inverse Compton scattering and nonlinear Breit-Wheeler electron-positron pair creation. Through reduced modeling and kinetic numerical simulations, we show that this single-beam setup can achieve interaction conditions similar to those envisioned in beam-beam collisions, but in a simpler and more controllable way owing to the automatic overlap of the beam and driving fields. This scheme thus eases the way to precision studies of SFQED and is also a promising milestone towards laserless studies of nonperturbative SFQED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Design, Manufacturing, and Installation of DIII-D Lower Hybrid Current Drive Outer Wall Waveguides and Expandable Vacuum Interface

High field side lower hybrid current drive (HFS LHCD) is promising method for efficient, off axis current drive. To operate this system, development of an expandable combination radio frequency (RF) and vacuum interface was necessary to allow the long (~2 m) waveguides to thermally expand with respect to the DIII-D 0 degree R-1 port extension without damage to the waveguide or excessive force on the vacuum seals which was previously observed when baking the waveguides during operations that were installed in 2024. Design of a custom bellows assembly with RF feedthrough that allows >20mm of waveguide thermal expansion is presented. Fabrication of this bellows assembly required simultaneously vacuum brazing multiple 304 stainless steel flanges to an oxygen free high conductivity copper waveguide section which required significant manufacturing development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Elucidating plasma dynamics in Hasegawa–Wakatani turbulence by information geometry

In this work, the impact of adiabatic electrons on drift-wave turbulence, modeled by the Hasegawa–Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system and represents the total number of statistically different states that a system evolves through in time. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ are computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behavior is identified by a rapid change in the information length. The effects of $A$ are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗