Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distributions”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
DOE OSTI2020

Simulation of Alfvénic avalanche onset in NSTX

In some tokamak discharges, a number of Alfvén eigenmodes (AEs) have been observed to cause a large-scale collapse of the high energy particle distribution, a phenomenon referred to as an avalanche. We examine the necessary conditions for an avalanche using the available experimental information from NSTX on the equilibrium and mode properties for two cases, one with a measurable but benign AE activity and one with an AE activity leading up to an avalanche. To produce an avalanche, the modes present in the discharge must possess resonances that can overlap with a modest increase in instability magnitude, providing a path to global particle stochastic motion. In this paper, we find that the modes present in the avalanche-free discharge do not provide such a path even at a very large amplitude. During the discharge which subsequently produces an avalanche, the high energy population is growing and the Alfvén frequency is dropping due to increasing density, and we find that both these changes, producing a small increased drive or an increased resonance width for the Alfvén modes, can lead in this case to uncontrolled mode growth and large-scale beam particle loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the collisional damping of plasma velocity space instabilities

For plasma velocity space instabilities driven by particle distributions significantly deviated from a Maxwellian, weak collisions can damp the instabilities by an amount that is significantly beyond the collisional rate itself. This is attributed to the dual role of collisions that tend to relax the plasma distribution toward a Maxwellian and to suppress the linearly perturbed distribution function. The former effect can dominate in cases where the unstable non-Maxwellian distribution is driven by collisionless transport on a timescale much shorter than that of collisions, and the growth rate of the ideal instability has a sensitive dependence on the distribution function. The whistler instability driven by electrostatically trapped electrons is used as an example to elucidate such a strong collisional damping effect of plasma velocity space instabilities, which is confirmed by first-principles kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Post-Newtonian kinetic theory

A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling Solution Drying by Moving a Liquid-Vapor Interface: Method and Applications

A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe “polymer-on-top” stratification, similar to that found previously with an explicit solvent model. Furthermore, “polymer-on-top” is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the “small-on-outside” stratification mechanism at fast evaporation rates. “Large-on-outside” stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pseudorapidity distributions of charged particles in pp($\bar p$), p(d)A and AA collisions using Tsallis thermodynamics

The pseudorapidity distributions of charged particles measured in p + p($\overline{\text{p}}$) collisions for energies ranging from $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=23.6$ GeV to 13 TeV and A + A collisions at RHIC and LHC are investigated in the fireball model with Tsallis thermodynamics. We assume that the rapidity axis is populated with fireballs following q-Gaussian distribution and the charged particles follow the Tsallis distribution in the fireball. We also extend the fireball model to asymmetric collision systems, i.e. d + Au collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=200$ GeV and p + Pb collisions at $\sqrt{{s}_{\mathrm{N}\mathrm{N}}}=5.02$ TeV, by taking into account the asymmetric geometry configuration. The model can fit the experimental data well for all the collision systems and centralities investigated. The collision energy and centrality dependence of the model parameters for the symmetric (asymmetric) collision system, i.e. the central position y 0 (y 0a , y 0A ) and its width σ (σ a , σ A ) of the fireball distribution, are also investigated and discussed. Furthermore, our results suggest that the fireball model with Tsallis thermodynamics can be used as a universal framework for the pseudorapidity distributions of charged particles in high energy collisions at RHIC and LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Thermal Stability of Silica for Application in Thermal Energy Storage

Thermal energy storage (TES) systems have enabled concentrating solar power (CSP) to remain competitive in the modern energy mix by providing economical load shifting grid services and firming up intermittent solar resource. Free from siting constraints, TES also shows promise as an economical alternative to traditional pumped-storage hydropower (PSH) and compressed air energy storage (CAES). As potential thermal energy storage media, some solid particles demonstrate stability over wide temperature ranges which allows for increased sensible energy storage density and is essential in achieving low-cost storage. Silica sand, in the form of a-quartz, is one such candidate. This work presents a brief review of relevant silica thermophysical properties and further investigates the thermal stability of silica particles as a candidate TES media by subjecting them to two different thermal tests: (1) a 500-hour thermal treatment at 1200 degrees C under varied atmospheres; and (2) cycling 25, 50, and 100 times between 300 degrees C and 1200 degrees C. For both tests, particle stability is examined by means of pre- and post-treatment Mie scattering. An additional XRD analysis is conducted for the 500-hour treatment in air. Results indicate limited changes in both particle distribution and crystallographic structure which is promising for the application as solid particle media for thermal energy storage.

