Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma theory”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Improving the stellarator through advances in plasma theory

Improvements to the stellarator concept can be realized through advancements in theoretical and computational plasma physics. Herein, recent advances are reported in the topical areas of: (1) improved energetic ion confinement, (2) the impact of three-dimensional (3D) shaping on turbulent transport, (3) reducing coil complexity, (4) novel optimization and design methods, and (5) computational magnetohydrodynamic tools. Furthermore, these advances enable the development of new stellarator configurations with improved confinement properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2020

Dissertation Proposal: Innovations in Pulsed Power and Plasma Science Theory Simulation and Experiment

Pulsed power and plasma physics are topics of great study at both Sandia National Laboratories (SNL or Sandia) and the University of New Mexico (UNM). The goal of this research is to further knowledge and understanding of these fields using the resources of both SNL and UNM in three ways. The first way is through the comprehension, application, and testing of theory. Reading and analytically deriving theoretical solutions of problems both real-world and simplified will allow for a fresh perspective and the furthering of the theory. One such theory is Ottinger's generalized theory for voltage measurement in magnetically insulated transmission lines (MITLs). By working through the math, a deeper understanding of the theory is gained from which one may add more physically accurate and/or more detailed physics into the theory. Additionally, understanding the theory lays a good foundation from which one can analyze, test, and compare results to the theory in the following two ways that will advance the fields of pulsed power and plasma physics. The second way is through the modeling and simulation of real-world and simplified problems that utilize and test the afore mentioned theories. Theory can be applied to a simulation domain by using the unstructured time-domain electromagnetic (UTDEM) codes EMPHASIS and EMPIRE as well as the physical modeling software CUBIT, all of which were developed at SNL. Problems such as the modeling and design of the extended MITL on HERMES III, the understanding of space-charge-limited emission from vacuum cathodes, and the interaction between a relativistic electron beam and an ideal gas can all be modeled, simulated, and analyzed with this set of codes. Here the advantage is three-fold. Firstly, theory that describes our understanding of these problems can be put to the test and advanced through iterative simulation and analysis. Secondly, the understanding of these problems will have a positive impact on national security through the advancement of the technological capability of the United States of America. Thirdly, and not unrelated to the prior advantage, is the validation and verification of EMPIRE and EMPHASIS. This segues into the third way, which is through experiment and the comparison of experiment to simulated and theoretical results. Performing experimental comparisons completes the scientific method and grounds all of the work in reality. Being able to physically test theory and simulation is necessary for any real conclusions to be drawn. Another advantage for carrying out experimental work is to advance the physical testing capabilities of SNL. Several systems will be developed and tested through the course of this work that positively impact technological advancement of Sandia National Labs. Lastly, all of the above work will converge to yield a well-rounded perspective that ties the three categories of research together.

42 ENGINEERING↗
DOE OSTI2025

Closure and transport theory for plasmas with multiple ion species

This report summarizes a DOE-funded research project on closure and transport theory for plasmas with multiple ion species. The work developed accurate moment-based closure models across the full range of collisionality, incorporating both kinetic and collisional effects. Key contributions include high-fidelity closure relations, multi-temperature models, and transport coefficients applicable to fusion, space, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Small-amplitude Compressible Magnetohydrodynamic Turbulence Modulated by Collisionless Damping in Earth’s Magnetosheath: Observation Matches Theory

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales and affects energetic particle transport. Recent advances in the understanding of compressible magnetohydrodynamic (MHD) turbulence demonstrate the important role of damping in shaping energy distributions on small scales, yet its observational evidence is still lacking. This study provides the first observational evidence of substantial collisionless damping (CD) modulation on the small-amplitude compressible MHD turbulence cascade in Earth’s magnetosheath using four Cluster spacecraft. Based on an improved compressible MHD decomposition algorithm, turbulence is decomposed into three eigenmodes: incompressible Alfvén modes and compressible slow and fast (magnetosonic) modes. Our observations demonstrate that CD enhances the anisotropy of compressible MHD modes because CD has a strong dependence on wave propagation angle. The wavenumber distributions of slow modes are mainly stretched perpendicular to the background magnetic field (B 0 ) and weakly modulated by CD. In contrast, fast modes are subjected to a more significant CD modulation. Fast modes exhibit a weak, scale-independent anisotropy above the CD truncation scale. Below the CD truncation scale, the anisotropy of fast modes enhances as wavenumbers increase. As a result, fast-mode fractions in the total energy of compressible modes decrease with the increase of perpendicular wavenumber (to B 0 ) or wave propagation angle. Our findings reveal how the turbulence cascade is shaped by CD and its consequences for anisotropies in the space environment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quasilinear theory of general electromagnetic fluctuations including discrete particle effects for magnetized plasmas: General analysis

