Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distribution functions”

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 127 records · Page 7

Microscopic calculation of fission product yields with particle-number projection

Fission fragments' charge and mass distribution is an important input to applications ranging from basic science to energy production or nuclear nonproliferation. In simulations of nucleosynthesis or calculations of superheavy elements, these quantities must be computed from models, as they are needed in nuclei where no experimental information is available. Until now, standard techniques to estimate these distributions were not capable of accounting for fine-structure effects, such as the odd-even staggering of the charge distributions. In this work, we combine a fully microscopic collective model of fission dynamics with a recent extension of the particle number projection formalism to provide the highest-fidelity prediction of the primary fission fragment distributions for the neutron-induced fission of 235 U and 239 Pu. Here, we show that particle-number projection is an essential ingredient to reproduce odd-even staggering in the charge yields and benchmark the performance of various empirical probability laws that could simulate its effect. This new approach also enables for the first time the realistic determination of two-dimensional isotopic yields within nuclear density functional theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic Energy Dissipation and γ-Ray Emission in Energetic Pulsars

Abstract Some of the most energetic pulsars exhibit rotation-modulated γ -ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power ( γ -ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower γ -ray efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗

Tracking of Marangoni driven motion during laser powder bed fusion

Laser Powder Bed Fusion (LPBF) can be used to create Functionally Graded Materials (FGMs) by selectively distributing property augmenting particles throughout the build matrix. With this method parts can be created with nonhomogeneous properties tuned to spatially varied, pre-determined design constraints. To unlock this capability, the fluid motion that occurs in the molten liquid phase must be thoroughly understood. More importantly, the actual tracking of motion during the LPBF process is needed. Further, a numerical model of heat transfer and fluid flow in LPBF process was created. Using results from the model the motion of particles starting at many different locations was tracked. The average change in position for these particles based on starting location was calculated. The effects of non-dimensional parameters on melt pool size were examined. These results were validated experimentally and compared to what is available in the literature.

42 ENGINEERING↗

Investigating the broadening phenomenon in two-particle correlations induced by gluon saturation

It has been found that the gluon density inside the proton grows rapidly at small momentum fractions. Quantum chromodynamics (QCD) predicts that this growth can be regulated by nonlinear effects, ultimately leading to gluon saturation. Within the color glass condensate framework, nonlinear QCD effects are predicted to suppress and broaden back-to-back angular correlations in collisions involving heavy nuclei. While suppression has been observed in various experiments in 𝑑/𝑝 + A collisions compared to 𝑝 + 𝑝 collisions, the predicted broadening remains unobserved. This study investigates the contributions of intrinsic transverse momentum (𝑘 𝑇 ), which is associated with saturation physics, as well as parton showers and transverse motion from fragmentation (𝑝$^{frag}_{𝑇}$), which are not saturation dependent, to the width of the correlation function. Our findings show that the non-saturation dependent effects, especially the initial-state parton shower and 𝑝$^{frag}_{𝑇}$, which occur independently of the collision system, smear the back-to-back correlation more than gluon saturation does, making the broadening phenomenon difficult to observe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Particle-in-cell simulations and passive bulk model of collisional capacitive discharge

Intermediate pressure (0.2–6.0 Torr), radio frequency capacitive discharges are of increasing importance in the thin film processing industry. However, existing low pressure (<0.1 Torr) discharge models are often used beyond their applicability in this pressure range. In this work, the authors performed one-dimensional particle-in-cell (PIC) simulations of a capacitive 2.5 cm gap argon discharge driven at 13.56 MHz with current density amplitude J 0 50 A/m 2 over a wide range of gas pressures p g = 0.04–20 Torr. For intermediate pressures, the PIC results showed that most of the ionization occurs near the plasma-sheath edges with very little occurring within the central bulk region. This led us to develop a “passive bulk” model of the discharge in which bulk electron heating does not contribute to ionization, and all the ionization is due solely to the sheath-heated electrons. The model assumes a constant ion-neutral mean free path λ i in the plasma bulk and either constant λ i or constant ion mobility μ i in the sheaths. The model is solved using an electron temperature iteration scheme. The PIC scalings of the discharge equilibrium parameters (e.g., electron density and temperature, sheath width, and voltage) with the input parameters p g and J 0 showed good agreement with those of the passive bulk model. A constant μ i sheath model showed the best agreement with the PIC results in the intermediate pressure regime of interest, but a constant λ i model is simpler to implement and showed reasonable agreement over a wide pressure range.

36 MATERIALS SCIENCE↗

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗

3D Printing of Transparent Silicone Elastomers

Here in this paper, inks for transparent elastomers that are formulated by matching the refractive index of silica and polysiloxanes are described. The transparent inks transmit up to ≈90% of 700 nm light through 1 cm and remain transparent when solidified. The inks and solidified materials exhibit a thermochromic effect. This thermochromic effect can be controlled by the refractive index mismatch. Transparency may increase or decrease as temperature increases, depending on the refractive index mismatch of the base polysiloxane and silica. It is found that the rheological properties of the ink depend on the distribution of silica particles, which is dictated by silica functionality, weight content, and processing. Siloxane precursors that vary in chemical functionality are introduced to tailor the mechanical properties of the printed elastomers, which obtain stretchability >65% along with a tensile modulus of 1.9 MPa. After optimizing siloxane chemistry and ink processing, the authors are able to print transparent elastomers. Potential applications are demonstrated for printed structures by printing encapsulation structures for light-emitting diodes, semitransparent dye-filled structures, a microfluidic mixing device, and a multimaterial structure that exhibits temperature-dependent camouflage.

36 MATERIALS SCIENCE↗

Computational predictions of porosities, pore size distributions, and conductivities of aerosol deposited particulate films

The structure and resulting physical properties of aerosol deposited particulate films are governed by the size distribution and morphology of the depositing particles as well as the physics governing particle deposition. While particulate film deposition processes are qualitatively understood, the link between particle characteristics and deposition physics to film properties has not been probed systematically. Here we apply a combination of Langevin dynamics deposition simulations, Monte Carlo pore size distribution calculations, and predictions of thermal conductivities to better establish such process parameter-film property relationships. Furthermore, we account for partial coalescence, polydispersity, and aggregate deposition. We establish that the deposition of polydisperse and aggregated particles leads to broadening of the pore size distribution function and an increase in the mode pore size. Via a non-continuum gas conductivity model, we show that aerosol deposited films can achieve porosities and thermal conductivities similar to conventional aerogels.

36 MATERIALS SCIENCE↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiplicity and rapidity dependence of $\textrm{K}^{*}(\mathrm {{\textbf {892}}})^{0}$ and $\mathrm {\phi ({\textbf {1020}})}$ production in p–Pb collisions at $\sqrt{s_{\textrm{NN}}} =$ 5.02 TeV

The transverse-momentum (p T ) spectra of K*(892) 0 and Φ(1020) measured with the ALICE detector up to p T = 16 GeV/c in the rapidity range -1.2 < y < 0.3, in p–Pb collisions at the center-of-mass energy per nucleon–nucleon collision $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured p T distributions show a dependence on both multiplicity and rapidity at low p T whereas no significant dependence is observed at high p T . A rapidity dependence is observed in the p T -integrated yield (dN/dy), whereas the mean transverse momentum ($\langle$p T $\rangle$) shows a flat behavior as a function of rapidity. The rapidity asymmetry (Y asym ) at low p T (< 5 GeV/c) is more significant for higher multiplicity classes. At high p T , no significant rapidity asymmetry is observed in any of the multiplicity classes. Both K*(892) 0 and Φ(1020) show similar Y asym . The nuclear modification factor (Q CP ) as a function of p T shows a Cronin-like enhancement at intermediate p T , which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high p T (> 5 GeV/c), the Q CP values are greater than unity and no significant rapidity dependence is observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Particle Acceleration and Pitch-angle Evolution in Relativistic Turbulence

Synchrotron radiation detected from relativistic astrophysical objects such as pulsar-wind nebulae and jets from active galactic nuclei depends on the magnetic fields and the distribution functions of energetic electrons in these systems. Relativistic magnetically dominated turbulence has been recognized as an efficient mechanism for structure formation and nonthermal particle acceleration in these environments. Recent numerical simulations of relativistic turbulence have provided insights into the energy distribution functions of accelerated electrons. Much less is currently understood about their pitch angle distributions, which are crucial for accurately interpreting the spectra of synchrotron radiation. We perform a detailed case study of the pitch angle distributions formed during the process of turbulent acceleration for B 0 /δB 0 = 10 and $\tilde{σ}$ 0 ~40, where B0 is the uniform component of the magnetic field, δB0 is the fluctuating component, and $\tilde{σ}$ 0 is the plasma magnetization based on the magnetic fluctuations. We find that even minimal numerical noise can cause substantial pitch angle scattering, but we demonstrate techniques for overcoming the numerical challenges associated with the evolution of very small pitch angles. Our numerical results are consistent with the phenomenological considerations found in C. Vega et al. (2024a, 2025).

79 ASTRONOMY AND ASTROPHYSICS↗

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. In this work, we present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, M ν ≲ 0.6eV, and applicable to redshifts z ≲ 2. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in ~ 25 steps, with a starting redshift of z=99. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between ~ 5-20% for runs with 1024 3 CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Surface Sampling Techniques for the Canister Deposition Field Demonstration

This report describes plans for dust sampling and analysis for the multi-year Canister Deposition Field Demonstration. The demonstration will use three commercial 32PTH2 NUHOMS welded stainless steel storage canisters, which will be stored at an ISFSI site in Advanced Horizontal Storage Modules. One canister will be unheated; the other two will have heaters to achieve canister surface temperatures that match, to the degree possible, spent nuclear fuel (SNF) loaded canisters with heat loads of 10 kW and 40 kW. Surface sampling campaigns will take place on a yearly or bi-yearly basis. The goal of the planned dust sampling and analysis is to determine important environmental parameters that impact the potential occurrence of stress corrosion cracking on SNF dry storage canisters. Specifically, the size, morphology, and composition of the deposited dust and salt particles will be quantified, as well as the soluble salt load per unit area and the rate of deposition, as a function of canister surface temperature, location, time, and orientation. Sampling locations on the canister surface will nominally include 25 locations, corresponding to 5 circumferential locations at each of the 5 longitudinal locations. At each sampling location, a 2x2 sampling grid (containing 4 sample cells) will be painted onto the metal surface. During each sampling campaign, two samples at each sampling location will be collected, in a specific routine to measure both periodic (yearly or bi-yearly) and cumulative deposition rates. For each sample, a wet and a dry sample will be collected. Wet samples will be analyzed to determine the composition of the soluble salt fraction and to estimate salt loading per unit area. Dry samples will be analyzed to assess particle size, morphology, mineralogy, and identity (e.g. for floral/faunal fragments). The data generated by this proposed sampling plan will provide detailed information on dust and salt aerosol deposits on spent nuclear fuel canister surfaces. The anticipated results include information regarding particle compositions, size distributions, and morphologies, in addition to particle deposition rates as a function of canister surface location, orientation, time, and temperature. The information gathered during the Canister Deposition Field Demonstration is critical for ongoing efforts to develop a detailed understanding of the potential for stress corrosion cracking on SNF dry storage canisters

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