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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

IRAD

IRAD (Ion Irradiation and Radiation Damage) is an open-source GUI for SRIM-like vacancy production and implanted ion profile calculations. Using updated databases, IRAD is accessible at https://code.ornl.gov/liny/irad for Windows systems. It performs energy-corrected Iradina calculations, automatically determining dpa and implanted ion concentration, considering total ion fluence. Additional results like ion trajectories, final positions, stopping power distribution, etc., are available. Advanced settings permit customization of parameters such as stopping power database, incident ion angle, and 3D target simulations.

Lin, Yan-Ru [Oak Ridge National Laboratory (ORNL),↗

NASA NTRS2005

Ground testing of the Li-ion batteries in support of JPL's 2003 Mars Exploration Rover Mission

In early 2004, JPL successfully landed two Rovers, named Spirit and Opportunity, on the surface of Mars after traveling > 300 million miles over a 6-7 month period. In order to operate for extended duration on the surface of Mars, both Rovers are equipped with rechargeable Lithium-ion batteries, which were designed to aid in the launch, correct anomalies during cruise, and support surface operations in conjunction with a triple-junction deployable solar arrays. The requirements of the Lithium-ion battery include the ability to provide power at least 90 sols on the surface of Mars, operate over a wide temperature range (-20(deg)C to +4O(deg)C), withstand long storage periods (e.g., cruise period), operate in an inverted position, and support high currents (e.g., firing pyro events). In order to determine the viability of meeting these requirements, ground testing was performed on a Rover Battery Assembly Unit (RBAU), consisting of two 8-cell 8 Ah lithium-ion batteries connected in parallel. The RBAU upon which the performance testing was performed is nearly identical to the batteries incorporated into the two Rovers currently on Mars. The testing performed includes, (a) performing initial characterization tests (discharge capacity at different temperatures), (b) simulating the launch conditions, (c) simulating the cruise phase conditions (including trajectory corrections), (d) simulating the entry, decent, and landing pulse load profile (if required to support the pyros) (e) simulating the Mars surface operation mission simulation conditions, as well as, (f) assessing performance capacity loss and impedance characteristics as a function of temperature and life. As will be discussed, the lithium-ion batteries (fabricated by LithiodYardney, Inc.) were demonstrated to far exceed the requirements defined by the mission, and are projected to support an extended mission (> 2 years) with margin to spare.

Li-ion batteries↗

DOE OSTI2025

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗

DOE OSTI2024

Heavy ion beam probe for Wendelstein 7-X measurement capabilities as projected through its design

A heavy ion beam probe (HIBP) diagnostic is being developed for studies of plasma equilibrium and turbulence in the optimized Wendelstein 7-X (W7-X) stellarator. Operation of W7-X has experimentally demonstrated that its optimized magnetic field results in improved neoclassical particle confinement and, as a result, turbulence is the predominant cause of energy transport. The HIBP will have the unique ability to provide experimental data needed to complement models of both neoclassical and turbulent transport. It will acquire direct measurements in the W7-X plasma interior of the electric potential (needed for understanding ambipolar particle flux) and fluctuations of electron density and potential (needed for understanding turbulence). The HIBP for W7-X will inject singly charged ion beams with energies of up to 2 MeV and is designed to access the upper cross section of the W7-X plasma. We use trajectory simulations to illustrate the plasma coverage that the diagnostic can achieve in the reference magnetic configurations of W7-X. We calculate beam signal levels, discuss anticipated measurement sensitivity of broadband fluctuations of electron density and plasma potential, and show how they depend on plasma density. We also discuss the diagnostic sensitivity to equilibrium plasma potential.

Instruments & Instrumentation↗

NASA NTRS1994

Particle/fluid simulations of an eruptive flare: Identifying the field-aligned currents responsible for the hard x-rays

While magnetohydrodynamics (MHD) can provide a reasonable description of the overall magnetic reconnection that is believed to drive flares, additional, and often separate processes have to be envoked to in order to explain the electron acceleration that is responsible for many of the observed flare emissions. A new model that incorporates the dynamic coronal current sheets, the reconnection site, and possible electron acceleration processes is developed through the use of two-dimensional particle and modified two-fluid simulations. The specific example of an eruptive flare driven by the coalescence of flux tubes supported by prescribed photospheric current elements is evaluated. It is shown that the electrons and ions have differential trajectories through the coronal current sheet which leads to the development of additonal plasma currents that flow around the surface of the current sheet. These surface currents are explicitly neglected in MHD but they are vital to the flare dynamics because they divert current from the coronal current sheet into the chromosphere, producing an effective resistivity that aids the development of fast reconnection. Because the surface currents are in the plane of the magnetic field, electrons in them experience strong acceleration and can account for the observed hard X-ray emissions. Model predictions are compared with observed time profiles of hard X-ray emissions and Doppler shifts seen in soft X-ray line emissions and are able to account for such features as (1) the asymmetry in the rise and decay time of the hard X-rays, (2) the apparent delay between the largest Doppler shifts and the hard X-ray peak, and (3) the relatively low intensity of the blue-shifted component. The use of particle and fluid simulations is important because it provides different, but complementary treatments of the electron acceleration, the global magnetic morphology, and the flare current system.

Winglee, R. M.↗

DOE OSTI2025

Structure–Activity Relationships in Ether-Functionalized Solid-State Metal–Organic Framework Electrolytes

The structure–property relationships of metal–organic framework (MOF)-based solid-state electrolytes are not well understood. Herein, a systematic investigation of 12 Zr(IV)-based UiO-66 MOFs with varying ether-chain functional groups was carried out to elucidate the critical microscopic interactions that facilitate improved solid-state electrolyte performance. Enhanced sampling molecular dynamics (MD) simulations were employed and revealed a three-tier ion hopping mechanism: linker–linker hopping, linker-counterion hopping, and counterion-counterion hopping. Detailed structural analysis of the MD trajectories revealed that the chemistry and morphology of the linker groups affect the relative stability and population distribution of the electrolyte components, such that crown-ether-based linker groups enhance the probability of extended, low-barrier ion percolation pathways. As a result, we were able to tune the ionic conductivities by rationally manipulating the counterion distributions, linker binding strengths, and the configurational entropy (multivariability of the linkers). The resulting performance of these MOF-based solid-state electrolytes was significantly enhanced, with a methoxy-functionalized framework (UiO-66-L1 100 ) achieving high ionic conductivities of 2.32 × 10 –4 S/cm and 2.07 × 10 –3 S/cm at 30 °C and 90 °C, respectively, an order of magnitude greater than other all-solid-state MOF electrolyte systems. The electrolyte stability was evaluated with LiIn|LPSCl|MOF:LiTFSI|LPSCl|LiIn symmetric cells, showing excellent Li plating/stripping processes for over 2 months.

36 MATERIALS SCIENCE↗

NASA NTRS2006

An Update of the Ground Testing of the Li-ion Batteries in Support of JPL's 2003 Mars Exploration Rover Mission

In early 2004, JPL successfully landed two Rovers, named Spirit and Opportunity, on the surface of Mars after traveling > 300 million miles over a 6-7 month period. In order to operate for extended duration (>9 months), both Rovers are equipped with rechargeable Lithium-ion batteries, which have enabled operation for over 854 and 834 Sols of operation, respectively, to date. Given that the batteries were required to support the mission for 90 Sols of operation by design, it is significant that the batteries have demonstrated over a nine fold increase in life over mission objectives. In addition to supporting the surface operations in conjunction with a triple-junction deployable solar arrays, the batteries were designed to aid in the launch and the EDL pyros, and allow for anomalies during cruise. In summary, the requirements of the Lithium-ion battery include the ability to provide power at least 90 sols on the surface of Mars, operate over a wide temperature range (-20 C to +30 C), withstand long storage periods (e.g., cruise period), operate in an inverted orientation, and support high current pulses (e.g., firing pyro events). In order to determine the viability of meeting these requirements, ground testing was performed on a Rover Battery Assembly Unit (RBAU), consisting of two 8-cell 10 Ah lithium-ion batteries connected in parallel. The RBAU upon which the performance testing was performed is nearly identical to the batteries incorporated into the two Rovers currently on Mars. The testing includes, (a) performing initial characterization tests (discharge capacity at different temperatures), (b) simulating the launch conditions, (c) simulating the cruise phase conditions (including trajectory correction maneuvers), (d) simulating the entry, decent, and landing (EDL) pulse load profile (required to support the pyros) (e) simulating the Mars surface operation mission simulation conditions, as well as, (f) assessing capacity loss and impedance characteristics as a function of temperature and life. This paper provides further detail to previously reported results1 of the RBAU testing program, especially with regard to the life characteristics. To date, the lithium-ion batteries (fabricated by Lithion/Yardney, Inc.) have been demonstrated to far exceed the requirements defined by the mission, both on Mars and on the ground, and are projected to support an extended mission (> 4 years).

Mars Exploration Rovers (MER) Mission↗

DOE OSTI2023

Dynamic Monkey Bar Mechanism of Superionic Li‐ion Transport in LiTaCl 6

Abstract The LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm −1 ). However, the mechanism has been unclear. We train machine‐learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on‐the‐fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li‐ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large‐amplitude vibrations of Cl − ions in TaCl 6 − enable the fast Li‐ion transport by allowing dynamic breaking and reforming of Li−Cl bonds across the space in between the TaCl 6 − octahedra. We term this process the dynamic‐monkey‐bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all‐solid‐state lithium batteries.

