SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

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 NTRS1983

Energy transfer in the quiet and disturbed magnetosphere

Progress on magnetospheric energy transfer reported during the 1979-1982 quadrennium is reviewed. A discussion of the macroscopic control of magnetospheric phenomena by solar-wind electric and magnetic fields and kinetic pressure is followed by sections on energy transfer processes at the dayside magnetopause, in the magnetotail, near geosynchronous orbit, and at altitudes of about 1 earth radius above the auroral oval. Rapid progress in these areas has been stimulated primarily by ISEE and S3-3 observations and by the development of techniques for the large-scale computer simulation of plasma phenomena.

Burch, J. L.↗

NASA NTRS1992

Analysis of cyclotron harmonic emissions at the outer planets

The flyby missions of Voyagers 1 and 2 at Jupiter, Saturn and Uranus revealed intense waves above the electron gyrofrequency. Observation of waves at the upper hybrid frequency is often accompanied by power at adjacent electron Bernstein harmonics, and the relative power in these modes depends both on the density and temperature ratios of the cold background electron population and the hot magnetospheric electrons which drive the instability. A model of electron distributions which is consistent with observations is used for analysis of the excited waves, their dependence upon plasma parameters, and the time scales of the saturation processes. It is shown that in the presence of two-temperature electron distributions the linear excitation is due to a fluidlike coupling of two eigenmodes for perpendicular propagation and to kinetic destalization of oblique modes. The dependence of linear growth rates on propagation angle is presented, along with results from particle simulations. A quasi-linear diffusion time for relaxing the hot electron loss cone is calculated and compared with simulation results. This time scale is faster than for local saturation by heating the cold population, and also the convective amplification time scale, suggesting that the waves saturate at quasi-linear levels, while being convectively localized to the equatorial regions of the outer planetary magnetospheres.

Roth, I.↗

NASA NTRS2006

Electrodynamics in a Very Thin Current Sheet Leading to Magnetic Reconnection

We study the formation of a very thin current sheet (CS) and associated plasma electrodynamics using three-dimensional (3-D) particle-in-cell simulations with ion to electron mass ratio M/m=1836. The CS is driven by imposed anti-parallel magnetic fields. The noteworthy features of the temporal evolution of the CS are the following: (i) Steepening of the magnetic field profile B,(z) in the central part of the CS, (ii) Generation of three-peak current distribution with the largest peak in the CS center as B,(z) steepens, (iii) Generation of converging electric fields forming a potential well in the CS center in which ions are accelerated. (iv) Electron and ion heating in the central part of the CS by current-driven instabilities (CDI). (v) Re-broadening of the CS due to increased kinetic plasma pressure in the CS center. (vi) Generation of electron temperature anisotropy with temperature perpendicular to the magnetic field being larger than the parallel one. (vii) Current disruption by electron trapping in an explosively growing electrostatic instability (EGEI) and electron tearing instability (ETI). (viii)The onset of EGEI coincides with an increase in the electron temperature above the temperature of the initially hot ions as well as the appearance of new shear in the electron drift velocity. (ix) Bifurcation of the central CS by the current disruption. (x) Magnetic reconnection (MR) beginning near the null in B, and spreading outward. (xi) Generation of highly energized electrons reaching relativistic speeds and having isotropic pitch-angle distribution in the region of reconnected magnetic fields. We compare some of these features of the current sheet with results from laboratory and space experiments.

Singh, Nagendra↗

NASA NTRS1995

Wave and ion evolution downstream of quasi-perpendicular bow shocks

Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Mckean, M. E.↗

NASA NTRS2016

MMS Observations of Large Guide Field Symmetric Reconnection Between Colliding Reconnection Jets at the Center of a Magnetic Flux Rope at the Magnetopause

We report evidence for reconnection between colliding reconnection jets in a compressed current sheet at the center of a magnetic flux rope at Earth's magnetopause. The reconnection involved nearly symmetric Inflow boundary conditions with a strong guide field of two. The thin (2.5 ion-skin depth (d(sub i) width) current sheet (at approximately 12 d(sub i) downstream of the X line) was well resolved by MMS, which revealed large asymmetries in plasma and field structures in the exhaust. Ion perpendicular heating, electron parallel heating, and density compression occurred on one side of the exhaust, while ion parallel heating and density depression were shifted to the other side. The normal electric field and double out-of-plane (bifurcated) currents spanned almost the entire exhaust. These observations are in good agreement with a kinetic simulation for similar boundary conditions, demonstrating in new detail that the structure of large guide field symmetric reconnection is distinctly different from antiparallel reconnection.

