Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle-in-cell 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Simulation of auroral double layers

Some basic properties of plasma double layers are deduced from a particle-in-cell computer simulation and related to parallel electric-field structures above the auroral regions. The simulation results on the processes leading to double-layer formation are examined, particularly in relation to the transient stage and double-layer structure and stability. It is concluded that: (1) a large potential difference applied to a finite-length plasma will be concentrated in a shocklike localized region instead of occurring over the entire length of the system; (2) the initial stage in double-layer formation is dominated by a large-potential pulse propagating in the direction of the induced electrostatic drift; (3) the entire potential is dropped over a specific scale length once the double layer has formed; and (4) this scale length is expected to be of the order of 1 km for a double layer above a discrete auroral arc with a potential of 10 kV and the electric-field vector parallel to the magnetic-field vector.

Hubbard, R. F.↗
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 NTRS1985

Self-consistent simulation of plasma interactions with secondary-emitting insulators

A cylindrical particle-in-cell (PIC) plasma simulation code applicable to plasma densities encountered in low Earth orbit (LEO) is described. The simulated geometries include that of a plain disk and a disk surrounded by a dielectric. Both configurations are mounted upon a ground plate in contact with a plasma environment. Techniques allowing simulation of dielectric charging using PIC time scales are discussed. Current versus voltage characteristic curves are calculated and the results are compared to experimental data.

Brandon, S. T.↗
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 NTRS1987

Plasma simulation using the massively parallel processor

Two dimensional electrostatic simulation codes using the particle-in-cell model are developed on the Massively Parallel Processor (MPP). The conventional plasma simulation procedure that computes electric fields at particle positions by means of a gridded system is found inefficient on the MPP. The MPP simulation code is thus based on the gridless system in which particles are assigned to processing elements and electric fields are computed directly via Discrete Fourier Transform. Currently, the gridless model on the MPP in two dimensions is about nine times slower that the gridded system on the CRAY X-MP without considering I/O time. However, the gridless system on the MPP can be improved by incorporating a faster I/O between the staging memory and Array Unit and a more efficient procedure for taking floating point sums over processing elements. The initial results suggest that the parallel processors have the potential for performing large scale plasma simulations.

Lin, C. S.↗
NASA NTRS1985

Charging effects in the cometary environment of Halley

Electrostatic charging of the Giotto spacecraft in different impact induced charged particle environments was studied with 3D numerical particle-in-cell models. The simulation results are assessed according to first experimental results of Giotto and Vega instruments. Decreasing and even negative spacecraft potentials measured on Vega near closest approach to the comet suggest the influence of the cometary plasma, which is also confirmed by corresponding numerical simulations.

Thiemann, H.↗
NASA NTRS1987

Particle simulation of auroral double layers

Work on the simulation of auroral double layers (DLs) with realistic particle-in-cell models is presented. An early model simulated weak DLs formed in a self-consistent circuit but under conditions subject to the ion-acoustic instability. Recent work has focused on strong DLs formed when currentless jets are injected into a dipole magnetic field.

Smith, Bruce L.↗
NASA NTRS1986

Studies on counterstreaming plasma expansion

Recent studies on counterstreaming plasma expansions are summarized. The basic phenomenon of plasma expansion is reviewed, and results from one-dimensional simulations of counterstreaming plasma expansion are discussed. Results from simulations based on an electrostatic particle-in-cell code, in which the dynamics of both the electrons and ions are exactly followed, are discussed. The formation of electrostatic shocks is addressed. Finally, results are presented on the ionospheric plasma expansion along the geomagnetic flux tubes by solving the hydrodynamic equations.

Singh, N.↗
NASA NTRS1985

Formation of the prelunar accretion disk

Considerations related to the angular momentum of the earth-moon system have led Ward and Cameron (1978) to the suggestion of a collisional origin of the moon. The required projectile would be about the mass of Mars (i.e., about 0.1 earth masses). The possibility has been considered that both the proto-earth and the projectile would be remnants of the formation of giant gaseous protoplanets. Processes of lunar formation studied take into account the viscous dissipation of an accretion disk. The present paper is concerned with the formation of the prelunar accretion disk. The hydrodynamics of a compressible viscous gas are simulated with the aid of the particle-in-cell method. The obtained results imply that the projectile is in an elliptical earth-crossing orbit.

Cameron, A. G. W.↗
NASA NTRS1983

