Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC simulation”

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 55 records · Page 3

Studies on equatorial shock formation during plasmaspheric refilling

During the grant period starting August 1, 1992, our major effort has been on examining the presence of equatorially trapped hot plasma on plasmaspheric refilling. We performed one-dimensional PIC simulations of cold plasmas expanding into a hot plasma, consisting of hot anisotropic ions and warm isotropic electrons, trapped in a region of minimum magnetic field. Simulations showed that the electric potential barrier built up by the anisotropy of the hot ion population facilitates in the formation of electrostatic shocks when the cold ion beams begin to come into contact with the hot plasma. The shock formation occurs even when the cold ion beams are highly supersonic with respect to the ion-acoustic speed. This finding is interesting because equatorial shock formation during the early stage of plasmaspheric refilling has been debated over about two decades. In the past ion-ion instability has been invoked as the main mechanism for the coupling between the cold ion beams approaching the equator from the conjugate ionspheres. This coupling occurs when the beams are sufficiently slow; the beam velocity being less than three times the ion-acoustic speed. In the presence of hot plasma, the beams slow down by the potential barrier. The slowing down and the reflection process lead to the formation of the electrostatic shock even for highly supersonic ion beams. The mixing of hot and cold plasma was also studied.

Singh, Nagendra↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the POLAR satellite is studied by using a 3-dimensional Particle-In-Cell (PIC) code. When the satellite passes through the region of low density plasma, the satellite charges to positive potentials as high as 4050Volts due to the photoelectrons emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, an ion-rich Xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at about 2Volts. Accordingly, in our 3-dimensional PIC simulation, we considered that the potential of the satellite is 2Volts as a fixed bias. Considering the relatively high density of the Xenon plasma in the sheath (approx. 10 - 10(exp 3)/cc), the ambient plasma of low density (less than 1/cc) is neglected. In the simulations, the electric fields and plasma dynamics are calculated self-consistently. We found that an "Apple" shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission, a high positive potential hill develops. Near the Thermal Ion Detection Experiment (TIDE) detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations, it takes about a hundred electron gyroperiods for the sheath to reach a quasi-steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. Using the steady state sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Leung, Wing C.↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗

Three-Dimensional Modeling of Guide-Field Magnetic Reconnection

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 analyze three-dimensional PIC simulations of guide-field magnetic reconnection. Specific emphasis will be on the question whether thermal-inertia processes, mediated by the electron pressure tensor, remain a viable dissipation mechanism in fully three-dimensional systems.

Hesse, Michael↗

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↗

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↗

The Dissipation Mechanism 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. 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 PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

Particle Acceleration, Magnetic Field Generation, and Associated Emission in Collisionless Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron)jets show that acceleration occurs within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

Particle Acceleration, Magnetic Field Generation and Associated Emission in Collisionless Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K. I.↗

The Electron Diffusion Region: Forces and Currents

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. 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 PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertia1 effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

Radiation from Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

What Turns Reconnection On and Off

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. 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 PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density. which matches key features of electron distribution functions.

Hesse, Michael↗

Microscopic Processes in Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

Radiation from Relativistic Shocks with Turbulent Magnetic Fields

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs at shocked region. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the shock. The "jitter" radiation from deflected electrons in turbulent magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. New recent calculation of spectra with various different Lorentz factors of jets (two electrons) and initial magnetic fields. New spectra based on small simulations will be presented.

Nishkawa, K.↗

Microscopic Processes On Radiation from Accelerated Particles in Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The jitter'' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

The Electron Diffusion Region: Forces and Currents

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. 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 PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

Numerical Studies of Hall Thruster Acceleration Region Electron Transport

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

Ortega, Alejandro Lopez↗

Reconnection in a Pinch

A recently published analysis of current sheets has updated the classic Harris 1D static solution by considering multiple classes of charged particle trajectories in a generalized and dynamic current sheet. It uses a 1D PIC simulation to describe dynamic pinching and bifurcation of the current sheet. These 1D results strongly suggest that properties of the inflowing plasma, including the plasma beta, have an important effect on the equilibrium thickness of the pinched current sheet. Since 1D studies cannot describe magnetic reconnection, the time appears right to carry such 1D studies over to 2D or 3D simulations to explore current sheet thickness effects on reconnection. The Magnetospheric Multiscale Mission (MMS), with its well-resolved multipoint measurements of collisionless plasma and fields, has found that collisionless reconnection is accompanied by non-adiabatic motions of electrons that only occur in magnetic structures with a narrow scale comparable to electron inertial lengths (d e ). The recent 1D studies suggest that a plasma pinch to such scales may only occur for inflowing magnetized plasmas with relatively low plasma beta. We conclude that a parametric exploration of simulated and observed reconnection inflow conditions, particularly plasma beta, should shed light on the enablement of reconnection in collisionless plasmas.

reconnection↗