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 181 records · Page 10

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

Nongyrotropic Electrons in Guide Field Reconnection

We apply a scalar measure of nongyrotropy to the electron pressure tensor in a 2D particle-in-cell simulation of guide field reconnection and assess the corresponding electron distributions and the forces that account for the nongyrotropy. The scalar measure reveals that the nongyrotropy lies in bands that straddle the electron diffusion region and the separatrices, in the same regions where there are parallel electric fields. Analysis of electron distributions and fields shows that the nongyrotropy along the inflow and outflow separatrices emerges as a result of multiple populations of electrons influenced differently by large and small-scale parallel electric fields and by gradients in the electric field. The relevant parallel electric fields include large-scale potential ramps emanating from the x-line and sub-ion inertial scale bipolar electron holes. Gradients in the perpendicular electric field modify electrons differently depending on their phase, thus producing nongyrotropy. Magnetic flux violation occurs along portions of the separatrices that coincide with the parallel electric fields. An inductive electric field in the electron EB drift frame thus develops, which has the effect of enhancing nongyrotropies already produced by other mechanisms and under certain conditions producing their own nongyrotropy. Particle tracing of electrons from nongyrotropic populations along the inflows and outflows shows that the striated structure of nongyrotropy corresponds to electrons arriving from different source regions. We also show that the relevant parallel electric fields receive important contributions not only from the nongyrotropic portion of the electron pressure tensor but from electron spatial and temporal inertial terms as well.

Parallel Electric Fields↗

Two-Scale Ion Meandering Caused by the Polarization Electric Field During Asymmetric Reconnection

Ion velocity distribution functions (VDFs) from a particle-in-cell simulation of asymmetric reconnection are investigated to reveal a two-scale structure of the ion diffusion region (IDR). Ions bouncing in the inner IDR are trapped mainly by the electric field normal to the current sheet (N direction), while those reaching the outer IDR are turned back mainly by the magnetic force. The resulting inner layer VDFs have counter-streaming populations along N with decreasing counter-streaming speeds away from the midplane while maintaining the out-of-plane speed, and the outer layer VDFs exhibit crescent shapes toward the out-of-plane direction. Observations of the above VDF features and the normal electric fields provide evidence for the two-scale meandering motion.

Wang, Shan↗

Fusion-Enabled Pluto Orbiter and Lander

The Pluto orbiter mission proposed here is credible and exciting. The benefits to this and all outer-planet and interstellar-probe missions are difficult to overstate. The enabling technology, Direct Fusion Drive, is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor under development at the Princeton Plasma Physics Laboratory. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. During this Phase I effort, we made great strides in modeling the engine efficiency, thrust, and specific impulse and analyzing feasible trajectories. Based on 2D fluid modeling of the fusion reactors outer stratum, its scrape-off-layer (SOL), we estimate achieving 2.5 to 5 N of thrust for each megawatt of fusion power, reaching a specific impulse, Isp, of about 10,000 s. Supporting this model are particle-in-cell calculations of energy transfer from the fusion products to the SOL electrons. Subsequently, this energy is transferred to the ions as they expand through the magnetic nozzle and beyond. Our point solution for the Pluto mission now delivers 1000 kg of payload to Pluto orbit in 3.75 years using 7.5 N constant thrust. This could potentially be achieved with a single 1 MW engine. The departure spiral from Earth orbit and insertion spiral to Pluto orbit require only a small portion of the total delta-V. Departing from low Earth orbit reduces mission cost while increasing available mission mass. The payload includes a lander, which utilizes a standard green propellant engine for the landing sequence. The lander has about 4 square meters of solar panels mounted on a gimbal that allows it to track the orbiter, which beams 30 to 50 kW of power using a 1080 nm laser. Optical communication provides dramatically high data rates back to Earth. Our mass modeling investigations revealed that if current high-temperature superconductors are utilized at liquid nitrogen temperatures, they drive the mass of the engine, partly because of the shielding required to maintain their critical temperature. Second generation materials are thinner but the superconductor is a very thin layer deposited on a substrate with additional layers of metallic classing. Tremendous research is being performed on a variety of these superconducting materials, and new irradiation data is now available. This raises the possibility of operating nearfuture high-temperature superconductors at a moderately low temperature to dramatically reduce the amount of shielding required. At the same time, a first generation space engine may require low-temperature superconductors, which are higher TRL and have been designed for space coils before (AMS-02 experiment for the ISS). We performed detailed analysis of the startup system and thermal conversion system components. The ideal working fluid was determined to be a blend of Helium and Xenon. No significant problems were identified with these subsystems. For the RF system, we conceived of a new, more efficient design using state-of-the-art switch amplifiers, which have the potential for 100% efficiency. This report presents details of our engine and trajectory analyses, mass modeling efforts, and updated vehicle designs.

temperature↗

Ion-Scale Secondary Flux Ropes Generated by Magnetopause Reconnection as Resolved by MMS

