Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle tracing”

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 NTRS1990

Shock-drift particle acceleration in superluminal shocks - A model for hot spots in extragalactic radio sources

Shock-drift acceleration at relativistic shock fronts is investigated using a fully relativistic treatment of both the microphysics of the shock-drift acceleration and the macrophysics of the shock front. By explicitly tracing particle trajectories across shocks, it is shown how the adiabatic invariance of a particle's magnetic moment breaks down as the upstream shock speed becomes relativistic, and is recovered at subrelativistic velocities. These calculations enable the mean increase in energy of a particle which encounters the shock with a given pitch angle to be calculated. The results are used to construct the downstream electron distribution function in terms of the incident distribution function and the bulk properties of the shock. The synchrotron emissivity of the transmitted distribution is calculated, and it is demonstrated that amplification factors are easily obtained which are more than adequate to explain the observed constrasts in surface brightness between jets and hot spots.

Begelman, Mitchell C.↗
NASA NTRS1994

Energy density of ionospheric and solar wind origin ions in the near-Earth magnetotail during substorms

Comprehensive energy density studies provide an important measure of the participation of various sources in energization processes and have been relatively rare in the literature. We present a statistical study of the energy density of the near-Earth magnetotail major ions (H(+), O(+), He(++), He(+)) during substorm expansion phase and discuss its implications for the solar wind/magnetosphere/ionosphere coupling. Our aim is to examine the relation between auroral activity and the particle energization during substorms through the correlation between the AE indices and the energy density of the major magnetospheric ions. The data we used here were collected by the charge-energy-mass (CHEM) spectrometer on board the Active Magnetospheric Particle Trace Explorer (AMPTE)/Charge Composition Explorer (CCE) satellite in the near-equatorial nightside magnetosphere, at geocentric distances approximately 7 to 9 R(sub E). CHEM provided the opportunity to conduct the first statistical study of energy density in the near-Earth magnetotail with multispecies particle data extending into the higher energy range (greater than or equal to 20 keV/E). the use of 1-min AE indices in this study should be emphasized, as the use (in previous statistical studies) of the (3-hour) Kp index or of long-time averages of AE indices essentially smoothed out all the information on substorms. Most distinct feature of our study is the excellent correlation of O(+) energy density with the AE index, in contrast with the remarkably poor He(++) energy density - AE index correlation. Furthermore, we examined the relation of the ion energy density to the electrojet activity during substorm growth phase. The O(+) energy density is strongly correlated with the pre-onset AU index, that is the eastward electrojet intensity, which represents the growth phase current system. Our investigation shows that the near-Earth magnetotail is increasingly fed with energetic ionospheric ions during periods of enhanced dissipation of auroral currents. The participation of the ionosphere in the substorm energization processes seems to be closely, although not solely, associated with the solar wind/magnetosphere coupling. That is, the ionosphere influences actively the substorm energization processes by responding to the increased solar wind/magnetosphere coupling as well as to the unloading dissipation of stored energy, with the increased feeding of new material into the magnetosphere.

Daglis, Loannis A.↗
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 OSTI2025

Distortions in charged-particle images of laser direct-drive inertial confinement fusion implosions

Energetic charged particles generated by inertial confinement fusion (ICF) implosions encode information about the spatial morphology of the hotspot and dense fuel during the time of peak fusion reactions. The knock-on deuteron imager (KoDI) was developed at the Omega Laser Facility to image these particles in order to diagnose low-mode asymmetries in the hotspot and dense fuel layer of cryogenic deuterium–tritium ICF implosions. However, the images collected are distorted in several ways that prevent reconstruction of the deuteron source. In this paper, we describe these distortions and a series of attempts to mitigate or compensate for them. We present several potential mechanisms for the distortions, including a new model for scattering of charged particles in filamentary electric or magnetic fields surrounding the implosion. Particle-tracing is used to create synthetic KoDI data based on the filamentary field model that reproduces the main experimentally observed image distortions. We conclude that the filamentary scattering model best matches the observed image distortions. Finally, we discuss potential impacts of filamentary fields on other charged-particle diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Geomagnetic transmission disturbances and heavy-ion fluences observed in low Earth orbit during the solar energetic particle events of October 1989

