Search NASA⌕ Search

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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Graph Neural Network for Object Reconstruction in Liquid Argon Time Projection Chambers

The Exa.TrkX project presents a graph neural network (GNN) technique for low-level reconstruction of neutrino interactions in a Liquid Argon Time Projection Chamber (LArTPC). GNNs are still a relatively novel technique, and have shown great promise for similar reconstruction tasks in the Large Hadron Collider (LHC). Graphs describing particle interactions are formed by treating each detector hit as a node, with edges describing the relationships between hits. We utilise a multi-head attention message passing network which performs graph convolutions in order to label each node with a particle type.We present an updated variant of our GNN architecture, with several improvements. After testing the model on more realistic simulation with regions of unresponsive wires, the target was modified from edge classification to node classification in order to increase robustness. Removing edges as a classification target opens up a broader possibility space for edge-forming techniques; we explore the model’s performance across a variety of approaches, such as Delaunay triangulation, kNN, and radius-based methods. We also extend this model to the 3D context, sharing information between detector views. By using reconstructed 3D spacepoints to map detector hits from each wire plane, the model naively constructs 2D representations that are independent yet fully consistent.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

COZMIC. III. Cosmological Zoom-in Simulations of Self-interacting Dark Matter with Suppressed Initial Conditions

We present eight cosmological dark matter (DM)-only zoom-in simulations of a Milky Way–like system that include suppression of the linear matter power spectrum P(k), and/or velocity-dependent DM self-interactions, as the third installment of the COZMIC suite. We consider a model featuring a massive dark photon that mediates DM self-interactions and decays into massless dark fermions. The dark photon and dark fermions suppress linear matter perturbations, resulting in dark acoustic oscillations in P(k), which ultimately affect dwarf galaxy scales. The model also features a velocity-dependent elastic self-interaction between DM particles (SIDM), with a cross section that can alleviate small-scale structure anomalies. For the first time, our simulations test the impact of P(k) suppression on gravothermal evolution in an SIDM scenario that leads to core collapse in (sub)halos with present-day virial masses below ≈10 9 M ⊙ . In simulations with P(k) suppression and self-interactions, the lack of low-mass (sub)halos and the delayed growth of structure reduce the fraction of core-collapsed systems relative to SIDM simulations without P(k) suppression. In particular, P(k) suppression that saturates current warm DM constraints almost entirely erases core collapse in isolated halos. Models with less extreme P(k) suppression produce core collapse in ≈20% of subhalos and ≈5% of isolated halos above 10 8 M ⊙ , and also increase the abundance of extremely low-concentration isolated low-mass halos relative to SIDM. These results reveal a complex interplay between early and late-Universe DM physics, revealing new discovery scenarios in the context of upcoming small-scale structure measurements.

dark matter↗
DOE OSTI2025

Forecasting Dark Matter Subhalo Constraints from Stellar Streams using Implicit Likelihood Inference

The evidence for dark matter (DM) remains compelling, although attempts to understand its particle nature remain inconclusive. One promising method to study DM is detecting DM subhalos through their gravitational interactions with stellar streams. In this study, we apply Neural Posterior Estimation (NPE) to constrain subhalo interaction parameters, including mass, scale radius, velocity, and encounter geometry, from stellar stream kinematics. We generate particle spray simulations based on the Lagrange Cloud stripping technique, focusing on the ATLAS-Aliqa Uma stream as a test case. We train multiple NPE models across multiple observational scenarios, quantifying how kinematic completeness affects inference and forecasting constraints from upcoming surveys including LSST, 4MOST, and 10-year Gaia data. Our results demonstrate that NPE can produce accurate and well-calibrated posteriors. In the idealized case with full 6D coordinates, we achieve subhalo mass uncertainties of 15-20% for a $10^7 \, \mathrm{M_\odot}$ subhalo, with 5D coordinates (excluding radial velocities) achieving similar performance. Under realistic observational conditions, mass uncertainties range from 50% (present-day) to 20-40% (future scenarios), with comparable performance between the photometric-only LSST sample and a smaller sample that includes Gaia proper motions and 4MOST radial velocities. Most notably, we find that velocity bimodality emerges when phase space is poorly sampled, whether due to missing kinematic information or limited stellar tracers. Combining large photometric samples with targeted spectroscopic follow-up can effectively resolves this degeneracy. These results demonstrate the power of implicit likelihood inference for optimizing stellar stream observational strategies and forecasting DM subhalo constraints from upcoming surveys.

