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 343 records · Page 19

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗

A CUPID Li2100MoO4 scintillating bolometer tested in the CROSS underground facility

A scintillating bolometer based on a large cubic Li$_{2}$$^{100}$MoO$_4$ crystal (45 mm side) and a Ge wafer (scintillation detector) has been operated in the CROSS cryogenic facility at the Canfranc underground laboratory in Spain. The dual-readout detector is a prototype of the technology that will be used in the next-generation $0\nu2\beta$ experiment CUPID. The measurements were performed at 18 and 12 mK temperature in a pulse tube dilution refrigerator. This setup utilizes the same technology as the CUORE cryostat that will host CUPID and so represents an accurate estimation of the expected performance. The Li$_{2}$$^{100}$MoO$_4$ bolometer shows a high energy resolution of 6 keV FWHM at the 2615 keV $\gamma$ line. The detection of scintillation light for each event triggered by the Li$_{2}$$^{100}$MoO$_4$ bolometer allowed for a full separation ($\sim$8$\sigma$) between $\gamma$($\beta$) and $\alpha$ events above 2 MeV. The Li$_{2}$$^{100}$MoO$_4$ crystal also shows a high internal radiopurity with $^{228}$Th and $^{226}$Ra activities of less than 3 and 8 $\mu$Bq/kg, respectively. Taking also into account the advantage of a more compact and massive detector array, which can be made of cubic-shaped crystals (compared to the cylindrical ones), this test demonstrates the great potential of cubic Li$_{2}$$^{100}$MoO$_4$ scintillating bolometers for high-sensitivity searches for the $^{100}$Mo $0\nu2\beta$ decay in CROSS and CUPID projects.

47 OTHER INSTRUMENTATION↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PASSAT: particle accelerator helioScopes for Slim Axion-like-particle deTection

We propose a novel method to search for axion-like particles (ALPs) at particle accelerator experiments. ALPs produced at the target via the Primakoff effect subsequently enter a region with a magnetic field, where they are converted to photons that are then detected. Dubbed Particle Accelerator helioScopes for Slim Axion-like-particle deTection (PASSAT), our proposal uses the principle of the axion helioscope but replaces ALPs produced in the Sun with those produced in a target material. Since we rely on ALP-photon conversions , our proposal probes light (slim) ALPs that are otherwise inaccessible to laboratory-based experiments which rely on ALP decay, and complements astrophysical probes that are more model-dependent. As a first application, we reinterpret existing data from the NOMAD experiment in light of PASSAT, and constrain the parameter space for ALPs lighter than $$\sim 100~\mathrm{eV}$$ ~ 100 eV and ALP-photon coupling larger than $$\sim 10^{-4}~\mathrm{GeV}^{-1}$$ ~ 10 - 4 GeV - 1 . As benchmarks of feasible low-cost experiments improving over the NOMAD limits, we study the possibility of re-using the magnets of the CAST and the proposed BabyIAXO experiments and placing them at the proposed BDF facility at CERN, together with some new detectors. We find that these realizations of PASSAT allow for a direct probe of the parameter space for ALPs lighter than $$\sim 100~\mathrm{eV}$$ ~ 100 eV and ALP-photon coupling larger than $$\sim 4\times 10^{-6}~\mathrm{GeV}^{-1}$$ ~ 4 × 10 - 6 GeV - 1 , which are regions that have not been probed yet by experiments with laboratory-produced ALPs. In contrast to other proposals aiming at detecting single or two-photon only events in hadronic beam dump environments, that rely heavily on Monte Carlo simulations, the background in our proposal can be directly measured in-situ , its suppression optimized, and the irreducible background statistically subtracted. Sensitivity evaluations with other beams will be the subject of a future paper. The measurements suggested in this paper represent an additional physics case for the BDF at CERN beyond those already proposed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring the Regimes of Particle Behavior upon Impact via the Discrete Element Method

Discrete element method simulations are conducted to probe the various regimes of post-impact behavior of particles with solid surfaces. The impacting particles are described as spherical agglomerates consisting of smaller constituent (or primary) particles held together via surface adhesion. Under the influence of a wide range of impact velocities and particle surface energies, five distinct behavioral regimes—rebounding, vibration, fragmentation, pancaking, and shattering—are identified, and force transmission patterns are linked to post-impact behavior. In the rebounding regime, the coefficient of restitution decreases linearly as impact velocity increases and the particle agglomerate experiences compaction. In the fragmentation regime, rebound velocity generally decreases with increasing fragment size. The rebound velocity of fragments decreases with time except for the smallest fragments, which can increase in velocity due to collisions with other fragments of high velocity. Particle breakage in the pancaking regime does not follow common mechanistic models of breakage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compositions including nano-particles and a nano-structured support matrix and methods of preparation as reversible high capacity anodes in energy storage systems

The present invention relates to compositions including nano-particles and a nano-structured support matrix, methods of their preparation and applications thereof. The compositions of the present invention are particularly suitable for use as anode material for lithium-ion rechargeable batteries. The nano-structured support matrix can include nanotubes, nanowires, nanorods, and mixtures thereof. The composition can further include a substrate on which the nano-structured support matrix is formed. The substrate can include a current collector material.

Kumta, Prashant Nagesh↗

Titanium-group nano-whiskers and method of production

Disclosed herein are structures comprising a titanium, zirconium, or hafnium powder particle with titanium carbide, zirconium carbide, or hafnium carbide (respectively) nano-whiskers grown directly from and anchored to the powder particle. Also disclosed are methods for fabrication of such structures, involving heating the powder particles and exposing the particles to an organic gas.

Seals, Roland D.↗

