Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle dynamics simulation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Advances in the simulation of space-charge effects

The nonlinear space-charge effects are an important topic in high intensity accelerator beam dynamics and have been extensively studied using macroparticle tracking simulations. Here, we report on recent advances in the simulation of space-charge effects using a symplectic multiparticle space-charge model. The transverse space-charge limit was explored for a periodic focusing and defocusing lattice. The artificial numerical emittance growth in the macroparticle space-charge simulation was analyzed using a one-dimensional model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Redesign of the Jefferson Lab -300 kV DC Photo-Gun for High Bunch Charge Operations

Production of high bunch charge beams for the Electron-Ion Collider (EIC) is a challenging task. High bunch charge (a few nC) electron beam studies at Jefferson Lab using an inverted insulator DC high voltage photo-gun showed evidence of space charge limitations starting at 0.3 nC, limiting the maximum delivered bunch charge to 0.7 nC for beam at -225 kV, 75 ps (FWHM) pulse width, and 1.64 mm (rms) laser spot size. The low extracted charge is due to the modest longitudinal electric field (Ez) at the photocathode leading to beam loss at the anode and downstream beam pipe. To reach the few nC high bunch charge goal, and to correct the beam deflection exerted by the non-symmetric nature of the inverted insulator photo-gun the existing photo-gun was modified. This contribution discusses the electrostatic design of the modified photo-gun obtained using CST Studio Suite’s electromagnetic field solver. Beam dynamics simulations performed using General Particle Tracer (GPT) with the resulting electrostatic field map obtained from the modified electrodes confirmed the validity of the new design.

Wijethunga, S.A.K.↗
DOE OSTI2023

Applying molecular dynamics simulations to the study of micro cold spray of yttria

Micro cold spray (MCS) is a novel manufacturing process with proven success depositing dense, thick films of technical ceramics without the need for elevated processing temperatures or chemical precursors. The process transforms a loose sub-micrometer powder into an adherent film with mechanical properties approaching that of the bulk material through a phenomenon known as room temperature impact consolidation. While this phenomenon has been demonstrated for a wide variety of ceramic materials and film substrates, yttria stands out as a material that is both well-characterized in experiments and of significant interest for industrial applications due to its excellent refractory properties. This dissertation seeks to investigate the fundamental mechanisms that enable room temperature impact consolidation by using molecular dynamics (MD) simulations of yttria particle impacts and characterizing the deformation behavior of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A versatile approach to stabilize liquid-liquid interfaces using surfactant self-assembly

Stabilizing liquid–liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. In this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid–liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Status and future plans for C 3 R&D

C 3 is an opportunity to realize an e + e - collider for the study of the Higgs boson at √s = 250 GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C 3 is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C 3 is the main linac. This paper presents the staged approach towards a facility to demonstrate C 3 technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The primary goal of the C 3 Demonstration R&D Plan is to reduce technical and cost risk by building and operating the key components of C 3 at an adequate scale. This R&D plan starts with the engineering design, and demonstration of one cryomodule and will culminate in the construction of a 3 cryomodule linac with pre-production prototypes. This R&D program would also demonstrate the linac rf fundamentals including achievable gradient and gradient stability over a full electron bunch train and breakdown rates. It will also investigate beam dynamics including energy spread, wakefields, and emittance growth. This work will be critical to confirm the suitability of the C 3 beam parameters for the physics reach and detector performance in preparation for a Conceptual Design Report (CDR), as well as for follow-on technology development and industrialization. The C 3 Demonstration R&D Plan will open up significant new scientific and technical opportunities based on development of high-gradient and high-efficiency accelerator technology. It will push this technology to operate both at the GeV scale and mature the technology to be reliable and provide high-brightness electron beams. The timeline for progressing with C 3 technology development will be governed by practical limitations on both the technical progress and resource availability. It consists of four stages: Stage 0) Ongoing fundamental R&D on structure prototypes, damping and vibrations. Stage 1) Advancing the engineering maturity of the design and developing start-to-end simulations including space-charge and wakefield effects. This stage will include testing of strucutres operating at cryogenic temperatures. Beam tests would be performed with high beam current to test full beam loading. Stage 2) Production and testing of the first cryomodule at cryogenic temperatures. This would provide sufficient experimental data to compile a CDR and it is anticipated for Stage 2 to last 3 years and to culminate with the transport of photo-electrons through the first cryomodule. Stage 3) Updates to the engineering design of the cryomodules, production of the second and third cryomodule and their installation. Lower charge and lower emittance beams will be used to investigate emittance growth. The successful full demonstration of the 3 cryomodules to deliver up to a 3 GeV beam and achieve the C 3 five gradient will allow a comprehensive and robust evaluation of the technical design of C 3 as well as mitigate technical, schedule, and cost risks required to proceed with a Technical Design Report (TDR).

radiation hardened magnets↗
DOE OSTI2023

Molecular dynamics analysis of particle number fluctuations in the mixed phase of a first-order phase transition

