Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle method”

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,315 records · Page 73

On the measurement of shape: With applications to lunar regolith

With the renewed commitment from NASA and other commercial entities for a presence on the Moon, the importance of understanding the characteristics of lunar regolith and how to utilize it have become the target of increasing scrutiny. Much of what is known about lunar regolith was collected during and immediately after the Apollo program, however, analytical techniques and instrumentation have advanced in leaps and bounds in the subsequent decades. Specifically, dynamic image analysis systems have advanced to the point that millions of particles can have morphological characteristics automatically determined in relatively short time frames. Particle morphological data was collected on several lunar samples and a pair of widely used regolith simulants to ascertain the accuracy of these simulants and to explore statistical analysis methods of these large datasets. It is found that these morphology datasets can vary widely depending on the particle size of the particles, and simple averaging of the data skews the results heavily towards the numerically abundant size ranges, the fines. Different reporting methods are suggested to ameliorate these problems. When applied to the lunar regolith, the particles are noted to be less morphologically complex than initially suspected. Compared to the lunar material, the simulants are found to contain some more morphological variability and angular grains. Such difference is likely due to the wildly different comminution processes that these different powder systems are subjected to.

2D shape↗

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding the Mesoscale Degradation in Nickel-Rich Cathode Materials through Machine-Learning-Revealed Strain–Redox Decoupling

The degradation of nickel-rich cathode materials for lithium-ion batteries upon prolonged electrochemical cycling features a complicated interplay among electronic structure, lattice configuration, and micro-morphology. The underlying mechanism for such an entanglement of different material properties at nano- to mesoscales is fundamental to the battery performance but not well-understood yet. Here we investigate the correlation between the local redox reaction and lattice mismatch through a nano-resolution synchrotron spectro-microscopy study of LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM 811) cathode particles. With assistance from a machine-learning-based data classification method, we identify local regions that demonstrate a strain–redox decoupling effect, which can be attributed to different side reactions. Overall, our results highlight the mesoscale reaction heterogeneity in the battery cathode and suggest that particle structure engineering could be a viable approach to mitigate the chemomechanical degradation of cathode materials.

25 ENERGY STORAGE↗

Particle Sensitivity Analysis

We propose to develop a computational sensitivity analysis capability for Monte Carlo sampling-based particle simulation relevant to Aleph, Cheetah-MC, Empire, Emphasis, ITS, SPARTA, and LAMMPS codes. These software tools model plasmas, radiation transport, low-density fluids, and molecular motion. Our report demonstrates how adjoint optimization methods can be combined with Monte Carlo sampling-based adjoint particle simulation. Our goal is to develop a sensitivity analysis to drive robust design-based optimization for Monte Carlo sampling-based particle simulation - a currently unavailable capability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advancing Multiscale Simulation of Plasma-Surface Interfaces

We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TRISO PARTICLE FUEL POLISHING PROCESS COMPARISON

The investigation of the TRISO coated particles is important in terms of verification the quality of the nuclear fuel samples before placing them in the reactor core. To determine the degree of damage to the coated particles, the samples must be properly polished in order to uncover the inner layers. The poster presents a comparison of two polishing methods. Mechanical polishing method, which is used frequently for sample preparations generates additional mechanical damages in the studied fuel particle, thus directly affecting the experimental results. Due to that, a novel way to polish the samples is presented - an ionic polishing method. Based on the obtained results, it was concluded that the ion polishing method is better because the level of interference with the structures of the individual layers of the tested samples is much lower than with the mechanical method.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Radial particle-based terrestrial thermocline for high temperature thermal storage

A system or method for thermal storage includes a recess or containment unit having a first storage layer and a second storage layer comprising a permeable filler material. An intermediate layer is disposed between the storage layers. A primary well traverses the layer in the recess. The primary well is in thermal communication with the first permeable filler material and the second permeable filler material. A heat source is provided for heating an inlet fluid. An input pump is in fluid communication with the primary well and the heat source. The primary well receives heated inlet fluid from the inlet pump and injects the fluid into the second layers. The heated inlet fluid transfers heat to the respective permeable filler material radially from the primary well toward an outer periphery of the thermocline recess.

Ho, Clifford K.↗

Tunnel ionization within a one-dimensional, undriven plasma sheath

