Search NASA⌕ Search

SEARCH · Search NASA

Results for “test particle simulations”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Self-consistent current sheet structures in the quiet-time magnetotail

The structure of the quiet-time magnetotail is studied using a test particle simulation. Vlasov equilibria are obtained in the regime where v(D) = E(y) c/B(z) is much less than the ion thermal velocity and are self-consistent in that the current and magnetic field satisfy Ampere's law. Force balance between the plasma and magnetic field is satisfied everywhere. The global structure of the current sheet is found to be critically dependent on the source distribution function. The pressure tensor is nondiagonal in the current sheet with anisotropic temperature. A kinetic mechanism is proposed whereby changes in the source distribution results in a thinning of the current sheet.

Holland, Daniel L.↗
NASA NTRS1994

A proposed neutral line signature

An identifying signature is proposed for the existence and location of the neutral line in the magnetotail. The signature, abrupt density and temperature changes in the Earth-tail direction, was first discovered in test particle simulations. Such temperature variations have been observed in International Sun Earth Explorer (ISEE) data (Huang et al., 1992), but their connection to the possible existence of a neutral line in the tail has not yet been established. The proposed signature develops earlier than the ion velocity space ridge of Martin and Speiser (1988) but can only be seen by spacecraft in the vicinity of the neutral line, while the latter can locate a neutral line remotely.

Doxas, I.↗
NASA NTRS1996

Excitation of Plasma Waves in Aurora by Electron Beams

In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

daSilva, C. E.↗
NASA NTRS1998

Dynamics of the Trans-Neptune Region: Apsidal Waves in the Kuiper Belt

The role of apsidal density waves propagating in a primordial trans-Neptune disk (i.e., Kuiper belt) is investigated. It is shown that Neptune launches apsidal waves at its secular resonance near 40 AU that propagate radially outward, deeper into the particle disk. The wavelength of apsidal waves is considerably longer than waves that might be launched at Lindblad resonances, because the pattern speed, g(sub s), resulting from the apsis precession of Neptune is much slower than its mean motion, Omega(sub s). If the early Kuiper belt had a sufficient surface density, sigma, the disk's wave response to Neptune's secular perturbation would have spread the disturbing torque radially over a collective scale lambda(sub *) approx. = r(2(mu)(sub d)Omega/ absolute value of r dg/dr)(sup 1/2), where mu(sub d)equivalent pi(sigma)r(exp 2)/(1 solar mass) and Omega(r) and g(r) are respectively the mean motion and precession frequency of the disk particles. This results in considerably smaller eccentricities at resonance than had the disk particles been treated as noninteracting test particles. Consequently, particles are less apt to be excited into planet-crossing orbits, implying that the erosion timescales reported by earlier test-particle simulations of the Kuiper belt may be underestimated. It is also shown that the torque the disk exerts upon the planet (due to its gravitational attraction for the disk's spiral wave pattern) damps the planet's eccentricity and further inhibits the planet's ability to erode the disk. Key words: celestial mechanics, stellar dynamics - comets: general minor planets, asteroids

Ward, William R.↗
NASA NTRS2011

Large-Amplitude Transmitter-Associated and Lightning-Associated Whistler Waves in the Earth's Inner Plasmasphere at L less than 2

We report observations of very large amplitude whistler mode waves in the Earth fs nightside inner radiation belt enabled by the STEREO Time Domain Sampler. Amplitudes range from 30.110 mV/m (zero ]peak), 2 to 3 orders of magnitude larger than previously observed in this region. Measurements from the peak electric field detector (TDSMax) indicate that these large ]amplitude waves are prevalent throughout the plasmasphere. A detailed examination of high time resolution electric field waveforms is undertaken on a subset of these whistlers at L < 2, associated with pump waves from lightning flashes and the naval transmitter NPM in Hawaii, that become unstable after propagation through the ionosphere and grow to large amplitudes. Many of the waveforms undergo periodic polarization reversals near the lower hybrid and NPM naval transmitter frequencies. The reversals may be related to finite plasma temperature and gradients in density induced by ion cyclotron heating of the plasma at 200 Hz, the modulation frequency of the continuous ]mode NPM naval transmitter signal. Test particle simulations using the amplitudes and durations of the waves observed herein suggest that they can interact strongly with high ]energy (>100 keV) electrons on a time scale of <1 s and thus may be an important previously unaccounted for source of energization or pitch ]angle scattering in the inner radiation belt.

