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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Digital Distortion Caused by Traveling- Wave-Tube Amplifiers Simulated

Future NASA missions demand increased data rates in satellite communications for near real-time transmission of large volumes of remote data. Increased data rates necessitate higher order digital modulation schemes and larger system bandwidth, which place stricter requirements on the allowable distortion caused by the high-power amplifier, or the traveling-wave-tube amplifier (TWTA). In particular, intersymbol interference caused by the TWTA becomes a major consideration for accurate data detection at the receiver. Experimentally investigating the effects of the physical TWTA on intersymbol interference would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to acquiring the required digital hardware. Thus, an accurate computational model is essential to predict the effects of the TWTA on system-level performance when a communication system is being designed with adequate digital integrity for high data rates. A fully three-dimensional, time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell code MAFIA (Solution of Maxwell's equations by the Finite-Integration-Algorithm). It comprehensively takes into account the effects of frequency-dependent AM (amplitude modulation)/AM and AM/PM (phase modulation) conversion, gain and phase ripple due to reflections, drive-induced oscillations, harmonic generation, intermodulation products, and backward waves. This physics-based TWT model can be used to give a direct description of the effects of the nonlinear TWT on the operational signal as a function of the physical device. Users can define arbitrary excitation functions so that higher order modulated digital signals can be used as input and that computations can directly correlate intersymbol interference with TWT parameters. Standard practice involves using communication-system-level software packages, such as SPW, to predict if adequate signal detection will be achieved. These models use a nonlinear, black-box model to represent the TWTA. The models vary in complexity, but most make several assumptions regarding the operation of the high-power amplifier. When the MAFIA TWT interaction model was used, these assumptions were found to be in significant error. In addition, digital signal performance, including intersymbol interference, was compared using direct data input into the MAFIA model and using the system-level analysis tool SPW for several higher order modulation schemes. Results show significant differences in predicted degradation between SPW and MAFIA simulations, demonstrating the significance of the TWTA approximations made in the SPW model on digital signal performance. For example, a comparison of the SPW and MAFIA output constellation diagrams for a 16-ary quadrature amplitude modulation (16-QAM) signal (data shown only for second and fourth quadrants) is shown. The upper-bound degradation was calculated from the corresponding eye diagrams. In comparison to SPW simulations, the MAFIA data resulted in a 3.6-dB larger degradation.

Kory, Carol L.↗
NASA NTRS2006

Particle Demagnetization in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS1996

Numerical Simulations of Virtual Anodes in Ion Beams Emissions in Space

Standard theories have shown that when the current density in a one-dimensional electron or ion beam exceeds a critical value, a potential hump will show up. This hump will behave as a virtual electrode to the beam particles. Measurements have demonstrated the existence of virtual cathodes and anodes. Results are presented of particle-in-cell simulations of ion beam emissions in space.

Ion Beam Plasma Ionization Particles↗
NASA NTRS1996

3-D Simulations of NSTAR Ion Thruster Plasma Interactions

Described is a Particle-in-Cell with Monte Carlo Collision code developed to perform detailed three-dimensional ion thruster simulations. To capture the full kinetic behavior of ion thruster plumes, both the electrons and ions are treated as test particles. Simulation results are given of the NSTAR ion thruster under ground test and in space conditions. Numerical results are compared.

NSTAR Ion Thruster Plasma Ionization Particles↗
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 NTRS2006

Electrodynamics in a Very Thin Current Sheet Leading to Magnetic Reconnection

We study the formation of a very thin current sheet (CS) and associated plasma electrodynamics using three-dimensional (3-D) particle-in-cell simulations with ion to electron mass ratio M/m=1836. The CS is driven by imposed anti-parallel magnetic fields. The noteworthy features of the temporal evolution of the CS are the following: (i) Steepening of the magnetic field profile B,(z) in the central part of the CS, (ii) Generation of three-peak current distribution with the largest peak in the CS center as B,(z) steepens, (iii) Generation of converging electric fields forming a potential well in the CS center in which ions are accelerated. (iv) Electron and ion heating in the central part of the CS by current-driven instabilities (CDI). (v) Re-broadening of the CS due to increased kinetic plasma pressure in the CS center. (vi) Generation of electron temperature anisotropy with temperature perpendicular to the magnetic field being larger than the parallel one. (vii) Current disruption by electron trapping in an explosively growing electrostatic instability (EGEI) and electron tearing instability (ETI). (viii)The onset of EGEI coincides with an increase in the electron temperature above the temperature of the initially hot ions as well as the appearance of new shear in the electron drift velocity. (ix) Bifurcation of the central CS by the current disruption. (x) Magnetic reconnection (MR) beginning near the null in B, and spreading outward. (xi) Generation of highly energized electrons reaching relativistic speeds and having isotropic pitch-angle distribution in the region of reconnected magnetic fields. We compare some of these features of the current sheet with results from laboratory and space experiments.

Singh, Nagendra↗
NASA NTRS2006

High-Frequency Electrostatic Wave Generation and Transverse Ion Acceleration by Low Alfvenic Wave Components of BBELF Turbulence

Satellite observations in the auroral plasma have revealed that extremely low frequency (ELF) waves play a dominant role in the acceleration of electrons and ions in the auroral plasma. The electromagnetic components of the ELF (EMELF) waves are the electromagnetic ion cyclotron (EMIC) waves below the cyclotron frequency of the lightest ion species in a multi-ion plasma. Shear Alfv6n waves (SAWS) constitute the lowest frequency components of the ELF waves below the ion cyclotron frequency of the heaviest ion. The -2 mechanism for the transfer of energy from such EMELF waves to ions affecting transverse ion heating still remains a matter of debate. A very ubiquitous fe8ture of ELF waves now observed in several rocket and satellite experiments is that they occur in conjunction with high-frequency electrostatic waves. The frequency spectrum of the composite wave turbulence extends from the low frequency of the Alfvenic waves to the high frequency of proton plasma frequency and/or the lower hybrid frequency. The spectrum does not show any feature organized by the ion cyclotron frequencies and their harmonics. Such broadband waves consisting of both the EM and ES waves are now popularly referred as BBELF waves. We present results here from 2.5-D particle-in-cell simulations showing that the ES components are directly generated by cross- field plasma instabilities driven by the drifts of the ions and electrons in the EM component of the BBELF waves.

Singh, Nagendra↗
NASA NTRS2006

Particle Acceleration, Magnetic Field Generation and Emission from Relativistic Jets and Supernova Remnants

We performed numerical simulations of particle acceleration, magnetic field generation, and emission from shocks in order to understand the observed emission from relativistic jets and supernova remnants. The investigation involves the study of collisionless shocks, where the Weibel instability is responsible for particle acceleration as well as magnetic field generation. A 3-D relativistic particle-in-cell (RPIC) code has been used to investigate the shock processes in electron-positron plasmas. The evolution of theWeibe1 instability and its associated magnetic field generation and particle acceleration are studied with two different jet velocities (0 = 2,5 - slow, fast) corresponding to either outflows in supernova remnants or relativistic jets, such as those found in AGNs and microquasars. Slow jets have intrinsically different structures in both the generated magnetic fields and the accelerated particle spectrum. In particular, the jet head has a very weak magnetic field and the ambient electrons are strongly accelerated and dragged by the jet particles. The simulation results exhibit jitter radiation from inhomogeneous magnetic fields, generated by the Weibel instability, which has different spectral properties than standard synchrotron emission in a homogeneous magnetic field.

Nishikawa, K.-I.↗
NASA NTRS2006

How Does Collisionless Magnetic Reconnection Work in the Presence of a Guide Magnetic Field?

The dissipation mechanism of guide field magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide-field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS2007

Electrostatic Wave Generation and Transverse Ion Acceleration by Alfvenic Wave Components of BBELF Turbulence

We present results here from 2.5-D particle-in-cell simulations showing that the electrostatic (ES) components of broadband extremely low frequency (BBELF) waves could possibly be generated by cross-field plasma instabilities driven by the relative drifts between the heavy and light ion species in the electromagnetic (EM) Alfvenic component of the BBELF waves in a multi-ion plasma. The ES components consist of ion cyclotron as well as lower hybrid modes. We also demonstrate that the ES wave generation is directly involved in the transverse acceleration of ions (TAI) as commonly measured with the BBELF wave events. The heating is affected by ion cyclotron resonance in the cyclotron modes and Landau resonance in the lower hybrid waves. In the simulation we drive the plasma by the transverse electric field, E(sub y), of the EM waves; the frequency of E(sub y), omega(sub d), is varied from a frequency below the heavy ion cyclotron frequency, OMEGA(sub h), to below the light ion cyclotron frequency, OMEGA(sub i). We have also performed simulations for E(sub y) having a continuous spectrum given by a power law, namely, |Ey| approx. omega(sub d) (exp -alpha), where the exponent alpha = _, 1, and 2 in three different simulations. The driving electric field generates polarization and ExB drifts of the ions and electrons. When the interspecies relative drifts are sufficiently large, they drive electrostatic waves, which cause perpendicular heating of both light and heavy ions. The transverse ion heating found here is discussed in relation to observations from Cluster, FAST and Freja.

Singh, Nagendra↗
NASA NTRS2006

e+/- Pair Loading and the Origin of the Upstream Field in GRB Shocks

We investigate here the effects of plasma instabilities driven by rapid e(sup plus or minus) pair cascades, which arise in the environment of GRB sources as a result of back-scattering of a seed fraction of their original spectrum. The injection of e(sup plus or minus) pairs induces strong streaming motions in the ambient medium. One therefore expects the pair-enriched medium ahead of the forward shock to be strongly sheared on length scales comparable to the radiation front thickness. Using three-dimensional particle-in-cell simulations, we show that plasma instabilities driven by these streaming e(sup plus or minus) pairs are responsible for the excitation of near-equipartition, turbulent magnetic fields. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective coupling between e(sup plus or minus) pairs and ions, and may help explain the origin of large upstream fields in GRB shocks.

Ramirez-Ruiz, Enrico↗
NASA NTRS2007

Simulation Study of Magnetic Fields generated by the Electromagnetic Filamentation Instability driven by Pair Loading

Using a 3-D relativistic particle-in-cell (RPIC) code, we have investigated particle acceleration associated with a relativistic electron-positron (cold) jet propagating into ambient electron-positron and electron-ion plasmas without initial magnetic fields in order to investigate the nonlinear stage of the Weibel instability. We have also performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which are assumed to be created by the photon annihilation. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions and ambient plasma. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-ion ambient plasma are larger than those in the electron-positron ambient plasma. We have shown that plasma instabilities driven by these streaming electron-positron pairs are responsible for the excitation of near-equipartition, turbulent magnetic fields. These fields maintain a strong saturated level on timescales much longer than the electron skin depth at least for the duration of the simulations. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective coupling between electron-positron pairs and ions, and may help explain the origin of large upstream fields in GRB shock.

Nishikawa, K.-I.↗
NASA NTRS2007

GRMHD and GRPIC Simulations

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous code. The simulation results show the jet formations from a geometrically thin accretion disk near a nonrotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field configuration including issues for future research. A General Relativistic Particle-in-Cell Code (GRPIC) has been developed using the Kerr-Schild metric. The code includes kinetic effects, and is in accordance with GRMHD code. Since the gravitational force acting on particles is extreme near black holes, there are some difficulties in numerically describing these processes. The preliminary code consists of an accretion disk and free-falling corona. Results indicate that particles are ejected from the black hole. These results are consistent with other GRMHD simulations. The GRPIC simulation results will be presented, along with some remarks and future improvements. The emission is calculated from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by GRMHD simulations considering thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We would like to extend this research using GRPIC simulations and examine a possible new mechanism for certain X-ray quasi-periodic oscillations (QPOs) observed in blackhole X-ray binaries.

Nishikawa, K.-I.↗
NASA NTRS2007

High Performance Simulations of Accretion Disk Dynamics and Jet Formations Around Kerr Black Holes

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

Nishikawa, Ken-Ichi↗
NASA NTRS2006

e(sup +/-) Pair Loading and the Origin of the Upstream Magnetic Field in GRB Shocks

We investigate here the effects of plasma instabilities driven by rapid e(sup +/-) pair cascades, which arise in the environment of GRB sources as a result of back-scattering of a seed fraction of their original spectrum. The injection of e(sup +/-) pairs induces strong streaming motions in the ambient medium. One therefore expects the pair-enriched medium ahead of the forward shock to be strongly sheared on length scales comparable to the radiation front thickness. Using three-dimensional particle-in-cell simulations, we show that plasma instabilities driven by these streaming e(sup +/-) pairs are responsible for the excitation of near-equipartition, turbulent magnetic fields. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective coupling between e(sup +/-) pairs and ions, and may help explain the origin of large upstream fields in GRB shocks.

Ramirez-Ruiz, Enrico↗
NASA NTRS2000

Effect of Geometric Azimuthal Asymmetrics of PPM Stack on Electron Beam Characteristics

A three-dimensional (3D) beam optics model has been developed using the electromagnetic particle-in-cell (PIC) code MAFIA. The model includes an electron beam with initial transverse velocity distribution focused by a periodic permanent magnet (PPM) stack. All components of the model are simulated in three dimensions allowing several azimuthally asymmetric traveling wave tube (TWT) characteristics to be investigated for the first time. These include C-magnets, shunts and magnet misalignment and their effects on electron beam behavior. The development of the model is presented and 3D TWT electron beam characteristics are compared in the absence of and under the influence of the azimuthally asymmetric characteristics described.

Kory, Carol L.↗
NASA NTRS2003

UCLA parallel PIC framework

The UCLA Parallel PIC Framework (UPIC) has been developed to provide trusted components for the rapid construction of new, parallel Particle-in-Cell (PIC) codes. The Framework uses object-based ideas in Fortran95, and is designed to provide support for various kinds of PIC codes on various kinds of hardware. The focus is on student programmers. The Framework supports multiple numerical methods, different physics approximations, different numerical optimizations and implementations for different hardware. It is designed with "defensive" programming in mind, meaning that it contains many error checks and debugging helps. Above all, it is designed to hide the complexity of parallel processing. It is currently being used in a number of new Parallel PIC codes.

Norton, Charles D.↗