Search NASA⌕ Search

SEARCH · Search NASA

Results for “TEST PARTICLES”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Time-dependent evolution of cosmic-ray-modified shock structure: Transition to steady state

Steady state solutions to the two-fluid equations of cosmic-ray-modified shock structure were investigated first by Drury and Volk (1981). Their analysis revealed, among other properties, that there exist regions of upstream parameter space where the equations possess three different downstream solutions for a given upstream state. In this paper we investigate whether or not all these solutions can occur as time-asymptotic states in a physically realistic evolution. To do this, we investigate the time-dependent evolution of the two-fluid cosmic-ray equations in going from a specified initial condition to a steady state. Our results indicate that the time-asymptotic solution is strictly single-valued, and it undergoes a transition from weakly to strongly cosmic-ray-modified at a critical value of the upstream cosmic ray energy density. The expansion of supernova remnant shocks is considered as an example, and it is shown that the strong to weak transition is in fact more likely. The third intermediate solution is shown to influence the time-dependent evolution of the shock, but it is not found to be a stable time-asymptotic state. Timescales for convergence to these states and their implications for the efficiency of shock acceleration are considered. We also investigate the effects of a recently introduced model for the injection of seed particles into the shock accelerated cosmic-ray population. The injection is found to result in a more strongly cosmic-ray-dominated shock, which supports our conclusion that for most classes of intermediate and strong cosmic-ray-modified shocks, the downstream cosmic-ray pressure component is at least as large as the thermal gas pressure, independent of the upstream state. As a result, cosmic rays almost always play a significant role in determining the shock structure and dissipation and they cannot be regarded as test particles.

Donohue, D. J.↗
NASA NTRS2020

3D Propagation of Relativistic Solar Protons through Interplanetary Space

Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the position of the observer and is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, resulting in a rollover in the spectrum.

Solar energetic particles↗
NASA NTRS2020

3D propagation of relativistic solar protons through interplanetary space

Context.Solar energetic particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce ground level enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation which is only parallel to the magnetic field lines of interplanetary space, that is, spatially 1D propagation. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the∼100 MeV–1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons using a 3D approach to assess the effect of the heliospheric current sheet(HCS) and drifts associated to the gradient and curvature of the Parker spiral. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used for comparison purposes with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings per particle are strongly influenced by 3D effects, with significant differences between periods of A+ and A−polarities. The decay time constant of 1 AU intensity profiles varies depending on the position of the observer and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, resulting in a rollover in the spectrum.

S. Dalla↗
NASA NTRS1991

Numerical exploration of planetary arc dynamics

After analytically characterizing the response of a test particle to the combined effects of corotation and Lindblad resonances, due to the satellites of Neptune, as two coupled dynamical systems whose coupling depends on the collective effects in the arc, a formula is obtained for the energy imparted by the Lindblad resonance. This energy confines the arc material around the corotation resonance radius. The fact that the gradient of the torque density around the arc must be negative confirms the constraints defined by Lin et al. (1987) on the relative positions of the corotation and the Lindblad resonances if the arc is to be stable. These results are tested with a direct numerical simulation for two different arc configurations.

Sicardy, Bruno↗
NASA NTRS1987

A scaling law for accretion zone sizes

Current theories of runaway planetary accretion require small random velocities of the accreted particles. Two body gravitational accretion cross sections which ignore tidal perturbations of the Sun are not valid for the slow encounters which occur at low relative velocities. Wetherill and Cox have studied accretion cross sections for rocky protoplanets orbiting at 1 AU. Using analytic methods based on Hill's lunar theory, one can scale these results for protoplanets that occupy the same fraction of their Hill sphere as does a rocky body at 1 AU. Generalization to bodies of different sizes is achieved here by numerical integrations of the three-body problem. Starting at initial positions far from the accreting body, test particles are allowed to encounter the body once, and the cross section is computed. A power law is found relating the cross section to the radius of the accreting body (of fixed mass).

Greenzweig, Yuval↗
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.↗
DOE OSTI2026

Improving ProtoDUNE pion cross-section measurements with NuGraph Michel-electron tagging

