Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Using Artificial Soiling to Rank Anti-Soiling Coatings for Arid Climates

A simple method to evaluate (e.g., screen and rank) the performance of coatings that are fully or partially anti-soiling (AS) is lacking within the PV industry. Artificial soiling may be used as a rapid and economical assessment approach, offering an efficient alternative to time-intensive, site-specific field testing. In this study, we present an artificial soiling method to replicate the anti-soiling performance rankings of two groups of coated glass samples supplied by two manufacturers (group C with CA, CP, and CS coatings; group A with AC coating). These are compared to field aging in two climates: semi-arid Lemoore, California over 4 months, and the hot desert in Mesa, Arizona over 7 months. Both field and indoor performance ranking utilize the transmittance ratio, defined as the optical transmittance between a coated sample and an uncoated reference in each group, as an evaluation metric to assess the effectiveness of the artificial soiling approach in replicating field soiling. It is critical to mimic the dominant field meteorological conditions associated with soiling-prone days and vulnerable times of day during the soiling season. Our results reveal that the use of artificial soiling for field performance ranking among coatings strongly depends on the soil type (composition and particle distribution), dust surface density, and prevalent environmental factors (wetting saturation by humidity, dew condensation extent, and prior weathering history of the coating). The similar rank order relative to the field suggests that the artificial soiling approach presented may be used for down-selecting anti-soiling coatings prior to prolonged field validation.

14 SOLAR ENERGY↗
DOE OSTI2022

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Microextraction-Single Particle-Inductively Coupled Plasma-Mass Spectrometry for the Direct Analysis of Nanoparticles on Surfaces

A novel employment of single particle-inductively coupled plasma-mass spectrometry (SP-ICP-MS) was developed, where a microextraction (ME) probe is used to sample nanoparticles from a surface and analyze them in a single analytical step. The effects of several parameters on the performance of ME-SP-ICP-MS were investigated, including the flow rate, choice of carrier solution, particle size, and the design of the microextraction probe head itself. The optimized ME-SP-ICP-MS technique was used to compare the extraction efficiency (EE, defined as the ratio of particles measured to particles deposited on the surface) of the commercial probe head to a newly designed SP polyether ether ketone (PEEK) probe head. The SP PEEK probe head was found to have increased EE compared to the commercial probe head (8.5 ± 3% vs 3.9 ± 3%, respectively). Increasing the carrier solution flow rate was found to decrease the total analysis time at the cost of decreasing EE. Extraction efficiencies for ME-SP-ICP-MS were typically 4–10%, which is similar to transport efficiencies (1–10%) for conventional SP-ICP-MS. Lastly, ME-SP-ICP-MS was employed for the analysis of nano- and microparticles. The sizes of gold nanoparticles, 30 ± 3 and 51 ± 1.9 nm (certified sizes), and iron-based microparticles, 1000 ± 50 nm (certified size), were accurately determined to be 32.2 ± 2.5, 50.8 ± 3.4, and 1030 ± 57 nm, respectively, by ME-SP-ICP-MS. Further, this work demonstrates the potential of ME-SP-ICP-MS for the direct analysis of particles on common collection surfaces (GSR tabs, carbon planchettes, etc.) while retaining spatial information on particle distribution across the surface.

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

Valorization of consolidated bioprocessing residues for bioplastics

This study demonstrates an organic solvent-free processing strategy to valorize consolidated bioprocessing (CBP) residues, from switchgrass and poplar biomass, into functional poly(butylene succinate) (PBS)-based biocomposites using high-shear homogenization (HSH). HSH transformed the switchgrass and poplar CBP residues (CBP-R) into fine, uniformly distributed particles and microfibers. The composites of PBS with homogenized switchgrass residues (H-CBP-R-SG) or homogenized poplar residues (H-CBP-R-P) at a 70/30 weight ratio exhibited improved processability and mechanical integrity, with the Young's modulus for the PBS/H-CBP-R-SG and PBS/H-CBP-R-P nearly doubling to 0.66 ± 0.07 GPa and 0.65 ± 0.04 GPa, respectively, compared to neat PBS (0.36 ± 0.02 GPa). Dynamic Mechanical Analysis (DMA) reveals a significant suppression of the tan δ peak magnitude, indicating that HSH-mediated physical activation facilitates stress transfer in composites typical of covalent chemical grafting systems. While the transition to a stiffness-dominated profile reduces ductility, the resulting composites exhibit the dimensional stability and resistance to thermal warping required for high-fidelity FDM 3D printing and injection molding. Beyond material performance, comprehensive techno-economic analysis (TEA) and life cycle assessment (LCA) confirmed that diverting CBP residues into composite products can improve the economic viability of the biorefinery without substantially increasing biorefinery global warming potential (GWP). At a 30 wt% blend ratio, incorporating residuals into PBS yielded a minimum selling price for the composite of $\$4.07$ per kg compared to the conventional bioplastic price of $\$5.00$ per kg. This approach aligns with circular bioeconomy principles by converting waste streams into value-added products. Furthermore, this innovative strategy addresses key challenges in bioplastic development, including cost, compatibility, and performance, while simultaneously advancing waste minimization strategies for sustainable manufacturing systems.

