Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle size distribution”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Production of Carbon Nanomaterials and Sorbents from Domestic U.S. Coal (Final Report)

The main goal of this project was to produce high-value carbon nanomaterials and carbon sorbents from domestic coal resources in a cost-effective manner. Four types of domestic coal samples were processed through a combination of deashing, devolatilization, oxidation, reduction, and activation treatments to produce graphene oxide (GO), reduced graphene oxide (RGO), and activated carbon (AC). The precursors and developed materials were extensively characterized by various methods to investigate the impact of the coal feedstock type on the yield and quality of each product. A commercial graphite-based GO sample was included in the experimental work as the baseline material for comparison with the coal-based materials developed in this work. The performed work also included a technoeconomic analysis and cost estimation for a plant processing 20 tons/day coal, a market evaluation for the graphene materials, and a technology gap analysis. A simple process by concentrated nitric acid oxidation is used to oxidize coal precursors for the production of GO. Fine oxidized coal particles that were separated from larger oxidized coal particles had significantly higher oxygen contents and were identified as coal-based GO samples. Coarse oxidized coal particles were used as precursors for production of AC. Based on the Raman spectroscopy results, coal-based GO samples exhibited G and D bands similar to those of graphite-based GO samples. X-ray photoelectron spectroscopy revealed that coal-based GO samples had surface oxygen contents of ~26-35% that were higher than the oxygen contents of graphite-based GO samples. Larger particles of oxidized coal samples were activated under different conditions to produce high surface area functionalized AC. Prepared materials had surface areas exceeding 1,500 m 2 /g and pore volumes more than 1 cm 3 /g with different pore size distributions. Reduced graphene oxide samples were prepared by thermal reduction of both coal-based and graphite-based GO samples at 170-2800 ºC. Coal-based and graphite-based RGO samples exhibited similar carbon contents, Raman spectra, and XRD profiles. Heat treatment above 1500 ºC shifted RGO to synthetic graphite. Among anthracite, bituminous, subbituminous, and lignite coals tested, anthracite was the best precursor to produce carbon nanomaterials exhibiting properties similar to those of graphite-based materials. Anthracite-based carbon nanomaterials also had the highest production yields. Technoeconomic analysis estimated the production cost of GO and RGO for anthracite-based samples at about 2,600 and 4,200 $/ton, respectively, which is about two orders of magnitude lower than the current estimated price of graphite-based materials. Several gaps to further develop the proposed technology were identified and discussed that include process and equipment optimization, process integration, need for additional bench- and pilot-scale experiments, and other items to reduce the scale up risk. Market analysis reports suggested a compound annual growth rate of 40% for graphene and related materials. Short-term applications include composites, inks, and coatings. However, energy storage applications appear to be the dominant potential long-term applications. Different applications of coal-based GO and RGO need to be explored and clear metrics and standards for each application need to be developed.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electrochemical-mechanical coupling failure mechanism of composite cathode in all-solid-state batteries

Composite cathode composed of active particles and solid electrolytes (SEs) can considerably enlarge the particle-SE contact areas and achieve high areal loadings in all-solid-state batteries (ASSBs). However, the challenging interfacial instability and particle damage problems remain unsolved. Herein, we establish a 3D electrochemical-mechanical coupled model to investigate the underlying failure mechanism by considering the governing electrochemical and physics processes. Micro-scale heterogeneous primary particles with random crystallographic orientation and size inside the LiNi 1/3 Co 1/3 Mn 1/3 O 2 (NCM111) secondary particle of the model result in the anisotropic Li diffusion and volume variation within the secondary particle, leading to significant nonuniformity of the Li concentration, and GPa-level stress distributions at primary particle boundaries, and finally causing the particle internal cracks. The particle volume shrinkage under the constraint of stiff Li 7 La 3 Zr 2 O 12 (LLZO) SE triggers the interface debonding (gap>50 nm) with increased interfacial impedance to degrade cell capacity. Higher C-rates result in larger residual stress (~100 MPa)/strain/debonding gap at dis-charging end, more likely to deteriorate the cell performance. Increasing the interfacial strength between the particle and SE can suppress the interface debonding but induces high stress (up to 10 GPa). In conclusion, results reveal the underlying mechanism of the electrochemical-mechanical coupling failure mechanism for composite cathode and provide promising guidance on the further improvement of a more robust composite cathode for ASSBs.