Numerically simulated two-dimensional auroral double layers

A magnetized 2 1/2-dimensional particle-in-cell system which is periodic in one direction and bounded by reservoirs of Maxwellian plasma in the other is used to numerically simulate electrostatic plasma double layers. For the cases of both oblique and two-dimensional double layers, the present results indicate periodic instability, Debye length rather than gyroradii scaling, and low frequency electrostatic turbulence together with electron beam-excited electrostatatic electron-cyclotron waves. Estimates are given for the thickness of auroral doule layers, as well as the separations within multiple auroral arcs. Attention is given to the temporal modulation of accelerated beams, and the possibilities for ion precipitation and ion conic production by the double layer are hypothesized. Simulations which include the atmospheric backscattering of electrons imply the action of an ionospheric sheath which accelerates ionospheric ions upward.

Borovsky, J. E.↗
NASA NTRS1987

Simulations of auroral plasma processes - Electric fields, waves and particles

Plasma processes driven by current sheets of finite thicknesses in an ambient magnetized plasma are studied using a 2 1/2 dimensional particle-in-cell code, and similarities are found between simulated plasma processes and those observed in the auroral plasma. Current sheets are shown to be bounded by large perpendicular electric fields occurring near their edges above the conducting boundary. Shaped potential structures form when the current sheets are narrow, and when the current sheets are wide, potential structures develop a significant parallel potential drop such that the electrons are accelerated upwards. Downward parallel electric fields of variable strength are noted in the downward current region, and double layer formation is seen in both narrow and wide current sheets. High frequency oscillations near the electron plasma frequency and its harmonic are seen, and low frequency waves are observed.

Singh, Nagendra↗
NASA NTRS1981

Simulation of a perpendicular bow shock

Simulations of a high-Mach-number shock with parameters typical of the earth's bow shock are performed. The simulations rely on a hybrid code in which the ions are treated kinetically using standard particle-in-cell techniques and the electrons are treated as a massless, charge neutralizing fluid. One spatial dimension and all velocity and field components are included in the calculation. The simulations reproduce the observed ion reflection and overshoots in the magnetic field and density, features that are shown to be closely associated with ion gyration. It is noted that gyrating ions play a crucial role in building up and maintaining overshoots in the potential, density, and magnetic field.

Leroy, M. M.↗
NASA NTRS1984

Numerical simulations of positively-biased probes and dielectric-conductor disks in a plasma

Plasma densities in the low-earth orbit range may be sufficient to cause difficulties for spacecraft operating at high voltages in this environment. The present investigation is concerned with the results of a continuing effort to develop a particle-in-cell (PIC), cylindrically-symmetric, 2-1/2 dimensional, self-consistent numerical simulation code. The simulation has the objective to explore the interactions of an ambient plasma with a conducting disk, which may be partially covered by a dielectric material. The disk and the surrounding dielectric material represent a hole in an insulator covering a conductor. Attention is given to a review of the simulation model, plain disk calculations, the disk and dielectric configuration, and the 'pinhole' effect.

Brandon, S. T.↗
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 NTRS1982

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗
NASA NTRS1978

Numerical simulation of plasma double layers

Numerical simulation results are presented for a plasma double layer, the computer model being a finite one-dimensional particle-in-cell plasma with specified potential difference across the system. A single pulse is formed which crosses the system with constant velocity; this is followed by the formation of a potential drop across a limited region of the plasma. An approximate expression relating the spatial extent of the double layer and the potential drop is presented. Electron and ion beams are generated which tend to lead to instabilities in the upstream and downstream regions.

Joyce, G.↗
NASA NTRS1987

Simulation of the electron acoustic instability for a finite-size electron beam system

Simulations of the electron acoustic instability for a finite-size electron beam system are performed with a particle-in-cell code to investigate the heating phenomena associated with the instability and the width of the heating region. The results show that the beam radiates electrostatic electron acoustic waves whose decay time outside the beam agrees with the spatial decay length derived from the linear dispersion equation. The ambient cold electrons in a diffusion region surrounding the beam are heated to a higher temperature by absorbing the radiated electron acoustic waves, with the heating occurring mainly in the parallel direction. In the heat diffusion region, the temperature of the cold electrons decreases with distance from the beam with a temperature gradient length smaller than the decay length of the wave energy. The results are discussed with respect to the DE 1 plasma and wave observations in the polar cusp region.

Lin, C. S.↗