The general quasilinear Fokker–Planck kinetic equation for the gyrophase-averaged plasma particle distribution functions in magnetized plasmas is derived, making no restrictions on the energy of the particles and on the frequency of the electromagnetic fluctuations and avoiding the often made Coulomb approximation of the electromagnetic interactions. Here, the inclusion of discrete particle effects breaks the dichotomy of nonlinear kinetic plasma theory divided into the test particle and the test fluctuation approximation because it provides expression of both the non-collective and collective electromagnetic fluctuation spectra in terms of the plasma particle distribution functions. Within the validity of the quasilinear approach, the resulting full quasilinear transport equation can be regarded as a determining nonlinear equation for the time evolution of the plasma particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transport Barriers in magnetized plasmas- general theory with dynamical constraints

Abstract A fundamental dynamical constraint—that fluctuation induced charge-weighted particle flux must vanish- can prevent instabilities from accessing the free energy in the strong gradients characteristic of Transport Barriers (TBs). Density gradients, when large enough, lead to a violation of the constraint and hence preclude unstable modes and turbulent transport. This mechanism, then, broadens the class of configurations (in magnetized plasmas) where these high confinement states can be formed and sustained. The need for velocity shear, the conventional agent for TB formation, is obviated. The most important ramifications of the constraint is to permit a charting out of the domains conducive to TB formation and hence to optimally confined fusion worthy states; the detailed investigation is conducted through new analytic methods and extensive gyrokinetic simulations.

Physics↗
DOE OSTI2023

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Magnetic fusion experiments have the largest gradients in the world that are in steady state. These gradients lead to instabilities, and thus, most of plasma theory over the past four decades have been devoted to the development of instability theory and their use in interpreting plasma results. Instabilities in tokamaks lead to 3 different phenomenological outcomes: coherent saturation, turbulent saturation, and sudden relaxation (which may include disruptions). For long-wavelength instabilities, the NIMROD extended MHD code has emerged as an important tool for understanding tokamak instabilities. Because of the long history of NIMROD, a new version of it has been started at Tech-X to be able to address multi-species capabilities as well as exploit modern GPU systems. The new version builds on the previous version and improves the workflow by enabling new equilibria from experiments to be resolved. The fundamental formulation of the equations underlying NIMROD has also been explored in multiples ways. First, understanding instabilities in tokamaks through numerical methods is aided by being able to understand the rich history of analytic studies. A paper was published which aid in understanding the literature by simplifying some of the analytic machinery inherent in these studies. Extended MHD, like gyrokinetics, are quasineutral models. From a theoretical point of view, plasma quasineutrality approximation is best understood as the same as the magneto-quasistatic approximation of the Maxwell equations. Finally, the standard model of tokamak theory is that of instability theory. A simple dynamical systems model has been developed to better illustrate the strengths and weaknesses of this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗
DOE OSTI2024

Overview of T and D–T results in JET with ITER-like wall

In 2021 JET exploited its unique capabilities to operate with T and D–T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D–T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements—new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors—as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D–T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D–T plasmas to date and expanding our understanding of isotopes and D–T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D–T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D–T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D–T operations is presented in (King et al Nucl. Fusion submitted).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Resonance cones in cold plasma: Origin, singularities, and power flow

In magnetized tenuous plasma, typical at the plasma edge of fusion devices, a nearly electrostatic wave mode with relatively enhanced electric field can propagate along a specific angle with the magnetic field. For this characteristic, it is known as a “resonance cone.” For instance, these waves can be excited by radio frequency antennas in the ion-cyclotron and lower-hybrid range of frequencies. We consider the resonance cones emitted by idealized spatially extended sources. In 2D, we use a novel geometric construction which generalizes the d'Alembert solution to curved boundaries/moving sources, and show, for the first time, that singular electric fields arise under these conditions, thereby bringing the resonance cones in line with the other resonances of the cold plasma theory. Still in 2D, we give an expression for the amount of power radiated by resonance cones in terms of surface quantities on the source, which is finite despite the singular electric field. In conclusion, we generalize the conclusions regarding the presence and location of singular electric fields to the 3D electromagnetic case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Integrated System and Application Continuous Performance Monitoring and Analysis Capability

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia's Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia's ElectroMagnetic Plasma In Realistic Environments ( EMPIRE ) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia's HPC users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Integrated System and Application Continuous Performance Monitoring and Analysis Capability (Final)

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia’s Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia’s ElectroMagnetic Plasma In Realistic Environments (EMPIRE) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia’s HPC users.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