SEARCH · Engineering Papers

Results for “accelerator modelling and 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗

DOE OSTI2021

Evaluation of Portable Acceleration Solutions for LArTPC Simulation Using Wire-Cell Toolkit

The Liquid Argon Time Projection Chamber (LArTPC) technology plays an essential role in many current and future neutrino experiments. Accurate and fast simulation is critical to developing efficient analysis algorithms and precise physics model projections. The speed of simulation becomes more important as Deep Learning algorithms are getting more widely used in LArTPC analysis and their training requires a large simulated dataset. Heterogeneous computing is an efficient way to delegate computationally intensive tasks to specialized hardware. However, as the landscape of compute accelerators quickly evolves, it becomes increasingly difficult to manually adapt the code to the latest hardware or software environments. A solution which is portable to multiple hardware architectures without substantially compromising performance would thus be very beneficial, especially for long-term projects such as the LArTPC simulations. In search of a portable, scalable and maintainable software solution for LArTPC simulations, we have started to explore high-level portable programming frameworks that support several hardware backends. In this paper, we present our experience porting the LArTPC simulation code in the Wire-Cell Toolkit to NVIDIA GPUs, first with the CUDA programming model and then with a portable library called Kokkos. Preliminary performance results on NVIDIA V100 GPUs and multi-core CPUs are presented, followed by a discussion of the factors affiecting the performance and plans for future improvements.

Yu, Haiwang↗

DOE OSTI2023

Earthquake response of head-mounted equipment in advanced nuclear reactors

The seismic response of safety-related equipment mounted on the head of an advanced reactor, including pumps, control rod drive mechanisms, and reactor monitoring devices, will affect the design and layout of many advanced reactors. High earthquake-induced accelerations in such equipment may challenge their seismic qualification and trigger the need for additional support framing on the reactor head. Base isolation is a design solution that can drastically reduce seismic demands on equipment. This article describes a set of earthquake-simulator experiments conducted on a scale-model of a base-isolated reactor vessel including four representations of head-mounted equipment, with frequencies spanning from 4.5 to 27 Hz. Dynamic responses of the head-mounted equipment, including displacements, accelerations, and strains, were measured in the experiments for three support conditions: conventional, and seismically isolated using single concave Friction Pendulum (SFP) bearings and triple Friction Pendulum (TFP) bearings. Seismic isolation was effective at reducing equipment responses (accelerations, displacements, and strains) with respect to those in the conventionally supported vessel across a range of seismic inputs. Companion numerical studies highlight the accuracy to be expected in the calculation of different response quantities for lightly damped equipment. The importance of characterizing damping in head-mounted, safety-related equipment through physical experiments to support design and risk assessment is made clear through the numerical simulations.

Engineering↗

DOE OSTI2026

The Role of Three-dimensional Effects on Ion Injection and Acceleration in Perpendicular Shocks

Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.

