Engineering Papers⌕ Search

Engineering topics

Roytershteyn, Vadim

Publications and source records attributed to Roytershteyn, Vadim.

10 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Oblique instability of quasi-parallel whistler waves in the presence of cold and warm electron populations

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effective viscosity, resistivity, and Reynolds number in weakly collisional plasma turbulence

ABSTRACT We examine dissipation and energy conversion in weakly collisional plasma turbulence, employing in situ observations from the Magnetospheric Multiscale mission and kinetic particle-in-cell simulations of proton–electron plasma. A previous result indicated the presence of viscous-like and resistive-like scaling of average energy conversion rates – analogous to scalings characteristic of collisional systems. This allows for extraction of collisional-like coefficients of effective viscosity and resistivity, and thus also determination of effective Reynolds numbers based on these coefficients. The effective Reynolds number, as a measure of the available bandwidth for turbulence to populate various scales, links turbulence macroscale properties with kinetic plasma properties in a novel way.

Astronomy & Astrophysics↗
DOE OSTI2023

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Using Kinetic Entropy to Understand Dissipation in Re-connection and Turbulence

This report summarizes results obtained at Space Science Institute (SSI) under a collaborative project funded by the NSF-DOE partnership in Basic Plasma Physics. The overall goal of the project was to perform a study of kinetic entropy as a diagnostic of dissipation in several archetypal models, such as magnetic reconnection and decaying plasma turbulence. The specific task of the SSI group was to perform collisionless and weakly collisional Particle-in-Cell (PIC) simulations modeling these systems. Overall, the project has demonstrated that, as expected, kinetic entropy carries an imprint of weakly collisional processes commonly thought to be associated with dissipation and thus provides an important diagnostic of dissipation. Specific entropy-related metrics have been proposed that could be used as quantitative measures of dissipation. These metrics have been implemented as diagnostics in PIC simulations and their utility, as well as relation to other, previously proposed, measures of dissipation have been demonstrated. Further, PIC simulations were used to investigate a general question of how energy is dissipated in magnetized weakly collisional turbulence in both non-relativistic and relativistic plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Scaling of Electron Heating by Magnetization During Reconnection and Applications to Dipolarization Fronts and Super-Hot Solar Flares

Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields (Shuster et al., 2014, https://doi. org/10.1002/2014GL060608). Here, we present a theory of the dependence of the major and minor radii of the ring distributions solely in terms of upstream (lobe) plasma conditions, thereby allowing a prediction of the associated temperature and temperature anisotropy of the rings in terms of upstream parameters. We test the validity of the prediction using 2.5-dimensional particle-in-cell (PIC) simulations with varying upstream plasma density and temperature, finding excellent agreement between the predicted and simulated values. We confirm the Shuster et al. suggestion for the cause of the ring distributions, and also find that the ring distributions are located in a region marked by a plateau, or shoulder, in the reconnected magnetic field profile. The predictions of the temperature are consistent with observed electron temperatures in dipolarization fronts, and may provide an explanation for the generation of plasma with temperatures in the 10s of MK in super-hot solar flares. A possible extension of the model to dayside reconnection is discussed. Since ring distributions are known to excite whistler waves, the present results should be useful for quantifying the generation of whistler waves in reconnection exhausts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Plasma Dynamics in Low-Electron-Beta Environments

Recent in situ measurements by the MMS and Parker Solar Probe missions bring interest to small-scale plasma dynamics (waves, turbulence, magnetic reconnection) in regions where the electron thermal energy is smaller than the magnetic one. Examples of such regions are the Earth’s magnetosheath and the vicinity of the solar corona, and they are also encountered in other astrophysical systems. In this brief review, we consider simple physical models describing plasma dynamics in such low-electron-beta regimes, discuss their conservation laws and their limits of applicability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nonlinear coupling of whistler waves to oblique electrostatic turbulence enabled by cold plasma

Kinetic simulations and theory demonstrate that whistler waves can excite oblique, short-wavelength fluctuations through secondary drift instabilities if a population of sufficiently cold plasma is present. The excited modes lead to heating of the cold populations and damping of the primary whistler waves. The instability threshold depends on the density and temperature of the cold population and can be relatively small if the temperature of the cold population is sufficiently low. Furthermore, this mechanism may thus play a significant role in controlling amplitude of whistlers in the regions of the Earth's magnetosphere where cold background plasma of sufficient density is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