Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma boundary 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.

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

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2022

Advances in prediction of tokamak experiments with theory-based models

The successful validation of theory-based models of transport, MHD stability, heating and current drive, with tokamak measurements over the last 20 years, has laid the foundation for a new era where these models can be routinely used in a "predict first" approach to design and predict the outcomes of experiments on tokamaks today. The capability to predict the plasma confinement and core profiles with a quantified uncertainty, based on a multi-machine, international, database of experience, will provide confidence that a proposed discharge will remain within the operational limits of the tokamak. Developing this predictive capability for the first generation of burning plasma devices, beginning with ITER, and progressing to tokamak demonstration reactors, is a critical mission of fusion energy research. Major advances have been made implementing this predict first methodology on today's tokamaks. In this work, an overview of several of these recent advances will be presented, providing the integrated modeling foundations of the experimental successes. The first steps to include boundary plasmas, and tokamak control systems, have been made. A commitment to predicting experiments as part of the planning process is needed in order to collect predictive accuracy data and evolve the models and software into a robust whole discharge pulse design simulator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Pushing the Boundaries of Relativistic Fluid Dynamics: From Particle Colliders to Neutron Star Mergers

This project aimed at establishing connections between fluid dynamics, kinetic theory, and string theory to develop new tools to characterize the far-from-equilibrium properties of relativistic systems, which are instrumental in understanding the emergence of liquid-like behavior of the quark-gluon plasma formed in heavy-ion collisions. This project investigated the fundamental properties of viscous fluids in the full nonlinear regime to pave the way for the systematic inclusion of viscous effects in future numerical simulations of neutron star mergers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Theoretical Investigation of the High-Altitude Cusp Region using Observations from Interball and ISTP Spacecraft

A fundamental goal of magnetospheric physics is to understand the transport of plasma through the solar wind-magnetosphere-ionosphere system. To attain such an understanding, we must determine the sources of the plasma, the trajectories of the particles through the magnetospheric electric and magnetic fields to the point of observation, and the acceleration processes they undergo enroute. This study employed plasma distributions observed in the near-Earth plasma sheet by Interball and Geotail spacecraft together with theoretical techniques to investigate the ion sources and the transport of plasma. We used ion trajectory calculations in magnetic and electric fields from a global Magnetohydrodynamics (MHD) simulation to investigate the transport and to identify common ion sources for ions observed in the near-Earth magnetotail by the Interball and Geotail spacecraft. Our first step was to examine a number of distribution functions and identify distinct boundaries in both configuration and phase space that are indicative of different plasma sources and transport mechanisms. We examined events from October 26, 1995, November 29-30, 1996, and December 22, 1996. During the first event Interball and Geotail were separated by approximately 10 R(sub E) in z, and during the second event the spacecraft were separated by approximately 4(sub RE). Both of these events had a strong IMF By component pointing toward the dawnside. On October 26, 1995, the IMF B(sub Z) component was northward, and on November 1-9-30, 1996, the IMF B sub Z) component was near 0. During the first event, Geotail was located near the equator on the dawn flank, while Interball was for the most part in the lobe region. The distribution function from the Coral instrument on Interball showed less structure and resembled a drifting Maxwellian. The observed distribution on Geotail, on the other hand, included a great number of structures at both low and high energies. During the third event (December 22, 1996) both spacecraft were in the plasma sheet and were separated bY approximately 20 R(sub E) in the y direction. During this event the IMF was southward.

Ashour-Abdalla, Maha↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Model predictive control of KSTAR equilibrium parameters enabled by TRANSP

Due to the complex behavior of tokamak plasmas and the importance of optimizing performance while avoiding instabilities and machine limits, plasma control algorithms are becoming increasingly dependent on sophisticated model-based control approaches. It is anticipated that the use of integrated modeling codes in the model-based control design process will reduce the amount of experimental time needed to implement new control algorithms by facilitating development of control-oriented models and enabling higher-fidelity closed-loop simulations. In this work, a reduced model is developed from a series of TRANSP simulations and is used to develop a model predictive control (MPC) algorithm for controlling important equilibrium parameters in KSTAR [1] discharges. The control algorithm uses the KSTAR neutral beam injection system and the target plasma current and plasma boundary as actuators, and optimizes the plasma stored energy, loop voltage, and internal inductance while avoiding constraints that could lead to disruptions. Higher fidelity testing of the control algorithm is performed using a flexible framework for enabling external processes to actively control plasma parameters in TRANSP simulations. Furthermore, closed-loop simulations demonstrate the ability of the control algorithm to respond to disturbances in density and confinement, handle actuator failures, and move the discharge to high non-inductive fraction conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dissipation of electron-beam-driven plasma wakes