New Magnetospheric Multiscale (MMS) observations of small-scale (approx. 7 ion inertial length radius) flux transfer events (FTEs) at the dayside magnetopause are reported. The 1O km MMS tetrahedron size enables their structure and properties to be calculated using a variety of multispacecraft techniques, allowing them to be identified as flux ropes, whose flux content is small (approx. 22 kWb).The current density, calculated using plasma and magnetic field measurements independently, is found to be filamentary. lntercomparison of the plasma moments with electric and magnetic field measurements reveals structured non-frozen-in ion behavior. The data are further compared with a particle-in-cell simulation. It is concluded that these small-scale flux ropes, which are not seen to be growing, represent a distinct class of FTE which is generated on the magnetopause by secondary reconnection.

Eastwood, J. P.↗

Temporal Evolution of the Plasma Sheath Surrounding Solar Cells in Low Earth Orbit

High voltage solar array interactions with the space environment can have a significant impact on array performance and spacecraft charging. Over the past 10 years, data from the International Space Station has allowed for detailed observations of these interactions over long periods of time. Some of the surprising observations have been floating potential transients, which were not expected and are not reproduced by existing models. In order to understand the underlying processes producing these transients, the temporal evolution of the plasma sheath surrounding the solar cells in low Earth orbit is being investigated. This study includes lumped element modeling and particle-in-cell simulation methods. This presentation will focus on recent results from the on-going investigations.

Willis, Emily M.↗

Temporal Evolution of the Plasma Sheath Surrounding Solar Cells in Low Earth Orbit

High voltage solar array interactions with the space environment can have a significant impact on array performance and spacecraft charging. Over the past 10 years, data from the International Space Station has allowed for detailed observations of these interactions over long periods of time. Some of the surprising observations have been floating potential transients, which were not expected and are not reproduced by existing models. In order to understand the underlying processes producing these transients, the temporal evolution of the plasma sheath surrounding the solar cells in low Earth orbit is being investigated. This study includes lumped element modeling and particle-in-cell simulation methods. This presentation will focus on recent results from the on-going investigations.

Willis, Emily 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.↗

3D Simulations of Ion Thruster Accelerator Grid Erosion Accounting for Charge Exchange Ion Space Charge

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

Chaplin, Vernon H.↗

Generation of Electron Whistler Waves at the Mirror Mode Magnetic Holes: MMS Observations and PIC Simulation

The Magnetospheric Multiscale mission has observed electron whistler waves at the center and at the edges of magnetic holes in the dayside magnetosheath. The magnetic holes are nonlinear mirror structures since their magnitude is anticorrelated with particle density. In this article, we examine the growth mechanisms of these whistler waves and their interaction with the host magnetic hole. In the observations, as magnetic holes develop and get deeper, an electron population gets trapped and develops a temperature anisotropy favorable for whistler waves to be generated. In addition, the decrease in magnetic field magnitude and the increase in density reduce the electron resonance energy, which promotes the electron cyclotron resonance. To investigate this process, we used expanding box particle-in-cell simulations to produce the mirror instability, which then evolve into magnetic holes. The simulation shows that whistler waves can be generated at the center and edges of magnetic holes, which reproduces the primary features of the MMS observations. The simulation shows that the electron temperature anisotropy develops in the center of the magnetic hole once the mirror instability reaches its nonlinear stage of evolution. The plasma is then unstable to whistler waves at the minimum of the magnetic field structures. In the saturation regime of mirror instability, when magnetic holes are developed, the electron temperature anisotropy appears at the edges of the holes and electron distributions become more isotropic at the magnetic field minimum. At the edges, the expansion of magnetic holes decelerates the electrons, which leads to temperature anisotropies.

Magnetospheric Multiscale↗

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↗

Energetic Particles of Cosmic Accelerators II: Active Galactic Nuclei and Gamma-ray Bursts

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Our pursuit of more than a century to uncover the origins and fate of these cosmic energetic particles has given rise to some of the most interesting and challenging questions in astrophysics. Within our own galaxy, we have seen that energetic particles engage in a complex interplay with the galactic environment and even drive many of its key characteristics (for more information, see the first white paper in this series). On cosmological scales, the energetic particles supplied by the jets of active galactic nuclei (AGN) are an important source of energy for the intracluster and intergalactic media, providing a mechanism for regulating star formation and black hole growth and cultivating galaxy evolution (AGN feedback). Gamma-ray burst (GRB) afterglows encode information about their circumburst environment, which has implications for massive stellar winds during previous epochs over the stellar lifecycle. As such, GRB afterglows provide a means for studying very high-redshift galaxies since GRBs can be detected even if their host galaxy cannot. It has even been suggest that GRB could be used to measure cosmological distance scales if they could be shown to be standard candles. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. Despite the enormous progress in particle-in-cell and magnetohydrodynamic simulations, we have yet to pinpoint the processes involved in jet formation and collimation and the conditions under which they can occur. For that matter, we have yet to identify the mechanism(s) through which the jet accelerates energetic particles – is it the commonly invoked diffusive shock acceleration process or is another mechanism, such as magnetic reconnection, required? Do AGNs and GRBs accelerate hadrons, and if so, do they accelerate them to ultra-high energies and are there high-energy neutrinos associated with them? MeV gamma-ray astronomy, enabled by technological advances that will be realized in the coming decade, will provide a unique and indispensable perspective on the persistent mysteries of the energetic universe. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. Specifically, MeV astronomy will: 1. Determine whether AGNs accelerate CRs to ultra-high energies; 2. Provide the missing pieces for the physics of the GRB prompt emission; 3. Measure magnetization in cosmic accelerators and search for acceleration via reconnection.

