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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Backstreaming ions from oblique earth bow shocks

One-dimensional, hybrid (particle ions, fluid electrons) numerical simulations of oblique shock structure are used to analyze the mechanism by which solar wind ions can be reflected off the earth's bow shock front and escape upstream. It is shown that no backstreaming ions are produced in the simulation over 50-90 degrees of theta(nB), the angle between the interplanetary magnetic field and the shock normal. Backstreaming ions are emitted for theta(nB) of 45 degrees or less, with characteristics in reasonable agreement with those of observed reflected ion streams. The simulated shock structure becomes rather turbulent as theta(nB) decreases from 45 to 30 degrees. Analysis of the backstreaming ion trajectories shows that all backstreaming ions are originally specularly reflected. The oblique shock structure and effects such as whistler formation and ion heating are discussed.

Leroy, M. M.↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Plasma electron hole kinematics. Hole tracking Particle-In-Cell simulation

The kinematics of a 1-D electron hole is studied using a novel Particle-In-Cell simulation code. A hole tracking technique enables us to follow the trajectory of a fast-moving solitary hole and study quantitatively hole acceleration and coupling to ions. We observe a transient at the initial stage of hole formation when the hole accelerates to several times the cold-ion sound speed. Artificially imposing slow ion speed changes on a fully formed hole causes its velocity to change even when the ion stream speed in the hole frame greatly exceeds the ion thermal speed, so there are no reflected ions. The behavior that we observe in numerical simulations agrees very well with our analytic theory of hole momentum conservation and the effects of “jetting.”

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Mapping of electrical potential distributions with charged particle beams

Methods for measuring electrostatic potentials on and near dielectric surfaces charged to several kilovolts are studied. Secondary emission from those charged dielectrics is measured. Candidates for potential measurement include the induced charge, from which potential is calculated; the trajectory endpoints of either high or low energy particles traversing the region near the surface; trajectory impact on the surface; and creating ions at points of interest near the surface. Some of the methods require computer simulations and iterative calculation if potential maps are to be generated. Several approaches are described and compared. A method using a half-cylinder as a test chamber and low-energy probing beams is adapted for the measurement of seconary emission.

Robinson, J. W.↗
NASA NTRS1994

The population of the magnetosphere by solar winds ions when the interplanetary magnetic field is northward

We have examined some possible entry mechanisms of solar wind ions into the magnetosphere by calculating the trajectories of thousands of non-interacting ions in the magnetic and electric fields from a three dimensional global magnetohydrodynamic (MHD) simulation of the magnetosphere and the magnetosheath, under northward interplanetary magnetic field (IMF) conditions. Particles, launched in the solar wind, entered the magnetosphere and formed the low latitude boundary layer (LLBL), plasma sheet and a region of trapped particles near the Earth. The densities and temperatures we obtained in these regions were realistic, with the exception of trapped particle densities. The dominant entry mechanism was convection into the magnetosphere on reconnecting field lines.

Richard, Robert L.↗
NASA NTRS2008

Simulating the Fate of an Ionospheric Mass Ejection

We report global ion kinetic (GIK) simulations of the 24-25 Sep 1988 storm, with all relevant ionospheric outflows including polar, auroral, and plasmaspheric winds. This storm included substantial periods of northward interplanetary magnetic field, but did develop a Dst of -200nT at its peak. The solar disturbance resulted from a coronal mass ejection that reached a peak dynamic pressure at the magnetosphere of 6.2 nPa, and produced a substantial enhancement of auroral wind oxygen outflow from the dayside, which has been termed an "ionospheric mass ejection" in an earlier paper. We use the LFM global simulation model to produce electric and magnetic fields in the outer magnetosphere, the Strangeway-Zheng outflow scalings with Delcourt ion trajectories to include ionospheric outflows, and the Fok-Obner inner magnetospheric model for the plasmaspheric and ring current response to all particle populations. We assess the combined contributions of heliospheric and geospheric plasmas to the ring current for this event.

Moore, Thomas E.↗
NASA NTRS1998

A Study of Transport in the Near-Earth Plasma Sheet During A Substorm Using Time-Dependent Large Scale Kinetics

In this study we investigate the transport of H+ ions that made up the complex ion distribution function observed by the Geotail spacecraft at 0740 UT on November 24, 1996. This ion distribution function, observed by Geotail at approximately 20 R(sub E) downtail, was used to initialize a time-dependent large-scale kinetic (LSK) calculation of the trajectories of 75,000 ions forward in time. Time-dependent magnetic and electric fields were obtained from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interaction with the solar wind and the interplanetary magnetic field (IMF) as observed during the interval of the observation of the distribution function. Our calculations indicate that the particles observed by Geotail were scattered across the equatorial plane by the multiple interactions with the current sheet and then convected sunward. They were energized by the dawn-dusk electric field during their transport from Geotail location and ultimately were lost at the ionospheric boundary or into the magnetopause.

