Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROSTATIC WAVE”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Growth and Saturation of the Electron Drift Instability in a Crossed Field Plasma

The linear growth and nonlinear energy transfer of the electron drift instability (EDI) are experimentally measured in the plume of a low-temperature, Hall effect discharge. A frequency-based bispectral analysis technique applied to fast ion density fluctuation measurements shows a growth rate function that is qualitatively similar to predictions from the linear instability dispersion relation, but an order of magnitude smaller. Calculation of the nonlinear transfer function indicates multiple three-wave interactions between high-frequency resonances of the instability in addition to an inverse energy cascade toward lower-frequency modes. Furthermore, these results are discussed in the context of recent theoretical, numerical, and experimental efforts on the EDI in Hall effect discharges and how the EDI may impact anomalous cross field transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Anomalous cross-field transport in a Hall thruster inferred from direct measurement of instability growth rates

The contribution of the electron drift instability to anomalous electron transport is experimentally assessed in a Hall effect discharge. The transport is represented by an anomalous collision frequency, which is related through quasilinear theory to the energy and growth rate of the instability. The wave energy is measured directly with ion saturation probes, while estimates of the growth rate are employed based on both linearized theory and direct measurement. The latter measurement is performed with a bispectral analysis method. The wave-driven collision frequency is compared to measurements of the actual collision frequency inferred from a method based on laser- induced fluorescence. It is found that estimates for transport using linearized theory for the growth differ by over an order of magnitude from the actual anomalous collision frequency in the plasma. The wave-driven anomalous collision frequency with measured growth, however, is shown to agree with the electron collision frequency in magnitude and capture aspects of the trends in spatial variation. This result demonstrates experimentally that wave-driven effects ultimately can explain the observed cross-field transport in these devices. As a result, the implications of this finding are discussed in the context of the key lengthscales that drive the transport as well as the implications identifying reduced fidelity models that could be used to predict anomalous collision frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

EMPIRE-PIC: A Performance Portable Unstructured Particle-in-Cell Code

In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
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 OSTI2026

Thomson Scattering with Gain

Thomson-scattering signals can be significantly modified by convective gains associated with the stimulated Raman scattering and stimulated Brillouin scattering instabilities as the scattered light traverses the probe beam. Gain results in amplification and narrowing of the red-shifted (i.e., Stokes) scattered light, whereas blue-shifted (anti-Stokes) features are depleted and broadened, making the relative prominence of these features a key signature of gain. Not accounting for instability growth has been a ubiquitous source of error in previous Thomson-scattering measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

jaxhps: An elliptic PDE solver built with machine learning in mind

Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. 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 to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Magnetogenesis by Wave-driven Momentum Exchange

When multiple species interact with an electrostatic ion acoustic wave, they can exchange momentum, despite the lack of momentum in the field itself. Furthermore, the resulting force on the electrons can have a curl, and thus give rise to compensating electric fields with curl on magnetohydrodynamic timescales. As a result, a magnetic field can be generated. Surprisingly, in some astrophysical settings, this mechanism can seed magnetic fields with growth rates even larger than through the traditional Biermann battery.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Localized Excitation of Electromagnetic Ion Cyclotron Waves From Anisotropic Protons Filtered by Magnetic Dips

The excitation of electrostatic and/or electromagnetic waves in the plasma universe is often associated with anisotropic velocity distributions of charged particles. In Earth's inner magnetosphere, this anisotropy can gradually develop as particles injected from the magnetotail drift around the Earth at different speeds depending on their energy and pitch angle. Here, we show that the perpendicular-moving and bouncing ions can be separated more abruptly near the injection front. These pitch-angle filters are localized magnetic dip structures formed by the diamagnetic behavior of the injected particles, which can trap perpendicular-moving ions and allow bouncing ions to overtake. The resulting ion anisotropy facilitates the rapid generation of electromagnetic ion cyclotron (EMIC) waves, which in turn can largely reshape the Van Allen radiation belts. This scenario is examined by case and statistical observations, together with numerical simulations that reproduce most of the observational signatures, to support the causal relationship between magnetic dips, anisotropic ion distributions, and localized excitation of EMIC waves. Finally, our study highlights the important roles of magnetic dips in the inner magnetospheric dynamics, as pitch-angle filters of the injected ions and traveling hotspots of EMIC wave activities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Whistler Waves Associated With Electron Beams in Magnetopause Reconnection Diffusion Regions

