Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle 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 NTRS1995

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

Parallel implementation of a particle simulation for modeling rarefied gas dynamic flow

When the conditions of flow are rarefied and hypersonic, a more suitable alternative to the use of the Navier-Stokes equations for developing a numerical solution is the Direct Simulation Monte Carlo method (DSMC), a method of simulation which employs a large number of particles in modeling a rarefied gas. The performance of a parallel DSMC code developed for the Intel iPSC/860 Touchstone Gamma prototype computer is studied and the scaleup is found to be very nearly over the range of 16-128 processors.

Fallavollita, M. A.↗
NASA NTRS1992

A hybrid Zakharov particle simulation of ionospheric heating

A 1D hybrid simulation model incorporating one of the Zakharov equations for the high-frequency waves and a particle-in-cell simulation of the ions is described and applied to ionospheric heating using realistic parameters. Results from the hybrid simulation are compared with a Zakharov simulation that incorporates a phenomenological model of ion damping. Both the hybrid and Zakharov simulations predict the formation of solitonlike waves. The early time behavior of the two simulations is different due to the high noise level in the hybrid simulation, but the late time behavior is quite similar.

Clark, K. L.↗
NASA NTRS1987

Particle simulation of auroral double layers

Work on the simulation of auroral double layers (DLs) with realistic particle-in-cell models is presented. An early model simulated weak DLs formed in a self-consistent circuit but under conditions subject to the ion-acoustic instability. Recent work has focused on strong DLs formed when currentless jets are injected into a dipole magnetic field.

Smith, Bruce L.↗
NASA NTRS2020

Dynamical Evolution of Simulated Particles Ejected from Asteroid Bennu

In early 2019, the OSIRIS-REx spacecraft discovered small particles being ejected from the surface of the near-Earth asteroid Bennu. Although they were seen to be ejected at slow speeds, on the order of tens of cm/s, a number of particles were surprisingly seen to orbit for multiple revolutions and days, which requires a dynamical mechanism to quickly and substantially modify the orbit to prevent re-impact upon their first periapse passage. This paper demonstrates that, based on simulations constrained by the conditions of the observed events, the combined effects of gravity, solar radiation pressure, and thermal radiation pressure from Bennu can produce many sustained orbits for ejected particles. Furthermore, the simulated populations exhibit two interesting phenomena that could play an important role in the geophysical evolution of bodies such as Bennu. First, small particles (< 1 cm radius) are preferentially removed from the system, which could lead to a deficit of such particles on the surface. Second, re-impacting particles preferentially land near or on the equatorial bulge of Bennu. Over time, this can lead to crater in-filling and growth of the equatorial radius without requiring landslides.

Jay W. McMahon↗
NASA NTRS1993

Density-conserving shape factors for particle simulations in cylindrical and spherical coordinates

It is established that a particle-to-grid weighting which is based on both conventional particle-in-cell (PIC) and cloud-in-cell (CIC) shape factors will yield nonuniform grid densities, even in the case of uniform particle distributions, in both cylindrical and spherical coordinates. In view of this, alternative density-conserving weighting schemes are discussed which encompass modified PIC and CIC weighting. These alternative shape factors are mathematically no more complicated than the conventional nondensity-conserving expressions.

Ruyten, Wilhelmus M.↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Buffer-IPyC separation process in TRISO fuel particles simulated with Bison code

During High Temperature Gas-cooled Reactor (HTGR) operation, tristructural isotropic (TRISO) coated-particle fuel undergoes irradiation-induced changes in morphology and thermomechanical properties. Experimental results from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program show, among other things, the mechanism of gap formation between the buffer and inner pyrolytic carbon (IPyC) layers, which could be explored further via computational simulations using the Bison code. Two simulation models were developed, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted, and the debonding enabled model, where the gap between those layers is created. The inputs of the simulated models are based on the irradiation conditions from the AGR-1 experiment. The research included simulations on spherical and aspherical fuel types. Under the specific temperatures and fluences of the AGR-1 irradiation experiment, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. The results obtained for spherical and aspherical fuel are nearly identical. Finally, performed simulations match conclusions observed from the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the TRISO fuel. Future simulations will include the validation of performed research and comparison of the results between Bison and PARFUME codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1990

Plasma particle simulation of electrostatic ion thrusters

Charge exchange collisons between beam ions and neutral propellant gas can result in erosion of the accelerator grid surfaces of an ion engine. A particle in cell (PIC) is developed along with a Monte Carlo method to simulate the ion dynamics and charge exchange processes in the grid region of an ion thruster. The simulation is two-dimensional axisymmetric and uses three velocity components (2d3v) to investigate the influence of charge exchange collisions on the ion sputtering of the accelerator grid surfaces. An example calculation has been performed for an ion thruster operated on xenon propellant. The simulation shows that the greatest sputtering occurs on the downstream surface of the grid, but some sputtering can also occur on the upstream surface as well as on the interior of the grid aperture.

