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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Small Particle Response to Fluid Motion using Tethered Particles to Simulate Microgravity

This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used for finding the position with an accuracy of a few microns. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. The electronic holograms can be downlinked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Trolinger, J. D.↗
NASA NTRS2002

Small Particle Response to Fluid Motion using Tethered Particles to Simulate Microgravity

This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used for finding the position with an accuracy of a few microns. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. The electronic holograms can be downlinked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Trolinger, J. D.↗
DOE OSTI2024

Skyrmion soliton motion on periodic substrates by atomistic and particle-based simulations

Abstract We compare the dynamical behavior of magnetic skyrmions interacting with square and triangular defect arrays just above commensuration using both atomistic and particle-based models. Under applied drives, the initial motion is a kink traveling through the pinned skyrmion lattice. For the square defect array, both models agree well and show a regime in which the soliton motion is locked along 45°. The atomistic model also produces locking of a soliton along 30°, which is absent in the particle-based model. For the triangular defect array, the atomistic model exhibits 30° soliton motion over a wide region of external current values. In contrast, the particle-based model gives 45° soliton motion over a small range of external driving force values. The difference arises because the nondeforming particle model facilitates meandering skyrmion orbits while the deformable atomistic model enables stronger skyrmion-skyrmion interactions that reduce the meandering. Our results indicate that soliton motion through pinned skyrmion lattices on a periodic substrate is a robust effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Relativistic electron scattering by magnetosonic waves: Effects of discrete wave emission and high wave amplitudes

In this paper, we study relativistic electron scattering by fast magnetosonic waves. We compare results of test particle simulations and the quasi-linear theory for different spectra of waves to investigate how a fine structure of the wave emission can influence electron resonant scattering. We show that for a realistically wide distribution of wave normal angles θ (i.e., when the dispersion δθ≥0.5{sup °}), relativistic electron scattering is similar for a wide wave spectrum and for a spectrum consisting in well-separated ion cyclotron harmonics. Comparisons of test particle simulations with quasi-linear theory show that for δθ>0.5{sup °}, the quasi-linear approximation describes resonant scattering correctly for a large enough plasma frequency. For a very narrow θ distribution (when δθ∼0.05{sup °}), however, the effect of a fine structure in the wave spectrum becomes important. In this case, quasi-linear theory clearly fails in describing accurately electron scattering by fast magnetosonic waves. We also study the effect of high wave amplitudes on relativistic electron scattering. For typical conditions in the earth's radiation belts, the quasi-linear approximation cannot accurately describe electron scattering for waves with averaged amplitudes >300 pT. We discuss various applications of the obtained results for modeling electron dynamics in the radiation belts and in the Earth's magnetotail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Resolution of differences between collision number definitions in particle and continuum simulations

The use of the same symbol, Z, representing a 'collision number' for thermal relaxation, has led to confusion regarding its definition in the context of both continuum and particle simulations. Examination of the relaxation mechanics employed in particle simulations demonstrates that these definitions differ by a numerical factor that depends upon the intermolecular potential. Particle and continuum simulations employing appropriate definitions of Z lead to identical results during isothermal and adiabatic stationary relaxation.

Lumpkin, Forrest E., III↗
DOE OSTI2019

Vectorization of Random Number Generation and Reproducibility of Concurrent Particle Transport Simulation

Efficient random number generation with high quality statistical properties and exact reproducibility of Monte Carlo simulations are important requirements in many areas of computational science. VecRNG is a package providing pseudo-random number generation (pRNG) in the context of a new library VecMath. This library bundles up several general-purpose mathematical utilities, data structures, and algorithms having both SIMD and SIMT (GPUs) support based on VecCore.~Several state-of-the-art RNG algorithms are implemented as kernels supporting parallel generation of random numbers in scalar, vector, and Cuda workflows. In this report, we will present design considerations, implementation details, and computing performance of parallel pRNG engines on both CPU and GPU. Reproducibility of propagating multiple particles in parallel for HEP event simulation is demonstrated, using GeantV-based examples, for both sequential and fine-grain track-level concurrent simulation wor kflows. Strategies for efficient uses of vectorized pRNG and non-overlapping streams of random number sequences in concurrent computing environments is discussed as well.

Jun, S. Y.↗
DOE OSTI2013

SciDAC Center for Plasma Edge Simulation