09 BIOMASS FUELS↗
DOE OSTI2020

Modeling q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of a drive cylinder for low collisional experiments on magnetic reconnection

For most laboratory plasma experiments, Coulomb collisions between the particle species are sufficiently frequent that the particle distribution functions are relaxed to a near-Maxwellian form. This hampers the applicability of such experiments to phenomena observed in tenuous and near-collisionless space plasma. The Terrestrial Reconnection EXperiment (TREX) at the Wisconsin Plasma Physics Laboratory aims to study collisionless reconnection for parameters relevant to the Earth’s magnetosphere. To reduce the role of collisional effects, a reconnection Drive Cylinder has been developed, which increases both the effective system size of the TREX configuration and the rate at which reconnection can be driven. These two effects now permit TREX to reach a kinetic reconnection regime where collisional effects are minimized. The Drive Cylinder is comprised of 12 single loop drive-coils connected in parallel to a 10 kV capacitor bank. Insulated sheets of aluminum are applied to smooth the magnetic fields and enhance the drive efficiency. Following is a description of the technical details and performance of the Drive Cylinder.

Instruments & Instrumentation↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2024

Lower hybrid drift instability in nonthermal plasmas

Lower hybrid drift instability (LHDI) is driven by the cross-field current and operates in the vicinity of the lower-hybrid frequency, between the ion- and electron-gyro frequencies, and with wavelengths between the electron and ion thermal gyro radii. The free energy source that drives this instability resides in the density gradient associated with an inhomogeneous plasma. The existing literature on LHDI assumes that the charged particle distribution function is given by a Maxwellian form, but the space plasma is pervasively observed to feature nonthermal characteristics. This paper extends the theory of LHDI to nonthermal plasmas. The generalized theory of LHDI is, thus, applicable to various space plasma environments characterized by nonthermal plasma velocity distribution functions.

Arya, Neetasha (ORCID:0009000236075605)↗
DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Imaging nuclei by smashing them at high energies: how are their shapes revealed after destruction?

High-energy nuclear collisions have recently emerged as a promising ‘imaging-by-smashing’ approach to reveal the intrinsic shapes of atomic nuclei. Here, I outline a conceptual framework for this technique, explaining how nuclear shapes are encoded during quark–gluon plasma (QGP) formation and evolution, and how they can be decoded from final-state particle distributions. Furthermore, I highlight the method’s potential to advance our understanding of both nuclear structure and QGP physics.

imaging by smashing↗
DOE OSTI2020

Hydrodynamic generators in relativistic kinetic theory

Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Benchmark eA generator for leptoproduction in high-energy lepton-nucleus collisions

The upcoming electron-ion collider (EIC) will address several outstanding puzzles in modern nuclear physics. Topics such as the partonic structure of nucleons and nuclei, the origin of their mass and spin, among others, can be understood via the study of high-energy electron-proton (ep) and electron-nucleus (eA) collisions. Achieving the scientific goals of the EIC will require a novel electron-hadron collider and detectors capable to perform high-precision measurements but also dedicated tools to model and interpret the data. To aid in the latter, we present a general-purpose e A Monte Carlo generator—BeAGLE. In this paper, we provide a general description of the models integrated into BeAGLE, applications of BeAGLE in eA physics, implications for detector requirements at the EIC, and the tuning of the parameters in BeAGLE based on available experimental data. Specifically, we focus on a selection of model and data comparisons in particle production in both ep and eA collisions, where baseline particle distributions provide essential information to characterize the event. In addition, we investigate the collision geometry determination in eA collisions, which could be used as an experimental tool for varying the nuclear density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quasiclassical evaluation of gluon saturation induced helicity effects