Peng, Xiaohang↗
NASA NTRS1995

Dawn-dusk asymmetries in the low-latitude boundary layer arising from the Kelvin-Helmholtz instability: A particle simulation

Along the low-latitude boundary layer (LLBL), the Kelvin-Helmholtz (K-H) instability can provide a means for injection of solar wind plasma across closed field lines of the magnetosphere. A fully electromagnetic, two-dimensional (three-velocity) particle code is used to investigate dawn-dusk asymmetries that can arise from velocity differences due to gradient drifts and electric field gradients, and ion acceleration across the narrow field transition layers in the dawn and dusk flanks. The model includes the dawn and dusk sides of the LLBL simultaneously in a slab geometry, incorporating seperate populations of ions and electrons for the magnetosphere and the magnetosheath plasmas. We report several effects: (1) asymmetries in the observed morphology of the turbulent structures, with familiar fluidlike vortex formation on the duskside, but tongues of magnetosheath plasma penetrating into the magnetosphere on the dawnside; (2) the formation of discrete current layers, characterized by strong currents and sharp gradients in the magnetic field, and discrete charge layers, having net charges and constant, weaker currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial boundary layer width; and (4) enhancement of the dawn-to-dusk electric field as magnetosheath particles and momenta are transported across the magnetopause.

Wilber, M.↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Particle simulations in magnetospheric plasmas

In view of the recent remarkable advancement of computer technology and simulation software, simulation studies are one of the most powerful academic tools for establishment of quantitative space physics and modelling of our space environment. The complex nature encountered in space plasma physics has motivated considerable development in computer simulations, which have played an essential role in the development of space plasma theory. This report describes research undertaken to understand physical processes involved in plasma waves observed in the magnetospheric plasmas, and associated nonlinear phenomena such as heating, diffusion, and acceleration of particles due to excited waves. The research explains and clarifies the observational data both qualitatively and quantitatively.

Nishikawa, Ken-Ichi↗
NASA NTRS2020

Lunar Dust Simulant Particle Adhesion on Copolyimide Alkyl Ethers

Mitigation of lunar dust contamination is one of the greatest challenges to be overcome to realize a sustained lunar presence. Likely solutions will integrate active mitigation strategies, requiring the input of external energy, and passive materials, exhibiting an intrinsic resistance to lunar dust adhesion. In this work, a series of copolyimide alkyl ethers containing per fluorinated side-chains were generated to evaluate the influence surface modification agents have on surface chemical, topographical, and mechanical properties. An expanded testing protocol to characterize the adhesion interaction between lunar dust simulant and the copolyimide substrate was carried out. The interfacial adhesion strength was in-situ measured by a custom-built particulate adhesion instrument, utilizing a sonic wand. Surface mechanical properties were characterized by nanoindentation, utilizing the continuous stiffness measurement approach. A nominal presence of surface modifying agents, 1 wt%, resulted in a six-fold reduction in adhesion strength of the interface. A strong inverse correlation between the adhesion strength and Young’s modulus of the substrate was identified. The reduction was attributed to a synergistic interaction between the surface energy, surface roughness, and modulus of the copolyimide alkyl ethers film.

Lunar Dust, Particulate Adhesion, Surface Migratio↗
DOE OSTI2022

Lagrangian particle simulation of hydrogen pellets and SPI into runaway electron beam in ITER

Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [Samulyak et al., Nucl. Fusion, 61(4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single-fragment ablation quantifies the influence of various factors, in particular, the impact ionization by runaway electrons and cross-field transport models, on the dynamics of ablated plasma and its penetration into the runaway beam. Simulations of SPI performed using different numbers of pellet fragments study the formation and evolution of the ablation clouds and their large-scale dynamics in ITER. In conclusion, the penetration depth of the ablation clouds is found to be of the order of 50 cm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification of gyrokinetic particle simulations of neoclassical tearing modes in fusion plasmas

Abstract The ability to simulate neoclassical tearing modes (NTMs) in the gyrokinetic toroidal code (GTC) has been developed and verified, in which ions are treated with a gyrokinetic model and electrons are treated as a resistive fluid. The simulation results demonstrate that the neoclassical bootstrap current effect can destabilize an otherwise stable classical tearing mode. In the cylindrical geometry, GTC simulations in the magnetohydrodynamic limit show quantitative agreement with the modified Rutherford theory, both in terms of the scaling law in the small island limit and in terms of the saturation level and pressure flattening effect in the large island limit. The toroidal effects are slightly destabilizing for the NTM, while the kinetic effects of thermal ions are stabilizing for the NTM and increase its excitation threshold.

Physics↗