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 631 records · Page 35

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗

Sparse Convolutional Neural Networks for particle classification in ProtoDUNE-SP events

Deep Learning (DL) methods and Computer Vision are becoming important tools for event reconstruction in particle physics detectors. In this work, we report on the use of submanifold sparse convolutional neural networks (SparseNets) for the classification of track and shower hits from a DUNE prototype liquid-argon detector at CERN (ProtoDUNE-SP). By taking advantage of the three-dimensional nature of the problem we use a set of nine input features to classify sparse and locally dense hits associated to track or shower particles. The SparseNet has been trained on a test sample and shows promising results: efficiencies and purities greater than 90%. This has also been achieved with a considerable speedup and substantially less resource utilization with respect to other DL networks such as graph neural networks. This method offers great scalability advantages for future large neutrino detectors such as the planned DUNE experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of scattering and absorption properties of submillimeter bauxite and silica particles

Submillimeter solid particles have been considered as thermal storage media for concentrated solar power applications. Knowledge of the scattering and absorption properties of individual particles is crucial for modeling the radiative heat transfer of the particle bed. In this work, a laser scatterometer is used to measure the single-particle scattering properties at a wavelength of 635 nm by using two configurations: (1) a falling particle curtain and (2) a taped particle layer. Because a one-particle nominal thickness is formed with area fractions of 5-55% depending on the configuration, multiple scattering is minimized and hence the single scattering phase function, averaged over all illuminated particles, is directly measured. Bauxite-based ceramic particles that are strongly absorbing in the solar spectrum and silica particles that are nonabsorbing in the visible and near-infrared are investigated. The directional-hemispherical reflectance and transmittance of the taped particles are also measured to deduce the forward and backward scattering efficiency factors and the absorption efficiency factors. Only weak wavelength dependence is observed in the measured region from 380 nm to 1020 nm. Furthermore, the scattering phase functions of all bauxite-based particles with varying sizes and compositions are very similar and can be fitted to a Henyey-Greenstein phase function with an asymmetry factor g = -0.20. For the silica particles, forward scattering dominates and g = 0.45 yields the best fit. A Monte Carlo method is developed to model the particle scattering characteristics, and reasonable agreements between the modeling and experimental results are observed by introducing a specularity parameter.

42 ENGINEERING↗

Particle-based thermal energy storage systems

Methods and devices for long-duration electricity storage using low-cost thermal energy storage and high-efficiency power cycle, are disclosed. In some embodiments it has the potential for superior long-duration, low-cost energy storage.

Ma, Zhiwen↗

Particle-based thermal energy storage systems

Methods and devices for long-duration electricity storage using low-cost thermal energy storage and high-efficiency power cycle, are disclosed. In some embodiments it has the potential for superior long-duration, low-cost energy storage.

Ma, Zhiwen↗

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of hydrothermal synthesis conditions on up-conversion luminescence intensity of β-NaYF{sub 4} : Er{sup 3+}, Yb{sup 3+} submicron particles

The differences in the luminescence intensities of up-conversion β-NaYF{sub 4} : Er{sup 3+}, Yb{sup 3+} particles synthesised by the hydrothermal method under various synthesis conditions are studied. The results of the study lead to the conclusion that in order to achieve the maximum luminescence intensity in such particles, it is necessary to use ammonium fluoride and a medium with pH = 3. In this case, the length of the particles increases, up to the formation of rod-shaped particles. Based on the data on the size of the coherent scattering region and on microstresses, we can assume that the particles are polycrystals. At the same time, limiting the size of the coherent scattering region is possible due to the defective structure. When the nanoparticles are synthesised in a medium with pH = 3, hydrolysed regions containing OH groups are formed on the crystallite surface. The presence of these groups does not affect the intensity of up-conversion luminescence of submicron-size particles. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds

Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.

42 ENGINEERING↗

Particle hit clustering and identification using point set transformers in liquid argon time projection chambers

