Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle methods”

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

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
DOE OSTI2024

Lie–Poisson Neural Networks (LPNets): Data-based computing of Hamiltonian systems with symmetries

An accurate data-based prediction of the long-term evolution of Hamiltonian systems requires a network that preserves the appropriate structure under each time step. Every Hamiltonian system contains two essential ingredients: the Poisson bracket and the Hamiltonian. Hamiltonian systems with symmetries, whose paradigm examples are the Lie–Poisson systems, have been shown to describe a broad category of physical phenomena, from satellite motion to underwater vehicles, fluids, geophysical applications, complex fluids, and plasma physics. The Poisson bracket in these systems comes from the symmetries, while the Hamiltonian comes from the underlying physics. We view the symmetry of the system as primary, hence the Lie–Poisson bracket is known exactly, whereas the Hamiltonian is regarded as coming from physics and is considered not known, or known approximately. Using this approach, we develop a network based on transformations that exactly preserve the Poisson bracket and the special functions of the Lie–Poisson systems (Casimirs) to machine precision. We present two flavors of such systems: one, where the parameters of transformations are computed from data using a dense neural network (LPNets), and another, where the composition of transformations is used as building blocks (G-LPNets). We also show how to adapt these methods to a larger class of Poisson brackets. We apply the resulting methods to several examples, such as rigid body (satellite) motion, underwater vehicles, a particle in a magnetic field, and others. The methods developed in this paper are important for the construction of accurate data-based methods for simulating the long-term dynamics of physical systems.

Data-based modeling↗
DOE OSTI2020

Pileup mitigation at CMS in 13 TeV data

With increasing instantaneous luminosity at the LHC come additional reconstruction challenges. At high luminosity, many collisions occur simultaneously within one proton-proton bunch crossing. The isolation of an interesting collision from the additional “pileup” collisions is needed for effective physics performance. In the CMS Collaboration, several techniques capable of mitigating the impact of these pileup collisions have been developed. Such methods include charged-hadron subtraction, pileup jet identification, isospin-based neutral particle “δβ” correction, and, most recently, pileup per particle identification. This paper surveys the performance of these techniques for jet and missing transverse momentum reconstruction, as well as muon isolation. The analysis makes use of data corresponding to 35.9 fb−1 collected with the CMS experiment in 2016 at a center-of-mass energy of 13 TeV. The performance of each algorithm is discussed for up to 70 simultaneous collisions per bunch crossing. Significant improvements are found in the identification of pileup jets, the jet energy, mass, and angular resolution, missing transverse momentum resolution, and muon isolation when using pileup per particle identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability of the particle-hole Pfaffian state and the 5 2 -fractional quantum Hall effect

Here we present a method for the exact construction of the fully particle-hole symmetric Pfaffian (PH-Pfaffian) ground state and its charged excitations on a sphere. We adopt the Moore-Read state, but with a nonholomorphic pairing component as in previous studies, and project it to the lowest Landau level. We study the energetics as well as other properties of these states and find that in a pure system interacting with the Coulomb forces the PH-Pfaffian cannot compete with either the Moore-Read state or its particle-hole conjugate, the anti-Pfaffian state, as an explanation for the 5/2 effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Ammonium-coordinated exchanger (ACE) functionalized silica sorbents for recovering/removing aqueous anionic contaminants

There are limited studies of functionalized silica anion exchange sorbents used for critical/heavy metal recovery/removal relative to polymeric and other inorganic materials. This work features ammonium-coordinated exchanger (ACE) anion exchange particle sorbents prepared by either acid-washing epoxy-crosslinked polyethylenimine (PEI) hydrogen bonded within/to a silica particle sorbent (two-step method) or reacting a di-chlorinated crosslinker, α,α-dichloro-p-xylene (DPX), with PEI within silica (single-step method). Energy dispersive X-ray spectroscopy (EDS) and infrared spectroscopy confirmed the presence of -NH 2 + ···Cl - and -NH 3 + ···Cl - groups, which removed oxyanionic species –arsenate, selenate, chromate, sulfate, phosphate, and nitrate– plus bromide from ideal solutions, authentic acid mine drainage (AMD), and authentic flue gas desulfurization (FGD) wastewater. Affinity of the anions for ACE varied across single- and mixed-element solutions. However, affinity for CrO 4 2- was among the highest in both cases. Total anion uptake by the optimized ACE, PEI-E3-HCl_1.1, reached 1.2 mmol anion/g-sorb. (0.56 mmol CrO 4 2- /g), or ∼2.3 mmol negative charge/g-sorb. This was close to the 0.52 mmol CrO 4 2- /g of a commercial anion exchange resin. Near-consistent removal of 20–80 % of each anion from FGD during an eight-cycle adsorption-desorption (1 M NaCl) test predicted good ACE viability for testing under practical conditions at larger scale.

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