Breneman, A.↗
NASA NTRS2013

Energetic Particles Dynamics in Mercury's Magnetosphere

We investigate the drift paths of energetic particles in Mercury's magnetosphere by tracing their motion through a model magnetic field. Test particle simulations solving the full Lorentz force show a quasi-trapped energetic particle population that gradient and curvature drift around the planet via "Shabansky" orbits, passing though high latitudes in the compressed dayside by equatorial latitudes on the nightside. Due to their large gyroradii, energetic H+ and Na+ ions will typically collide with the planet or the magnetopause and will not be able to complete a full drift orbit. These simulations provide direct comparison for recent spacecraft measurements from MESSENGER. Mercury's offset dipole results in an asymmetric loss cone and therefore an asymmetry in particle precipitation with more particles precipitating in the southern hemisphere. Since the planet lacks an atmosphere, precipitating particles will collide directly with the surface of the planet. The incident charged particles can kick up neutrals from the surface and have implications for the formation of the exosphere and weathering of the surface

Energetic particles↗
NASA NTRS2016

Outer Radiation Belt Dropout Dynamics Following the Arrival of Two Interplanetary Coronal Mass Ejections

Magnetopause shadowing and wave-particle interactions are recognized as the two primary mechanisms for losses of electrons from the outer radiation belt. We investigate these mechanisms, sing satellite observations both in interplanetary space and within the magnetosphere and particle drift modeling. Two interplanetary shocks sheaths impinged upon the magnetopause causing a relativistic electron flux dropout. The magnetic cloud (C) and interplanetary structure sunward of the MC had primarily northward magnetic field, perhaps leading to a concomitant lack of substorm activity and a 10 day long quiescent period. The arrival of two shocks caused an unusual electron flux dropout. Test-particle simulations have shown 2 to 5 MeV energy, equatorially mirroring electrons with initial values of L 5.5can be lost to the magnetosheath via magnetopause shadowing alone. For electron losses at lower L-shells, coherent chorus wave-driven pitch angle scattering and ULF wave-driven radial transport have been shownto be viable mechanisms.

Alves, L. R.↗
NASA NTRS2017

Transverse eV Ion Heating by Random Electric Field Fluctuations in the Plasmasphere

Charged particle acceleration in the Earth inner magnetosphere is believed to be mainly due to the local resonant wave-particle interaction or particle transport processes. However, the Van Allen Probes have recently provided interesting evidence of a relatively slow transverse heating of eV ions at distances about 2-3 Earth radii during quiet times. Waves that are able to resonantly interact with such very cold ions are generally rare in this region of space, called the plasmasphere. Thus, non-resonant wave-particle interactions are expected to play an important role in the observed ion heating. We demonstrate that stochastic heating by random transverse electric field fluctuations of whistler (and possibly electromagnetic ion cyclotron) waves could explain this weak and slow transverse heating of H+ and O+ ions in the inner magnetosphere. The essential element of the proposed model of ion heating is the presence of trains of random whistler (hiss) wave packets, with significant amplitude modulations produced by strong wave damping, rapid wave growth, or a superposition of wave packets of different frequencies, phases, and amplitudes. Such characteristics correspond to measured characteristics of hiss waves in this region. Using test particle simulations with typical wave and plasma parameters, we demonstrate that the corresponding stochastic transverse ion heating reaches 0.07-0.2 eV/h for protons and 0.007-0.015 eV/h for O+ ions. This global temperature increase of the Maxwellian ion population from an initial Ti approx. 0.3 eV could potentially explain the observations.

Artemyev, A. V.↗
NASA NTRS2017

Impact of Substorm Time O+ Outflow on Ring Current Enhancement

Energetic O+ ions (tens of keV) rapidly increase in the inner magnetosphere and contribute significantly to the ring current during substorms. Previously, two source regions of the energetic O+ ions have been proposed. The first one is the dayside polar region. Ions from the dayside polar region are transported to the lobe; then they are injected to the nightside plasma sheet during substorm expansion phase. The second one is the nightside aurora region. After the substorm onset, energetic O+ ions are extracted from the ionosphere with the auroral acceleration processes, and the O+ ions are directly supplied to the nightside plasma sheet. We investigated the relative importance of these two regions on supplying the energetic O+ ions in the inner magnetosphere. We performed a test particle simulation in global MHD electromagnetic fields. We obtained the following results. (1) During the substorm growth phase, O+ ions at tens of eV are extracted from the dayside polar region, resulting in the enhancement of the warm O+ ions (hundreds of eV) in the lobe. After the substorm onset, the warm O+ ions are nonadiabatically accelerated to tens of keV and injected to the inner magnetosphere. These O+ ions contribute to most of the O+ ring current. (2) O+ ions at less than a few keVs are supplied from the nightside aurora region to the plasma sheet. However, their contribution to the O+ ring current remains small. From the results, we concluded that the main source of the energetic O+ ions is the dayside polar region.

