SEARCH · Engineering Papers

Results for “ion trajectory 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Investigation of post-breakup Coulomb acceleration using a trajectory model

Intermediate mass fragments ejected during the deexcitation of excited projectilelike fragments may promptly decay following ejection; the daughter particles that are subsequently produced are subject to interactions with the residual nucleus that affect final-state observables, a process herein referred to as post-breakup Coulomb acceleration. A simple classical Coulomb interaction model was used to study modification of 8 Be (2 + ), 5 Li (3/2 – ), 7 Li (7/2 – ), 7 Be (7/2 – ), and states in 12 B emitted following heavy-ion collisions of 28 Si + 12 C at 35 MeV/nucleon. Here, in contrast to previous work studying 8 Be (2 + ), excellent agreement between simulation and experiment was obtained using only Coulomb forces when either a Lorentzian or R-matrix line shape was used to describe the initial relative energy rather than a Gaussian. In consideration of the obtained results, improvements to the model and evaluation of experimental data are discussed as future directions, but it was concluded that the effects observed in the present data can be accurately described using only elements of classical mechanics and that the process is largely understood for a wide range of state lifetimes and mass and charge (a)symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1994

Simulations of phase space distributions of storm time proton ring current

We use results of guiding-center simulations of ion transport to map phase space densities of the stormtime proton ring current. We model a storm as a sequence of substorm-associated enhancements in the convection electric field. Our pre-storm phase space distribution is an analytical solution to a steady-state transport model in which quiet-time radial diffusion balances charge exchange. This pre-storm phase space spectra at L approximately 2 to 4 reproduce many of the features found in observed quiet-time spectra. Using results from simulations of ion transport during model storms having main phases of 3, 6, and 12 hr, we map phase space distributions from the pre-storm distribution in accordance with Liouville's theorem. We find stormtime enhancements in the phase space densities at energies E approximately 30-160 keV for L approximately 2.5 to 4. These enhancements agree well with the observed stormtime ring current. For storms with shorter main phases (approximately 3 hr), the enhancements are caused mainly by the trapping of ions injected from open night side trajectories, and diffusive transport of higher-energy (greater than or approximately 160 keV) ions contributes little to the stormtime ring current. However, the stormtime ring current is augmented also by the diffusive transport of higher-energy ions (E greater than or approximately 160 keV) durinng stroms having longer main phases (greater than or approximately 6 hr). In order to account for the increase in Dst associated with the formation of the stormtime ring current, we estimate the enhancement in particle-energy content that results from stormtime ion transport in the equatorial magnetosphere. We find that transport alone cannot account for the entire increase in absolute value of Dst typical of a major storm. However, we can account for the entire increase in absolute value of Dst by realistically increasing the stormtime outer boundary value of the phase space density relative to the quiet-time value. We compute the magnetic field produced by the ring current itself and find that radial profiles of the magnetic field depression resemble those obtained from observational data.

Chen, Margaret W.↗

NASA NTRS1986

Three dimensional simulation of the operation of a hollow cathode electron emitter on the Shuttle orbiter

Several researchers have suggested using hollow cathodes as plasma contactors for electrodynamic tethers, particularly to prevent the shuttle orbiter from charging to large negative potentials. Previous studies have shown that fluid models with anomalous scattering can describe the electron transport in hollow cathode generated plasmas. An improved theory of the hollow cathode plasmas is developed and computational results using the theory are compared with laboratory experiment. Numerical predictions for a hollow cathode plasma source of the type considered for use on the shuttle are presented as are three-dimensional NASCAP/LEO calculations of the emitted ion trajectories and the resulting potentials in the vicinity of the orbiter. The computer calculations show that the hollow cathode plasma source makes vastly superior contact with the ionospheric plasma compared with either an electron gun or passive ion collection by the orbiter.

Davis, V. A.↗

DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗

NASA NTRS1992

The formation and structure of plasma wakes behind large high-voltage space platforms in ionosphere

Theory and particle simulation results are presented for ionospheric plasma flow over space platforms in the large-dimension and high-voltage range. Both the transient formation of the space-charge wake and its steady state structure are studied. The wake-side ion impact and current collection are obtained. It is found that the wake behind a high-voltage plate is characterized by two ion-rich sheaths embedded in a quasi-neutral background wake. The embedded sheath is formed by the ions passing through the sheath around the plate edge and serves as their trajectory path. Depending on the surface potential, the plate dimension, and the angle of attack, the embedded sheath may either extend downstream or curve back to the plate causing a high, localized ion flux density at the location it strikes.

Wang, J.↗

NASA NTRS2023