concentrating solar power↗
DOE OSTI2020

Thermal Stability of Silica for Application in Thermal Energy Storage

Free from siting constraints, thermal energy storage (TES) shows promise as an economical alternative to traditional pumped-storage hydropower (PSH) and compressed air energy storage (CAES). As potential thermal energy storage media, many solid particles demonstrate stability over wide temperature ranges which allows for increased sensible energy storage density and is essential in achieving low-cost storage. Silica sand, in the form of a-quartz, is one such candidate. This work presents a brief review of relevant silica thermophysical properties and further investigates the thermal stability of silica particles as a candidate TES media by subjecting them to two different thermal campaigns: (1) a 500-hour thermal treatment at 1200°C under varied atmospheres; and (2) cycling 25, 50, and 100 times between 300°C and 1200°C. For both campaigns, particle stability is examined by means of pre- and post-treatment Mie scattering. An additional XRD analysis is conducted for the 500-hour treatment in air. Results indicate limited changes in both particle distribution and crystallographic structure which is promising for the application as solid particle media for thermal energy storage.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2020

Thermal Stability of Silica for Application in Thermal Energy Storage: Preprint

Free from siting constraints, thermal energy storage (TES) shows promise as an economical alternative to traditional pumped-storage hydropower (PSH) and compressed air energy storage (CAES). As potential thermal energy storage media, many solid particles demonstrate stability over wide temperature ranges which allows for increased sensible energy storage density and is essential in achieving low-cost storage. Silica sand, in the form of a-quartz, is one such candidate. This work presents a brief review of relevant silica thermophysical properties and further investigates the thermal stability of silica particles as a candidate TES media by subjecting them to two different thermal campaigns: (1) a 500-hour thermal treatment at 1200°C under varied atmospheres; and (2) cycling 25, 50, and 100 times between 300°C and 1200°C. For both campaigns, particle stability is examined by means of pre- and post-treatment Mie scattering. An additional XRD analysis is conducted for the 500-hour treatment in air. Results indicate limited changes in both particle distribution and crystallographic structure which is promising for the application as solid particle media for thermal energy storage.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2026

QCD Theory and Transverse Momentum Dependent Physics

A major challenge for QCD theorists is to understand how to different categories of QCD physics contribute to a single observable. This is the goal of QCD factorization theory. Today, QCD physicists are tasked with describing increasingly detailed properties of transverse momentum dependent particle distributions, including their shapes and normalizations, their variation with energy and other kinematical variables, the identity of particles involved, correlations with spin, and the relationship between transverse momentum distributions in different processes. Embedded within that broad effort are specific fundamental questions at the heart of QCD theory. Namely, how does one give precise, quantitative meaning to quark and gluon concepts in hard scattering events, and what are the limits of those concepts? In the past two decades, there has been a renaissance in TMD QCD, with a broad spectrum of motivations driving the study of TMD observables. TMD QCD studies are in a transformative stage; subfields that traditionally progressed nearly independently of one other, such as hadronic structure studies and the study of perturbative QCD gluon radiation in high-energy collisions, are becoming intertwined. It is thus becoming increasingly clear that unified approaches -- those that incorporate all the relevant phenomena into a single formalism rooted in QCD fundamentals -- will be critical to the mission to expand the current understanding of QCD nuclear/hadronic dynamics. TMD QCD factorization theorems and their generalizations will be central to these efforts. The factorization theorems, in addition, provide the practical prescription for applying specific calculations to phenomenology. The goal of this project was to unify and push the boundaries of traditional TMD research and prepare it for use in future phenomenological applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Friction stir processing on a strontium modified, thin-wall, vacuum-assisted high-pressure die-cast Aural-5 alloy to improve tensile and fatigue performance