Abstract Whistler waves are often observed in magnetopause reconnection associated with electron beams. We analyze seven MMS crossings surrounding the electron diffusion region (EDR) to study the role of electron beams in whistler excitation. Waves have two major types: (a) Narrow‐band waves with high ellipticities and (b) broad‐band waves that are more electrostatic with significant variations in ellipticities and wave normal angles. While both types of waves are associated with electron beams, the key difference is the anisotropy of the background population, with perpendicular and parallel anisotropies, respectively. The linear instability analysis suggests that the first type of wave is mainly due to the background anisotropy, with the beam contributing additional cyclotron resonance to enhance the wave growth. The second type of broadband waves are excited via Landau resonance, and as seen in one event, the beam anisotropy induces an additional cyclotron mode. The results are supported by particle‐in‐cell simulations. We infer that the first type occurs downstream of the central EDR, where background electrons experience Betatron acceleration to form the perpendicular anisotropy; the second type occurs in the central EDR of guide field reconnection. A parametric study is conducted with linear instability analysis. A beam anisotropy alone of above ∼3 likely excites the cyclotron mode waves. Large beam drifts cause Doppler shifts and may lead to left‐hand polarizations in the ion frame. Future studies are needed to determine whether the observation covers a broader parameter regime and to understand the competition between whistler and other instabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Periodic Narrowband Radio Wave Emissions and Inward Plasma Transport at Saturn's Magnetosphere

The abrupt brightening of an Energetic Neutral Atom (ENA) blob or cloud has been interpreted as plasma injection in the Kronian magnetosphere (termed ENA injection herein). Morphologically, there appears to be two types of abrupt ENA cloud brightening: (1) the brightening of a large cloud usually seen at distances >10–12 R {sub s} (R {sub s} ∼ 60,268 km) in the midnight or postmidnight region; (2) the brightening of a smaller cloud usually seen at distances <10–12 R {sub s} around 21–03 magnetic local time. Among many radio waves observed at Saturn, type 2 ENA injections correlate best with the 5 kHz narrowband (NB) waves. Using Cassini Ion and Neutral Camera (INCA) and Radio and Plasma Wave Science (RPWS) data, we examine the periodicities of the type 2 ENA injections and the 5 kHz NB emissions as well as their cross-correlations, which have been previously used to measure the lag times or phase differences. Because correlational analysis can only establish linear relationships, we also use mutual information to establish linear and nonlinear relationships. On average, the peak of the 5 kHz NB emission lags those of the type 2 ENA injection by about a few minutes to 2 hr. The injection of hot plasma to the inner magnetosphere can lead to temperature anisotropy, which can generate electrostatic upper hybrid waves, which upon encountering the high-density gradient at the outer edge of the Enceladus density torus, can mode convert to Z mode and then to O mode. The 5 kHz NB waves commonly propagate in the O mode.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Density Irregularities on Radio Frequency Wave Propagation in Ionospheric Plasmas

Density irregularities play a vital role in determining how radio frequency (RF) waves travel through plasmas. In the Earth’s ionosphere, these density irregularities also impact radio communication. In this study, we conduct a detailed numerical analysis of RF wave propagation in small-scale ionospheric density irregularities using the advanced Petra-M code. We focus specifically on high-frequency (HF) waves, ranging from 3 to 30 MHz, which are essential for military, amateur radio operators, and emergency communications. By introducing density structures, such as equatorial plasma bubbles derived from fluid simulations, we demonstrate that HF waves can scatter in multiple directions when they encounter these irregularities. Additionally, we observe significant mode conversion, where incoming electromagnetic waves transform into electrostatic modes within the density gradient layer. This shows that smaller density irregularities can greatly weaken signals or cause complete signal loss for receivers, emphasizing the need for increased awareness and innovative solutions in radio communication transmission.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