Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic simulation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection (RR): ▪ Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in RR. ▪ The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research. ▪ RR plays a key role in global models of high-energy sources, in terms of both global-scale layers and reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools.

acceleration of particles↗
DOE OSTI2024

Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Effects of Gasoline Composition and Operating Parameters on the Response of End-Gas Autoignition to Nitric Oxide in a Lean-Burn DI-SI Engine

Lean operation of spark-ignition engines can lead to engine thermal efficiency gains and lower NOx emissions due to reduced combustion temperatures. Yet, lean operation could still face challenges in end-gas autoignition and knock generation due to higher intake pressures and trapped NO in the residual gas. Here, this study evaluates the impact of NO on end-gas autoignition for two gasoline fuels with similar octane rating but different composition: high cycloalkane fuel (HCA) and high olefin fuel (HO). Experiments were performed at stoichiometric and lean (λ = 2) conditions and at two engine speeds of 1400 rpm and 2000 rpm. Accompanying chemical kinetics simulations in CHEMKIN revealed that the mechanisms controlling the effect of NO on autoignition are similar λ = 2 and λ = 1, with NO + HO 2 = NO 2 + OH being the main pathway for enhancing reactivity by promoting low-temperature heat release (LTHR). The compositionally different fuels reacted differently to NO seeding and engine speed, and differences were augmented at λ = 2 compared to λ = 1 as the end-gas autoignition shifted to the low temperature regime. HO, which has inherent low temperature chemistry, was strongly impacted by engine speed at low NO seeding levels, with no noticeable peak of LTHR detected at 2000 rpm. On the other hand, LTHR of HCA was marginally affected by shortened residence time at higher engine speed as NO + HO 2 reaction was not greatly affected by shorter time scales, since HO 2 production was sustained even at 2000 rpm to support OH generation from NO + HO 2 . Contrary to HO, HCA exhibited greater sensitivity to NO seeding, as the increased OH production at higher NO concentrations offset the OH-quenching effect of cyclopentane, which accounts for 28.6% of HCA’s composition. Consequently, a sensitivity analysis revealed that fuels with weak inherent low-temperature chemistry, like HCA, are likely to be more sensitive to variations in NO concentration and charge temperature, whereas fuels with strong low temperature chemistry are more sensitive to variations in end-gas λ and intake pressure.

Energy↗
NASA NTRS1986

Nonadiabatic electron heating at high-Mach-number perpendicular shocks

Fully kinetic simulations of high-Mach-number (HMN) perpendicular collisionless shocks are described. It is shown that electron acceleration in the cross-shock electron field can produce downstream electron temperature significantly higher than those expected for adiabatic compression. The momentum space for test electrons at Mach 6 is illustrated.

Tokar, R. L.↗
NASA NTRS1986

Interaction of two collisionless shocks

Kinetic simulations of the interaction between two collisionless shocks are presented. During the collision of two perpendicular shocks, the shock electromagnetic field structures pass through each other, while the previously shocked ions are kept separate by the electric field arising in the collision. When two supercritical shocks collide, a fraction of ions are accelerated up to an order of magnitude in energy by first being reflected at one shock, then interacting with the electric fields of the other shock.

Cargill, P. J.↗
NASA NTRS1992

Effects of equatorially trapped ions on refilling of the plasmasphere

The kinetic aspects of shock formation in response to the equatorial heating of ions are reported by means of small-scale kinetic simulations of countersteaming plasma flows along a spatially varying magnetic field having the features of the dipolar magnetic field lines. Attention is given to the trapping of the field-aligned flowing ions in response to the simulated equatorial heating of ions and the associated self-consistent structures in the electric potentials and the f-plasma flow pattern. Trapping is found to lead to the formation of electrostatic shocks, which propagate away from the 'equator' as seen in the hydrodynamic model of Singh and Torr (1990). The shocks form near the mirror points of the 'equatorially' heated ions. The second stage of the refilling with equatorially trapped ions is reported. A variety of structures in density and temperature distributions as the refilling proceeds slowly is found.

Singh, Nagendra↗
NASA NTRS1992

The formation of the wall region - Consequences in the near earth magnetotail

This paper discusses important new findings obtained from global kinetic simulations of magnetotail plasma. A region of strongly nonadiabatic ion acceleration (known as the wall region) exists in the near earth tail and demarcates two very different regimes of ion motion: adiabatic and quasi-adiabatic. A strong enhancement of the cross-tail current occurs on the tailward side of the wall. A comparison of numerical and adiabatic pressure profiles indicates that nonadiabatic processes operating in this region may contribute significantly to a pressure balance relief in the course of quasi-steady magnetospheric convection.

