Engineering Papers⌕ Search

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.

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

At least 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shock Propagation and Associated Particle Acceleration in the Presence of Ambient Solar-Wind Turbulence

The topic of this review paper is on the influence of solar wind turbulence on shock propagation and its consequence on the acceleration and transport of energetic particles at shocks. As the interplanetary shocks sweep through the turbulent solar wind, the shock surfaces fluctuate and ripple in a range of different scales. We discuss particle acceleration at rippled shocks in the presence of ambient solar-wind turbulence. This strongly affects particle acceleration and transport of energetic particles (both ions and electrons) at shock fronts. In particular, we point out that the effects of upstream turbulence is critical for understanding the variability of energetic particles at shocks. Moreover, the presence of pre-existing upstream turbulence significantly enhances the trapping near the shock of low-energy charged particles, including those near the thermal energy of the incident plasma, even when the shock propagates normal to the average magnetic field. Pre-existing turbulence, always present in space plasmas, provides a means for the efficient acceleration of low-energy particles and overcoming the well known injection problem at shocks.

, solar wind↗
DOE OSTI2025

PREFACE: A Search for Long-Lived Particles at the Large Hadron Collider

The Standard Model (SM) fails to explain many problems (neutrino masses, dark matter, and matter–antimatter asymmetry, among others) that may be resolved with new particles beyond the SM. No observation of such new particles may be explained either by their exceptionally high mass or by considerably small coupling to SM particles. The latter case implies relatively long lifetimes. Such long-lived particles (LLPs) then to have signatures different from those of SM particles. Searches in the “central region” are covered by the LHC general purpose experiments. The forward small angle region far from the interaction point (IP) is unexplored. Such particles are expected to have the energy as large as E = O(1 TeV) and Lorentz time dilation factor γ = E/m ≈ 10 2 –10 3 (with m the particle mass) hence long enough decay distances. A new class of specialized LHC detectors dedicated to LLP searches has been proposed for the forward regions. Among these experiments, FASER is already operational, and FACET is under consideration at a location 100 m from the LHC IP5 (the CMS detector intersection). However, some features of FACET require a specially enlarged beam pipe, which cannot be implemented for LHC Run 4. In this study, we explore a simplified version of the proposed detector PREFACE compatible with the standard LHC beam pipe in the HL-LHC Run 4. Realistic G EANT 4 simulations are performed and the background is evaluated. An initial analysis of the physics potential with the PREFACE geometry indicates that several significant channels could be accessible with sensitivities comparable to FACET and other LLP searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Investigating Particle Size‐Dependent Redox Kinetics and Charge Distribution in Disordered Rocksalt Cathodes

Abstract Understanding how various redox activities evolve and distribute in disordered rocksalt oxides (DRX) can advance insights into manipulating materials properties for achieving stable, high‐energy batteries. Herein, the authors present how the reaction kinetics and spatial distribution of redox activities are governed by the particle size of DRX materials. The size‐dependent electrochemical performance is attributed to the distinct cationic and anionic reaction kinetics at different sizes, which can be tailored to achieve optimal capacity and stability. Overall, the local charged domains in DRX particles display random heterogeneity caused by the isotropic delithiation pathways. Owing to the kinetic limitation, the micron‐sized particles exhibit a holistic “core‐shell” charge distribution, whereas sub‐micron particles show more uniform redox reactions throughout the particles and ensembles. Sub‐micron DRX particles exhibit increasing anionic redox activities yet inferior cycling stability. In summary, engineering particle size can effectively modulate how cationic and anionic redox activities evolve and distribute in DRX materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Acoustic‐Modulated, Selective Particle Deposition for Muti‐Material Additive Manufacturing

