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.

SourceResearch connectionFrom yearThrough year

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

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Predicting mesoscale microstructural evolution in electron beam welding

Using the kinetic Monte Carlo simulator, Stochastic Parallel PARticle Kinetic Simulator, from Sandia National Laboratories, a user routine has been developed to simulate mesoscale predictions of a grain structure near a moving heat source. Here, we demonstrate the use of this user routine to produce voxelized, synthetic, three-dimensional microstructures for electron-beam welding by comparing them with experimentally produced microstructures. When simulation input parameters are matched to experimental process parameters, qualitative and quantitative agreement for both grain size and grain morphology are achieved. The method is capable of simulating both single- and multipass welds. As a result, the simulations provide an opportunity for not only accelerated design but also the integration of simulation and experiments in design such that simulations can receive parameter bounds from experiments and, in turn, provide predictions of a resultant microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

First Estimation of Model Parameters for Neutrino-Induced Nucleon Knockout Using Simulation-Based Inference

To enable an accurate determination of oscillation parameters, accelerator-based neutrino experiments require detailed simulations of nuclear interaction physics in the GeV regime. While substantial effort from both theory and experiment is currently being invested to improve the fidelity of these simulations, their present deficiencies typically oblige experimental collaborations to resort to empirical tuning of simulation model parameters. As the precision requirements of the field continue to become more stringent, machine learning techniques may provide a powerful means of handling corresponding growth in the complexity of future neutrino interaction model tuning exercises. To study the suitability of simulation-based inference (SBI) for this physics application, in this paper we revisit a tuned configuration of the GENIE neutrino event generator that was originally developed by the MicroBooNE collaboration. Despite closely reproducing the adopted values of four physics parameters when confronted with the tuned cross-section predictions as input, we find that our trained SBI algorithm prefers modestly different values (within MicroBooNE's assigned uncertainties) and achieves slightly better goodness-of-fit when inference is run on the experimental data set originally used by MicroBooNE. We also find that our trained algorithm can create a fair approximation of an alternative neutrino scattering simulation, NuWro, that shares only a subset of its physics model parameters with GENIE.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗

DOE OSTI2013

PIXE simulation: Models, methods and technologies

The simulation of PIXE (Particle Induced X-ray Emission) is discussed in the context of general-purpose Monte Carlo systems for particle transport. Dedicated PIXE codes are mainly concerned with the application of the technique to elemental analysis, but they lack the capability of dealing with complex experimental configurations. General-purpose Monte Carlo codes provide powerful tools to model the experimental environment in great detail, but so far they have provided limited functionality for PIXE simulation. This paper reviews recent developments that have endowed the Geant4 simulation toolkit with advanced capabilities for PIXE simulation, and related efforts for quantitative validation of cross sections and other physical parameters relevant to PIXE simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2012

DIFFUSIVE ACCELERATION OF PARTICLES AT OBLIQUE, RELATIVISTIC, MAGNETOHYDRODYNAMIC SHOCKS

Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic nuclei and gamma-ray bursts. These sources remain good candidate sites for the generation of ultrahigh energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are germane to the production of relativistic electrons and ions are outlined. The technique employed to identify these characteristics is a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, magnetohydrodynamic (MHD) shocks. Using a compact prescription for diffusion of charges in MHD turbulence, this approach generates particle angular and momentum distributions at any position upstream or downstream of the shock. Simulation output is presented for both small angle and large angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffmann-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat-spectrum generation is addressed; we provide evidence for it being due to shock drift acceleration, a phenomenon well known in heliospheric shock studies. The impact of these theoretical results on blazar science is outlined. Specifically, Fermi Large Area Telescope gamma-ray observations of these relativistic jet sources are providing significant constraints on important environmental quantities for relativistic shocks, namely, the field obliquity, the frequency of scattering, and the level of field turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Fission Product Transport in TRISO Particles and Pebbles

This document demonstrates completion of the goals described in the technical narrative of the FOA project titled: ”Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing” regarding fission product transport in the Kairos-proposed fuel pebble by INL and Kairos Power. Showcased in this report are code developments and simulations in BISON that extend the state of the art in computation and understanding of fission product transport in a TRISO fuel particle and pebble. These enhancements lay the foundation for making predictions of fission product transport that can be used as input in the fuel licensing process. This was achieved by installing existing fuel material models originally used in PARFUME, developing a new failure probability method that is efficient and multi-dimensional, employing material homogenization, and expanding verification and validation simulations to demonstrate the efficacy of the work. All this work is leveraged to spotlight the main deliverable; a three-dimensional model and corresponding demonstration simulation of a pebble, which will serve as the starting point for models used to predict fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Fission Product Transport in TRISO Particles and Pebbles

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

Progress on Beam-Plasma Effect Simulations in Muon Ionization Cooling Lattices

New computational tools are essential for accurate modeling and simulation of the next generation of muon-based accelerators. One of the crucial physics processes specific to muon accelerators that has not yet been simulated in detail is beam-induced plasma effect in liquid, solid, and gaseous absorbers. We report here on the progress of developing the required simulation tools and applying them to study the properties of plasma and its effects on the beam in muon ionization cooling channels.

43 PARTICLE ACCELERATORS↗

DOE OSTI2011

ACCELERATION OF PARTICLES AT THE TERMINATION SHOCK OF A RELATIVISTIC STRIPED WIND

The relativistic wind of obliquely rotating pulsars consists of toroidal stripes of opposite magnetic field polarity, separated by current sheets of hot plasma. By means of two- and three-dimensional particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the termination shock of a relativistic striped wind. At the shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength {lambda} or the wind magnetization {sigma} (in the regime {sigma} >> 1 of magnetically dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of {sigma} with respect to the pre-shock value. The shape of the post-shock spectrum depends primarily on the ratio {lambda}/(r{sub L} {sigma}), where r{sub L} is the relativistic Larmor radius in the wind. The spectrum becomes broader as the value of {lambda}/(r{sub L} {sigma}) increases, passing from a relativistic Maxwellian to a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. Close to the equatorial plane of the wind, where the stripes are symmetric, the highest energy particles resulting from magnetic reconnection can escape ahead of the shock, and be injected into a Fermi-like acceleration process. In the post-shock spectrum, they populate a power-law tail with slope around -2.5, which extends beyond the flat component produced by reconnection. Our study suggests that the spectral break between the radio and the optical band in Pulsar Wind Nebulae can be a natural consequence of particle acceleration at the termination shock of striped pulsar winds.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Maturation of biomass‐to‐biofuels conversion technology pathways for rapid expansion of biofuels production: a system dynamics perspective

Abstract The Biomass Scenario Model ( BSM ) is a system‐dynamics simulation model intended to explore the potential for rapid expansion of the biofuels industry. The model is not predictive — it uses scenario assumptions based on various types of data to simulate industry development, emphasizing how incentives and technological learning‐by‐doing might accelerate industry growth. The BSM simulates major sectors of the biofuels industry, including feedstock production and logistics, conversion, distribution, and end uses, as well as interactions among sectors. The model represents conversion of biomass to biofuels as a set of technology pathways, each of which has allowable feedstocks, capital and operating costs, allowable products, and other defined characteristics. This study and the BSM address bioenergy modeling analytic needs that were identified in recent literature reviews. Simulations indicate that investments are most effective at expanding biofuels production through learning‐by‐doing when they are coordinated with respect to timing, pathway, and target sector within the biofuels industry. Effectiveness metrics include timing and magnitude of increased production, incentive cost and cost effectiveness, and avoidance of windfall profits. Investment costs and optimal investment targets have inherent risks and uncertainties, such as the relative value of investment in more‐mature versus less mature pathways. These can be explored through scenarios, but cannot be precisely predicted. Dynamic competition, including competition for cellulosic feedstocks and ethanol market shares, intensifies during times of rapid growth. Ethanol production increases rapidly, even up to Renewable Fuel Standards‐targeted volumes of biofuel, in simulations that allow higher blending proportions of ethanol in gasoline‐fueled vehicles. Published 2014. This article is a U.S. Government work and is in the public domain in the USA. Biofuels, Bioproducts, Biorefining published by John Wiley & Sons, Ltd on behalf of Society of Chemical Industry.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2021

How does vertical wind shear influence entrainment in squall lines?

The influence of vertical wind shear on updraft entrainment in squall lines is not well understood. To address this knowledge gap, a suite of high-resolution idealized numerical model simulations of squall lines were run in various vertical wind shear (hereafter “shear”) environments to study the effects of shear on entrainment in deep convective updrafts. Low-level horizontal mass flux into the leading edge of the cold pool was strongest in the simulations with the strongest low-level shear. These simulations consequently displayed wider updrafts, less entrainment-driven dilution, and larger buoyancy than the simulations with comparatively weak low-level shear. An analysis of vertical accelerations along trajectories that passed through updrafts showed larger net accelerations from buoyancy in the simulations with stronger low-level shear, which demonstrates how less entrainment-driven dilution equated to stronger updrafts. Additionally, the effects of upper-level shear on entrainment and updraft vertical velocities were generally less pronounced than the effects of low-level shear. We argue that in addition to the outflow boundary-shear interactions and their effect on updraft tilt established by previous authors, decreased entrainment-driven dilution is yet another beneficial effect of strong low-level shear on squall line updraft intensity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Distributions of Particles Accelerated by Strong Alfvénic Turbulence

This work presents a model for generating nonthermal power-law tails of particles’ energy probability density functions in turbulent collisionless plasmas, applicable to both nonrelativistic and relativistic scenarios. We propose that strong Alfvénic turbulence energizes plasma particles through curvature acceleration, particularly for particles with Larmor radii comparable to the scales of turbulence. When the energy density of the energized particles increases, the efficiency of the energy exchange process diminishes. As a result, the acceleration process saturates, leading to power-law distributions of particle momentum and energy. In the nonrelativistic case, the momentum probability density function scales as f(p)dp ∝ p −3 dp, while in the ultrarelativistic case, the energy probability density function scales as f(γ)dγ ∝ γ −3 dγ, where γ is the Lorentz factor. This model provides a unified framework for understanding particle acceleration in both energy regimes, complementing existing analytical approaches. The predicted scalings are consistent with available observations of energetic ion distributions in the heliosphere and with the findings from numerical simulations of ultrarelativistic particle acceleration in magnetically dominated plasma turbulence.

Alfven waves↗

DOE OSTI2015

Investigation of differential surface removal due to electropolishing at JLab

Surface chemistry carried out for Superconducting Radio Frequency (SRF) cavities such as Buffered Chemical Polishing (BCP) and Electropolishing (EP) aims to uniformly remove the internal surface of a cavity along the entire structure and within each cell from equator to iris in order to obtain an equally etched surface. A uniform removal, however, is not readily achievable due to the complex fluid flow and varying temperatures of the acid mixture, which can lead to differential etching. This needs to be considered when envisaging a certain surface damage removal throughout the interior. The process-specific differential etching influences the target frequency set at the manufacturing stage as well as the field flatness and length of the as-built cavity. We report on analyses of JLab's present EP system using experimental data for six nine-cell cavities that have been processed recently in the frame of the LCLS-II high-Q development plan. In conjunction with numerical simulations, the differential etching and the impact on field flatness is assessed.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Simulation study of high-current 7 Li 3+ beam acceleration with laser ion source and RFQ linac using direct plasma injection

Accelerator-based neutron sources (ABNS) utilizing the 7 Li(p,n) 7 Be reaction in inverse kinematics offer strong forward-directed neutron emission, making them attractive for compact and clean neutron source with low background neutrons and unwanted radiation. However, practical use of such systems requires lithium ion beam currents that exceed the capability of conventional ion accelerators by two orders of magnitude. In this study, we numerically designed and evaluated a high-current acceleration system based on a direct plasma injection scheme (DPIS), combining a laser ion source (LIS) and a radio-frequency quadrupole (RFQ) linac. The ion extraction optics and RFQ entrance were optimized using IGUN, OPERA, and GPT, demonstrating that over 1.2 A-class beam current can be injected into the RFQ. The RFQ structure was designed under realistic engineering constraints including surface field strength (Kilpatrick factor) and RF power. GPT simulations showed that a 370 mA 7 Li 3+ beam can be successfully accelerated within a ±10 % energy spread. To validate downstream compatibility, we also designed a simplified medium-energy beam transport (MEBT) section and an interdigital H-type (IH) linac, confirming successful acceleration of a 320 mA beam to the total energy of 14 MeV. These results support the feasibility of the DPIS and RFQ approach as a promising solution for compact neutron sources and other high-intensity ion beam applications.

43 PARTICLE ACCELERATORS↗

NASA NTRS1984

Particle-acceleration by turbulent magnetohydrodynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.

Matthaeus, W. H.↗

DOE OSTI2020

OpenABLext: An automatic code generation framework for agent-based simulations on CPU-GPU-FPGA heterogeneous platforms

The execution of agent-based simulations (ABSs) on hardware accelerator devices such as graphics processing units (GPUs) has been shown to offer great performance potentials. However, in heterogeneous hardware environments, it can become increasingly difficult to find viable partitions of the simulation and provide implementations for different hardware devices. To automate this process, we present OpenABLext, an extension to OpenABL, a model specification language for ABSs. By providing a device-aware OpenCL backend, OpenABLext enables the co-execution of ABS on heterogeneous hardware platforms consisting of central processing units, GPUs, and field programmable gate arrays (FPGAs).We present a novel online dispatching method that efficiently profiles partitions of the simulation during run-time to optimize the hardware assignment while using the profiling results to advance the simulation itself. In addition, OpenABLext features automated conflict resolution based on user-specified rules, supports graph-based simulation spaces, and utilizes an efficient neighbor search algorithm. We show the improved performance of OpenABLext and demonstrate the potential of FPGAs in the context of ABS. We illustrate how co-execution can be used to further lower execution times. OpenABLext can be seen as an enabler to tap the computing power of heterogeneous hardware platforms for ABS.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

The production of ion conics by oblique double layers

Magnetized test ions are subjected to acceleration through a numerically simulated oblique double layer in order to determine whether they emerge with velocity vectors aligned with or oblique to the ambient magnetic field. A criterion for oblique alignment, depending on the double-layer parameters and on the external magnetization, is obtained. When it is applied to observed and theoretical auroral double layers, this criterion predicts that accelerated heavy ions will be substantially less magnetic field aligned than will accelerated hydrogen ions, thus suggesting auroral double layers as a source of high-energy ion conics. Test particle simulations are also used to investigate the perpendicular heating of ions at low altitudes by the electric fields associated with moving auroral arcs. The rapid motion of small-scale structures in the arcs is suggested as a source of low-energy conical ion distributions, and the slow drifts of the entire arc forms are inferred to heat ionospheric ions.

Borovsky, J. E.↗