Nguyen, Tri [Northwestern U. (main); SkAI, Chicago↗
DOE OSTI2022

Learning to simulate high energy particle collisions from unlabeled data

In many scientific fields which rely on statistical inference, simulations are often used to map from theoretical models to experimental data, allowing scientists to test model predictions against experimental results. Experimental data is often reconstructed from indirect measurements causing the aggregate transformation from theoretical models to experimental data to be poorly-described analytically. Instead, numerical simulations are used at great computational cost. We introduce Optimal-Transport-based Unfolding and Simulation (OTUS), a fast simulator based on unsupervised machine-learning that is capable of predicting experimental data from theoretical models. Without the aid of current simulation information, OTUS trains a probabilistic autoencoder to transform directly between theoretical models and experimental data. Identifying the probabilistic autoencoder’s latent space with the space of theoretical models causes the decoder network to become a fast, predictive simulator with the potential to replace current, computationally-costly simulators. Here, we provide proof-of-principle results on two particle physics examples, Z-boson and top-quark decays, but stress that OTUS can be widely applied to other fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

SPH modeling of biomass granular flow: Theoretical implementation and experimental validation

The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.

09 BIOMASS FUELS↗
DOE OSTI2024

Method for Real-Time Salt Aerosol Concentration and Size Measurements for Molten Salt Reactor Safety Assessments

Licensing a new nuclear reactor involves evaluating credible accident scenarios by developing models that simulate accident progression to predict possible outcomes and their impact on safety. The formation of radionuclide-bearing aerosols in the respirable size range can significantly affect dose consequence and threaten human health and is therefore a focus of nuclear reactor safety assessments. Recent reviews of the literature identified a lack of existing experimental data describing the mechanisms of formation and properties (size, concentration, and composition) of radionuclide-bearing aerosols that may be produced from molten salt reactor (MSR) facilities during postulated accident scenarios. Experiments on aerosol formation from molten salt systems are a high priority need that will support MSR licensing by indicating the radiological consequences of aerosol formation and providing the data required for model development and validation. The evolving conditions that occur during a MSR accident, such as a spill of molten fuel salt, may affect aerosol formation mechanisms and aerosol properties over time. There is a need to conduct experiments that simulate credible MSR accidents in a laboratory to generate aerosols with realistic characteristics and behaviors and a need to measure these aerosols in real time under accident-relevant conditions. This report describes the development of a method that can be used to quantify the size and concentration of salt aerosol particles that form from molten salt systems in real time. This method will be employed in future integral effects tests that are conducted at an engineering scale to simulate realistic MSR accidents and in future separate effects tests that will provide mechanistic insight into aerosol formation and properties to support process model development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Fast dynamic aperture optimization with forward-reversal integration

A fast dynamic aperture (DA) optimization method for storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symmetry, a numerical forward integration differs from its reversal. For a chaotic trajectory, cumulative round-off errors are scaled, which results in an exponential growth on the difference. The exponential effect is a generic chaos indicator which represents the sensitivity of the chaotic motion to its initial condition. The chaos indicator of the charged particle motion can be obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. “backward integrations.” The indicator is observable even through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. Finally, the DA of the National Synchrotron Light Source II storage ring, and another test diffraction-limited light source ring, were optimized using this method for the purpose of demonstration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Integrated Process Testing of MSR Salt Spill Accidents

Part of the licensing process for new nuclear reactors requires vendors to assess the potential consequences of identified accident scenarios using accident progression modeling. The accident scenario that will likely be evaluated by all molten salt reactor (MSR) developers is a spill of radionuclide-bearing fuel salt onto the reactor containment floor (i.e., a salt spill accident). The development of accident progression models requires experimental data to inform which processes to incorporate, to enable the calculation of parameters to model these processes, and to validate the model predictions. The data should quantify the sensitivities of key processes (e.g., molten salt spreading, heat transfer, containment structure corrosion, radionuclide vaporization, and aerosol generation) towards the initial conditions of the spill, the ambient environment, and the features of the containment. In addition, results from integrated process tests that quantify coupled processes are required to validate systems-level models. This report documents results from integrated process tests conducted on simulated molten salt spill accidents. The generated experimental data simultaneously quantify the heat transfer behavior of the spilled salt, compositional changes to the bulk salt, and the release of surrogate fission products from the spilled salt as aerosol particles. All tests that were conducted used FLiNaK doped with surrogate fission products, and the variables that were evaluated included the initial salt temperature and the concentration of surrogate fission products present in the salt. The major accomplishments of this work include identifying surrogate fuel salt compositions that provide insight into the dispersal behavior of radionuclides of potential significance to the source term, employing previously developed methods and measurement techniques to simultaneously measure key processes, generating data on the coupled processes of molten salt heat transfer and surrogate fission product release as aerosol particles, demonstrating new test methods for real-time monitoring of the flow rate of the spill and aerosol size quantification in an argon atmosphere, and developing a mass transfer model for cesium and iodine release from molten FLiNaK to provide insight into aerosol formation by vapor condensation. The same methodology applied herein can be employed to study different salt compositions of interest to MSR developers, different environmental conditions, and other variables that are relevant to postulated accident scenarios. The insights gained from these integrated process tests conducted at a laboratory scale will be incorporated into future integral effects tests conducted at an engineering scale.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗
DOE OSTI2020

Design and Analysis of the Staggered Double Grating Slow Wave Circuit for 263 GHz Sheet Beam TWT

The complete design process of a high-frequency circuit for a 263 GHz staggered double grating, sheet beam traveling wave tube is presented in this article. Here, this device, for use in a pulsed electron paramagnetic resonance spectrometer instrument, requires >30 dB gain and a 20 GHz bandwidth centered at 263 GHz. In addition to the standard Pierce gain analysis and point by point particle in cell (PIC) simulation, a novel fast pulse response analysis technique was utilized to determine the accurate synchronization voltage and test the ``hot” feature of the high-frequency circuit. Since this new method can evaluate the entire bandwidth in a single simulation, it can significantly reduce the computational resources normally required for the conventional PIC analysis while also illustrating the potential instabilities near the band edge. Good agreement between the fast pulse response and PIC simulation results validate the accuracy of this novel method. According to the analysis results and given the manufacturing limits of nano/micro-computer numerical control (CNC) machining, an appropriate SWS, coupler, and sever designs were chosen to complement the simulated design. Finally, a prototype circuit was manufactured and cold tested to verify the design and inspect the machining tolerances.

42 ENGINEERING↗
DOE OSTI2021

Simulation-based methodology to assess the lattice defects creation as energy storing process

In nuclear industry, damages caused to in-core structural materials by high-energy neutrons have been studied for years. Recently, the possibility of storing energy through the creation of defects in crystalline materials via heavy particle bombardment was investigated. To evaluate the potential of the Wigner effect as energy-storing process, a graphite-moderated Molten Salt Reactor was adopted as a test case. Neutronics simulations performed with OpenMC Monte Carlo code allowed estimating the amount of energy released during fission events, the fraction that is deposited in graphite, and the fraction that has the potential to cause radiation damages. Additionally, molecular dynamics simulations were then performed to study the evolution of cascades triggered by high-energy neutrons. As a major outcome, a more accurate estimate of the radiation damage with respect to the traditional approaches was obtained, and the fraction of the deposited neutron energy that can be stored through long-lived defects was quantified.

25 ENERGY STORAGE↗
DOE OSTI2024

Macroscopic approach to the radar echo scatter from high-energy particle cascades

To probe the cosmic particle flux at the highest energies, large volumes of dense material like ice have to be monitored. This can be achieved by exploiting the radio signal. In this work, we provide a macroscopic model to predict the radar echo signatures found when a radio signal is reflected from a cosmic-ray or neutrino-induced particle cascade propagating in a dense medium like ice. Its macroscopic nature allows for an energy independent run-time, taking less than 10 s for simulating a single scatter event. As a first application, we discuss basic signal properties and simulate the expected signal for the T-576 beam-test experiment at the Stanford Linear Accelerator Center. Finally, we find good signal strength agreement with the only observed radar echo from a high-energy particle cascade to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GPU-accelerated kinetic Hall thruster simulations in WarpX

Abstract Two-dimensional (axial-azimuthal) simulations of a Hall thruster are performed using the open-source particle-in-cell code WarpX. The simulation conditions are chosen to match those of the axial-azimuthal benchmark first reported by Charoy et al. in 2019. A range of numerical and solver parameters is investigated in order to find those which yield the best performance. It is found that WarpX completes the benchmark case in 3.8 days on an Nvidia V100 GPU, and in as low as 1.5 days on a more recent Nvidia H100 GPU. Of the numerical parameters investigated, it is determined that the field-solve tolerance and particle resampling thresholds have the largest effect on the simulation wall time and that particle resampling may artificially widen electron velocity distribution functions, leading to unphysical heating. A semi-implicit scheme for the electrostatic field solve is tested and is found to produce results consistent to within 10% of the benchmark in less than twelve hours. The scaling properties of the electrostatic solver to multiple GPUs are also assessed on a uniform plasma test problem. The results of this work are discussed in the context of advancements in GPU hardware and the suitability of kinetic Hall thruster simulations for engineering applications.

Marks, Thomas A.↗
DOE OSTI2023

Cladding and Butt-Joining Dissimilar Aluminum Alloys Simultaneously via Shear Assisted Processing and Extrusion

Shear assisted processing and extrusion has been used to fabricate multi-metallic tubing that incorporates regions with clad and butt-joined structures. In this work, aluminum alloys 6061, 7075, and 1100 are combined to create a monolithic extruded tubing with weldable ends (6061), a high strength center (7075), and a corrosion-resistant cladding (1100). Optical microscopy shows a gradual, defect-free, transition between alloys at the butt-joint interfaces and a sharp transition between alloys at the clad interface. Tensile testing proved the mechanical coherency of the joint, with failure occurring outside the tri-metallic structure. In conclusion, smoothed particle hydrodynamic simulations identified a spiral pattern of the cladding alloy that leads to inconsistent cladding thickness, and identified a simplified 2-alloy combination of aluminum alloys 6082 and 2014 that accomplished the design objectives of the tri-metallic structure.

Coatings↗
DOE OSTI2022

The influence of multiple groups of biological ice nucleating particles on microphysical properties of mixed-phase clouds observed during MC3E

Abstract. A new empirical parameterization (EP) for multiple groups of primary biological aerosol particles (PBAPs) is implemented in the aerosol–cloud model (AC) to investigate their roles as ice nucleating particles (INPs). The EP describes the heterogeneous ice nucleation by (1) fungal spores, (2) bacteria, (3) pollen, (4) detritus of plants, animals, and viruses, and (5) algae. Each group includes fragments from the originally emitted particles. A high-resolution simulation of a midlatitude mesoscale squall line by AC is validated against airborne and ground observations. Sensitivity tests are carried out by varying the initial vertical profiles of the loadings of individual PBAP groups. The resulting changes in warm and ice cloud microphysical parameters are investigated. The changes in warm microphysical parameters, including liquid water content and cloud droplet number concentration, are minimal (<10 %). Overall, PBAPs have little effect on the ice number concentration (<6 %) in the convective region. In the stratiform region, increasing the initial PBAP loadings by a factor of 1000 resulted in less than 40 % change in ice number concentrations. The total ice concentration is mostly controlled by various mechanisms of secondary ice production (SIP). However, when SIP is intentionally shut down in sensitivity tests, increasing the PBAP loading by a factor of 100 has an effect of less than 3 % on the ice phase. Further sensitivity tests revealed that PBAPs have little effect on surface precipitation and on the shortwave and longwave flux (<4 %) for a 100-fold perturbation in PBAPs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Background studies and normalization of signal events in the Mu2e experiment

The Mu2e experiment is currently being constructed at Fermilab to search for the neutrino-less conversion of negative muons into electrons in the field of an aluminum nucleus. The experiment aims at a sensitivity of four orders of magnitude higher than previous related experiments, which implies highly demanding accuracy requirements both in the design and during the operation. To achieve such a goal, two important tasks should be accomplished. First, it is essential to estimate precisely the particle yields and all the backgrounds that could mimic the monoenergetic conversion electron signal. Second, it is necessary to normalize the signal events accurately. The normalization of the signal events is planned to be done using a detector system made of an HPGe detector and a Lanthanum Bromide detector, which will measure the rate of muons stopped on the aluminum target by looking at the emitted characteristic X-and γ-rays of energies up to 1809 keV. Therefore, it is essential to evaluate the detector system's performance before the start of the actual experiment. In this study, the first task was addressed by an extensive campaign of Monte Carlo simulations to investigate the relevant parameters and their impact on the experiment's sensitivity. The second task was handled by taking advantage of the Helmholtz-Zentrum Dresden-Rossendorf (HZDR) pulsed Bremsstrahlung photon beam at the ELBE facility. The detector system was tested at the ELBE facility under timing and background conditions similar to the ones expected at the Mu2e experiment. The study presents and discusses the simulation results and the detector system testing campaign.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

HIGH BURNUP FUEL-COOLANT INTERACTION ANALYSIS SUPPORTING FUEL SAFETY TESTING AT IDAHO NATIONAL LABORATORY

In the near future, experiments on HBu fuel under loss-of-coolant accident (LOCA) and reactivity-initiated accident (RIA) conditions will be performed within the Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL). These experiments will be performed using the Transient Water Irradiation System for TREAT (TWIST) experiment vehicle. To support these experiments, analysis of fuel-coolant interaction (FCI) energetics is underway. This paper discusses FCIs in the context of light water reactor (LWR) safety, differentiating between the severe accident focus of commercial reactors and experimental RIA test programs where FCIs have occurred. However, it is highlighted that as the nuclear industry aims for increased burnup limits, the FCI events observed in RIA test programs may become relevant to commercial LWR safety analysis. The paper then presents developments to the UW-FCI computer program to enable simulation of FCIs initiated by solid fuel particles dispersing into the coolant during RIAs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