Venters, Tonia M.↗

Energy Conversion and Partition in the Asymmetric Reconnection Diffusion Region

We investigate the energy conversion and partition in the asymmetric reconnection diffusion region using two-dimensional particle-in-cell simulations and Magnetosphere Multiscale (MMS) mission observations. Under an upstream condition with equal temperatures in the two inflow regions, the simulation analysis indicates that the energy partition between ions and electrons depends on the distance from the X-line. Within the central electron diffusion region (EDR), nearly all dissipated electromagnetic field energies are converted to electrons. From the EDR to the ion diffusion region (IDR) scales, the rate of the electron energy gain decreases to be lower than that of ions. A magnetopause reconnection event inside the IDR observed by MMS shows comparable ion and electron energy gains, consistent with the simulation result in the transition region from EDR to IDR. At the EDR scale, the electron energization is mainly by the reconnection electric field (E(sub r)); in-plane electric fields (E(sub xz)) provide additional positive contributions near the X-line and do negative work on electrons beyond the EDR. The guide field reduces the electron energization by both E(sub r) and E(sub xz) in the EDR. For ion energization, E(sub r) and E(sub xz) have comparable contributions near the time of the peak reconnection rate, while E(sub xz) dominants at later time. At the IDR scale, the guide field causes asymmetry in the amount of the energy gain and energization mechanisms between two exhausts but does not have significant effects on energy partition. Our study advances understanding of ion and electron energization in asymmetric reconnect IDRs.

Wang, Shan↗

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth's bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (di), and the other has a wavelength less than 1 di. Electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

N. Bessho↗

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

Plasmaspheric Plume Turbulence: Signature of an Electrostatic Corotation Convection Shear-Layer Instability

As demonstrated by IMAGE-EUV observations, plasmaspheric drainage plumes are a common feature emanating from the plasmasphere during periods of enhanced convection. Concurrent IMAGE-EUV and LANL-MPA observations reveal that plasmaspheric plumes are not a uniform density structure that extends out to geosynchronous orbit and beyond. Instead, in–situ observations reveal a large spatial-scale structure populated by small-scale density structures suggestive of the presence of a turbulent process. Here, we investigate the occurrence and nature of plasmaspheric drainage plumes as observed by concurrent EUV-MPA observations, as well as explore the possibility that the in–situ observed small-scale density structures are the signature of an instability produced by the sheared velocity found within the plasmaspheric layers separating the corotation of the main plasmasphere from the convection-driven flow generating the plume. Particle-in-cell simulations indicate that the generation of shear-flow driven instabilities is possible under plasmaspheric plasma conditions. The initial results of these simulations indicate that the plasmaspheric instabilities generated are electrostatic in nature, consist with in–situ plume observations.

Mark L Adrian↗

Modeling the Lunar Wake Response to a CME Using a Hybrid PIC Model

In the solar wind, a low-density wake region forms downstream of the nightside lunar surface. In this study, we use a series of 3D hybrid particle-in-cell simulations to model the response of the lunar wake to a passing coronal mass ejection (CME). Average plasma parameters are derived from the Wind spacecraft located at 1 au during three distinct phases of a passing halo (Earth-directed) CME on 2015 June 22. Each set of plasma parameters, representing the shock/plasma sheath, a magnetic cloud, and plasma conditions we call the mid-CME phase, are used as the time-static upstream boundary conditions for three separate simulations. These simulation results are then compared with results that use nominal solar wind conditions. Results show a shortened plasma void compared to nominal conditions and a distinctive rarefaction cone originating from the terminator during the CME's plasma sheath phase, while a highly elongated plasma void reforms during the magnetic cloud and mid-CME phases. Developments of electric and magnetic field intensification are also observed during the plasma sheath phase along the central wake, while electrostatic turbulence dominates along the plasma void boundaries and 2–3 lunar radii RM downstream in the central wake during the magnetic cloud and mid-CME phases. The simulations demonstrate that the lunar wake responds in a dynamic way with the changes in the upstream solar wind during a CME.

Anthony P. Rasca↗