In high density, high temperature plasmas, the plasma sheath that develops can result in extremely high electric fields, on the order of tens to hundreds of V/nm. Under the right conditions, these electric fields can reach magnitudes that can increase the probability of electron tunneling ionization to occur, resulting in one or more electron-ion pairs. The presence of tunneling ionization can then modify the development of the plasma sheath, as well as properties such as the ion and electron densities and plasma potential. The tunnel ionization process for hydrogen atoms is demonstrated, in this work, as implemented in a Sandia National Laboratories, particle-in-cell code Aleph. Results are presented for the application of the tunnel ionization process to a one-dimensional, undriven plasma sheath. Additional results for cases that consider warm ions and neutrals, the inclusion of electron–neutral collisions, and the injection of neutral particles, as well as the application to various plasma devices, will be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gas Stopper Developments for Improved Purity and Intensity of Low-Energy, Rare Isotope Ion Beams (Final Technical Report)

This final technical report summarizes the work of the Michigan State University (MSU) team supported by grant # DE-SC0021423 awarded by the Office of Nuclear Physics, Department of Energy. Objectives: The successful fulfillment of the FRIB science mission hinges on ensuring the availability of fast, stopped, and reaccelerated beams consisting of rare isotopes. This project's research and development focus was dedicated to supporting the advancement and creation of a cutting-edge linear gas stopper. The primary aim is to efficiently convert the high-intensity fast beams of rare isotopes provided by FRIB into high-quality, low-energy beams. These beams are essential for conducting stopped beam experiments or for subsequent reacceleration. The overarching goal is to advance technology, aiming to increase the beam rate capability of the linear gas stopper for medium-to-heavy-mass rare isotopes by more than tenfold compared to the currently most effective gas stopper in operation, and to improve the purity of the delivered beams. Project Description: The existing technology employed in gas stopping devices designed for low-energy, rare-isotope beams presents limitations in both the purity of extracted beams and the intensities of injected beams. These limitations are incompatible with the requirements of the recently commissioned rare isotope beam facility, FRIB. Our research and development efforts, aligned with the previously outlined objectives, focused on addressing the most critical aspects for enhancing beam-rate capability and purity. Specifically, advanced particle-in-cell simulations were developed and integrated into a simulation pipeline to explore the efficacy of multi-layer RF carpets on increasing ion transport efficiency with high incoming beam rates that generate space charge fields which can limit it. We also explored the possibility of using a collision-induced-dissociation (CID) gas cell to break up molecular contaminant ions that are generated during the stopping process. A prototype CID gas cell was constructed and tested with beams from an offline ion source, validating the concept with the successful demonstration of breaking of molecular ions. The outcome of this research enabled the formulation a conceptual design for a next-generation linear gas stopping device specifically tailored for FRIB. This device is envisioned to deliver rare-isotope-ion beams at a rate of 10 8 particles per second or higher, accompanied by advancements in purity. Methods employed: This project leverages advancements in technologies initially designed for the Advanced Cryogenic Gas Stopper (ACGS), the current state-of-the-art linear gas stopper, through the use of new simulations and beam purification via collision-induced-dissociation. The methods include: 1. Development of a prototype low-energy, low-pressure CID gas stopper. This prototype features a thin, approximately 20 nm, Si 3 N 4 entrance window designed for dissociating stable and rare-isotope molecular ions. The goal is to enhance beam purification and overall efficiency. 2. Creation of Particle-in-cell (PIC) simulations to assess the advantages of multi-layer RF carpets and multi-point extraction for ion transport efficiency. These simulations rely on the 3DCylPIC package, specifically designed for studying devices of this nature. The goal is to quantify and mitigate ion transport losses due to space charge generated in the stopping process of large numbers of ions. 3. Perform ion transport simulations across an RF carpet using an 8-phase travelling wave and evaluate its performance. Compared to the 4-phase RF carpets used in ACGS, the 8-phase carpets will double the wavelength of the generated traveling wave allowing for larger maximum RF amplitudes that could result in improved ion transport efficiency for high-intensity incoming beams when large space charge fields are present. Impact: Tackling the primary challenges associated with transforming high-energy projectile fragment beams into low-energy beams—specifically, addressing efficiency and purity—holds significant promise for advancing FRIB science. This advancement will particularly impact precision mass measurements, laser spectroscopy of short-lived nuclei, and studies in astrophysics and nuclear reactions using reaccelerated beams. These domains play a crucial role in addressing key questions outlined in the 2023 NSAC long-range plan, spanning nuclear structure, nuclear astrophysics, and fundamental symmetries. Additionally, they contribute to addressing 10 out of the 17 benchmarks identified by the NRC RISAC. The development of a next-generation gas stopping device capable of delivering low-energy, rare-isotope beams at a rate of 10 8 particles per second, or more, with high purity holds the potential to unlock experiments that would otherwise be unfeasible. Furthermore, it is expected to reduce the time required for experiments at FRIB, thereby maximizing scientific output. The research and development activities performed as part of this project bolstered essential competencies at FRIB in beam physics and ion source technologies, provided valuable training opportunities for junior scientists.

