SEARCH · Engineering Papers

Results for “Particles”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗

DOE OSTI2020

An Assessment of the Drag Models in the Case of a Shock Interacting With a Fixed Bed of Point Particles

Abstract In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.

Engineering↗

DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗

DOE OSTI2015

Excitation of XPS spectra from nanoscaled particles by local generation of x-rays

In preliminary work, the authors have shown that use of an aluminum substrate to support a distribution of copper particles enables their characteristic photoelectrons to be observed within the Auger electron spectrum generated by an incident electron beam. This observation raises the possibility of the use of chemical shifts and the corresponding Auger parameter to identify the chemical states present on the surface of individual submicrometer particles within a mixture. In this context, the technique has an advantage in that, unlike conventional Auger electron spectroscopy, the electron beam does not dwell on the particle but on the substrate adjacent to it. Given the importance, for both medical and toxicological reasons, of the surface composition of such particles, the authors have continued to explore the potential of this development. In this contribution, the authors show that proximal excitation of x-rays is equally successful with magnesium substrates. In some regions of the x-ray photoelectron spectrum, the much larger Auger peaks generated by the electron beam can cause inconvenient clustering of Auger and photoelectron peaks. As in conventional x-ray photoelectron spectroscopy, the ability to switch between Al and Mg sources is useful in such situations. In this context, the authors have extended the studies to iron particles where the authors show that use of Al or Mg substrates, as necessary, can make a contribution to clear identification of individual components in the Fe 2p peaks. For this development in electron spectroscopy to achieve its full potential, it is necessary to optimize the beam conditions used to generate the local x-ray to give good selectivity of a given particle. Measurements made in support of this will be given. Of greater concern is a possible problem of local heating associated with x-ray generation. The authors continue to explore this problem and report some progress in minimizing heating of the particle while maintaining the particle selectivity that is central to this exciting development.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗

DOE OSTI2019

Building Advanced Materials via Particle Aggregation and Molecular Self-Assembly

Hierarchical and other advanced materials have attracted increasing attention due to their unique physical and chemical properties, which strongly depend on morphology and size.1-2 These materials have been applied in important technological fields such as energy, catalysis, optical devices, water purification, pollutant removal, CO2 sequestration, and biomedicine.3-7 Particle-based crystallization and self-assembly of molecules are important pathways to synthesize advanced materials of complex structures.8-11 Unlike monomer-by-monomer addition or Ostwald ripening, particle-based crystallization occurs via particle-by-particle addition, to form larger crystals or clusters.8, 12 To date, numerous advanced materials have been synthesized in the lab using particle-based crystallization. Examples include metals such as Pt, Pd, Au, Ag, and Cu;13 alloys such as Pt-Ni, Pt-Cu, Pt-Fe, and Au-Ag;14 metal oxides such as ZnO, TiO2, CuO, and a-Fe2O3, Fe3O4;15-19 and metal sulfides such as PbS, PbSe, ZnS, and CdS.20-21 In addition, particle aggregation-based crystallization has been observed in nature, such as in various geological and biological minerals including calcite, collagen, and others.22-24 Different from the particle-based crystallization, self-assembly of molecules has also been used to build advanced materials such as molecular clusters and nanoparticles. For instance, advanced luminescent materials have been prepared by aggregation-induced emission (AIE) of intrinsically non-emissive molecules.25-26 One of the challenges facing this fast-growing field of advanced materials is to develop a fundamental understanding of the interactions between particles or molecules in a growth medium and the resulting response dynamics.

Zhang, Xin↗

DOE OSTI2022

EPMA-Based Mass Balance Method for Quantitative Fission Product Distribution Comparison between TRISO Particles

Two irradiated AGR-2 TRISO particles were chosen to demonstrate a recently developed mass balance technique in which EPMA-generated concentration data was used to determine fission product mass on a layer-by-layer basis in TRISO particles. EPMA-calculated fission product masses for most fission products in the two particles were within +/- 20% of their ORIGEN-modelled masses. Results show that Sr, Ba, and Eu accumulate preferentially in the carbon-rich kernel periphery on the particles’ side that lacks a gap between the buffer and IPyC. In addition, the more mobile elements--Cs, Sr, and Pd, accumulate in greater quantity in the outer layers of particle AGR2-223-RS34 compared to particle AGR2-223-RS06, which has relatively more of those elements’ mass located in the kernel and kernel periphery, suggesting enhanced fission product transport in particle AGR2-223-RS34. This model can be used better understand and test fission product transport in TRISO particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2012

Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (Final Report)

