Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Plasma Wakes from the Baffle Scanning Mechanism on a LEO Mass Spectrometer

The Modular Spectrometer for Atmosphere and Ionosphere Characterization (MoSAIC) is a Quadrupole Mass Spectrometer (QMS) which will fly on NASA's Geospace Dynamics Constellation (GDC) mission in low Earth orbit (LEO), providing measurements of neutral and ionized gas density, temperature, composition, and wind/drift. MoSAIC includes a Baffle Scanning Mechanism (BSM) which periodically sweeps in front of the QMS aperture, enabling determination of the ion temperature and cross-track velocity. Accurate knowledge of the plasma wake formed by the BSM is critical for accurate determination of these quantities. We conduct particle-in-cell (PIC) simulations modeling the interactions between the baffle and plasma across the full range of ion densities, temperatures, and flow speeds GDC is expected to encounter. We discuss the wake physics and present an empirical model which can be implemented in the MoSAIC data analysis process to determine the relevant ion quantities.

mass spectrometer↗
NASA NTRS2003

Kinetic Simulation of Inhomogeneous Plasma with a Variable Sized Grid System

Space plasmas are usually inhomogeneous and irregular. Particle in cell (PIC) plasma simulations, however, often consider the plasma to be homogeneous and periodic. Here non-periodic PIC simulations with inhomogeneities in the density, temperature and background magnetic field are considered. Boundary conditions and a variable sized grid system are discussed, along with an application to a magnetospheric plasma physics problem in the auroral zone.

Schriver, David↗
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 NTRS1996

3-D Simulations of NSTAR Ion Thruster Plasma Environment

Full three-dimensional numerical models based on partical-in-cell with Monte Carlo collisions (PIC-MCC) simulations are developed to study ion thruster plasma interactions. 3-D simulation results for the NSTAR ion thruster plasma environment are presented.

NSTAR↗
NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗
NASA NTRS2014

Open Boundary Particle-in-Cell Simulation of Dipolarization Front Propagation

First results are presented from an ongoing open boundary 2-1/2D particle-in-cell simulation study of dipolarization front (DF) propagation in Earth's magnetotail. At this stage, this study is focused on the compression, or pileup, region preceding the DF current sheet. We find that the earthward acceleration of the plasma in this region is in general agreement with a recent DF force balance model. A gyrophase bunched reflected ion population at the leading edge of the pileup region is reflected by a normal electric field in the pileup region itself, rather than through an interaction with the current sheet. We discuss plasma wave activity at the leading edge of the pileup region that may be driven by gradients, or by reflected ions, or both; the mode has not been identified. The waves oscillate near but above the ion cyclotron frequency with wavelength several ion inertial lengths. We show that the waves oscillate primarily in the perpendicular magnetic field components, do not propagate along the background magnetic field, are right handed elliptically (close to circularly) polarized, exist in a region of high electron and ion beta, and are stationary in the plasma frame moving earthward. We discuss the possibility that the waves are present in plasma sheet data, but have not, thus far, been discovered.

dipolarization front↗
NASA NTRS1999

Generation of Electron Suprathermal Tails in the Upper Solar Atmosphere: Implications for Coronal Heating

We present a mechanism for the generation of non-Maxwellian electron distribution function in the upper regions of the solar atmosphere in the presence of collisional damping. It is suggested that finite amplitude, low frequency, obliquely propagating electromagnetic waves can carry a substantial electric field component parallel to the mean magnetic field that can be significantly larger than the Dreicer electric field. This long wavelength electric fluctuation is capable of generating high frequency electron plasma oscillations and low frequency ion acoustic-like waves. The analysis has been performed using 1-1/2D Vlasov and PIC numerical simulations in which both electrons and ions are treated kinetically and self consistently. The simulation results indicate that high frequency electron plasma oscillations and low frequency ion acoustic-like waves are generated. The high frequency electron plasma oscillation drives electron plasma turbulence, which subsequently is damped out by the background electrons. The turbulence damping results in electron acceleration and the generation of non-Maxwellian suprathermal tails on time scales short compared to collisional damping. Bulk heating also occurs if the fluctuating parallel electric field is strong enough. This study suggests that finite amplitude, low frequency, obliquely propagating, electromagnetic waves can play a significant role in the acceleration and heating of the solar corona electrons and in the coupling of medium and small-scale phenomena.

Vinas, Adolfo F.↗
NASA NTRS1995

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗
NASA NTRS2021

An Encounter with the Ion and Electron Diffusion Regions at a Flapping and Twisted Tail Current Sheet

