Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator 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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Particle acceleration in relativistic magnetic reconnection with strong inverse-Compton cooling in pair plasmas

Particle-in-cell (PIC) simulations have shown that relativistic collisionless magnetic reconnection drives non-thermal particle acceleration (NTPA), potentially explaining high-energy (X-ray/γ-ray) synchrotron and/or inverse Compton (IC) radiation observed from various astrophysical sources. The radiation back-reaction force on radiating particles has been neglected in most of these simulations, even though radiative cooling considerably alters particle dynamics in many astrophysical environments where reconnection may be important. We present a radiative PIC study examining the effects of external IC cooling on the basic dynamics, NTPA, and radiative signatures of relativistic reconnection in pair plasmas. We find that, while the reconnection rate and overall dynamics are basically unchanged, IC cooling significantly influences NTPA: the particle spectra still show a hard power law (index ≥ -2) as in non-radiative reconnection, but transition to a steeper power law that extends to a cooling-dependent cut-off. The steep power-law index fluctuates in time between roughly -3 and -5. The time-integrated photon spectra display corresponding power laws with indices ≈-0.5 and ≈-1.1, similar to those observed in hard X-ray spectra of accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Expanding the representation of aerosol, cloud, and precipitation processes with graph network-based simulators

We explored a novel framework for simulating the small-scale processes that drive the evolution of aerosol, cloud, and precipitation particles, which are a critical gap in the predictive understanding of weather and climate. Particle-based methods have emerged as an effective tool for modeling aerosol-cloud-precipitation interactions, but existing particle-based models are computationally too expensive to simulate the large domains relevant for the atmosphere or to represent the full suite of relevant processes. The lack of a comprehensive and efficient reference model is a critical bottleneck in our understanding of cloud and precipitation processes and our ability to parameterize these processes for regional- and global-scale simulations. To address this need, we explored an approach to accelerate and expand particle-based models using a new machine learning approach, graph network-based simulators (GNS). Rather than modeling the evolution of the system by numerically integrating continuity equations, the GNS represents dynamics through learned message passing. Our aim was to develop fast and accurate surrogate models for particle-based simulations. We explored applying GNS to simulate cloud droplet transport, growth, and evaporation under turbulent conditions, but we found the GNS over-smoothed the simulations. We then applied the GNS to simulate aerosol dynamics through gas condensation and found the GNS was able to reproduce the benchmark, physics-based simulation with high accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Infrastructure improvements in the National Reactor Innovation Center Virtual Test Bed