Electron Acceleration and Heating during Magnetic Reconnection in the Earth's Quasi-parallel Bow Shock

We perform a 2.5-dimensional particle-in-cell simulation of a quasi-parallel shock, using parameters for the Earth's bow shock, to examine electron acceleration and heating due to magnetic reconnection. The shock transition region evolves from the ion-coupled reconnection dominant stage to the electron-only reconnection dominant stage, as time elapses. The electron temperature enhances locally in each reconnection site, and ion-scale magnetic islands generated by ion-coupled reconnection show the most significant enhancement of the electron temperature. The electron energy spectrum shows a power law, with a power-law index around 6. We perform electron trajectory tracing to understand how they are energized. Some electrons interact with multiple electron-only reconnection sties, and Fermi acceleration occurs during multiple reflections. Electrons trapped in ion-scale magnetic islands can be accelerated in another mechanism. Islands move in the shock transition region, and electrons can obtain larger energy from the in-plane electric field than the electric potential in those islands. These newly found energization mechanisms in magnetic islands in the shock can accelerate electrons to energies larger than the achievable energies by the conventional energization due to the parallel electric field and shock drift acceleration. This study based on the selected particle analysis indicates that the maximum energy in the nonthermal electrons is achieved through acceleration in ion-scale islands, and electron-only reconnection accounts for no more than half of the maximum energy, as the lifetime of sub-ion-scale islands produced by electron-only reconnection is several times shorter than that of ion-scale islands.

Solar magnetic reconnection↗

DOE OSTI2025

String-based axial and helicity-flip GPDs: A comparison to lattice QCD

We construct an analytic, string-based representation of the nucleon’s axial and helicity-flip conformal moments of generalized parton distributions (GPDs) that holds for skewness and for both the quark and gluon channels. The starting point is the Mellin-Barnes resummation of the conformal partial-wave expansion, where the moments are parametrized by open- (Reggeon) and closed-string (Pomeron) trajectories with slopes determined by experimental form factors and meson/glueball spectroscopy. The forward limits are fixed by the empirical unpolarized and polarized parton distributions. Polynomiality, crossing symmetry, and support are satisfied by construction. After next-to-leading order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi/Efremov-Radyushkin-Brodsky-Lepage evolution to μ = 2 GeV our analytic framework (i) reproduces some of the currently available lattice moments of E and H˜ in the nonsinglet sector, (ii) predicts sea-quark and gluon polarized moments that will be testable by forthcoming simulations and experiments at Jefferson Lab and the future Electron-Ion Collider, and (iii) yields axial and helicity-flip GPDs in x-space in reasonable agreement with lattice QCD.

Gauge-gravity dualities↗

DOE OSTI2025

Light control of intramolecular nuclear dynamics by vortex electron localization

In strong-field ionization of molecules, intense light pulses are thought to have a negligible direct influence on atomic nuclei. Molecular dissociation is thus expected to be determined by the geometrical configuration of the molecular ion at the ionization instant. Contrary to this picture, we observe a counterintuitive electron-proton angular correlation and the formation of proton vortices following strong-field ionization of H 2 molecules by bicircularly polarized two-color laser fields. We explain this phenomenon by the pathway interference and localization of the residual H 2 ⁡+ electron in different angular-momentum states formed in the tail of the driving laser pulse. We validate this interpretation by combining a quantum-mechanical numerical simulation of the field-driven coupled electronic-nuclear dynamics and a semiclassical-trajectory model for the phase accumulation of the laser-driven electronic-nuclear wave packet. Our joint experimental and theoretical study reveals a general picture of vortex electron localization which can be used for controlling molecular-bond breaking with circularly polarized laser fields.

Atomic & molecular processes in external fields↗

NASA NTRS1995

Test-particle simulations in increasingly strong turbulence

Quasi-linear theory supposes that the energy in resonant fluctuations is small compared to that in the mean magnetic field. This is evident in the fact that the zeroth-order particle trajectories are helices about a mean field B(sub o) that is spatially uniform over many correlation lengths. However, in the solar wind it is often the case that the fluctuating part of the field is comparable in magnitude to the mean part. It is generally expected that quasi-linear theory remains viable for particles that are in resonance with a region of the fluctuation spectrum having only small energy density, but even so, care must be taken when comparing simulations to theoretical predictions. We have performed a series of test-particle simulations to explore the evolution of ion distributions in turbulent situations with varying levels of magnetic fluctuations. As delta-B/B(sub o) is increased the distinctions among absolute pitch angle (defined relative to B(sub o)), local pitch angle (defined relative to B(x)), and magnetic moment become important, some of them exhibiting periodic sloshing unrelated to the nonadiabatic processes of interest. Comparing and contrasting the various runs illustrates the phenomena that must be considered when the premise underlying quasi-linear theory are relaxed.

