Optimal iteration and its application to some problems in aerosol science and particle dynamics
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, the effective diffusion coefficient of the fission product Ag through the grain boundary (GB) network is calculated using a combination of atomistic and phase-field methods. Atomistic calculations of Ag diffusivity in SiC bulk and GBs are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Since GBs serve as pathways for Ag diffusion, Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.
Abstract Synthetic self‐propelled nano and microparticles have a growing appeal for targeted drug delivery, collective functionality, and manipulation at the nanoscale. However, it is challenging to control their positions and orientations under confinement, e.g., in microchannels, nozzles, and microcapillaries. This study reports on the synergistic effect of acoustic and flow‐induced focusing in microfluidic nozzles. In a microchannel with a nozzle, the balance between the acoustophoretic forces and the fluid drag due to streaming flows generated by the acoustic field controls the microparticle's dynamics. This study manipulates the positions and orientations of dispersed particles and dense clusters inside the channel at a fixed frequency by tuning the acoustic intensity. The main findings are: first, this study successfully manipulates the positions and orientations of individual particles and dense clusters inside the channel at a fixed frequency by tuning the acoustic intensity. Second, when an external flow is applied, the acoustic field separates and selectively extrudes shape‐anisotropic passive particles and self‐propelled active nanorods. Finally, the observed phenomena are explained by multiphysics finite‐element modeling. The results shed light on the control and extrusion of active particles in confined geometries and enable applications for acoustic cargo (e.g., drug) delivery, particle injection, and additive manufacturing via printed self‐propelled active particles.
Several experimental disciplines require simultaneous neutron flux and temperature monitoring. This paper provides the groundwork on a diamond sensor concept that could be used in harsh environments where the limits of current technologies are exceeded. Here, we present initial results on the functionality of a multimodal diamond sensor, which can simultaneously sense the presence of alpha particles at a few counts per second at temperatures up to 494 K. Such a sensor may also have applications where both charged particles or neutron detection and temperature measurement are simultaneously desired.
This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.
Here, we perform a comprehensive study of cosmological constraints on nonstandard neutrino self-interactions using cosmic microwave background and baryon acoustic oscillation data. We consider different scenarios for neutrino self-interactions distinguished by the fraction of neutrino states allowed to participate in self-interactions and how the relativistic energy density, N eff , is allowed to vary. Specifically, we study cases in which all neutrino states self-interact and N eff varies; two species free-stream, which we show alleviates tension with laboratory constraints, while the energy in the additional interacting states varies; and a variable fraction of neutrinos self-interact with either the total N eff fixed to the Standard Model value or allowed to vary. In no case do we find compelling evidence for new neutrino interactions or nonstandard values of N eff . In several cases, we find additional modes with neutrino decoupling occurring at lower redshifts z dec ~ 10 3–4 . We do a careful analysis to examine whether new neutrino self-interactions solve or alleviate the so-called H 0 tension and find that, when all Planck 2018 CMB temperature and polarization data are included, none of these examples eases the tension more than allowing a variable N eff comprised of free-streaming particles. Although we focus on neutrino interactions, these constraints are applicable to any light relic particle.
As a second generation of oxy-fuel combustion technologies, pressurized oxy-combustion has a potential to improve the process efficiency and economy. With the modern decarbonization needs, biomass combustion is considered to be nearly CO2 neutral, and thus co-firing biomass with coal in an oxy-combustion process can achieve a negative CO2 balance. The size of the biomass fuel particles is a key parameter in an entrained flow combustor because large biomass particles, exceeding some threshold size, may only partially be burnt due to the restrictions on the heat transport, thereby diminishing the combustion efficiency. While such a threshold particle size has been scrutinized at an atmospheric pressure, the present study extends the analysis to elevated pressures. Specifically, two scenarios were set up, for experimental-scale and full-scale pressurized oxy-combustors, with the threshold sizes of biomass particles appearing 0.5 mm and 5 mm, respectively, without temperature gradient in biomass particles. With temperature gradient, the parametric study included the particles in the range from 0.15 mm to 5 mm. It is shown that the thermal-thin model is applicable to spherical biomass particles not exceeding 0.15 mm. The impact of the heating rate on biomass particle devolatilization was also investigated and the biomass devolatilization could happen in less than 1 second; the same amount of heat energy could induce different volatile release rates.
While many aspects of single-particle electron cryo-microscopy (cryo-EM) of biological macromolecules have reached a sophisticated level of development, this is not yet the case when it comes to preparing thin samples on specimen grids. As a result, there currently is considerable interest in achieving better control of both the sample thickness and the amount of area that is useful, but this is only one aspect in which improvement is needed. This Perspective addresses the further need to prevent the macromolecular particles from making contact with the air-water interface, something that can result in preferential orientation and even structural disruption of macromolecular particles. This unwanted contact can occur either as the result of free diffusion of particles during the interval between application, thinning and vitrification of the remaining buffer, or—when particles have been immobilized—by the film of buffer becoming too thin prior to vitrification. An opportunity now exists to apply theoretical and practical insights from the fields of thin-film physical chemistry and interfacial science, in an effort to bring cryo-EM sample preparation to a level of sophistication that is comparable to that of current data collection and analysis.
GRAPPA is a Geant4 application that simulates the interaction of an incoming particle beam (typically electron or photons) with an High-Z material. Particles are generated at a certain energy and interact with a solid target. The product of this interaction is detected and saved on file. The application implements a basic target that is just a simple box, and it accepts complex target descriptions via GDML files. The user has the freedom to set detectors to capture the desired particles. This application enables studying processes such as antiparticle production, muon generation, and atomic decay.
In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.
Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.
Increasing penetration of variable renewable energy resources requires the deployment of energy storage at a range of durations. Long-duration energy storage (LDES) technologies will fulfill the need to firm variable renewable energy resource output year round; lithium-ion batteries are uneconomical at these durations. Thermal energy storage (TES) is one promising technology for LDES applications because of its siting flexibility and ease of scaling. Particle-based TES systems use low-cost solid particles that have higher temperature limits than the molten salts used in traditional concentrated solar power systems. A key component in particle-based TES systems is the containment silo for the high-temperature (>1100 degrees C) particles. This study combined experimental testing and computational modeling methods to design and characterize the performance of a particle containment silo for LDES applications. A laboratory-scale silo prototype was built and validated the congruent transient finite element analysis (FEA) model. The performance of a commercial-scale silo was then characterized using the validated model. The commercial-scale model predicted a storage efficiency above 95% after 5 days of storage with a design storage temperature of 1200 degrees C. Insulation material and concrete temperature limits were considered as well. The validation of the methodology means the FEA model can simulate a range of scenarios for future applications. This work supports the development of a promising LDES technology with implications for grid-scale electrical energy storage, but also for thermal energy storage for industrial process heating applications.
A semiconductor detection material, Cadmium Zinc Telluride, was investigated for potential use as a particle detector for space flight applications. Our set of detectors was determined to be flawed for particle detection methods, likely due to charge transport mechanisms within the detector material, or due to reduced voltage required to avoid dielectric breakdown. It was observed that the method of crystal growth may play a role in the performance, but performance may also vary between manufacturers using the same growth method. It is recommended to perform further study at a larger proton beam facility.
The detonation of high explosives produces a wide variety of particulate matter (PM) with distinct properties, not all of which are traditionally studied for chemical composition, formation processes, and forensic applications. We report particle-resolved measurements of Composition B detonation soot using soot particle aerosol mass spectrometry (SP-AMS), identifying carbonaceous species and metals not previously characterized on a single-particle basis. Results are combined with scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) to enhance source-dependent signatures. Black carbon, including graphitic carbon and detonation nanodiamonds, contributed 50.5-71.4% of PM 2.5 mass, while 22.5-43.4% was non-refractory organic carbon, a previously overlooked component that exhibited a complex and varying composition. Detonations were performed with and without PMMA confinement and under steady and overdriven conditions. Distinct particles enriched with polycyclic aromatic hydrocarbons (PAHs) were observed in experiments utilizing confinement, with quantities dependent on manufacturing method. SEM-EDS validated SP-AMS findings of metals internally mixed with carbonaceous species and extended the particle size range to 100 µm. This work makes detecting non-refractory organics using single-particle techniques more feasible for detonation forensics and understanding high-energy soot formation. While this analysis uses offline sample collection, SP-AMS could be deployed for in-situ measurements of detonation PM 2.5 transported in the atmosphere.
In the presence of large size disparities, single-grid neighbor search algorithms lead to inflated neighbor lists that significantly degrade the performance of Lagrangian particle solvers. If Eulerian–Lagrangian (EL) frameworks are to remain performant when simulating realistic systems, improved neighbor detection approaches must be adopted. To this end, we consider the application of a multigrid neighbor search (MGNS) algorithm in the mfix-exa software package, an exascale EL solver built upon the AMReX library. Here, details regarding the implementation and verification of MGNS are provided along with speedup curves for a bidisperse mixing layer. MGNS is shown to yield up to 15$\times$ speedup on CPU and 6$\times$ speedup on GPU for the problems considered here. The mfix-exa software is then validated for a variety of polydisperse flows. Finally, a brief discussion is given for how dynamic MGNS may be completed, with application to spatially varying particle size distributions.
The creation of optically powered self-assembling nano-to-meso-scale machines that do work is a long-standing goal in photonics. We demonstrate an optical matter (OM) machine that converts the spin angular momentum (SAM) of light into orbital angular momentum (OAM) to do mechanical work. The specific OM machine we study is based on a sixfold symmetric hexagonally ordered nanoparticle array that operates as an OM “gear” that is assembled and made to rotate in a circularly polarized Gaussian beam. The rotational symmetry of the OM gear leads to a selection rule for the allowed scattering modes based on their angular momentum. Electrodynamics calculations show that the collective scattering modes with the largest angular momentum scatter strongly in the transverse direction. Simulations and experiments show that the angular momentum that accompanies the scattered light causes a “negative torque” response on the OM gear and drives a “probe” particle placed outside the OM gear around the gear in an asymmetric force field analogously to Brownian ratchets. The gear–probe OM machine concept can be expanded to applications in nanofluidics and particle sorting.
Herein we present the cost-effective production of superconducting radio frequency (SRF) cavities made of medium grain (MG) niobium (Nb) discs directly sliced from forged and annealed billet. This production method provides clean surface conditions and reliable mechanical characteristics with sub-millimeter average grain size resulting in stable SRF cavity production. We propose to apply this material to particle accelerator in the science and industrial applications. The science applications require high field gradients(≥ 30 MV/m) particularly in pulsed mode. The industrial applications require high $Q_0$ values with moderate gradients (~ 20 MV/m) in CW mode operation. This paper describes the MG Nb disc production recently demonstrated and discusses future prospects for application in advanced particle accelerators in the science and industrial applications.
Accurately measuring the translations of objects between images is essential in many fields, including biology, medicine, chemistry, and physics. One important application is tracking one or more particles by measuring their apparent displacements in a series of images. Popular methods, such as the center of mass, often require idealized scenarios to reach the shot noise limit of particle tracking and, therefore, are not generally applicable to multiple image types. More general methods, such as maximum likelihood estimation, reliably approach the shot noise limit, but are too computationally intense for use in real-time applications. These limitations are significant, as real-time, shot-noise-limited particle tracking is of paramount importance for feedback control systems. To fill this gap, we introduce a new cross-correlation-based algorithm that approaches shot-noise-limited displacement detection and a graphics processing unit-based implementation for real-time image analysis of a single particle.