Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle interactions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2025

Moment-preserving Monte-Carlo Coulomb collision method for particle codes

Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. In conclusion, the efficacy of the model is illustrated with various test problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Relativistic approach to manipulating angular distribution of charged particles via kinetic equations

Deflection angles of charged particles interacting with materials play a critical role in various plasma applications. The development of a mathematically well-posed kinetic collision operator that accounts for deflection angles of strong Coulomb interactions remains a fundamental open problem. This paper presents a relativistic method for modifying the electromagnetic field in an anisotropic and adjustable manner to manipulate a system of charged particles, specifically by the transfer of angular momentum from a superluminal wave source to particles at specific times and locations. The method provides a mechanism to influence the scattering outcomes of strong interactions by manipulating the angular distribution of particles, and thus the deflection angles of their interactions with a solid surface, without requiring detailed knowledge of the kinetic collision operator. To this end, we demonstrate how a specific type of singularity, generated by Maxwell's equations for a superluminal wave source at the boundary of the plasma, can modify the electromagnetic field in a highly directional manner. The proposed method can lead to the development of novel approaches for controlling interactions of charged particles with a material in plasma systems. Published by the American Physical Society 2025

Moini, Nima (ORCID:0009000929568824)↗
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2025

Detailed characterization of runaway electron driven whistler waves in low-density DIII-D discharges

RE-driven whistler waves during quiescent DIII-D shots have been investigated further. The waves are confirmed to be mostly perpendicularly propagating and are observed for the first time with frequencies up to 700 MHz. Phase-spectral analysis has been used to infer their toroidal mode numbers, n, which are expected to scale with the wavenumber, k, of the mode. Though we derive a theoretical scaling of k ≈ 4n, the measured mode numbers are found to exhibit a very weak dependence on k. In addition, increases in synchrotron emission have been found to consistently lag whistler wave bursts by roughly 3–5 ms, suggesting the waves are causing pitch-angle scattering, since the emitted synchrotron radiation is a strong function of the REs' perpendicular energy. The stronger the wave bursts, the greater the subsequent increase in synchrotron emission. A predator-prey model is used to describe these nonlinear wave-particle interactions, from which the wave damping rates and the loss parameter can be inferred. The damping rates are found to be of the order of (1.6 ± 0.8) × 10 4 /s, and the unitless loss parameter is found to be approximately 2, suggesting that the loss mechanism is diffusive. These observations will serve to validate models of RE-driven waves in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

High‐Speed Embedded Ink Writing of Anatomic‐Size Organ Constructs

Embedded ink writing (EIW) is an emerging 3D printing technique that fabricates complex 3D structures from various biomaterial inks but is limited to a printing speed of ∼10 mm s −1 due to suboptimal rheological properties of particulate-dominated yield-stress fluids when used as liquid baths. In this work, a particle-hydrogel interactive system to design advanced baths with enhanced yield stress and extended thixotropic response time for realizing high-speed EIW is developed. In this system, the interactions between particle additive and three representative polymeric hydrogels enable the resulting nanocomposites to demonstrate different rheological behaviors. Accordingly, the interaction models for the nanocomposites are established, which are subsequently validated by macroscale rheological measurements and advanced microstructure characterization techniques. Filament formation mechanisms in the particle-hydrogel interactive baths are comprehensively investigated at high printing speeds. To demonstrate the effectiveness of the proposed high-speed EIW method, an anatomic-size human kidney construct is successfully printed at 110 mm s −1 , which only takes ∼4 h. This work breaks the printing speed barrier in current EIW and propels the maximum printing speed by at least 10 times, providing an efficient and promising solution for organ reconstruction in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Towards a Public Event Display for DUNE

The Deep Underground Neutrino Experiment (DUNE) is a next generation long baseline neutrino experiment based at Fermilab, with a near detector near the beam target and a Far Detector (FD) in South Dakota. As the experiment prepares for its first data runs, creating pathways for public engagement and data transparency is essential. We present the first-ever DUNE event display designed for public outreach and education. Developed using data from the ProtoDUNE detectors at the CERN Neutrino Platform, this tool provides an intuitive and interactive interface that allows non-experts to visualise and explore particle interactions in a Liquid Argon Time Projection Chamber (LArTPC). By translating raw experimental data into a browser-accessible format, we establish the essential infrastructure for DUNE’s pathway to open data. This talk will detail the technical development of the display, its current implementation with ProtoDUNE data, and the strategic roadmap for integrating it into DUNE’s long-term open-access framework.