25 ENERGY STORAGE↗
DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Functionalized Silicon Particles for Enhanced Half- and Full-Cell Cycling of Si-Based Li-Ion Batteries

Vinylene carbonate (VC) and polyethylene oxide (PEO) have been investigated as functional agents that mimic the solid electrolyte interphase (SEI) chemistry of silicon (Si). VC and PEO are known to contribute to the stability of Si-based lithium-ion batteries as an electrolyte additive and as a SEI component, respectively. In this work, covalent surface functionalization was achieved via a facile route, which involves ball-milling the Si particles with sacrificial VC and PEO. Thermogravimetric analysis (TGA), X-ray photoelectron spectroscopy (XPS), and magic angle spinning nuclear magnetic resonance (MAS NMR) spectroscopy indicate that the additives are strongly bound to Si. In particular, MAS NMR shows Si–R or Si–O–R groups, which confirm functionalization of the Si after milling in VC or PEO. Particle size analysis by dynamic light scattering reveals that the additives facilitate particle size reduction and that the functionalized particles result in more stable dispersions based on zeta potential measurements. Raman mapping of the electrodes fabricated from the VC and PEO-coated active material with a polyacrylic acid (PAA) binder reveals a more homogenous distribution of Si and the carbon conductive additive compared to the electrodes prepared from the neat Si. Furthermore, the VC-milled Si strikingly exhibited the highest capacity in both half- and full-cell configurations, with more than 200 mAh g –1 measured capacity compared to the neat Si in the half-cell format. This is linked to an improved electrode processing based on the Raman and zeta potential measurements as well as a thinner SEI (with more organic components for the functionalized Si relative to the neat Si) based on XPS analysis of the cycled electrodes. In conclusion, the effect of binder was also investigated by comparing PAA with P84 (polyimide type), where an increased capacity is observed in the latter case.

25 ENERGY STORAGE↗
DOE OSTI2020

3D Nanotomography of calcium silicate hydrates by transmission electron microscopy

Calcium silicate hydrate (C-S-H), is the principal hydration product of Portland cement that mainly contributes to the physical and mechanical properties of concrete. This paper aims to investigate the three-dimensional structure of C-S-H with Ca/Si ratios of 1.0 and 1.6 at the nanoscale using electron tomography. The 3D reconstructions and selected region of interest analysis confirm that the morphology of both C-S-H materials are foil-like structures. Additionally, the difference between the two materials is the density of elongated structures. C-S-H with Ca/Si ratio 1.6 is clearly composed of denser particles compared to the other C-S-H material due to overlapping of the foil-like structure. Pore analysis shows that C-S-H 1.0 and C-S-H 1.6 have porosities 69.2% and 49.8% respectively. Pore size distribution also reveals that C-S-H 1.0 has pore size range between 0-250 nm and C-S-H 1.6 between 0-100 nm. The pore network's size of C-S-H 1.0 is significantly larger than 1.6. This study illustrates the capability of using electron tomography to determine the 3D nanoscale structure of cementitious products and to distinguish between C-S-H 1.0 and 1.6.

42 ENGINEERING↗
DOE OSTI2024

Smart Pixels: In-pixel AI for on-sensor data filtering

