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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Production of Magnetic Turbulence by Cosmic Rays Drifting Upstream of Supernova Remnant Shocks

I will present results of our recent two- and three-dimensional Particle-In-Cell simulations of magnetic-turbulence production by cosmic-ray ions drifting upstream of supernova remnant shocks. These studies' aim is twofold: test recent predictions of strong amplification in short wavelength, non-resonant wave modes, and study the subsequent evolution of the magnetic turbulence, including its backreaction on cosmic-ray trajectories. We confirm that the drifting cosmic rays give rise to a turbulent magnetic field, but show that an oblique filamentary mode grows more rapidly than the non-resonant parallel modes found in analytical theory. The field perturbations grow more slowly than estimated using a quasi-linear analytical approach for the parallel plane-wave mode, and saturate in amplitude at deltaB/B approximately equal to 1. The backreaction of the magnetic turbulence on the particles leads to an alignment of the bulk-flow velocities of the cosmic rays and the background medium. This is an essential characteristic of cosmic ray-modified shocks: the upstream flow speed is continuously changed by the cosmic rays. The reduction of relative drift between cosmic rays and background medium accounts for the saturation of the instability at only moderate magnetic-field amplitudes. It is possible that the prolonged magnetic field growth observed in recent MHD simulations results from a cosmic-ray current assumed to be constant and thus immune to the backreaction from the turbulent field. We speculate that the parallel plane-wave mode found in analytical treatments very quickly leads co filamentation, which we observe in our PIC modeling and is also apparent in the MHD simulations.

Stroman, Thomas↗

NASA NTRS2004

Numerical Simulation of Two-grid Ion Optics Using a 3D Code

A three-dimensional ion optics code has been developed under NASA's Project Prometheus to model two grid ion optics systems. The code computes the flow of positive ions from the discharge chamber through the ion optics and into the beam downstream of the thruster. The rate at which beam ions interact with background neutral gas to form charge exchange ions is also computed. Charge exchange ion trajectories are computed to determine where they strike the ion optics grid surfaces and to determine the extent of sputter erosion they cause. The code has been used to compute predictions of the erosion pattern and wear rate on the NSTAR ion optics system; the code predicts the shape of the eroded pattern but overestimates the initial wear rate by about 50%. An example of use of the code to estimate the NEXIS thruster accelerator grid life is also presented.

ion optics↗

DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗

NASA NTRSDate not supplied

Phase-Free Orbital Element Model Designed to Enable Rapid Assessment of Eclipse and Radiation Profiles for Low-Thrust Spiral Transfers Around the Earth

The Gateway Power and Propulsion Element (PPE) will be the first low-thrust solar electric ion propulsion mission to transfer from a highly elliptical Earth-bound orbit to a southern near-rectilinear halo orbit (NRHO) at the Earth-Moon L2 point. Due to the low thrust nature of the transfer orbit, it is desirable to locate viable trajectories that minimize the time spent in the Van Allen radiation belts to reduce solar array degradation and keep radiation dosages below design limits. In addition to radiation, low thrust trajectories which spiral around the Earth will pass through at least one or more seasons of Earth eclipses. The solar electric propulsion system relies on sunlight to generate the power necessary to operate the ion thrusters. Therefore, it is advantageous to find trajectories which also minimize the amount of time spent in eclipse and avoid excessive battery draw-down periods. We present an analytical method which rapidly approximates low-thrust spiral trajectories, fit to high-fidelity simulated data and parameterized to allow for changes in vehicle thrust characteristics, that is post-processed to determine eclipse profiles and time spent in the belts using a novel approach that estimates the geometry of the belts using a first-order dipole approximation of the Earth’s magnetic field. This method can be used to find satisfactory launch dates and orbit orientations that can serve as initial guesses when optimizing such missions in high-fidelity software.

low thrust trajectory design↗

NASA NTRSDate not supplied

Phase-Free Orbital Element Model Designed to Enable Rapid Assessment of Eclipse and Radiation Profiles for Low-Thrust Spiral Transfers Around the Earth

low thrust trajectory design↗

DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

NASA NTRS2006

21st Century Lunar Exploration: Advanced Radiation Exposure Assessment

On January 14, 2004 President George W Bush outlined a new vision for NASA that has humans venturing back to the moon by 2020. With this ambitious goal, new tools and models have been developed to help define and predict the amount of space radiation astronauts will be exposed to during transit and habitation on the moon. A representative scenario is used that includes a trajectory from LEO to a Lunar Base, and simplified CAD models for the transit and habitat structures. For this study galactic cosmic rays, solar proton events, and trapped electron and proton environments are simulated using new dynamic environment models to generate energetic electron, and light and heavy ion fluences. Detailed calculations are presented to assess the human exposure for transit segments and surface stays.