The large solar energetic particle (SEP) events and simultaneous large geomagnetic disturbances observed during October 1989 posed a significant, rapidly evolving space radiation hazard. Using data from the GOES-7, NOAA-10, IMP-8 and LDEF satellites, we determined the geomagnetic transmission, heavy ion fluences, mean Fe ionic charge state, and effective radiation hazard observed in low Earth orbit (LEO) for these SEPs. We modeled the geomagneitc transmission by tracing particles through the combination to the internal International Geomagnetic Reference Field (IGRF) and the Tsyganenko (1989) magnetospheric field models, extending the modeling to large geomagnetic disturbances. We used our results to assess the radiation hazard such very large SEP events would pose in the anticipated 52 deg inclination space station orbit.

Boberg, P. R.↗
DOE OSTI2023

Loss of energetic ions due to n = 1 internal kink instability in HL-2M

Effects of three-dimensional (3D) perturbations due to an unstable n = 1 (n is the toroidal mode number) internal kink (IK) on the energetic particles (EPs) are systematically investigated for the HL-2M tokamak, utilizing the MARS-F/K code and a recently developed test particle tracing module. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code, is chosen for this study. In general, the 3D perturbation associated with an unstable IK is found to affect the EP drift orbit, confinement, and loss in HL-2M. The instability-induced EP loss fraction is found to be typically less than 10%, without counting for the prompt orbit loss associated with the 2D equilibrium field for counter-current particles. The latter reaches about 16% in HL-2M. For co-current EPs, a 100 G 3D magnetic field (inside the plasma) due to the IK does not induce any EP loss assuming a static perturbation. A sawtooth-like time-varying perturbation field, with the peak amplitude reaching 1000 G, can however produce about 30% loss for the co-current EPs in HL-2M. The majority of lost EPs tend to strike the lower divertor region, with a small fraction of particles striking the low-field side mid-plane region of the limiting surface.

36 MATERIALS SCIENCE↗
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 OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Influence of Photoemission Geometry on Timing and Efficiency in 4D Ultrafast Electron Microscopy

Broader adoption of 4D ultrafast electron microscopy (UEM) for the study of chemical, materials, and quantum systems is being driven by development of new instruments as well as continuous improvement and characterization of existing technologies. Perhaps owing to the still-high barrier to entry, the full range of capabilities of laser-driven 4D UEM instruments has yet to be established, particularly when operated at extremely low beam currents (~fA). Accordingly, with an eye on beam stability, we have conducted particle tracing simulations of unconventional off-axis photoemission geometries in a UEM equipped with a thermionic-emission gun. Specifically, we have explored the impact of experimentally adjustable parameters on the time-of-flight (TOF), the collection efficiency (CE), and the temporal width of ultrashort photoelectron packets. The adjustable parameters include the Wehnelt aperture diameter (DW), the cathode set-back position (Ztip), and the position of the femtosecond laser on the Wehnelt aperture surface relative to the optic axis (Rphoto). Notable findings include significant sensitivity of TOF to DW and Ztip, as well as non-intuitive responses of CE and temporal width to varying Rphoto. As a means to improve accessibility, practical implications and recommendations are emphasized wherever possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Delineation of the impact on temporal behaviors of off-axis photoemission in an ultrafast electron microscope

Efforts to push the spatiotemporal imaging-resolution limits of femtosecond laser-driven ultrafast electron microscopes (UEMs) to the combined angstrom–fs range will benefit from stable sources capable of generating high bunch charges. Recent demonstrations of unconventional off-axis photoemitting geometries are promising, but connections to the observed onset of structural dynamics are yet to be established. Here we use the in-situ photoexcitation of coherent phonons to quantify the relative time-of-flight (r-TOF) of photoelectron packets generated from the Ni Wehnelt aperture and from a Ta cathode set-back from the aperture plane. We further support the UEM experiments with particle-tracing simulations of the precise electron-gun architecture and photoemitting geometries. In this way, we measure discernible shifts in electron-packet TOF of tens of picoseconds for the two photoemitting surfaces. Furthermore, these shifts arise from the impact that the Wehnelt-aperture off-axis orientation has on the electron-momentum distribution, which modifies both the collection efficiency and the temporal-packet distribution relative to on-axis emission. Future needs are identified; we expect this and other developments in UEM electron-gun configuration to expand the range of material phenomena that can be directly imaged on scales commensurate with fundamental structural dynamics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2026

Persistence of deuterium and tritium nuclear spin-polarization in presence of high-frequency plasma waves