We present a smart pixel prototype readout integrated circuit (ROIC) designed in CMOS 28 nm bulk process, with in-pixel implementation of an artificial intelligence (AI) / machine learning (ML) based data filtering algorithm designed as proof-of-principle for a Phase III upgrade at the Large Hadron Collider (LHC) pixel detector. The first version of the ROIC consists of two matrices of 256 smart pixels, each 25$\times$25 µm\textsuperscript{2} in size. Each pixel consists of a charge-sensitive preamplifier with leakage current compensation and three auto-zero comparators for a 2-bit flash-type ADC. The frontend is capable of synchronously digitizing the sensor charge within 25 ns. Measurement results show an equivalent noise charge (ENC) of $\sim$30e\textsuperscript{-} and a total dispersion of $\sim$100e\textsuperscript{-} The second version of the ROIC uses a fully connected two-layer neural network (NN) to process information from a cluster of 256 pixels to determine if the pattern corresponds to highly desirable high-momentum particle tracks for selection and readout. The digital NN is embedded in-between analog signal processing regions of the 256 pixels without increasing the pixel size and is implemented as fully combinatorial digital logic to minimize power consumption and eliminate clock distribution, and is active only in the presence of an input signal. The total power consumption of the neural network is $\sim$ 300 $\mu$W. The NN performs momentum classification based on the generated cluster patterns and even with a modest momentum threshold, it is capable of 54.4\% – 75.4\% total data rejection, opening the possibility of using the pixel information at 40MHz for the trigger. The total power consumption of analog and digital functions per pixel is $\sim$ 6 $\mu$W per pixel, which corresponds to $\sim$ 1 W/cm\textsuperscript{2} staying within the experimental constraints.

Parpillon, Benjamin↗
DOE OSTI2024

Smart Pixels: In-pixel AI for on-sensor data filtering

We present a smart pixel prototype readout integrated circuit (ROIC) designed in CMOS 28 nm bulk process, with in-pixel implementation of an artificial intelligence (AI) / machine learning (ML) based data filtering algorithm designed as proof-of-principle for a Phase III upgrade at the Large Hadron Collider (LHC) pixel detector. The first version of the ROIC consists of two matrices of 256 smart pixels, each 25$\times$25 $\mu$m$^2$ in size. Each pixel consists of a charge-sensitive preamplifier with leakage current compensation and three auto-zero comparators for a 2-bit flash-type ADC. The frontend is capable of synchronously digitizing the sensor charge within 25 ns. Measurement results show an equivalent noise charge (ENC) of $\sim$30e$^-$ and a total dispersion of $\sim$100e$^-$ The second version of the ROIC uses a fully connected two-layer neural network (NN) to process information from a cluster of 256 pixels to determine if the pattern corresponds to highly desirable high-momentum particle tracks for selection and readout. The digital NN is embedded in-between analog signal processing regions of the 256 pixels without increasing the pixel size and is implemented as fully combinatorial digital logic to minimize power consumption and eliminate clock distribution, and is active only in the presence of an input signal. The total power consumption of the neural network is $\sim$ 300 $\mu$W. The NN performs momentum classification based on the generated cluster patterns and even with a modest momentum threshold, it is capable of 54.4% - 75.4% total data rejection, opening the possibility of using the pixel information at 40MHz for the trigger. The total power consumption of analog and digital functions per pixel is $\sim$ 6 $\mu$W per pixel, which corresponds to $\sim$ 1 W/cm$^2$ staying within the experimental constraints.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Kinetic modeling of solitary wave dynamics in a neutralizing ion beam

In this work, we characterize the formation and evolution of electrostatic solitary waves (ESWs) in the space-charge neutralization of ion beams using particle-in-cell simulations. These waves become excited when the electrons emitted from an external filament source initiate a two-stream instability in the beam. We show that such electrostatic waves become excited in both two-dimensional (2D) and three-dimensional (3D) beams with different shapes and sizes. Through a 1D Bernstein–Greene–Kruskal (BGK) analysis of the 2D beam, we find that the non-Maxwellian nature of the beam electrons gives rise to large-sized ESWs that are not predicted by BGK theory since it assumes a Maxwellian electron velocity distribution in the beam. Finally, we show that a 1D BGK theory is inadequate to describe ESWs in 3D beams because of complex electron trajectories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗
DOE OSTI2023

