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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Apparatus for continuous, fast, and precise measurements of position and velocity of a small spherical particle

The position and velocity of optically levitated glass spheres (radii 10-20 microns) movng in a gas are measured accurately, rapidly, and continuously using a high-speed rotating polygon mirror. The experimental technique developed here has repeatable position accuracies better than 20 microns. Each measurement takes less than 1 microsec and can be repeated every 100 microsec. The position of the levitated glass spheres can be manipulated accurately by modulating the laser power with an acoustic optic modulator. The technique provides a fast and accurate method to study general particle dynamics in a fluid.

Venkataraman, T. S.↗
NASA NTRS2011

Spectral Analyses and Radiation Exposures from Several Ground-Level Enhancement (GLE) Solar Proton Events: A Comparison of Methodologies

Several methods for analyzing the particle spectra from extremely large solar proton events, called Ground-Level Enhancements (GLEs), have been developed and utilized by the scientific community to describe the solar proton energy spectra and have been further applied to ascertain the radiation exposures to humans and radio-sensitive systems, namely electronics. In this paper 12 GLEs dating back to 1956 are discussed, and the three methods for describing the solar proton energy spectra are reviewed. The three spectral fitting methodologies are EXP [an exponential in proton rigidity (R)], WEIB [Weibull fit: an exponential in proton energy], and the Band function (BAND) [a double power law in proton rigidity]. The EXP and WEIB methods use low energy (MeV) GLE solar proton data and make extrapolations out to approx.1 GeV. On the other hand, the BAND method utilizes low- and medium-energy satellite solar proton data combined with high-energy solar proton data deduced from high-latitude neutron monitoring stations. Thus, the BAND method completely describes the entire proton energy spectrum based on actual solar proton observations out to ~10 GeV. Using the differential spectra produced from each of the 12 selected GLEs for each of the three methods, radiation exposures are presented and discussed in detail. These radiation exposures are then compared with the current 30-day and annual crew exposure limits and the radiation effects to electronics.

Atwell, William↗
DOE OSTI2020

Singularity variables for missing energy event kinematics

We discuss singularity variables which are properly suited for analyzing the kinematics of events with missing transverse energy at the LHC. We consider six of the simplest event topologies encountered in studies of leptonic W-bosons and top quarks, as well as in SUSY-like searches for new physics with dark matter particles. In each case, we illustrate the general prescription for finding the relevant singularity variable, which in turn helps delineate the visible parameter subspace on which the singularities are located. Our results can be used in two different ways — first, as a guide for targeting the signal-rich regions of parameter space during the stage of discovery, and second, as a sensitive focus point method for measuring the particle mass spectrum after the initial discovery.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2020

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The ABACUS cosmological N -body code

We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct 1/r 2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10 -5 . We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Learning-based approach to plasticity in athermal sheared amorphous packings: Improving softness

The plasticity of amorphous solids undergoing shear is characterized by quasi-localized rearrangements of particles. While many models of plasticity exist, the precise relationship between the plastic dynamics and the structure of a particle’s local environment remains an open question. Previously, machine learning was used to identify a structural predictor of rearrangements called “softness.” Although softness has been shown to predict which particles will rearrange with high accuracy, the method can be difficult to implement in experiments where data are limited and the combinations of descriptors it identifies are often difficult to interpret physically. Here, we address both of these weaknesses, presenting two major improvements to the standard softness method. First, we present a natural representation of each particle’s observed mobility, allowing for the use of statistical models that are both simpler and provide greater accuracy in limited datasets. Second, we employ persistent homology as a systematic means of identifying simple, topologically informed, structural quantities that are easy to interpret and measure experimentally. We test our methods on two-dimensional athermal packings of soft spheres under quasi-static shear. We find that the same structural information that predicts small variations in the response is also predictive of where plastic events will localize. We also find that an excellent accuracy is achieved in athermal sheared packings using simply a particle’s species and the number of nearest neighbor contacts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

The binary Monte Carlo (MC) collision algorithm is a standard and robust method to include binary Coulomb collision effects in particle-in-cell (PIC) simulations of plasmas. Here we show that the coupling between PIC and MC algorithms can give rise to (nonphysical) numerical heating of the system that significantly exceeds that observed when these algorithms operate independently. We argue that this deleterious effect results from an inconsistency between the particle motion associated with MC collisions and the work performed by the collective electromagnetic field on the PIC grid. This inconsistency manifests as the (artificial) stochastic production of electromagnetic energy, which ultimately heats the plasma particles. Here, the MC-induced numerical heating can significantly impact the evolution of the simulated system for long simulation times (≳10 3 collision periods, for typical numerical parameters). We describe the source of the MC-induced numerical heating analytically and discuss strategies to minimize it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Deep neural network uncertainty quantification for LArTPC reconstruction