Ashour-Abdalla, Maha↗
NASA NTRS1993

Plasma irregularities caused by cycloid bunching of the CRRES G-2 barium release

The Combined Release and Radiation Effects Satellite (CRRES) spacecraft carried a number of barium thermite canisters for release into the upper atmosphere. The barium release labeled G-2 showed evidence of curved irregularities not aligned with the ambient magnetic field B. The newly discovered curved structures can be explained by a process called cycloid bunching. Cycloid bunching occurs when plasma is created by photoionization of a neutral cloud injected at high velocity perpendicular to B. If the injection velocity is much larger than the expansion speed of the cloud, the ion trail will form a cycloid that has irregularities spaced by the product of the perpendicular injection speed and the ion gyroperiod, Images of the solar-illuminated barium ions are compared with the results of a three-dimensional kinetic simulation. Cycloid bunching is shown to be responsible for the rapid generation of both curved and field-aligned irregularities in the CRRES G-2 experiment.

Bernhardt, P. A.↗
NASA NTRS1993

Comment on 'Comparison of observed and calculated implanted ion distributions outside Comet Halley's bow shock' by T. I. Gombosi, M. Neugebauer, A. D. Johnstone, A. J. Coates, and D. E. Huddleston

It is argued that the formulation used by Gombosi (1991) does not distinguish between two fundamental types of resonant interactions (those between cometary pickup ions and LF electromagnetic waves). Two problems arise from using their overly simplified expressions. The wave polarization of the solar wind Alfven wave is not purely right-handed, and values from alpha sub com and alpha sub sw are not well known, nor are the 'polarizations' of the 'backward'-going waves in either case. Gombosi et al. discuss the antisunward vs the sunward propagating solar wind Alfven wave, and assume a reasonable value of 0.2 for alpha sub sw. However, alpha sub com is totally unknown. It is concluded that a more correct analysis of wave polarization can lead to significant changes in the formulation and results from Gombosi et al. In his reply Gombosi argues that his simplified wave-particle interaction model is consistent with Giotto observations and plasma kinetic simulations and it leads to a reasonable approximation of the implanted ion transport coefficients. The more complicated approach suggested by Tsurutani & Thorne (1993) requires the introduction of several ad hoc parameters and it does not lead to any change of the energy diffusion coefficient.

Tsurutani, Bruce T.↗
NASA NTRS1993

Nonintrusive, multipoint velocity measurements in high-pressure combustion flows

A combined experimental and analytical effort was conducted to demonstrate the applicability of OH Doppler-shifted fluorescence imaging of velocity distributions in supersonic combustion gases. The experiments were conducted in the underexpanded exhaust flow from a 6.8 atm, 2400 K, H2-O2-N2 burner exhausting into the atmosphere. In order to quantify the effects of in-plane variations of the gas thermodynamic properties on the measurement accuracy, a set of detailed measurements of the OH (1,0) band collisional broadening and shifting in H2-air gases was produced. The effect of pulse-to-pulse variations in the dye laser bandshape was also examined in detail and a modification was developed which increased in the single pulse bandwidth, thereby increasing the intraimage velocity dynamic range as well as reducing the sensitivity of the velocity measurement to the gas property variations. Single point and imaging measurements of the velocity field in the exhaust flowfield were compared with 2D, finite-rate kinetics simulations of the flowfield. Relative velocity accuracies of +/- 50 m/s out of 1600 m/s were achieved in time-averaged imaging measurements of the flow over an order of magnitude variation in pressure and a factor of two variation in temperature.

Allen, M.↗
NASA NTRS1992

Effects of near-Earth stochastic acceleration and reflections of magnetotail ions on the formation of auroral arcs

Findings obtained from global kinetic simulations of magnetotail plasma are discussed. A region of strongly nonadiabatic ion acceleration (known as the 'wall' region) exists in the near earth tail and demarcates two very different regimes of ion motion: adiabatic and quasi-adiabatic. After convection through the wall, ion distributions rapidly become isotropized and thermalized. A strong enhancement of the cross tail current occurs on the tailward side of the wall. Comparison of numerical and adiabatic pressure profiles indicates that nonadiabatic processes operating in this region may contribute significantly to a pressure balance relief in the course of quasi-steady magnetospheric convection.

Ashour-Abdalla, Maha↗
NASA NTRS1994

On the cross-field diffusion of ions in one- and two-dimensional hybrid simulations of collisionless shocks