Powder‐based additive manufacturing (AM) is an important technology in fabricating complex 3D products. In powder‐based AM, a powder deposition system plays a key role in determining the resolution, density, and properties of printed parts. An acoustic‐actuated particle deposition approach, which utilizes acoustic waves to actuate a glass nozzle, enables the deposition of particles in a cost‐effective way. However, the capability to programmably deposit different particles using acoustic waves has yet to be fully demonstrated. Herein, a selective particle deposition approach via acoustic modulation for muti‐material AM is demonstrated. Sequential and parallel depositions of two types of particles by modulating the signal inputs and resonant frequencies of piezo transducers are demonstrated. An AM process is demonstrated by integrating this acoustic‐actuated particle deposition approach with a heat‐curing mechanism. Several products from this system are produced as well. It is expected that this particle deposition approach can be used in a wide spectrum of applications, such as direct manufacturing of multi‐functional materials/devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-Reciprocity, Metastability, and Dynamic Reconfiguration in Co-Assembly of Active and Passive Particles

Living organisms often exhibit non-reciprocal interactions where the forces acting on the objects are not equal in magnitude or opposite in direction. The combination of reciprocal and non-reciprocal interactions between synthetic building blocks remains largely unexplored. Here, out-of-equilibrium assemblies of non-motile isotropic passive and metal-patched motile active particles are formed by overlapping bulk interactions with directed self-propulsion. An external alternating current (AC) electric field generates concurrent dipolar and induced-charge electrophoretic forces between the particles which are evaluated using microscopy. The interaction force measurements allow to determine the degree of reciprocity in interactions, which is tunable by designing the active particle and its trajectory. While linearly-propelled active particles evade assembly with passive particles, helically propelled active particles form active-passive clusters with dynamic reconfiguration and long-lived metastability. Large clusters display programmable fluctuations and reconfigurability by controlling the fraction of active particles. The study establishes principles of integrating reciprocal and non-reciprocal interactions in guided colloidal assembly of reconfigurable metastable structures.

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

Energy flow networks: deep sets for particle jets

A key question for machine learning approaches in particle physics is how to best represent and learn from collider events. As an event is intrinsically a variable-length unordered set of particles, we build upon recent machine learning efforts to learn directly from sets of features or “point clouds”. Adapting and specializing the “Deep Sets” framework to particle physics, we introduce Energy Flow Networks, which respect infrared and collinear safety by construction. We also develop Particle Flow Networks, which allow for general energy dependence and the inclusion of additional particle-level information such as charge and flavor. These networks feature a per-particle internal (latent) representation, and summing over all particles yields an overall event-level latent representation. We show how this latent space decomposition unifies existing event representations based on detector images and radiation moments. To demonstrate the power and simplicity of this set-based approach, we apply these networks to the collider task of discriminating quark jets from gluon jets, finding similar or improved performance compared to existing methods. We also show how the learned event representation can be directly visualized, providing insight into the inner workings of the model. These architectures lend themselves to efficiently processing and analyzing events for a wide variety of tasks at the Large Hadron Collider. Implementations and examples of our architectures are available online in our EnergyFlow package.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Settling of two-way momentum and energy coupled particles subject to Boussinesq and non-Boussinesq heating

This work establishes a procedure to accurately compute heat transfer between an Eulerian fluid and Lagrangian point-particles. Recent work has focused on accurately computing momentum transfer between fluid and particles. The coupling term for momentum involves the undisturbed fluid velocity at the particle location which is not directly accessible in the simulation. Analogously, in the context of thermal coupling, the undisturbed fluid temperature at the particle location is not directly accessible in simulations and must be estimated. In this paper, we develop a scheme to accurately estimate the undisturbed fluid temperature of a point-particle exchanging thermal energy with a surrounding fluid. Furthermore, the temperature disturbance is correlated with the enhanced temperature curvature in the vicinity of the particle and is formally valid in the low heating, low convection limit. We conduct extensive verification of the correction procedure for a settling particle subject to radiation. This setup allows the simultaneous testing of thermal and momentum corrections. By considering equations of drag and Nusselt number extended to finite Péclet and Boussinesq numbers, we establish a large range over which the correction procedure can be applied.

