Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synchronization of dust acoustic waves in a forced Korteweg–de Vries–Burgers model

The synchronization of dust acoustic waves to an external periodic source is studied in the framework of a driven Korteweg–de Vries–Burgers equation that takes into account the appropriate nonlinear and dispersive nature of low-frequency waves in a dusty plasma medium. For a spatiotemporally varying source term, the system is shown to demonstrate harmonic (1:1) and superharmonic (1:2) synchronized states. Here, the existence domains of these states are delineated in the form of Arnold tongue diagrams in the parametric space of the forcing amplitude and forcing frequency and their resemblance to some past experimental results is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

In this paper, we consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Parasitic current collection by PASP Plus solar arrays

Solar cells at potentials positive with respect to a surrounding plasma collect electrons. Current is collected by the exposed high voltage surfaces: the interconnects and the sides of the solar cells. This current is a drain on the array power that can be significant for high-power arrays. In addition, this current influences the current balance that determines the floating potential of the spacecraft. One of the objectives of the Air Force (PL/GPS) PASP Plus (Photovoltaic Array Space Power Plus Diagnostics) experiment is an improved understanding fo parasitic current collection. We have done computer modeling of parasitic current collection and have examined current collection flight data from the first year of operations. Prior to the flight we did computer modeling to improve our understanding of the physical processes that control parasitic current collection. At high potentials, the current rapidly rises due to a phenomenon called snapover. Under snapover conditions, the equilibrium potential distribution across the dielectric surface is such that part of the area is at potentials greater than the first crossover of the secondary yield curve. Therefore, each incident electron generates more than one secondary electron. The net effect is that the high potential area and the collecting area increase. We did two-dimensional calculations for the various geometries to be flown. The calculations span the space of anticipated plasma conditions, applied potential, and material parameters. We used the calculations and early flight data to develop an analytic formula for the dependence of the current on the primary problem variables. The analytic formula was incorporated into the EPSAT computer code. EPSAT allows us to easily extend the results to other conditions. PASP Plus is the principal experiment integrated onto the Advanced Photovoltaic and Electronics Experiments (APEX) satellite bus. The experiment is testing twelve different solar array designs. Parasitic current collection is being measured for eight of the designs under various operational and environment conditions. We examined the current collected as a function of the various parameters for the six non-concentrator designs. The results are similar to those obtained in previous experiments and predicted by the calculations. We are using the flight data to validate the analytic formula developed. The formula can be used to quantify the parasitic current collected. Anticipating the parasitic current value allows the spacecraft designer to include this interaction when developing the design.

Davis, Victoria Ann↗
DOE OSTI2022

Ultralight millicharged dark matter via misalignment

We explore the cosmology and phenomenology of millicharged and millicharge-like dark matter with masses from 1 eV to 10 keV and charges of 10 -18 to 10 -14 . Dark matter in this mass range cannot be thermally produced, but can arise from non-thermal mechanisms. We propose a concrete model employing a spontaneously broken approximate global symmetry, in which millicharged dark matter is produced via the misalignment mechanism. We show that this production mechanism is cosmologically consistent and compatible with the observed dark matter abundance. This model can be implemented using either fundamental scalars or hidden-sector quarks, and coupled either to the Standard Model photon or to a hidden photon. We then consider the phenomenology of light millicharged dark matter, regardless of its cosmological origin, and determine the parameter space consistent with existing experiments and observations. A significant part of the new parameter space we consider may be accessible in the near future through direct deflection experiments, measurements of the cosmic microwave background blackbody spectrum, and future constraints on plasma instabilities due to dark matter self-interaction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1980

Plasma physics analysis of SERT-2 operation

An analysis of the major plasma processes involved in the SERT 2 spacecraft experiments was conducted to aid in the interpretation of recent data. A plume penetration model was developed for neutralization electron conduction to the ion beam and showed qualitative agreement with flight data. In the SERT 2 configuration conduction of neutralization electrons between thrusters was experimentally demonstrated in space. The analysis of this configuration suggests that the relative orientation of the two magnetic fields was an important factor in the observed results. Specifically, the opposed field orientation appeared to provide a high conductivity channel between thrusters and a barrier to the ambient low energy electrons in space. The SERT 2 neutralizer currents with negative neutralizer biases were up to about twice the theoretical prediction for electron collection by the ground screen. An explanation for the higher experimental values was a possible conductive path from the neutralizer plume to a nearby part of the ground screen. Plasma probe measurements of SERT 2 gave the clearest indication of plasma electron temperature, with normal operation being near 5 eV and discharge only operation near 2 eV.

Kaufman, H. R.↗
NASA NTRS1981

