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

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Solid state electrolytes

A product includes a solid state electrolyte particle coated with a coating consisting essentially of metal oxide. A method includes depositing a coating on a solid state electrolyte particle by atomic layer deposition. A method includes fabricating a product using a plurality of solid state electrolyte particles coated with a coating consisting essentially of metal oxide.

Ye, Jianchao↗
DOE OSTI2025

ASCR Workshop Position Paper: Challenges and Opportunities in High Energy Physics

High energy particle physics and cosmology concern themselves with estimating fundamental parameters of nature, such as the masses and interactions of fundamental particles like the Higgs boson and the rate of expansion of the universe. In doing so, they analyze exabyte-scale datasets, some of the largest in all of science, and face many challenges in subsequent data analysis. These challenges are shared between the two disciplines, but we focus on particle physics to highlight one specific domain. In particle physics, the standard method for estimating parameters involves performing Monte Carlo (MC) integration as a function of both parameters of interest and nuisance parameters using an expensive simulator, counting the number of observed collision events (i.i.d. samples) from an experiment in the corresponding integration domains, and forming a Poisson likelihood function. This likelihood function is then used in a Frequentist manner to construct a maximum likelihood point estimate (MLE) and confidence set for the parameters. To sufficiently populate the high-dimensional integration domains, simulators consume billions of CPU-hours annually and produce hundreds of petabytes of intermediate output data. Several techniques have been developed to: optimize definitions of the integration domains so as to be maximally sensitive to a particular subset of parameters, efficiently estimate the integrals, and build robust surrogate models by interpolating between integral evaluations at different parameter points. One can view this whole endeavor as classical Simulation-Based Inference (SBI).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Development of interatomic potential and effect of ordering on defect properties in CrMnV

Developing materials that can withstand extreme environments, such as high radiation doses and elevated temperatures, is crucial for next-generation particle accelerators, including the 2.4 MW Long-Baseline Neutrino Facility. High-Entropy Alloys have emerged as promising candidates for beam window materials due to their superior mechanical strength, corrosion resistance, and radiation tolerance. In this study, we focus on the Cr–Mn–V alloy system, developing and employing machine-learning interatomic potentials (MLIPs) to investigate the formation of an ordered phase and its influence on defect properties. Using hybrid Monte Carlo-Molecular Dynamics simulations, we observe the formation of a B2-ordered phase at lower temperatures, consistent with Density Functional Theory (DFT) predictions. Ordered structures display a bimodal distribution of migration energies and reduced mean square displacement values, indicating suppressed vacancy diffusion. Our results also show that the migration energy barrier varies based on the atomic species, with Mn and V exhibiting the highest and lowest average barriers, respectively. These findings suggest that atomic ordering inhibits defect mobility, potentially enhancing the radiation resistance of CrMnV alloys. The validated MLIP provides a reliable framework for simulations that are faster than traditional DFT while maintaining the accuracy required to study defect and ordering properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploration of LIBS as a novel and rapid elemental mapping technique of nuclear fuels in the form of surrogate TRISO particles

Laser-induced breakdown spectroscopy (LIBS) was employed to characterize coatings on surrogate fuel particles. Tri-structural isotropic (TRISO) particles are a proposed nuclear fuel alternative for high temperature reactors. These particles are constructed of a ZrO 2 kernel (as a surrogate to uranium), surrounded by an inner pyrolytic carbon layer and are surrounded by an outer carbide layer (ZrC, presented here) to act as a barrier to fission products generated during nuclear reactions. These particles are embedded within a graphite compact and housed within the reactor core. Simply put, due to their robust nature, performing elemental analysis of these particles poses a challenge. Presented here, LIBS is explored as a method for characterizing elemental constituents of these particles, with the focus being on rapid elemental mapping and depth profiling. Different from traditional elemental analysis techniques (e.g., inductively coupled plasma – based methods), LIBS is advantageous because it can directly analyze the sample surface and can detect light elements such as C and O, making it a viable technique for the analysis of small, multilayered particles as spatial elemental information is warranted in the production of these particles. In the work presented here, LIBS was successfully used for discerning small layers (30–50 μm), detecting the location of carbon and oxygen layers, providing fast 2-D mapping (<5 min per particle) and rapid depth profiling (10 s per particle).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Measurement of forward charged hadron flow harmonics in peripheral PbPb collisions at s N N = 5.02 TeV with the LHCb detector