We evaluate uncertainty quantification (UQ) methods for deep learning applied to liquid argon time projection chamber (LArTPC) physics analysis tasks. As deep learning applications enter widespread usage among physics data analysis, neural networks with reliable estimates of prediction uncertainty and robust performance against overconfidence and out-of-distribution (OOD) samples are critical for their full deployment in analyzing experimental data. While numerous UQ methods have been tested on simple datasets, performance evaluations for more complex tasks and datasets are scarce. Here we assess the application of selected deep learning UQ methods on the task of particle classification using the PiLArNet monte carlo 3D LArTPC point cloud dataset. We observe that UQ methods not only allow for better rejection of prediction mistakes and OOD detection, but also generally achieve higher overall accuracy across different task settings. We assess the precision of uncertainty quantification using different evaluation metrics, such as distributional separation of prediction entropy across correctly and incorrectly identified samples, receiver operating characteristic curves (ROCs), and expected calibration error from observed empirical accuracy. We conclude that ensembling methods can obtain well calibrated classification probabilities and generally perform better than other existing methods in deep learning UQ literature.

47 OTHER INSTRUMENTATION↗
NASA NTRS1983

Settling of fixed erythrocyte suspension droplets

It is pointed out that when particles behave collectively rather than individually, the fractionation of micron-size particles on the basis of size, density, and surface characteristics by centrifugation and electrophoresis is hindered. The formation and sedimentation of droplets containing particles represent an extreme example of collective behavior and pose a major problem for these separation methods when large quantities of particles need to be fractionated. Experiments are described that measure droplet sizes and settling rates for a variety of particles and droplets. Expressions relating the particle concentration in a drop to measurable quantities of the fluids and particles are developed. The number of particles in each droplet is then estimated, together with the effective droplet density. Red blood cells from different animals fixed in glutaraldehyde provide model particle groups.

Omenyi, S. N.↗
DOE OSTI2021

Self-Stabilized LiNi 0.8 Mn 0.1 Co 0.1 O 2 in thiophosphate-based all-solid-state batteries through extra LiOH

Nickle-rich LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC 811) cathode material exhibits engaging properties in high energy density and low cost, making it great potential for the next generation high-energy all-solid-state lithium batteries (ASSLBs). However, NMC 811 suffers from severe surface electrochemical, chemical, and voltage incompatibility towards solid-state electrolytes (SSE), especially thiophosphate-based electrolytes like Li 6 PS 5 Cl. Although diverse coating methods have been made to overcome this issue, they are typically cumbersome and expensive. A coating strategy that satisfied all the requirements of cost-efficiency, stability, uniformity, scalability, and easy-achieving is still challenging. Here, we developed a LiOH-based surface stabilization strategy that provides a ~10 nm stable permeable layer on NMC 811. After one-step sintering of NMC 811 precursor mixed with LiOH, which is commonly used for NMC 811 lithiation process, excessive LiOH simultaneously distributes on NMC 811 particles. Unlike other reported methods, this coating method can be easily controlled and fabricated without additional complicated processes. By simply controlling the thickness of LiOH layer, which protects the Li 6 PS 5 Cl solid electrolyte materials from being oxidized, optimized cycling stability can be obtained for 600 cycles with capacity of 130 mAh g – 1 on average at a wide electrochemical window of 2.50–4.20 V (vs. Li-In).

25 ENERGY STORAGE↗
DOE OSTI2023

Pulse shape measurements for neutron/gamma discrimination using the TOFPET2 ASIC

Many highly pixelated organic scintillator detection systems would benefit from independent readout of each scintillator pixel. Recent advances in Silicon Photomultiplier (SiPM) technology makes this goal feasible, however the data acquisition from potentially hundreds or thousands of channels requires a low-cost and compact solution. For pixelated neutron detection with organic scintillators, the capability to distinguish between neutron and gamma interactions using Pulse Shape Discrimination (PSD) is required along with pulse charge and time of arrival. The TOFPET2 ASIC from PETsys Electronics is a 64-channel readout chip providing pulse time and charge integration measurements from SiPMs, and is specifically designed for time-of-flight positron-emission tomography. Using an 8 × 8 array of 6 mm × 6 mm J-series SiPMs from SensL/OnSemi (ArrayJ-60035-64P-PCB), we have studied the energy and PSD performance of the TOFPET2 ASIC using a 4 × 4 array of 6 mm × 6 mm × 30 mm trans -Stilbene crystals from Inrad Optics and a custom SiPM routing board from PETsys Electronics. In conclusion, using a time-over-threshold method, we measure a maximum PSD figure-of-merit of approximately 1.2 at 478 keV (the Compton edge of 662 keV) for a J-series SiPM operating at an over-voltage of 3V.

47 OTHER INSTRUMENTATION↗
NASA NTRS2001

Monte Carlo track structure for radiation biology and space applications

