Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasma PIC”

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

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.↗

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↗

The Heliopause Electrostatic Rapid Transit System (HERTS) Design, Trades, and Analyses Performed in the First Year of a Two Year Investigation

The Heliopause Electrostatic Rapid Transit System (HERTS)1 was one of the seven total Phase II NASA Innovative Advanced Concepts (NIAC) that was down-selected in 2015 for continued funding and research. In Phase I we learned that a spacecraft propelled by an Electric Sail (E-Sail) can travel great astronomical distances, such as to the Heliopause region of the solar system (approx.100 to 120 AU) in approximately one quarter of the time (10 years) versus the time it took the Voyager spacecraft launched in 1977 (36 years). The current work within the Phase II NIAC effort builds upon the work that was done in the Phase I NIAC and is focused on: 1) Testing of plasma interaction with a charged wire in a unique MSFC test chamber, 2) Development of a Particle-in-Cell (PIC) models that are validated in the plasma testing and used to extrapolate to the E-Sail propulsion system design. 3) Further down select of a wire deployment and control approach from those narrowed down in the Phase I effort. This paper will document the findings to date (June, 2016) of the above focused areas.

Wiegmann, Bruce M.↗

The Heliopause Electrostatic Rapid Transit System (HERTS) Design, Trades, and Analyses Performed in a Two Year NASA Investigation of Electric Sail Propulsion Systems

The Heliopause Electrostatic Rapid Transit System (HERTS) was one of the seven total Phase II NASA Innovative Advanced Concepts (NIAC) that was down-selected in 2015 for continued funding and research. In Phase I our team learned that a spacecraft propelled by an Electric Sail (E-Sail) can travel great astronomical distances, such as to the Heliopause region of the solar system (approx. 100 to 120 AU) in approximately one quarter of the time (10 years) versus the time it took the Voyager spacecraft launched in 1977 (36 years). The completed work within the Phase II NIAC funded effort builds upon the work that was done in the Phase I NIAC and is focused on: 1) Testing of plasma interaction with a charged wire in a MSFC simulated solar environment vacuum test chamber. 2) Development of a Particle-in-Cell (PIC) models that are validated in the plasma testing and used to extrapolate to the E-Sail propulsion system design. 3) Conceptual design of a Technology Demonstration Mission (TDM) spacecraft developed to showcase E-Sail propulsion systems. 4) Down selection of both: a) Materials for a multi km length conductor and, b) Best configuration of the proposed conductor deployment subsystem. This paper will document the findings to date (June, 2017) of the above focused areas.

Wiegmann, Bruce M.↗

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↗

Concurrent Algorithm For Particle-In-Cell Simulations

Separate decompositions used for particle-motion and field calculations. General Concurrent Particle-in-Cell (GCPIC) algorithm used to implement motions of individual plasma particles (ions and electrons) under influence of particle-in-cell (PIC) computer codes on concurrent processors. Simulates motions of individual plasma particles under influence of electromagnetic fields generated by particles themselves. Performed to study variety of nonlinear problems in plasma physics, including magnetic and inertial fusion, plasmas in outer space, propagation of electron and ion beams, free-electron lasers, and particle accelerators.

Liewer, Paulett C.↗

Electron cyclotron thruster new modeling results preparation for initial experiments

The following topics are discussed: a whistler-based electron cyclotron resonance heating (ECRH) thruster; cross-field coupling in the helicon approximation; wave propagation; wave structure; plasma density; wave absorption; the electron distribution function; isothermal and adiabatic plasma flow; ECRH thruster modeling; a PIC code model; electron temperature; electron energy; and initial experimental tests. The discussion is presented in vugraph form.

Hooper, E. Bickford↗

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↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the POLAR satellite is studied by using a 3-dimensional Particle-In-Cell (PIC) code. When the satellite passes through the region of low density plasma, the satellite charges to positive potentials as high as 4050Volts due to the photoelectrons emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, an ion-rich Xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at about 2Volts. Accordingly, in our 3-dimensional PIC simulation, we considered that the potential of the satellite is 2Volts as a fixed bias. Considering the relatively high density of the Xenon plasma in the sheath (approx. 10 - 10(exp 3)/cc), the ambient plasma of low density (less than 1/cc) is neglected. In the simulations, the electric fields and plasma dynamics are calculated self-consistently. We found that an "Apple" shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission, a high positive potential hill develops. Near the Thermal Ion Detection Experiment (TIDE) detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations, it takes about a hundred electron gyroperiods for the sheath to reach a quasi-steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. Using the steady state sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Leung, Wing C.↗