Pontius, D. H., Jr.↗

DOE OSTI2021

Strong-field ionization of water. II. Electronic and nuclear dynamics en route to double ionization

We investigate the role of nuclear motion and strong-field-induced electronic couplings during the double ionization of deuterated water using momentum-resolved coincidence spectroscopy. By examining the three-body dicationic dissociation channel, D$^{+}$/D$^{+}$/O, for both few- and multi-cycle laser pulses, strong evidence for intra-pulse dynamics is observed. The extracted angle- and energy-resolved double ionization yields are compared to classical trajectory simulations of the dissociation dynamics occurring from different electronic states of the dication. In contrast with measurements of single photon double ionization, pronounced departure from the expectations for vertical ionization is observed, even for pulses as short as 10~fs in duration. We outline numerous mechanisms by which the strong laser field can modify the nuclear wavefunction en-route to final states of the dication where molecular fragmentation occurs. Specifically, we consider the possibility of a coordinate-dependence to the strong-field ionization rate, intermediate nuclear motion in monocation states prior to double ionization, and near-resonant laser-induced dipole couplings in the ion. These results highlight the fact that, for small and light molecules such as D$_2$O, a vertical-transition treatment of the ionization dynamics is not sufficient to reproduce the features seen experimentally in the strong field coincidence double-ionization data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2022

Exploring the Inorganic Composition of the Asian Tropopause Aerosol Layer Using Medium-Duration Balloon Flights

Satellite observations have revealed an enhanced aerosol layer near the tropopause over Asia during the summer monsoon, called the “Asian Tropopause Aerosol Layer” (ATAL). In this work, aerosol particles in the ATAL were collected with a balloon-borne impactor near the tropopause region over India, using extended-duration balloon flights, in summer 2017 and winter 2018. The chemical composition of these particles was further investigated by quantitative analysis using offline ion chromatography. Nitrate (NO_3 ^1) and nitrite (NO_2 ^-) were found to be the dominant ions in the collected aerosols with values ranging between 87 and 343 ng m^(−3) at STP (standard temperature and pressure) during the summer campaign. In contrast, sulfate (SO_4 ^(2-)) levels were found to be above the detection limit (>10 ng m^(−3) at STP) only in winter. In addition, we determined the origin of the air masses sampled during the flights using the analysis of back trajectories as well as a convective proxy from cloud-top temperature fields derived from a geostationary satellite. The results obtained from this analysis were put into the context of large-scale transport and aerosol distribution using GEOS-Chem chemical transport model simulations. The first flight in summer 2017 which sampled an air mass within the Asian monsoon anticyclone (AMA), influenced by convection over Western China, was associated with particle size diameters from 0.05 to 0.15 µm. In contrast, the second flight sampled air masses at the edge of the AMA associated with a larger particle size radius (>2 µm) with a higher NO_2 ^- concentration. The sampled air masses in winter 2018 were likely affected by smoke from the Pacific Northwest fire event in Canada, which occurred 7 months before our campaign, associated with concentration enhancements of SO_4 (^2-) and Ca^(2+). Overall, our results suggest that nitrogen-containing particles represent a large fraction of cloud-free and in-cloud aerosols populating the ATAL, which is partially in agreement with the results from aircraft measurements during the StratoClim (Stratospheric and upper tropospheric processes for better climate predictions) campaign. The exact nature of those particles is still unknown, but their coincidences with subvisible cirrus clouds and their sizes suggest nitric acid trihydrate (NAT) as a possible candidate, as NAT has already been observed in the tropical upper troposphere and lower stratosphere in other studies. Furthermore, GEOS-Chem model simulations indicate that lightning NOx emissions could have significantly impacted the production of nitrate aerosols sampled during the summer of 2017.

Hazel Vernier↗

NASA NTRS2005

Design and Optimization of Low-thrust Orbit Transfers Using Q-law and Evolutionary Algorithms

Future space missions will depend more on low-thrust propulsion (such as ion engines) thanks to its high specific impulse. Yet, the design of low-thrust trajectories is complex and challenging. Third-body perturbations often dominate the thrust, and a significant change to the orbit requires a long duration of thrust. In order to guide the early design phases, we have developed an efficient and efficacious method to obtain approximate propellant and flight-time requirements (i.e., the Pareto front) for orbit transfers. A search for the Pareto-optimal trajectories is done in two levels: optimal thrust angles and locations are determined by Q-law, while the Q-law is optimized with two evolutionary algorithms: a genetic algorithm and a simulated-annealing-related algorithm. The examples considered are several types of orbit transfers around the Earth and the asteroid Vesta.