42 ENGINEERING↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Particle Acceleration by Magnetic Reconnection in Geospace

Particles are accelerated to very high, non-thermal energies during explosive energy-release phenomena in space, solar, and astrophysical plasma environments. While it has been established that magnetic reconnection plays an important role in the dynamics of Earth’s magnetosphere, it remains unclear how magnetic reconnection can further explain particle acceleration to non-thermal energies. Here we review recent progress in our understanding of particle acceleration by magnetic reconnection in Earth’s magnetosphere. With improved resolutions, recent spacecraft missions have enabled detailed studies of particle acceleration at various structures such as the diffusion region, separatrix, jets, magnetic islands (flux ropes), and dipolarization front. With the guiding-center approximation of particle motion, many studies have discussed the relative importance of the parallel electric field as well as the Fermi and betatron effects. However, in order to fully understand the particle acceleration mechanism and further compare with particle acceleration in solar and astrophysical plasma environments, there is a need for further investigation of, for example, energy partition and the precise role of turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2019

Predicting effective thermal and elastic properties of cementitious composites containing polydispersed hollow and core-shell micro-particles

Micro- and nano-sized hollow and core-shell particles (CSPs) have attracted tremendous interests in developing functional cementitious composites and concretes. In this paper, a general method is developed to predict the thermal and elastic properties of cementitious composites containing core-shell and hollow micro-particles. The model follows a two-stage homogenization process – CSP inclusions together with their surrounding interfacial transition zone (ITZ) are first treated as equivalent solid particles, and then the homogenized properties of the composite system are obtained using numerical (i.e., finite element) or analytical (Mori-Tanaka) approaches. Compared with other homogenization methods, the proposed model is easy-to-use and versatile. The numerical model, validated by experimental results, is then used to guide the design of functional cementitious composites. The results show that CSP particle size (or relative shell thickness), volume fraction, shell property, and the size of the ITZ have significant impacts on the effective thermal and elastic properties of cementitious composites; whereas the particle size distribution pattern has relatively minor influence. In addition, the effective property of the composite is more sensitive to the particle property change when CSP inclusions are much stiffer or more thermally conductive than the matrix. Lastly, a scoping study is conducted to elucidate the relationship between the effective thermal conductivity and effective elastic moduli for cementitious composites containing CSP additives with different shell materials, particle sizes, and volume concentrations etc.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Kinetic modeling of the formation and growth of inorganic nano-particles during pulverized coal char combustion in O 2 /N 2 and O 2 /CO 2 atmospheres

In this formation of nano-particles during coal char combustion, the vaporization of inorganic components in char and the subsequent homogeneous particle nucleation, heterogeneous condensation, coagulation, and coalescence play decisive roles. Furthermore, conventional measurements cannot provide detailed information on the dynamics of nano-particle formation and evolution, In this study, a sophisticated intrinsic char kinetics model that considers ash effects (including ash film formation, ash dilution, and ash vaporization acting in tandem), both oxidation and gasification by CO 2 and H 2 O, homogeneous particle nucleation, heterogeneous vapor condensation, coagulation, and and coalescence mechanisms is developed and used to compare the temporal evolution of the number and size of nano-particles during coal char particle combustion as a function of char particle size, ash content, and oxygen content in O 2 /N 2 and O 2 /CO 2 atmospheres .

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

In situ observations of reacting Al/Fe 2 O 3 thermite: Relating dynamic particle size to macroscopic burn time

