Engineering Papers⌕ Search

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.

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

At least 73 records · Page 4

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).

NASA Center ARC↗

Monte Carlo simulation of charged particle impact on the satellites of Jupiter and Saturn

Ensemble-averaged lifetimes of Jovian and Saturnian trapped charged particles against loss by impact on natural satellites are presented. Precise values of these loss rates can be obtained from the known geometry of the satellite orbits, magnetic fields, and particle trajectory. These loss rates are needed as an ingredient of transport theory which is usually formulated in terms of the first and second adiabatic invariants and radial distance. The effects of satellite orbital inclination and eccentricity, ion charge to mass ratio and energy, equatorial pitch, and distance from the planet on the particle lifetime are discussed utilizing a tilted, but not offset, dipolar field geometry resulting in lifetime values expressed as functions of the first two adiabatic invariants and the equatorial radial distance from the planet.

Bell, E. V., II↗

The SPACE Computer Code for Analyzing the International Space Station Electrical Power System: Past, Present, and Future

The System Power Analysis for Capability Evaluation (SPACE) computer code was initially developed by NASA in 1988 to assess the Space Station Freedom electric power system and later adapted to support contractor electrical power system capability analyses for the International Space Station (ISS). Over time, the code has supported many efforts such as ISS redesign activities in the early 1990s, assessment of time-phased loads against power system operating limits for future ISS assembly flights (including Certification of Flight Readiness reviews by the ISS program office), and determining the optimum solar array gimbal positions while respecting keep-out zones which minimize both solar array contamination and structural loads. The code has been validated by comparisons with ISS on-orbit data in multiple validation episodes. Recent updates to the code include the incorporation of a Lithium-Ion battery model in addition to the Nickel Hydrogen battery model and modifications to the solar array degradation model to better match on-orbit test results. SPACE has also been extended beyond the ISS to include modeling of the Orion Multi-Purpose Crew Vehicle electrical power system (SPACE-MPCV) and Mars Surface Electrical Power Systems (MSEPS). Portions of SPACE were integrated with a trajectory code to form a Solar Electric Propulsion Simulation (SEPSim), which can be used for analyzing solar electric propulsion missions. In addition, SPACE methods and subroutines have been adapted to a multitude of other projects. This paper summarizes the initial code development and subsequent code utilization in the context of the overall ISS program development and on-orbit operations. Recent updates and results from the code are discussed, including preliminary analyses for the Orion power system.

Electrical Power System↗

New Modeling Approaches to Investigate Cell Signaling in Radiation Response

Ionizing radiation damages individual cells and tissues leading to harmful biological effects. Among many radiation-induced lesions, DNA double-strand breaks (DSB) are considered the key precursors of most early and late effects [1] leading to direct mutation or aberrant signal transduction processes. In response to damage, a flow of information is communicated to cells not directly hit by the radiation through signal transduction pathways [2]. Non-targeted effects (NTE), which includes bystander effects and genomic instability in the progeny of irradiated cells and tissues, may be particularly important for space radiation risk assessment [1], because astronauts are exposed to a low fluence of heavy ions and only a small fraction of cells are traversed by an ion. NTE may also have important consequences clinical radiotherapy [3]. In the recent years, new simulation tools and modeling approaches have become available to study the tissue response to radiation. The simulation of signal transduction pathways require many elements such as detailed track structure calculations, a tissue or cell culture model, knowledge of biochemical pathways and Brownian Dynamics (BD) propagators of the signaling molecules in their micro-environment. Recently, the Monte-Carlo simulation code of radiation track structure RITRACKS was used for micro and nano-dosimetry calculations [4]. RITRACKS will be used to calculate the fraction of cells traversed by an ion and delta-rays and the energy deposited in cells in a tissue model. RITRACKS also simulates the formation of chemical species by the radiolysis of water [5], notably the .OH radical. This molecule is implicated in DNA damage and in the activation of the transforming growth factor beta (TGF), a signaling molecule involved in NTE. BD algorithms for a particle near a membrane comprising receptors were also developed and will be used to simulate trajectories of signaling molecules in the micro-environment and characterize autocrine and paracrine cell communication and signal transduction.

Plante, Ianik↗

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 three-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 40-50 V, owing to the photoelectron 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, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<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 Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of 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. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady 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.

Singh, N.↗

Particle acceleration due to shocks in the interplanetary field: High time resolution data and simulation results

Data were examined from two experiments aboard the Explorer 50 (IMP 8) spacecraft. The Johns Hopkins University/Applied Lab Charged Particle Measurement Experiment (CPME) provides 10.12 second resolution ion and electron count rates as well as 5.5 minute or longer averages of the same, with data sampled in the ecliptic plane. The high time resolution of the data allows for an explicit, point by point, merging of the magnetic field and particle data and thus a close examination of the pre- and post-shock conditions and particle fluxes associated with large angle oblique shocks in the interplanetary field. A computer simulation has been developed wherein sample particle trajectories, taken from observed fluxes, are allowed to interact with a planar shock either forward or backward in time. One event, the 1974 Day 312 shock, is examined in detail.