Anderson, Brooke↗

DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2000

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

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

Leung, Wing C.↗

DOE OSTI2025

Jet Drift in Heavy Ion Collisions

We introduce a sub-eikonal anisotropic contribution to jetbroadening, “jet drift”, that couples to the flow of the medium, showing that this effect results in a deflection of hard partons, and thus jets, in the direction of the medium flow. Next, we study this effect in a full-fledged hybrid transport simulation of √s = 5.02 TeV PbPb collisions at the LHC, tracking trajectories of hard partons with perturbative energy loss and drift. We show that sub-eikonal anisotropic effects, including flow-mediated jet drift, are sensitive to properties of the medium that traditional eikonal isotropic effects are insensitive to, demonstrating that including these effects leads to modifications to jet and hard particle observables that survive averaging over events. We show that jet drift leads to an enhancement of the elliptic flow (v2) of hard particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Effects of convection electric field on upwelling and escape of ionospheric O(+)

A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.

Cladis, J. B.↗

DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

NASA NTRS2012

An Efficient Deterministic Approach to Model-based Prediction Uncertainty Estimation

Prognostics deals with the prediction of the end of life (EOL) of a system. EOL is a random variable, due to the presence of process noise and uncertainty in the future inputs to the system. Prognostics algorithm must account for this inherent uncertainty. In addition, these algorithms never know exactly the state of the system at the desired time of prediction, or the exact model describing the future evolution of the system, accumulating additional uncertainty into the predicted EOL. Prediction algorithms that do not account for these sources of uncertainty are misrepresenting the EOL and can lead to poor decisions based on their results. In this paper, we explore the impact of uncertainty in the prediction problem. We develop a general model-based prediction algorithm that incorporates these sources of uncertainty, and propose a novel approach to efficiently handle uncertainty in the future input trajectories of a system by using the unscented transformation. Using this approach, we are not only able to reduce the computational load but also estimate the bounds of uncertainty in a deterministic manner, which can be useful to consider during decision-making. Using a lithium-ion battery as a case study, we perform several simulation-based experiments to explore these issues, and validate the overall approach using experimental data from a battery testbed.

Daigle, Matthew J.↗

DOE OSTI2025

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Waveform-dependent air fluorescence from neutral and ionic nitrogen molecules

Laser-induced air fluorescence in the ultraviolet regime is primarily attributed to transitions between the C and B states in excited neutral N 2 molecules and between the B and X states in N$^+_2$ ions. However, the mechanism underlying the former remains contentious, as direct population to the C state by light fields is forbidden by electron spin constraints. In this work, we investigate the mechanism of air fluorescence from excited neutral N 2 molecules by carrier-envelope phase–stabilized sub–4 femtosecond pulses. Our results show that fluorescence from N$^+_2$ ions reaches a maximum with cosine-like pulses, while fluorescence from excited neutral N 2 molecules peaks with sine-like pulses. In addition, by scanning the chirp of the driving pulse, we find that ionic fluorescence is maximized with chirp-free pulses, whereas neutral fluorescence favors negatively chirped pulses. These observations, supported by classical trajectory Monte Carlo simulations, support the mechanism of intersystem crossing from excited spin-singlet states, which are populated via recollision-induced strong-field excitation.

Science & Technology - Other Topics↗

NASA NTRS1998

Wall effects observed in tissue-equivalent proportional counters from 1.05 GeV/nucleon iron-56 particles

Tissue-equivalent proportional counters (TEPCs) have been used to measure energy deposition in simulated volumes of tissue ranging in diameter from 0.1 to 10 microm. There has been some concern that the wall used to define the volume of interest could influence energy deposition within the sensitive volume because it has a density significantly greater than that of the cavity gas. These effects become important for high-velocity heavy ions. Measurements of energy deposition were made for 1 GeV/nucleon iron particles in a TEPC simulating a 1-microm-diameter sphere of tissue. The TEPC was nested within a particle spectrometer that provided identification and flight path of individual particles. Energy deposition was studied as a function of pathlength through the TEPC. Approximately 30% of the energy transfer along trajectories through the center of the detector escapes the sensitive volume. The response of the TEPC, for trajectories through the detector, is always larger than calculations for energy loss in a homogeneous medium. This enhancement is greatest for trajectories near the cavity/wall interface. An integration of the response indicates that charged-particle equilibrium is essentially achieved for a wall thickness of 2.54 mm. However, estimates of the linear energy transfer for the incident particles are influenced by these wall effects.

Non-NASA Center↗