Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 451 records · Page 25

A Semi-Analytical Approach for State-Space Electromagnetic Transient Simulation

Here, this paper proposes a semi-analytical approach for efficient and accurate electromagnetic transient (EMT) simulation of a power grid. The approach first derives a high-order semi-analytical solution (SAS) of the grid’s state-space EMT model using the differential transformation (DT), and then evaluates the solution over enlarged, variable time steps to significantly accelerate the simulations while maintaining its high accuracy on detailed fast EMT dynamics. The approach also addresses switches during large time steps by using a limit violation detection algorithm with a binary search-enhanced quadratic interpolation. Case studies are conducted on EMT models of the IEEE 39-bus system and large-scale systems to demonstrate the merits of the new simulation approach against traditional numerical methods.

electromagnetic transient↗

Cascade Models in Simulation of Extended Heavy Targets Irradiated by Accelerated Proton and Deuteron Beams

The paper presents a survey of the main numerical models used for simulation of interaction of accelerated particle beams with target nuclei. These models form the core of the software for simulation of various experiments and experimental facilities both for scientific and applied purposes. The beam and target parameters considered in detail in this study (protons and deuterons with energies from 0.66 to 4 AGeV and bulk U targets) cover the range of interest in development of new concepts of nuclear power production aided by accelerated particle beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Autoencoder Neural Network for chemically reacting systems

Incorporating detailed chemical kinetic models is critical for accurate simulations of reacting flows. However, detailed models involve a large number of thermochemical (TC) state variables. Solving the governing equations to evolve these TC variables becomes impractical for real-world applications. In this work, we propose an autoencoder (AE) neural network (NN)-based reduced model to accelerate such simulations. The AE NN is first trained to find a low-dimensional latent representation of the TC states. Then, the evolving state of a chemical system can be tracked by solving the equations of the latent variables instead of the original TC equations. We demonstrate the reduced model in a syngas CO/H 2 combustion system, using training data collected from canonical perfectly stirred reactors (PSRs). It is found that the AE model can reduce the dimension of the combustion system from 12 to 2 while maintaining low reconstruction error and excellent elemental mass conservation for the test dataset. In the a posteriori test, the combustion states obtained from solving the two latent equations are compared to those from solving the 12 equations of the full model. The AE reduced method is found to be able to capture the diverse combustion states on the top two branches of the S-curve well including the extinction turning point, but with higher prediction errors for states near the ignition turning point.

97 MATHEMATICS AND COMPUTING↗

Advanced Computing Annual Report 2025 [Slides]

In Fiscal Year (FY) 2025, the National Laboratory of the Rockies (NLR) continued to advance computing as a cornerstone of energy innovation, expanding the Kestrel high-performance computing (HPC) system to 56 peak petaflops. This growth strengthened Kestrel's role as a national asset for applied energy research, enabling larger, more complex simulations and accelerating the integration of artificial intelligence (AI) methods across the laboratory's computing portfolio. In FY 2025, AI was a component of most projects running on Kestrel, underscoring its central role in modern energy science and engineering. Kestrel supported a broad and diverse set of 507 modeling and simulation projects, engaging 855 researchers across the U.S. Department of Energy's (DOE's) Office of Critical Minerals and Energy Innovation (CMEI) portfolio and other offices, as well as partners from industry, academia, and utilities. These efforts span critical materials discovery, energy systems modeling, grid modernization, advanced manufacturing, and other areas essential to strengthening U.S. energy security and competitiveness. Together, these collaborations produced 708 technical outputs, including 293 peer-reviewed publications, reflecting both the depth and impact of the science enabled by NLR's computing capabilities. This year's report highlights the growing importance and benefit of AI throughout NLR's research programs and features work by early career researchers who are helping shape the future of computing-enabled energy innovation. Explore these sections and the many project successes captured in the pages that follow.

97 MATHEMATICS AND COMPUTING↗

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods

Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrating Experiments, Simulations, and Artificial Intelligence to Accelerate the Discovery of High-Performance Green Composites

The imperative for incorporating greener materials into the aerospace industry necessitates addressing significant challenges associated with the microstructural variability exhibited by recycled and sustainable feedstocks. In this study, we propose an integrated methodology that combines experimental investigations, finite element analysis, and artificial intelligence to develop sustainable composites with consistent properties. Our approach utilizes a pipeline comprising an automated mechanical tester, a finite element method simulator, and a convolutional neural network predictor to identify and optimize fabrication parameters for achieving desired mechanical characteristics in composites. By employing a nested-loop pipeline, our methodology improves sample efficiency, accuracy, and effectively bridges the gap between simulations and real-world performance. This unique methodology offers a promising avenue for facilitating the adoption of aerospace-appropriate green composites.

Athanasiou, Christos↗

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bidimensional Particle-In-Cell simulations for laser-driven proton acceleration using ultra-short, ultra-high contrast laser

