Implementation of an Electromagnetic 3D Plasma Particle-in-Cell Code on a MIMD Parallel Computer
Explore the source record for details and available documents.
PIC GCPOC algorithm electrostatic field↗
Engineering topics
Liewer, P. C.
Publications and source records attributed to Liewer, P. 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↗
Hybrid simulations of the effects of interstellar pickup hydrogen on the solar wind termination shock
Hybrid (kinetic ions/fluid electrons) plasma simulations are used to study the effects of a population of energetic interstellar pickup hydrogen ions on the solar wind termination shock. The pickup hydrogen is treated as a second ion species in the simulations, and thus the effects of the pick-ups on the shock, as well as the effects of the shock on the pickups, are treated in a fully self-consistent manner. For quasi-perpendicular shocks with 10-20 percent pickup hydrogen the pickup ions manifest themselves in a small foot ahead of the shock ramp caused by pickup ion reflection. For oblique shocks with smaller angles between the field and the shock normal, a large fraction of the pickup ions are reflected and move back upstream where they excite large amplitude magnetosonic waves which steepen into shocklets. These backstreaming pickup ions may provide advance warning of a spacecraft encounter with the termination shock.
Liewer, P. C.↗
Two dimensional PIC simulations of plasma heating by the dissipation of Alfven waves
Two dimensional plasma particle simulations of the evolution of large amplitude circularly polarized Alfven waves propagating parallel to the magnetic field show that the waves decay via both one- and two- dimensional parametric decay instabilities. For parameters studied, one-dimensional processes dominate the simulations, but two-dimensional decay processes, including the recently predicted filamentation instability are also observed. The daughter waves generated by the parametric decay are primarily damped by the ions, leading to ion heating. The parametric decay processes efficiently convert the ordered fluid ion motion in the Alfven wave into ion thermal energy. These processes may be important for the dissipation of Alfven waves in the solar wind, the corona and other space plasma environments. The computations were performed on the Intel Touchstone parallel supercomputer.
Liewer, P. C.↗
Numerical studies of electron dynamics in oblique quasi-perpendicular collisionless shock waves
Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.
Liewer, P. C.↗
Hypercube-Computer Analysis Of Electromagnetic Scattering
Capabilities of hypercube and parallel processing demonstrated. Report describes use of Mark III Hypercube computer to analyze scattering of electromagnetic waves. Purpose of study to assess utility of parallel computing in such computation-intensive problems as large-scale electromagnetic scattering. Two electromagnetic codes based on different algorithms converted to run on Mark III Hypercube. First code implements finite-difference, time-domain solution of Maxwell's curl equations. Second code is Numerical Electromagnetics Code (NEC-2) which embodies frequency-domain method and developed to analyze electromagnetic responses of antennas and other metallic structures. On Mark III Hypercube with 32 active nodes, largest lattice contains about 2,048,000 unit cells.
Patterson, J. E.↗