Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗
DOE OSTI2020

High-energy ballistic electrons in low-pressure radio-frequency plasmas

This study demonstrates the presence of a small number of high-energy ballistic electrons (HEBEs) that originate from secondary electrons in low-pressure radio-frequency (rf) plasmas. The kinetic behaviors of the HEBEs are illustrated through electron energy probability functions from the fully kinetic particle-in-cell simulations, showing two wavy high-energy tails and two bifurcations during one rf cycle. Test-particle simulations and a semi-analytical method associated with nonlocal electron kinetics are performed to characterize the HEBE trajectories, which reveal the ballistic nature of the HEBEs and their typical bouncing features between the rf sheaths. Parameter dependence of the HEBEs on the discharge conditions (e.g., gas pressure, gap distance, and rf frequency) are identified, which is relevant to the plasma collisionality. With a pronounced presence of HEBEs, the overall impacts of the secondary electron emission on discharge parameters, such as electron power absorption and ionization rate, are also illustrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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

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

42 ENGINEERING↗
DOE OSTI2023

Mapping the performance of a versatile water-based condensation particle counter (vWCPC) with numerical simulation and experimental study

Accurate airborne aerosol instrumentation is required to determine the spatial distribution of ambient aerosol particles, particularly when dealing with the complex vertical profiles and horizontal variations of atmospheric aerosols. A versatile water-based condensation particle counter (vWCPC) has been developed to provide aerosol concentration measurements under various environments with the advantage of reducing the health and safety concerns associated with using butanol or other chemicals as the working fluid. However, the airborne deployment of vWCPCs is relatively limited due to the lack of characterization of vWCPC performance at reduced pressures. Given the complex combinations of operating parameters in vWCPCs, modeling studies have advantages in mapping vWCPC performance. In this work, we thoroughly investigated the performance of a laminar-flow vWCPC using COMSOL Multiphysics® simulation coupled with MATLAB™. We compared it against a modified vWCPC (vWCPC model 3789, TSI, Shoreview, MN, USA). Our simulation determined the performance of particle activation and droplet growth in the vWCPC growth tube, including the supersaturation, $D$ p,kel,0 (smallest size of particle that can be activated), $D$ p,kel,50 (particle size activated with 50 % efficiency) profile, and final growth particle size D d under wide operating temperatures, inlet pressures $P$ (30–101 kPa), and growth tube geometry (diameter $D$ and initiator length $L$ ini ). The effect of inlet pressure and conditioner temperature on vWCPC 3789 performance was also examined and compared with laboratory experiments. The COMSOL simulation result showed that increasing the temperature difference (Δ$T$) between conditioner temperature $T$ con and initiator $T$ ini will reduce $D$ p,kel,0 and the cut-off size $D$ p,kel,50 of the vWCPC. In addition, lowering the temperature midpoint ($T$ mid = $\frac{T_{con}+T_{ini}}{2}$ increases the supersaturation and slightly decreases the $D$ p,kel . The droplet size at the end of the growth tube is not significantly dependent on raising or lowering the temperature midpoint but significantly decreases at reduced inlet pressure, which indirectly alters the vWCPC empirical cut-off size. Our study shows that the current simulated growth tube geometry ($D$=6.3 mm and $L$ ini =30 mm) is an optimized choice for current vWCPC flow and temperature settings. The current simulation can more realistically represent the $D$ p,kel for 7 nm vWCPC and also achieved good agreement with the 2 nm setting. Using the new simulation approach, we provide an optimized operation setting for the 7 nm setting. This study will guide further vWCPC performance optimization for applications requiring precise particle detection and atmospheric aerosol monitoring.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Sustained Oxygen Spectral Gaps and Their Dynamic Evolution in the Inner Magnetosphere

Van Allen Probes observations of ion spectra often show a sustained gap within a very narrow energy range throughout the full orbit. To understand their formation mechanism, we statistically investigate the characteristics of the narrow gaps for oxygen ions and find that they are most frequently observed near the noon sector with a peak occurrence rate of over 30%. The magnetic moment (μ) of the oxygen ions in the gap shows a strong dependence on magnetic local time (MLT), with higher and lower μ values in the morning and afternoon sectors, respectively. Moreover, we find through superposed epoch analysis that the gap formation also depends on geomagnetic conditions. Those gaps formed at lower magnetic moments (μ < 3,000 keV/G) are associated with stable convection electric fields, which enable magnetospheric ions to follow a steady drift pattern that facilitates the gap formation by corotational drift resonance. On the other hand, gaps with higher μ values are statistically preceded by a gradual increase of geomagnetic activity. Here, we suggest that ions within the gap were originally located inside the Alfven layer following closed drift paths, before they were transitioned into open drift paths as the convection electric field was enhanced. The sunward drift of these ions, with very low fluxes, forms a drainage void in the dayside magnetosphere manifested as the sustained gap in the oxygen spectrum. This scenario is supported by particle-tracing simulations, which reproduce most of the observed characteristics and therefore provide new insights into inner magnetospheric dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Integrating HPC simulations and physical experiments to characterize the effects of gamma radiation on seismic protective devices

Seismic protective systems, composed of seismic isolators and dampers, can substantially reduce the effects of earthquake shaking on nuclear power plants and components therein. To enable the use of these devices to protect equipment inside a plant and close to a source of radiation, the U.S. Department of Energy (DOE) funded a project at the Idaho National Laboratory (INL) and the University at Buffalo to characterize the effects of absorbed gamma dose on their mechanical properties. An early task in the project was to determine the exposure time required in the INL Foss Therapy Services (FTS) 60 Co gamma irradiator to achieve a target absorbed dose in the materials used to construct isolators and dampers, including fluids, polymers, composites, and metals. This task required the novel integration of high-performance computing (HPC), Monte Carlo N-Particle (MCNP) simulations, and irradiation experiments using Fricke dosimetry. An MCNP model of the FTS irradiator at INL was developed and validated using Fricke dosimetry. Simulations of three experiments in the irradiator, two with Fricke vials only and one with Fricke vials and a large-size isolator, predicted the Fricke-measured absorbed dose rate to within 15% in all three cases, providing high confidence in the calculation of the gamma dose absorbed in the materials comprising the seismic protective devices. The simulations demonstrated that the effects of photon scattering on absorbed dose rate in the FTS irradiator are negligible for test articles installed close to the cobalt sources and near the rear of the irradiator. The validated MCNP model of the FTS irradiator is being used to support ongoing DOE-funded experiments on seismic protective devices and could be applied to future, non-seismic-related experiments. In conclusion, the novel validation process successfully deployed for the FTS irradiator at INL could be applied to other irradiators, requiring new MCNP models and simulations, and irradiation experiments using dosimeters.

42 - ENGINEERING↗
DOE OSTI2023

Reconstruction of interactions in the ProtoDUNE-SP detector with Pandora

The Pandora Software Development Kit and algorithm libraries provide pattern-recognition logic essential to the reconstruction of particle interactions in liquid argon time projection chamber detectors. Pandora is the primary event reconstruction software used at ProtoDUNE-SP, a prototype for the Deep Underground Neutrino Experiment far detector. ProtoDUNE-SP, located at CERN, is exposed to a charged-particle test beam. This paper gives an overview of the Pandora reconstruction algorithms and how they have been tailored for use at ProtoDUNE-SP. In complex events with numerous cosmic-ray and beam background particles, the simulated reconstruction and identification efficiency for triggered test-beam particles is above 80% for the majority of particle type and beam momentum combinations. Specifically, simulated 1 GeV/c charged pions and protons are correctly reconstructed and identified with efficiencies of 86.1$\pm 0.6$% and 84.1$\pm 0.6$%, respectively. The efficiencies measured for test-beam data are shown to be within 5% of those predicted by the simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Improving microstructures segmentation via pretraining with synthetic data

Image analysis of material microstructures through microscopy is an integral capability in the field of materials science. The topological and chemical information obtained through microscopy allow us to draw vital connections between material microstructures, properties, and processing. While scanning electron microscopy (SEM) is able to yield a considerable wealth of information interpretable by the intuition of experts, there has been considerable interest in using machine learning, convolutional neural networks (CNNs) in particular, for such image analysis task. Training CNNs for an image analysis task requires a large annotated dataset. However, in many materials science applications, obtaining a large annotated dataset is cost and labor intensive. In this work, we study the use of synthetic data to enlarge the available annotated experimental data of uranium oxide. We utilize a modified Potts model to simulate uranium oxide particles with morphologies similar to those observed experimentally. We then leverage an image-to-image translation model to synthesize the simulated particles as if they are acquired with SEM. Through this process, we obtain pairs of particle images and their corresponding SEM representations, which corresponds to pairs of annotations and images. Unlike previous works, we leverage synthetic data for pretraining a CNN model prior, and finetune that model further with experimental data. We experimentally demonstrate that using synthetic data as incremental learning process benefits the overall performance compared to training a model on combined synthetic and experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Three-Dimensional Particle-in-Cell Analysis of Mechanical Design Tolerances of the Scorpius Injector With Warp

Here, this article examines the mechanical design tolerances for the Scorpius injector using a 3-D particle-in-cell (PIC) code. In order to meet its objective as a high-fidelity, multipulse flash X-ray source, the accelerator has requirements for the electron beam in order to minimize the spot size of the beam focused onto a foil target to create high fluence X-rays at energies greater than 10 MeV. The requirements for alignment of the cathode, anode beam pipe, focusing solenoids, and transport solenoids specified using Trak and AMBER are independently verified using a 3-D electrostatic/magnetostatic model of the Scorpius injector, including transverse magnetic field errors that arise from mechanical alignment errors. It is predicted that the injector will meet the performance requirements for radiographic operation and has a design margin to correct for potentially unforeseen operating conditions.

42 ENGINEERING↗