It can be demonstrated analytically that under certain geometries used in numerical simulations of collisionless shocks in which there is at least one ignorable spatial coordinate, the transport of particles across the magnetic field is essentially zero. This notion is tested using one- and two-dimensional hybrid simulations (kinetic ions/fluid electrons). We find, as the theorem predicts, the particles treated kinetically are tied to the same field line on which they start.

Giacalone, Joe↗
NASA NTRS1993

The interaction with the lower ionosphere of electromagnetic pulses from lightning: Excitation of optical emissions

A self consistent and fully kinetic simulation of the interaction of lightning radiated electromagnetic (EM) pulses with the nighttime lower ionosphere indicates that optical emissions observable with conventional instruments would be excited. For example, emissions of the 1st and 2nd positive bands of N2 occur at rates reaching 7 x 10(exp 7) and 10(exp 7) cu cm/s respectively at 92 km altitude for a lightning discharge with an electric field E(sub 100) = 20 V/m (normalized to a 100 km distance). The maximum height integrated intensities of these emissions are 4 x 10(exp 7) and 6 x 10(exp 6) R respectively, lasting for approx. 50 micrometers.

Taranenko, Y. N.↗
NASA NTRS1994

Observations of nonadiabatic acceleration of ions in Earth's magnetotail

We present observations of the three-dimensional velocity distributions of protons in the energy range 20 eV to 52 keV at locations within and near the current sheet of Earth's magnetotail at geocentric radial distances 35 to 87 R(sub E). These measurements were acquired on December 8, 1990, with a set of electrostatic analyzers on board the Galileo spacecraft during its approach to Earth in order to obtain one of its gravitational assists to Jupiter. It is found that the velocity distributions are inadequately described as quasi-Maxwellian distributions such as those found in the central plasma sheet at positions nearer to Earth. Instead the proton velocity distributions can be categorized into two major types. The first type is the 'lima bean' shaped distribution with high-speed bulk flows and high temperatures that are similar to those found nearer to Earth in the plasma sheet boundary layer. The second type consists of colder protons with considerably lesser bulk flow speeds. Examples of velocity distributions are given for the plasma mantle, a region near the magnetic neutral line, positions earthward and tailward of the neutral line, and the plasma sheet boundary layer. At positions near the neutral line, only complex velocity distributions consisting of the colder protons are found, whereas both of the above types of distributions are found in and near the current sheet at earthward and tailward locations. Bulk flows are directed generally earthward and tailward at positions earthward and tailward of the neutral line, respectively. Only the high-speed, hot distribution is present in the plasma sheet boundary layer. The observations are interpreted in terms of the nonadiabatic acceleration of protons that flow into the current sheet from the plasma mantle. For this interpretation the hot, 'lima bean' shaped distributions are associated with meandering, or Speiser, orbits in the current sheet. It is suggested that the colder, lower-speed proton velocity distributions are the result of fractional or few gyromotions before ejection out of the current sheet, but this speculation must be further investigated with appropriate kinetic simulation of trajectories.

Frank, L. A.↗
NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗
NASA NTRS1995

Heating and scattering of ring-beam distributions by turbulence

Pickup ions in the solar wind are initially are born in ring-beam distributions, i.e. f(v) varies as delta(v(sub perpendicular) - V(sub sw)sin(Theta)) delta(v(sub parallel) - V(sub sw)cos(Theta)), where Theta is the angle between the solar wind velocity and the IMF(Interplanetary Magnetic Field), and V(sub sw) is the solar wind speed. Often the distribution has been presumed to relax to a distribution that is isotropic in Theta and essentially mono-energetic, a shell or a 'bi-spherical distribution.' However solar wind turbulence is capable of heating the ring distribution on the timescale of a few tens of gyroperiods, a timescale not greatly distinct from that required for pitch angle scattering to a shell. To describe this effect, we have performed test-particle studies of the heating/scattering of the ring beam distribution by MHD turbulence, adopting various models for the MHD fluctuations, including slab and fully dynamic 2D and 3D incompressible turbulence. Furthermore, a system composed of a cold ion ring and a background plasma is unstable to several kinetic plasma instabilities. We carried out kinetic simulations of the ring beam distribution, showing that plasma instabilities also rapidly energize and scatter particles. Results will be presented comparing relaxation and heating rates of the ring-beam distribution by the various mechanisms.

Gray, P. C.↗
NASA NTRS1998

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic-simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost supercomputers, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enables us to make a number of important discoveries predictions, which need to be looked for in satellite data. Examples include the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH), (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs), and (3) this finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗
NASA NTRS1998

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost desktop superconductors, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enabled us to make a number of important discoveries/predictions, which need to be looked for in satellite data. Examples include: (1) the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH;) (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs;) and (3) the finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