Understanding hadron-argon interactions is essential for precise neutrino energy reconstruction and final-state interaction modeling in liquid-argon time projection chamber (LArTPC) experiments such as DUNE. In particular, pion absorption and charge-exchange processes constitute significant sources of systematic uncertainty in neutrino oscillation measurements. ProtoDUNE-SP, a large-scale LArTPC prototype operated at the CERN Neutrino Platform and exposed to charged-particle test beams in the few-GeV range, enables direct measurements of these processes. This work focuses on the measurement of differential cross sections for pion absorption and charge exchange using the 2 GeV/c pion beam data from the ProtoDUNE-SP run. A key component of this analysis is the identification of Michel electrons from $\pi \rightarrow \mu \rightarrow e$ decay chains, which helps separate different interaction topologies and improves background rejection. Michel electron identification will also assist in reliably calibrating the electromagnetic response in ProtoDUNE-SP data and for the future DUNE detectors. In this analysis, we apply NuGraph to identify Michel electrons. NuGraph is a graph neural network that models detector hits as nodes connected by spatial and temporal edges for particle and topology classification in LArTPC detectors. We first benchmark NuGraph’s Michel electron classification performance using ICEBERG data, a small-scale LArTPC prototype used for DUNE electronics and reconstruction development, and then transfer the approach to ProtoDUNE-SP. This poster presents the analysis strategy, NuGraph-based classification studies, and discusses how these developments are expected to improve the pion cross-section measurement.

