SEARCH · Search NASA
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.
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.
A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.
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.
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.
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.
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.
The particle acceleration processes here studied are driven by the relaxation of unstable ion ring distributions; these produce strong wave activity at the lower hybrid resonance frequency which collapses, and forms energetic electron and ion tails. The results obtained are applied to the problem posed by the production of energetic particles by solar flares. The numerical simulation results thus obtained by a 2 1/2-dimensional particle-in-cell code show a simultaneous acceleration of electrons to 10-500 keV energies, and of ions to as much as the 1 MeV range; the energy of the latter is still insufficient to account for gamma-ray emission in the 4-6 MeV range, but furnish a seed population for further acceleration.
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.
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.
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.
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.
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.
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.
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.
Results of 2D plasma simulations are presented and related to auroral observations. The formation of V-shaped potentials is studied with a 2 1/2 dimensional electrostatic particle-in-cell code for a magnetized plasma. It is shown that amplitudes for perpendicular electric fields are larger than for parallel electric fields, and for Te less than 100 eV, the amplitudes are comparable to the electric fields associated with the electrostatic shocks observed from the S3-3 satellite. The excitation of electrostatic ion-cyclotron EIC waves which occurs in the region below the parallel potential drop is discussed. In auroral plasmas EIC waves are observed above the V-shaped double layers in association with ion beams and field-aligned currents. The results also show that oppositely directed electric fields in the center and at the edges of the simulation region produce oppositely directed currents. Precipitating auroral ions in association with electron inverted-V events are seen by the DMSP-F6 satellite.
A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.
A Fourier-Fourier transformed version of the splitting algorithm for simulating solutions of the Vlasov-Poisson system of equations is introduced. It is shown that with the inclusion of filamentation filtration in this transformed algorithm it is both faster and more stable than the standard splitting algorithm. It is further shown that in a scalar computer environment this new algorithm is approximately equal in speed and far less noisy than its particle-in-cell counterpart. It is conjectured that in a multiprocessor environment the filtered splitting algorithm would be faster while producing more precise results.
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.