Oieroset, M.↗

NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

NASA NTRS1994

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗

NASA NTRS2021

Magnetospheric Multiscale Observations of an Expanding Oxygen Wave in Magnetic Reconnection

Heavier plasma species such as oxygen ions can have a large impact on the magnetic reconnection process. It has been hypothesized that the acceleration of demagnetized oxygen ions by the Hall electric field will lead to the formation of an oxygen wave that expands into the exhaust. By comparing data from NASA's Magnetospheric Multiscale mission to a fully kinetic particle-in-cell simulation, we can for the first time provide observational evidence of such an expanding oxygen wave. The wave is characterized by an oxygen jet consisting of cold ions directed toward the neutral sheet associated with a density cavity. This density cavity forms as the EO are subject to collective acceleration by the Hall electric field leaving behind a region of low-density oxygen ions. Our results are important for the understanding of the role and effect of oxygen ions in magnetic reconnection.

Hakon Midthun Kolsto↗

NASA NTRS2011

Effect of Inductive Coil Geometry and Current Sheet Trajectory of a Conical Theta Pinch Pulsed Inductive Plasma Accelerator

Results are presented demonstrating the e ect of inductive coil geometry and current sheet trajectory on the exhaust velocity of propellant in conical theta pinch pulsed induc- tive plasma accelerators. The electromagnetic coupling between the inductive coil of the accelerator and a plasma current sheet is simulated, substituting a conical copper frustum for the plasma. The variation of system inductance as a function of plasma position is obtained by displacing the simulated current sheet from the coil while measuring the total inductance of the coil. Four coils of differing geometries were employed, and the total inductance of each coil was measured as a function of the axial displacement of two sep- arate copper frusta both having the same cone angle and length as the coil but with one compressed to a smaller size relative to the coil. The measured relationship between total coil inductance and current sheet position closes a dynamical circuit model that is used to calculate the resulting current sheet velocity for various coil and current sheet con gura- tions. The results of this model, which neglects the pinching contribution to thrust, radial propellant con nement, and plume divergence, indicate that in a conical theta pinch ge- ometry current sheet pinching is detrimental to thruster performance, reducing the kinetic energy of the exhausting propellant by up to 50% (at the upper bound for the parameter range of the study). The decrease in exhaust velocity was larger for coils and simulated current sheets of smaller half cone angles. An upper bound for the pinching contribution to thrust is estimated for typical operating parameters. Measurements of coil inductance for three di erent current sheet pinching conditions are used to estimate the magnetic pressure as a function of current sheet radial compression. The gas-dynamic contribution to axial acceleration is also estimated and shown to not compensate for the decrease in axial electromagnetic acceleration that accompanies the radial compression of the plasma in conical theta pinches.

Hallock, Ashley K.↗

NASA NTRS1989

Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.

Max, C. E.↗

NASA NTRS1993

Solar wind-magnetosphere interaction as simulated by a 3D, EM particle code

The results of simulating the solar wind-magnetosphere interaction with a three dimensional, electromagnetic (EM) particle code are presented. Hitherto such global simulations were done with magnetohydrodynamic (MHD) codes while lower dimensional particle or hybrid codes served to account for microscopic processes and such transport parameters as have to be introduced ad hoc in MHD. The kinetic model combines macroscopic and microscopic tasks. It relies only on the Maxwell curl equations and the Lorentz equation for particles. The preliminary results are for an unmagnetized solar wind plasma streaming past a dipolar magnetic field. The results show the formation of a bow shock and a magnetotail, the penetration of energetic particles into cusp and radiation belt regions, and dawn to dusk asymmetries.

Buneman, O.↗

NASA NTRS2002

Prebiotic Processing induced by Comet and Meteor Impact

In their study of organic synthesis from impact shocks using the laser-induced-plasma (LIP) technique, McKay and Borucki(l) found that organic synthesis preferentially occurred in a reducing gas mixture rich in methane, and not in a mixture rich in carbon dioxide. This result means chemical models based on the thermodynamical equilibrium approach do not apply to shock chemistry. In this study, we employ the technique of reacting flow, i.e., chemical kinetics in a fluid flow, to simulate the chemistry occurring in LIP and in the wake region from comet or meteor impact. Three different air compositions have been used: (1) 1/3 CO2 and 2/3 H2, (2) pure CH4, and (3) 1/4 CH4, 1/4 CO2, and 1/2 H2O. The stoichiometric ratio of gas mixtures (1) and (3) are kept the same. For (1) we obtain equal mole fractions of CO and H2O as the major products and for (2) C2H2 is the major product. In both cases our results are in agreement with Ref. (1). For (3) we find an interesting case where the nature of chemicals produced to be critically dependent on the flowfield temperature. At the higher temperature part of the wake region, CO and H2O are the dominant products, whereas in the cooler region C2H2 is the dominant product. Further studies of these reactions, as well as for the gas mixture including N2, are being pursued.