Metre-scale plasma wakefield accelerators have imparted energy gain approaching 10 gigaelectronvolts to single nano-Coulomb electron bunches. To reach useful average currents, however, the enormous energy density that the driver deposits into the wake must be removed efficiently between shots. Yet mechanisms by which wakes dissipate their energy into surrounding plasma remain poorly understood. Here, we report picosecond-time-resolved, grazing-angle optical shadowgraphic measurements and large-scale particle-in-cell simulations of ion channels emerging from broken wakes that electron bunches from the SLAC linac generate in tenuous lithium plasma. Measurements show the channel boundary expands radially at 1 million metres-per-second for over a nanosecond. Simulations show that ions and electrons that the original wake propels outward, carrying 90 percent of its energy, drive this expansion by impact-ionizing surrounding neutral lithium. The results provide a basis for understanding global thermodynamics of multi-GeV plasma accelerators, which underlie their viability for applications demanding high average beam current.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

NSTX-U liquid metal core-edge facility (LMCE)

NSTX-U/LMCE will provide a unique and world-leading research facility to address the primary challenge to delivering economic and timely magnetic fusion energy, namely the need to develop a power and particle exhaust and first-wall system that can withstand very high edge heat fluxes, maximize energy confinement, and avoid the production of large masses of solid eroded first-wall material. The NSTX-U/LMCE facility will assess the ability of liquid metals (LMs) – especially liquid lithium – to provide a new boundary condition for magnetic fusion systems, to extend the lifetime of the plasma facing components (PFCs) and improve core plasma confinement. Such capability is needed to establish the basis for next-step fusion facilities including fusion pilot plants, and to maintain U.S. world leadership in core-edge integration research. NSTX-U/LMCE will leverage the ability to generate very high divertor perpendicular heat flux q⊥ ~ 100MW/m 2 , extensive diagnostics, and liquid-metal-applicable infrastructure of NSTX-U. NSTX-U/LMCE will provide access to a high-confinement plasma core with majority self-driven plasma current, the flexibility to test a range of liquid metal divertor concepts, access to a range of separatrix collisionalities (from high to very low), and the ability to controllably vary the first-wall temperature to vary the plasma- wall interaction physics on liquid lithium components. Further, NSTX-U/LMCE will utilize more reactor-relevant high-Z refractory-metal PFC substrates. With these capabilities the NSTX-U/LMCE facility will explore the full continuum of core-edge solutions ranging from high core radiated power, to conditions with radiative losses concentrated in the scrape-off layer (SOL), and ultimately low recycling conditions. The low collisionality SOL that may be accessible in the low recycling regime is relatively unexplored and will require a kinetic treatment of the edge, which can be addressed theoretically, and with experiments in LTX-β. Additional smaller-scale preparatory R&D facilities will be required to reduce the risk of premature technical/engineering failure of liquid metal systems implemented in NSTX-U. The NSTX-U/LMCE facility aligns very well with recommendations in the FESAC Long-Range Plan and NASEM Pilot Plant reports and the Bold Decadal Vision, will be unique in the world program throughout the next decade, and is garnering private company interest in utilizing NSTX-U/LMCE for development of LM PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Real-time plasma equilibrium reconstruction and shape control for the MAST Upgrade tokamak