Sabater, Eva [U. Sussex (main)] (ORCID:00090001748↗
DOE OSTI2024

DarkSide-20k sensitivity to light dark matter particles

The dual-phase liquid argon time projection chamber is presently one of the leading technologies to search for dark matter particles with masses below 10 GeV c$^{−2}$. This was demonstrated by the DarkSide-50 experiment with approximately 50 kg of low-radioactivity liquid argon as target material. The next generation experiment DarkSide-20k, currently under construction, will use 1,000 times more argon and is expected to start operation in 2027. Based on the DarkSide-50 experience, here we assess the DarkSide-20k sensitivity to models predicting light dark matter particles, including Weakly Interacting Massive Particles (WIMPs) and sub-GeV c$^{−2}$ particles interacting with electrons in argon atoms. With one year of data, a sensitivity improvement to dark matter interaction cross-sections by at least one order of magnitude with respect to DarkSide-50 is expected for all these models. A sensitivity to WIMP–nucleon interaction cross-sections below 1 × 10$^{−42}$ cm$^{2}$ is achievable for WIMP masses above 800 MeV c$^{−2}$. With 10 years exposure, the neutrino fog can be reached for WIMP masses around 5 GeV c$^{−2}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Overview of Geant4 Physics

Geant4 is a powerful toolkit for simulating particle interactions across a wide energy range (0 to ~TeV). It offers comprehensive physics processes for various particle types, allowing users to define custom processes using the G4VProcess interface. Major categories include electromagnetic, hadronic, decay, parameterized processes, and transportation. A critical feature is the PhysicsList, which specifies particles and processes for simulations. Users can choose between G4VUserPhysicsList for simpler environments or G4VModularPhysicsList for more complex scenarios, leveraging pre-defined or custom G4VPhysicsConstructors. Additionally, Geant4 provides reference physics lists and the G4PhysicsListFactory for well-established applications, facilitating simulation development in diverse domains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Physics-tailored machine learning reveals unexpected physics in dusty plasmas

Dusty plasma is a mixture of ions, electrons, and macroscopic charged particles that is commonly found in space and planetary environments. The particles interact through Coulomb forces mediated by the surrounding plasma, and as a result, the effective forces between particles can be nonconservative and nonreciprocal. Machine learning (ML) models are a promising route to learn these complex forces, yet their structure should match the underlying physical constraints to provide useful insight. Here, we demonstrate and experimentally validate an ML approach that incorporates physical intuition to infer force laws in a laboratory dusty plasma. Trained on 3D particle trajectories, the model accounts for inherent symmetries, nonidentical particles, and learns the effective nonreciprocal forces between particles with exquisite accuracy (R 2 > 0.99). We validate the model by inferring particle masses in two independent yet consistent ways. The model’s accuracy enables precise measurements of particle charge and screening length, identifying large deviations from common theoretical assumptions. Our ability to identify unknown physics from experimental data demonstrates how ML-powered approaches can guide new routes of scientific discovery in many-body systems. Furthermore, we anticipate our ML approach to be a starting point for inferring laws from dynamics in a wide range of many-body systems, from colloids to living organisms.

Science & Technology - Other Topics↗
DOE OSTI2025

Selection Algorithm for electron neutrino charged current interactions in SBND

The Short Baseline Neutrino (SBN) program at Fermilab is a joint proposal by three experimental collaborations primarily for the investigation of the cause behind the low-energy electron-like event excess observed by the MiniBoone experiment. This dissertation focuses on the near detector of the project, named the Short Baseline Near Detector (SBND), a Liquid Argon Time Projection Chamber apparatus which will conduct searches for sterile neutrinos in the mass range of 1 ${eV}^2/{c}^4$, as well as provide cross-section measurements for neutrino interactions in argon and perform other beyond the standard model studies.\\ \indent As is the case for all detectors in the program, the SBND will use the Booster Neutrino Beam as its source, which will provide it with both muon and electron neutrinos. Given that the ability to discern between the neutrino flavors will be crucial to the fulfillment of the detector's physics goals, the objective of this work is to provide the collaboration with a tool capable of doing so. As such, we here present the development process for an inclusive selection algorithm for the identification of electron neutrino charged current (CC) events regardless of their interaction channel. This is done through a combination of traditional techniques, such as the implementation of cuts on the reconstructed interaction properties, with the use of the Convolutional Visual Network, a machine learning algorithm capable of classifying particle interactions through the analysis of the topology of their final states. With this approach, we have developed a selection process that is capable of identifying $\nu_e$ CC interactions across a wide range of topologies with 34.4\% efficiency, as well as a purity of 91.2\%, making it especially promising for use in cross section studies.

Freire, Hector Moya [ABC Federal U.] (ORCID:000900↗
DOE OSTI2025

Boron Nitride-Driven Strengthening of Aluminum Composites via Friction Stir Processing

Friction stir welding and processing (FSW/P) has emerged as an effective solid-state joining technique for fabricating metal matrix composites (MMCs), offering improved mechanical properties through refined microstructural evolution. In this study, an aluminum-boron nitride nanoparticle (Al-BNNP) composite was synthesized via FSW, and its indentation-based mechanical properties were systematically evaluated. Microhardness mapping across the weld cross-section revealed a progressive increase in hardness toward the stir zone (SZ), attributed to severe plastic deformation, dynamic recrystallization (DRX), and the reinforcing effect of BNNPs. Profilometry-based indentation plastometry (PIP) inferred yield strength (YS) demonstrates a 47.8% increase compared to the base metal (BM) and a 75% improvement compared to FSP pure aluminum reported in literature. This enhancement is attributed to strengthening mechanisms, including grain boundary pinning, load transfer, and increased dislocation density. The strain rate sensitivity (SRS) measurements at the nanoscale demonstrated a substantial decrease in the SZ, correlated with ultrafine grain structures and strong BNNP-matrix interactions. Activation volume analysis revealed a significant reduction in the SZ, suggesting that dislocation motion is increasingly restricted by dislocation-dislocation and dislocation-particle interactions. These findings suggest that incorporating BNNPs in FSW/P enables tailoring the microstructure without thermal degradation of the secondary particles, thereby significantly enhancing the mechanical performance of aluminum composites, particularly for structural applications in aerospace and automotive industries.

Aluminum↗
DOE OSTI2026

Precision spectroscopy of nuclear decays using quantum optomechanical sensors

This project developed and demonstrated a fundamentally new way to detect individual nuclear decays: rather than capturing the energy that decay products deposit in a detector, we measure the tiny mechanical recoil of the entire micron- or nanometer-sized particle in which the decaying nucleus is embedded. Because momentum is conserved, this approach is sensitive even to neutral, weakly interacting particles, including neutrinos, that escape conventional detectors. During the award, the Yale group reported the first-ever mechanical detection of single nuclear decays, a result featured widely in the scientific press, and pushed the sensitivity of smaller levitated nanoparticles into the quantum measurement regime, reaching an impulse resolution within a factor of five of the fundamental Standard Quantum Limit (SQL), good enough in principle to detect the momentum kicked to the particle by a single emitted neutrino. In parallel, the LBNL group developed the theory of quantum-enhanced (sub-SQL) readout tailored to this experiment, showing how squeezed light can push the sensitivity below the SQL. Together these results establish levitated optomechanical sensors as a new tool for precision nuclear decay spectroscopy, with further applications in neutrino physics, quantum metrology, and nuclear forensics and safeguards.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Investigating the Impacts of Aqueous-Phase Processing on Organic Aerosol Chemical Climatology Using ARM and ASR Observations

This project improved understanding of how atmospheric aerosol particles form and evolve, with a focus on the role of water-driven (aqueous-phase) chemical reactions in the atmosphere. These processes occur in clouds, fog, and humid air and can significantly change the composition and properties of airborne particles, known as aerosols, which influence air quality and climate. By combining field measurements, laboratory experiments, and advanced analytical techniques, the project identified key chemical signatures that allow scientists to distinguish particles formed through aqueous processes from those formed in the gas phase. Observations from multiple environments, including wildfire smoke, urban regions, and cloud-influenced areas, show that aqueous chemistry is an important pathway for particle formation and aging. The project also developed new measurement approaches using uncrewed aerial systems (UAS) to capture how aerosol composition varies with altitude, providing critical insights into how particles interact with clouds. In addition, new data analysis frameworks were created to better interpret long-term aerosol measurements and improve characterization of particle sources and transformations. These results have been integrated into a global database of aerosol measurements and used to support atmospheric modeling efforts. Overall, the project provides important tools and knowledge to improve predictions of how aerosols affect climate and air quality, particularly through their interactions with radiation and clouds.

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

Decoherence and Brownian motion of a polarizable particle near a medium

Optically levitated nanoparticles are ideal experimental testbeds for investigating macroscopic superpositions and microscopic thermodynamics. Integrating such levitated nanoparticles with photonic structures can enable strong coupling between their center-of-mass motion and guided photonic modes, facilitating enhanced control and probing of their motion. When coupling a particle to a photonic structure, such as a waveguide, the effects of fluctuations become prominent at nanoscales. In this work, we analyze the classical and quantized center-of-mass motion of a polarizable particle interacting with the fluctuations of the electromagnetic field in the presence of a medium. We derive a position localization master equation for the particle's quantized center of mass, and examine its classical center-of-mass momentum diffusion, elucidating correspondences between classical and quantum Brownian motion of polarizable particles near media. We study the decoherence rate of the particle in the presence of a planar surface as a function of temperature and distance from the surface, comparing it to common sources of decoherence. Our results are pertinent to experiments aimed at preparing levitated nanospheres in macroscopic quantum states and investigating their Brownian dynamics.

Brownian motion↗