Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle simulations”

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 163 records · Page 9

Particle-in-cell simulations of stochastic electron acceleration

The results of a series of particle-in-cell simulations of stochastic wave-particle interaction are presented. The threshold for stochasticity was confirmed. The simulations demonstrate that in a strong magnetic field plasma waves with quiver velocities much less than the speed of light but above a certain threshold can stochastically accelerate electrons to energies far greater than 1 MeV. Moreover, self-consistency effects drive return currents and produce energetic runaway electrons that violate an invariant of motion.

Akimoto, K.↗

Global hybrid simulation of unmagnetized planets - Comparison of Venus and Mars

Results from three-dimensional hybrid particle simulations of the solar wind interaction with the planets Mars and Venus are presented. The simulations produce shocks and magnetic barriers which are asymmetric. These results are qualitatively in agreement with data. In the absence of an ionosphere the subsolar shock standoff distance was found to agree with the observations if the Hall current is limited. It was also found that the solar wind interaction with Mars and Venus was substantially different. The interaction with Venus can be generally viewed as a magnetized interaction. The Mars interaction is very kinetic in nature and appears not to have a shock in the classic sense.

Brecht, Stephen H.↗

Simulation of magnetic field line stochasticity at the magnetopause

We have conducted a three-dimensional particle simulation to study the magnetic field line stochasticity at the magnetopause current layer. Our results show that the magnetic field lines become stochastic due to the overlap of the destabilized multiple tearing mode islands, which agrees with the percolation model suggested by Galeev et al. (1986). After the field lines become stochastic, these tearing modes grow even 2-3 times faster than in the linear stage and saturate at an amplitude level 3-4 times bigger than the single tearing mode without mode-mode coupling. The field line stochasticity also causes a strong particle diffusion across the current layer. The diffusion coefficient reaches to 10(exp 9) sq m/s for typical magnetopause parameters. Associated with the particle diffusion, the current layer becomes broader in width. As a result, the magnetic energy is dissipated into particle energy by heating parallel to the local magnetic field. The particle energy increases by 60%, while the magnetic helicity, which has always been regarded as a good invariant, changes by 20%.

Wang, Zhi↗

A simulation study of cold electron heating by loss cone instabilities

Results are presented of a computer simulation study of electron cyclotron harmonic waves; an electrostatic finite size particle simulation code is used. The initial electron velocity distribution function is modeled as a hot ring perpendicular to the ambient magnetic field and a cold Maxwellian component. Nonlinear cyclotron resonance is found to be the saturation mechanism of the unstable waves and the heating mechanism of the cold electrons. The results can be related to certain ionospheric phenomena.

Ashour-Abdalla, M.↗

Simulations and observations of heating of auroral ion beams

Two-dimensional three-velocity electrostatic particle simulations were used to determine the nonlinear evolution of the distributions of auroral ion beams and thereby to determine quantitatively signatures in the ion distributions produced by the ion-ion instability for a variety of plasma conditions in the auroral zone. The signatures determined from these simulations were compared with observations from DE 1, making it possible to characterize semiquantitatively the heating of the ionospheric ions, and to investigate the causes of variability seen in the observations of Reiff et al. (1988). A comparison of the simulation with observations showed features consistent with heating via the ion-ion instability including perpendicular heating in the supersonic regime and parallel heating in the subsonic regime, and a change in the heating between these regimes as the ratio of the H(+) beam speed to the local sound speed decreases.

Winglee, R. M.↗

Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability

Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).

Dum, C. T.↗

Launch Environment Water Flow Simulations Using Smoothed Particle Hydrodynamics

This paper describes the use of Smoothed Particle Hydrodynamics (SPH) to simulate the water flow from the rainbird nozzle system used in the sound suppression system during pad abort and nominal launch. The simulations help determine if water from rainbird nozzles will impinge on the rocket nozzles and other sensitive ground support elements.

Smoothed Particle Hyrodynamics↗

Simulation of ion-cyclotron-like modes in a magnetoplasma with transverse inhomogeneous electric field

The new kinetic ion-cyclotronlike (IC-like) instability driven by a localized electric field perpendicular to the ambient magnetic field is investigated using particle simulation. The existence of the new mode is demonstrated and the mode is distinguished from the well-known Kelvin-Helmholtz branch and the kinetic modification of this branch. The issue of the initial loading of particles in phase space in the simulation is clarified. The simulation results show that the growing ion waves are associated with the small vortices in the linear state, which evolve into a nonlinear stage dominated by large vortices with lower frequencies.

Nishikawa, K.-I.↗

Numerical studies of multiple-encounter shock drift acceleration

Test particle simulation is performed to examine the mechanism previously proposed by Chiueh for particle acceleration by a nearly perpendicular shock where only downstream MHD turbulence is present. A nonnegligible fraction of low-energy particles is observed gaining energy by more than three orders of magnitude in a few hundred cyclotron periods. Criteria for effective acceleration given previously are examined and confirmed. The numerical results suggest that in the nonrelativistic energy range, a power-law particle energy spectrum, with the power index determined by the degree to which the criteria are satisfied, may form.

Chiueh, Tzihong↗

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗

Simulation of the current-driven electrostatic ion cyclotron instability

Results are presented of a self-consistent particle simulation of the growth and saturation of electrostatic ion cyclotron (EIC) waves along auroral field lines. The driving mechanism for the waves is an initial drifting Maxwellian distribution for the electrons. The magnetic field is taken to lie in the simulation plane, so that modes such as the ion acoustic instability and the EIC instability can be excited, thus making it possible to study the behavior of EIC waves in the presence of ion acoustic turbulence. A dc electric field, which tends to keep the current constant along the field lines, is included. It is shown that the EIC instability grows to a moderate level before saturating due to plateau formation. The perpendicular ion temperature increases by typically 40%, with an increase arising from bulk heating only. The anomalous resistivity, associated with the EIC instability is determined and compared with previous calculations and with that due to the ion acoustic instability. Implications of the results for observation of ion cyclotron waves on auroral field lines are discussed.