Lei, Ming↗

DOE OSTI2022

All-optical nanoscopic spatial control of molecular reaction yields on nanoparticles

Molecular adsorbate reactions on nanoparticles play a fundamental role in areas such as nano-photocatalysis, atmospheric, and astrochemistry. They can be induced, enhanced, and controlled by field localization and enhancement on the nanoparticle surface. In particular, the ability to perform highly controlled near-field-mediated reactions is key to deepening our understanding of surface photoactivity on nanosystems. Here, using reaction nanoscopy, we experimentally demonstrate all-optical nanoscopic control of surface reaction yields by tailoring the near fields on nanoparticles with waveform-controlled linear and bicircular two-color laser pulses, respectively. We observe site-selective proton emission from the dissociative ionization of adsorbate molecules on Si O 2 nanoparticles as a function of the polarization and relative phase of the two-color pulses. The angularly resolved close-to-uniform mapping between the surface reaction yields and the measured ion momentum enables the observation and spatial control of molecular reactions on the nanoparticle surface with nanoscopic resolution. The experimental results are modeled and reproduced qualitatively by classical trajectory Monte Carlo simulations. Our work paves the way toward reliable all-optical control of photocatalytic chemical reactions on nanoscale surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Self‐consistent transport simulation of boron dust particle injection in the peripheral plasma in Large Helical Device

