Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell algorithm”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

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

A general concurrent algorithm for plasma particle-in-cell simulation codes

The general concurrent particle-in-cell (GCPIC) algorithm has been used to implement an electrostatic particle-in-cell code on a 32-node hypercube parallel computer. The GCPIC algorithm decomposes the PIC code by dividing the particle simulation physical domain into subdomains that are equal in number to the number of processors; all subdomains will accordingly possess approximately equal numbers of particles. The portion of the code which updates particle positions and velocities is nearly 100 percent efficient when the number of particles increases linearly with that of hypercube processors.

Liewer, Paulett C.↗
NASA NTRS1992

Electrostatic Particle-In-Cell Code For Hypercube Computer

Code simulates two-dimensional motions of plasma particles in self-consistent electrostatic field and externally applied magnetic field developed for execution on Mark IIIfp hypercube computer. Based on generalization of one-dimensional particle-in-cell algorithm applicable to many different parallel computing architectures. Intermediate product of continuing effort to speed particle-in-cell computations by taking advantage of distributed-memory parallel computers like those of hypercube class.

Ferraro, Robert D.↗
NASA NTRS1990

A parallel particle-in-cell model for the massively parallel processor

The availability of the nearest-neighbor communication-incorporating Massively Parallel Processor has prompted the development of a two-dimensional, particle-in-cell algorithm which loads particles in a cell randomly onto a row of processors, filling only half of them with particles. Due to the simplification of communications among processors achieved in a row by the vacant processors and the random-particle sequence, the algorithm efficiently sorts particles and performs gather/scatter procedures for collecting charge density according to their cells. The algorithm calculates electric fields at the cells by FFT.

Lin, C. S.↗
NASA NTRS1997

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗
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 NTRS1995

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗
NASA NTRS1994

Extending HPF for advanced data parallel applications

The stated goal of High Performance Fortran (HPF) was to 'address the problems of writing data parallel programs where the distribution of data affects performance'. After examining the current version of the language we are led to the conclusion that HPF has not fully achieved this goal. While the basic distribution functions offered by the language - regular block, cyclic, and block cyclic distributions - can support regular numerical algorithms, advanced applications such as particle-in-cell codes or unstructured mesh solvers cannot be expressed adequately. We believe that this is a major weakness of HPF, significantly reducing its chances of becoming accepted in the numeric community. The paper discusses the data distribution and alignment issues in detail, points out some flaws in the basic language, and outlines possible future paths of development. Furthermore, we briefly deal with the issue of task parallelism and its integration with the data parallel paradigm of HPF.

Chapman, Barbara↗
NASA NTRS1994

A splitting algorithm for Vlasov simulation with filamentation filtration

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.

Klimas, A. J.↗
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 NTRS2001

Accurate Time-Dependent Traveling-Wave Tube Model Developed for Computational Bit-Error-Rate Testing

The phenomenal growth of the satellite communications industry has created a large demand for traveling-wave tubes (TWT's) operating with unprecedented specifications requiring the design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies on computational modeling. However, the TWT industry's computational modeling capabilities need to be improved because there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This limits the analysis and design of novel devices or TWT's with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. To address these concerns, a fully three-dimensional, time-dependent, helical TWT interaction model was developed using the electromagnetic particle-in-cell code MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm). The model includes a short section of helical slow-wave circuit with excitation fed by radiofrequency input/output couplers, and an electron beam contained by periodic permanent magnet focusing. A cutaway view of several turns of the three-dimensional helical slow-wave circuit with input/output couplers is shown. This has been shown to be more accurate than conventionally used two-dimensional models. The growth of the communications industry has also imposed a demand for increased data rates for the transmission of large volumes of data. To achieve increased data rates, complex modulation and multiple access techniques are employed requiring minimum distortion of the signal as it is passed through the TWT. Thus, intersymbol interference (ISI) becomes a major consideration, as well as suspected causes such as reflections within the TWT. To experimentally investigate effects of the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers, in addition to acquiring the required digital hardware. As an alternative, the time-domain TWT interaction model developed here provides the capability to establish a computational test bench where ISI or bit error rate can be simulated as a function of TWT operating parameters and component geometries. Intermodulation products, harmonic generation, and backward waves can also be monitored with the model for similar correlations. The advancements in computational capabilities and corresponding potential improvements in TWT performance may prove to be the enabling technologies for realizing unprecedented data rates for near real time transmission of the increasingly larger volumes of data demanded by planned commercial and Government satellite communications applications. This work is in support of the Cross Enterprise Technology Development Program in Headquarters' Advanced Technology & Mission Studies Division and the Air Force Office of Scientific Research Small Business Technology Transfer programs.

Kory, Carol L.↗
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 NTRS2001

Traveling-Wave Tube Amplifier Model to Predict High-Order Modulation Intersymbol Interference

Demands for increased data rates in satellite communications necessitate higher order modulation schemes, larger system bandwidth, and minimum distortion of the modulated signal as it is passed through the traveling wave tube amplifier (TWTA). One type of distortion that the TWTA contributes to is intersymbol interference (ISI), and this becomes particularly disruptive with wide-band, complex modulation schemes. It is suspected that in addition to the dispersion of the TWT, frequency dependent reflections due to mismatches within the TWT are a significant contributor to ISI. To experimentally investigate the effect of these mismatches within the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to the acquisition of the required digital hardware. In an attempt to develop a more accurate model to correlate IS1 with the TWTA and the operational signal, a fully three-dimensional (3D), time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell (PIC) code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm). The model includes a user defined slow-wave circuit with a spatially tapered region of loss to implement a sever, and spatially varied geometry (such as helical pitch) to implement a phase velocity taper. The model also includes user defined input/output coupling and an electron beam contained by solenoidal, electrostatic, or periodic permanent magnet (PPM) focusing allowing standard or novel TWTs to be investigated. This model comprehensively takes into account the effects of frequency dependent nonlinear distortions (MAM and AMPM); gain ripple due to frequency dependent reflections at the input/output coupling, severs, and mismatches from dynamic pitch variations; drive induced oscillations; harmonic generation; intermodulation products; and backward waves.