Pritchett, P. L.↗

Nonlinear longitudinal resonance interaction of energetic charged particles and VLF waves in the magnetosphere

The longitudinal resonance of waves and energetic electrons in the Earth's magnetosphere, and the possible role this resonance may play in generating various magnetospheric phenomena are studied. The derivation of time-averaged nonlinear equations of motion for energetic particles longitudinally resonant with a whistler mode wave propagating with nonzero wave normal is considered. It is shown that the wave magnetic forces can be neglected at lower particle pitch angles, while they become equal to or larger than the wave electric forces for alpha 20 deg. The time-averaged equations of motion were used in test particle simulation which were done for a wide range of wave amplitudes, wave normals, particle pitch angles, particle parallel velocities, and in an inhomogeneous medium such as the magnetosphere. It was found that there are two classes of particles, trapped and untrapped, and that the scattering and energy exchange for those two groups exhibit significantly different behavior.

Tkalcevic, S.↗

Current Collection in a Magnetic Field

It is found that the upper-bound limit for current collection in the case of strong magnetic field from the current is close to that given by the Parker-Murphy formula. This conclusion is consistent with the results obtained in laboratory experiments. This limit weakly depends on the shape of the wire. The adiabatic limit in this case will be easily surpassed due to strong magnetic field gradients near the separatrix. The calculations can be done using the kinetic equation in the drift approximation. Analytical results are obtained for the region where the Earth's magnetic field is dominant. The current collection can be calculated (neglecting scattering) using a particle simulation code. Dr. Singh has agreed to collaborate, allowing the use of his particle code. The code can be adapted for the case when the current magnetic field is strong. The needed dm for these modifications is 3-4 months. The analytical description and essential part of the program is prepared for the calculation of the current in the region where the adiabatic description can be used. This was completed with the collaboration of Drs. Khazanov and Liemohn. A scheme of measuring the end body position is also proposed. The scheme was discussed in the laboratory (with Dr. Stone) and it was concluded that it can be proposed for engineering analysis.

Krivorutsky, E. N.↗

Simulations of beam excited minor species gyroharmonics in the Porcupine experiment

An active experiment for the study of wave particle interaction in the presence of a perpendicular beam type distribution is considered. The experiment involves the injection of a beam perpendicular to the magnetic field. Beam parameters can then be varied while local plasma parameters and associated wave phenomena are measured. Such an experiment was conducted as part of the Porcupine project. The present investigation has the objective to explain theoretically and by means of plasma particle simulations, some of the wave observations obtained. In the experiment, a main rocket payload and four subpayloads were employed. The subpayloads were ejected radially in different directions at an altitude of 240 km. One of the spinning subpayloads contained a plasma gun which emitted a beam of 200 eV xenon ions toward the main payload. The data recorded on the main payload revealed the existence of electrostatic waves with peaks at hydrogen gyroharmonics. A series of computer simulations approximating the xenon experiment was performed. The interpretation of the data is discussed.

Roth, I.↗

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗

Simulations of active ion injection experiments on ARCS 3

The ARCS 3 sounding rocket experiment used two argon guns oriented parallel and perpendicular to the geomagnetic field. The parameters of this experiment have been used in a series of particle simulations in order to interpret ion beam interaction with the multispecies ionospheric plasma. Along with the lower hybrid mode between the electron and the hydrogen gyrofrequencies, ion-ion hybrid modes between the hydrogen and helium gyrofrequencies and between the helium and oxygen gyrofrequencies are identified in the simulations. Besides these cold plasma modes, warm Bernstein modes of the respective ion species are identified. The parallel argon beam injections excite lower hybrid waves perpendicular to the magnetic field because of the finite thermal spread of the parallel beam. Perpendicular ion heating of background hydrogen, helium, and oxygen is comparable in the simulations. Comparisons are made with the wave and particle data from the ARCS 3 rocket experiment.

Hudson, M. K.↗

Windblown sand on Venus - Preliminary results of laboratory simulations

Small particles and winds of sufficient strength to move them have been detected from Venera and Pioneer-Venus data and suggest the existence of aeolian processes on Venus. The Venus wind tunnel (VWT) was fabricated in order to investigate the behavior of windblown particles in a simulated Venusian environment. Preliminary results show that sand-size material is readily entrained at the wind speeds detected on Venus and that saltating grains achieve velocities closely matching those of the wind. Measurements of saltation threshold and particle flux for various particle sizes have been compared with theoretical models which were developed by extrapolation of findings from Martian and terrestrial simulations. Results are in general agreement with theory, although certain discrepancies are apparent which may be attributed to experimental and/or theoretical-modeling procedures. Present findings enable a better understanding of Venusian surface processes and suggest that aeolian processes are important in the geological evolution of Venus.

Greeley, R.↗

Simulation studies of plasma waves in the electron foreshock - The generation of Langmuir waves by a gentle bump-on-tail electron distribution

Particle simulation experiments were used to study the basic physical ingredients needed for building a global model of foreshock wave phenomena. In particular, the generation of Langmuir waves by a gentle bump-on-tail electron distribution is analyzed. It is shown that, with appropriately designed simulations experiments, quasi-linear theory can be quantitatively verified for parameters corresponding to the electron foreshock.

Dum, C. T.↗