Abstract The trajectories and the ablation positions of boron dust particles dropped from an impurity powder dropper in the peripheral plasma in the Large Helical Device (LHD) were calculated using a three‐dimensional edge plasma simulation code (EMC3‐EIRENE) and a dust transport simulation code (DUSTT). The simulation shows that the trajectory of the boron dust particles is deflected at the upper divertor leg due to the effect of the hydrogen plasma flow, and the ablation positions of the dust particles in an ergodic layer change toward the outboard side of the torus for higher plasma densities. The effect of the boron ion flow in the divertor leg on the deflection is investigated by coupling the two codes self‐consistently. The simulation predicts that the boron ions in the divertor leg, which are produced by sputtering on the divertor plates, which do not affect the change in the ablation positions. It also shows that the ablation positions move toward the inboard side and approach the Last Closed Flux Surface (LCFS) in case of increased boron dust drop rates, which is caused by the lowered plasma flow in the upper divertor leg due to the lowered electron temperature by radiation cooling by the dropped dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Dynamic Monkey Bar Mechanism of Superionic Li–ion Transport in LiTaCl 6

Here, the LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm –1 ). However, the mechanism has been unclear. We train machine-learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on-the-fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li-ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large-amplitude vibrations of Cl – ions in TaCl 6 – enable the fast Li-ion transport by allowing dynamic breaking and reforming of Li–Cl bonds across the space in between the TaCl 6 – octahedra. We term this process the dynamic-monkey-bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all-solid-state lithium batteries.

25 ENERGY STORAGE↗

DOE OSTI2023

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

Evolutionary Computing for Low-thrust Navigation

The development of new mission concepts requires efficient methodologies to analyze, design and simulate the concepts before implementation. New mission concepts are increasingly considering the use of ion thrusters for fuel-efficient navigation in deep space. This paper presents parallel, evolutionary computing methods to design trajectories of spacecraft propelled by ion thrusters and to assess the trade-off between delivered payload mass and required flight time. The developed methods utilize a distributed computing environment in order to speed up computation, and use evolutionary algorithms to find globally Pareto-optimal solutions. The methods are coupled with two main traditional trajectory design approaches, which are called direct and indirect. In the direct approach, thrust control is discretized in either arc time or arc length, and the resulting discrete thrust vectors are optimized. In the indirect approach, a thrust control problem is transformed into a costate control problem, and the initial values of the costate vector are optimized. The developed methods are applied to two problems: 1) an orbit transfer around the Earth and 2) a transfer between two distance retrograde orbits around Europa, the closest to Jupiter of the icy Galilean moons. The optimal solutions found with the present methods are comparable to other state-of-the-art trajectory optimizers and to analytical approximations for optimal transfers, while the required computational time is several orders of magnitude shorter than other optimizers thanks to an intelligent design of control vector discretization, advanced algorithmic parameterization, and parallel computing.