At Colorado University-Boulder the primary task is to extend our gyrokinetic Particle-in-Cell simulation of tokamak micro-turbulence and transport to the area of energetic particle physics. We have implemented a gyrokinetic ion/massless fluid electron hybrid model in the global δƒ-PIC code GEM, and benchmarked the code with analytic results on the thermal ion radiative damping rate of Toroidal Alfvén Eigenmodes (TAE) and with mode frequency and spatial structure from eigenmode analysis. We also performed nonlinear simulations of both a single-n mode (n is the toroidal mode number) and multiple-n modes, and in the case of single-n, benchmarked the code on the saturation amplitude vs. particle collision rate with analytical theory. Most simulations use the f method for both ions species, but we have explored the full-ƒ method for energetic particles in cases where the burst amplitude of the excited instabilities is large as to cause significant re-distribution or loss of the energetic particles. We used the hybrid model to study the stability of high-n TAEs in ITER. Our simulations show that the most unstable modes in ITER lie in the rage of 10 < n < 20. Thermal ion pressure effect and alpha particles non-perturbative effect are important in determining the mode radial location and stability threshold. The thermal ion Landau damping rate and radiative damping rate from the simulations are compared with analytical estimates. The thermal ion Landau damping is the dominant damping mechanism. Plasma elongation has a strong stabilizing effect on the alpha driven TAEs. The central alpha particle pressure threshold for the most unstable n=15 mode is about β α (0) = 0.7% for the fully shaped ITER equilibrium. We also carried nonlinear simulations of the most unstable n = 15 mode and found that the saturation amplitude for the nominal ITER discharge is too low to cause large redistribution or loss of alpha particles. To include kinetic electron effects in the hybrid model we have studied a kinetic electron closure scheme for the fluid electron model. The most important element of the closure scheme is a complete Ohm’s law for the parallel electric field E || , derived by combining the quasi-neutrality condition, the Ampere’s equation and the v || moment of the gyrokinetic equations. A discretization method for the closure scheme is studied in detail for a three-dimensional shear-less slab plasma. It is found that for long-wavelength shear Alfvén waves the kinetic closure scheme is both more accurate and robust than the previous GEM algorithm using the split-scheme, whereas for the ion-gradient-driven instability the previous algorithm is more efficient. This kinetic electron closure scheme will be implemented in GEM in the future. We have studied the beam driven Reverse Shear Alfven Eigenmodes (RSAE) observed in DIII-D discharge 142111. For this purpose a new scheme for obtaining the electric potential is implemented, i.e., by solving the gyrokinetic moment (GKM) equation, which is essentially the equation for $∂\phi/∂t$ used in GEM’s split-weight scheme, and then integrating in time. Due to charge-neutrality the ExB motions of the equilibrium densities of all species cancel each other and do not cause charge separation if there is no finite Larmor radius effect. The advantage of solving the GKM equation is that this lowest-order cancellation can be made explicit. The GKM approach is found to be more accurate and robust. GEM simulations have reproduced many features of RSAE seen in the experiment, such as frequency chirping and the chirping range. It has been reported by other simulation codes that the shearing direction of the mode structure in the poloidal plane disagrees with observation. We found that the mode structure, including the shearing in the poloidal plane, is in general sensitive to the beam distribution. Using the same beam density profile as in other codes but with a slowing-down distribution in velocity, GEM simulations reproduce the mode shearing direction seen in the experiment. We have carried out extensive nonlinear simulations of RSAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Effects of volatile coatings on the morphology and optical detection of combustion-generated black carbon particles.

We have measured time-resolved laser-induced incandescence (LII) from combustion-generated mature soot extracted from a burner and (1) coated with oleic acid or (2) coated with oleic acid and then thermally denuded using a thermodenuder. The soot samples were size selected using a differential mobility analyser and characterized with a scanning mobility particle sizer, centrifugal particle mass analyser, and transmission electron microscope. The results demonstrate a strong influence of coatings particle morphology and on the magnitude and temporal evolution of the LII signal. For coated particles higher laser fluences are required to reach LII signal levels comparable to those of uncoated particles. This effect is predominantly attributable to the additional energy needed to vaporize the coating while heating the particle. LII signals are higher and signal decay rates are significantly slower for thermally denuded particles relative to coated or uncoated particles, particularly at low and intermediate laser fluences.

Bambha, Ray.↗

DOE OSTI2022

Design and Technoeconomic Analysis of High-Temperature Particle Conveyance Components for a 100 MW e Concentrating Solar Power Plant