This project with a total funding of $592,998 for six years has partially supported four postdoctoral researchers at the University of California, Irvine (UCI). The UCI team has formulated electrostatic and electromagnetic global gyrokinetic particle simulation models with kinetic electrons, implemented these models in the edge code XGC1, performed benchmark between GTC and XGC1, developed computational tools for gyrokinetic particle simulation in tokamak edge geometry, and initiated preparatory study of edge turbulence using GTC code. The research results has been published in 12 papers and presented at many international and national conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Comparison of dust charging between orbital-motion-limited theory and particle-in-cell simulations

The Orbital-Motion-Limited (OML) theory has been modified to predict the dust charge and the results were contrasted with the Whipple approximation [X. Z. Tang and G. L. Delzanno, Phys. Plasmas 21, 123708 (2014)]. To further establish its regime of applicability, here, the OML predictions (for a non-electron-emitting, spherical dust grain at rest in a collisionless, unmagnetized plasma) are compared with particle-in-cell simulations that retain the absorption radius effect. It is found that for large dust grain radius r d relative to the plasma Debye length λ D , the revised OML theory remains a very good approximation as, for the parameters considered (r d /λ D ≤ 10, equal electron and ion temperatures), it yields the dust charge to within 20% accuracy. This is a substantial improvement over the Whipple approximation. The dust collected currents and energy fluxes, which remain the same in the revised and standard OML theories, are accurate to within 15%–30%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

A fast integrated mobility spectrometer for rapid measurement of sub-micrometer aerosol size distribution, Part II: Experimental characterization

A Fast Integrated Mobility Spectrometer (FIMS) with a wide dynamic size range has been developed for rapid aerosol size distribution measurements. The design and model evaluation of the FIMS are presented in the preceding paper (Paper I), and this paper focuses on the experimental characterization of the FIMS. Monodisperse aerosol with diameter ranging from 8 to 600 nm was generated using Differential Mobility Analyzer (DMA), and was measured by the FIMS in parallel with a Condensation Particle Counter (CPC). The mean particle diameter measured by the FIMS is in good agreement with the DMA centroid diameter. Comparison of the particle concentrations measured by the FIMS and CPC indicates the FIMS detection efficiency is essentially 100% for particles with diameters of 8 nm or larger. For particles smaller than 20 nm or larger than 200 nm, FIMS transfer function and resolution can be well represented by the calculated ones based on simulated particle trajectories in the FIMS. For particles between 20 and 200 nm, the FIMS transfer function is boarder than the calculated, likely due to non-ideality of the electric field, including edge effects near the end of the electrode, which are not represented by the 2-D electric field used to simulate particle trajectories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Implementation of a Particle Accelerator Beam Dynamics Code on Multi-Node GPUs

Particle accelerators play an important role in a wide range of scientific discoveries and industrial applications. The self-consistent multi-particle simulation based on the particle-in-cell (PIC) method has been used to study charged particle beam dynamics inside those accelerators. However, the PIC simulation is time-consuming and needs to use modern parallel computers for high-resolution applications. In this paper, we implemented a parallel beam dynamics PIC code on multi-node hybrid architecture computers with multiple Graphics Processing Units (GPUs). We used two methods to parallelize the PIC code on multiple GPUs and observed that the replication method is a better choice for moderate problem size and current computer hardware while the domain decomposition method might be a better choice for large problem size and more advanced computer hardware that allows direct communications among multiple GPUs. Using the multi-node hybrid architectures at Oak Ridge Leadership Computing Facility (OLCF), the optimized GPU PIC code achieves a reasonable parallel performance and scales up to 64 GPUs with 16 million particles.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
DOE OSTI2015

Analytical approximation of transit time scattering due to magnetosonic waves

Abstract Recent test particle simulations have shown that energetic electrons traveling through fast magnetosonic (MS) wave packets can experience an effect which is specifically associated with the tight equatorial confinement of these waves, known as transit time scattering. However, such test particle simulations can be computationally cumbersome and offer limited insight into the dominant physical processes controlling the wave‐particle interactions, that is, in determining the effects of the various wave parameters and equatorial confinement on the particle scattering. In this paper, we show that such nonresonant effects can be effectively captured with a straightforward analytical treatment that is made possible with a set of reasonable, simplifying assumptions. It is shown that the effect of the wave confinement, which is not captured by the standard quasi‐linear theory approach, acts in such a way as to broaden the range of particle energies and pitch angles that can effectively resonate with the wave. The resulting diffusion coefficients can be readily incorporated into global diffusion models in order to test the effects of transit time scattering on the dynamical evolution of radiation belt fluxes.