Multiplicity dependence of charged-particle production in pp, p–Pb, Xe–Xe and Pb–Pb collisions at the LHC

Multiplicity (N ch ) distributions and transverse momentum (p T ) spectra of inclusive primary charged particles in the kinematic range of |η| < 0.8 and 0.15 GeV/c < p T < 10 GeV/c are reported for pp, p–Pb, Xe–Xe and Pb–Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{NN}}$ = 2.76 TeV up to 13 TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the p T spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗
DOE OSTI2022

Dose Coefficients for Internal Dose Assessments for Exposure to Radioactive Fallout

This paper presents values as well as the bases for calculating internal dose coefficients suitable for estimating organ doses from the exposure to radioactive fallout that could result from the detonation of a nuclear fission device. The 34 radionuclides discussed are the same as those given in a priority list of radionuclides for fallout dose assessments presented in a companion overview paper. The radionuclides discussed are those that are believed to account for a preponderance of the organ doses that might be received by intake by persons of all ages (including in utero and via breast feeding for infants) following exposure to radioactive fallout. The presented dose coefficients for ingestion account for age and include modifications for variations in solubility with distance as discussed previously in the literature, and those for inhalation similarly account for age, solubility, and particle sizes that would be relevant at various distances of exposure as discussed in a companion paper on ingestion dose methods. The proposed modifications peculiar to radioactive fallout account for systematic changes in solubility and particle sizes with distance from the site of detonation, termed here as the region of “local fallout” and the region “beyond local fallout.” Brief definitions of these regions are provided here with more detailed discussion in a companion paper on estimating deposition of fallout radionuclides. This paper provides the dose coefficients for ingestion and inhalation (for particle sizes of 1 μm, 5 μm, 10 μm, and 20 μm) for the region “local fallout.” These dose coefficients for “local fallout” are specific for particles formed in a nuclear explosion that can be large and have radionuclides, particularly the more refractory ones, distributed throughout the volume where the radionuclide has reduced solubility. The dose coefficients for the region “beyond local fallout” are assumed to be the ones published by the International Commission on Radiological Protection (ICRP) in 1995. Comparisons of the presented dose coefficients are made with values published by the ICRP.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗
DOE OSTI2025

Going Off Grid: A Comparative Study of the Lagrangian and Eulerian Perspectives of New Particle Formation Events

New particle formation and growth (NPF&G) is the process by which ultrafine particles are formed from gas-phase precursors. NPF&G is the dominant source of global aerosol number with important influences on climate. Most observations of NPF&G events are conducted at stationary sites; however, NPF&G observed from stationary sites is influenced by gradual or rapid changes in the air masses passing over the site, complicating NPF&G analysis. In this work, we use observations and a 3D aerosol model to compare aerosol size distributions at a stationary site (Southern Great Plains [SGP] observatory, Oklahoma, USA) and along Lagrangian trajectories crossing the site. The model simulates the NPF&G events reasonably well at SGP. Using the model to compare the Lagrangian and stationary perspectives, we can explain previously unanalyzable days with some evidence of NPF&G as either non-event or analyzable NPF&G days. We find most of the unanalyzable NPF&G days are due to isolated and inhomogeneous NPF&G occurring upwind of the stationary site, often in the outflow of urban regions. Finally, we compare formation rates of 3 nm particles, growth rates, and the survival probability of 3 nm particles growing to 25 nm between the stationary and Lagrangian perspectives. Because of the much larger number of analyzable days along the Lagrangian trajectories, this perspective potentially provides more robust statistics and better characterization of NPF&G event extremes. Our method for extracting chemical/physical properties along Lagrangian trajectories from 3D models can be applied to a wide range of science questions.

O’Donnell, Samuel E. [Colorado State Univ., Fort C↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