We analyze data returned by the Magnetospheric Multiscale mission (MMS) constellation during a rapid (~1.5 s) traversal of a flapping and reconnecting current sheet (CS) in the near-Earth magnetotail (X ~ -20 RE). The CS was highly tilted, with its normal pointing strongly duskward. Its extreme thinness was confirmed by a curvature analysis of the magnetic field lines. The event was associated with a guide field of 8% of the reconnecting components. From the pitch angle distributions of low-energy electrons we infer that the crossing occurred earthward of the X-line. Traveling practically normal to the CS, MMS encountered an ion diffusion region (IDR) in which was embedded an electron diffusion region (EDR). IDR signatures included breaking of the ion frozen-in condition in the presence of Hall B and E fields. EDR signatures included a strong out of-plane current associated with a super-Alfvenic electron jet, evidence of positive energy transfer, and a temperature anisotropy (TeII > Te1) which disappeared at the field reversal. Derived scale sizes normal to the CS are: ~6.9 de (EDR) and ~0.4 di (IDR) (40 and 100 km). We estimate the average dimensionless reconnection rate as 0.077 + 0.050. The observations and inferences are supported by PIC numerical simulations. We find very good agreement in the reconnection rates. We also discuss the effects of asymmetries in the density, temperature and magnetic field strength on the Hall fields and length of the outflow jets. The event is associated with a substorm onset which began 7 min after the MMS observations.

C J Farrugia↗
NASA NTRS2021

The Micro‐Macro Coupling of Mass‐Loading in Symmetric Magnetic Reconnection With Cold Ions

We investigate how magnetic reconnection is influenced by an inflow of a dense cold ion population. We compare two 2.5D Particle-In-Cell simulations, one containing the cold population and one without. We find that the cold population influences the reconnection process on both global and kinetic scales, and that the dominant contribution can be explained through mass-loading. We provide an analysis of how these multiscale changes are related through kinetic processes in the ion diffusion region, the so-called micro-macro coupling of mass-loading. The inertia of the cold ion population is found to be the significant link that connects the changes on different scales. The cold and warm populations exhibit counter streaming behavior when and after the ion diffusion region reorganizes itself in response to the arrival of the cold population. This signature of the cold population should be observable by spacecraft observatories such as MMS.

Magnetic reconnection↗
NASA NTRSDate not supplied

Development of Candidate Airspace Procedures for Urban Air Mobility in the Dallas Area

As the concept of on-demand electric air taxis gains popularity and begins to take shape, many research efforts are underway to address the challenges of integrating this new class of passenger-carrying air vehicles into the already congested U.S. National Airspace System. Through collaborations between the National Aeronautics and Space Administration (NASA), the Federal Aviation Administration, and industry, this concept, commonly referred to as Urban Air Mobility (UAM), envisions a safe, reliable, and efficient mode of transportation traversing metropolitan and urban areas. Among the many challenges being tackled, airspace procedures and information requirements are critical areas of research that need to be addressed. In a joint effort between NASA and Joby Aviation, a human-in-the-loop study was conducted at the NASA Ames Research Center to evaluate initial and midterm operations with air traffic controllers and on-board UAM pilots in the Dallas area. This area was chosen due to its complex Class Bravo Airspace that extends to the surface over a relatively large area. To achieve the goals and objectives of the study, multiple preceding efforts were conducted to develop the candidate procedures and information requirements. This paper provides an overview of the operational concept used for the study, the process that was followed, and the findings from the two tabletop meetings and a shakedown activity. The combination of these led to a set of airspace procedures, letters of agreement, and information requirements that were evaluated in the human-inthe-loop study.

Urban Air Mobility↗
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 NTRS2011

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Using our new 3-D relativistic particle-in-cell (PIC) code, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic compression. Behind the bow shock, in the jet shock, strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to go beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations and calculated radiation based on first principles. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the technique to calculate emission from electrons based on simulations with a small system. We obtain spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system may generate a jitter/synchrotron spectrum.

Nishikawa, K.-I.↗
NASA NTRS1993

(abstract) 3D Electromagnetic Plasma Particle Simulations

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.

electromagnetic plasma particles simulations 3D al↗
NASA NTRS2016

Modeling an Iodine Hall Thruster Plume in the Iodine Satellite (ISAT)

An iodine-operated 200-W Hall thruster plume has been simulated using a hybrid-PIC model to predict the spacecraft surface-plume interaction for spacecraft integration purposes. For validation of the model, the plasma potential, electron temperature, ion current flux, and ion number density of xenon propellant were compared with available measurement data at the nominal operating condition. To simulate iodine plasma, various collision cross sections were found and used in the model. While time-varying atomic iodine species (i.e., I, I+, I2+) information is provided by HP Hall simulation at the discharge channel exit, the molecular iodine species (i.e., I2, I2+) are introduced as Maxwellian particles at the channel exit. Simulation results show that xenon and iodine plasma plumes appear to be very similar under the assumptions of the model. Assuming a sticking coefficient of unity, iodine deposition rate is estimated.

Iodine↗