Flow harmonic coefficients, v n , which are the key to studying the hydrodynamics of the quark-gluon plasma (QGP) created in heavy-ion collisions, have been measured in various collision systems and kinematic regions and using various particle species. The study of flow harmonics in a wide pseudorapidity range is particularly valuable to understand the temperature dependence of the shear viscosity to entropy density ratio of the QGP. This paper presents the first LHCb results of the second- and the third-order flow harmonic coefficients of charged hadrons as a function of transverse momentum in the forward region, corresponding to pseudorapidities between 2.0 and 4.9, using the data collected from PbPb collisions in 2018 at a center-of-mass energy of 5.02 TeV . The coefficients measured using the two-particle angular correlation analysis method are smaller than the central-pseudorapidity measurements at ALICE and ATLAS from the same collision system but share similar features. ©2024 CERN, for the LHCb Collaboration 2024 CERN

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nanolipoprotein particle (NLP) vaccine confers protection against Yersinia pestis aerosol challenge in a BALB/c mouse model

Introduction: Yersinia pestis is the etiological agent of plague, a disease that remains a concern as demonstrated by recent outbreaks in Madagascar. Infection with Y. pestis results in a rapidly progressing illness that can only be successfully treated with antibiotics given shortly after symptom onset. Live attenuated or whole cell inactivated vaccines confer protection against bubonic plague, but pneumonic plague has been more difficult to prevent. Novel effective subunit vaccine formulations may circumvent some of these shortfalls. Here, we compare the immunogenicity generated by an advanced subunit vaccine (F1V fusion protein) and a nanolipoprotein particle (NLP)-based vaccine. Methods: The NLP, a high-density lipoprotein mimetic, provides a nanoscale delivery platform for recombinant Y. pestis antigens LcrV (V) and F1. BALB/c mice were immunized via subcutaneous injection twice, three or four weeks apart. Four weeks later, splenocytes and sera were collected for immune profiling, and mice were challenged with aerosolized Y. pestis CO92. Results: Both formulations induced a strong IgG response against the F1 and V proteins, along with a robust memory B cell response and a balanced cell-mediated immune response as evidenced by both Th1- and Th2-related cytokines. The NLP-based vaccine induced a stronger cytokine response against F1, V, and F1V proteins relative to the F1V vaccine. As with F1V, the inclusion of Alhydrogel (Alu) in NLP vaccine formulations was critical for enhanced immunogenicity and protective efficacy. Mice that received two doses of F1:V:NLP + Alu and CpG were completely protected from a challenge with approximately eight median lethal doses of aerosolized Y. pestis CO92 and this protection confirmed the well-documented synergy between the F1 and V antigens in context of pneumonic plague. The NLPs have defined regions of polarity that facilitates the incorporation of a wide range of adjuvants and antigens with distinct physicochemical properties and are an excellent candidate platform for the development of multi-antigen vaccines.

F1↗
DOE OSTI2026

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗
DOE OSTI2022

On the perturbation effect and LET dependence of beam quality correction factors in carbon ion beams