Kessel, R. L.↗

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

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

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

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↗

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

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↗

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↗

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

Features of Afterbody Radiative Heating for Earth Entry

Radiative heating is identified as a major contributor to afterbody heating for Earth entry capsules at velocities above 10 km/s. Because of rate-limited electron-ion recombination processes, a large fraction of the electronically-excited N and O atoms produced in the high temperature/pressure forebody remain as they expand into the afterbody region, which results in significant afterbody radiation. Large radiative heating sensitivities to electron-impact ionization rates and escape factors are identified. Ablation products from a forebody ablator are shown to increase the afterbody radiation by as much as 40%. The tangent-slab radiation transport approach is shown to over-predict the radiative flux by as much as 40% in the afterbody, therefore making the more computationally expensive ray-tracing approach necessary for accurate radiative flux predictions. For the Stardust entry, the afterbody radiation is predicted to be nearly twice as large as the convective heating during the peak heating phase of the trajectory. Comparisons between simulations and the Stardust Echelle observation measurements, which are shown to be dominated by afterbody emission, indicate agreement within 20% for various N and O lines. Similarly, calorimeter measurements from the Fire II experiment are identified as a source of validation data for afterbody radiation. For the afterbody calorimeter measurement closest to the forebody, which experiences the largest afterbody radiative heating component, the convective heating alone is shown to under-predict the measurement, even for the fullycatalytic assumption. Agreement with the measurements is improved with the addition of afterbody radiation. These comparisons with Stardust and Fire II measurements provide validation that the significant afterbody radiation values proposed in this work are legitimate.

Johnston, Christopher O.↗

Stormtime transport of ring current and radiation belt ions

A dynamical guiding-center simulation model is used to study the stormtime ion transport which leads to the formation of the ring current and diffusion in the radiation belts. Representative ions guiding-center motion in response to model storm-associated impulses in the convection electric field is traced for a range of ion mu values. The present numerical results are compared with previously formulated limiting idealization of particle transport in order to assess the limits of validity of these approximations. For ions having drift periods that exceed the duration of the main phase of the storm, their inward transport to form the stormtime ring current is appropriately described as direct convective access. For ions having drift periods comparable to the duration of the main phase of the storm, there is a transition between direct convective access and transport that resembles radial diffusion. Lower-energy ring-current ions at L of about 3 are freshly injected there from open adiabatic trajectories, whereas the higher-energy ring-current population consists of a mixture of freshly injected and previously trapped ions.

Chen, Margaret W.↗

Modeling of Orbital and Attitude Dynamics of Nanosatellites Controlled via Active Electrostatic Charging

The large-scale exploration of airless bodies, such as asteroids and moons, is gaining interest, however it is limited by mobility issues: the lack of atmosphere, low gravity, and unknown soil properties pose difficult challenges for many forms of traditional locomotion. The environment in proximity of these bodies is also electrically charged due to interactions with solar wind and UV radiation. The EGlider (Electrostatic Glider) concept would be able to overcome these mobility issues by leveraging the natural environment, allowing operations in close proximity of the surface, while enabling long duration missions by minimizing propellant consumption. The EGlider is an advanced concept for small satellite mobility and propulsion, which relies on the electric fields naturally present around airless bodies in order to generate forces and torques useful for maneuvering. It does so by extending electrically charged appendages, which enable it to electrostatically soar above the surface. By differentially charging its electrodes it can also produce torques to control its attitude. The charges are maintained by continuous active ion or electron emission from the spacecraft, in order to cancel out the neutralizing influx of charges from surrounding plasma. An investigation of the spacecraft-plasma interaction was carried out. This included studying the effect of electrode geometry and calculating the charge-to-mass ratios required to enable several mission scenarios. Long, thin wire electrodes were identified to be the most power-efficient and would allow power-to-weight ratios achievable with current nanosatellite technologies. High electrode potential represents the main limiting factor for the system design. In order to test the feasibility of active control by means of differential charging, a simple 2D interaction model was developed, and a feedback controller to stabilize the vehicle was tested in a simulation environment. The results confirmed that good performance could be obtained for both position and attitude control. Finally, a dedicated software was developed for future simulation and testing of control strategies for the EGlider. This software allows to study the trajectory and attitude of an arbitrarily configured spacecraft in the proximity of an arbitrarily defined main airless body. The spacecraft can be assembled from basic parts, each with specified electrical, mass and optical properties. Efficient models allow to calculate gravitational and electrical interactions with the rotating main body and the local plasma field, as well as solar radiation pressure effects. Control models can be implemented as simple plug-in functions and easily tested. The preliminary validation campaign showed good matching with the reference cases that have been analyzed.

Corpino, Sabrina↗

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