Titanium-group nano-whiskers and method of production

Disclosed herein are structures comprising a titanium, zirconium, or hafnium powder particle with titanium carbide, zirconium carbide, or hafnium carbide (respectively) nano-whiskers grown directly from and anchored to the powder particle. Also disclosed are methods for fabrication of such structures, involving heating the powder particles and exposing the particles to an organic gas.

Menchhofer, Paul A.↗

Quantifying the impact of operating temperature on cracking in battery electrodes, using super-resolution of microscopy images and stereology

There are numerous factors that can have an impact on the degradation behavior of batteries, such as the number of recharge cycles or the charge rate. Here, we investigate the influence of operating temperature on the structural degradation of the microstructure in lithium-ion positive electrodes. For that purpose, the microstructure is characterized for cathodes which have been cycled for 200 cycles under 6C (10-minute) charging at different operating temperatures, namely, 20°C, 30°C, 40°C, and 50°C. For each operating condition scanning electron microscopy (SEM) images of cross-sectioned Li x Ni 0.5 Mn 0.3 Co 0.2 O 2 (NMC532) electrodes have been analyzed, to determine structural descriptors such as global particle porosity, crack size/length/width distribution, and porosity and specific surface area distribution of individual particles. Additionally, a stereological method has been deployed to investigate the local particle porosity as a function of distance to the particle center. Results show that particle porosity increases with increasing cycling temperature. Particle porosity is greatest at the particle center and decreases along the particle radius to the exterior. Particle surface area is similar across the four cycling-temperature aging conditions.

25 ENERGY STORAGE↗

Timing Performance of the CMS High Granularity Calorimeter Prototype

This paper describes the experience with the calibration,reconstruction and evaluation of the timing capabilities of the CMSHGCAL prototype in the beam tests in 2018. The calibrationprocedure includes multiple steps and corrections ranging from tensof nanoseconds to a few hundred picoseconds. The timing performanceis studied using signals from positron beam particles with energiesbetween 20 GeV and 300 GeV. The performance is studied as afunction of particle energy against an external timing reference aswell as standalone by comparing the two different halves of theprototype. The timing resolution is found to be 60 ps forsingle-channel measurements and better than 20 ps for full showersat the highest energies, setting excellent perspectives for theHGCAL calorimeter performance at the HL-LHC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Noble Gas Management: SBMOF 1 vs. NUCON Carbon

Pacific Northwest National Laboratory (PNNL), in collaboration with Flibe Energy, Inc., demonstrated the fabrication of metal organic frameworks (MOFs) into engineered beads to manage off-gases released from molten salt reactors. Among the many MOF options, SBMOF 1 was selected as an ideal sorbent to manage off-gases released from a molten salt reactor because its pore size is optimal. The SBMOF 1 was synthesized and fabricated into engineered particles using three different methods, with and without polymers. All the engineered particles were exposed to Co-60 radiation (1,000 kGy) to demonstrate their radiation stability. Single column breakthrough experiments were performed at different temperatures, activation conditions, flow rates, and off-gas compositions to optimize the Xe capacity at the breakthrough point and at saturation. The Xe loading in SBMOF 1 was compared with that of standard solid sorbent (Nucon carbon) under identical conditions. Based on the single column experiments, SBMOF 1 outperforms the Nucon carbon at all the temperatures and reduces the mass (to as little as 55%) of sorbent needed to achieve the same performance as carbon. The bed volume is also reduced, down to 48% of that of Nucon carbon, when the MOF is used. These results clearly demonstrate improvement in the amount of adsorbent needed and show a reduction in bed volume compared to traditional sorbent material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗

Progressive Tree-Based Compression of Large-Scale Particle Data

Scientific simulations and observations using particles have been creating large datasets that require effective and efficient data reduction to store, transfer, and analyze. However, current approaches either compress only small data well while being inefficient for large data, or handle large data but with insufficient compression. Toward effective and scalable compression/decompression of particle positions, we introduce new kinds of particle hierarchies and corresponding traversal orders that quickly reduce reconstruction error while being fast and low in memory footprint. Our solution to compression of large-scale particle data is a flexible block-based hierarchy that supports progressive, random-access, and error-driven decoding, where error estimation heuristics can be supplied by the user. For low-level node encoding, we introduce new schemes that effectively compress both uniform and densely structured particle distributions. Our proposed methods thus target all three phases of a tree-based particle compression pipeline, namely tree construction, tree traversal, and node encoding. In conclusion, the improved efficacy and flexibility of these methods over existing compressors are demonstrated through extensive experimentation, using a wide range of scientific particle datasets.

97 MATHEMATICS AND COMPUTING↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Image registration for in situ X-ray nano-imaging of a composite battery cathode with deformation

The structural and chemical evolution of battery electrodes at the nanoscale plays an important role in affecting the cell performance. Nano-resolution X-ray microscopy has been demonstrated as a powerful technique for characterizing the evolution of battery electrodes under operating conditions with sensitivity to their morphology, compositional distribution and redox heterogeneity. In real-world batteries, the electrode could deform upon battery operation, causing challenges for the image registration which is necessary for several experimental modalities, e.g. XANES imaging. To address this challenge, this work develops a deep-learning-based method for automatic particle identification and tracking. This approach was not only able to facilitate image registration with good robustness but also allowed quantification of the degree of sample deformation. The effectiveness of the method was first demonstrated using synthetic datasets with known ground truth. The method was then applied to an experimental dataset collected on an operating lithium battery cell, revealing a high degree of intra- and interparticle chemical complexity in operating batteries.

25 ENERGY STORAGE↗

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Furthermore, our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

97 MATHEMATICS AND COMPUTING↗