In this paper, we report on bi-dimensional Particle-In-Cell simulations performed in order to reproduce the laser-driven proton acceleration obtained when a commercial 200 TW Ti:Sa Laser hits a solid target. The laser-to prepulse contrast was enhanced using plasma mirrors yielding to a main-to-prepulse contrast of ∼10{sup 12}. We varied the pulse duration from 30 fs to 500 fs and the target thickness from 30 nm to several tens of μm. The on-target laser energy was up to 1.8 J leading to an intensity in excess of 10{sup 20 }W cm{sup −2}. A comparison between numerical and existing experimental data [S. Fourmaux et al., Phys. Plasmas 20, 013110 (2013)] is performed, showing a good agreement between experimental results and simulations which confirms that for ultra-thin targets there is an optimum expansion regime. This regime depends on the target thickness and on the laser intensity: if the target is too expanded, the laser travels through the target without being able to deposit its energy within the target. If the target is not sufficiently expanded, the laser energy is reflected by the target. It is important to note that maximum proton energies are reached at longer pulse durations (in the 100 fs regime) than what is currently the best compression pulse length for this type of lasers (typically 20–30 fs). This duration, around 50–100 fs, can be considered a minimum energy transfer time between hot electrons to ions during the considered acceleration process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum Computing for Quantum Simulation for NNSA Mission

Beginning with Feynman in 1982, it was proposed that highly controllable quantum devices, later to become known as quantum computers, would be especially good at the simulation of other quantum systems. Around the turn of the 21st century, this notion was later formalized (by Lloyd and coworkers at MIT and Laflamme and coworkers at Los Alamos) to show how in particular instances, we expect an exponential speedup in the solution of the Schroedinger equation for chemical systems. Subsequently, quantum-accelerated quantum simulation for chemistry techniques have been developed by IBM, Google, Harvard, etc. This effort has led to the OpenFermion (openfermion.org, OF) community, an open source effort for compiling and analyzing quantum algorithms to simulate fermionic systems, including quantum chemistry. Among other functionalities, the current version of OF features data structures and tools for obtaining and manipulating representations of fermionic and qubit Hamiltonians. In October, 2017, the cover of Chemical and Engineering News posited “Quantum computers: Chemistry could be the machine’s killer app” for first principles understanding of biological systems and complex materials. Prof. Aspuru-Guzik (Harvard) predicts “I think that before 2035 QC will be beating classical computers at chemistry.” The algorithms at the heart of quantum chemistry codes (eigensolvers, fast fourier transforms) are the computation bottleneck.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation of cardiovascular response to acceleration stress following weightless exposure

Physiological adjustments taking place during space flight tend to reduce the tolerance of the crew to headward (+Gz) acceleration experienced during the reentry phase of the flight. This reduced tolerance to acceleration stress apparently arises from an adaptation to the microgravity environment of space, including a decrease in the total circulating blood volume. Countermeasures such as anti-g garments have long been known to improve the tolerance to headward g-force, but their effectiveness in space flight has not been fully evaluated. The simulation study presented in this paper is concerned with the response of the cardiovascular system to g-stress following cardiovascular deconditioning, resulting from exposure to weightlessness, or any of its ground-based experimental analogs. The results serve to demonstrate the utility of mathematical modeling and computer simulation for studying the causes of orthostatic intolerance and the remedial measures to lessen it.

Srinivasan, R.↗

SU-F-P-51: Similarity Analysis of the Linear Accelerator Machines Based On Clinical Simulation

Purpose: To evaluate the clinical rationale for Truebeam and Trilogy Linac machines from Varian Medical System as exchangeable treatment modalities in the same radiation oncology department. Methods: Intensity Modulated Radiotherapy (IMRT) plans for different diseases were selected for this study. These disease sites included brain, head and neck, breast, lung, and prostate. The parameters selected for this study were the energy amount, Monitor Unit (MU); dose coverage of target reflected by prescription isodose volume(PIV); dose spillage described by the volume of 50% isodoseline of the prescription; and dose homogeneities represented by the maximum dose (MaxD) and the minimum dose (MinD) of target volume (TV) and critical structure (CS). Each percentage difference between the values of these parameters formed an element of a matrix, which was called Similarity Comparison Matrix(SCM). The elements of the SCM were then simplified by dimensional conversion algorithm, which was used to determine clinical similarity between two machines through a single value. Results: For the selected clinical cases in this study, the average percentage differences between Trilogy and Truebeam in MU was 0.28% with standard deviation(SD) 0.66%, PIV was 0.23% with SD 0.20%, Volume at 50% prescription dose was 0.31% with SD at 0.78%, MaxD at TV is 0.26% with SD 0.35%, MinD at TV is −0.04% with SD 0.51%, MaxD in CS is −0.53% with SD 0.92%, and MinD in CS 3.31%, with SD at 2.89%. The sum, product, geometric and harmonic mean for the matrix elements were 19.0%, 0.00%, 0.19%, and 0.00%. Conclusion: A method to compare two machines in clinical level was developed and some reference values were calculated for decision-making in clinical practice, and this strategy could be expanded to different clinical applications.

60 APPLIED LIFE SCIENCES↗