Highly Durable and Active Pt/Sb-Doped SnO 2 Oxygen Reduction Reaction Electrocatalysts Produced by Atomic Layer Deposition

Platinum supported on mixed-metal oxides (MMOs) are a class of active and durable cathode catalysts for proton exchange membrane fuel cell (PEMFC) due to a combination of the high oxidative stability of the supports and strong-metal-support interactions (SMSI) that enable them to exceed the activity of Pt/C. Herein, we solve a significant remaining challenge with Pt/MMO systems, namely the relatively low surface area and porosity. This is achieved by dispersing nearly-uniform Pt clusters using atomic layer deposition (ALD) on highly conductive (6.2 S/cm) and stable antimony doped tin dioxide (ATO) support. ALD-Pt/ATO exhibited significantly higher electrochemically active surface area (ECSA) (74 m 2 /g) and oxygen reduction reaction (ORR) catalytic activity (102 mA/mg Pt at 0.9 V vs. RHE) compared to Pt/ATO synthesized using ethylene glycol (ECSA=31 m 2 /g Pt , mass activity=52 mA/mg Pt at 0.9 V vs. RHE) and formic acid reduction methods (ECSA=28 m 2 /g Pt , mass activity=46 mA/mg Pt at 0.9 V vs. RHE). Further characterization showed that wet chemical methods resulted in poorer Pt particle dispersion, poor control over Pt particle size distribution and chemical degradation of the support (during Pt deposition). Given the near-ideal Pt particle size distribution of the ALD-Pt/ATO, particle size growth and loss of ECSA was found to be minimal over the course of rigorous potential cycling. Thus, after 10,000 potential cycles between 1V and 1.5V vs. RHE, ALD-Pt/ATO and other Pt/ATOs were found to retain 100% of their initial ECSA compared to 57.6% retention for Pt/C. Upon testing in a H 2 /air PEMFC, following 1,000 potential cycles, the change in ALD-Pt/ATO performance was negligible while Pt/C exhibited a 68.2% loss of initial peak power density. Thus, ALD-Pt/ATO is an active and highly durable ORR electrocatalyst in PEMFCs under start-up-shut down conditions.

25 ENERGY STORAGE↗
DOE OSTI2024

Tutorial: The equations of electron emission and their evaluation

Electron emission and transport through and over potential barriers is an essential process requiring modeling and simulation to meet the design needs and characterization of an exceedingly broad range of technologically important devices and processes. The simulation and description of thermal, field, and photoemission, and the related concerns of space–charge affected electron flow, often make use of specialized formulations developed in the early days of quantum mechanics. Advancements in the utilization of electron sources and particularly the simulation of devices and applications using advanced particle-in-cell and trajectory methods for beam optics codes create a strong need for a pedagogical account of the emission models to ensure correct numerical evaluation of their equations. This Tutorial starts from simple phenomenological accounts and progressively builds to comprehensive models emphasizing straightforward and often rapid calculation. It recommends formulations to supplant the canonical Richardson–Laue–Dushman (thermal), Fowler–Nordheim (field), Fowler–DuBridge (photo), and Baroody (secondary) equations and provides a useful formulation of space–charge affected flow commonly described by the Child–Langmuir relation that takes into account cathode dependence on surface field.

Physics↗
DOE OSTI2024

Cooper minima in high-Z atoms: effects of correlation and relativity on np photoionization

Abstract Photoionization dipole transition matrix elements pass through a zero or attain a minimum that leaves imprints on photoionization parameters like the cross-section, angular distribution asymmetry parameter, phase shift, and photoionization time delay. This minimum is commonly known as the ‘Cooper minimum’ (CM). The CM, in general, is strongly affected by relativistic and correlation effects. Previous works investigated CM in the 6pand 5psubshell photoionization up toZ= 100 using the single-particle Dirac-Slater (DS) method. The present work extends the earlier work toZup to 120 using more accurate methods; Dirac–Hartree–Fock (DHF) which includes the relativistic effects and exchange correlations, and the relativistic random phase approximation (RRPA) which includes both initial and final state electron-electron correlations along with relativistic effects. In addition to the study of photoionization from the 6pand 5psubshells, the 4psubshell has also been investigated in the present work. To demonstrate the prominent effects in the high-Z atoms, Rn (Z= 86), Ra (Z= 88), No (Z= 102), Cn (Z = 112), Og (Z= 118), and Ubn (Z= 120) are investigated.

Optics↗
DOE OSTI2025

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗
DOE OSTI2022

Variational Adiabatic Gauge Transformation on Real Quantum Hardware for Effective Low-Energy Hamiltonians and Accurate Diagonalization