We present first-principles numerical calculations of the depolarization rate of spin-polarized deuterium and tritium nuclei in realistic tokamak plasmas, driven by resonant interactions with plasma waves. Backed up by first-of-a-kind linear and nonlinear simulations, we find that alpha particle-driven Alfvénic modes cause only negligible depolarization, which is contrary to expectations in prior literature. Other Alfvénic instabilities can in principle degrade polarization, but only under conditions unlikely to be realized on transport timescales. By combining full-orbit particle tracing with a dedicated depolarization solver, we demonstrate that wave-driven depolarization is surprisingly weak in SPARC and ITER-scale devices. These results provide strong evidence that spin-polarized fuel can maintain its polarization long enough to boost fusion reactivity, opening a viable path toward substantially enhanced performance in magnetic confinement fusion power plants.

Cook, James William Setchfield [UK Atomic Energy A↗
DOE OSTI2025

Accelerating discoveries at DIII-D with the Integrated Research Infrastructure

DIII-D research is being accelerated by leveraging high performance computing (HPC) and data resources available through the National Energy Research Scientific Computing Center (NERSC) Superfacility initiative. As part of this initiative, a high-resolution, fully automated, whole discharge kinetic equilibrium reconstruction workflow was developed that runs at the NERSC for most DIII-D shots in under 20 min. This has eliminated a long-standing research barrier and opened the door to more sophisticated analyses, including plasma transport and stability. These capabilities would benefit from being automated and executed within the larger Department of Energy Advanced Scientific Computing Research program’s Integrated Research Infrastructure (IRI) framework. The goal of IRI is to empower researchers to meld DOE’s world-class research tools, infrastructure, and user facilities seamlessly and securely in novel ways to radically accelerate discovery and innovation. For transport, we are looking at producing flux matched profiles and also using particle tracing to predict fast ion heat deposition from neutral beam injection before a shot takes place. Our starting point for evaluating plasma stability focuses on the pedestal limits that must be navigated to achieve better confinement. This information is meant to help operators run more effective experiments, so it needs to be available rapidly inside the DIII-D control room. So far this has been achieved by ensuring the data is available with existing tools, but as more novel results are produced new visualization tools must be developed. In addition, all of the high-quality data we have generated has been collected into databases that can unlock even deeper insights. This has already been leveraged for model and code validation studies as well as for developing AI/ML surrogates. The workflows developed for this project are intended to serve as prototypes that can be replicated on other experiments and can be run to provide timely and essential information for ITER, as well as next stage fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Increasing the isotropic properties of a discharge spark chamber with a large gap

The recording of particle tracks in the form of sparks or streamers running at different angles toward the electric field of the chamber was studied. The investigations were performed with two discharge spark chambers consisting of boxes filled with neon at 1 atm. The spark chambers were powered by a pulse with an amplitude of 300 kV applied from a pulse generator. The recording of particle tracks in the form of sparks took place at track angles smaller than 30 deg with respect to the electrical field. At angle larger than 40 deg the tracks are recorded in the form of streamers. Particle traces with respect to the angle are given.

Nurgozhin, N. N.↗
NASA NTRS1985

Computation of three-dimensional shock wave and boundary-layer interactions

Computations of the impingement of an oblique shock wave on a cylinder and a supersonic flow past a blunt fin mounted on a plate are used to study three dimensional shock wave and boundary layer interaction. In the impingement case, the problem of imposing a planar impinging shock as an outer boundary condition is discussed and the details of particle traces in windward and leeward symmetry planes and near the body surface are presented. In the blunt fin case, differences between two dimensional and three dimensional separation are discussed, and the existence of an unique high speed, low pressure region under the separated spiral vortex core is demonstrated. The accessibility of three dimensional separation is discussed.

Hung, C. M.↗
NASA NTRS1985

Graphics and flow visualization in computational fluid dynamics

Techniques for displaying two- and three-dimensional flowfield solutions are described. Several methods of illustrating flow structure are addressed including particle tracing, simulated oil flow, and shock finding. These are incorporated into an interactive graphics program for CFD flowfields, called PLOT3D. Emphasis is made on the difficulty in visualizing three-dimensional flow features, and the importance of color, fast 3D image manipulation, and dynamic movie play-back in displaying such flows. The need for advanced algorithms to identify shock waves, vortices, and separation lines is pointed out. It is likely that the supercomputer will be needed for this process because of the size of 3D and/or unsteady CFD databases.

Buning, P. G.↗