43 PARTICLE ACCELERATORS↗

Automating Bayesian inference and design to quantify acoustic particle levitation

Self-propulsion of micro- and nanoparticles powered by ultrasound provides an attractive strategy for the remote manipulation of colloidal matter using biocompatible energy inputs. Quantitative understanding of particle motion and its dependence on size, shape, and composition requires accurate characterization of the acoustic field, which depends sensitively on the experimental setup. Here, we show how automated experiments based on Bayesian inference and design can accurately and efficiently characterize the acoustic field within resonant chambers used to propel acoustic nanomotors. Repeated cycles of observation, inference, and design (OID) are guided by a physical model that describes the rate at which levitating particles approach the nodal plane. Using video microscopy, we observe the relaxation of tracer particles to this plane following the application of the acoustic field. We use sequential Monte Carlo methods to infer model parameters such as the amplitude and frequency of the resonant chamber while accounting for particle-level measurement noise and population-level heterogeneity in the field. Guided by simulated outcomes, we select the optimal design for the next experiment as to maximize the information gain in the relevant parameters. We show how this iterative process serves to discriminate between competing hypotheses and efficiently converges to accurate parameter estimates using only few automated experiments. We discuss the need for model criticism to ensure the validity of the guiding model throughout automated cycles of observation, inference, and design. Furthermore, this work demonstrates how Bayesian methods can learn the parameters of nonlinear, hierarchical models used to describe video microscopy data of active colloids.

42 ENGINEERING↗

Investigating late-stage particle production in pp collisions with balance functions

Balance functions have been regarded in the past as a method of investigating the late-stage hadronization found in the presence of a strongly-coupled medium. They are also used to constrain mechanisms of particle production in large and small collision systems. Charge balance functions for inclusive and identified particle pairs are reported as a function of charged particle multiplicity in proton–proton collisions simulated with the PYTHIA8 and the EPOS4 models. The charge balance functions of inclusive, pion, kaon, and proton pairs exhibit amplitudes and shapes that depend on particle species and differ significantly in the two models due to the different particle production mechanisms implemented in PYTHIA and EPOS. The shapes and amplitudes also evolve with multiplicity in both models. In addition, the evolution of the longitudinal rms width and that of balance functions integrals with multiplicity (and average transverse momentum) feature significant differences in the two models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Formation and Detriments of Residual Alkaline Compounds on High-Nickel Layered Oxide Cathodes

High-nickel layered oxides LiNi x M 1-x O 2 (x ≥ 0.9) have emerged as promising cathode materials for automotive batteries due to their high energy density and lower cost. However, the formation and accumulation of surface alkaline compounds during storage hinder their mass production and commercialization. Here, in this study, a validated chemical method is employed to deconvolute and quantify the evolution of each residual lithium compound in four representative cathodes during ambient-air storage, viz., LiNiO 2 (LNO), LiNi 0.95 Co 0.05 O 2 (NC), LiNi 0.95 Mn 0.05 O 2 (NM), and LiNi 0.95 Al 0.05 O 2 (NA). Furthermore, the activation energy of the reaction between water and the cathode is determined by measuring the leached LiOH concentration at various temperatures. While residual lithium and time-of-flight secondary-ion mass spectrometry measurements collectively reveal that the air stability overall follows the trend of NM > NA ≈ NC > LNO, the aged NM exhibits the highest charge-transfer resistance and the worst electrochemical performance among the cathodes. In situ, X-ray diffraction and scanning transmission electron microscopy unveil that the aged NM is plagued by a large area of resistive spinel-like M 3–x Li x O 4 phases, leading to aggravated particle reaction heterogeneity. Finally, a one-step recalcination method is demonstrated effective in fully restoring the degraded cathodes. This work provides insights into overcoming air sensitivity issues of high-Ni cathodes.

(S)TEM↗

Neutrino Transport with Monte Carlo Method. I. Toward Fully Consistent Implementation of Nucleon Recoils in Core-collapse Supernova Simulations

The small energy exchange via nucleon recoils in neutrino–nucleon scattering is now supposed to be one of the important factors for successful explosions of core-collapse supernovae (CCSNe), as they can change neutrino spectra through accumulation of a large number of scatterings. In deterministic methods employed for neutrino transport in CCSN simulations, we normally cannot afford to deploy a large enough number of energy bins needed to resolve this small energy exchange, and subgrid techniques are employed one way or another. In this paper, we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by ~15% for heavy-lepton neutrinos and much smaller amounts for other types of neutrinos in a typical postbounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the deterministic methods. In order to mimic the latter calculation, we redistribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected, depending on the energy resolution. We also discuss the possible implications of our results for the deterministic methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Determining g A / g V with High-Resolution Spectral Measurements Using a LiInSe 2 Bolometer