Dateo, Christopher E.↗

NASA NTRS1998

A Study of Transport in the Near-Earth Plasma Sheet During A Substorm Using Time-Dependent Large Scale Kinetics

In this study we investigate the transport of H+ ions that made up the complex ion distribution function observed by the Geotail spacecraft at 0740 UT on November 24, 1996. This ion distribution function, observed by Geotail at approximately 20 R(sub E) downtail, was used to initialize a time-dependent large-scale kinetic (LSK) calculation of the trajectories of 75,000 ions forward in time. Time-dependent magnetic and electric fields were obtained from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interaction with the solar wind and the interplanetary magnetic field (IMF) as observed during the interval of the observation of the distribution function. Our calculations indicate that the particles observed by Geotail were scattered across the equatorial plane by the multiple interactions with the current sheet and then convected sunward. They were energized by the dawn-dusk electric field during their transport from Geotail location and ultimately were lost at the ionospheric boundary or into the magnetopause.

El-Alaoui, M.↗

NASA NTRS1983

Monte Carlo simulation of steady state shock structure including cosmic ray mediation and particle escape

Both hydrodynamic calculations (Drury and Volk, 1981, and Axford et al., 1982) and kinetic simulations imply the existence of thermal subshocks in high-Mach-number cosmic-ray-mediated shocks. The injection efficiency of particles from the thermal background into the diffusive shock-acceleration process is determined in part by the sharpness and compression ratio of these subshocks. Results are reported for a Monte Carlo simulation that includes both the back reaction of accelerated particles on the inflowing plasma, producing a smoothing of the shock transition, and the free escape of particles allowing arbitrarily large overall compression ratios in high-Mach-number steady-state shocks. Energy spectra and estimates of the proportion of thermal ions accelerated to high energy are obtained.

Ellison, D. C.↗

NASA NTRS2016

Global Effects of Transmitted Shock Wave Propagation Through the Earth's Inner Magnetosphere: First Results from 3-D Hybrid Kinetic Modeling

We use a new hybrid kinetic model to simulate the response of ring current, outer radiation belt, and plasmaspheric particle populations to impulsive interplanetary shocks. Since particle distributions attending the interplanetary shock waves and in the ring current and radiation belts are non-Maxwellian, waveparticle interactions play a crucial role in energy transport within the inner magnetosphere. Finite gyroradius effects become important in mass loading the shock waves with the background plasma in the presence of higher energy ring current and radiation belt ions and electrons. Initial results show that shocks cause strong deformations in the global structure of the ring current, radiation belt, and plasmasphere. The ion velocity distribution functions at the shock front, in the ring current, and in the radiation belt help us determine energy transport through the Earth's inner magnetosphere.

Lipatov, A. S.↗

NASA NTRS2014

Bidirectional Energy Cascades and the Origin of Kinetic Alfvenic and Whistler Turbulence in the Solar Wind

The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore how the free energy released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfvénic and whistler turbulence are excited that evolve through inverse and forward magnetic energy cascades.

plasmas↗

NASA NTRS1981

Effects of auroral-particle anisotropies and mirror forces on high-latitude electric fields

It is noted that, for most of the mechanisms for the strong electric fields that characterize the narrow regions in which there is acceleration and precipitation of ring current and/or plasma-sheet plasma, certain effects must be taken into account in simulations of auroral electric fields. The effects are those of auroral particle anisotropy, of mirror forces due to the inhomogeneous geomagnetic field, of auroral electron backscatter by the atmosphere, and of electron trapping by the combination of magnetic mirroring and electrostatic forces. What is more, the effects of the very strong perpendicular electric field must also be taken into account in a kinetic description of the Poisson equation in order to achieve a unified theory of the auroral electrostatic structure. Progress in these areas during the past few years is reviewed. It is shown that particle anisotropies and mirror forces can account for some basic electrostatic features of the quiet arc, while additional effects may be occurring in strong events in which the parallel potential drop is more than about 10 kV.

Chiu, Y. T.↗