El-Alaoui, M.↗
NASA NTRS2014

Titan's Plasma Environment : 3D Hybrid Kinetic Modeling of the TA Flyby and Comparison with CAPS-ELS and RPWS LP Observations

In this report we discuss the global plasma environment of the TA flyby from the perspective of 3D hybrid modeling. In our model the background, pickup, and ionospheric ions are considered as particles, whereas the electrons are described as a fluid. In homogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. Our modeling shows that mass loading of the background plasma (H(+), O(+)) by pick up ions H2(+), CH4(+) and N2(+) differs from theT9 encounter simulations when O(+) ions are not introduced into the background plasma. In our hybrid-modeling we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M=28 amu ions that were generated inside the ionosphere. Titan's interior is considered as a weakly conducting body. Special attention has been paid to comparing the simulated pickup ion density distribution with CAPS-ELS and with RPWS LP observations by the Cassini-Huygens spacecraft along theta trajectory. Our modeling shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures.

induced magnetospheres↗
NASA NTRS1983

Voyager 2 observations of energetic particle variations in the Ganymede wake region - A possible acceleration mechanism

Voyager 2's passage through the downstream corotation region of Ganymede found disturbances in the field and particle environment. Large fluctuations in the intensities and energy spectra of ions in the 0.1-to 4.0-MeV interval were also observed with the low energy charged particle experiment. All ion species were evidently affected, up through medium (C, N, O) nuclei. Electrons intensities, while disturbed, did not behave similarly to the ions. These effects are probably associated with the wake produced by the absorption of Jovian magnetospheric corotating plasma by Ganymede. A two-dimensional theoretical model for perturbed field in the wake region is proposed and it is shown by direct numerical simulation of exact particle trajectories that such perturbations are capable of modulating the charged particle energies to about the observed amount for the observed magnetic and plasma parameters. Analysis shows that the energization process is dependent on the pitch angles of the energetic ions and the process itself is nonadiabatic (the Larmor radius is approximately the size of the disturbance region).

Tariq, G.F.↗
DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.

dust↗
NASA NTRS2011

Modeling the Superstorm in November 2003

The superstorm on 20.21 November 2003 was the largest geomagnetic storm in solar cycle 23 as measured by Dst, which attained a minimum value of .422 nT. We have simulated this storm to understand how particles originating from the solar wind and ionosphere get access to the magnetosphere and how the subsequent transport and energization processes contribute to the buildup of the ring current. The global electromagnetic configuration and the solar wind H+ distribution are specified by the Lyon-Fedder-Mobarry (LFM) magnetohydrodynamics model. The outflow of H+ and O+ ions from the ionosphere are also considered. Their trajectories in the magnetosphere are followed by a test-particle code. The particle distributions at the inner plasma sheet established by the LFM model and test-particle calculations are then used as boundary conditions for a ring current model. Our simulations reproduce the rapid decrease of Dst during the storm main phase and the fast initial phase of recovery. Shielding in the inner magnetosphere is established at early main phase. This shielding field lasts several hours and then breaks down at late main phase. At the peak of the storm, strong penetration of ions earthward to L shell of 1.5 is revealed in the simulation. It is surprising that O+ is significant but not the dominant species in the ring current in our calculation for this major storm. It is very likely that substorm effects are not well represented in the models and O+ energization is underestimated. Ring current simulation with O+ energy density at the boundary set comparable to Geotail observations produces excellent agreement with the observed symH. As expected in superstorms, ring current O+ is the dominant species over H+ during the main to mid-recovery phase of the storm.

Fok, Mei-Ching↗
NASA NTRS2008

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗
DOE OSTI2021

HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017) (and subsequently titled) Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020) (Final Technical Report)