Razafinime, Soamasina Herilala [Cincinnati U.] (OR↗
DOE OSTI2023

Understanding hadronic interactions using advanced reconstruction techniques

ProtoDUNE-SP served as the prototype of the future Deep Underground Neutrino Experiment (DUNE). With a total liquid argon mass of 0.77 kt, it stood as the most extensive monolithic single-phase Liquid Argon Time Projection Chamber (LArTPC) ever constructed. Strategically located and operated at CERN, ProtoDUNE-SP benefited from a specialized charged-particle test beam, offering a momentum range of 0.3 - 7 GeV/c. This setup provided a unique opportunity to conduct in-depth studies on hadronic interactions within argon. Hadronic interactions play a fundamental role in neutrino physics, particularly in the intricacies of neutrino detection. We depend on analyzing the outcomes of these interactions to accurately determine both the flavor and energy of the interacting neutrinos. Pions, as one of the primary byproducts of neutrino interactions, are crucial for understanding the overall dynamics and kinematics of these processes. Accurately characterizing and understanding the interaction of pions with argon can significantly enhance the precision of neutrino simulations and measurements.Being a Liquid Argon Time Projection Chamber (LArTPC), ProtoDUNE-SP is distinguished by its ability to produce high-definition images of charged particles as they traverse through the detector's active volume. Yet, accurately reconstructing these particle interactions and determining their kinematic attributes remains a challenging task. The successes of deep learning in diverse domains, especially in image recognition, provide a promising approach for addressing this challenge in ProtoDUNE. This thesis discusses a measurement of the Piplus - argon inelastic cross section in the energy range of 2400 to 3000 MeV, using data taken by protoDUNE during the fall of 2018. Additionally, panoptic segmentation, a machine learning technique, is introduced and validated through its application in the reconstruction of the neutral pion rest mass. The aim is to showcase how these advanced methods can enhance the quality of event reconstruction.

Sarasty Segura, Carlos Eduardo↗
NASA NTRS1998

Radio to Gamma-Ray Emission from Shell-Type Supernova Remnants: Predictions from Non-Linear Shock Acceleration Models

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays, produced by diffusive shock acceleration in the environs of the remnant's expanding blast wave. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. The recently reported observation of TeV gamma-rays from SN1006 by the CANGAROO Collaboration, combined with the fact that several unidentified EGRET sources have been associated with known radio/optical/X-ray-emitting remnants, provides powerful motivation for studying gamma-ray emission from SNRs. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency considerations and impact photon intensities and spectral shapes at all energies, producing GeV/TeV intensity ratios that are quite different from test particle predictions.

Baring, Matthew G.↗
DOE OSTI2023

Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector

Measurements of electrons from 𝜈𝑒 interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectra is derived, and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50 MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Heating and acceleration of ions with Kappa distribution functions by low‐frequency Alfvén wave

Abstract Heating and acceleration of ions with Kappa distribution functions (with parameter ) in low‐beta plasmas, by a low‐frequency Alfvén wave, is investigated using test‐particle simulations, yielding interesting new results. As long as the Alfvén wave amplitude is sufficiently large, the computed net heating energy of ions becomes independent of the wave frequency and amplitude, always approaching the same value of . The eventual energy of ions is dictated only by the initial ion energy and the ratio of the magnetic field energy density to the plasma density. The heating effect of the Kappa ions increases with . During the heating process, the ions are picked up by the Alfvén wave and pitch angle scattered, forming a quasi‐isotropic spherical shell velocity distribution. The Kappa ions are accelerated in the parallel direction, reaching a bulk flow speed roughly equal to the local Alfvén speed. Higher ‐value in the initial Kappa distribution leads to faster saturation. The above results may explain certain features of the ion heating and acceleration in the solar wind and corona.

Li, Kehua↗
DOE OSTI2020

Quantifying the Effect of Magnetic Field Line Curvature Scattering on the Loss of Ring Current Ions

During geomagnetic storms, the ring current ions sometimes exhibit rapid loss as suggested by the fast recovery of the Dst index on a time scale of a few hours. Here, the effects of magnetic field line curvature (FLC) scattering on the loss of ring current ions, which have not been well quantified, are studied here by test particle simulations under the T89c magnetic field model. Our simulation results show that the prediction of ion loss based on a single–value cutoff of the κ parameter or maximum of δμ/μ of a single FLC scattering is not accurate. Instead, the e–folding lifetime (τ) for the loss of ring current ions due to cumulative FLC scattering has been calculated for different initial ion energies, equatorial pitch angles, and L shells under different geomagnetic conditions. The results show that in general the FLC scattering loss is faster for ions of higher energy, higher mass, smaller pitch angle, higher L, and at high Kp level. Specifically, we find that at Kp = 6 the lifetime can be <10 h at L > 5 for 100 s keV protons and at L > 4 for 100 s keV O + , which demonstrates that FLC scattering can be an important mechanism for the observed fast loss (τ < 10 h) of ring current ions during geomagnetic storms. Furthermore, we formulate an empirical formula for τ as a function of ion energy, pitch angle, position, species, and Kp. The empirical formula can be directly included in ring current models to account for the FLC scattering effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Observations and Simulations of Dropout Events and Flux Decays in October 2013: Comparing MEO Equatorial With LEO Polar Orbit

We compare ESA PROBA-V observations of electron flux at LEO with those from the NASA Van Allen Probes mostly at MEO for October 2013. Dropouts are visible at all energy during four storms from both satellites. Equatorially trapped electron fluxes are higher than at LEO by 10 2 (<1 MeV) to 10 5 (>2.5 MeV). We observe a quite isotropic structure of the outer belt during quiet times, contrary to the inner belt, and pitch angle dependence of high energy injection. We find a very good overlap of the outer belt at MEO and LEO at ~0.5 MeV. We use test-particle simulations of the energetic electrons trapped in the terrestrial magnetic field to study the outer radiation belt electron flux changes during geomagnetic storms. We show that the Dst (Disturbance storm time) effect during the main phase of a geomagnetic storm results in a betatron mechanism causing outward radial drift and a deceleration of the electrons. This outward drift motion is energy independent, pitch angle-dependent, and represents a significant distance (~1 L-shell at L = 5 for moderate storms). At fixed L-shell, this causes a decay of the LEO precipitating flux (adiabatic outward motion), followed by a return to the normal state (adiabatic inward motion) during main and recovery phases. Dst effect, associated with magnetopause shadowing and radial diffusion can explain the main characteristics of outer radiation belt electron dropouts in October 2013. We also use Fokker-Planck simulations with event-driven diffusion coefficients at high temporal resolution, to distinguish instantaneous loss from the gradual scattering that depopulates the slot region and the outer belt after storms. Simulations reproduce the slot formation and the gradual loss in the outer belt. The typical energy dependence of these losses leads to the absence of scattering for relativistic and ultra-relativistic electrons in the outer belt, oppositely to dropouts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Acceleration of Oxygen Ions In Dipolarization Events: 2. PSBL Distributions

Abstract This paper represents the second part of an investigation of the acceleration of energetic oxygen ions from encounters with a dipolarization front (DF), based on test particle tracing in the fields of an MHD simulation. In this paper, we focus on distributions in the plasma sheet boundary layer (PSBL). O + beams close to the plasma sheet boundary are found to be less pronounced and/or delayed against the H + beams. The reason is that these particles are accelerated by nonadiabatic motion in the duskward electric field such that O + ions gain the same amount of energy, but only 1/4 of the speed of protons. This causes a delay and larger equatorward displacement by the E × B drift. In contrast, the O + beams somewhat deeper inside the plasma sheet, where previously multiple proton beams were found, are accelerated at an earthward propagating DF just like H + , forming a field‐aligned beam at a similar speed as the lowest‐energy H + beam. We found that the source location depends on the adiabaticity of the orbit. For larger adiabaticity, the beam ions originate initially from the outer plasma sheet, but later from the opposite PSBL or lobe, but for low adiabaticity, sources are well inside the plasma sheet. The energy gained from a single encounter of a DF is comparable to the kinetic energy associated with the front speed. Assuming maximum speeds of 500–1,000 km/s, this yields a mass dependent acceleration of about 1–5 keV for protons and 20–80 keV for oxygen ions, independent of their charge state.

Birn, J.↗
DOE OSTI2021

Acceleration of Oxygen Ions in Dipolarization Events: 1. CPS Distributions

Abstract Using an MHD simulation of near tail reconnection associated with a flow burst and the collapse (dipolarization) of the inner tail in combination with test particle tracing we study the acceleration and flux increases of energetic oxygen ions (O + ). The characteristic orbits, distributions, and acceleration mechanisms are governed by the dimensionless parameter σ = ω ci t n , where ω ci is the ion gyro frequency and t n a characteristic Alfvén time of the MHD simulation. For σ < 1, central plasma sheet (CPS) populations after the passage of the dipolarization front are found to resemble half‐shells in velocity space oriented toward dusk. They originate from within the CPS and are energized typically by a single encounter of the region of enhanced cross‐tail electric field associated with the flow burst. For larger σ values ( σ > 1) the O + distributions resemble more closely those of protons, consisting of two counter‐streaming field‐aligned beams and an, albeit more tenuous and irregular, ring population perpendicular to the magnetic field. The existence of the beams, however, depends on suitable earthward moving source populations in the plasma sheet boundary layer or the adjacent lobes. The acceleration to higher energies is found to indicate a charge dependence, consistent with a dominance of more highly charged ions at energies of a few hundred keV. As in earlier simulations, the simulated fluxes show large anisotropies and nongyrotropic effects, phase bunching, and spatially and temporally localized beams.

Birn, J.↗
DOE OSTI2026

Revealing the Formation of the <20 MeV Inner Proton Radiation Belt at L ∼ 2 During the 10–11 May 2024 Superstorm

A new proton radiation belt was identified during the geomagnetic superstorm of 10–11 May 2024. To investigate its origin, we use an MHD-test particle simulation to model solar energetic proton (SEP) trapping and the evolution of the initial trapped proton population during the storm. The simulation weights the injected SEP population with interplanetary proton measurements and the initial trapped proton radial profile with energy-resolved Weather System Follow-on-Microwave (WSF-M) proton data, enabling quantitative comparison with observations. Results show that SEP contributions are minimal; the new ∼2–18 MeV belt at L ∼ 2 forms mainly through redistribution and energization of the initial trapped population within hours of the CME shock, driven by subsequent electric field impulses, rather than the initial shock alone. These results provide quantitative evidence based on realistic, plasma-dependent field dynamics and significantly advance understanding of inner radiation belt formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Computational investigation of hysteresis and phase equilibria of n-alkanes in a metal-organic framework with both micropores and mesopores

Abstract Adsorption hysteresis is a phenomenon related to phase transitions that can impact applications such as gas storage and separations in porous materials. Computational approaches can greatly facilitate the understanding of phase transitions and phase equilibria in porous materials. In this work, adsorption isotherms for methane, ethane, propane, and n-hexane were calculated from atomistic grand canonical Monte Carlo (GCMC) simulations in a metal-organic framework having both micropores and mesopores to better understand hysteresis and phase equilibria between connected pores of different size and the external bulk fluid. At low temperatures, the calculated isotherms exhibit sharp steps accompanied by hysteresis. As a complementary simulation method, canonical (NVT) ensemble simulations with Widom test particle insertions are demonstrated to provide additional information about these systems. The NVT+Widom simulations provide the full van der Waals loop associated with the sharp steps and hysteresis, including the locations of the spinodal points and points within the metastable and unstable regions that are inaccessible to GCMC simulations. The simulations provide molecular-level insight into pore filling and equilibria between high- and low-density states within individual pores. The effect of framework flexibility on adsorption hysteresis is also investigated for methane in IRMOF-1.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of Coulomb coupling on friction in strongly magnetized plasmas

The friction force on a test particle traveling through a plasma that is both strongly coupled and strongly magnetized is studied using molecular dynamics simulations. In addition to the usual stopping power component aligned antiparallel to the velocity, a transverse component that is perpendicular to both the velocity and Lorentz force is observed. This component, which was recently discovered in weakly coupled plasmas, is found to increase in both absolute and relative magnitude in the strongly coupled regime. Strong coupling is also observed to induce a third component of the friction force in the direction of the Lorentz force. These first-principles simulations reveal novel physics associated with collisions in strongly coupled, strongly magnetized, plasmas that are not predicted by existing kinetic theories. The effect is expected to influence macroscopic transport in a number of laboratory experiments and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