Here, this study explores the application of friction stir processing (FSP) to enhance the material properties of Sr-modified Aural-5 alloy, with a focus on improved tensile and fatigue properties. Aural-5 is a well-known vacuum-assisted high-pressure die-cast (HPDC) Al-Si7-Mg alloy used in the automotive industry to reduce vehicle weight, enhance fuel efficiency, and lower carbon emissions. This alloy modifies its material chemistry with Sr for fine fibrous networks of eutectic silicon and manganese (Mn) to reduce die soldering. It has significantly less iron (Fe) content resulting in the elimination of detrimental needle-shaped Fe-bearing ß-phase intermetallic and improving ductility. The initial microstructure of as-received HPDC Aural-5 exhibits shrinkage porosity in the middle section, a dendritic microstructure with fibrous Al-Si eutectic colonies, a shear-band structure beneath the die-wall, large dendritic externally solidified crystals (ESCs), needle-shaped Mg 2 Si phase and significant second-phase particulates. Some of those microstructural features, such as porosity, ESCs, needle-shaped Mg 2 Si phase, and large second-phase particles, serve as initiation sites for cracks under mechanical loading, resulting in adverse effects on tensile properties, particularly ductility. FSP effectively transforms the microstructure into a wrought configuration with uniform particle distribution by eliminating porosity and disintegrating dendrites, eutectic colonies, ESCs, second-phase particles, and shear-band structures. FSP-driven microstructure modification enhances yield strength and tensile ductility by ~30% and ~35%, respectively. The fatigue life of the material in a bending mode configuration (stress ratio R = 0.1) after FSP exhibits enhancements ranging from 2.0 to 3.9 times that of the original HPDC Aural-5 alloy, depending on the applied stress level.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Particle resonances in stellarators

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the particle distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through an avalanche. The existence of a mode particle resonance depends on the properties of the equilibrium and particle parameters, and their number, location, and density can vary with device design. Recently, the advent of more powerful computing capabilities and advanced theoretical understanding has led to the design of non-axisymmetric devices or stellarators, which could prove to be more advantageous than tokamaks. Stellarators have the advantage of being immune to major disruptions because of the very low plasma current. One of the problems shared by both types of devices is the existence of resonances in particle orbits, which can lead to large amplitude high frequency instabilities and subsequent induced particle loss. We examine the number of resonances, their location, and dependence on particle energy for some stellarator designs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Skopi: a simulation package for diffractive imaging of noncrystalline biomolecules

X-ray free-electron lasers (XFELs) have the ability to produce ultra-bright femtosecond X-ray pulses for coherent diffraction imaging of biomolecules. While the development of methods and algorithms for macromolecular crystallography is now mature, XFEL experiments involving aerosolized or solvated biomolecular samples offer new challenges in terms of both experimental design and data processing. Skopi is a simulation package that can generate single-hit diffraction images for reconstruction algorithms, multi-hit diffraction images of aggregated particles for training machine learning classifiers using labeled data, diffraction images of randomly distributed particles for fluctuation X-ray scattering algorithms, and diffraction images of reference and target particles for holographic reconstruction algorithms. Skopi is a resource to aid feasibility studies and advance the development of algorithms for noncrystalline experiments at XFEL facilities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Pseudo-rapidity distributions of charged particles in asymmetric collisions using Tsallis thermodynamics

Abstract The pseudo-rapidity distributions of the charged particles produced in the asymmetric collision systems p+Al, p+Au and$$^3$$ 3 He+Au at$$\sqrt{s_\text {NN}} = 200$$ s NN = 200 GeV are evaluated in the framework of a fireball model with Tsallis thermodynamics. The fireball model assumes that the experimentally measured particles are produced by fireballs following the Tsallis distribution and it can effectively describe the experimental data. Our results as well as previous results for d+Au collisions at$$\sqrt{s_\text {NN}} = 200$$ s NN = 200 GeV and p+Pb collisions at$$\sqrt{s_\text {NN}} = 5.02$$ s NN = 5.02 TeV validate that the fireball model based on Tsallis thermodynamics can provide a universal framework for pseudo-rapidity distribution of the charged particles produced in asymmetric collision systems. We predict the centrality dependence of the total charged particle multiplicity in the p+Al, p+Au and$$^3$$ 3 He+Au collisions. Additionally, the dependences of the fireball model parameters ($$y_{0a}$$ y 0 a ,$$y_{0A}$$ y 0 A ,$$\sigma _{a}$$ σ a and$$\sigma _{A}$$ σ A ) on the centrality and system size are studied.

Nuclear Science & Technology↗