Effective low-energy theories represent powerful theoretical tools to reduce the complexity in modeling interacting quantum many-particle systems. However, common theoretical methods rely on perturbation theory, which limits their applicability to weak interactions. Here we introduce the Variational Adiabatic Gauge Transformation (VAGT), a nonperturbative hybrid quantum algorithm that can use nowadays quantum computers to learn the variational parameters of the unitary circuit that brings the Hamiltonian to either its block-diagonal or full-diagonal form. If a Hamiltonian can be diagonalized via a shallow quantum circuit, then VAGT can learn the optimal parameters using a polynomial number of runs. The accuracy of VAGT is tested through numerical simulations, as well as simulations on Rigetti and IonQ quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bioaerosols are the dominant source of warm-temperature immersion-mode INPs and drive uncertainties in INP predictability

Ice-nucleating particles (INPs) are rare atmospheric aerosols that initiate primary ice formation, but accurately simulating their concentrations and variability in large-scale climate models remains a challenge. Doing so requires both simulating major particle sources and parameterizing their ice nucleation (IN) efficiency. Validating and improving model predictions of INP concentrations requires measuring their concentrations delineated by particle type. We present a method to speciate INP concentrations into contributions from dust, sea spray aerosol (SSA), and bioaerosol. Field campaign data from Bodega Bay, California, showed that bioaerosols were the primary source of INPs between –12° and –20°C, while dust was a minor source and SSA had little impact. We found that recent parameterizations for dust and SSA accurately predicted ambient INP concentrations. However, the model did not skillfully simulate bioaerosol INPs, suggesting a need for further research to identify major factors controlling their emissions and INP efficiency for improved representation in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Measurements of multiplicity fluctuations of identified hadrons in inelastic proton–proton interactions at the CERN Super Proton Synchrotron

Measurements of multiplicity fluctuations of identified hadrons produced in inelastic p+p interactions at 31, 40, 80, and 158 $\text {Ge}\text {V}/c$ beam momentum are presented. Three different measures of multiplicity fluctuations are used: the scaled variance $\omega $ and strongly intensive measures $\Sigma $ and $\Delta $. These fluctuation measures involve second and first moments of joint multiplicity distributions. Data analysis is preformed using the Identity method which corrects for incomplete particle identification. Strongly intensive quantities are calculated in order to allow for a direct comparison to corresponding results on nucleus–nucleus collisions. The results for different hadron types are shown as a function of collision energy. A comparison with predictions of string-resonance Monte-Carlo models: Epos, Smash and Venus, is also presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2020

Hydrometeor Shape Variability in Snowfall as Retrieved from Polarimetric Radar Measurements

A polarimetric radar–based method for retrieving atmospheric ice particle shapes is applied to snowfall measurements by a scanning Ka-band radar deployed at Oliktok Point, Alaska (70.495°N, 149.883°W). The mean aspect ratio, which is defined by the hydrometeor minor-to-major dimension ratio for a spheroidal particle model, is retrieved as a particle shape parameter. The radar variables used for aspect ratio profile retrievals include reflectivity, differential reflectivity, and the copolar correlation coefficient. The retrievals indicate that hydrometeors with mean aspect ratios below 0.2–0.3 are usually present in regions with air temperatures warmer than approximately from –17° to –15°C, corresponding to a regime that has been shown to be favorable for growth of pristine ice crystals of planar habits. Radar reflectivities corresponding to the lowest mean aspect ratios are generally between –10 and 10 dBZ. For colder temperatures, mean aspect ratios are typically in a range between 0.3 and 0.8. There is a tendency for hydrometeor aspect ratios to increase as particles transition from altitudes in the temperature range from –17° to –15°C toward the ground. Furthermore, this increase is believed to result from aggregation and riming processes that cause particles to become more spherical and is associated with areas demonstrating differential reflectivity decreases with increasing reflectivity. Aspect ratio retrievals at the lowest altitudes are consistent with in situ measurements obtained using a surface-based multiangle snowflake camera. Pronounced gradients in particle aspect ratio profiles are observed at altitudes at which there is a change in the dominant hydrometeor species, as inferred by spectral measurements from a vertically pointing Doppler radar.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Geant4 Simulations on Faraday Cup Design for PIP-II Laser Wire Scanner System

The Proton Improvement Plan-II (PIP-II) accelerator upgrade at Fermilab represents a groundbreaking leap forward in high-energy physics research. This ambitious initiative involves enhancing Fermilab's accelerator complexby replacing the current linear accelerator (linac) with awarm front end (WFE) capable of accelerating H- beamsup to 2.1 MeV. Subsequently, a superconducting linac further accelerates these beams up to 800 MeV. To preciselymeasure the transverse beam profile, a combination of traditional wire scanners at the WFE section and Laser wirescanners along the superconducting linac are planned forimplementation. This investigation is centered on refiningthe Faraday cup design for the PIP-II Laser wire scannersby utilizing GEANT4, a Monte Carlo simulation toolkit.Leveraging this method enables a comprehensive analysisof particle trajectories, energy deposition, secondary particle emission, backscattering, etc., facilitating optimizationthrough adjustments to cup geometries, materials, andplacement to maximize its efficacy in beam diagnostics.

43 PARTICLE ACCELERATORS↗