The National Reactor Innovation Center’s (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. The VTB is being developed in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished via three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: • A documentation website detailing the models: https://mooseframework.inl.gov/virtual_test_bed • A GitHub repository that hosts the corresponding files: https://github.com/idaholab/virtual_test_bed. This paper presents the infrastructure added to the VTB in the last two fiscal years. Additional information about the VTB can be found in various publications [1, 2, 3, 4, 5, 6, 7, 8].

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Status of New Models Hosted on the Virtual Test Bed (VTB) in 2024

The National Reactor Innovation Center (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. It is in development in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished by three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: 1. A documentation website detailing the models (https://mooseframework.inl.gov/virtual_test_bed). 2. A GitHub repository that hosts the corresponding files (https://github.com/idaholab/virtual_test_bed). Previous documentation on the models hosted in the VTB can be found in [1,2,3,4]. These references also include additional background information on the various NEAMS codes showcased in the VTB (which is omitted here for brevity). This paper primarily provides a status update of the most recent additions to the repository.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Mechanism of the G-protein mimetic nanobody binding to a muscarinic G-protein-coupled receptor

Protein–protein binding is key in cellular signaling processes. Molecular dynamics (MD) simulations of protein–protein binding, however, are challenging due to limited timescales. In particular, binding of the medically important G-protein-coupled receptors (GPCRs) with intracellular signaling proteins has not been simulated with MD to date. Here, we report a successful simulation of the binding of a G-protein mimetic nanobody to the M 2 muscarinic GPCR using the robust Gaussian accelerated MD (GaMD) method. Through long-timescale GaMD simulations over 4,500 ns, the nanobody was observed to bind the receptor intracellular G-protein-coupling site, with a minimum rmsd of 2.48 Å in the nanobody core domain compared with the X-ray structure. Binding of the nanobody allosterically closed the orthosteric ligand-binding pocket, being consistent with the recent experimental finding. In the absence of nanobody binding, the receptor orthosteric pocket sampled open and fully open conformations. The GaMD simulations revealed two low-energy intermediate states during nanobody binding to the M 2 receptor. The flexible receptor intracellular loops contribute remarkable electrostatic, polar, and hydrophobic residue interactions in recognition and binding of the nanobody. As a result, these simulations provided important insights into the mechanism of GPCR–nanobody binding and demonstrated the applicability of GaMD in modeling dynamic protein–protein interactions.

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

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

MULTI-WAVELENGTH EMISSION FROM THE FERMI BUBBLE. III. STOCHASTIC (FERMI) RE-ACCELERATION OF RELATIVISTIC ELECTRONS EMITTED BY SNRs

We analyze the model of stochastic re-acceleration of electrons that are emitted by supernova remnants (SNRs) in the Galactic Disk and then propagate into the Galactic Halo in order to explain the origin of nonthermal (radio and gamma-ray) emission from Fermi bubbles (FB). We assume that the energy for re-acceleration in the Halo is supplied by shocks generated by processes of star accretion onto the central black hole. Numerical simulations show that regions with strong turbulence (places for electron re-acceleration) are located high up in the Galactic Halo several kpc above the disk. The energy of the SNR electrons that reach these regions does not exceed several GeV due to synchrotron and inverse Compton energy losses. At appropriate parameters of re-acceleration these electrons can be re-accelerated up to an energy of 10{sup 12} eV, which explains in this model the origin of the observed radio and gamma-ray emission from the FB. However, although the model gamma-ray spectrum is consistent with the Fermi results, the model radio spectrum is steeper than that observed by WMAP and Planck. If adiabatic losses due to plasma outflows from the Galactic central regions are taken into account, then the re-acceleration model nicely reproduces the Planck data points.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Toward fusion plasma scenario planning for NSTX-U using machine-learning-accelerated models

One of the most promising devices for realizing power production through nuclear fusion is the tokamak. To maximize performance, it is preferable that tokamak reactors achieve advanced operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic (MHD) stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high \emph{fusion gain} --- the ratio of produced fusion power to the external power provided through the plasma boundary. Precise and robust control of the evolution of the plasma boundary shape as well as the spatial distribution of the plasma current, density, temperature, and rotation will be essential to achieving and maintaining such scenarios. The complexity of the evolution of tokamak plasmas, arising due to nonlinearities and coupling between various parameters, motivates the use of model-based control algorithms that can account for the system dynamics. In this work, a learning-based accelerated model trained on data from the National Spherical Torus Experiment Upgrade (NSTX-U) is employed to develop planning and control strategies for regulating the density and temperature profile evolution around desired trajectories. The proposed model combines empirical scaling laws developed across multiple devices with neural networks trained on empirical data from NSTX-U and a database of first-principles-based computationally intensive simulations. The reduced execution time of the accelerated model will enable practical application of optimization algorithms and reinforcement learning approaches for scenario planning and control development. An initial demonstration of applying optimization approaches to the learning-based model is presented, including a strategy for mitigating the effect of leaving the finite validity range of the accelerated model. The approach shows promise for actuator planning between experiments and in real-time.

machine learning↗
DOE OSTI2014

An improved Monte Carlo (MC) dose simulation for charged particle cancer therapy

Heavy-particle therapy such as carbon ion therapy are more popular nowadays because of the nature characteristics of charged particle and almost no side effect to patients. An effective treatment is achieved with high precision of dose calculation, in this research work, Geant4 based Monte Carlo simulation method has been used to calculate the radiation transport and dose distribution. The simulation have the same setting with the treatment room in Heavy Ion Medical Accelerator, HIMAC. The carbon ion beam at the isocentric gantry nozzle for the therapeutic energy of 290 MeV/u was simulated, experimental work was carried out in National Institute of Radiological Sciences, NIRS, Chiba, Japan by using the HIMAC to confirm the accuracy and qualities dose distribution by MC methods. The Geant4 based simulated dose distribution were verified with measurements for Bragg peak and spread out Bragg peak (SOBP) respectively. The verification of results shows that the Bragg peak depth-dose and SOBP distributions in simulation has good agreement with measurements. In overall, the study showed that Geant4 based can be fully applied in the heavy-ion therapy field for simulation, further works need to be carry on to refine and improve the Geant4 MC simulations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

COLA with scale-dependent growth: applications to screened modified gravity models

We present a general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers all the main examples of screening mechanisms in the literature. We test the code by comparing it to full simulations of two popular modified gravity models, namely f ( R ) gravity and nDGP, and find good agreement in the modified gravity boost-factors relative to ΛCDM even when using a fairly small number of COLA time steps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Novel Materials R&D for Next-Generation Accelerator Target Facilities

High-Entropy Alloys and Electrospun Nanofiber materials are two classes of novel materials that can offer improved resistance to beam-induced radiation damage and thermal shock. Research to develop these new materials specifically for multi-megawatt accelerator target applications, such as beam windows and particle-production targets, are ongoing at Fermilab within the scope of a DOE Early Career Research Program. The research program combines in-beam experiments with complementary simulations to tailor the microstructures of these novel materials for use in next-generation accelerator target facilities. Iterative simulations to optimize the material composition, physics performance, and beam-induced thermomechanical response will guide the material design and fabrication processes based on established figures of merit. This will be followed by material irradiation experiments using low-energy ions and prototypic high-energy protons with extensive post-irradiation material characterization to assess and qualify the selected novel materials. This talk will describe the alloy design and synthesis, microstructural pre-characterization of the alloys, and plans for the eventual down selection following low-energy ion irradiation studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Numerical Analysis of Parasitic Crossing Compensation with Wires in DA$\Phi$NE

Current-bearing wire compensators were successfully used in the 2005-2006 run of the DAΦNE collider to mitigate the detrimental effects of parasitic beam-beam interactions. A marked improvement of the positron beam lifetime was observed in machine operation with the KLOE detector. In view of the possible application of wire beam-beam compensators for the High Luminosity LHC upgrade, we revisit the DAΦNE experiments. We use an improved model of the accelerator with the goal to validate the modern simulation tools and provide valuable input for the LHC upgrade project.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Physics of the Aurora: Laboratory Measurements of Electron Acceleration by Inertial Alfven Waves (Final Technical Report)

The physics of the aurora is one of the foremost unsolved problems of space physics. The mechanisms responsible for accelerating electrons that precipitate onto the ionosphere are not fully understood. For more than three decades, particle interactions with Alfven waves have been proposed as a possible means for accelerating electrons and generating auroras. Due to the limitations of spacecraft conjunction studies and other multi-spacecraft approaches, it is unlikely that it will ever be possible, through spacecraft observations alone, to confirm definitively the proposed electron acceleration mechanism by making simultaneous measurements of both the accelerated electrons and the inertial Alfven wave responsible for the acceleration. Here, we show laboratory measurements of the resonant transfer of energy from Alfven waves to electrons under conditions relevant to the auroral zone. Experiments are performed by launching Alfven waves and simultaneously recording the electron velocity distribution. Numerical simulations and analytical theory support that the measured energy transfer process produces accelerated electrons capable of reaching auroral energies. The experiments, theory, and simulations demonstrate a clear causal relationship between Alfven waves and accelerated electrons that directly cause auroras.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of radiation damage via alpha decay in BFS:PC grouts using 4 He 2+ ion acceleration

The impact of alpha radiation on cements used to encapsulate intermediate-level waste (ILW) is not well understood. ILW wastes can contain high levels of alpha-emitting radionuclides, meaning that the grouts used to encapsulate them are exposed to significant ionising radiation. Thus, a damaged region could develop in the grout adjacent to the alpha-emitting species. This work attempted to recreate this behaviour through nonradioactive 4 He 2+ ion-accelerator experiments, which have not previously been applied to common encapsulation grouts. The influence of this irradiation on a slag-Portland cement was investigated at different ages via transmission electron microscopy energy-dispersive x-ray spectroscopy (TEM-EDX) and supporting techniques, to assess whether 4 He 2+ irradiation caused textural or chemical zonation. No significant changes in hydrate phases or textures were observed, other than minor variations associated with carbonation. This paper provides a proof of concept for using ion acceleration techniques on cements and furthers knowledge on their radiation response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