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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Kinetic Simulations of Current-Sheet Formation and Reconnection at a Magnetic X Line

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.

Black, C.↗
NASA NTRS2011

The Inner Structure of Collisionless Magnetic Reconnection: The Electron-Frame Dissipation Measure and Hall Fields

It was recently proposed that the electron-frame dissipation measure, the energy transfer from the electromagnetic field to plasmas in the electron s rest frame, identifies the dissipation region of collisionless magnetic reconnection [Zenitani et al., Phys. Rev. Lett. 106, 195003 (2011)]. The measure is further applied to the electron-scale structures of antiparallel reconnection, by using two-dimensional particle-in-cell simulations. The size of the central dissipation region is controlled by the electron-ion mass ratio, suggesting that electron physics is essential. A narrow electron jet extends along the outflow direction until it reaches an electron shock. The jet region appears to be anti-dissipative. At the shock, electron heating is relevant to a magnetic cavity signature. The results are summarized to a unified picture of the single dissipation region in a Hall magnetic geometry.

Zenitani, Seiji↗
NASA NTRS2012

New Measure of the Dissipation Region in Collisionless Magnetic Reconnection

A new measure to identify a small-scale dissipation region in collisionless magnetic reconnection is proposed. The energy transfer from the electromagnetic field to plasmas in the electron s rest frame is formulated as a Lorentz-invariant scalar quantity. The measure is tested by two-dimensional particle-in-cell simulations in typical configurations: symmetric and asymmetric reconnection, with and without the guide field. The innermost region surrounding the reconnection site is accurately located in all cases. We further discuss implications for nonideal MHD dissipation.

Zenitani, Seiji↗
NASA NTRS2011

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

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
NASA NTRS2014

Direct Energy Conversion for Nuclear Propulsion at Low Specific Mass

The project will continue the FY13 JSC IR&D (October-2012 to September-2013) effort in Travelling Wave Direct Energy Conversion (TWDEC) in order to demonstrate its potential as the core of a high potential, game-changing, in-space propulsion technology. The TWDEC concept converts particle beam energy into radio frequency (RF) alternating current electrical power, such as can be used to heat the propellant in a plasma thruster. In a more advanced concept (explored in the Phase 1 NIAC project), the TWDEC could also be utilized to condition the particle beam such that it may transfer directed kinetic energy to a target propellant plasma for the purpose of increasing thrust and optimizing the specific impulse. The overall scope of the FY13 first-year effort was to build on both the 2012 Phase 1 NIAC research and the analysis and test results produced by Japanese researchers over the past twenty years to assess the potential for spacecraft propulsion applications. The primary objective of the FY13 effort was to create particle-in-cell computer simulations of a TWDEC. Other objectives included construction of a breadboard TWDEC test article, preliminary test calibration of the simulations, and construction of first order power system models to feed into mission architecture analyses with COPERNICUS tools. Due to funding cuts resulting from the FY13 sequestration, only the computer simulations and assembly of the breadboard test article were completed. The simulations, however, are of unprecedented flexibility and precision and were presented at the 2013 AIAA Joint Propulsion Conference. Also, the assembled test article will provide an ion current density two orders of magnitude above that available in previous Japanese experiments, thus enabling the first direct measurements of power generation from a TWDEC for FY14. The proposed FY14 effort will use the test article for experimental validation of the computer simulations and thus complete to a greater fidelity the mission analysis products originally conceived for FY13.

Scott, John H.↗
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS2014

Bidirectional Energy Cascades and the Origin of Kinetic Alfvenic and Whistler Turbulence in the Solar Wind

The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore how the free energy released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfvénic and whistler turbulence are excited that evolve through inverse and forward magnetic energy cascades.

plasmas↗
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 NTRS2013

Aspects of Collisionless Magnetic Reconnection in Asymmetric Systems

Asymmetric reconnection is being investigated by means of particle-in-cell simulations. The research has two foci: The direction of the reconnection line in configurations with nonvanishing magnetic fields; and the question why reconnection can be faster if a guide field is added to an otherwise unchanged asymmetric configuration. We find that reconnection prefers a direction, which maximizes the available magnetic energy, and show that this direction coincides with the bisection of the angle between the asymptotic magnetic fields. Regarding the difference in reconnection rates between planar and guide field models, we demonstrate that a guide field can provide essential confinement for particles in the reconnection region, which the weaker magnetic field in one of the inflow directions cannot necessarily provide.

Hesse, Michael↗
NASA NTRS2013

Aspects of Collisionless Magnetic Reconnection in Asymmetric Systems

Asymmetric reconnection is being investigated by means of particle-in-cell simulations. The research has two foci: the direction of the reconnection line in configurations with non-vanishing magnetic fields; and the question why reconnection can be faster if a guide field is added to an otherwise unchanged asymmetric configuration. We find that reconnection prefers a direction, which maximizes the available magnetic energy, and show that this direction coincides with the bisection of the angle between the asymptotic magnetic fields. Regarding the difference in reconnection rates between planar and guide field models, we demonstrate that a guide field can provide essential confinement for particles in the reconnection region, which the weaker magnetic field in one of the inflow directions cannot necessarily provide.

Magnetic Reconnection↗
NASA NTRS2014

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗
NASA NTRS2012

Multiscale Studies of the Three-Dimensional Dayside X-Line

We review recent experience from the Cluster, Double Star, and THEMIS missions for lessons that apply to the upcoming Magnetospheric Multiscale Mission (MMS) being developed for launch in 2014. On global scales, simulation and statistical studies lead to mean congurations of dayside reconnection, implying specific relative alignments of the inflow magnetic fields and X-line, with implications for MMS operations designed to maximize the number of close encounters with the diffusion region. At intermediate MHD-to-ion scales, reconstruction of features created by one or two X-lines have developed to the point where data from a cluster of spacecraft can determine their temporal trends and the approximate three-dimensional X-line structure. Recent petascale particle-in-cell (PIC) simulations of reconnection encompass three spatial dimensions with excellent resolution, and make striking predictions of electron scale physics that creates complex interacting flux ropes under component reconnection. High time resolution measurements from MMS will determine the detailed electron scale kinetics embedded within the global and MHDion scale contexts. These developments will lead to the refinement of our three-dimensional multiscale picture of reconnection, yielding improved understanding of the global, MHD, and local physics controlling the onset or quenching, variability, and mean rate of reconnection. This in turn will enable improved predictability of the structural features created by transient reconnection, and their space weather consequences.

Reconnection↗
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.↗
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 Particl-in-Cell (PIC) codes.

particle-in-cell↗