Bortnik, J.↗
DOE OSTI2017

Numerical Test of Analytical Theories for Perpendicular Diffusion in Small Kubo Number Turbulence

In the literature, one can find various analytical theories for perpendicular diffusion of energetic particles interacting with magnetic turbulence. Besides quasi-linear theory, there are different versions of the nonlinear guiding center (NLGC) theory and the unified nonlinear transport (UNLT) theory. For turbulence with high Kubo numbers, such as two-dimensional turbulence or noisy reduced magnetohydrodynamic turbulence, the aforementioned nonlinear theories provide similar results. For slab and small Kubo number turbulence, however, this is not the case. In the current paper, we compare different linear and nonlinear theories with each other and test-particle simulations for a noisy slab model corresponding to small Kubo number turbulence. We show that UNLT theory agrees very well with all performed test-particle simulations. In the limit of long parallel mean free paths, the perpendicular mean free path approaches asymptotically the quasi-linear limit as predicted by the UNLT theory. For short parallel mean free paths we find a Rechester and Rosenbluth type of scaling as predicted by UNLT theory as well. The original NLGC theory disagrees with all performed simulations regardless what the parallel mean free path is. The random ballistic interpretation of the NLGC theory agrees much better with the simulations, but compared to UNLT theory the agreement is inferior. We conclude that for this type of small Kubo number turbulence, only the latter theory allows for an accurate description of perpendicular diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Properties of toroidal Alfvén eigenmode in DIII-D plasma

Here, linear properties of the toroidal Alfvén eigenmode (TAE) excited by energetic particles (EP) in a DIII-D tokamak experiment have been examined in global gyrokinetic particle simulations treating self-consistently kinetic effects of EP, thermal ions, and electrons. Simulation results of the TAE frequency and mode structure agree very well with the experimental measurements. The non-perturbative EP contribution induces a radial localization of the TAE mode structure, a break-down of mode radial symmetry, as well as a frequency dependence on the toroidal mode number. The simulations further exhibit the dependence of the growth rate and mode structure on EP pressure gradients. The in-out asymmetry of the mode structure and the experimental identification of the poloidal harmonics have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Simulation of particle velocity in a laser-produced tin plasma extreme ultraviolet source

In connection with fast heating in a laser produced plasma (LPP) extreme ultraviolet (EUV) source, the superheating behavior of bulk tin (Sn) at high heating rates is investigated. A constant temperature and pressure molecular dynamics simulation using modified Lennard-Jones and Coulomb potentials suitable for studying the liquid structure of Sn is employed in order to derive the caloric curves of the solid and liquid phases. The results have shown transient effects on the phase transitions. Superheating is observed during the melting and vaporizing processes. The velocity distribution of Sn particles against typical laser fluence in a LPP EUV light source has been numerically investigated using a simplified method including a one-dimensional, two-temperature, molecular dynamics, and steady-state ionization model. In the framework of our model, it was found that ejected Sn particles have a maximum velocity on the order of 10 to 40 km/s in plasma created using a nanosecond pre-pulse neodymium-doped yttrium aluminum garnet (Nd:YAG, 1.06 {mu}m) laser in EUV lithography experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Particle Tracking Simulation of Collective Modes. Parametric Landau Damping Off Coupling Resonance.

Employing Synergia simulations with the DMD method we investigate the Landau damping of space charge modes in bunched beams. The simulations reveal two instances of the parametric damping mechanism in bunched beams. The first example occurs in the proximity of coupling resonance and is due to the oscillation of particles’ amplitudes in the transverse plane. This oscillation modulates the modeparticle coupling with particle dependent trapping frequency. The second example is due to the modulation of the modeparticle coupling in one transverse plane by the oscillatory motion in the other plane.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Effects of trapped electrons on ion reflection in an oblique shock wave

A magnetosonic shock wave propagating obliquely to an external magnetic field can trap electrons and accelerate them to ultrarelativistic energies. The trapped electrons excite two-dimensional (2D) electromagnetic fluctuations with finite wavenumbers along the shock front. We study effects of the trapped electrons on ion motions through the 2D fluctuations. It is analytically shown that the fraction of ions reflected from the shock front is enhanced by the 2D fluctuations. This is confirmed by 2D (two space coordinates and three velocities) relativistic, electromagnetic particle simulations with full ion and electron dynamics and calculation of test ions in the electromagnetic fields averaged along the shock front. A comparison between 2D and one-dimensional electromagnetic particle simulations is also shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