optimization↗

NASA NTRS2004

L2 Earth Atmosphere Observatory: Formation Guidance, Metrology and Control Synthesis

The Earth Observatory Formation at L2, a Lagrange libration point, is a unique large aperture (25 m diameter) space telescope concept that will improve the knowledge and understanding of dynamic, chemical and radiative mechanisms that cause changes in the atmosphere, and can lead to the development of models and techniques to predict short and long-term climate changes. The results of this concept definition study show that the telescope concept is feasible, and can have technology readiness in the 2020 time frame. Further advanced development in several subsystems is needed, such as higher efficiency Xenon ion thrusters with throttling, and optical quality large membrane mirror with active shape control. It presents an analysis and solution of guidance, sensing, control, and propulsion problems for a formation of two spacecraft on the Sun-Earth line in the neighborhood of the Sun-Earth L2 point, that observes Earth s atmosphere during continuous solar occultation by the Earth. A system architecture is described for the observatory, and its components that include unique mission specific metrology. The formation must follow a powered trajectory with strictly limited fuel use to observe solar occultation. A configuration of ion thrusters and reaction wheels for translation and attitude control is designed along with algorithms for orbit following and formation control. Simulation results of the orbital and formation dynamics are presented that verify performance of the control systems.

Earth Atmosphere Observatory Formation↗

DOE OSTI2022

Li-Ion Diffusion Correlations in LiAlGeO 4 : Quasielastic Neutron Scattering and Ab Initio Simulation

Here, we investigated the impact of Li stoichiometry and host flexibility on Li + diffusion processes in LiAlGeO 4 at the microscopic level using quasielastic neutron scattering (QENS) and ab initio molecular dynamics (AIMD) simulations. Using sufficiently long AIMD trajectories, we could simulate the observed QENS signal and identify the localized dynamics of Li in crystalline LiAlGeO 4 . Such information is vital to identify the bottleneck of diffusion processes and design materials for battery application. Our AIMD simulations in LiAlGeO 4 reveal that the Li + conductivity can be significantly improved by manipulating the Li stoichiometry and/or host flexibility via amorphization. We determined that excess Li stoichiometry enhances the Coulomb repulsion of neighboring Li sites and softens the host structure to enable faster Li + diffusion along the hexagonal c-axis. In the amorphous structure, random orientations of AlO 4 and GeO 4 polyhedral units create a wide distribution of intersite distances and significantly soften the host structure, greatly enhancing the Li + diffusion. The simulations are used to understand the nature of diffusion, especially the role of the host structure, the possible hopping pathways, and the diffusion behavior in various structures of LiAlGeO 4 .

localized dynamics↗

DOE OSTI2025

Structural Changes in Metal Chalcogenide Nanoclusters Associated with Single Heteroatom Incorporation

Atomically precise nanoclusters (NCs) are promising building blocks for designing materials and interfaces with unique properties. By incorporating heteroatoms into the core, the electronic and magnetic properties of NCs can be precisely tuned. To accurately predict these properties, density functional theory (DFT) is often employed, making the rigorous benchmarking of DFT results particularly important. In this study, we present a benchmarking approach based on metal chalcogenide NCs as a model system. We synthesized a series of bimetallic, iron-cobalt chalcogenide NCs [Co 6-x Fe x S 8 (PEt 3 ) 6 ] + (x = 0-6) (PEt = triethyl phosphine) and investigated the effect of heteroatoms in the octahedral metal chalcogenide core on their size and electronic properties. Using ion mobility-mass spectrometry (IM-MS), we observed a gradual increase in the collision cross section (CCS) with an increase in the number of Fe atoms in the core. DFT calculations combined with trajectory method CCS simulations successfully reproduced this trend, revealing that the increase in cluster size is primarily due to changes in metal-ligand bond lengths, while the electronic properties of the core remain largely unchanged. Moreover, this method allowed us to exclude certain multiplicity states of the NCs, as their CCS values were significantly different from those predicted for the lowest-energy structures. Here, this study demonstrates that gas-phase IM-MS is a powerful technique for detecting subtle size differences in atomically precise NCs, which are often challenging to observe using conventional NC characterization methods. Accurate CCS measurements are established as a benchmark for comparison with theoretical calculations. The excellent correspondence between experimental data and theoretical predictions establishes a robust foundation for investigating structural changes of transition metal NCs of interest to a broad range of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2018

3D Simulations of Ion Thruster Accelerator Grid Erosion Accounting for Charge Exchange Ion Space Charge

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

Chaplin, Vernon H.↗