Background: In a recent study, we reported beam quality correction factors, f Q , in carbon ion beams using Monte Carlo (MC) methods for a cylindrical and a parallel-plate ionization chamber (IC). A non-negligible perturbation effect was observed; however, the magnitude of the perturbation correction due to the specific IC subcomponents was not included. Furthermore, the stopping power data presented in the International Commission on Radiation Units and Measurements (ICRU) report 73 were used, whereas the latest stopping power data have been reported in the ICRU report 90. Purpose: The aim of this study was to extend our previous work by computing f Q correction factors using the ICRU 90 stopping power data and by reporting IC-specific perturbation correction factors. Possible energy or linear energy transfer (LET) dependence of the f Q correction factor was investigated by simulating both pristine beams and spread-out Bragg peaks (SOBPs). Methods: The TOol for PArticle Simulation (TOPAS)/GEANT4 MC code was used in this study. A 30 × 30 × 50 cm 3 water phantom was simulated with a uniform 10 × 10 cm 2 parallel beam incident on the surface. A Farmer-type cylindrical IC (Exradin A12) and two parallel-plate ICs (Exradin P11 and A11) were simulated in TOPAS using the manufacturer-provided geometrical drawings. The f Q correction factor was calculated in pristine carbon ion beams in the 150–450 MeV/u energy range at 2 cm depth and in the middle of the flat region of four SOBPs. The k Q correction factor was calculated by simulating the f Qo correction factor in a 60 Co beam at 5 cm depth. The perturbation correction factors due to the presence of the individual IC subcomponents, such as the displacement effect in the air cavity, collecting electrode, chamber wall, and chamber stem, were calculated at 2 cm depth for monoenergetic beams only. Additionally, the mean dose-averaged and track-averaged LET was calculated at the depths at which the f Q was calculated. Results: The ICRU 90 f Q correction factors were reported. The p dis correction factor was found to be significant for the cylindrical IC with magnitudes up to 1.70%. The individual perturbation corrections for the parallel-plate ICs were <1.0% except for the A11 p cel correction at the lowest energy. The f Q correction for the P11 IC exhibited an energy dependence of >1.00% and displayed differences up to 0.87% between pristine beams and SOBPs. Conversely, the f Q for A11 and A12 displayed a minimal energy dependence of <0.50%. The energy dependence was found to manifest in the LET dependence for the P11 IC. A statistically significant LET dependence was found only for the P11 IC in pristine beams only with a magnitude of <1.10%. Conclusions: The perturbation and k Q correction factor should be calculated for the specific IC to be used in carbon ion beam reference dosimetry as a function of beam quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Solubility controls on plutonium and americium release in subsurface environments exposed to acidic processing wastes

To identify the role of waste composition and sediment interactions in controlling Pu and Am mobility in contaminated sediments at Hanford, a legacy nuclear site, Pu and Am concentrations in solutions equilibrated with contaminated sediments from beneath the 216-Z-9 (Z-9) Trench were compared to the solubilities of PuO 2 materials, synthesized by methods representative of the disposed wastes, in the absence of sediments. Furthermore, this work shows that the solubilities of PuO 2 materials synthesized by different methods, and with varying particle sizes, agree with PuO 2 (am,hyd), although dissolution kinetics differed between materials. According to saturation index (SI) calculations, PuO 2 (am,hyd) is likely also controlling Pu release from sediments under conditions where phosphate concentrations are low. However, both Pu-phosphate and Am-phosphate phases, identified in SI calculations and by high resolution transmission electron microscopy, play roles in controlling release in low pH, high phosphate, shallow sediments just below the Z-9 Trench. The elevated phosphate is likely due to decomposition of tributyl phosphate from waste solutions over time. Sediments from deeper in the subsurface beneath the Z-9 Trench are less acidic and contain less phosphate, with Pu solubility likely controlled by PuO 2 (am,hyd) that precipitated following neutralization of the acidic waste stream. Controls on Am concentrations in deeper sediments are more complex and potentially involve sediment adsorption and/or release from Pu (1-x) Am x O 2 following Am in-growth. The concentrations of both Pu and Am were elevated in the colloidal fraction associated with shallow sediments, but not in PuO 2 experiments, suggesting the presence of Pu/Am pseudocolloids (e.g., Pu/Am associated with mineral colloids). However, Pu and Am association with the colloidal size fraction was not observed in deeper sediments, suggesting transport of Pu and Am to these depths beneath the Z-9 Trench was not due to colloidal transport.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Novel imaging technique for thermal neutrons using a fast optical camera

A novel imaging technique for thermal neutrons using a fast optical camera is presented. Thermal neutrons are reacted with 6 Lithium to produce a pair of 2.73 MeV tritium and 2.05 MeV alpha particles, which in turn interact in a thin layer of LYSO crystal scintillator to produce a localized flash of light. These photons are directed by a pair of lenses to a micro-channel plate intensifier, and its output is connected to the optical camera, TPX3CAM. The results from the camera are reconstructed through a custom algorithm. Various cutting parameters were found through data analysis to eliminate the background, and they were shown effective in matching the simulated rate of the neutron source. The system is fast with 40 ns decay time and allows free-space light collection, both vastly enhances flexibility of neutron detection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Importance of non-flow background on the chiral magnetic wave search