Levelized costs of electricity (LCOE) approaching the U.S. Department of Energy Solar Energy Technologies Office 2030 goal of 0.05 $/kWh may be achievable using Brayton power cycles that use supercritical CO2 as the working fluid and flowing solid particles with temperatures >700° C as the heat transfer media. The handling and conveyance of bulk solid particles at these temperatures in an insulated environment is a critical technical challenge that must be solved for this approach to be used. A design study was conducted at the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories in Albuquerque, NM, with the objective of identifying the technical readiness level, performance limits, capital and O&M costs, and expected thermal losses of particle handling and conveyance components in a particle-based CSP plant. Key findings indicated that chutes can be a low-cost option for particle handling but uncertainties in tower costs make it difficult to know whether they can be cost effective in areas above the receiver if tower heights must then be increased. Skips and high temperature particle conveyance technology are available for moving particles up to 640° C. This limits the use of mechanical conveyance above the heat exchanger and suggests vertical integration of the hot storage bin and heat exchanger to facilitate direct gravity fed handling of particles.

14 SOLAR ENERGY↗

DOE OSTI2020

Reducing Exposures to Airborne Particles Through Improved Filtration: A High-Level Modeling Analysis

Indoor airborne particulates are well-known health hazards and filtration is one common method of reducing exposures. Based on our previously developed Regional Shelter Analysis methodology and parameters that characterize the existing US building stock, we perform a high-level assessment of the potential benefits of upgrading existing filters in furnace and in heating, ventilation, and air conditioning systems using off-the-shelf filters. We use three metrics to assess the improvement: Building Transmission Factor (a measure of protection against outdoor airborne particles), Indoor Normalized Time and Space Integrated Air Concentration (a measure of indoor exposure to indoor-origin airborne particles), and Building Exit Fraction (fraction of indoor airborne particles that are released to the outdoor atmosphere). We also discuss the potential reduction in regional exposures due to particles exiting the building and exposing downwind building occupants. Our modeling indicates that while buildings provide their occupants passive protection against airborne particulate hazards, including but not limited to PM 2.5 , PM 10 , and wildfire smoke; improving particle filtration efficiency may further improve this protection. The degree of improvement varies with particle size and building type. Of the building types studied, apartments are predicted to benefit most, with greater than a factor of 2 improvement (≥50% reduction in exposures) for 1 µm particle exposures when using MERV 7 to 12 rated filters. Non-residential buildings were notably less responsive to improved filtration but had the highest Building Exit Fractions with 30% to 40% of indoor airborne particles released to the outdoor atmosphere (apartment buildings only released 6% to 9%). Improvements predicted for single family homes were intermediate between apartments and non-residential buildings. Improvements in the Regional Exposure metric are larger, ranging from a factor of 2.5x to 10x for residences (when using MERV 7 to 12 rated filters) and up to 25x for large apartments with MERV 14 or 15 rated filters. The results of our modeling analysis are broadly consistent with the limited experimental data and modeling results available in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗

DOE OSTI2026

SILICON CARBIDE THERMOMETRY USING RAMAN SPECTROSCOPY ON IRRADIATED TRISO PARTICLES

Silicon carbide (SiC) passive thermometry has emerged as a promising post-irradiation examination (PIE) technique for estimating irradiation temperature near the end of irradiation. While dilatometry techniques have been traditionally used to analyze prismatic samples after irradiation, Raman spectroscopy has recently been shown to provide comparable results by analyzing Raman-active phonon modes. In this work, Raman spectroscopy has been applied to the SiC layer of cross-sectioned irradiation tristructural isotropic (TRISO) particles from the AGR-5/6/7 experiment to evaluate the feasibility of particle-scale passive thermometry. Two-dimensional Raman mapping was used to measure the position of the SiC longitudinal optical (LO) phonon, which was then converted to an apparent irradiation temperature using a previously established empirical correlation. Particles from AGR-5/6/7 Compacts 2-2-1 and 5-1-3 were selected as their calculated time-averaged, volume-averaged (TAVA) temperatures (828°C and 706°C, respectively) fall within the range of the sensitivity of the experimental approach. The use of SiC thermometry was anticipated to confirm or highlight potential deviations from calculated end-of-life TAVA temperature across compacts. Four particles from Compact 2-2-1 and two particles from Compact 5-1-3 were selected based on 110mAg inventory (either some measurable activity or below the minimum detection limit) which is commonly used as an indicator of in-pile temperature variation of particles within the same compact. Across every particle selected it was determined that the average LO peak position was located around 968 cm-1 to 969 cm-1. Using the previously determined empirical correlation, this corresponds to an irradiation temperature around 950°C to 966°C, which does not align with the reported TAVA temperature values. This discrepancy in apparent irradiation temperatures likely reflects a combination of uncertainties in the calculated particle temperatures and differences in irradiation history between the present specimens and those used to establish the empirical Raman calibration such as neutron flux (damage rate) and SiC microstructure (as fabricated and irradiated).