Time-resolved x-ray imaging is used to capture high-resolution images of the Al/Fe 2 O 3 thermite reaction as it proceeds in an extended burn tube. The observed thermite reaction products are generally large, multi-phase composite particles. Image processing algorithms are developed to extract a continuous record of composite particle size as a function of time and distance from the ignition point. Small particles are observed immediately after ignition, followed by a rapid increase in particle size and a gradual decline toward late times. Correlating the dynamic particle size with visible light emission from the reaction, we propose a two-stage mechanism where the combination of early ignition and delayed expansion causes the material closest to the ignition point to undergo substantial coalescence before it begins moving down the tube. We also identify two mechanisms that may contribute to particle size refinement at late times, both based on the collapse of large gas bubbles that are observed preferentially in particles at earlier times. A new approach to extracting burn times from luminous intensity data is explored, and we find that a more nuanced analysis of burn time and luminous intensity may help to improve reactivity assessments when data from in situ studies is unavailable.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗
DOE OSTI2018

Investigation and quantification of flow unsteadiness in shock-particle cloud interaction

This work aims to study the interaction of a shock wave with a cloud of particles to quantify flow unsteadiness and velocity fluctuations using particle-resolved direct numerical simulation (PR-DNS). Three cases are studied, with each case revealing one aspect of the intricate flow phenomena involved in this interaction. The unsteady interaction of a shock wave with a transverse array of particles reveals the origin of unsteadiness under the effect of mutual wave-wave and wave-wake interactions between the particles. In the second case, the interaction of a shock with a particle cloud is studied, with a focus on the interaction of the complex wave system with the vortical structure. A budget analysis of the vorticity equation reveals the sources of strong unsteadiness in the particle cloud. A detailed analysis of the velocity fluctuation and kinetic energy in the fluctuating motion is performed to ascertain the importance of the velocity fluctuations that arise from the strong unsteadiness. Here, an analogous analysis is presented, in the third case, for a gradually-induced flow on the same particle cloud along with a comparison to the shock induced case to assess the impulsive effect of shock on intensity of the fluctuating field statistics.

42 ENGINEERING↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A second-order accurate immersed boundary-lattice Boltzmann method for particle-laden flows

A new immersed boundary-lattice Boltzmann method (IB-LBM) is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The immersed boundary method (IBM) recently developed by Breugem (2012) [19] is adopted in the present method, development including the retraction technique, the multi-direct forcing method and the direct account of the inertia of the fluid contained within the particles. The present IB-LBM is, however, formulated with further improvement with the implementation of the high-order Runge–Kutta schemes in the coupled fluid–particle interaction. The major challenge to implement high-order Runge–Kutta schemes in the LBM is that the flow information such as density and velocity cannot be directly obtained at a fractional time step from the LBM since the LBM only provides the flow information at an integer time step. This challenge can be, however, overcome as given in the present IB-LBM by extrapolating the flow field around particles from the known flow field at the previous integer time step. The newly calculated fluid–particle interactions from the previous fractional time steps of the current integer time step are also accounted for in the extrapolation. The IB-LBM with high-order Runge–Kutta schemes developed in this study is validated by several benchmark applications. It is demonstrated, for the first time, that the IB-LBM has the capacity to resolve the translational and rotational motion of particles with the second-order accuracy. The optimal retraction distances for spheres and tubes that help the method achieve the second-order accuracy are found to be around 0.30 and −0.47 times of the lattice spacing, respectively. Simulations of the Stokes flow through a simple cubic lattice of rotational spheres indicate that the lift force produced by the Magnus effect can be very significant in view of the magnitude of the drag force when the practical rotating speed of the spheres is encountered. This finding may lead to more comprehensive studies of the effect of the particle rotation on fluid–solid drag laws. It is also demonstrated that, when the third-order or the fourth-order Runge–Kutta scheme is used, the numerical stability of the present IB-LBM is better than that of all methods in the literature, including the previous IB-LBMs and also the methods with the combination of the IBM and the traditional incompressible Navier–Stokes solver. - Highlights: • The IBM is embedded in the LBM using Runge–Kutta time schemes. • The effectiveness of the present IB-LBM is validated by benchmark applications. • For the first time, the IB-LBM achieves the second-order accuracy. • The numerical stability of the present IB-LBM is better than previous methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