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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1992

The effect of induced charges on low-energy particle trajectories near conducting and semiconducting plates

The effect of the induced charge was found on particles less than 1 eV as they passed through simulated parallel, grounded channels that are comparable in dimension to those that are presently in space plasma instruments which measure the flux of low-energy ions. Applications were made to both conducting and semiconducting channels that ranged in length from 0.1 to 50 mm and in aspect ratio from 1 to 100. The effect of the induced charge on particle trajectories from simple straight lines. Several configurations of channel aspect ratio and detector locations are considered. The effect is important only at very low energies with small dimensions.

Coffey, Victoria N.↗

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↗

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.↗

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↗

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↗

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.↗

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↗

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↗

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↗

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.↗

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.↗

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.↗

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↗

NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗

NASA NTRS1989

A three-dimensional numerical model of ionospheric plasma in the magnetosphere

A three-dimensional particle trajectory tracing in empirical models of the geoelectric and geomagnetic fields is used to study the ionospheric contribution to magnetospheric plasma. Various ionospheric outflows are examined and results on ion transport are presented in terms of density, composition, and energy. Results are presented for two opposite magnetospheric configurations, ground state and storm phases. An estimate of the contribution of ionospheric O(+) to the hot plasma sheet is given. The simulation results are compared with observational data.

Delcourt, D. C.↗