This is the Final Technical Report for the DoE Measurement Innovation (previously Diagnostic Development) grant DE-SC0006077 titled HIBP Applications to Advance Understanding of Plasma Transport Physics (2011-2017), and subsequently titled Development of Beam Measurements to Advance Understanding of Transients and Improve Validation (2017-2020). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 May 2011–17 July 2020. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that will improve understanding of transport, magnetic equilibrium, and electric fields in magnetically confined plasmas. We have developed novel hardware and advanced techniques to expand beam-based diagnostic capabilities, and extend measurements to new operating regimes. Some have been realized using a Heavy Ion Beam Probe (HIBP) having traditional features, while other applications have joined our innovations with subsets of advantageous HIBP features. This work has established smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key Measurement Innovations include: Computer Models and Simulations of HIBP Applications-Computer modeling and simulating of HIBP applications is central to predicting interaction of beam particles with plasmas, designing diagnostic systems, and anticipating measurement characteristics. We have extended the capability and improved measurement fidelity of the HIBP diagnostic on the Madison Symmetric Torus (MST) reversed field pinch (RFP); simulated feasibility of HIBP operation in the Helically Symmetric eXperiment (HSX) optimized stellarator; and investigated the plausibility of HIBP measurements in the ASDEX Upgrade tokamak. We have also developed a new technique that uses ion optics to model a finite phase-space beam. It yields more realistic (than traditional method) estimates of sample volume characteristics, which influence measurement resolution and sensitivity. Extension of Diagnostic Capabilities and Measurements on the MST RFP-Capabilities have been extended using the first and only HIBP installation on an RFP. Challenges associated with operation of the diagnostic on MST include the three-dimensional nature of particle trajectories, temporal topology and amplitude changes in the plasma equilibrium, and strong particle and radiative emission from the plasma. We have addressed these challenges through development of detailed and higher precision diagnostic simulations, calibrated subtraction of noise to resolve secondary ion signals, and implementation of hardware that enables higher fidelity measurement of plasma fluctuations. Development of New Hardware and Measurement Techniques - We have advanced beam-based measurement capabilities via development of new detectors and techniques to determine the poloidal magnetic flux ψ in the plasma, and enable placement of detectors close to the plasma and thus allow a smaller and less expensive diagnostic. We have demonstrated measurement of ψ, an ability that was made possible by development of highly effective noise reduction techniques. We also modeled the effect of non-ideal and instrumentation effects by developing a virtual (simulated) diagnostic. Application of beam-based diagnostics on various devices and multiple magnetic configurations (e.g. stellarator, tokamak, and reversed field pinch) enables investigation of critical physics issues and measurements resulting in the broad parameter space data needed to test and validate theory & modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Use of the thin sheath approximation for obtaining ion temperatures from the ISEE 1 limited aperture RPA

A procedure for analyzing low-energy (less than approximately 100 eV) ion data from the plasma composition experiment on ISEE 1 is set forth. The method is based on a derived analytic expression for particle flux to a limited aperture retarding potential analyzer (RPA) in the thin sheath approximation, which makes allowance for some effects of a charged spacecraft on plasma particle trajectories. Calculations using simulated data are employed in testing the efficacy and accuracy of the technique. On the basis of an analysis of these calculation results and the mathematical model, the method is seen as being able to provide accurate ion temperatures from all good plasmaspheric RPA data. It is noted that corresponding densities and spacecraft potentials should be accurate when spacecraft potentials are negative but that they are subject to error for positive spacecraft potentials, particularly when ion Mach numbers are much less than 1. An analysis of data from a representative ISEE 1 pass produces a plasmasphere temperature profile that is consistent in overall structure with previous observations.

Comfort, R. H.↗
NASA NTRS1997

Ion Sources and Acceleration Mechanisms Inferred from Local Distribution Functions

This study investigates the sources of the ions up the complex and nonisotropic H(+) velocity distribution functions observed by the Geotail spacecraft on May 23, 1995, in the near-Earth magnetotail region and recently reported by Frank et al. [1996]. A distribution function observed by Geotail at -10 R(sub E) downtail is used as input for the large scale kinetic (LSK) technique to follow the trajectories of approximately 90,000 H(+) ions backward in time. Time-dependent magnetic and electric fields are taken from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interactions with appropriate solar wind and IMF conditions. The ion population described by the Geotail distribution function was found to consist of a mixture of particles originating from three distinct sources: the ionosphere, the low latitude boundary layer (LLBL), and the high latitude plasma mantle. Ionospheric particles had direct access along field lines to Geotail, and LLBL ions convected adiabatically to the Geotail location. Plasma mantle ions, on the other hand, exhibited two distinct types of behavior. Most near-Earth mantle ions reached Geotail on adiabatic orbits, while distant mantle ions interacted with the current sheet tailward of Geotail and had mostly nonadiabatic orbits. Ions from the ionosphere, the LLBL, and the near-Earth mantle were directly responsible for the well-separated, low energy structures easily discernible in the observed and modeled distribution functions. Distant mantle ions formed the higher energy portion of the Geotail distribution. Thus, we have been successful in extracting useful information about particle sources, their relative contribution to the measured distribution and the acceleration processes that affected particle transport during this time.

Ashour-Abdalla, M.↗
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.↗
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