Vawdrey, Josh [ORNL]↗

DOE OSTI2017

PARTICLE SCATTERING OFF OF RIGHT-HANDED DISPERSIVE WAVES

Resonant scattering of fast particles off low frequency plasma waves is a major process determining transport characteristics of energetic particles in the heliosphere and contributing to their acceleration. Usually, only Alfvén waves are considered for this process, although dispersive waves are also present throughout the heliosphere. We investigate resonant interaction of energetic electrons with dispersive, right-handed waves. For the interaction of particles and a single wave a variable transformation into the rest frame of the wave can be performed. Here, well-established analytic models derived in the framework of magnetostatic quasi-linear theory can be used as a reference to validate simulation results. However, this approach fails as soon as several dispersive waves are involved. Based on analytic solutions modeling the scattering amplitude in the magnetostatic limit, we present an approach to modify these equations for use in the plasma frame. Thereby we aim at a description of particle scattering in the presence of several waves. A particle-in-cell code is employed to study wave–particle scattering on a micro-physically correct level and to test the modified model equations. We investigate the interactions of electrons at different energies (from 1 keV to 1 MeV) and right-handed waves with various amplitudes. Differences between model and simulation arise in the case of high amplitudes or several waves. Analyzing the trajectories of single particles we find no microscopic diffusion in the case of a single plasma wave, although a broadening of the particle distribution can be observed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Modeling Particle Acceleration and Release from Solar Eruptions

Determining the relative contribution of solar flares versus coronal mass ejections in large solar energetic particle (SEP) events is a long-standing problem. Flare-accelerated particles may travel through complex magnetic fields in the eruption region and escape into interplanetary space, thereby contributing to large SEP events. The process by which flare accelerated particles are released into the heliosphere is poorly understood and yet is critical to advancing our understanding of SEPs. In this work, we address the release problem by solving the focused transport equation in the context of a 2.5D ARMS magnetohydrodynamic simulation of a breakout coronal mass ejection (CME)/flare event. We find that particles accelerated by flare reconnection can be released into interplanetary space through interchange reconnection between closed and open field lines. These particles can contribute directly to SEP events and may become an important seed population for further acceleration by CME-driven shocks. Additionally, we find that the energetic particle fluxes in the inner heliosphere remain elevated for an extended period, allowing them to contribute to SEP acceleration by subsequent CMEs. This study represents the first direct particle modeling of how flare-accelerated particles can contribute to major SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Solid Particle Number and Mass Emissions from Lean and Stoichiometric Gasoline Direct Injection Engine Operation

In this work, engine-out particle mass (PM) and particle number (PN) emissions were experimentally examined from a gasoline direct injection (GDI) engine operating in two lean combustion modes and one stoichiometric mode with a fuel of known properties. Ten steady state operating points, two constant speed load steps, and an engine cold start were examined. Results showed that solid particles emitted from the engine under steady state stoichiometric conditions had a uniquely broad size distribution that was relatively flat between the diameters of 10 and 100 nm. In most operating conditions, lean homogenous modes can achieve lower particle emissions than stoichiometric modes while improving engine thermal efficiency. Alternatively, lean stratified operating modes resulted in significantly higher PN and PM emissions than both lean homogeneous and stoichiometric modes with increased efficiency only at low engine load. Stoichiometric load steps showed minimal soot emissions while ash-mode emissions spike dramatically due to oil consumption caused by piston ring adjustment. Correlation of PN to PM for steady state stoichiometric cases was in good agreement with that reported from multiple prior studies for both diesel and stoichiometric GDI engines. However, the lean cases resulted in higher PN to PM ratios indicating more small particles per unit mass. High ash particle concentration, especially in lean operation illustrates that oil control is important for mitigating impacts on downstream gasoline particulate filters (GPF) from which ash particles are cannot be removed during filter regeneration. Further research is necessary to elucidate the origin and composition of sub-23 nm particles from GDI engines and their effect on aftertreatment technology if lean homogeneous modes are to be employed.

Bock, Noah↗

DOE OSTI2026

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

FOS: Physical sciences↗