Waves in Space Plasmas Program

The Waves in Space Plasmas (WISP) program is a joint international effort involving instrumentation to be designed and fabricated by funding from NASA and the National Research Council of Canada. The instrumentation, with a tentatively planned payload for 1986, can be used to perturb the plasma with radio waves to solve problems in ionospheric, atmospheric, magnetospheric, and plasma physics. Among the ionospheric and plasma phenomena to be investigated using WISP instrumentation are VLF wave-particle interactions; ELF/VLF propagation; traveling ionospheric disturbances and gravity wave coupling; equatorial plasma bubble phenomena; plasma wave physics such as mode-coupling, dispersion, and instabilities; and plasma physics of the antenna-plasma interactions.

Fredricks, R. W.↗
NASA NTRS2012

The Magnetic Origins of Solar Activity

The defining physical property of the Sun's corona is that the magnetic field dominates the plasma. This property is the genesis for all solar activity ranging from quasi-steady coronal loops to the giant magnetic explosions observed as coronal mass ejections/eruptive flares. The coronal magnetic field is also the fundamental driver of all space weather; consequently, understanding the structure and dynamics of the field, especially its free energy, has long been a central objective in Heliophysics. The main obstacle to achieving this understanding has been the lack of accurate direct measurements of the coronal field. Most attempts to determine the magnetic free energy have relied on extrapolation of photospheric measurements, a notoriously unreliable procedure. In this presentation I will discuss what measurements of the coronal field would be most effective for understanding solar activity. Not surprisingly, the key process for driving solar activity is magnetic reconnection. I will discuss, therefore, how next-generation measurements of the coronal field will allow us to understand not only the origins of space weather, but also one of the most important fundamental processes in cosmic and laboratory plasmas.

Antiochos, S. K.↗
DOE OSTI2025

Search for Quasiparticle Scattering in the Quark-Gluon Plasma with Jet Splittings in 𝑝⁢𝑝 and Pb-Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration reports measurements of the large relative transverse momentum (𝑘 𝑇 ) component of jet substructure in 𝑝⁢𝑝 and Pb-Pb collisions at center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV. Enhancement in the yield of such large-𝑘 T emissions in head-on Pb-Pb collisions is predicted to arise from partonic scattering with quasiparticles of the quark-gluon plasma. The analysis utilizes charged-particle jets reconstructed by the anti-𝑘 T algorithm with resolution parameter 𝑅 = 0.2 in the transverse-momentum interval 60 < 𝑝 T,ch,jet < 80 GeV/𝑐. The soft drop and dynamical grooming algorithms are used to identify high transverse momentum splittings in the jet shower. Comparison of measurements in Pb-Pb and 𝑝⁢𝑝 collisions shows medium-induced narrowing, corresponding to yield suppression of high-𝑘 𝑇 splittings, in contrast to the expectation of yield enhancement due to quasiparticle scattering. The measurements are compared to theoretical model calculations incorporating jet modification due to jet-medium interactions (“jet quenching”), both with and without quasiparticle scattering effects. These measurements provide new insight into the underlying mechanisms and theoretical modeling of jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Particlization in fluid dynamical simulations of heavy-ion collisions: The IS3D module

The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamic Phase Alignment in Inertial Alfvén Turbulence

In weakly collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin depth, k ⊥ d e ≳ 1. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as cos α k ∝ $k$$^{–1}_{⊥}$. The kinetic and magnetic energy spectra scale as ∝ $k$$^{–5/3}_{⊥}$ and ∝ $k$$^{–11/3}_{⊥}$, respectively. As a result of the dual direct cascade, the generalized helicity spectrum scales as ∝ $k$$^{–5/3}_{⊥}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the k ⊥ d e >> 1 range. Turbulent eddies exhibit a phase-space anisotropy k ∥ ∝ $k$$^{5/3}_{⊥}$, consistent with critically balanced inertial Alfvén fluctuations. Furthermore, our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth’s magnetosheath and ionosphere, solar corona, and nonrelativistic pair plasmas, as well as to strongly rotating nonionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Solar-terrestrial physics.

The work is concerned with solar influences on the earth's upper atmosphere and electromagnetic environment partly through solar ultraviolet and X rays, but particularly through corpuscular radiations. The sun and interplanetary space, the internal structure and magnetic field of the earth, the terrestrial atmosphere and photochemistry, and the dynamics of the upper atmosphere and dynamo action are considered. The formation of the magnetosphere, energetic particles, plasma, and electromagnetic waves in the magnetosphere, solar storms and their extension into interplanetary space, and magnetospheric storms are treated.

Akasofu, S.-I.↗
DOE OSTI2025

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