Liquid argon time projection chambers are often used in neutrino physics and dark-matter searches because of their high spatial resolution. The images generated by these detectors are extremely sparse, as the energy values detected by most of the detector are equal to 0, meaning that despite their high resolution, most of the detector is unused in a particular interaction. Instead of representing all of the empty detections, the interaction is usually stored as a sparse matrix, a list of detection locations paired with their energy values. Traditional machine learning methods that have been applied to particle reconstruction such as convolutional neural networks (CNNs), however, cannot operate over data stored in this way and therefore must have the matrix fully instantiated as a dense matrix. Operating on dense matrices requires a lot of memory and computation time, in contrast to directly operating on the sparse matrix. We propose a machine learning model using a point set neural network that operates over a sparse matrix, greatly improving both processing speed and accuracy over methods that instantiate the dense matrix, as well as over other methods that operate over sparse matrices. Compared to competing state-of-the-art methods, our method improves classification performance by 14%, segmentation performance by more than 22%, while taking 80% less time and using 66% less memory. Compared to state-of-the-art CNN methods, our method improves classification performance by more than 86%, segmentation performance by more than 71%, while reducing runtime by 91% and reducing memory usage by 61%.

calibration and fitting methods↗

Correlation of optical properties with particle size, morphology, and polymorph of fine- and nano-particle formulations of titanium dioxide powders

Titanium dioxide (TiO 2 ) particulates are known to exhibit different visible and infrared optical properties compared to the bulk material, showing strong dependence on particle size, crystal structure, and morphology. In this study, the optical properties, sizes, and morphologies of TiO 2 particles from two different sources (nano and fine powders) having a) nominally different particle sizes and b) various crystal polymorph mixture fractions are compared using a combination of single particle mass spectrometry, optical spectroscopies, and aerosol characterization methods. The nano sample was found to be largely particles of the anatase polymorph (88% by mass), while the fine sample was found to consist largely of rutile particles (95% by mass). Two distinct particle morphologies (fractal and compact) were found in each powder sample and could be identified and separated in-situ based on particle aerodynamic properties. The attenuation of near-infrared, visible and ultraviolet light by TiO 2 particles shows strong dependence on particle morphology. Furthermore, while the fine particles were found to have larger near-infrared (675–800 nm) extinction coefficients by mass than the nanoparticles, the reverse was true in the ultraviolet and visible regions (370–675 nm). However, for polydisperse particles with different sizes and shapes, the optical behaviors are not straightforward to directly correlate to a combination of physical parameters.

Lockwood, Schuyler P. [Pacific Northwest National ↗

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗

An implicit particle code with exact energy and charge conservation for studies of dense plasmas in axisymmetric geometries

A collisional particle code based on implicit energy- and charge-conserving methods in axisymmetric geometries is presented. A new particle pusher for axisymmetric systems is introduced that is compatible with exact energy and charge conservation and yields improved accuracy compared to other methods. How to appropriately treat all aspects of the algorithm near the r = 0 axis of symmetry is described in detail. Here, the axisymmetric model is verified by simulating the free expansion of a plasma sphere in 2D cylindrical and 1D spherical geometries. The algorithm's ability to study the dynamic compression of a dense plasma is illustrated by simulating the dynamic Z-pinch in 1D cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cryogenic Correlative Single-Particle Photoluminescence Spectroscopy and Electron Tomography for Investigation of Nanomaterials

Cryogenic single-particle photoluminescence (PL) spectroscopy has been used with great success to directly observe the heterogeneous photophysical states present in a population of luminescent particles. Cryogenic electron tomography provides complementary nanometer scale structural information to PL spectroscopy, but the two techniques have not been correlated due to technical challenges. Here, we present a method for correlating single-particle information from these two powerful microscopy modalities. We simultaneously observe PL brightness, emission spectrum, and in-plane excitation dipole orientation of CdSSe/ZnS quantum dots suspended in vitreous ice. Stable and fluctuating emitters were observed, as well as a surprising splitting of the PL spectrum into two bands with an average energy separation of 80 meV. In some cases, the onset of the splitting corresponded to changes in the in-plane excitation dipole orientation. These dynamics were assigned to structures of individual quantum dots and the excitation dipoles were visualized in the context of structural features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