Molecular dynamics simulations are performed for a finite nonrelativistic system of particles with Lennard-Jones potential. We study the effect of liquid-gas mixed phase on particle number fluctuations in coordinate subspace. A metastable region of the mixed phase, the so-called nucleation region, is analyzed in terms of a noninteracting cluster model. Large fluctuations due to spinodal decomposition are observed. They arise due to the interplay between the size of the acceptance region and that of the liquid phase. These effects are studied with a simple geometric model. The model results for the scaled variance of particle number distribution are compared with those obtained from the direct molecular dynamic simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2023

Quantum simulation of fundamental particles and forces

Key static and dynamic properties of matter — from creation in the Big Bang to evolution into subatomic and astrophysical environments — arise from the underlying fundamental quantum fields of the standard model and their effective descriptions. However, the simulation of these properties lies beyond the capabilities of classical computation alone. Advances in quantum technologies have improved control over quantum entanglement and coherence to the point at which robust simulations of quantum fields are anticipated in the foreseeable future. In this Perspective article, we discuss the emerging area of quantum simulations of standard-model physics, outlining the challenges and opportunities for progress in the context of nuclear and high-energy physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Effect of an oxide layer on high velocity impact of tantalum particles characterized using molecular dynamics

Micro-cold spray, also commonly referred to as the aerosol deposition method or vacuum kinetic spraying, is a process for depositing nanostructured films of metals and ceramics by impacting solid nanoparticles at high velocity onto a substrate at room temperature and at low pressures. Tantalum films deposited by micro-cold spray are desirable because tantalum is non-reactive to certain molten metals and can therefore serve as a non-reactive barrier coating. However, unlike most other metals deposited by micro-cold spray, tantalum quickly forms an oxide layer when exposed to air. In the deposition of films using traditional cold spray with much larger particles, native oxides must fracture for particle–substrate bonding to occur. Furthermore, the fraction of the particle that is oxidized in micro-cold spray and the strain rates experienced by the particles upon impact are both much larger than in cold spray, suggesting that the role of the oxide on film deposition may be distinct. The effect of this oxide layer on particle deformation and adherence to a substrate under high strain rates is studied using molecular dynamics simulations. The effect of oxide layer thickness, particle diameter, and impact velocity are examined.

42 ENGINEERING↗
DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

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↗
DOE OSTI2024

Enhanced beam-beam modeling to include longitudinal variation during weak-strong simulation

Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multiscale CFD simulation of biomass fast pyrolysis with a machine learning derived intra-particle model and detailed pyrolysis kinetics

Coupling particle and reactor scale models is as essential as reactor fluid dynamics and particle motion for accurate Computational Fluid Dynamic (CFD) simulations of biomass fast pyrolysis reactors due to intraparticle heat transfer and chemical reactions controlling conversion time and product distributions. Direct online coupling of a particle model with a reactor model is computationally expensive, while offline coupling is case-dependent. In this research, solutions from a series of particle pyrolysis simulations were regressed with Artificial Neural Network (ANN). Furthermoer, this machine learning-derived model predicted the same temperature and conversion profiles compared with particle resolved simulation while the isothermal approach overpredicted the temperature by 130 K and underpredicted the conversion time by 30 s. The ANN model was then integrated into CFD simulations of fluidized bed biomass fast pyrolysis with varied feedstocks via coupling PyTorch and MFiX. The averaged error of simulation predicted bio-oil yields with four feedstocks is 6.4%. This multi-scale approach provides an efficient tool for the coupled particle and reactor scale simulations of biomass pyrolysis.

09 BIOMASS FUELS↗
DOE OSTI2022

Single electron in a storage ring: a probe into the fundamental properties of synchrotron radiation and a powerful diagnostic tool

An experimental study of undulator radiation generated by a single electron was carried out at the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. Photons were detected with a single-photon avalanche diode (SPAD) at an average rate of 20 kHz, corresponding to one detection every few hundred revolutions. Events were continuously recorded by a picosecond event timer for as long as 5 minutes at a time. The collected data were used to directly test the predictions of classical Poisson photostatistics. The study was motivated in part by a previous observation of sub-Poissonian statistics in a similar experiment with synchrotron radiation. In addition, we show that the time series of recorded events can be used to study the synchrotron motion of a single electron, yielding dynamical parameters such as the phase jitter of the radiofrequency cavity and the dependence of the synchrotron period on amplitude. Research plans for further experiments are also discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Simulation of Electron-Proton Scattering Events by a Feature-Augmented and Transformed Generative Adversarial Network (FAT-GAN)

We apply generative adversarial network (GAN) technology to build an event generator that simulates particle production in electron-proton scattering that is free of theoretical assumptions about underlying particle dynamics. The difficulty of efficiently training a GAN event simulator lies in learning the complicated patterns of the distributions of the particles physical properties. We develop a GAN that selects a set of transformed features from particle momenta that can be generated easily by the generator, and uses these to produce a set of augmented features that improve the sensitivity of the discriminator. The new Feature-Augmented and Transformed GAN (FAT-GAN) is able to faithfully reproduce the distribution of final state electron momenta in inclusive electron scattering, without the need for input derived from domain-based theoretical assumptions. The developed technology can play a significant role in boosting the science of existing and future accelerator facilities, such as the Electron-Ion Collider.

Alanazi, Yasir↗
DOE OSTI2024

Modeling Flow and Particle Deposition in a Spent Nuclear Fuel Assembly

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. Here, we perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