Multi-modality deep learning for pulse prediction in homogeneous nonlinear systems via parametric conversion

In this Letter, we introduce FusionNet, a multi-modality deep learning framework designed to predict and analyze output pulses in high-power rare-earth-doped laser systems driving parametric conversion in homogeneous guided nonlinear media. FusionNet integrates temporal, spectral, and physical experimental conditions to model ultrafast nonlinear phenomena, including parametric nonlinear frequency conversion, self-phase modulation, and cross-phase modulation in homogeneous guided systems such as gas-filled hollow-core fibers. These systems bridge physical models with experimental data, advancing our understanding of light-guiding principles and nonlinear interactions while expediting the design and optimization of on-demand high-power, high-brightness systems. Our results demonstrate a 73% reduction in prediction error and an 83% improvement in computational efficiency compared to conventional neural networks. This work establishes a new paradigm for accelerating parametric simulations and optimizing experimental designs in high-power laser systems, with further implications for high-precision spectroscopy, quantum information science, and distributed entangled interconnects.

47 OTHER INSTRUMENTATION↗

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

Advanced Simulation Tools for Muon-Based Accelerators

The project developed new software tools incorporating the most accurate theoretical calculations and experimental data available for crucial and not-yet-considered physics processes specific to muon accelerators. These new tools were applied to modeling muon ionization cooling channels for a muon collider or a neutrino factory, the facilities that have the potential to discover and explore new exciting fundamental physics but require the development of demanding technologies and innovative concepts. G4beamline, COSY Infinity, and WARP computer codes are particularly well-suited for addressing muon accelerator modeling problems. G4beamline was chosen for its flexibility and modularity, and for being one of the de facto standard codes used by the muon accelerator community. COSY was selected for speed, ability to take into account high-order nonlinear effects essential for the large beams under consideration, and availability of built-in optimization tools. WARP was chosen for its space charge calculation capabilities, speed and flexibility of implementing new lattice elements such as RF cavities. This final technical report summarizes the efforts on implementing advanced new tools in these computer codes addressing the shortcomings in modeling muon physics processes and using these tools to model muon ionization cooling channels.

43 PARTICLE ACCELERATORS↗

Lagrangian simulations of spray interaction with a surface using a stochastic multi-regime impingement model

Modeling of spray impingement upon the piston or cylinder surface is important for hydrocarbon emission and soot predictions in high-fidelity simulation of internal combustion engines with direct injections. Due to the complexity of the physical process of spray impingement, current computational models heavily rely on the empirical studies which involve large uncertainties. A computational model is presented for spray impingement on solid walls in general computational fluid dynamics (CFD) simulations. The model adopts the Marmanis-Thoroddsen theory for determining the number of secondary droplets and a stochastic method for calculating the velocities of the secondary droplets. The numerical methods and parallel implementation of the impingement model are presented in detail. The model is implemented in Grit, which is a recently developed CPU/GPU performance-portable Lagrangian particle library, for accelerating spray simulations. The model and its implementation are tested for examples of spray impingements on flat and smooth surfaces with different impingement angles.

Ge, Wenjun↗

Characteristics of Escaping Magnetospheric Ions Associated with Magnetic Field Fluctuations

The four Magnetospheric Multiscale (MMS) spacecraft observed a hot ow anomaly (HFA) at the boundary of the quasi-parallel and quasi-perpendicular bow shock and energetic (E > 50 keV) ion bursts exhibitinginverse dispersion in the magnetosheath on 28 December 2015. We consider the possibility that the ions seen both upstream in the foreshock and down-stream in the magnetosheath originate from the magnetosphere. The ion composition ratios, flux levels, and the spectral slopes of the energetic ion energy spectra observed in the region upstream from the bow shock and in themagnetosheath are similar to those in the outer magnetosphere but significantly differ from those seen further upstream from the bow shock at ACE. We can exclude an explanation of the particle source in terms of HFA acceleration. A simulation shows that escaping magnetospheric ions can be scattered and transported across the magnetosheath. Ground magnetometer observations indicate that a solar wind pressure decrease subsequently allows the bow shock to move outward past the MMS spacecraft placing them in the magnetosheath. This was followed by a pressure increase in the magnetosphere. We posit that the enhanced pressure applied to the magnetosphere accelerates the escaping magnetospheric ions (the betatron acceleration), resulting in the inverse dispersion.

S. H. Lee↗

Ultrashort Laser Pulse Shaping and Characterization for Tailored Electron Bunch Generation

Temporally shaped laser pulses are desirable in various applications including emittance reduction and beam-driven acceleration. Pulse shaping techniques enable flexible controls over the longitudinal distribution of electron bunches emitted from the photocathode. While direct manipulation and measurement of an ultrashort pulse can be challenging in the time domain, both actions can be performed in the frequency domain. In this paper, we report the study and development of laser shaper and diagnostics at Argonne Wakefield Accelerator (AWA). Simulations of the shaping process for several sought-after shapes are presented along with the temporal diagnosis. Status of the experiment at the AWA facility is also discussed.

43 PARTICLE ACCELERATORS↗