Kory, Carol L.↗
NASA NTRS2000

Effects of Geometric Azimuthal Asymmetries of the PPM Stack on Electron Beam Characteristics

The effects of geometric azimuthally asymmetric properties of a periodic permanent magnet (PPM) focusing stack on electron beam characteristics obtained using a fully three dimensional (3D) particle-in-cell (PIC) code will be presented. The simulation model, using MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm), incorporates 3D behavior of the beam immersed in static fields calculated directly from the exact geometry and material properties of the 3D magnetic focusing structure. The Hughes 8916H, 18-40 GHz helical TWT for the millimeter-wave power module (MMPM) was used as a prototype. Firstly, the effects of C-magnets used at the input and output of the TWT to allow for coupling of the RF signal into and out of the tube are considered. The 8916H input and output C-magnets differ because coaxial couplers are used at the input and waveguide couplers are used at the output The repositioning of the beam from its central axis due to the inclusion of the output C-magnet was found to be most significant. The modeled output C-magnet and its orientation in the Cartesian coordinate system is shown, and a two-dimensional beam profile including the output C-magnet is also shown. A table presents the shift of the beam center off the central axis relative to the average radius of the beam at the longitudinal points A, B and C designated on an enclosed figure. Secondly, the addition of shunts, or rectangular iron pieces applied manually by a skilled technician in order to improve beam transmission, is considered. The shunts are applied to the top of the tube; thus, azimuthal symmetry of the focusing stack is interrupted. Although shunts are typically added during RF focusing, they are also typically added at the input section of the tube where RF forces are minimal, making an electron optics analysis meaningful. Because several shunts are usually applied to one pole piece, the simulations have been simplified by modeling a half washer with the same radius and longitudinal length as a shunt over the entire x, positive-y half of the transverse plane. A modeled pole piece and shunt as described are shown. Lastly, in order to study the effects of magnet misalignments, a magnet in the PPM stack was arbitrarily chosen and adjusted so that its central axis was shifted both 0.7 percent and 1.0 percent of the magnet outer diameter in the positive-y direction. In practice, positioning the magnets so that their central axis is accurately aligned with the central axis of the tube is challenging. Thus, it is a strong possibility that one or more magnets will be misaligned relative to the tube central axis.

Kory, Carol L.↗
NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗
NASA NTRS2000

Three-Dimensional Electron Optics Model Developed for Traveling-Wave Tubes

A three-dimensional traveling-wave tube (TWT) electron beam optics model including periodic permanent magnet (PPM) focusing has been developed at the NASA Glenn Research Center at Lewis Field. This accurate model allows a TWT designer to develop a focusing structure while reducing the expensive and time-consuming task of building the TWT and hot-testing it (with the electron beam). In addition, the model allows, for the first time, an investigation of the effect on TWT operation of the important azimuthally asymmetric features of the focusing stack. The TWT is a vacuum device that amplifies signals by transferring energy from an electron beam to a radiofrequency (RF) signal. A critically important component is the focusing structure, which keeps the electron beam from diverging and intercepting the RF slow wave circuit. Such an interception can result in excessive circuit heating and decreased efficiency, whereas excessive growth in the beam diameter can lead to backward wave oscillations and premature saturation, indicating a serious reduction in tube performance. The most commonly used focusing structure is the PPM stack, which consists of a sequence of cylindrical iron pole pieces and opposite-polarity magnets. Typically, two-dimensional electron optics codes are used in the design of magnetic focusing devices. In general, these codes track the beam from the gun downstream by solving equations of motion for the electron beam in static-electric and magnetic fields in an azimuthally symmetric structure. Because these two-dimensional codes cannot adequately simulate a number of important effects, the simulation code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm) was used at Glenn to develop a three-dimensional electron optics model. First, a PPM stack was modeled in three dimensions. Then, the fields obtained using the magnetostatic solver were loaded into a particle-in-cell solver where the fully three-dimensional behavior of the beam was simulated in the magnetic focusing field. For the first time, the effects of azimuthally asymmetric designs and critical azimuthally asymmetric characteristics of the focusing stack (such as shunts, C-magnets, or magnet misalignment) on electron beam behavior have been investigated. A cutaway portion of a simulated electron beam focused by a PPM stack is illustrated.

Kory, Carol L.↗
NASA NTRS1998

Three-Dimensional Deformable Grid Electromagnetic Particle-in-cell for Parallel Computers

We describe a new parallel, non-orthogonal grid, three-dimensional electromagnetic particle-in-cell (EMPIC) code based on a finite-volume formulation. This code uses a logically Cartesian grid of deformable hexahedral cells, a discrete surface integral (DSI) algorithm to calculate the electromagnetic field, and a hybrid logical-physical space algorithm to push particles.

Cartesian grid Grid Electromagnetic electromagneti↗