Real-time magnetic control has been developed to deliver precise control of multiple plasma shape parameters for advanced divertor configurations, including double-null, Super-X, X-point target and X-divertor for the first time on the MAST Upgrade (MAST-U) spherical tokamak. Successful real-time magnetic equilibrium control of different plasma shape variables has been accomplished in the 2022–2023 MAST-U experimental campaign for the advanced MAST-U divertor configurations. Application of the MAST-U boundary reconstruction algorithm, LEMUR, is described and compared with off-line equilibrium reconstruction and diagnostic measurements. The process of design and verification of the axisymmetric magnetic control schemes using a suite of control analysis tools (known collectively as TokSys) is also described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

A Review of Low Frequency Electromagnetic Wave Phenomena Related to Tropospheric-Ionospheric Coupling Mechanisms

Investigation of coupling mechanisms between the troposphere and the ionosphere requires a multidisciplinary approach involving several branches of atmospheric sciences, from meteorology, atmospheric chemistry, and fulminology to aeronomy, plasma physics, and space weather. In this work, we review low frequency electromagnetic wave propagation in the Earth-ionosphere cavity from a troposphere-ionosphere coupling perspective. We discuss electromagnetic wave generation, propagation, and resonance phenomena, considering atmospheric, ionospheric and magnetospheric sources, from lightning and transient luminous events at low altitude to Alfven waves and particle precipitation related to solar and magnetospheric processes. We review in situ ionospheric processes as well as surface and space weather phenomena that drive troposphere-ionosphere dynamics. Effects of aerosols, water vapor distribution, thermodynamic parameters, and cloud charge separation and electrification processes on atmospheric electricity and electromagnetic waves are reviewed. We also briefly revisit ionospheric irregularities such as spread-F and explosive spread-F, sporadic-E, traveling ionospheric disturbances, Trimpi effect, and hiss and plasma turbulence. Regarding the role of the lower boundary of the cavity, we review transient surface phenomena, including seismic activity, earthquakes, volcanic processes and dust electrification. The role of surface and atmospheric gravity waves in ionospheric dynamics is also briefly addressed. We summarize analytical and numerical tools and techniques to model low frequency electromagnetic wave propagation and solving inverse problems and summarize in a final section a few challenging subjects that are important for a better understanding of tropospheric-ionospheric coupling mechanisms.

Simoes, Fernando↗
DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Anticathode effect on electron kinetics in electron beam generated E × B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E x B) fields are promising for low damage processing of materials with applications to microelectronics and quantum information systems. In cylindrical e-beam E x B plasmas, radial confinement of electrons and ions is achieved by an axial magnetic field and radial electric field, respectively. To control the axial confinement of electrons, such e-beam generated plasma sources may incorporate a conducting boundary known as an anticathode, which is placed on the axially opposite side of the plasma from the cathode. In this work, it is shown that varying the anticathode voltage bias can control the degree to which the anticathode collects or repels incident electrons, allowing control of warm electron (electron energies in 10–30 eV range) and beam electron population confinement. It is suggested that the effect of the anticathode bias on the formation of these distinct electron populations is also associated with the transition between weak turbulence and strong Langmuir turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Three-dimensional cross-field flows at the plasma-material interface in an oblique magnetic field

This article describes experimental evidence that the magnetic presheath is a fully three-dimensional structure modified by ion–neutral collisions. Velocity distributions of both ions and neutrals, obtained via laser-induced fluorescence, show that cross field ion drifts do not result from entrainment of ions in a flowing neutral background. Ion flows parallel to E x B arise and accelerate to as much as 0.2c s within several ion gyroradii of the boundary surface, where c s is the sound speed. Within measurement resolution, the onset of the E x B aligned flow occurs at the same distance to the surface that ions begin to deflect from travel along magnetic field lines. Collisional fluid and particle-in-cell simulations of the boundary region are compared to the experimental measurements. We find that, in contrast to the classical collisionless Chodura model, collisional effects between the ions and the non-flowing neutral population are essential to quantitatively predict the observed ion drift velocities. No momentum coupling between ions and neutrals, separable from noise and other effects, is observed in either signal. Here, we discuss several explanations and implications of this observation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2026

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. II. Interaction via Landau resonance

We report that conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes driven by the Landau resonance with super-Alfvénic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite ω/ω ci and |k ∥ /k ⊥ | is retained and found to be responsible for the existence of the GAE instability via this resonance. Finally, encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