(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↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗

Merging for Particle-Mesh Complex Particle Kinetic Modeling of the Multiple Plasma Beams

We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.

Lipatov, Alexander S.↗

PIC code modeling of spacecraft charging potential during electron beam injection into a background of neutral gas and plasma, part 1

Injections of nonrelativistic electron beams from an isolated equipotential conductor into a uniform background of plasma and neutral gas were simulated using a 2-D electrostatic particle code. The ionization effects on spacecraft charging are examined by including interactions of electrons with neutral gas. The simulations show that the conductor charging potential decreases with increasing neutral background density due to the production of secondary electrons near the conductor surface. In the spacecraft wake, the background electrons accelerated towards the charged spacecraft produce an enhancement of secondary electrons and ions. Simulations run for longer times indicate that the spacecraft potential is further reduced and short wavelength beam-plasma oscillations appear. The results are applied to explain the spacecraft charging potential measured during the SEPAC experiments from Spacelab 1.

Koga, J. K.↗

NASA's Electric Sail Propulsion System Investigations over the Past Three Years

Personnel from NASA's MSFC have been investigating the feasibility of an advanced propulsion system known as the Electric Sail for future scientific missions of exploration. This team initially won a NASA Space Technology Mission Directorate (STMD) Phase I NASA Innovative Advanced Concept (NIAC) award and then a two year follow-on Phase II NIAC award. This paper documents the findings from this three year investigation. An Electric sail propulsion system is a propellant-less and extremely fast propulsion system that takes advantage of the ions that are present in the solar wind to provide very rapid transit speeds whether to deep space or to the inner solar system. Scientific spacecraft could arrive to Pluto in ~5 years, to the boundary of the solar system in ten to twelve years vs. thirty five plus years it took the Voyager spacecraft. The team's recent focused activities are: 1) Developing a Particle in Cell (PIC) numeric engineering model from the experimental data collected at MSFC's Solar Wind Facility on the interaction between simulated solar wind interaction with a charged bare wire that can be applied to a variety of missions, 2) The development of the necessary tether deployers/tethers to enable successful deployment of multiple, multi km length bare tethers, 3) Determining the different missions that can be captured from this revolutionary propulsion system 4) Conceptual designs of spacecraft to reach various destinations whether to the edge of the solar system, or as Heliophysics sentinels around the sun, or to trips to examine a multitude of asteroids These above activities, once demonstrated analytically, will require a technology demonstration mission (~2021 to 2023) to demonstrate that all systems work together seamlessly before a Heliophysics Electrostatic Rapid Transit System (HERTS) could be given the go-ahead. The proposed demonstration mission will require that a small spacecraft must first travel to cis-lunar space as the Electric Sail must be outside of Earth's Magnetic fields to produce thrust. The paper will outline what was done over the past three years from performing various plasma chamber tests to obtain data for the PIC model development, investigation of tether material trades, and conceptual designs of proposed spacecraft.

Wiegmann, Bruce M.↗

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗

Current collection in a flowing magnetoplasma

Effects of plasma drift on the current collection by a long conducting cylinder in a magnetized plasma is studied by means of a 2 1/2 dimensional PIC code. It is found that for the drift velocity (V sub 0) perpendicular to the magnetic field (B sub 0), the electron current collected by a positive cylinder is considerably enhanced depending on the drift velocity. The distributions of plasma and the potential structure around the cylinder for several relative orientations between V sub 0 and B are presented along with the comparisons of current with and without the magnetic field. Simulations with the magnetic field in the simulation plane show that the potential structures around the cylinder are two-dimensional double layers with dimension (L) perpendicular to B much smaller than the dimension (L) parallel to B. In fact, L perpendicular is found to be approximately determined by the current limiting radius given by the Parker-Murphy model. However, it is found that the collected currents in the simulations are generally higher than those given by this model.

Singh, Nagendra↗

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.↗