Orusa, Luca [Princeton Univ., NJ (United States); ↗

DOE OSTI2022

Transforming Energy Through Computational Excellence: Bringing Low Mach Number Reactive Flow Simulations at the Exascale

PeleLMeX's unique capabilities are allowing for reactive flow modeling at unprecedented scales and a reasonable time and cost. The code is currently being extended to tackle more practical, design-oriented simulations by implementing Large Eddy Simulation and data-driven chemical models, providing a fast but accurate tool for engineers considering the emergence of GPU-accelerated platforms. These extensions are critical for enabling the physical insight required to design the next generation of combustion devices as a key component of a renewable energy future.

MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗

DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

INTEGRATED WORKFLOW MANAGEMENT FOR PARTICLE ACCELERATOR SIMULATION

Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.

Tourtellott, John↗

DOE OSTI2023

Stellar feedback-regulated black hole growth: driving factors from nuclear to halo scales

ABSTRACT Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are undermassive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multichannel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fuelling: constant thresholds in these properties are typically crossed within ∼0.1 Hubble time of accelerated BH fuelling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold $\Sigma^\star _{1\,\rm kpc}\approx 10^{9.5} \, {\rm M_{\odot }}\,{\rm kpc}^{-2}$, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived thin gas discs, as well as with virialization of the inner circumgalactic medium. The halo mass Mhalo ∼ 1012 M⊙ and stellar mass M* ∼ 1010.5 M⊙ at which BH growth accelerates correspond to ∼L⋆ galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above L⋆.

Astronomy & Astrophysics↗

DOE OSTI2022

Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays

Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Computational Design of Alloys for Energy Technologies

Advanced materials that maintain their mechanical performance under elevated temperatures, corrosive environments, and a range of static and evolving stresses are needed to improve the efficiency and reduce the environmental impact of a wide spectrum of energy technologies. For instance, cost-efficient alloys that can withstand high temperatures (e.g., 700 °C) have a critical role in improving the efficiency and economics of power generation to support decarbonization of the energy sector; such is true of both the nuclear and fossil energy sectors. Considering both the threats of the energy crisis, namely soaring costs of greenhouse gas emission-producing energy and climate change, it is essential to increase the pace of material discovery and enable rapid paths for material qualification to advance clean energy technologies. Conventionally, alloy development has followed a slow Edisonian process that uses repeated cycles of making, characterizing, and modifying to arrive at optimum composition and processing conditions to achieve the desired component performance. This optimization is followed by the necessary stepwise materials qualification. Furthermore, the increasing adoption of sound data management and physics-informed machine learning represents the next step in the acceleration of materials design and development. In the integrated computational materials engineering (ICME) approach, computational modeling and simulation data from different length and time scales can be combined with complex microstructural details from multimodal experimental characterization and selective property testing to close the design loop for rapid alloy development.

Computational Design Of Materials↗

DOE OSTI2024

A real-time energy and cost efficient vehicle route assignment neural recommender system

Here, this paper presents a neural network recommender system algorithm for assigning vehicles to routes based on energy and cost criteria. In this work, we applied this new approach to efficiently identify the most cost-effective medium and heavy duty truck (MDHDT) powertrain technology, from a total cost of ownership (TCO) perspective, for given trips. We employ a machine learning based approach to efficiently estimate the energy consumption of various candidate vehicles over given routes, defined as sequences of links (road segments), with little information known about internal dynamics, i.e. using high level macroscopic route information. A complete recommendation logic is then developed to allow for real-time optimum assignment for each route, subject to the operational constraints of the fleet. We show how this framework can be used to (1) efficiently provide a single trip recommendation with a top-k vehicles star ranking system, and (2) engage in more general assignment problems where n vehicles need to be deployed over m (m ≤ n) trips. This new assignment system has been deployed and integrated into the POLARIS. Transportation System Simulation Tool for use in research conducted by the Department of Energy's Systems and Modeling for Accelerated Research in Transportation (SMART) Mobility Consortium (SMART, 2024).

Energy consumption↗

DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Modeling ablator grain structure impacts in ICF implosions

High-density carbon is a leading ablator material for inertial confinement fusion (ICF). This and some other ablator materials have grain structure which is believed to introduce very small-scale (~nm) density inhomogeneity. In principle, such inhomogeneity can affect key ICF metrics like fuel compression and yield, by, for example, acting as a seed for instabilities and inducing mix between ablator and fuel. However, assessments of such effects are uncertain due to the difficulty of modeling this small-scale structure in ICF simulations, typically requiring reduced-resolution modeling that scales these features. We present a grain model and show both the impact of de-resolving grains and the complex mixing dynamics such structures can induce. We find that different methods for de-resolving grains can yield both different total deposition of kinetic energy perturbations and different fuel–ablator mixing. We then show a simple-to-implement approach for approximately conserving the deposition of perturbed kinetic energy and demonstrate that, for the present grain model and test cases, this approach yields a reasonably matched time history of mix width between less and more resolved grain models. The simulations here also demonstrate the complex interaction history between grain-induced mixing and instability around the fuel–ablator interface, showing, for example, that the grain-induced perturbations typically trigger instability of conduction-driven density gradients in the DT fuel, enhancing mix penetration early in the acceleration of the shell. Furthermore simulating both microscale and nanoscale grains, we find initial evidence for larger mixing in the microscale case of the present model, despite smaller deposited kinetic energy perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