Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell”

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
NASA NTRS1990

Modeling solar flare conduction fronts. I - Homogeneous plasmas and ion-acoustic turbulence. II - Inhomogeneous plasmas and ambipolar electric fields

A one-dimensional, electrostatic, particle-in-cell simulation is used here to model the expansion of a heated electron population in a coronal loop during a solar flare and the characteristics of the associated X-ray emissions. The hot electrons expand outward from the localized region, creating an ambipolar electric field which accelerates a return current of cooler, ambient electrons. Ion-acoustic waves are generated by the return currents as proposed by Brown et al. (1979), but they play little or no role in containing energetic electrons and the conduction front proposed by Brown et al. does not form. The X-ray emission efficiency of the electrons is too low in the corona for them to be the source of hard X-ray bursts. The particle dynamics changes dramatically if the heated plasma is at low altitudes and expands upward into the more tenuous plasma at higher altitudes. Two important applications of this finding are the radio-frequency heating of the corona and the collisional heating of the chromosphere by precipitating energetic electrons. In both cases, the overlying plasma has a density that is too low to supply a balancing return current to the expanding hot electrons. As a result, an ambipolar electric field develops that tends to confine the energetic electrons behind a front that propagate outward at about the speed of sound.

Mckean, M. E.↗
NASA NTRS1991

Monte Carlo simulation of ion-neutral charge exchange collisions and grid erosion in an ion thruster

A combined particle-in-cell (PIC)/Monte Carlo simulation model has been developed in which the PIC method is used to simulate the charge exchange collisions. It is noted that a number of features were reproduced correctly by this code, but that its assumption of two-dimensional axisymmetry for a single set of grid apertures precluded the reproduction of the most characteristic feature of actual test data; namely, the concentrated grid erosion at the geometric center of the hexagonal aperture array. The first results of a three-dimensional code, which takes into account the hexagonal symmetry of the grid, are presented. It is shown that, with this code, the experimentally observed erosion patterns are reproduced correctly, demonstrating explicitly the concentration of sputtering between apertures.

Peng, Xiaohang↗
NASA NTRS1990

Simulation studies of electron acceleration by ion ring distributions in solar flares

A 2.5-dimensional fully relativistic EM, particle-in-cell code (PIC) is used to investigate a potential electron acceleration mechanism in solar flares. The free energy is provided by ions which have a ring velocity distribution about the magnetic-field direction. Ion rings may be produced by perpendicular shocks, which could in turn be generated by the super-Alfvenic motion of magnetic flux tubes emerging from the photosphere or by coronal mass ejections. Such ion distributions are known to be unstable to the generation of lower hybrid waves, which have phase velocities in excess of the electron thermal speed parallel to the field and can, therefore, resonantly accelerate electrons in that direction. The simulations show the transfer of perpendicular ion energy to energetic electrons via lower hybrid wave turbulence. With plausible ion ring velocities, the process can account for the observationally inferred fluxes and energies of non-thermal electrons during the impulsive phase of flares.

Mcclements, K. G.↗
NASA NTRS1991

Preferential positron heating and acceleration by synchrotron maser instabilities in relativistic positron-electron-proton plasmas

A new process of the preferential strong heating of positrons through the ion synchrotron maser instability in positron-electron-proton magnetized plasmas is investigated using particle-in-cell simulations. It is shown that the positrons form a nonthermal power-law-like energy distribution via their gyroresonant interaction with the extraordinary modes emitted by the ions. It is noted that this process may be of significance in connection with the shock excitation of nonthermal synchrotron radiation from astrophysical systems powered by relativistic outflows from compact central objects, e.g., supernova remnants powered by pulsars and jets from active galactic nuclei.

Hoshino, Masahiro↗
NASA NTRS1991

Computer simulation of astrophysical plasmas

The role of sophisticated numerical models and simulations in the field of plasma astrophysics is discussed. The need for an iteration between microphysics and macrophysics in order for astrophysical plasma physics to produce quantitative results that can be related to astronomical data is stressed. A discussion on computational requirements for simulations of astrophysical plasmas contrasts microscopic plasma simulations with macroscopic system models. An overview of particle-in-cell simulations (PICS) is given and two examples of PICS of astrophysical plasma are discussed including particle acceleration by collisionless shocks in relativistic plasmas and magnetic field reconnection in astrophysical plasmas.

Max, Claire E.↗
NASA NTRS1991

Comparison of 2-D and 3-D models of grid erosion in an ion thruster

Numerical results of particle-in-cell/Monte Carlo calculations of accelerator grid erosion in an ion thruster are presented. Specifically, it is shown that a three-dimensional model is required to account for the experimentally observed pitting of the accelerator grid between grid apertures. Some comparisons with earlier two-dimensional, axisymmetric model are made, and it is shown that, for identical operating conditions of the thruster, the wear-through time in the three-dimensional model is about two to three times higher than that obtained previously with the two-dimensional model, namely on the order of 10,000 hours for sample calculation.

Peng, Xiaohang↗
NASA NTRS1991

The expansion of polarization charge layers into magnetized vacuum - Theory and computer simulations