We report an observable sensitive to the chiral magnetic wave (CMW) is the charge asymmetry dependence of the π - and π + anisotropic flow difference, Δv n (A ch ). We show that, due to non-flow correlations, the flow measurements by the Q-cumulant method using all charged particles as reference introduce a trivial linear term to Δv n (A ch ). The trivial slope contribution to the triangle flow difference Δv 3 (A ch ) can be negative if the non-flow is dominated by back-to-back pairs. This can explain the observed negative Δv 3 (A ch ) slope in the preliminary STAR data. We further find that the non-flow correlations give rise to additional backgrounds to the slope of Δv 2 (A ch ) from the competition among different pion sources and from the larger multiplicity dilution to π + (π - ) at positive (negative) A ch .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Alpha Shape Analysis (ASA) Framework for Post- Clustering Property Determination in Atom Probe Tomographic Data

While application of clustering algorithms to atom probe tomography data have enabled quantification of solute clusters in terms of number density, size, and subcomposition there exist other properties (e.g., volume, surface area, and composition) that are better determined by defining an interface between the cluster and the surrounding matrix. The limitation in composition results from an ion selection step where the expected matrix ion types are omitted from the cluster search algorithm to enhance the contrast between the matrix and cluster and to reduce the complexity of the search. Previously, composition determination within solute clusters has utilized a secondary envelopment and erosion step on top of conventional methods such as maximum separation. In this work, we present a novel stochastic method that combines the particle identification fidelity of a conventional clustering algorithm with the analytical flexibility of mesh-based approaches through the generation of alpha shapes for each identified cluster. Furthermore, the corresponding mesh accounts for concave components of the clusters and determines the volume and surface area of the clusters; additionally, the mesh boundary is utilized to update the total composition according to the internal ions.

42 ENGINEERING↗
DOE OSTI2020

Development of Density Functional Tight-Binding Parameters Using Relative Energy Fitting and Particle Swarm Optimization

We provide a strategy to optimize density functional tight-binding (DFTB) parameterization for the calculation of the structures and properties of organic molecules consisting of hydrogen, carbon, nitrogen, and oxygen. We utilize an objective function based on similarity measurements and the Particle Swarm Optimization (PSO) method to find an optimal set of parameters. This objective function considers not only the common DFTB descriptors of binding energies and atomic forces but also incorporates relative energies of isomers into the fitting procedure for more chemistry-driven results. The quality in the description of the binding energies and atomic forces is measured based on the Ballester similarity index and relative energies through a similarity index induced by the Levenshtein edit distance to quantify the correct energetic order of isomers. Training and testing datasets were created to include all relevant chemical functional groups. Finally, the accuracy of this strategy is assessed, and its range of applicability is discussed by comparison against our previous parameterization. The improved performance of the new DFTB parameterization is validated with respect to the density functional theory large datasets QM-9 and ANI-1, where excellent agreement is found between the structures and properties available in these datasets, and the ones obtained with DFTB.

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

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗
DOE OSTI2020

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗
DOE OSTI2025

Direct evidence of the effect of a moderate external magnetic field on stimulated Raman scattering in the kinetic regime

We present results from an experiment carried out at the OMEGA-EP laser facility that investigated the effect of a perpendicular magnetic field on stimulated Raman scattering (SRS) and report the first direct measurement of magnetic mitigation of SRS. A 13-T magnetic field generated by pulsed-power coils was imposed on a gas jet plasma, and a novel three-picket interaction beam was used to explore SRS reflectivity for several plasma conditions within single shots. The time-resolved backscattered light shows that SRS was mitigated by the external 13-T magnetic field in the kinetic regime (kλ D ∼ 0.3, where k is the electron plasma wave's wavenumber and λ D is electron Debye length) and at a density of n e /n cr ∼ 0.10. On the other hand, we also measured an enhancement of SRS reflectivity at lower density (n e /n cr < 0.08). We discuss experimental results in the context of magnetohydrodynamic and particle-in-cell simulations that scan SRS dynamics for a variety of plasma conditions. While the experimental evidence of SRS mitigation validates prior work on the kinetic simulation of SRS in an external magnetic field, we also find that other mechanisms such as SRS rescatter can lead to an enhancement of measured SRS reflectivity in simulations of parameters relevant to our experiment.

Gas jet↗
DOE OSTI2024

Transfer functions for Q A /Q B international regulatory limits for the safe transport of radioactive materials

This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.

61 RADIATION PROTECTION AND DOSIMETRY↗