TESTING AND SIMULATIONS OF SPATIAL AND TEMPORAL TEMPERATURE VARIATIONS IN A PARTICLE-BASED THERMAL ENERGY STORAGE BIN.
Abstract not provided.
Sment, Jeremy Niko Iversen↗
SEARCH · Search NASA
Results for “test particle simulations”
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.
Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos
When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Generating multi-scale Li-ion battery cathode particles with radial grain architectures using stereological generative adversarial networks
Abstract Understanding structure-property relationships of Li-ion battery cathodes is crucial for optimizing rate-performance and cycle-life resilience. However, correlating the morphology of cathode particles, such as in LiNi0.8Mn0.1Co0.1O2 (NMC811), and their inner grain architecture with electrode performance is challenging, particularly, due to the significant length-scale difference between grain and particle sizes. Experimentally, it is not feasible to image such a high number of particles with full granular detail. A second challenge is that sufficiently high-resolution 3D imaging techniques remain expensive and are sparsely available at research institutions. Here, we present a stereological generative adversarial network-based model fitting approach to tackle this, that generates representative 3D information from 2D data, enabling characterization of materials in 3D using cost-effective 2D data. Once calibrated, this multi-scale model can rapidly generate virtual cathode particles that are statistically similar to experimental data, and thus is suitable for virtual characterization and materials testing through numerical simulations. A large dataset of simulated particles with inner grain architecture has been made publicly available.
25 ENERGY STORAGE↗
Particle‐In‐Cell Simulations of Starfish Prime
The Starfish Prime high altitude nuclear test created a transient diamagnetic cavity in the Earth's magnetic field above Johnston Island and launched an electromagnetic pulse (EMP) that was detected around the globe. For this work, we use the ion‐kinetic particle‐in‐cell code Topanga to simulate diamagnetic cavity evolution and the E3 EMP signal for Starfish Prime out to over a minute of physical time. The simulation domain has a longitudinal and latitudinal extent of 60° X 60° and a vertical extent of 2,000 km from the surface of the Earth. We compare our simulated results to magnetometer measurements taken in space and on the ground, finding good agreement in both cases. The diamagnetic cavity in the simulation forms in about a second, while the associated debris flux tube takes approximately 30 s to decay. The debris flux tube undergoes significant motion during this timeframe, rising upward. The measured E3 EMP signal on the ground consists of several components, all of which are present in our simulations. We discuss the physical origin of these components in relation to E3a (blast) and E3b (heave).
Belyaev, Mikhail A. [Lawrence Livermore National L↗
Selecting XFEL single-particle snapshots by geometric machine learning
A promising new route for structural biology is single-particle imaging with an X-ray Free-Electron Laser (XFEL). This method has the advantage that the samples do not require crystallization and can be examined at room temperature. However, high-resolution structures can only be obtained from a sufficiently large number of diffraction patterns of individual molecules, so-called single particles. Here, we present a method that allows for efficient identification of single particles in very large XFEL datasets, operates at low signal levels, and is tolerant to background. This method uses supervised Geometric Machine Learning (GML) to extract low-dimensional feature vectors from a training dataset, fuse test datasets into the feature space of training datasets, and separate the data into binary distributions of “single particles” and “non-single particles.” As a proof of principle, we tested simulated and experimental datasets of the Coliphage PR772 virus. We created a training dataset and classified three types of test datasets: First, a noise-free simulated test dataset, which gave near perfect separation. Second, simulated test datasets that were modified to reflect different levels of photon counts and background noise. These modified datasets were used to quantify the predictive limits of our approach. Third, an experimental dataset collected at the Stanford Linear Accelerator Center. The single-particle identification for this experimental dataset was compared with previously published results and it was found that GML covers a wide photon-count range, outperforming other single-particle identification methods. Moreover, a major advantage of GML is its ability to retrieve single particles in the presence of structural variability.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Simulation Studies on the Interactions of Electron Beam with Wastewater
The manufactured chemical pollutants, like 1,4 dioxane and PFAS (per- and polyfluroralkyl substances), found in the underground water and/or drinking water are challenging to be removed or biodegraded. Energetic electrons are capable of mediating and removing them. This paper utilizes FLUKA code to evaluate the beam-wastewater interaction effects with different energy, space and divergence distributions of the electron beam. With 8 MeV average energy, the electron beam exits from a 0.0127 cm thick titanium window, travels through a 4.3 cm distance air and a second 0.0127 cm thick stainless water container window with 2.43 cm radius, and finally is injected into the water area, where the volume of water is around 75 cubic cm. The distribution parameters of the electron beam are from the GPT (General Particle Tracer) simulations for UITF (Upgraded Injector Test Facility) in Jefferson lab. By varying the distributions, several measurements including the dose (or energy deposition) distribution, electron fluence, photon fluence are scored and compared. Taking the comparisons into consideration, this paper is aiming to find better electron beams for the wastewater irradiation.
Li, X.↗
Simulation Studies on the Interactions of Electron Beam with Wastewater
The manufactured chemical pollutants, like 1,4 dioxane and PFAS (per- and polyfluroralkyl substances), found in the underground water and/or drinking water are challenging to be removed or biodegraded. Energetic electrons are capable of mediating and removing them. This paper utilizes FLUKA code to evaluate the beam-wastewater interaction effects with different energy, space and divergence distributions of the electron beam. With 8 MeV average energy, the electron beam exits from a 0.0127 cm thick titanium window, travels through a 4.3 cm distance air and a second 0.0127 cm thick stainless water container window with 2.43 cm radius, and finally is injected into the water area, where the volume of water is around 75 cubic cm. The distribution parameters of the electron beam are from the GPT (General Particle Tracer) simulations for UITF (Upgraded Injector Test Facility) in Jefferson lab. By varying the distributions, several measurements including the dose (or energy deposition) distribution, electron fluence, photon fluence are scored and compared. Taking the comparisons into consideration, this paper is aiming to find better electron beams for the wastewater irradiation.
Li, X.↗
Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks
Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
High-Temperature Particle Heat Exchanger for sCO2 Power Cycles [Award 30342]
This report describes the design, development, and testing of a prototype 100 kWt particle-to-supercritical CO 2 (sCO2) heat exchanger. An analytic hierarchy process was implemented to compare and evaluate alternative heat-exchanger designs (fluidized bed, shell-and-plate moving packed bed, and shell-and-tube moving packed bed) that could meet the high pressure (≥ 20 MPa) and high temperature (≥ 700 °C) operational requirements associated with sCO2 power cycles. Cost, heat-transfer coefficient, structural reliability, manufacturability, parasitics and heat losses, scalability, compatibility, erosion and corrosion, transient operation, and inspection ease were considered in the evaluation. A 100 kWt shell-and-plate design was selected for construction and integration with Sandia’s falling particle receiver system that heats the particles using concentrated sunlight. Sandia worked with industry to design and construct the moving packed-bed shell-and-plate heat exchanger. Tests were performed to evaluate its performance using both electrical heating and concentrated sunlight to heat the particles. Overall heat transfer coefficients at off-design conditions (reduced operating temperatures and only three stainless steel banks in the counter-crossflow heat exchanger) were measured to be approximately ~25 - 70 W/m 2 -K, significantly lower than simulated values of >100 W/m 2 -K. Tests using the falling particle receiver to heat the particles with concentrated sunlight yielded overall heat transfer coefficients of ~35 – 80 W/m 2 -K with four banks (including a nickel-alloy bank above the three stainless steel banks). The overall heat transfer coefficient was observed to decrease with increasing particle inlet temperatures, which contrasted the results of simulations that showed an increase in heat transfer coefficient with temperature due to increased effective particle-bed thermal conductivity from radiation. The likely cause of the discrepancy was particle-flow maldistributions and funnel flow within the heat exchanger caused by internal ledges and cross-bracing, which could have been exacerbated by increased particle-wall friction at higher temperatures. Additional heat loss at higher temperatures may also contribute to a lower overall heat-transfer coefficient. Design challenges including pressure drop, particle and sCO2 flow maldistribution, and reduced heat transfer coefficient are discussed with approaches for mitigation in future designs. Lessons learned regarding instrumentation, performance characterization, and operation of particle components and sCO2 flow loops are also discussed. Finally, a 200 MWt commercial-scale shell-and-plate heat-exchanger design based on the concepts investigated in this report is proposed.
14 SOLAR ENERGY↗
Development of a Flotation Simulator That Can Predict Grade vs. Recovery Curves from Mineral Liberation Data
Flotation is a kinetic process designed to separate mineral fines by control of hydrophobicity. In general, both the recovery and selectivity of a flotation process increase with contact angle. It is difficult, however, to model flotation kinetics using the thermodynamic parameter. In the present work, we consider flotation as heterocoagulation so that the bubble-particle interactions occurring in the pulp phase of a flotation cell can be modeled using the hydrophobic force as a kinetic parameter that can represent the role of contact angles in flotation. Here, the pulp-phase model developed in this manner has been combined with a froth-phase model developed earlier by considering the effect of particles in film thinning kinetics and bubble coarsening (Park et al, 2018). The combined model has been used to develop a flotation simulator that can be used to construct grade vs. recovery curves from feed characteristics (e.g., mineral liberation and particle size) and various operating parameters. The simulator has been tested successfully against a set of bench-scale continuous flotation test results reported in the literature (dos Santos and Galery, 2018).
58 GEOSCIENCES↗
Post-irradiation Heating Tests of As-Irradiated AGR-3/4 TRISO Fuel Compacts
Four post-irradiation heating tests of fuel compacts from the U.S. Advanced Gas Reactor (AGR)-3/4 irradiation experiment were completed. In addition to tristructural isotropic (TRISO)-coated driver fuel, each compact contained designed-to-fail (DTF) particles with fuel kernels coated only in pyrocarbon so as to simulate exposed kernels. Tests at 1600/1700°C, 1400°C, and 1200°C were performed to measure fission product release as a function of time and temperature. Silver release was highest in the 1200°C test, supporting the observation that silver release rates are highest in the 1100–1300°C range. Compared to tests of AGR-1 compacts with no exposed kernels, the Cs-134 and Kr-85 releases were noticeably higher in AGR-3/4. The exposed kernels’ contributions to Eu and Sr release are inconclusive, due to the difficulty in distinguishing among the combined effects of higher irradiation temperatures in these particular AGR-3/4 compacts, the presence of the DTF particles, and the Fuel Accident Condition Simulator (FACS) test temperatures. These data can be used to make inferences about fission product retention in exposed kernels as a function of time and temperature.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
gRASPA
GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.
Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
When should PIC simulations be applied to atmospheric pressure plasmas? Impact of correlation heating
Abstract Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to artificial correlation heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w 2 / 3 , where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations
The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.
, Monte Carlo methods↗
2D radial-azimuthal particle-in-cell benchmark for E × B discharges
In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Long-term Dynamical Stability in the Outer Solar System. I. The Regular and Chaotic Evolution of the 34 Largest Trans-Neptunian Objects
We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly 100,000 test particles, as well as a statistical analysis of 200 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects “dwarf planets,” DPs, even if probably only the largest of them are true DPs). From the frequency analysis, we produced statistical diffusion maps for a wide region of the a–e phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs, making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories, including 17 stable, 11 unstable, and 6 resonant objects; we also found that, statistically, two objects from the sample will leave the trans-Neptunian region within the next gigayear, most likely being ejected from the solar system, but with a nonnegligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.
79 ASTRONOMY AND ASTROPHYSICS↗
Acceleration of Oxygen Ions In Dipolarization Events: 2. PSBL Distributions
Abstract This paper represents the second part of an investigation of the acceleration of energetic oxygen ions from encounters with a dipolarization front (DF), based on test particle tracing in the fields of an MHD simulation. In this paper, we focus on distributions in the plasma sheet boundary layer (PSBL). O + beams close to the plasma sheet boundary are found to be less pronounced and/or delayed against the H + beams. The reason is that these particles are accelerated by nonadiabatic motion in the duskward electric field such that O + ions gain the same amount of energy, but only 1/4 of the speed of protons. This causes a delay and larger equatorward displacement by the E × B drift. In contrast, the O + beams somewhat deeper inside the plasma sheet, where previously multiple proton beams were found, are accelerated at an earthward propagating DF just like H + , forming a field‐aligned beam at a similar speed as the lowest‐energy H + beam. We found that the source location depends on the adiabaticity of the orbit. For larger adiabaticity, the beam ions originate initially from the outer plasma sheet, but later from the opposite PSBL or lobe, but for low adiabaticity, sources are well inside the plasma sheet. The energy gained from a single encounter of a DF is comparable to the kinetic energy associated with the front speed. Assuming maximum speeds of 500–1,000 km/s, this yields a mass dependent acceleration of about 1–5 keV for protons and 20–80 keV for oxygen ions, independent of their charge state.
Birn, J.↗
The structure of dissipative dark matter halos
Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10–100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.
79 ASTRONOMY AND ASTROPHYSICS↗