low-thrust orbit transfer↗

DOE OSTI2024

Decoding the Pair Distribution Function of Uranium in Molten Fluoride Salts from X-Ray Absorption Spectroscopy Data by Machine Learning

Thermal properties of actinides in molten salts are linked to the strongly disordered local environment of actinide ions. Here, we illustrate both the limitations of the commonly used fitting method for analysis of extended X-ray absorption fine structure (EXAFS) spectra in molten UF 4 and a possible solution using an "objective neural network - EXAFS" (ONNE) method. ONNE provides both extraction of the pair distribution function, as validated by its application to the EXAFS spectra calculated on molecular dynamics trajectory, and the EXAFS data reconstruction. The ONNE analysis of the molten UF4 has revealed reduction of the first nearest neighbor U-F coordination number, expansion of the U-F bond length and smaller contribution to the second shell compared to its crystalline counterpart. This method is therefore an attractive alternative to conventional EXAFS analysis and molecular dynamics simulations for studies of disordered environment of actinides in molten salts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2023

Kinetic modeling of solitary wave dynamics in a neutralizing ion beam

In this work, we characterize the formation and evolution of electrostatic solitary waves (ESWs) in the space-charge neutralization of ion beams using particle-in-cell simulations. These waves become excited when the electrons emitted from an external filament source initiate a two-stream instability in the beam. We show that such electrostatic waves become excited in both two-dimensional (2D) and three-dimensional (3D) beams with different shapes and sizes. Through a 1D Bernstein–Greene–Kruskal (BGK) analysis of the 2D beam, we find that the non-Maxwellian nature of the beam electrons gives rise to large-sized ESWs that are not predicted by BGK theory since it assumes a Maxwellian electron velocity distribution in the beam. Finally, we show that a 1D BGK theory is inadequate to describe ESWs in 3D beams because of complex electron trajectories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1992

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Guiding-center simulations of stormtime transport of ring-current and radiation-belt ions having first adiabatic invariants mu is approximately greater than 15 MeV/G (E is approximately greater than 165 keV at L is approximately 3) are surprisingly well described (typically within a factor of approximately less than 4) by the quasilinear theory of radial diffusion. This holds even for the case of an individual model storm characterized by substorm-associated impulses in the convection electric field, provided that the actual spectrum of the electric field is incorporated in the quasilinear theory. Correction of the quasilinear diffusion coefficient D(sub LL)(sup ql) for drift-resonance broadening (so as to define D(sub LL)(sup ql)) reduced the typical discrepancy with the diffusion coefficients D(sub LL)(sup sim) deduced from guiding-center simulations of representative-particle trajectories to a factor of approximately 3. The typical discrepancy was reduced to a factor of approximately 1.4 by averaging D(sub LL)(sup sim), D(sub LL)(sup ql), and D(sub LL)(sup rb) over an ensemble of model storms characterized by different (but statistically equivalent) sets of substorm-onset times.

Chen, Margaret W.↗

DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Monte Carlo Simulation of Electron Interactions in an MeV-STEM for Thick Frozen Biological Sample Imaging

A variety of volume electron microscopy techniques have been developed to visualize thick biological samples. However, the resolution is limited by the sliced section thickness (>30–60 nm). To preserve biological samples in a hydrated state, cryo-focused ion beam scanning electron microscopy has been developed, providing nm resolutions. However, this method is time-consuming, requiring 15–20 h to image a 10 μm thick sample with an 8 nm slice thickness. There is a pressing need for a method that allows the rapid and efficient study of thick biological samples while maintaining nanoscale resolution. The remarkable ability of mega-electron-volt (MeV) electrons to penetrate thick biological samples, even exceeding 10 μm in thickness, while maintaining nanoscale resolution, positions MeV-STEM as a suitable microscopy tool for such applications. Our research delves into understanding the interactions between MeV electrons and frozen biological specimens through Monte Carlo simulations. Single elastic scattering, plural elastic scattering, single inelastic scattering, and plural inelastic scattering events have been simulated. The electron trajectories, the beam profile, and the intensity change of electrons in each category have been investigated. Additionally, the effects of the detector collection angle and the focal position of the electron beam were investigated. As electrons penetrated deeper into the specimen, single and plural elastic scattered electrons diminished, and plural inelastic scattered electrons became dominant, and the beam profile became wider. Even after 10 μm of the specimen, 42% of the MeV electrons were collected within 10 mrad. This confirms that MeV-STEM can be employed to study thick biological samples.

electron specimen interaction↗