At sub-eikonal order in the high-energy limit, helicity dependences are generally contained in the transverse components of the gluon field. In the gluon saturation regime, novel helicity effects arise from the nonlinear interaction between the eikonal-order longitudinal gluon field and the sub-eikonal-order transverse gluon field. We derive this saturation induced helicity-dependent field both by solving the classical Yang-Mills equations and through direct diagrammatic calculations. Our analysis shows that the saturation induced helicity effect is intrinsically a two-particle (or multiparticle) correlation effect, rather than a single-particle distribution effect. Furthermore, we evaluate this effect in the context of double-spin asymmetry for incoherent diffractive dijet production in longitudinally polarized electron-nucleus collisions, using a simplified helicity-extended McLerran-Venugopalan model. We find that the saturation induced helicity effect further suppresses the back-to-back peak in the dijet azimuthal angle correlation, with a magnitude comparable to that of the direct helicity effect. This gluon saturation induced helicity effect may offer a novel avenue to probe gluon saturation in polarized collisions. Published by the American Physical Society 2025

Li, Ming (ORCID:0000000195509691)↗
DOE OSTI2021

Adaptive Spatially Aware I/O for Multiresolution Particle Data Layouts

Large-scale simulations on nonuniform particle distributions that evolve over time are widely used in cosmology, molecular dynamics, and engineering. Such data are often saved in an unstructured format that neither preserves spatial locality nor provides metadata for accelerating spatial or attribute subset queries, leading to poor performance of visualization tasks. Furthermore, the parallel I/O strategy used typically writes a file per process or a single shared file, neither of which is portable or scalable across different HPC systems. We present a portable technique for scalable, spatially aware adaptive aggregation that preserves spatial locality in the output. We evaluate our approach on two supercomputers, Stampede2 and Summit, and demonstrate that it outperforms prior approaches at scale, achieving up to 2.5× faster writes and reads for nonuniform distributions. Furthermore, the layout written by our method is directly suitable for visual analytics, supporting low-latency reads and attribute-based filtering with little overhead.

Usher, Will↗
DOE OSTI2024

Jet substructure

Jet substructure observables hold the keys to identifying the inner working of the quark–gluon plasma through its imprints in jet modification patterns. Because of the multi-scale nature of jet evolution, the strategy has focused on developing jet reclustering-based and jet shape-based observables to specifically probe different stages of jet evolution, which provide multi-dimensional quantification of jet particle distribution phase spaces. The developments of jet substructure in the past few years allowed us to explore the correlations among jet observables. The richer information beyond single-observable modifications could give provides a crucial step beyond the energy loss paradigm. This review surveys the theory and phenomenology of modern jet substructure observables, which help constraining the energy scales of medium interactions and identifying heavy flavor quarks within jets. Also, the peculiar soft activities point to more refined studies of medium responses to jets.

Physics↗
DOE OSTI2024

Investigation of Proton Beam-Driven Fusion Reactions Generated by an Ultra-Short Petawatt-Scale Laser Pulse

We present results from a pitcher-catcher experiment utilizing a proton beam generated with nanostructured targets at a petawatt-class, short-pulse laser facility to induce proton-boron fusion reactions in a secondary target. A 45-fs laser pulse with either 400 nm wavelength and 7 J energy, or 800 nm and 14 J, and an intensity of up to 5 × 10 21 W/cm 2 was used to irradiate either thin foil targets or near-solid density, nanostructured targets made of boron nitride (BN) nanotubes. In particular, for 800 nm wavelength irradiation, a BN nanotube target created a proton beam with about five times higher maximum energy and about ten times more protons than a foil target. This proton beam was used to irradiate a thick plate made of boron nitride placed in close proximity to trigger 11 B (p, α) 2α fusion reactions. A suite of diagnostics consisting of Thomson parabola ion spectrometers, postshot nuclear activation measurements, neutron time-of-flight detectors, and differentially filtered solid-state nuclear track detectors were used to measure both the primary proton spectrum and the fusion products. From the primary proton spectrum, we calculated (p, n) and (α,n) reactions in the catcher and compare with our measurements. The nuclear activation results agree quantitatively and neutron signals agree qualitatively with the calculations, giving confidence that primary particle distributions can be obtained from such measurements. These results provide new insights for measuring the ion distributions inside of proton-boron fusion targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulating Boundary Region Plasma Instabilities

The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