Neutrinoless double beta decay (0νββ) processes sample a wide range of intermediate forbidden nuclear transitions, which may be impacted by quenching of the axial vector coupling constant (g A /g V ), the uncertainty of which plays a pivotal role in determining the sensitivity reach of 0νββ experiments. Here we present measurements performed on a high-resolution LiInSe 2 bolometer in a “source = detector” configuration to measure the spectral shape of the fourfold forbidden β decay of 115 In . The value of g A /g V is determined by comparing the spectral shape of theoretical predictions to the experimental β spectrum taking into account various simulated background components as well as a variety of detector effects. We find evidence of quenching of g A /g V at >5 σ with a model-dependent quenching factor of 0.655 ± 0.002 as compared to the free-nucleon value for the interacting shell model. We also measured the 115 In half-life to be [5.18 ± 0.06(stat) $^{+0.005}_{- 0.015}$(sys)] x 10 14 yr within the interacting shell model framework. This Letter demonstrates the power of the bolometeric technique to perform precision nuclear physics single-β decay measurements, which along with improved nuclear modeling can help reduce the uncertainties in the calculation of several decay nuclear matrix elements including those used in 0νββ sensitivity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ultra-Microtome for the Preparation of TEM Samples from Battery Cathodes

An ultra-microtome method is developed for the preparation of transmission electron microscopy (TEM) specimens from the cathodes of Li-ion battery. Results indicate that the samples prepared with the ultra-microtome method yield atomic resolution imaging from the core area of particles of hundreds to thousands of nanometers. Analysis of the TEM samples prepared with ultra-microtome indicates that no observable structural damage, such as microcracking, is introduced during the cutting process, and the cathodes cut by ultra-microtome exhibit unchanged micro-morphologies compared with the uncut ones. These results demonstrate that ultra-microtome is a promising method for the preparation of TEM specimens and revealing in-depth microstructures of fragile materials such as battery cathodes.

transmission electron microscopy, Sample Preparati↗

Dose Coefficients for Internal Dose Assessments for Exposure to Radioactive Fallout

This paper presents values as well as the bases for calculating internal dose coefficients suitable for estimating organ doses from the exposure to radioactive fallout that could result from the detonation of a nuclear fission device. The 34 radionuclides discussed are the same as those given in a priority list of radionuclides for fallout dose assessments presented in a companion overview paper. The radionuclides discussed are those that are believed to account for a preponderance of the organ doses that might be received by intake by persons of all ages (including in utero and via breast feeding for infants) following exposure to radioactive fallout. The presented dose coefficients for ingestion account for age and include modifications for variations in solubility with distance as discussed previously in the literature, and those for inhalation similarly account for age, solubility, and particle sizes that would be relevant at various distances of exposure as discussed in a companion paper on ingestion dose methods. The proposed modifications peculiar to radioactive fallout account for systematic changes in solubility and particle sizes with distance from the site of detonation, termed here as the region of “local fallout” and the region “beyond local fallout.” Brief definitions of these regions are provided here with more detailed discussion in a companion paper on estimating deposition of fallout radionuclides. This paper provides the dose coefficients for ingestion and inhalation (for particle sizes of 1 μm, 5 μm, 10 μm, and 20 μm) for the region “local fallout.” These dose coefficients for “local fallout” are specific for particles formed in a nuclear explosion that can be large and have radionuclides, particularly the more refractory ones, distributed throughout the volume where the radionuclide has reduced solubility. The dose coefficients for the region “beyond local fallout” are assumed to be the ones published by the International Commission on Radiological Protection (ICRP) in 1995. Comparisons of the presented dose coefficients are made with values published by the ICRP.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Symplectic Particle Tracking in a Thick Nonlinear McMillan Lens for the Fermilab Integrable Optics Test Accelerator (IOTA)

The McMillan system is a novel method to increase the tune spread of a beam without decreasing its dynamic aperture due to the systems integrability. While the ideal system is based on an infinitely thin kick, the physical design requires a thick electron lens, including a solenoid. Particle transport through the lens is difficult to simulate due to the nature of the force on the circulating beam. This paper demonstrates accurate simulation of a thick McMillan lens in a solenoid using symplectic integrators derived from Yoshidas method.

43 PARTICLE ACCELERATORS↗