The formation and evolution of polarization charge layers on cylindrical plasma streams moving in vacuum are investigated using analytic theory and 2D electrostatic particle-in-cell computer simulations. It is shown that the behavior of the electron charge layer goes through three stages. An early time expansion is driven by electrostatic repulsion of electrons in the charge layer. At the intermediate stage, the simulations show that the electron-charge-layer expansion is halted by the positively charged plasma stream. Electrons close to the stream are pulled back to the stream and a second electron expansion follows in time. At the late stage, the expansion of the ion charge layer along the magnetic field lines accompanies the electron expansion to form an ambipolar expansion. It is found that the velocities of these electron-ion expansions greatly exceed the velocities of ambipolar expansions which are driven by plasma temperatures.

Galvez, Miguel↗
NASA NTRS1991

A simulation study of interactions of space-shuttle generated electron beams with ambient plasma and neutral gas

The objective was to conduct large scale simulations of electron beams injected into space. The study of the active injection of electron beams from spacecraft is important, as it provides valuable insight into the plasma beam interactions and the development of current systems in the ionosphere. However, the beam injection itself is not simple, being constrained by the ability of the spacecraft to draw current from the ambient plasma. The generation of these return currents is dependent on several factors, including the density of the ambient plasma relative to the beam density, the presence of neutrals around the spacecraft, the configuration of the spacecraft, and the motion of the spacecraft through the plasma. Two dimensional (three velocity) particle simulations with collisional processes included are used to show how these different and often coupled processes can be used to enhance beam propagation from the spacecraft. To understand the radial expansion mechanism of an electron beam injected from a highly charged spacecraft, two dimensional particle-in-cell simulations were conducted for a high density electron beam injected parallel to magnetic fields from an isolated equipotential conductor into a cold background plasma. The simulations indicate that charge build-up at the beam stagnation point causes the beam to expand radially to the beam electron gyroradius.

Winglee, Robert M.↗
NASA NTRS1991

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1990

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1991

A model for the frequency fine structure of auroral kilometric radiation

A model for AKR frequency fine structure in which the fine structure is a natural consequence of the emission of electron cyclotron maser radiation in a nonuniform magnetic field is presented. The model is examined by using 1D electromagnetic particle-in-cell simulations. It is found that maser radiation is emitted as individual wave packets combine to form the drifting emission features that make up the fine structure. Rising frequency drifts are associated with wave packets emitted independently of one another. Wave packets that interact with one another produce features with falling frequency drift. The emission features produced have a large spread in both positive and negative frequency drift rates. The average positive drift rate of rapidly drifting features is found to be approximately proportional to the maser growth rate.

Mckean, M. E.↗
NASA NTRS1992

A hybrid Zakharov particle simulation of ionospheric heating

A 1D hybrid simulation model incorporating one of the Zakharov equations for the high-frequency waves and a particle-in-cell simulation of the ions is described and applied to ionospheric heating using realistic parameters. Results from the hybrid simulation are compared with a Zakharov simulation that incorporates a phenomenological model of ion damping. Both the hybrid and Zakharov simulations predict the formation of solitonlike waves. The early time behavior of the two simulations is different due to the high noise level in the hybrid simulation, but the late time behavior is quite similar.

Clark, K. L.↗
NASA NTRS1992

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗
NASA NTRS1992

Relativistic, perpendicular shocks in electron-positron plasmas

One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators.

Gallant, Yves A.↗
NASA NTRS1992

Effects of magnetospheric electrons on polar plasma outflow - A semikinetic model

The effect or hot magnetospheric electrons on the polar-plasma outflow was investigated, using a semikinetic model developed by Wilson et al. (1990) and Ho et al. (1991) to simulate the effect. The model is based on a hybrid particle-in-cell approach, in which the H(+) and O(+) ions are treated as adiabatic parallel-drifting gyrocenters injected as the upgoing portions of drifting bi-Maxwellian distributions at 1.6 R(E), while the electrons are treated as a massless neutralizing fluid. The results show that, in order to simulate the polar outflow under the influence of hot magnetospheric electrons, it is necessary to consider the effect of the electron temperature gradient.

Ho, C. W.↗
NASA NTRS1992

Further study of the effect of the downstream plasma condition on accelerator grid erosion in an ion thruster

Further numerical results are presented of earlier particle-in-cell/Monte Carlo calculations of accelerator grid erosion in an ion thruster. A comparison between numerical and experimental results suggests that the accelerator grid impingement is primarily due to ions created far downstream from the accelerator grid. In particular, for the same experimental conditions as those of Monheiser and Wilbur at Colorado State University, it is found that a downstream plasma density of 2 x 10 exp 14/cu m is required to give the same ratio of accelerator grid impingement current to beam current (5 percent). For this condition, a potential hill is found in the downstream region of 2.5 V.

Peng, Xiaohang↗
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗
NASA NTRS1991

A new approach to plasmasphere refilling: Anomalous plasma effects

During the last 10 months of the grant, both laminar and anomalous plasma processes occurring during the refilling of the outer plasmasphere after magnetic storms are investigated. Theoretical investigations were based on two types of models: (1) two-stream hydrodynamic model in which plasma flows from the conjugate ionospheres are treated as separate fluids and the ion temperature anisotropies are treated self-consistently; and (2) large-scale particle-in-cell code.

Singh, N.↗