outflow↗
NASA NTRS1986

Acceleration of charged particles in magnetic reconnection Solar flares, the magnetosphere, and solar wind

A possible source of free energy available for accelerating charged particles is conversion of magnetic energy to particle energy in reconnecting magnetic fields. Recent simulations using test particles suggests that reconnection may efficiently accelerate particles to the maximum energies that are observed in several astrophysical contexts. A simple analytic formula is used in conjunction with the simulation results to predict the maximum energy achievable in a particular plasma environment with the result that in solar flares reconnection is capable of accelerating particles to several GeV. In magnetospheric substorms the predicted maximum can reach several hundred keV, and near magnetic sector crossings in the solar wind the maximum energy can approach 100 keV.

Goldstein, M. L.↗
NASA NTRS1978

Particle impact tests

Particle impact tests were performed on three types of orbiter surface with a micrometeoroid facility. The test equipment electrostatically accelerated micron sized particles to high velocities simulating micrometeoroid impacts. Test particles were titanium diboride with typical velocities in the range 1 to 2.3 km x sec/1 and equivalent particle diameters in the range 4 to 16 microns. Impact angles to the material surface were 90, 60 and 30 degrees. The particle impact sites were located on the sample surfaces and craters were photographed with a magnification of 400X.

Komatsu, G. K.↗
DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Sustained Oxygen Spectral Gaps and Their Dynamic Evolution in the Inner Magnetosphere

Van Allen Probes observations of ion spectra often show a sustained gap within a very narrow energy range throughout the full orbit. To understand their formation mechanism, we statistically investigate the characteristics of the narrow gaps for oxygen ions and find that they are most frequently observed near the noon sector with a peak occurrence rate of over 30%. The magnetic moment (μ) of the oxygen ions in the gap shows a strong dependence on magnetic local time (MLT), with higher and lower μ values in the morning and afternoon sectors, respectively. Moreover, we find through superposed epoch analysis that the gap formation also depends on geomagnetic conditions. Those gaps formed at lower magnetic moments (μ < 3,000 keV/G) are associated with stable convection electric fields, which enable magnetospheric ions to follow a steady drift pattern that facilitates the gap formation by corotational drift resonance. On the other hand, gaps with higher μ values are statistically preceded by a gradual increase of geomagnetic activity. Here, we suggest that ions within the gap were originally located inside the Alfven layer following closed drift paths, before they were transitioned into open drift paths as the convection electric field was enhanced. The sunward drift of these ions, with very low fluxes, forms a drainage void in the dayside magnetosphere manifested as the sustained gap in the oxygen spectrum. This scenario is supported by particle-tracing simulations, which reproduce most of the observed characteristics and therefore provide new insights into inner magnetospheric dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Monte Carlo Analysis of the Performance of the ITER Diagnostic Residual Gas Analyzer

The composition of exhausted gas is a key parameter in long-pulse plasma fusion experiments, and its evolution shall be monitored at timescales relevant to plasma dynamics and plasma-wall interactions. A diagnostic residual gas analyzer (DRGA) is a multisensor instrument particularly suited to these studies, and ITER will adopt DRGAs in the equatorial and in the divertor tokamak regions. In this work, we have revisited the design of the ITER divertor DRGA through simple vacuum analytical considerations supported by simulations conducted with Molflow+, a test particle Monte Carlo (TPMC) simulation code commonly used in the particle accelerator community. Here, starting with recommendations on the manufacturing of the vacuum piping of the DRGA, this work is followed by a complete vacuum characterization of the diagnostic vacuum setup (pressure profiles at base pressure and during sampling, orifice diameter, and length optimization), and finally, the in-vessel residence time of the most important gas species is simulated. These studies have allowed us to give insights into some experimental results recently found on the prototype DRGA installed in the Wendelstein W7-X stellarator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