Over the past two decades event by event Monte Carlo track structure codes have increasingly been used for biophysical modelling and radiotherapy. Advent of these codes has helped to shed light on many aspects of microdosimetry and mechanism of damage by ionising radiation in the cell. These codes have continuously been modified to include new improved cross sections and computational techniques. This paper provides a summary of input data for ionizations, excitations and elastic scattering cross sections for event by event Monte Carlo track structure simulations for electrons and ions in the form of parametric equations, which makes it easy to reproduce the data. Stopping power and radial distribution of dose are presented for ions and compared with experimental data. A model is described for simulation of full slowing down of proton tracks in water in the range 1 keV to 1 MeV. Modelling and calculations are presented for the response of a TEPC proportional counter irradiated with 5 MeV alpha-particles. Distributions are presented for the wall and wall-less counters. Data shows contribution of indirect effects to the lineal energy distribution for the wall counters responses even at such a low ion energy.

NASA Center JSC↗
DOE OSTI2022

Particulate Fuel Modeling of MC 2 -3 using Iterative Local Spatial Self-shielding Method

We report a new spatial self-shielding method for particulate fuels has been developed based on disadvantage factors and implemented in the MC 2 -3 code. This method named the iterative local spatial self-shielding (ILSS) method considers the shadowing effect of randomly distributed particles on spatial self-shielding in particles through a homogenized composition region added outside the particle of interest at the center. The self-shielded cross sections of the central particle are determined iteratively since they are used in determining the cross sections of the homogenized composition region. The ILSS method was verified for infinite stochastic medium problems of single and multiple types of particles, VHTR unit cell problems, and HTTR assembly problems. The verification test results show that the ILSS method accurately predicts the stochastic particle shadowing effect and reaction rates in particles, whereas the regular array model and the stochastic collision probability method underpredict the particle shadowing effect and overestimate reaction rates in particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Quality Control Methods for Measurement of UCO Kernel Composition and SiC Microstructure

Quality control (QC) is critically important to tristructural-isotropic (TRISO) particle fuels owing to the complexity of and reliance on the fuel form to contain fission products during irradiation. Characterization methods for particle fuel QC have decades of history and have continued to develop as new insights into fuel performance inform revised fuel specifications and as advances in underlying technologies expand the possibilities of what may be characterized. Two relatively new methods for characterization of TRISO fuels have been published in open literature: optical microscopy image analysis for mixed uranium carbide/uranium oxide (UCO) kernel composition analysis and automated grain boundary detection in backscattered electron (BSE) images of the silicon carbide (SiC) layer in TRISO particles for grain size characterization. Suggestions and guidelines for the application of these methods to TRISO fuel qualification are provided in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Observation of Deep Convective Cloud-Top Height and Vertical Temperature Structure of Hurricane Using Hyperspectral Infrared Sounder and its Single-Field-View Retrieval Products

Hurricanes, severe tropical cyclones (TC), or typhoons are significant natural disasters that often result in substantial loss of life and property damage. Numerous studies have indicated that changes in TC intensity are closely linked to deep convective clouds (DCC), with stronger TCs typically exhibiting higher cloud top heights (CTH) compared to weaker TCs. The CTH can help determine if a tropical depression is at the onset of rapid intensification based on case studies. Therefore, accurate determination of TC CTH will be greatly helpful for monitoring TC development and studying TC dynamics. One traditional and most common method to derive CHT from satellite observations is using the thermal brightness temperature in atmospheric channels to match the sounding temperature profile. However, it was found that thermally derived CTH has a lower bias of approximately 1 km, and this bias tends to worsen for the tallest clouds. A new method using the hyperspectral infrared sounder will be presented. From the measurements of Cross-track Infrared Sounder (CrIS) on S-NPP and J-1, along with radiative transfer simulations, we identified the inverted-V spectral feature in the ozone (O3) band (near 9.6 μm) corresponding to high clouds. The depth of the inverted-V can be used to estimate the CTH. Since the depth is computed using the peak absorption O3 channel and the nearby most transparent O3 channel in this O3 band, the uncertainties associated with cloud emissivity and scattering by cloud particles in the traditional method can be ignored. From several hurricane case studies, we found that the CHT derived using this method can accurately capture the structure of the cloud tops in the eyewall, spiral rainbands, and surrounding regions. For example, Hurricane Dorian on September 2, 2019, showed a nicely outward-sloping and circular shape eye cloud in the early morning, but the circular shape of the eyewall cloud became distorted in the afternoon. For various hurricanes we examined, the distribution of CHT for the eyewall clouds differed significantly. To better study the thermodynamic structure of hurricane clouds, this research will analyze the vertical temperature profiles from a new single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP), derived using CrIS and the Advanced Technology Microwave Sounder (ATMS) onboard SNPP and JPSS-1. SiFSAP has a spatial resolution of 15 km at nadir, which surpasses most global weather and climate models and other current operational sounding products. The combined use of ATMS and CrIS allows for retrievals near hurricane eyewalls and spiral rainbands. Wind fields from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) and ERA5 will be used to characterize transport, and comparisons between the model temperature and water vapor profiles with the corresponding SiFSAP products will also be provided.

SiFSAP↗
NASA NTRS2010

Effects of Ice Particle Size Vertical Inhomogeneity on the Passive Remote Sensing of Ice Clouds

The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.

Zhang, Zhibo↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